#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	ekinMin_

std::vector< double >	ekinMinRegion_

int	endPrintTrackID_

const G4Track *	g4Track_ = nullptr

std::vector< G4Region * >	ptrRegion_

bool	saveCaloBoundaryInformation_

CMSSteppingVerbose *	steppingVerbose_

SimTrackManager *	trackManager_

TrackInformation *	trkInfo_ = nullptr

Detailed Description

Definition at line 19 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_, ekinMin_, ekinMinRegion_, AlCaHLTBitMon_ParallelJobs::p, ptrRegion_, 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")),
       ekinMin_(p.getParameter<double>("PersistencyEmin") * CLHEP::GeV),
       ekinMinRegion_(p.getParameter<std::vector<double>>("RegionEmin")) {
   double eth = p.getParameter<double>("EminFineTrack") * CLHEP::MeV;
   if (doFineCalo_ && eth < ekinMin_) {
     ekinMin_ = eth;
   }
   edm::LogVerbatim("SimG4CoreApplication") << "TrackingAction: boundary: " << saveCaloBoundaryInformation_
                                            << "; DoFineCalo: " << doFineCalo_ << "; ekinMin(MeV)=" << ekinMin_;
   if (!ekinMinRegion_.empty()) {
     ptrRegion_.resize(ekinMinRegion_.size(), nullptr);
   }
 }

◆ ~TrackingAction()

TrackingAction::~TrackingAction ( )

overridedefault

Member Function Documentation

◆ currentTrackWithHistory()

TrackWithHistory* TrackingAction::currentTrackWithHistory ( )

inline

Definition at line 27 of file TrackingAction.h.

References currentTrack_.

27 { return currentTrack_; }

TrackingAction::currentTrack_

TrackWithHistory * currentTrack_

Definition: TrackingAction.h:39

◆ geant4Track()

const G4Track* TrackingAction::geant4Track ( ) const

inline

Definition at line 28 of file TrackingAction.h.

References g4Track_.

Referenced by StackingAction::ClassifyNewTrack().

28 { return g4Track_; }

TrackingAction::g4Track_

const G4Track * g4Track_

Definition: TrackingAction.h:37

◆ getTrackManager()

G4TrackingManager* TrackingAction::getTrackManager ( )

inline

Definition at line 29 of file TrackingAction.h.

29 { return fpTrackingManager; }

◆ PostUserTrackingAction()

void TrackingAction::PostUserTrackingAction ( const G4Track * aTrack )

override

Definition at line 77 of file TrackingAction.cc.

References SimTrackManager::addTrack(), TrackInformation::crossedBoundary(), currentTrack_, doFineCalo_, l1tnanotables_cff::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().

                                                                  {
   // Tracks in history may be upgraded to stored secondary tracks,
   // which cross the boundary between Tracker and Calo
   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() << " proposed to be saved= " << ok
                                      << " 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 39 of file TrackingAction.cc.

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

                                                                 {
   g4Track_ = aTrack;
   currentTrack_ = new TrackWithHistory(aTrack, aTrack->GetParentID());
 
   BeginOfTrack bt(aTrack);
   m_beginOfTrackSignal(&bt);
 
   trkInfo_ = static_cast<TrackInformation*>(aTrack->GetUserInformation());
 
   // Always save primaries
   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 > ekinMin_) {
     // Each track with energy above the threshold should be saved
     trkInfo_->putInHistory();
   }
 }