iGeant4::Geant4TruthIncident Class Reference

#include <Geant4TruthIncident.h>

Inheritance diagram for iGeant4::Geant4TruthIncident:

Collaboration diagram for iGeant4::Geant4TruthIncident:

Public Member Functions
	Geant4TruthIncident (const G4Step *, const ISF::ISFParticle &baseISP, AtlasDetDescr::AtlasRegion geoID, AtlasG4EventUserInfo *atlasG4EvtUserInfo)
virtual	~Geant4TruthIncident ()
const HepMC::FourVector &	position () const override final
	Return HepMC position of the truth vertex. More...
int	physicsProcessCategory () const override final
	Return category of the physics process represented by the truth incident (eg hadronic, em, ..) More...
int	physicsProcessCode () const override final
	Return specific physics process code of the truth incident (eg ionisation, bremsstrahlung, ..) More...
double	parentP2 () const override final
	Return p^2 of the parent particle. More...
double	parentPt2 () const override final
	Return pT^2 of the parent particle. More...
double	parentEkin () const override final
	Return Ekin of the parent particle. More...
int	parentPdgCode () const override final
	Return the PDG Code of the parent particle. More...
int	parentBarcode () override final
	Return the barcode of the parent particle. More...
int	parentUniqueID () override final
	Return the unique ID of the parent particle. More...
int	parentStatus () override final
	Return the status of the parent particle. More...
bool	parentSurvivesIncident () const override final
	Return a boolean whether or not the parent particle survives the incident. More...
HepMC::GenParticlePtr	parentParticleAfterIncident (int newBC) override final
	Return the parent particle after the TruthIncident vertex (and give it a new barcode) More...
const G4ThreeVector	childP (unsigned short index) const
	Return p of the i-th child particle. More...
double	childP2 (unsigned short index) const override final
	Return p^2 of the i-th child particle. More...
double	childPt2 (unsigned short index) const override final
	Return pT^2 of the i-th child particle. More...
double	childEkin (unsigned short index) const override final
	Return Ekin of the i-th child particle. More...
int	childPdgCode (unsigned short index) const override final
	Return the PDG Code of the i-th child particle. More...
int	childBarcode (unsigned short index) const override final
	Return the barcode of the i-th child particle (if defined as part of the TruthIncident) otherwise return 0. More...
ISF::InteractionClass_t	interactionClassification () const override final
	The interaction classifications are described as follows: STD_VTX: interaction of a particle without a pre-defined decay; QS_SURV_VTX: a particle with a pre-defined decay under-going a non-destructive interaction; QS_DEST_VTX: a particle with a pre-defined decay under-going a destructive interaction other than its pre-defined decay; QS_PREDEF_VTX: a particle under-going its pre-defined decay. More...
HepMC::GenParticlePtr	parentParticle () override final
	Return the parent particle as a HepMC particle type. More...
HepMC::GenParticlePtr	childParticle (unsigned short index, int bc) override final
	Return the i-th child as a HepMC particle type and assign the given Barcode to the simulator particle. More...
AtlasDetDescr::AtlasRegion	geoID ()
	Return the SimGeoID corresponding to the vertex of the truth incident. More...
unsigned short	numberOfChildren () const
	Return total number of child particles. More...
bool	childrenP2Pass (double p2cut)
	Return true if at least one child particle passes the given p^2 cut (= at least one child with p^2 >= pt2cut) More...
bool	childrenPt2Pass (double pt2cut)
	Return true if at least one child particle passes the given pT^2 cut (= at least one child with pT^2 >= pt2cut) More...
bool	childrenEkinPass (double ekincut)
	Return true if at least one child particle passes the given Ekin cut (= at least one child with Ekin >= ekincut) More...
void	setChildPassedFilters (unsigned short index)
	Record that a particular child passed a check. More...
bool	childPassedFilters (unsigned short index) const
	Should a particular child be written out to the GenEvent. More...
void	setPassWholeVertices (bool passWholeVertex)
	Set whether this TruthIncident should pass the vertex as a whole or individual children. More...

Protected Attributes
int	m_numChildren
bool	m_passWholeVertex
std::vector< bool >	m_childPassedFilters

Private Member Functions
	Geant4TruthIncident ()
