#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 (EventAction ea, const edm::ParameterSet &ps, CMSSteppingVerbose )

	~TrackingAction () override

Public Attributes
SimActivityRegistry::BeginOfTrackSignal	m_beginOfTrackSignal

SimActivityRegistry::EndOfTrackSignal	m_endOfTrackSignal

Private Attributes
bool	checkTrack_

TrackWithHistory *	currentTrack_

bool	doFineCalo_

double	eMinFine_

EventAction *	eventAction_

TrackInformationExtractor	extractor_

const G4Track *	g4Track_

bool	saveCaloBoundaryInformation_

CMSSteppingVerbose *	steppingVerbose_

Detailed Description

Definition at line 16 of file TrackingAction.h.

Constructor & Destructor Documentation

◆ TrackingAction()

TrackingAction::TrackingAction	(	EventAction *	ea,
		const edm::ParameterSet &	ps,
		CMSSteppingVerbose *	sv
	)

explicit

Definition at line 19 of file TrackingAction.cc.

     : eventAction_(e),
       currentTrack_(nullptr),
       steppingVerbose_(sv),
       g4Track_(nullptr),
       checkTrack_(p.getUntrackedParameter<bool>("CheckTrack", false)),
       doFineCalo_(p.getParameter<bool>("DoFineCalo")),
       saveCaloBoundaryInformation_(p.getParameter<bool>("SaveCaloBoundaryInformation")),
       eMinFine_(p.getParameter<double>("EminFineTrack") * CLHEP::MeV) {}

◆ ~TrackingAction()

TrackingAction::~TrackingAction ( )

override

Definition at line 29 of file TrackingAction.cc.

29 {}

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:34

◆ geant4Track()

const G4Track* TrackingAction::geant4Track ( ) const

inline

Definition at line 25 of file TrackingAction.h.

References g4Track_.

Referenced by StackingAction::ClassifyNewTrack(), and LowEnergyFastSimModel::ModelTrigger().

25 { return g4Track_; }

TrackingAction::g4Track_

const G4Track * g4Track_

Definition: TrackingAction.h:36

◆ 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 68 of file TrackingAction.cc.

References EventAction::addTkCaloStateInfo(), EventAction::addTrack(), TrackWithHistory::checkAtEnd(), checkTrack_, TrackInformation::crossedBoundary(), currentTrack_, doFineCalo_, EgHLTOffHistBins_cfi::et, eventAction_, extractor_, TrackInformation::getMomentumAtBoundary(), TrackInformation::getPositionAtBoundary(), m_endOfTrackSignal, AlCaHLTBitMon_ParallelJobs::p, TrackWithHistory::save(), saveCaloBoundaryInformation_, TrackWithHistory::saved(), TrackWithHistory::setCrossedBoundaryPosMom(), steppingVerbose_, EventAction::trackContainer(), CMSSteppingVerbose::TrackEnded(), and TrackWithHistory::trackID().

Referenced by RunManagerMTWorker::abortEvent().

                                                                  {
   if (eventAction_->trackContainer() != nullptr) {
     uint32_t id = aTrack->GetTrackID();
     math::XYZVectorD pos(aTrack->GetStep()->GetPostStepPoint()->GetPosition().x(),
                          aTrack->GetStep()->GetPostStepPoint()->GetPosition().y(),
                          aTrack->GetStep()->GetPostStepPoint()->GetPosition().z());
     math::XYZTLorentzVectorD mom;
     std::pair<math::XYZVectorD, math::XYZTLorentzVectorD> p(pos, mom);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "PostUserTrackingAction:"
                                    << " aTrack->GetTrackID()=" << aTrack->GetTrackID()
                                    << " currentTrack_->saved()=" << currentTrack_->saved() << " PostStepPosition=("
                                    << pos.x() << "," << pos.y() << "," << pos.z() << ")";
 #endif
 
     if (doFineCalo_) {
       TrackInformation* trkInfo = (TrackInformation*)aTrack->GetUserInformation();
       // 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.
       eventAction_->addTkCaloStateInfo(id, p);
       if (trkInfo->crossedBoundary()) {
         currentTrack_->save();
         currentTrack_->setCrossedBoundaryPosMom(id, trkInfo->getPositionAtBoundary(), trkInfo->getMomentumAtBoundary());
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("DoFineCalo") << "PostUserTrackingAction:"
                                        << " Track " << aTrack->GetTrackID() << " crossed boundary; pos=("
                                        << trkInfo->getPositionAtBoundary().x() << ","
                                        << trkInfo->getPositionAtBoundary().y() << ","
                                        << trkInfo->getPositionAtBoundary().z() << ")"
                                        << " mom[GeV]=(" << trkInfo->getMomentumAtBoundary().x() << ","
                                        << trkInfo->getMomentumAtBoundary().y() << ","
                                        << trkInfo->getMomentumAtBoundary().z() << ","
                                        << trkInfo->getMomentumAtBoundary().e() << ")";
 #endif
       }
     }
 
     TrackInformation* trkInfo = (TrackInformation*)aTrack->GetUserInformation();
     if (extractor_(aTrack).storeTrack() || currentTrack_->saved() ||
         (saveCaloBoundaryInformation_ && trkInfo->crossedBoundary())) {
       if (trkInfo->crossedBoundary()) {
         currentTrack_->setCrossedBoundaryPosMom(id, trkInfo->getPositionAtBoundary(), trkInfo->getMomentumAtBoundary());
       }
       currentTrack_->save();
 
       eventAction_->addTkCaloStateInfo(id, p);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("SimTrackManager")
           << "TrackingAction addTkCaloStateInfo " << id << " of momentum " << mom << " at " << pos;
 #endif
     }
 
     bool withAncestor =
         ((extractor_(aTrack).getIDonCaloSurface() == aTrack->GetTrackID()) || (extractor_(aTrack).isAncestor()));
 
     if (extractor_(aTrack).isInHistory()) {
       // check with end-of-track information
       if (checkTrack_) {
         currentTrack_->checkAtEnd(aTrack);
       }
 
       eventAction_->addTrack(currentTrack_, true, withAncestor);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("SimTrackManager")
           << "TrackingAction addTrack " << currentTrack_->trackID() << " E(GeV)= " << aTrack->GetKineticEnergy() << "  "
           << aTrack->GetDefinition()->GetParticleName() << " added= " << withAncestor << " at "
           << aTrack->GetPosition();
       edm::LogVerbatim("SimTrackManager") << "TrackingAction addTrack " << currentTrack_->trackID() << " added with "
                                           << true << " and " << withAncestor << " at " << pos;
 #endif
 
     } else {
       eventAction_->addTrack(currentTrack_, false, false);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("SimTrackManager")
           << "TrackingAction addTrack " << currentTrack_->trackID() << " added with " << false << " and " << false;
 #endif
 
       delete currentTrack_;
     }
   }
   if (nullptr != steppingVerbose_) {
     steppingVerbose_->TrackEnded(aTrack);
   }
 
   EndOfTrack et(aTrack);
   m_endOfTrackSignal(&et);
 
   currentTrack_ = nullptr;  // reset for next track
 }

◆ PreUserTrackingAction()

void TrackingAction::PreUserTrackingAction ( const G4Track * aTrack )

override

Definition at line 31 of file TrackingAction.cc.

References currentTrack_, doFineCalo_, eventAction_, g4Track_, TrackInformation::isPrimary(), m_beginOfTrackSignal, EventAction::prepareForNewPrimary(), TrackInformation::putInHistory(), TrackWithHistory::save(), steppingVerbose_, and CMSSteppingVerbose::TrackStarted().

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