Merge pull request #129 from llebanowski/main_AmBe_Cf

Fix solid angle bug in AmBe and Cf generators
rat-pac · Apr 26, 2024 · 7d82469 · 7d82469
2 parents cffd522 + 7d9a21e
commit 7d82469
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/src/gen/src/AmBeSource.cc b/src/gen/src/AmBeSource.cc
@@ -100,7 +100,7 @@ AmBeSource::AmBeSource() {
     double energy = massNeutron + neutronKE;
     // 	  info << "energy = " << energy*CLHEP::MeV << newline;
     // Generate momentum direction uniformly in phi and cos(theta).
-    double phi = CLHEP::RandFlat::shoot(0., M_PI);
+    double phi = CLHEP::RandFlat::shoot(0., 2.0*M_PI);
     double cosTheta = CLHEP::RandFlat::shoot(-1., 1.);
     double sinTheta = sqrt(1. - cosTheta * cosTheta);
 
@@ -147,7 +147,7 @@ AmBeSource::AmBeSource() {
   for (int nn = 0; nn < Ngamma; nn++) {
     double energy = 4.43;  // from the C12 first excited state
     // Generate momentum direction uniformly in phi and cos(theta).
-    double phi = CLHEP::RandFlat::shoot(0., M_PI);
+    double phi = CLHEP::RandFlat::shoot(0., 2.0*M_PI);
     double cosTheta = CLHEP::RandFlat::shoot(-1., 1.);
     double sinTheta = sqrt(1. - cosTheta * cosTheta);
     double px = energy * sinTheta * cos(phi);

diff --git a/src/gen/src/CfSource.cc b/src/gen/src/CfSource.cc
@@ -132,7 +132,7 @@ CfSource::CfSource(int newIsotope) : Isotope(newIsotope) {
       double neutronKE = fGenerate->shoot() * (fhigh - flow) + flow;
       double energy = massNeutron + neutronKE;
       // Generate momentum direction uniformly in phi and cos(theta).
-      double phi = CLHEP::RandFlat::shoot(0., M_PI);
+      double phi = CLHEP::RandFlat::shoot(0., 2.0*M_PI);
       double cosTheta = CLHEP::RandFlat::shoot(-1., 1.);
       double sinTheta = sqrt(1. - cosTheta * cosTheta);
       double neutronP2 = std::max(0., energy * energy - massNeutron * massNeutron);
@@ -181,7 +181,7 @@ CfSource::CfSource(int newIsotope) : Isotope(newIsotope) {
     for (int nn = 0; nn < Ngamma; nn++) {
       double energy = gGenerate->shoot() * (ghigh - glow) + glow;
       // Generate momentum direction uniformly in phi and cos(theta).
-      double phi = CLHEP::RandFlat::shoot(0., M_PI);
+      double phi = CLHEP::RandFlat::shoot(0., 2.0*M_PI);
       double cosTheta = CLHEP::RandFlat::shoot(-1., 1.);
       double sinTheta = sqrt(1. - cosTheta * cosTheta);
       double px = energy * sinTheta * cos(phi);