void	prepareChildren ()
	prepare the child particles More...
bool	particleAlive (const G4Track *track) const
	check if the given G4Track represents a particle that is alive in ISF or ISF-G4 More...
HepMC::GenParticlePtr	convert (const G4Track *particle, const int barcode, const bool secondary) const

Private Attributes
bool	m_positionSet
HepMC::FourVector	m_position
const G4Step *	m_step {}
const ISF::ISFParticle &	m_baseISP
AtlasG4EventUserInfo *m_atlasG4EvtUserInfo	ATLAS_THREAD_SAFE {}
std::vector< const G4Track * >	m_children
HepMC::GenParticlePtr	m_parentParticleAfterIncident {}
AtlasDetDescr::AtlasRegion	m_geoID
	region that the TruthIncident is located in More...

Detailed Description

ISF_Geant4 specific implementation of the ISF::ITruthIncident

Author

Andre.nosp@m.as.S.nosp@m.chael.nosp@m.icke.nosp@m.@cern.nosp@m..ch

Definition at line 44 of file Geant4TruthIncident.h.

Constructor & Destructor Documentation

Geant4TruthIncident() [1/2]

iGeant4::Geant4TruthIncident::Geant4TruthIncident	(	const G4Step *	step,
		const ISF::ISFParticle &	baseISP,
		AtlasDetDescr::AtlasRegion	geoID,
		AtlasG4EventUserInfo *	atlasG4EvtUserInfo
	)

Definition at line 61 of file Geant4TruthIncident.cxx.

                                                                                         :
  ITruthIncident(geoID, step->GetSecondaryInCurrentStep()->size()), // switch to G4Step::GetNumberOfSecondariesInCurrentStep() once we're using G4 10.2 or later
  m_positionSet(false),
  m_position(),
  m_step(step),
  m_baseISP(baseISP),
  m_atlasG4EvtUserInfo(atlasG4EvtUserInfo)
{
  // prepare children:
  prepareChildren();
 
  // calculate position:
  const G4StepPoint *postStepPoint = m_step->GetPostStepPoint();
  const G4ThreeVector         &pos = postStepPoint->GetPosition();
  const G4double              time = postStepPoint->GetGlobalTime()*Gaudi::Units::c_light;
  m_position.set( pos.x(), pos.y(), pos.z(), time );
}

~Geant4TruthIncident()

virtual iGeant4::Geant4TruthIncident::~Geant4TruthIncident ( )

inlinevirtual

Definition at line 50 of file Geant4TruthIncident.h.

{};

Geant4TruthIncident() [2/2]

iGeant4::Geant4TruthIncident::Geant4TruthIncident ( )

private

Member Function Documentation

childBarcode()

int iGeant4::Geant4TruthIncident::childBarcode ( unsigned short  index ) const

finaloverridevirtual

Return the barcode of the i-th child particle (if defined as part of the TruthIncident) otherwise return 0.

Implements ISF::ITruthIncident.

Definition at line 211 of file Geant4TruthIncident.cxx.

                                                                        {
  // the G4Track instance for the current child particle
  const G4Track* track = m_children[index];
  // This should be a *secondary* track.  If it has a primary, it was a decay and
  //  we are running with quasi-stable particle simulation.  Note that if the primary
  //  track is passed in as a secondary that survived the interaction, then this was
  //  *not* a decay and we should not treat it in this way
  if (track->GetDynamicParticle() &&
      track->GetDynamicParticle()->GetPrimaryParticle() &&
      track->GetDynamicParticle()->GetPrimaryParticle()->GetUserInformation()){
    // Then the new particle should use the same barcode as the old one!!
    PrimaryParticleInformation* primaryPartInfo = dynamic_cast<PrimaryParticleInformation*>( track->GetDynamicParticle()->GetPrimaryParticle()->GetUserInformation() );
    return primaryPartInfo->GetParticleBarcode();
  }
  return 0;
}

childEkin()

double iGeant4::Geant4TruthIncident::childEkin ( unsigned short  index ) const

finaloverridevirtual

Return Ekin of the i-th child particle.

Implements ISF::ITruthIncident.

Definition at line 203 of file Geant4TruthIncident.cxx.

                                                                   {
  return (m_children[i]->GetKineticEnergy()/CLHEP::MeV);
}

