#include <TrackWithHistory.h>

Public Member Functions
void	checkAtEnd (const G4Track *)

const G4VProcess *	creatorProcess () const

bool	crossedBoundary () const

int	genParticleID () const

int	getIDAtBoundary () const

const math::XYZTLorentzVectorF &	getMomentumAtBoundary () const

const math::XYZTLorentzVectorF &	getPositionAtBoundary () const

double	globalTime () const

double	localTime () const

const math::XYZVectorD &	momentum () const

void	operator delete (void *TrackWithHistory)

void *	operator new (size_t)

int	parentID () const

int	particleID () const

double	properTime () const

void	save ()

bool	saved () const

void	setCrossedBoundaryPosMom (int id, const math::XYZTLorentzVectorF position, const math::XYZTLorentzVectorF momentum)

void	setGenParticleID (int i)

void	setParentID (int i)

void	setTrackID (int i)

bool	storeTrack () const

double	totalEnergy () const

unsigned int	trackID () const

	TrackWithHistory (const G4Track *g4track)

const math::XYZVectorD &	vertexPosition () const

double	weight () const

	~TrackWithHistory ()

Private Member Functions
int	extractGenID (const G4Track *gt) const

Private Attributes
const G4VProcess *	creatorProcess_

bool	crossedBoundary_

int	genParticleID_

double	globalTime_

int	idAtBoundary_

bool	isPrimary_

double	localTime_

math::XYZVectorD	momentum_

math::XYZTLorentzVectorF	momentumAtBoundary_

int	parentID_

int	particleID_

math::XYZTLorentzVectorF	positionAtBoundary_

double	properTime_

bool	saved_

bool	storeTrack_

double	totalEnergy_

unsigned int	trackID_

math::XYZVectorD	vertexPosition_

double	weight_

Detailed Description

The part of the information about a SimTrack that we need from a G4Track

Definition at line 16 of file TrackWithHistory.h.

Constructor & Destructor Documentation

◆ TrackWithHistory()

TrackWithHistory::TrackWithHistory ( const G4Track * g4track )

The constructor is called at PreUserTrackingAction time, when some of the information is not available yet.

Definition at line 17 of file TrackWithHistory.cc.

     : trackID_(0),
       particleID_(0),
       parentID_(0),
       momentum_(math::XYZVectorD(0., 0., 0.)),
       totalEnergy_(0),
       vertexPosition_(math::XYZVectorD(0., 0., 0.)),
       globalTime_(0),
       localTime_(0),
       properTime_(0),
       creatorProcess_(nullptr),
       weight_(0),
       storeTrack_(false),
       saved_(false) {
   if (g4trk != nullptr) {
     TrackInformationExtractor extractor;
     trackID_ = g4trk->GetTrackID();
     particleID_ = G4TrackToParticleID::particleID(g4trk);
     parentID_ = g4trk->GetParentID();
     momentum_ = math::XYZVectorD(g4trk->GetMomentum().x(), g4trk->GetMomentum().y(), g4trk->GetMomentum().z());
     totalEnergy_ = g4trk->GetTotalEnergy();
     vertexPosition_ = math::XYZVectorD(g4trk->GetPosition().x(), g4trk->GetPosition().y(), g4trk->GetPosition().z());
     globalTime_ = g4trk->GetGlobalTime();
     localTime_ = g4trk->GetLocalTime();
     properTime_ = g4trk->GetProperTime();
     creatorProcess_ = g4trk->GetCreatorProcess();
     storeTrack_ = extractor(g4trk).storeTrack();
     saved_ = false;
     crossedBoundary_ = false;
     genParticleID_ = extractGenID(g4trk);
     // V.I. weight is computed in the same way as before
     // without usage of G4Track::GetWeight()
     weight_ = 10000 * genParticleID_;
 #ifdef DEBUG
     LogDebug("TrackInformation") << " TrackWithHistory : created history for " << trackID_ << " with mother "
                                  << parentID_;
 #endif
   }
 }

