TrackWithHistory Class Reference

#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

(

const G4Track *

g4track

)

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

Definition at line 16 of file TrackWithHistory.cc.

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

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

TrackWithHistory::~TrackWithHistory ( )

inline

Definition at line 22 of file TrackWithHistory.h.

{}

Member Function Documentation

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 56 of file TrackWithHistory.cc.

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

Referenced by TrackingAction::PostUserTrackingAction().

                                                   {
  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)
    G4Exception("TrackWithHistory::checkAtEnd()", "mc001", FatalException, "check at track end failed");
}

const G4VProcess* TrackWithHistory::creatorProcess ( ) const

inline

Definition at line 38 of file TrackWithHistory.h.

References creatorProcess_.

Referenced by SimTrackManager::getOrCreateVertex().

{ return creatorProcess_; }

bool TrackWithHistory::crossedBoundary ( ) const

inline

Definition at line 55 of file TrackWithHistory.h.

References crossedBoundary_.

Referenced by G4SimTrack::copyCrossedBoundaryVars(), and CaloSD::findBoundaryCrossingParent().

{ return crossedBoundary_; }

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

private

Definition at line 80 of file TrackWithHistory.cc.

Referenced by TrackWithHistory().

                                                          {
  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();
}

int TrackWithHistory::genParticleID ( ) const

inline

Definition at line 31 of file TrackWithHistory.h.

References genParticleID_.

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

{ return genParticleID_; }

int TrackWithHistory::getIDAtBoundary ( ) const

inline

Definition at line 58 of file TrackWithHistory.h.

References idAtBoundary_.

Referenced by G4SimTrack::copyCrossedBoundaryVars().

{ return idAtBoundary_; }

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

inline

Definition at line 57 of file TrackWithHistory.h.

References momentumAtBoundary_.

Referenced by G4SimTrack::copyCrossedBoundaryVars().

{ return momentumAtBoundary_; }

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

inline

Definition at line 56 of file TrackWithHistory.h.

References positionAtBoundary_.

Referenced by G4SimTrack::copyCrossedBoundaryVars().

{ return positionAtBoundary_; }

double TrackWithHistory::globalTime ( ) const

inline

Definition at line 35 of file TrackWithHistory.h.

References globalTime_.

Referenced by SimTrackManager::getOrCreateVertex().

{ return globalTime_; }

double TrackWithHistory::localTime ( ) const

inline

Definition at line 36 of file TrackWithHistory.h.

References localTime_.

{ return localTime_; }

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

inline

Definition at line 32 of file TrackWithHistory.h.

References momentum_.

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

{ return momentum_; }

void TrackWithHistory::operator delete ( void *  TrackWithHistory )

inline

Definition at line 99 of file TrackWithHistory.h.

References fpTrackWithHistoryAllocator.

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

void * TrackWithHistory::operator new ( size_t  )

inline

Definition at line 93 of file TrackWithHistory.h.

References fpTrackWithHistoryAllocator.

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

int TrackWithHistory::parentID ( ) const

inline

Definition at line 30 of file TrackWithHistory.h.

References parentID_.

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

{ return parentID_; }

int TrackWithHistory::particleID ( ) const

inline

Definition at line 29 of file TrackWithHistory.h.

References particleID_.

Referenced by SimTrackManager::reallyStoreTracks().

{ return particleID_; }

double TrackWithHistory::properTime ( ) const

inline

Definition at line 37 of file TrackWithHistory.h.

References properTime_.

{ return properTime_; }

void TrackWithHistory::save ( )

inline

Definition at line 27 of file TrackWithHistory.h.

References saved_.

Referenced by CaloSD::createNewHit(), CaloSD::findBoundaryCrossingParent(), TrackingAction::PostUserTrackingAction(), and TrackingAction::PreUserTrackingAction().

{ saved_ = true; }

bool TrackWithHistory::saved ( ) const

inline

Definition at line 44 of file TrackWithHistory.h.

References saved_.

Referenced by SimTrackManager::cleanTracksWithHistory(), and TrackingAction::PostUserTrackingAction().

{ return saved_; }

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

inline

Definition at line 47 of file TrackWithHistory.h.

References crossedBoundary_, gpuClustering::id, idAtBoundary_, momentum(), momentumAtBoundary_, position, and positionAtBoundary_.

Referenced by TrackingAction::PostUserTrackingAction().

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

void TrackWithHistory::setGenParticleID ( int  i )

inline

Definition at line 42 of file TrackWithHistory.h.

References genParticleID_, and mps_fire::i.

Referenced by SimTrackManager::resetGenID().

{ genParticleID_ = i; }

void TrackWithHistory::setParentID ( int  i )

inline

Definition at line 41 of file TrackWithHistory.h.

References mps_fire::i, and parentID_.