childP()

const G4ThreeVector iGeant4::Geant4TruthIncident::childP

(

unsigned short

index

)

const

Return p of the i-th child particle.

Definition at line 193 of file Geant4TruthIncident.cxx.

                                                                             {
  const G4ThreeVector & mom= m_children[i]->GetMomentum();
  return mom;
}

childP2()

double iGeant4::Geant4TruthIncident::childP2 ( unsigned short  index ) const

finaloverridevirtual

Return p^2 of the i-th child particle.

Implements ISF::ITruthIncident.

Definition at line 188 of file Geant4TruthIncident.cxx.

                                                                 {
  const G4ThreeVector & mom= m_children[i]->GetMomentum();
  return mom.mag2();
}

childParticle()

HepMC::GenParticlePtr iGeant4::Geant4TruthIncident::childParticle ( unsigned short  index, int  bc  )

finaloverridevirtual

Return the i-th child as a HepMC particle type and assign the given Barcode to the simulator particle.

Implements ISF::ITruthIncident.

Definition at line 228 of file Geant4TruthIncident.cxx.

                                                                                  {
  // the G4Track instance for the current child particle
  const G4Track* thisChildTrack = m_children[i];
 
  // NB: NOT checking if secondary is actually alive. Even with zero momentum,
  //     secondary could decay right away and create further particles which pass the
  //     truth strategies.
 
  HepMC::GenParticlePtr hepParticle = convert( thisChildTrack, newBarcode, true );
  TrackHelper tHelper(thisChildTrack);
  TrackInformation *trackInfo = tHelper.GetTrackInformation();
 
  // needed to make AtlasG4 work with ISF TruthService
  if (trackInfo==nullptr) {
    trackInfo = new TrackInformation( hepParticle );
    thisChildTrack->SetUserInformation( trackInfo );
  }
 
  trackInfo->SetCurrentGenParticle(hepParticle);
  trackInfo->SetClassification(TrackInformation::RegisteredSecondary);
  trackInfo->SetRegenerationNr(0);
 
  return hepParticle;
}

childPassedFilters()

bool ISF::ITruthIncident::childPassedFilters ( unsigned short  index ) const

inlineinherited

Should a particular child be written out to the GenEvent.

Definition at line 186 of file ITruthIncident.h.

                                                                       {
    return m_childPassedFilters[index];
  }

childPdgCode()

int iGeant4::Geant4TruthIncident::childPdgCode ( unsigned short  index ) const

finaloverridevirtual

Return the PDG Code of the i-th child particle.

Implements ISF::ITruthIncident.

Definition at line 207 of file Geant4TruthIncident.cxx.

                                                                   {
  return m_children[i]->GetDefinition()->GetPDGEncoding();
}

childPt2()

double iGeant4::Geant4TruthIncident::childPt2 ( unsigned short  index ) const

finaloverridevirtual

Return pT^2 of the i-th child particle.

Implements ISF::ITruthIncident.

Definition at line 198 of file Geant4TruthIncident.cxx.

                                                                  {
  const G4ThreeVector & mom= m_children[i]->GetMomentum();
  return mom.perp2();
}

childrenEkinPass()

bool ISF::ITruthIncident::childrenEkinPass ( double  ekincut )

inlineinherited

Return true if at least one child particle passes the given Ekin cut (= at least one child with Ekin >= ekincut)

Definition at line 170 of file ITruthIncident.h.

                                                         {
    bool pass = false; // true if cut passed
    // as soon as at a particle passes the cut -> end loop and return true
    for ( unsigned short i=0; !(pass && m_passWholeVertex) && (i<m_numChildren); ++i) {
      bool thispassed = (childEkin(i) >= ekincut);
      if(thispassed) { setChildPassedFilters(i); }
      pass |= thispassed;
    }
    return pass;
  }

childrenP2Pass()

bool ISF::ITruthIncident::childrenP2Pass ( double  p2cut )

inlineinherited

Return true if at least one child particle passes the given p^2 cut (= at least one child with p^2 >= pt2cut)

