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
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	updateCurrentBaseISFParticle (ISF::ISFParticle *baseISFParticle)
	Set the base ISFParticle for 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 43 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 46 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = nullptr;
  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 212 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 52 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = 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 76 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;
        auto* trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(*aTrack);
        HepMC::GenParticlePtr currentGenParticle = trackInfo ? trackInfo->GetCurrentGenParticle() : nullptr;
        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
      if (secondaryVector && !secondaryVector->empty()) {
        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
      } // <-- end if secondary vector has entries
 
      return;
    }

newISFParticle()

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

private

Definition at line 244 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 224 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->GetGenerationZeroGenParticle();
 
      ISF::TruthBinding* tBinding = new ISF::TruthBinding(currentGenParticle, primaryGenParticle, generationZeroGenParticle);
 
      return tBinding;
    }

PostUserTrackingAction()

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

finaloverridevirtualinherited

Definition at line 236 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 102 of file TrackProcessorUserActionBase.cxx.

{
  // AdePT currently does not return the G4Secondary vector per G4Step. As the vector is required to attach
  // the G4VUserInfo below the G4VUserInfo is not available, which would lead to nullptr access crashes
#ifdef ATHSIMULATION_USE_ADEPT
  return;
#endif
 
  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 263 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 242 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 121 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();
 
  // Keep the legacy classification based on generation-zero identity; the
  // true primary is stored separately on TrackInformation by ISFG4Helper.
  HepMC::GenParticlePtr generationZeroGenParticle = truthBinding->getGenerationZeroGenParticle();
  HepMC::GenParticlePtr primaryGenParticle = generationZeroGenParticle;
  HepMC::GenParticlePtr currentGenParticle = truthBinding->getCurrentGenParticle();
 
  auto classification = classify(primaryGenParticle,
                                 generationZeroGenParticle,
                                 currentGenParticle,
                                 regenerationNr);
 
  auto* newTrackInfo = ::iGeant4::ISFG4Helper::attachTrackInfoToNewG4Track(aTrack,
                                                                 *baseISP,
                                                                 classification,
                                                                 generationZeroGenParticle );
  newTrackInfo->SetRegenerationNr(regenerationNr);
 
  updateCurrentBaseISFParticle(baseISP);
 
  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 193 of file TrackProcessorUserActionBase.cxx.

{
  auto* trackInfo = ::iGeant4::ISFG4Helper::getISFTrackInfo(aTrack);
 
  ISF::ISFParticle* baseISFParticle = trackInfo->GetBaseISFParticle();
 
  updateCurrentBaseISFParticle(baseISFParticle);
 
  return;
}

updateCurrentBaseISFParticle()

void G4UA::iGeant4::TrackProcessorUserActionBase::updateCurrentBaseISFParticle ( ISF::ISFParticle * baseISFParticle )

privateinherited

Set the base ISFParticle for the currently traced particle.

Definition at line 204 of file TrackProcessorUserActionBase.cxx.

{
  m_curBaseISP = baseISFParticle;
  return;
}

UserSteppingAction()

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

finaloverridevirtualinherited

Definition at line 58 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 );
 
  // AdePT currently does not return the G4Secondary vector per G4Step. As the vector is required to attach
  // the G4VUserInfo in the PreUserTrackingAction below, the G4VUserInfo is not available, which
  // would lead to nullptr access crashes
#ifdef ATHSIMULATION_USE_ADEPT
  return;
#endif
  //
  // 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();
  if (secondaryVector && !secondaryVector->empty()) {
    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_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 68 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 45 of file TrackProcessorUserActionBase.h.

---

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

• TrackProcessorUserActionPassBack.h
• TrackProcessorUserActionPassBack.cxx