Corrects the mc loop in GModelSpectralGauss.cpp, in GModelSpectralGauss::mc, as suggested

ellisowen
1 parent 79885c60
Showing 1 changed file with 3 additions and 4 deletions
src/model/GModelSpectralGauss.cpp
@@ -466,15 +466,14 @@ GEnergy GModelSpectralGauss::mc(const GEnergy&amp; emin,
     }
  
     // Get uniform value between 0 and 1
-    double u = ran.uniform();
  
     // Do ...
     do {
         double x1;
         double w;
     	do{
-    		x1        = 2.0 * u - 1.0;
-    		double x2 = 2.0 * u - 1.0;
+    		x1        = 2.0 * ran.uniform() - 1.0;
+    		double x2 = 2.0 * ran.uniform() - 1.0;
     		w         = x1 * x1 + x2 * x2;
     	} while (w >= 1.0);
  
@@ -485,7 +484,7 @@ GEnergy GModelSpectralGauss::mc(const GEnergy&amp; emin,
     	double val = x1 * w;
  
     	// Map into [emin,emax] range
-    	energy = (xmax - xmin) * val + xmin;
+    	energy = m_sigma.value() * val + m_mean.value();
  
     } while ((xmin <= energy) && (energy < xmax));