From 05ae584ff83afd5ae6caf6e605beb81046e4b4b4 Mon Sep 17 00:00:00 2001 From: Richard Jones Date: Sat, 13 Nov 2021 11:50:32 -0500 Subject: [PATCH] * forgot to set the lepton type to muons in handling the control.in option GENBEAM "BHmuons_A...G" [rtj] --- src/GlueXBeamConversionProcess.cc | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/src/GlueXBeamConversionProcess.cc b/src/GlueXBeamConversionProcess.cc index 2d4d3dd..672f5ea 100644 --- a/src/GlueXBeamConversionProcess.cc +++ b/src/GlueXBeamConversionProcess.cc @@ -272,6 +272,7 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name, genbeampars[1] == "bhmuons_a" )) { fStopBeamAfterTarget = 1; + fLeptonPairFamily = 1; if (genbeampars.find(2) != genbeampars.end()) { fBHpair_mass_min = std::atof(genbeampars[2].c_str()); } @@ -284,6 +285,7 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name, genbeampars[1] == "bhmuons_c" )) { fStopBeamAfterTarget = 1; + fLeptonPairFamily = 1; if (genbeampars.find(2) != genbeampars.end()) { fBHpair_mass_min = std::atof(genbeampars[2].c_str()); } @@ -296,6 +298,7 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name, genbeampars[1] == "bhmuons_d" )) { fStopBeamAfterTarget = 1; + fLeptonPairFamily = 1; if (genbeampars.find(2) != genbeampars.end()) { fBHpair_mass_min = std::atof(genbeampars[2].c_str()); } @@ -308,6 +311,7 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name, genbeampars[1] == "bhmuons_e" )) { fStopBeamAfterTarget = 1; + fLeptonPairFamily = 1; if (genbeampars.find(2) != genbeampars.end()) { fBHpair_mass_min = std::atof(genbeampars[2].c_str()); } @@ -320,6 +324,7 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name, genbeampars[1] == "bhmuons_f" )) { fStopBeamAfterTarget = 1; + fLeptonPairFamily = 1; if (genbeampars.find(2) != genbeampars.end()) { fBHpair_mass_min = std::atof(genbeampars[2].c_str()); } @@ -332,6 +337,7 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name, genbeampars[1] == "bhmuons_g" )) { fStopBeamAfterTarget = 1; + fLeptonPairFamily = 1; if (genbeampars.find(2) != genbeampars.end()) { fBHpair_mass_min = std::atof(genbeampars[2].c_str()); }