#include <KillSecondariesTrackAction.h>

Inheritance diagram for KillSecondariesTrackAction:

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

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

	~KillSecondariesTrackAction ()

Public Member Functions inherited from SimWatcher
	SimWatcher ()

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 ( edm::ParameterSet const & p )

Definition at line 12 of file KillSecondariesTrackAction.cc.

References edm::ParameterSet::getParameter(), and MeV.

                                                                                 {
 
   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 ( )

Definition at line 26 of file KillSecondariesTrackAction.cc.

26 {}

Member Function Documentation

void KillSecondariesTrackAction::update ( const BeginOfTrack * )

virtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 28 of file KillSecondariesTrackAction.cc.

References TrackInformation::isPrimary(), ke, and MeV.

Referenced by Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom().

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