Definition at line 146 of file ITruthIncident.h.

                                                     {
    bool pass = false; // true if cut passed
    // as soon as at a particle passes the cut -> end loop and return true
    for ( unsigned short i=0; !(pass && m_passWholeVertex) && (i<m_numChildren); ++i) {
      bool thispassed = (childP2(i) >= p2cut);
      if(thispassed) { setChildPassedFilters(i); }
      pass |= thispassed;
    }
    return pass;
  }

childrenPt2Pass()

bool ISF::ITruthIncident::childrenPt2Pass ( double  pt2cut )

inlineinherited

Return true if at least one child particle passes the given pT^2 cut (= at least one child with pT^2 >= pt2cut)

Definition at line 158 of file ITruthIncident.h.

                                                       {
    bool pass = false; // true if cut passed
    // as soon as at a particle passes the cut -> end loop and return true
    for ( unsigned short i=0; !(pass && m_passWholeVertex) && (i<m_numChildren); ++i) {
      bool thispassed = (childPt2(i) >= pt2cut);
      if(thispassed) { setChildPassedFilters(i); }
      pass |= thispassed;
    }
    return pass;
  }

convert()

HepMC::GenParticlePtr iGeant4::Geant4TruthIncident::convert ( const G4Track *  particle, const int  barcode, const bool  secondary  ) const

private

Definition at line 276 of file Geant4TruthIncident.cxx.

                                                                                                                           {
#endif
 
  const G4ThreeVector & mom =  track->GetMomentum();
  const double energy =  track->GetTotalEnergy();
  const int pdgCode = track->GetDefinition()->GetPDGEncoding();
  const HepMC::FourVector fourMomentum( mom.x(), mom.y(), mom.z(), energy);
 
  const HepMC::GenParticlePtr parent = m_atlasG4EvtUserInfo->GetCurrentGenParticle();
  int status = (secondary) ? 1 + HepMC::SIM_STATUS_THRESHOLD : parent->status() + HepMC::SIM_STATUS_INCREMENT;
  // Treat child particles of pre-defined decays differently
  if (this->interactionClassification() == ISF::QS_PREDEF_VTX) {
    const G4DynamicParticle* dynPart = track->GetDynamicParticle();
    bool hasPredefinedDecay = (dynPart && (nullptr!=(dynPart->GetPreAssignedDecayProducts())));
    status = (hasPredefinedDecay)? 2 : 1;
  }
  // Treat particles with a pre-defined decay which have survived an
  // interaction differently for now.
  if (this->interactionClassification() == ISF::QS_SURV_VTX) {
    status = status%HepMC::SIM_STATUS_INCREMENT;
    // Such particles were previously assigned barcodes below
    // HepMC::SIM_BARCODE_OFFSET and therefore
    // HepMC::is_simulation_particle(barcode) would return false. Not
    // applying HepMC::SIM_STATUS_INCREMENT means that
    // HepMC::is_simulation_particle(status) will also return false.
    // TODO This should be revisited in the future should we decide
    // that we want to be able to easily identify such particles.
  }
  HepMC::GenParticlePtr newParticle = HepMC::newGenParticlePtr(fourMomentum, pdgCode, status);
 
#ifndef HEPMC3
  // This should be a *secondary* track.  If it has a primary, it was a decay and
  //  we are running with quasi-stable particle simulation.  Note that if the primary
  //  track is passed in as a secondary that survived the interaction, then this was
  //  *not* a decay and we should not treat it in this way
  if (secondary &&
      track->GetDynamicParticle() &&
      track->GetDynamicParticle()->GetPrimaryParticle() &&
      track->GetDynamicParticle()->GetPrimaryParticle()->GetUserInformation()){
    // Then the new particle should use the same barcode as the old one!!
    PrimaryParticleInformation* primaryPartInfo = dynamic_cast<PrimaryParticleInformation*>( track->GetDynamicParticle()->GetPrimaryParticle()->GetUserInformation() );
    HepMC::suggest_barcode( newParticle, primaryPartInfo->GetParticleBarcode() );
  } else {
    HepMC::suggest_barcode( newParticle, barcode );
  }
#endif
 
  return newParticle;
}

geoID()

AtlasDetDescr::AtlasRegion ISF::ITruthIncident::geoID ( )

inlineinherited

Return the SimGeoID corresponding to the vertex of the truth incident.

Definition at line 53 of file ITruthIncident.h.

