G4UA::iGeant4::TrackProcessorUserActionPassBack Class Reference

#include <TrackProcessorUserActionPassBack.h>

Inheritance diagram for G4UA::iGeant4::TrackProcessorUserActionPassBack:

Collaboration diagram for G4UA::iGeant4::TrackProcessorUserActionPassBack:

Classes
struct	Config

Public Member Functions
	TrackProcessorUserActionPassBack (const Config &config)
virtual void	BeginOfEventAction (const G4Event *) override final
virtual void	EndOfEventAction (const G4Event *) override final
virtual void	PreUserTrackingAction (const G4Track *) override
virtual void	PostUserTrackingAction (const G4Track *) override final
virtual void	UserSteppingAction (const G4Step *) override final
ISF::ISFParticleContainer	ReturnSecondaries (ISF::ISFParticle const *parent)

Protected Attributes
AtlasG4EventUserInfo *	m_atlasG4EvtUserInfo {}
	event-global G4 UserInformation
ISF::ISFParticleContainer	m_storedSecondaries

Private Member Functions
void	ISFSteppingAction (const G4Step *, ISF::ISFParticle *curISP) override final
	Called by the base class after the G4Track->ISFParticle association has been established.
ISF::TruthBinding *	newTruthBinding (const G4Track *aTrack, HepMC::GenParticlePtr currentGenParticle) const
	create a new TruthBinding object for the given G4Track (may return 0 if unable)
ISF::ISFParticle *	newISFParticle (G4Track *aTrack, const ISF::ISFParticle *parent, HepMC::GenParticlePtr currentGenParticle, AtlasDetDescr::AtlasRegion nextGeoID)
void	returnParticleToISF (G4Track *aTrack, const ISF::ISFParticle *parentISP, HepMC::GenParticlePtr currentGenParticle, AtlasDetDescr::AtlasRegion nextGeoID)
	kills the given G4Track, converts it into an ISFParticle and returns it to the ISF particle broker
void	setupPrimary (G4Track &)
	Setup the given G4Track as the current primary particle which we'll process.
void	setupSecondary (const G4Track &)
	Setup the given G4Track as the current secondary particle which we'll process.
void	updateCachedParticleInfo (ISF::ISFParticle *baseISFParticle, HepMC::ConstGenParticlePtr primaryGenParticle, HepMC::GenParticlePtr currentGenParticle)
	Set the following information as the currently traced particle.
VTrackInformation::TrackClassification	classify (HepMC::ConstGenParticlePtr primaryGenParticle, HepMC::ConstGenParticlePtr generationZeroGenParticle, HepMC::ConstGenParticlePtr currentGenParticle, int regenerationNumber) const
	Classify the particle represented by the given set of truth links.

Private Attributes
Config	m_config
ISF::IParticleBroker *	m_particleBrokerQuick {}
	quick access avoiding gaudi overhead
ISF::IGeoIDSvc *	m_geoIDSvcQuick {}
	quick access avoiding gaudi overhead
ISF::ISFParticle *	m_curBaseISP {}
	The most recent ISFParticle ancestor that triggers the currently processed G4Track.

Detailed Description

Definition at line 34 of file TrackProcessorUserActionPassBack.h.

Constructor & Destructor Documentation

TrackProcessorUserActionPassBack()

G4UA::iGeant4::TrackProcessorUserActionPassBack::TrackProcessorUserActionPassBack

(

const Config &

config

)

Definition at line 44 of file TrackProcessorUserActionPassBack.cxx.

      : m_config(config)
    {
      if(4<m_config.verboseLevel)
        {
          G4cout << "Initializing TrackProcessorUserActionPassBack" << G4endl;
        }
 
      if ( !m_config.particleBroker.empty() ) {
        if (m_config.particleBroker.retrieve().isFailure()) {
          G4ExceptionDescription description;
          description << G4String("TrackProcessorUserActionPassBack: ") + "Could not retrieve ISF Particle Broker: " << m_config.particleBroker;
          G4Exception("iGeant4::TrackProcessorUserActionPassBack", "NoISFParticleBroker", FatalException, description);
          return; //The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
        }
        m_particleBrokerQuick = &(*m_config.particleBroker);
      }
 
      if ( !m_config.geoIDSvc.empty() ) {
        if (m_config.geoIDSvc.retrieve().isFailure()) {
          G4ExceptionDescription description;
          description << G4String("TrackProcessorUserActionPassBack: ") + "Could not retrieve ISF GeoID Svc: " << m_config.geoIDSvc;
          G4Exception("iGeant4::TrackProcessorUserActionPassBack", "NoISFGeoIDSvc", FatalException, description);
          return; //The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
        }
 
        m_geoIDSvcQuick = &(*m_config.geoIDSvc);
      }
 
      return;
 
    }

