small optimization for case of haveSecondaries=False, as pointed out …

…by Julien
CRPropa · May 21, 2024 · 14c7835 · 14c7835
1 parent 4b27792
commit 14c7835
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/module/EMInverseComptonScattering.cpp b/src/module/EMInverseComptonScattering.cpp
@@ -189,9 +189,9 @@ void EMInverseComptonScattering::performInteraction(Candidate *candidate) const
 	double Enew = distribution.sample(E, s);
 
 	// add up-scattered photon
-	double Esecondary = E - Enew;
-	double f = Enew / E;
 	if (havePhotons) {
+		double Esecondary = E - Enew;
+		double f = Enew / E;
 		if (random.rand() < pow(1 - f, thinning)) {
 			double w = 1. / pow(1 - f, thinning);
 			Vector3d pos = random.randomInterpolatedPosition(candidate->previous.getPosition(), candidate->current.getPosition());