{ return m_geoID; };

interactionClassification()

ISF::InteractionClass_t iGeant4::Geant4TruthIncident::interactionClassification ( ) const

finaloverridevirtual

The interaction classifications are described as follows: STD_VTX: interaction of a particle without a pre-defined decay; QS_SURV_VTX: a particle with a pre-defined decay under-going a non-destructive interaction; QS_DEST_VTX: a particle with a pre-defined decay under-going a destructive interaction other than its pre-defined decay; QS_PREDEF_VTX: a particle under-going its pre-defined decay.

Reimplemented from ISF::ITruthIncident.

Definition at line 344 of file Geant4TruthIncident.cxx.

                                                                                  {
  G4Track* track=m_step->GetTrack();
  const G4DynamicParticle* dynPart = track->GetDynamicParticle();
  const double preAsgnDecPropTime = dynPart->GetPreAssignedDecayProperTime();
  const double trackPropTime = track->GetProperTime();
  bool parentIsQuasiStable = (nullptr!=(dynPart->GetPreAssignedDecayProducts()));
  const G4VProcess *process = m_step->GetPostStepPoint()->GetProcessDefinedStep();
  const int processType = process->GetProcessType();
  ISF::InteractionClass_t classification(ISF::STD_VTX);
  if(parentIsQuasiStable) {
    if(this->parentSurvivesIncident()) {
      classification = ISF::QS_SURV_VTX;
    }
    else if(processType==6 && (trackPropTime-preAsgnDecPropTime) < -std::numeric_limits<double>::epsilon()) {
      // Particle decayed before its expected pre-defined decay time
      classification = ISF::QS_DEST_VTX;
    }
    else {
      classification = ISF::QS_PREDEF_VTX;
    }
  }
  return classification;
  // end of code to determine interaction type
}

numberOfChildren()

unsigned short ISF::ITruthIncident::numberOfChildren ( ) const

inlineinherited

Return total number of child particles.

Definition at line 139 of file ITruthIncident.h.

                                                           {
    return m_numChildren;
  }

parentBarcode()

int iGeant4::Geant4TruthIncident::parentBarcode ( )

finaloverridevirtual

Return the barcode of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 116 of file Geant4TruthIncident.cxx.

                                              { // TODO Remove this method
  auto parent = parentParticle();
 
  return (parent) ? HepMC::barcode(parent) : HepMC::UNDEFINED_ID;
}

parentEkin()

double iGeant4::Geant4TruthIncident::parentEkin ( ) const

finaloverridevirtual

Return Ekin of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 108 of file Geant4TruthIncident.cxx.

                                                    {
  return (m_step->GetPreStepPoint()->GetKineticEnergy()/CLHEP::MeV);
}

parentP2()

double iGeant4::Geant4TruthIncident::parentP2 ( ) const

finaloverridevirtual

Return p^2 of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 98 of file Geant4TruthIncident.cxx.

                                                  {
  const G4ThreeVector & mom= m_step->GetPreStepPoint()->GetMomentum();
  return mom.mag2();
}

parentParticle()

HepMC::GenParticlePtr iGeant4::Geant4TruthIncident::parentParticle ( )

finaloverridevirtual

Return the parent particle as a HepMC particle type.

Implements ISF::ITruthIncident.

Definition at line 134 of file Geant4TruthIncident.cxx.

                                                               {
  return m_atlasG4EvtUserInfo->GetCurrentGenParticle();
}

parentParticleAfterIncident()

HepMC::GenParticlePtr iGeant4::Geant4TruthIncident::parentParticleAfterIncident ( int  newBC )

finaloverridevirtual

Return the parent particle after the TruthIncident vertex (and give it a new barcode)

Implements ISF::ITruthIncident.