Member Function Documentation

BeginOfEventAction()

void G4UA::iGeant4::TrackProcessorUserActionBase::BeginOfEventAction ( const G4Event * )

finaloverridevirtualinherited

Definition at line 48 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = nullptr;
  m_atlasG4EvtUserInfo = ::iGeant4::ISFG4Helper::getAtlasG4EventUserInfo();
  return;
}

classify()

VTrackInformation::TrackClassification G4UA::iGeant4::TrackProcessorUserActionBase::classify ( HepMC::ConstGenParticlePtr primaryGenParticle, HepMC::ConstGenParticlePtr generationZeroGenParticle, HepMC::ConstGenParticlePtr currentGenParticle, int regenerationNumber ) const

privateinherited

Classify the particle represented by the given set of truth links.

Definition at line 208 of file TrackProcessorUserActionBase.cxx.

{
  // if particle points to a non-zero truth particle it can not just be a 'simple' Secondary
 
  if (currentGenParticle) {
    if (currentGenParticle == primaryGenParticle) {
      return VTrackInformation::Primary;
    }
    else if (generationZeroGenParticle == primaryGenParticle && regenerationNumber>0) {
      return VTrackInformation::RegeneratedPrimary;
    }
    else {
      return VTrackInformation::RegisteredSecondary;
    }
  }
 
  return VTrackInformation::Secondary;
}

EndOfEventAction()

void G4UA::iGeant4::TrackProcessorUserActionBase::EndOfEventAction ( const G4Event * )

finaloverridevirtualinherited

Definition at line 55 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = nullptr;
  m_atlasG4EvtUserInfo = nullptr;
  return;
}

ISFSteppingAction()

void G4UA::iGeant4::TrackProcessorUserActionPassBack::ISFSteppingAction ( const G4Step * aStep, ISF::ISFParticle * curISP )

finaloverrideprivatevirtual

Called by the base class after the G4Track->ISFParticle association has been established.

Implements G4UA::iGeant4::TrackProcessorUserActionBase.

