#include <TrackingAction.h>

Inheritance diagram for TrackingAction:

Public Member Functions
TrackWithHistory *	currentTrackWithHistory ()

const G4Track *	geant4Track () const

G4TrackingManager *	getTrackManager ()

void	PostUserTrackingAction (const G4Track *aTrack) override

void	PreUserTrackingAction (const G4Track *aTrack) override

	TrackingAction (SimTrackManager , CMSSteppingVerbose , const edm::ParameterSet &ps)

	~TrackingAction () override=default

Public Attributes
SimActivityRegistry::BeginOfTrackSignal	m_beginOfTrackSignal

SimActivityRegistry::EndOfTrackSignal	m_endOfTrackSignal

Private Attributes
bool	checkTrack_

TrackWithHistory *	currentTrack_ = nullptr

bool	doFineCalo_

double	eMinFine_

int	endPrintTrackID_

const G4Track *	g4Track_ = nullptr

bool	saveCaloBoundaryInformation_

CMSSteppingVerbose *	steppingVerbose_

SimTrackManager *	trackManager_

TrackInformation *	trkInfo_ = nullptr

Detailed Description

Definition at line 16 of file TrackingAction.h.

Constructor & Destructor Documentation

◆ TrackingAction()

TrackingAction::TrackingAction	(	SimTrackManager *	stm,
		CMSSteppingVerbose *	sv,
		const edm::ParameterSet &	ps
	)

explicit

Definition at line 19 of file TrackingAction.cc.

References doFineCalo_, eMinFine_, and saveCaloBoundaryInformation_.

     : trackManager_(stm),
       steppingVerbose_(sv),
       endPrintTrackID_(p.getParameter<int>("EndPrintTrackID")),
       checkTrack_(p.getUntrackedParameter<bool>("CheckTrack", false)),
       doFineCalo_(p.getParameter<bool>("DoFineCalo")),
       saveCaloBoundaryInformation_(p.getParameter<bool>("SaveCaloBoundaryInformation")),
       eMinFine_(p.getParameter<double>("EminFineTrack") * CLHEP::MeV) {
   if (!doFineCalo_) {
     eMinFine_ = DBL_MAX;
   }
   edm::LogVerbatim("SimG4CoreApplication") << "TrackingAction: boundary: " << saveCaloBoundaryInformation_
                                            << "; DoFineCalo: " << doFineCalo_ << "; EminFineTrack(MeV)=" << eMinFine_;
 }

◆ ~TrackingAction()

TrackingAction::~TrackingAction ( )

overridedefault

Member Function Documentation

◆ currentTrackWithHistory()

TrackWithHistory* TrackingAction::currentTrackWithHistory ( )

inline

Definition at line 24 of file TrackingAction.h.

References currentTrack_.

24 { return currentTrack_; }

TrackingAction::currentTrack_

TrackWithHistory * currentTrack_

Definition: TrackingAction.h:36

◆ geant4Track()

const G4Track* TrackingAction::geant4Track ( ) const

inline

Definition at line 25 of file TrackingAction.h.

References g4Track_.

Referenced by StackingAction::ClassifyNewTrack().

25 { return g4Track_; }

TrackingAction::g4Track_

const G4Track * g4Track_

Definition: TrackingAction.h:34

◆ getTrackManager()

G4TrackingManager* TrackingAction::getTrackManager ( )

inline

Definition at line 26 of file TrackingAction.h.

26 { return fpTrackingManager; }

◆ PostUserTrackingAction()

void TrackingAction::PostUserTrackingAction ( const G4Track * aTrack )

override

Definition at line 75 of file TrackingAction.cc.

References SimTrackManager::addTrack(), TrackInformation::crossedBoundary(), currentTrack_, doFineCalo_, EgHLTOffHistBins_cfi::et, TrackInformation::getIDonCaloSurface(), TrackInformation::getMomentumAtBoundary(), TrackInformation::getPositionAtBoundary(), TrackInformation::isAncestor(), TrackInformation::isInHistory(), m_endOfTrackSignal, convertSQLiteXML::ok, saveCaloBoundaryInformation_, TrackWithHistory::saved(), TrackWithHistory::setCrossedBoundaryPosMom(), TrackWithHistory::setToBeSaved(), TrackInformation::storeTrack(), trackManager_, and trkInfo_.