Definition at line 149 of file Geant4TruthIncident.cxx.

                                                                                          {
  const G4Track *track = m_step->GetTrack();
 
  // check if particle is a alive in G4 or in ISF
  if ( !parentSurvivesIncident() ) {
    return nullptr;
  }
 
  // internally cache the HepMC representation of the parent particle
  // (to allow multiple calls to this method on the same instance)
  if ( !m_parentParticleAfterIncident ) {
    // create new HepMC particle, using momentum and energy
    // from G4DynamicParticle (which should be equivalent to postStep)
    m_parentParticleAfterIncident = convert(track, newBarcode, false);
    
    m_atlasG4EvtUserInfo->SetCurrentGenParticle( m_parentParticleAfterIncident );
    
    // store (new) hepmc particle in track's UserInformation
    TrackHelper       tHelper(track);
    TrackInformation *tInfo = tHelper.GetTrackInformation();
    if (tInfo) {
      // do NOT update the TrackInformation for regenerated particles!
      // (most recent truth info is kept in AtlasG4EventUserInfo)
      //tInfo->SetCurrentGenParticle( m_parentParticleAfterIncident );
      int regenerationNr = tInfo->GetRegenerationNr();
      regenerationNr++;
      tInfo->SetRegenerationNr(regenerationNr);
      if ( tHelper.IsPrimary() ) { 
        tInfo->SetClassification(TrackInformation::RegeneratedPrimary);
      }
    }
 
  } // <-- end if m_parentParticleAfterIncident is not filled
 
 
  // return the HepMC particle representation of the parent
  return m_parentParticleAfterIncident;
}

parentPdgCode()

int iGeant4::Geant4TruthIncident::parentPdgCode ( ) const

finaloverridevirtual

Return the PDG Code of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 112 of file Geant4TruthIncident.cxx.

                                                    {
  return  m_step->GetTrack()->GetDefinition()->GetPDGEncoding();
}

parentPt2()

double iGeant4::Geant4TruthIncident::parentPt2 ( ) const

finaloverridevirtual

Return pT^2 of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 103 of file Geant4TruthIncident.cxx.

                                                   {
  const G4ThreeVector & mom= m_step->GetPreStepPoint()->GetMomentum();
  return mom.perp2();
}

parentStatus()

int iGeant4::Geant4TruthIncident::parentStatus ( )

finaloverridevirtual

Return the status of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 128 of file Geant4TruthIncident.cxx.

                                              {
 auto parent = parentParticle();
 
 return (parent) ? parent->status() : 0;
}

parentSurvivesIncident()

bool iGeant4::Geant4TruthIncident::parentSurvivesIncident ( ) const

finaloverridevirtual

Return a boolean whether or not the parent particle survives the incident.

Implements ISF::ITruthIncident.

Definition at line 138 of file Geant4TruthIncident.cxx.

                                                              { 
  const G4Track *track = m_step->GetTrack();
 
  // check if particle is a alive in G4 or in ISF
  if ( particleAlive(track) ) {
    return true;
  } else {
    return false;
  }
}

parentUniqueID()

int iGeant4::Geant4TruthIncident::parentUniqueID ( )

finaloverridevirtual

Return the unique ID of the parent particle.

Implements ISF::ITruthIncident.

Definition at line 122 of file Geant4TruthIncident.cxx.

                                               {
  auto parent = parentParticle();
 
  return (parent) ? HepMC::uniqueID(parent) : HepMC::UNDEFINED_ID;
}

particleAlive()

bool iGeant4::Geant4TruthIncident::particleAlive ( const G4Track *  track ) const

inlineprivate

check if the given G4Track represents a particle that is alive in ISF or ISF-G4

Definition at line 255 of file Geant4TruthIncident.cxx.

                                                                         {
  G4TrackStatus  trackStatus = track->GetTrackStatus();
 
  if ( trackStatus != fAlive && trackStatus != fStopButAlive ) {
    // parent does not exist in G4 anymore after this step
 
    // check whether the particle was returned to ISF
    auto*    trackInfo = ISFG4Helper::getISFTrackInfo( *track );
    bool returnedToISF = trackInfo ? trackInfo->GetReturnedToISF() : false;
    if ( !returnedToISF ) {
      // particle was not sent to ISF either
      return false;
    }
  }
 
  return true;
}

physicsProcessCategory()

int iGeant4::Geant4TruthIncident::physicsProcessCategory ( ) const

finaloverridevirtual

Return category of the physics process represented by the truth incident (eg hadronic, em, ..)

Implements ISF::ITruthIncident.

Definition at line 86 of file Geant4TruthIncident.cxx.

                                                             {
  const G4VProcess *process = m_step->GetPostStepPoint()->GetProcessDefinedStep();
  // TODO: need to check that G4ProcessSubTypes Match int
  return process->GetProcessType();
}