Definition at line 77 of file TrackProcessorUserActionPassBack.cxx.

    {
      G4Track*           aTrack   = aStep->GetTrack();
      G4TrackStatus aTrackStatus  = aTrack->GetTrackStatus();
 
      //const G4StepPoint *preStep  = aStep->GetPreStepPoint(); // Only used for DEBUG messages
      const G4StepPoint *postStep = aStep->GetPostStepPoint();
 
      // get geoID from parent
      AtlasDetDescr::AtlasRegion curGeoID = curISP->nextGeoID();
 
      //std::cout<<"retrieved isp "<<curISP<<" for trackID "<<curISP<<std::endl;
 
      // check geoID of postStep
      const G4ThreeVector     &postPos  = postStep->GetPosition();
      //const G4ThreeVector     &postMom  = postStep->GetMomentum();
      //AtlasDetDescr::AtlasRegion  nextGeoID = m_geoIDSvcQuick->identifyNextGeoID( postPos.x(),
      //                                                                            postPos.y(),
      //                                                                            postPos.z(),
      //                                                                            postMom.x(),
      //                                                                            postMom.y(),
      //                                                                            postMom.z() );
      AtlasDetDescr::AtlasRegion  nextGeoID = m_geoIDSvcQuick->identifyGeoID( postPos.x(),
                                                                              postPos.y(),
                                                                              postPos.z() );
 
      //ATH_MSG_DEBUG("PostStep point resolved to geoID = "<<nextGeoID);
 
      // return if particle did NOT cross boundary
      if ( nextGeoID==curGeoID ) {
        //ATH_MSG_DEBUG(" -> G4Track stays inside geoID = "<<curGeoID);
 
        // //
        // // for debugging:
        // //
        // if ( msgLvl(MSG::DEBUG) ) {
        //   const G4ThreeVector &prePos = preStep->GetPosition();
        //   AtlasDetDescr::AtlasRegion  preStepGeoID  = m_geoIDSvcQuick->identifyGeoID( prePos.x(),
        //                                                                               prePos.y(),
        //                                                                               prePos.z() );
        //   AtlasDetDescr::AtlasRegion  postStepGeoID = m_geoIDSvcQuick->identifyGeoID( postPos.x(),
        //                                                                               postPos.y(),
        //                                                                               postPos.z() );
 
        //   if( preStepGeoID!=postStepGeoID ) {
        //     const G4VPhysicalVolume     *preVol  = preStep->GetPhysicalVolume();
        //     const G4VPhysicalVolume     *postVol = postStep->GetPhysicalVolume();
        //     const G4ThreeVector         &preMom  = preStep->GetMomentum();
        //     const G4ThreeVector        &postMom  = postStep->GetMomentum();
        //     const G4TrackVector *secondaryVector = aStep->GetSecondary();
        //     const G4ThreeVector&      aTrack_pos = aTrack->GetPosition();
        //     const G4ThreeVector&      aTrack_mom = aTrack->GetMomentum();
        //     int pdgID=aTrack->GetDefinition()->GetPDGEncoding();
        //     int                         aTrackID = aTrack->GetTrackID();
        //     ATH_MSG_WARNING("pre "<<preVol->GetName()<<" x="<<prePos.x()<<" y="<<prePos.y()<<" z="<<prePos.z()<<" p="<<preMom.mag()<<" geoID="<<preStepGeoID<<"; post "<<postVol->GetName()<<" x="<<postPos.x()<<" y="<<postPos.y()<<" z="<<postPos.z()<<" p="<<postMom.mag()<<" geoID="<<nextGeoID<<"; length="<<aStep->GetStepLength()<<"; n2nd="<<secondaryVector->size()<<" track  x="<<aTrack_pos.x()<<" y="<<aTrack_pos.y()<<" z="<<aTrack_pos.z()<<" p="<<aTrack_mom.mag()<<" curgeoID="<<curGeoID<<" pdgid="<<pdgID<<" trackID="<<aTrackID<<" ISF="<<curISP<<"; ploss="<<(postMom-preMom).mag());
        //   }
        // }
 
        return;
      }
 
 
      //
      // this point is only reached if particle has crossed
      // a sub-det boundary in the non-Geant4-only mode
      //
 
      if ( aTrack->GetKineticEnergy() < m_config.passBackEkinThreshold ) {
        // kinetic energy of primary particle below threshold
        // ATH_MSG_DEBUG(" -> G4Track enters geoID = " << nextGeoID <<
        //               " but is below Ekin threshold. Not returned to ISF.");
        if ( m_config.killBoundaryParticlesBelowThreshold ) {
          aTrack->SetTrackStatus( fStopAndKill );
        } else {
          // TODO: link G4Track to ISF particle with the new GeoID
        }
      } else if ( aTrackStatus!=fAlive && aTrackStatus != fStopButAlive  ) {
        // particle is killed by G4 in this step
        // TODO: do we need to handle this case specifically?
        // ATH_MSG_DEBUG(" -> G4Track enters geoID = " << nextGeoID <<
        //               " but is destroyed in this step. Not returned to ISF.");
 
      } else {
        // particle is above kinetic energy threshold and alive after this step
        // -> push new ISFParticle back to ISF particle broker
        // ATH_MSG_DEBUG(" -> G4Track enters geoID = " << nextGeoID <<
        //               " and is returned to ISF.");
 
        const ISF::ISFParticle*    parent = curISP;
        HepMC::GenParticlePtr currentGenParticle = m_atlasG4EvtUserInfo->GetCurrentGenParticle();
        this->returnParticleToISF(aTrack, parent, currentGenParticle, nextGeoID); // TODO CHECK THIS LOGIC
      }
 
      //
      // handle secondaries that were created in this G4Step
      //
      const std::vector<const G4Track*> *secondaryVector = aStep->GetSecondaryInCurrentStep();
      // loop over new secondaries
      for ( auto* aConstTrack_2nd : *secondaryVector ) {
        // get a non-const G4Track for current secondary (nasty!)
        G4Track *aTrack_2nd ATLAS_THREAD_SAFE = const_cast<G4Track*>( aConstTrack_2nd ); // imposed by Geant4 interface
 
        // check if new secondary position is behind boundary
        const G4ThreeVector&             pos_2nd = aTrack_2nd->GetPosition();
        AtlasDetDescr::AtlasRegion nextGeoID_2nd = m_geoIDSvcQuick->identifyGeoID( pos_2nd.x(),
                                                                                   pos_2nd.y(),
                                                                                   pos_2nd.z() );
        if( nextGeoID_2nd!=curGeoID ) {
          // secondary was generated in this step and has
          // a different geoID than the currently tracked one
 
          if ( aTrack_2nd->GetKineticEnergy() < m_config.passBackEkinThreshold ) {
            // kinetic energy of secondary particle below threshold
            // ATH_MSG_DEBUG(" -> Secondary particle generated in this G4Step does not pass Ekin cut." <<
            //               " Not returned to ISF.");
            if ( m_config.killBoundaryParticlesBelowThreshold ) {
              // TODO: should we use fKillTrackAndSecondaries instead?
              aTrack_2nd->SetTrackStatus( fStopAndKill );
            } else {
              // TODO: link G4Track to ISF particle with the new GeoID
            }
          } else {
            // secondary particle is above kinetic energy threshold
            // -> return it to ISF
            // ATH_MSG_DEBUG(" -> Secondary particle generated in this G4Step is returned to ISF.");
 
            // attach TrackInformation instance to the new secondary G4Track
            ISF::ISFParticle *parent                  = curISP;
            HepMC::GenParticlePtr generationZeroGenParticle = nullptr;
            ::iGeant4::ISFG4Helper::attachTrackInfoToNewG4Track( *aTrack_2nd,
                                                       *parent,
                                                       VTrackInformation::Secondary,
                                                       generationZeroGenParticle );
 
            HepMC::GenParticlePtr currentGenParticle{};
            returnParticleToISF(aTrack_2nd, parent, currentGenParticle, nextGeoID_2nd); // TODO CHECK THIS LOGIC
          }
        }
 
      } // <-- end loop over new secondaries
 
      return;
    }