{ parentID_ = i; }

void TrackWithHistory::setTrackID ( int  i )

inline

Definition at line 40 of file TrackWithHistory.h.

References mps_fire::i, and trackID_.

{ trackID_ = i; }

bool TrackWithHistory::storeTrack ( ) const

inline

Definition at line 43 of file TrackWithHistory.h.

References storeTrack_.

{ return storeTrack_; }

double TrackWithHistory::totalEnergy ( ) const

inline

Definition at line 33 of file TrackWithHistory.h.

References totalEnergy_.

Referenced by SimTrackManager::reallyStoreTracks().

{ return totalEnergy_; }

unsigned int TrackWithHistory::trackID ( ) const

inline

Definition at line 28 of file TrackWithHistory.h.

References trackID_.

Referenced by SimTrackManager::addTrack(), SimTrackManager::cleanTracksWithHistory(), SimTrackManager::StrictWeakOrdering::operator()(), trkIDLess::operator()(), TrackingAction::PostUserTrackingAction(), SimTrackManager::reallyStoreTracks(), and CaloSD::update().

{ return trackID_; }

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

inline

Definition at line 34 of file TrackWithHistory.h.

References vertexPosition_.

Referenced by SimTrackManager::getOrCreateVertex().

{ return vertexPosition_; }

double TrackWithHistory::weight ( ) const

inline

Definition at line 39 of file TrackWithHistory.h.

References weight_.

{ return weight_; }

Member Data Documentation

const G4VProcess* TrackWithHistory::creatorProcess_

private

Definition at line 77 of file TrackWithHistory.h.

Referenced by creatorProcess(), and TrackWithHistory().

bool TrackWithHistory::crossedBoundary_

private

Definition at line 83 of file TrackWithHistory.h.

Referenced by crossedBoundary(), setCrossedBoundaryPosMom(), and TrackWithHistory().

int TrackWithHistory::genParticleID_

private

Definition at line 70 of file TrackWithHistory.h.

Referenced by genParticleID(), setGenParticleID(), and TrackWithHistory().

double TrackWithHistory::globalTime_

private

Definition at line 74 of file TrackWithHistory.h.

Referenced by globalTime(), and TrackWithHistory().

int TrackWithHistory::idAtBoundary_

private

Definition at line 84 of file TrackWithHistory.h.

Referenced by getIDAtBoundary(), and setCrossedBoundaryPosMom().

bool TrackWithHistory::isPrimary_

private

Definition at line 82 of file TrackWithHistory.h.

double TrackWithHistory::localTime_

private

Definition at line 75 of file TrackWithHistory.h.

Referenced by localTime(), and TrackWithHistory().

math::XYZVectorD TrackWithHistory::momentum_

private

Definition at line 71 of file TrackWithHistory.h.

Referenced by checkAtEnd(), momentum(), and TrackWithHistory().

math::XYZTLorentzVectorF TrackWithHistory::momentumAtBoundary_

private

Definition at line 86 of file TrackWithHistory.h.

Referenced by getMomentumAtBoundary(), and setCrossedBoundaryPosMom().

int TrackWithHistory::parentID_

private

Definition at line 69 of file TrackWithHistory.h.

Referenced by parentID(), setParentID(), and TrackWithHistory().

int TrackWithHistory::particleID_

private

Definition at line 68 of file TrackWithHistory.h.

Referenced by particleID(), and TrackWithHistory().

math::XYZTLorentzVectorF TrackWithHistory::positionAtBoundary_

private

Definition at line 85 of file TrackWithHistory.h.

Referenced by getPositionAtBoundary(), and setCrossedBoundaryPosMom().

double TrackWithHistory::properTime_

private

Definition at line 76 of file TrackWithHistory.h.

Referenced by properTime(), and TrackWithHistory().

bool TrackWithHistory::saved_

private

Definition at line 80 of file TrackWithHistory.h.

Referenced by save(), saved(), and TrackWithHistory().

bool TrackWithHistory::storeTrack_

private

Definition at line 79 of file TrackWithHistory.h.

Referenced by storeTrack(), and TrackWithHistory().

double TrackWithHistory::totalEnergy_

private

Definition at line 72 of file TrackWithHistory.h.

Referenced by totalEnergy(), and TrackWithHistory().

unsigned int TrackWithHistory::trackID_

private

Definition at line 67 of file TrackWithHistory.h.

Referenced by setTrackID(), trackID(), and TrackWithHistory().

math::XYZVectorD TrackWithHistory::vertexPosition_

private

Definition at line 73 of file TrackWithHistory.h.

Referenced by checkAtEnd(), TrackWithHistory(), and vertexPosition().

double TrackWithHistory::weight_

private

Definition at line 78 of file TrackWithHistory.h.

Referenced by TrackWithHistory(), and weight().