physicsProcessCode()

int iGeant4::Geant4TruthIncident::physicsProcessCode ( ) const

finaloverridevirtual

Return specific physics process code of the truth incident (eg ionisation, bremsstrahlung, ..)

Implements ISF::ITruthIncident.

Definition at line 92 of file Geant4TruthIncident.cxx.

                                                         {
  const G4VProcess *process = m_step->GetPostStepPoint()->GetProcessDefinedStep();
  // TODO: need to check that G4ProcessSubTypes Match int
  return process->GetProcessSubType();
}

position()

const HepMC::FourVector & iGeant4::Geant4TruthIncident::position ( ) const

finaloverridevirtual

Return HepMC position of the truth vertex.

Implements ISF::ITruthIncident.

Definition at line 82 of file Geant4TruthIncident.cxx.

                                                                  {
  return m_position;
}

prepareChildren()

void iGeant4::Geant4TruthIncident::prepareChildren ( )

inlineprivate

prepare the child particles

Definition at line 327 of file Geant4TruthIncident.cxx.

                                                 {
 
  const std::vector<const G4Track*>* tracks = m_step->GetSecondaryInCurrentStep();
  const int iSize = tracks->size();
  const int iLast = iSize - numberOfChildren() - 1; //NB can be -1.
  for(int i=iSize-1;i>iLast;i--) {
    m_children.push_back((*tracks)[i]);
  }
}

setChildPassedFilters()

void ISF::ITruthIncident::setChildPassedFilters ( unsigned short  index )

inlineinherited

Record that a particular child passed a check.

Definition at line 181 of file ITruthIncident.h.

                                                                    {
    m_childPassedFilters[index] = true;
    return;
  }

setPassWholeVertices()

void ISF::ITruthIncident::setPassWholeVertices ( bool  passWholeVertex )

inlineinherited

Set whether this TruthIncident should pass the vertex as a whole or individual children.

Definition at line 190 of file ITruthIncident.h.

                                                                   {
    m_passWholeVertex=passWholeVertex;
  }

Member Data Documentation

ATLAS_THREAD_SAFE

AtlasG4EventUserInfo* m_atlasG4EvtUserInfo iGeant4::Geant4TruthIncident::ATLAS_THREAD_SAFE {}

private

Definition at line 125 of file Geant4TruthIncident.h.

m_baseISP

const ISF::ISFParticle& iGeant4::Geant4TruthIncident::m_baseISP

private

Definition at line 123 of file Geant4TruthIncident.h.

m_childPassedFilters

std::vector<bool> ISF::ITruthIncident::m_childPassedFilters

protectedinherited

Definition at line 131 of file ITruthIncident.h.

m_children

std::vector<const G4Track*> iGeant4::Geant4TruthIncident::m_children

private

Definition at line 126 of file Geant4TruthIncident.h.

m_geoID

AtlasDetDescr::AtlasRegion ISF::ITruthIncident::m_geoID

privateinherited

region that the TruthIncident is located in

Definition at line 127 of file ITruthIncident.h.

m_numChildren

int ISF::ITruthIncident::m_numChildren

protectedinherited

Definition at line 129 of file ITruthIncident.h.

m_parentParticleAfterIncident

HepMC::GenParticlePtr iGeant4::Geant4TruthIncident::m_parentParticleAfterIncident {}

private

Definition at line 128 of file Geant4TruthIncident.h.

m_passWholeVertex

bool ISF::ITruthIncident::m_passWholeVertex

protectedinherited

Definition at line 130 of file ITruthIncident.h.

m_position

HepMC::FourVector iGeant4::Geant4TruthIncident::m_position

private

Definition at line 121 of file Geant4TruthIncident.h.

m_positionSet

bool iGeant4::Geant4TruthIncident::m_positionSet

private

Definition at line 120 of file Geant4TruthIncident.h.

m_step

const G4Step* iGeant4::Geant4TruthIncident::m_step {}

private

Definition at line 122 of file Geant4TruthIncident.h.

---

The documentation for this class was generated from the following files:

• Geant4TruthIncident.h
• Geant4TruthIncident.cxx