newISFParticle()

ISF::ISFParticle * G4UA::iGeant4::TrackProcessorUserActionPassBack::newISFParticle ( G4Track * aTrack, const ISF::ISFParticle * parent, HepMC::GenParticlePtr currentGenParticle, AtlasDetDescr::AtlasRegion nextGeoID )

private

Definition at line 242 of file TrackProcessorUserActionPassBack.cxx.

    {
      ISF::TruthBinding* tBinding = newTruthBinding(aTrack, currentGenParticle);
 
      ISF::ISFParticle* isp = ::iGeant4::ISFG4Helper::convertG4TrackToISFParticle( *aTrack,
                                                                         *parentISP,
                                                                         tBinding );
 
      if (nextGeoID!=AtlasDetDescr::fUndefinedAtlasRegion) {
        isp->setNextGeoID( AtlasDetDescr::AtlasRegion(nextGeoID) );
        isp->setNextSimID( parentISP->nextSimID() );
      }
 
      return isp;
    }

newTruthBinding()

ISF::TruthBinding * G4UA::iGeant4::TrackProcessorUserActionPassBack::newTruthBinding ( const G4Track * aTrack, HepMC::GenParticlePtr currentGenParticle ) const

private

create a new TruthBinding object for the given G4Track (may return 0 if unable)

