#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
virtual void	beginRun (edm::EventSetup const &)

bool	isMT () const

const SimWatcher &	operator= (const SimWatcher &)=delete

virtual void	registerConsumes (edm::ConsumesCollector)

	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
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

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.

References AlCaHLTBitMon_ParallelJobs::p.

                                                                                {
   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";
 }

◆ ~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.

References TrackInformation::isPrimary().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                                {
   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";
         }
       }
     }
   }
 }