#include <KillSecondariesTrackAction.h>

Inheritance diagram for KillSecondariesTrackAction:

Public Member Functions
	KillSecondariesTrackAction (edm::ParameterSet const &p)

void	update (const BeginOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...

	~KillSecondariesTrackAction () override

Public Member Functions inherited from SimWatcher
const SimWatcher &	operator= (const SimWatcher &)=delete

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

virtual	~SimWatcher ()

Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()

void	slotForUpdate (const BeginOfTrack * iT)

virtual	~Observer ()

Private Attributes
bool	killHeavy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

Additional Inherited Members

Detailed Description

Definition at line 10 of file KillSecondariesTrackAction.h.

Constructor & Destructor Documentation

◆ KillSecondariesTrackAction()

KillSecondariesTrackAction::KillSecondariesTrackAction ( edm::ParameterSet const & p )

Definition at line 12 of file KillSecondariesTrackAction.cc.

                                                                                {
   killHeavy = p.getParameter<bool>("KillHeavy");
   kmaxIon = p.getParameter<double>("IonThreshold") * MeV;
   kmaxProton = p.getParameter<double>("ProtonThreshold") * MeV;
   kmaxNeutron = p.getParameter<double>("NeutronThreshold") * MeV;
  
   edm::LogInfo("KillSecondaries") << "KillSecondariesTrackAction:: Killing"
                                   << " Flag " << killHeavy << " protons below " << kmaxProton << " MeV, neutrons below "
                                   << kmaxNeutron << " MeV and ions below " << kmaxIon << " MeV\n";
 }

References MeV, and AlCaHLTBitMon_ParallelJobs::p.

◆ ~KillSecondariesTrackAction()

KillSecondariesTrackAction::~KillSecondariesTrackAction ( )

override

Definition at line 23 of file KillSecondariesTrackAction.cc.

23 {}

Member Function Documentation

◆ update()

void KillSecondariesTrackAction::update ( const BeginOfTrack * )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 25 of file KillSecondariesTrackAction.cc.

                                                                {
   if (killHeavy) {
     G4Track *theTrack = (G4Track *)((*trk)());
     TrackInformation *trkInfo = (TrackInformation *)(theTrack->GetUserInformation());
     if (trkInfo) {
       int pdg = theTrack->GetDefinition()->GetPDGEncoding();
       if (!(trkInfo->isPrimary())) {  // Only secondary particles
         double ke = theTrack->GetKineticEnergy() / MeV;
         if ((((pdg / 1000000000 == 1 && ((pdg / 10000) % 100) > 0 && ((pdg / 10) % 100) > 0)) && (ke < kmaxIon)) ||
             ((pdg == 2212) && (ke < kmaxProton)) || ((pdg == 2112) && (ke < kmaxNeutron))) {
           theTrack->SetTrackStatus(fStopAndKill);
           edm::LogInfo("KillSecondaries")
               << "Kill Track " << theTrack->GetTrackID() << " Type " << theTrack->GetDefinition()->GetParticleName()
               << " Kinetic Energy " << ke << " MeV";
         }
       }
     }
   }
 }