Definition at line 222 of file TrackProcessorUserActionPassBack.cxx.

    {
      auto* trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(*aTrack);
      if (!trackInfo) {
        G4ExceptionDescription description;
        description << G4String("newTruthBinding: ") + "No TrackInformation associated with G4Track (trackID: "
                    << aTrack->GetTrackID() << ", track pos: "<<aTrack->GetPosition() << ", mom: "<<aTrack->GetMomentum()
                    << ", parentID " << aTrack->GetParentID() << ")";
        G4Exception("iGeant4::TrackProcessorUserActionPassBack", "NoTrackInformation", FatalException, description);
        return nullptr; //The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
      }
 
      HepMC::GenParticlePtr         primaryGenParticle = trackInfo->GetPrimaryGenParticle();
      HepMC::GenParticlePtr  generationZeroGenParticle = trackInfo->GetCurrentGenParticle(); // TODO CHECK THIS LOGIC
 
      ISF::TruthBinding* tBinding = new ISF::TruthBinding(currentGenParticle, primaryGenParticle, generationZeroGenParticle);
 
      return tBinding;
    }

PostUserTrackingAction()

void G4UA::iGeant4::TrackProcessorUserActionBase::PostUserTrackingAction ( const G4Track * )

finaloverridevirtualinherited

Definition at line 232 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = nullptr;
  return;
}

PreUserTrackingAction()

void G4UA::iGeant4::TrackProcessorUserActionBase::PreUserTrackingAction ( const G4Track * aTrack )

overridevirtualinherited

Reimplemented in G4UA::iGeant4::TrackProcessorUserActionFullG4.

Definition at line 98 of file TrackProcessorUserActionBase.cxx.

{
  bool isPrimary = ! aTrack->GetParentID();
  if (isPrimary) {
    G4Track* nonConstTrack ATLAS_THREAD_SAFE = const_cast<G4Track*> (aTrack); // imposed by Geant4 interface
    setupPrimary(*nonConstTrack);
  } else {
    setupSecondary(*aTrack);
  }
 
  return;
}

returnParticleToISF()

void G4UA::iGeant4::TrackProcessorUserActionPassBack::returnParticleToISF ( G4Track * aTrack, const ISF::ISFParticle * parentISP, HepMC::GenParticlePtr currentGenParticle, AtlasDetDescr::AtlasRegion nextGeoID )

private

kills the given G4Track, converts it into an ISFParticle and returns it to the ISF particle broker

Definition at line 261 of file TrackProcessorUserActionPassBack.cxx.

    {
      // kill track inside G4
      aTrack->SetTrackStatus( fStopAndKill );
 
      // create new ISFParticle and attach it to current G4Track
      ISF::ISFParticle *newISP = newISFParticle( aTrack, parentISP, currentGenParticle, nextGeoID );
 
      // update TrackInformation
      auto trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(*aTrack);
      trackInfo->SetReturnedToISF( true );
      trackInfo->SetBaseISFParticle( newISP );
 
      // push the particle back to ISF via the particle broker
      if ( m_particleBrokerQuick )
      {
        [&] ATLAS_NOT_THREAD_SAFE () {  // suppress checker warning, in MT mode there is no broker
          m_particleBrokerQuick->push(newISP, parentISP);
        }();
      }
      else
      {
        // store the particle for retrieval in MT mode
        //parentISP must be non-null by here, it has already been deeferenced
        if (!newISP->getTruthBinding()) newISP->setTruthBinding(new ISF::TruthBinding(*parentISP->getTruthBinding()));
        m_storedSecondaries.push_back( newISP );
      }
 
      return;
    }

ReturnSecondaries()

ISF::ISFParticleContainer G4UA::iGeant4::TrackProcessorUserActionBase::ReturnSecondaries ( ISF::ISFParticle const * parent )

inherited

Definition at line 238 of file TrackProcessorUserActionBase.cxx.

{
  // For now, just return all particles
  ISF::ISFParticleContainer result;
  std::swap( result, m_storedSecondaries );
  return result;
}

setupPrimary()

void G4UA::iGeant4::TrackProcessorUserActionBase::setupPrimary ( G4Track & aTrack )

privateinherited

Setup the given G4Track as the current primary particle which we'll process.

Definition at line 111 of file TrackProcessorUserActionBase.cxx.

