#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_ = nullptr

bool	doFineCalo_

double	eMinFine_

EventAction *	eventAction_

const G4Track *	g4Track_ = nullptr

bool	saveCaloBoundaryInformation_

CMSSteppingVerbose *	steppingVerbose_

TrackInformation *	trkInfo_ = nullptr

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.

References doFineCalo_, eMinFine_, and saveCaloBoundaryInformation_.

     : eventAction_(e),
       steppingVerbose_(sv),
       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 ( )

override

Definition at line 33 of file TrackingAction.cc.

33 {}

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

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

Referenced by RunManagerMTWorker::abortEvent().

                                                                  {
   int id = aTrack->GetTrackID();
   const auto& ppos = aTrack->GetStep()->GetPostStepPoint()->GetPosition();
   math::XYZVectorD pos(ppos.x(), ppos.y(), ppos.z());
   math::XYZTLorentzVectorD mom;
   std::pair<math::XYZVectorD, math::XYZTLorentzVectorD> p(pos, mom);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("DoFineCalo") << "PostUserTrackingAction:"
                                  << " aTrack->GetTrackID()=" << id
                                  << " currentTrack_->saved()=" << currentTrack_->saved();
 #endif
 
   bool boundary = false;
   if (doFineCalo_) {
     // 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);
     boundary = true;
   } else if (trkInfo_->storeTrack() || currentTrack_->saved() ||
              (saveCaloBoundaryInformation_ && trkInfo_->crossedBoundary())) {
     currentTrack_->save();
     eventAction_->addTkCaloStateInfo(id, p);
     boundary = true;
   }
   if (boundary && trkInfo_->crossedBoundary()) {
     currentTrack_->setCrossedBoundaryPosMom(id, trkInfo_->getPositionAtBoundary(), trkInfo_->getMomentumAtBoundary());
     currentTrack_->save();
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "PostUserTrackingAction:"
                                    << " Track " << id << " 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
   }
 
   bool withAncestor = (trkInfo_->getIDonCaloSurface() == id || trkInfo_->isAncestor());
 
   if (trkInfo_->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 " << id << "  "
                                         << aTrack->GetDefinition()->GetParticleName() << " added= " << withAncestor
                                         << " at " << aTrack->GetPosition();
 #endif
 
   } else {
     eventAction_->addTrack(currentTrack_, false, false);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("SimTrackManager") << "TrackingAction addTrack " << id << " added with " << false << " and "
                                         << false;
 #endif
     delete currentTrack_;
   }
   if (nullptr != steppingVerbose_) {
     steppingVerbose_->TrackEnded(aTrack);
   }
 
   EndOfTrack et(aTrack);
   m_endOfTrackSignal(&et);
 }

◆ PreUserTrackingAction()

void TrackingAction::PreUserTrackingAction ( const G4Track * aTrack )

override

Definition at line 35 of file TrackingAction.cc.

References currentTrack_, eMinFine_, eventAction_, g4Track_, TrackInformation::isPrimary(), m_beginOfTrackSignal, EventAction::prepareForNewPrimary(), TrackInformation::putInHistory(), TrackWithHistory::save(), steppingVerbose_, 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 (see NewTrackAction), but are not saved by
   // default. The primary is the earliest ancestor, and it must be saved.
   if (trkInfo_->isPrimary()) {
     eventAction_->prepareForNewPrimary();
     currentTrack_->save();
   }
   if (nullptr != steppingVerbose_) {
     steppingVerbose_->TrackStarted(aTrack, false);
   }
   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();
   }
 }