Merge pull request #151 from JeffersonLab/stop_but_not_kill_on_cuts_rtj

Stop but not kill on cuts rtj

Author: markito3

src/GlueXSpecialCuts.cc

@@ -93,7 +93,7 @@ GlueXSpecialCuts::PostStepDoIt( const G4Track& aTrack,
                                  const G4Step&  )
 {
 
-   // Kill the current particle, if requested by G4UserLimits
+   // Stop (but do not kill) the current particle, if requested by G4UserLimits
 
    if (verboseLevel > 1) {
       G4cout << "    " << aTrack.GetParticleDefinition()->GetParticleName()
@@ -107,7 +107,7 @@ GlueXSpecialCuts::PostStepDoIt( const G4Track& aTrack,
    aParticleChange.Initialize(aTrack);
    aParticleChange.ProposeEnergy(0.) ;
    aParticleChange.ProposeLocalEnergyDeposit(aTrack.GetKineticEnergy()) ;
-   aParticleChange.ProposeTrackStatus(fStopAndKill);
+   //aParticleChange.ProposeTrackStatus(fStopAndKill);
    return &aParticleChange;
 }
 

src/GlueXUserEventInformation.cc

@@ -76,6 +76,11 @@ std::map<G4String, G4String> process_4letter_abbrev = {
     {"OpBoundary", "boun"},
     {"Cerenkov", "ckov"},
     {"TPolBeamConversion", "trip"},
+    {"GlueX tracking limits for gammas", "cuts"},
+    {"GlueX tracking limits for electrons", "cuts"},
+    {"GlueX tracking limits for muons", "cuts"},
+    {"GlueX tracking limits for protons", "cuts"},
+    {"GlueX tracking limits for neutrons", "cuts"},
     {"User Limit", "user"}
 };