{
  //
  // Get PrimaryParticleInformation from G4PrimaryParticle (assigned by TransportTool::addPrimaryVertex)
  //
 
  auto* trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(aTrack);
  if ( trackInfo ) {
    G4ExceptionDescription description;
    description << G4String("PreUserTrackingAction: ")
                << "Started simulation of primary particle which already has a TrackInformation/TrackBarcodeInfo object attached (trackID: "
                << aTrack.GetTrackID() << ", track pos: "<<aTrack.GetPosition() << ", mom: "<<aTrack.GetMomentum()
                << ", parentID " << aTrack.GetParentID() << ")";
    G4Exception("iGeant4::TrackProcessorUserActionBase", "TrackInformationAlreadyExists", FatalException, description);
    return; // The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
  }
 
  auto* primaryPartInfo = dynamic_cast <PrimaryParticleInformation*> (aTrack.GetDynamicParticle()->GetPrimaryParticle()->GetUserInformation());
  if (!primaryPartInfo) {
    G4ExceptionDescription description;
    description << G4String("PreUserTrackingAction: ") + "NULL PrimaryParticleInformation pointer for current G4Step (trackID "
                << aTrack.GetTrackID() << ", track pos: "<<aTrack.GetPosition() << ", mom: "<<aTrack.GetMomentum()
                << ", parentID " << aTrack.GetParentID() << ")";
    G4Exception("iGeant4::TrackProcessorUserActionBase", "NoPPInfo", FatalException, description);
    return; // The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
  }
 
  // get base ISFParticle and link to TrackInformation
  auto* baseISP = primaryPartInfo->GetISFParticle();
  if (!baseISP) {
    G4ExceptionDescription description;
    description << G4String("PreUserTrackingAction: ") + "No ISFParticle associated with primary particle (trackID: "
                << aTrack.GetTrackID() << ", track pos: "<<aTrack.GetPosition() << ", mom: "<<aTrack.GetMomentum()
                << ", parentID " << aTrack.GetParentID() << ")";
    G4Exception("iGeant4::TrackProcessorUserActionBase", "NoISFParticle", FatalException, description);
    return; // The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
  }
 
  ISF::TruthBinding* truthBinding = baseISP->getTruthBinding();
  if (!truthBinding) {
    G4ExceptionDescription description;
    description << G4String("PreUserTrackingAction: ") + "No ISF::TruthBinding associated with primary particle (trackID: "
                << aTrack.GetTrackID() << ", track pos: "<<aTrack.GetPosition() << ", mom: "<<aTrack.GetMomentum()
                << ", parentID " << aTrack.GetParentID() << ")";
    G4Exception("iGeant4::TrackProcessorUserActionBase", "NoISFTruthBinding", FatalException, description);
    return; // The G4Exception call above should abort the job, but Coverity does not seem to pick this up.
  }
 
  int regenerationNr = primaryPartInfo->GetRegenerationNr();
 
  HepMC::GenParticlePtr primaryGenParticle = truthBinding->getGenerationZeroGenParticle();
  HepMC::GenParticlePtr generationZeroGenParticle = truthBinding->getGenerationZeroGenParticle();
  HepMC::GenParticlePtr currentGenParticle = truthBinding->getCurrentGenParticle();
 
  auto classification = classify(primaryGenParticle,
                                 generationZeroGenParticle,
                                 currentGenParticle,
                                 regenerationNr);
 
  auto* newTrackInfo = ::iGeant4::ISFG4Helper::attachTrackInfoToNewG4Track(aTrack,
                                                                 *baseISP,
                                                                 classification,
                                                                 generationZeroGenParticle );
  newTrackInfo->SetRegenerationNr(regenerationNr);
 
  updateCachedParticleInfo(baseISP,
                     primaryGenParticle,
                     currentGenParticle);
 
  return;
}

setupSecondary()

void G4UA::iGeant4::TrackProcessorUserActionBase::setupSecondary ( const G4Track & aTrack )

privateinherited

Setup the given G4Track as the current secondary particle which we'll process.

Definition at line 183 of file TrackProcessorUserActionBase.cxx.

{
  auto* trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(aTrack);
 
  HepMC::GenParticlePtr currentGenParticle = trackInfo->GetCurrentGenParticle();
  HepMC::GenParticlePtr primaryGenParticle = trackInfo->GetPrimaryGenParticle();
  ISF::ISFParticle* baseISFParticle = trackInfo->GetBaseISFParticle();
 
  updateCachedParticleInfo(baseISFParticle, primaryGenParticle, currentGenParticle);
 
  return;
}

updateCachedParticleInfo()