References creatorProcess_, crossedBoundary_, extractGenID(), genParticleID_, globalTime_, localTime_, LogDebug, momentum_, parentID_, G4TrackToParticleID::particleID(), particleID_, properTime_, saved_, storeTrack_, totalEnergy_, trackID_, vertexPosition_, and weight_.

◆ ~TrackWithHistory()

TrackWithHistory::~TrackWithHistory ( )

inline

Definition at line 22 of file TrackWithHistory.h.

22 {}

Member Function Documentation

◆ checkAtEnd()

void TrackWithHistory::checkAtEnd ( const G4Track * gt )

Internal consistency check (optional). Method called at PostUserTrackingAction time, to check if the information is consistent with that provided to the constructor.

Definition at line 57 of file TrackWithHistory.cc.

                                                    {
   math::XYZVectorD vposdir(gt->GetVertexPosition().x(), gt->GetVertexPosition().y(), gt->GetVertexPosition().z());
   math::XYZVectorD vmomdir(
       gt->GetVertexMomentumDirection().x(), gt->GetVertexMomentumDirection().y(), gt->GetVertexMomentumDirection().z());
   bool ok = true;
   double epsilon = 1.e-6;
   double eps2 = epsilon * epsilon;
   if ((vertexPosition_ - vposdir).Mag2() > eps2) {
     edm::LogWarning("TrackInformation") << "TrackWithHistory vertex position check failed"
                                         << "\nAt construction: " << vertexPosition_ << "\nAt end:          " << vposdir;
     ok = false;
   }
   math::XYZVectorD dirDiff = momentum_.Unit() - vmomdir;
   if (dirDiff.Mag2() > eps2 && momentum_.Unit().R() > eps2) {
     edm::LogWarning("TrackInformation") << "TrackWithHistory momentum direction check failed"
                                         << "\nAt construction: " << momentum_.Unit()
                                         << "\nAt end:          " << vmomdir;
     ok = false;
   }
   if (!ok)
     throw SimG4Exception("TrackWithHistory::checkAtEnd failed");
 }

References geometryDiff::epsilon, submitPVValidationJobs::gt, momentum_, convertSQLiteXML::ok, and vertexPosition_.

Referenced by TrackingAction::PostUserTrackingAction().

◆ creatorProcess()

const G4VProcess* TrackWithHistory::creatorProcess ( ) const

inline

Definition at line 38 of file TrackWithHistory.h.

38 { return creatorProcess_; }

References creatorProcess_.

Referenced by SimTrackManager::getOrCreateVertex().

◆ crossedBoundary()

bool TrackWithHistory::crossedBoundary ( ) const

inline

Definition at line 55 of file TrackWithHistory.h.

55 { return crossedBoundary_; }

References crossedBoundary_.

Referenced by CaloSD::hitBookkeepingFineCalo().

◆ extractGenID()

int TrackWithHistory::extractGenID ( const G4Track * gt ) const

private

Definition at line 80 of file TrackWithHistory.cc.

                                                           {
   void* vgprimary = gt->GetDynamicParticle()->GetPrimaryParticle();
   if (vgprimary == nullptr)
     return -1;
   // replace old-style cast with appropriate new-style cast...
   G4PrimaryParticle* gprimary = (G4PrimaryParticle*)vgprimary;
   GenParticleInfoExtractor ext;
   return ext(gprimary).id();
 }

References submitPVValidationJobs::gt.

Referenced by TrackWithHistory().

◆ genParticleID()

int TrackWithHistory::genParticleID ( ) const

inline

Definition at line 31 of file TrackWithHistory.h.

31 { return genParticleID_; }

References genParticleID_.

Referenced by SimTrackManager::reallyStoreTracks(), and SimTrackManager::resetGenID().

◆ getIDAtBoundary()

int TrackWithHistory::getIDAtBoundary ( ) const

inline

Definition at line 58 of file TrackWithHistory.h.

58 { return idAtBoundary_; }

References idAtBoundary_.

Referenced by CaloSD::hitBookkeepingFineCalo().

◆ getMomentumAtBoundary()

const math::XYZTLorentzVectorF& TrackWithHistory::getMomentumAtBoundary ( ) const