Referenced by RunManagerMTWorker::abortEvent().

                                                                  {
   // Add the post-step position for every track in history to the TrackManager.
   // Tracks in history may be upgraded to stored tracks, at which point
   // the post-step position is needed again.
   int id = aTrack->GetTrackID();
   bool ok = (trkInfo_->storeTrack() || currentTrack_->saved());
   if (trkInfo_->crossedBoundary()) {
     currentTrack_->setCrossedBoundaryPosMom(id, trkInfo_->getPositionAtBoundary(), trkInfo_->getMomentumAtBoundary());
     ok = (ok || saveCaloBoundaryInformation_ || doFineCalo_);
   }
   if (ok) {
     currentTrack_->setToBeSaved();
   }
 
   bool withAncestor = (trkInfo_->getIDonCaloSurface() == id || trkInfo_->isAncestor());
   bool isInHistory = trkInfo_->isInHistory();
 
   trackManager_->addTrack(currentTrack_, aTrack, isInHistory, withAncestor);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("TrackingAction") << "TrackingAction end track=" << id << "  "
                                      << aTrack->GetDefinition()->GetParticleName() << " ansestor= " << withAncestor
                                      << " saved= " << currentTrack_->saved() << " end point " << aTrack->GetPosition();
 #endif
   if (!isInHistory) {
     delete currentTrack_;
   }
 
   EndOfTrack et(aTrack);
   m_endOfTrackSignal(&et);
 }

◆ PreUserTrackingAction()

void TrackingAction::PreUserTrackingAction ( const G4Track * aTrack )

override

Definition at line 34 of file TrackingAction.cc.

References SimTrackManager::cleanTracksWithHistory(), currentTrack_, eMinFine_, endPrintTrackID_, g4Track_, TrackInformation::isPrimary(), m_beginOfTrackSignal, TrackInformation::putInHistory(), TrackWithHistory::setToBeSaved(), steppingVerbose_, CMSSteppingVerbose::stopEventPrint(), trackManager_, CMSSteppingVerbose::trackStarted(), and trkInfo_.

                                                                 {
   g4Track_ = aTrack;
   currentTrack_ = new TrackWithHistory(aTrack);
 
   BeginOfTrack bt(aTrack);
   m_beginOfTrackSignal(&bt);
 
   trkInfo_ = static_cast<TrackInformation*>(aTrack->GetUserInformation());
 
   // Always save primaries
   // Decays from primaries are marked as primaries, but are not saved by
   // default. The primary is the earliest ancestor, and it must be saved.
   if (trkInfo_->isPrimary()) {
     trackManager_->cleanTracksWithHistory();
     currentTrack_->setToBeSaved();
   }
   if (nullptr != steppingVerbose_) {
     steppingVerbose_->trackStarted(aTrack, false);
     if (aTrack->GetTrackID() == endPrintTrackID_) {
       steppingVerbose_->stopEventPrint();
     }
   }
   double ekin = aTrack->GetKineticEnergy();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("DoFineCalo") << "PreUserTrackingAction: Start processing track " << aTrack->GetTrackID()
                                  << " pdgid=" << aTrack->GetDefinition()->GetPDGEncoding()
                                  << " ekin[GeV]=" << ekin / CLHEP::GeV << " vertex[cm]=("
                                  << aTrack->GetVertexPosition().x() / CLHEP::cm << ","
                                  << aTrack->GetVertexPosition().y() / CLHEP::cm << ","
                                  << aTrack->GetVertexPosition().z() / CLHEP::cm << ")"
                                  << " parentid=" << aTrack->GetParentID();
 #endif
   if (ekin > eMinFine_) {
     // It is impossible to tell whether daughter tracks if this track may need to be saved at
     // this point; Therefore, every track above the threshold is put in history,
     // so that it can potentially be saved later.
     trkInfo_->putInHistory();
   }
 }