void G4UA::iGeant4::TrackProcessorUserActionBase::updateCachedParticleInfo ( ISF::ISFParticle * baseISFParticle, HepMC::ConstGenParticlePtr primaryGenParticle, HepMC::GenParticlePtr currentGenParticle )

privateinherited

Set the following information as the currently traced particle.

Definition at line 196 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = baseISFParticle;
  m_atlasG4EvtUserInfo->SetCurrentPrimaryGenParticle( primaryGenParticle );
  m_atlasG4EvtUserInfo->SetCurrentGenParticle( currentGenParticle );
  return;
}

UserSteppingAction()

void G4UA::iGeant4::TrackProcessorUserActionBase::UserSteppingAction ( const G4Step * aStep )

finaloverridevirtualinherited

Definition at line 62 of file TrackProcessorUserActionBase.cxx.

{
  // get geoID from parent
  //TODO ELLI AtlasDetDescr::AtlasRegion curGeoID = m_curBaseISP->nextGeoID();
  //TODO ELLI ATH_MSG_DEBUG( "Currently simulating TrackID = " << aStep->GetTrack()->GetTrackID() <<
  //TODO ELLI                " inside geoID = " << curGeoID );
 
  //
  // call the ISFSteppingAction method of the implementation
  //
  ISFSteppingAction( aStep, m_curBaseISP );
 
  //
  // propagate the current ISFParticle link to all secondaries
  //
  const std::vector<const G4Track*>  *secondaryVector = aStep->GetSecondaryInCurrentStep();
  for ( auto* aConstSecondaryTrack : *secondaryVector ) {
    // get a non-const G4Track for current secondary (nasty!)
    G4Track* aSecondaryTrack ATLAS_THREAD_SAFE = const_cast<G4Track*>( aConstSecondaryTrack ); // imposed by Geant4 interface
 
    auto *trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(*aSecondaryTrack);
 
    // G4Tracks aready returned to ISF will have a TrackInformation attached to them
    bool particleReturnedToISF = trackInfo && trackInfo->GetReturnedToISF();
    if (!particleReturnedToISF) {
      HepMC::GenParticlePtr generationZeroGenParticle{};
      ::iGeant4::ISFG4Helper::attachTrackInfoToNewG4Track( *aSecondaryTrack,
                                                *m_curBaseISP,
                                                VTrackInformation::Secondary,
                                                generationZeroGenParticle );
    }
  } // <- loop over secondaries from this step
 
  return;
}

Member Data Documentation

m_atlasG4EvtUserInfo

AtlasG4EventUserInfo* G4UA::iGeant4::TrackProcessorUserActionBase::m_atlasG4EvtUserInfo {}

protectedinherited

event-global G4 UserInformation

Definition at line 48 of file TrackProcessorUserActionBase.h.

{};   

m_config

Config G4UA::iGeant4::TrackProcessorUserActionPassBack::m_config

private

Definition at line 55 of file TrackProcessorUserActionPassBack.h.

m_curBaseISP

ISF::ISFParticle* G4UA::iGeant4::TrackProcessorUserActionBase::m_curBaseISP {}

privateinherited

The most recent ISFParticle ancestor that triggers the currently processed G4Track.

Definition at line 75 of file TrackProcessorUserActionBase.h.

{};

m_geoIDSvcQuick

ISF::IGeoIDSvc* G4UA::iGeant4::TrackProcessorUserActionPassBack::m_geoIDSvcQuick {}

private

quick access avoiding gaudi overhead

Definition at line 72 of file TrackProcessorUserActionPassBack.h.

{}; 

m_particleBrokerQuick

ISF::IParticleBroker* G4UA::iGeant4::TrackProcessorUserActionPassBack::m_particleBrokerQuick {}

private

quick access avoiding gaudi overhead

Definition at line 71 of file TrackProcessorUserActionPassBack.h.

{}; 

m_storedSecondaries

ISF::ISFParticleContainer G4UA::iGeant4::TrackProcessorUserActionBase::m_storedSecondaries

protectedinherited

Definition at line 50 of file TrackProcessorUserActionBase.h.

---

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

• TrackProcessorUserActionPassBack.h
• TrackProcessorUserActionPassBack.cxx