inline

Definition at line 57 of file TrackWithHistory.h.

57 { return momentumAtBoundary_; }

References momentumAtBoundary_.

Referenced by CaloSD::hitBookkeepingFineCalo().

◆ getPositionAtBoundary()

const math::XYZTLorentzVectorF& TrackWithHistory::getPositionAtBoundary ( ) const

inline

Definition at line 56 of file TrackWithHistory.h.

56 { return positionAtBoundary_; }

References positionAtBoundary_.

Referenced by CaloSD::hitBookkeepingFineCalo().

◆ globalTime()

double TrackWithHistory::globalTime ( ) const

inline

Definition at line 35 of file TrackWithHistory.h.

35 { return globalTime_; }

References globalTime_.

Referenced by SimTrackManager::getOrCreateVertex().

◆ localTime()

double TrackWithHistory::localTime ( ) const

inline

Definition at line 36 of file TrackWithHistory.h.

36 { return localTime_; }

References localTime_.

◆ momentum()

const math::XYZVectorD& TrackWithHistory::momentum ( ) const

inline

Definition at line 32 of file TrackWithHistory.h.

32 { return momentum_; }

References momentum_.

Referenced by CaloSD::createNewHit(), SimTrackManager::reallyStoreTracks(), and setCrossedBoundaryPosMom().

◆ operator delete()

void TrackWithHistory::operator delete ( void * TrackWithHistory )

inline

Definition at line 99 of file TrackWithHistory.h.

                                                         {
   fpTrackWithHistoryAllocator->FreeSingle((TrackWithHistory *)aTwH);
 }

References fpTrackWithHistoryAllocator.

◆ operator new()

void * TrackWithHistory::operator new ( size_t )

inline

Definition at line 93 of file TrackWithHistory.h.

                                                   {
   if (!fpTrackWithHistoryAllocator)
     fpTrackWithHistoryAllocator = new G4Allocator<TrackWithHistory>;
   return (void *)fpTrackWithHistoryAllocator->MallocSingle();
 }

References fpTrackWithHistoryAllocator.

◆ parentID()

int TrackWithHistory::parentID ( ) const

inline

Definition at line 30 of file TrackWithHistory.h.

30 { return parentID_; }

References parentID_.

Referenced by SimTrackManager::addTrack(), CaloSD::hitBookkeepingFineCalo(), SimTrackManager::reallyStoreTracks(), and SimTrackManager::saveTrackAndItsBranch().

◆ particleID()

int TrackWithHistory::particleID ( ) const

inline

Definition at line 29 of file TrackWithHistory.h.

29 { return particleID_; }

References particleID_.

Referenced by SimTrackManager::reallyStoreTracks().

◆ properTime()

double TrackWithHistory::properTime ( ) const

inline

Definition at line 37 of file TrackWithHistory.h.

37 { return properTime_; }

References properTime_.

◆ save()

void TrackWithHistory::save ( )

inline

Definition at line 27 of file TrackWithHistory.h.

27 { saved_ = true; }

References saved_.

Referenced by CaloSD::createNewHit(), CaloSD::hitBookkeepingFineCalo(), TrackingAction::PostUserTrackingAction(), TrackingAction::PreUserTrackingAction(), and SimTrackManager::saveTrackAndItsBranch().

◆ saved()

bool TrackWithHistory::saved ( ) const

inline

Definition at line 44 of file TrackWithHistory.h.

44 { return saved_; }

References saved_.

Referenced by CaloSD::hitBookkeepingFineCalo(), and TrackingAction::PostUserTrackingAction().

◆ setCrossedBoundaryPosMom()

void TrackWithHistory::setCrossedBoundaryPosMom	(	int	id,
		const math::XYZTLorentzVectorF	position,
		const math::XYZTLorentzVectorF	momentum
	)

inline

Definition at line 47 of file TrackWithHistory.h.

                                                                        {
     crossedBoundary_ = true;
     idAtBoundary_ = id;
     positionAtBoundary_ = position;
     momentumAtBoundary_ = momentum;
   }