#include <ActsFatrasSimTool.h>

Inheritance diagram for ISF::ActsFatrasSimTool:

Collaboration diagram for ISF::ActsFatrasSimTool:

Classes
struct	HitSurfaceSelector
	Simple struct to select surfaces where hits should be generated. More...

struct	SingleParticleSimulation
	Single particle simulation with fixed propagator, interactions, and decay. More...

Public Types
using	Generator = std::ranlux48

using	Navigator = Acts::Navigator

using	ChargedStepper = Acts::EigenStepper< Acts::EigenStepperDefaultExtension >

using	ChargedPropagator = Acts::Propagator< ChargedStepper, Navigator >

using	NeutralStepper = Acts::StraightLineStepper

using	NeutralPropagator = Acts::Propagator< NeutralStepper, Navigator >

using	ChargedSelector = ActsFatras::ChargedSelector

using	ChargedInteractions = ActsFatras::StandardChargedElectroMagneticInteractions

using	ChargedSimulation = SingleParticleSimulation< ChargedPropagator, ChargedInteractions, HitSurfaceSelector, ActsFatras::NoDecay >

using	NeutralSelector = ActsFatras::NeutralSelector

using	NeutralInteractions = ActsFatras::InteractionList< ActsFatras::PhotonConversion >

using	NeutralSimulation = SingleParticleSimulation< NeutralPropagator, NeutralInteractions, ActsFatras::NoSurface, ActsFatras::NoDecay >

using	Simulation = ActsFatras::Simulation< ChargedSelector, ChargedSimulation, NeutralSelector, NeutralSimulation >

Public Member Functions
	ActsFatrasSimTool (const std::string &type, const std::string &name, const IInterface *parent)

virtual	~ActsFatrasSimTool ()

virtual StatusCode	initialize () override

virtual StatusCode	simulate (const EventContext &ctx, ISFParticle &isp, ISFParticleContainer &, McEventCollection *) override

virtual StatusCode	simulateVector (const EventContext &ctx, const ISFParticleVector &particles, ISFParticleContainer &secondaries, McEventCollection mcEventCollection, McEventCollection shadowTruth=nullptr) override
	Simulation call for vectors of particles. More...

virtual StatusCode	setupEvent (const EventContext &) override
	Setup Event chain - in case of a begin-of event action is needed. More...

virtual StatusCode	releaseEvent (const EventContext &ctx) override
	Release Event chain - in case of an end-of event action is needed. More...

virtual ISF::SimulationFlavor	simFlavor () const override

virtual Acts::MagneticFieldContext	getMagneticFieldContext (const EventContext &) const

virtual StatusCode	sysInitialize () override
	Gaudi sysInitialize() methods. More...

virtual StatusCode	setupEventST () override
	Setup Event chain - in case of a begin-of event action is needed (called by ISimulationSvc) More...

virtual StatusCode	releaseEventST () override
	Release Event chain - in case of an end-of event action is needed (called by ISimulationSvc) More...

const ChronoEntity *	chronoStart (const IChronoSvc::ChronoTag &tag)
	wrapper call to start chrono with given tag More...

const ChronoEntity *	chronoStop (const IChronoSvc::ChronoTag &tag)
	wrapper call to stop chrono with given tag More...

Protected Attributes
ServiceHandle< IChronoStatSvc >	m_chrono {this, "ChronoStatService", "ChronoStatSvc"}
	The timing service for general usage. More...

Private Member Functions
template<class T >
StatusCode	retrieveTool (ToolHandle< T > &thandle)

int	getATLASProcessCode (ActsFatras::ProcessType actspt)

Private Attributes
SiHitCollection	m_pixelSiHits

SiHitCollection	m_sctSiHits

ServiceHandle< IAthRNGSvc >	m_rngSvc {this, "RNGService", "AthRNGSvc"}

ATHRNG::RNGWrapper *m_randomEngine	ATLAS_THREAD_SAFE {}

Gaudi::Property< std::string >	m_randomEngineName

PublicToolHandle< ActsTrk::ITrackingGeometryTool >	m_trackingGeometryTool

std::shared_ptr< const Acts::TrackingGeometry >	m_trackingGeometry

SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCacheCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}

std::shared_ptr< const Acts::Logger >	m_logger {nullptr}

PublicToolHandle< ISF::IParticleFilter >	m_particleFilter

ServiceHandle< ISF::ITruthSvc >	m_truthRecordSvc {this, "TruthRecordService", "ISF_TruthRecordSvc", ""}

ToolHandle< ActsFatrasWriteHandler >	m_ActsFatrasWriteHandler

Gaudi::Property< double >	m_interact_minPt

Gaudi::Property< bool >	m_meanEnergyLoss {this, "MeanEnergyLoss", true, "Toggle between mean and mode evaluation of energy loss"}

Gaudi::Property< bool >	m_includeGgradient {this, "IncludeGgradient", true, "Boolean flag for inclusion of d(dEds)d(q/p) into energy loss"}

Gaudi::Property< double >	m_momentumCutOff {this, "MomentumCutOff", 0., "Cut-off value for the momentum in SI units"}

Gaudi::Property< double >	m_maxStep

Gaudi::Property< double >	m_maxRungeKuttaStepTrials

Gaudi::Property< double >	m_maxStepSize

Gaudi::Property< double >	m_pathLimit

Gaudi::Property< bool >	m_loopProtection

Gaudi::Property< double >	m_loopFraction

Gaudi::Property< double >	m_tolerance

Gaudi::Property< double >	m_stepSizeCutOff

Gaudi::Property< std::map< int, int > >	m_processTypeMap

Detailed Description

Definition at line 70 of file ActsFatrasSimTool.h.

Member Typedef Documentation

◆ ChargedInteractions

using ISF::ActsFatrasSimTool::ChargedInteractions = ActsFatras::StandardChargedElectroMagneticInteractions

Definition at line 188 of file ActsFatrasSimTool.h.

◆ ChargedPropagator

using ISF::ActsFatrasSimTool::ChargedPropagator = Acts::Propagator<ChargedStepper, Navigator>

Definition at line 179 of file ActsFatrasSimTool.h.

◆ ChargedSelector

using ISF::ActsFatrasSimTool::ChargedSelector = ActsFatras::ChargedSelector

Definition at line 187 of file ActsFatrasSimTool.h.

◆ ChargedSimulation

using ISF::ActsFatrasSimTool::ChargedSimulation = SingleParticleSimulation< ChargedPropagator, ChargedInteractions, HitSurfaceSelector, ActsFatras::NoDecay>

Definition at line 190 of file ActsFatrasSimTool.h.

◆ ChargedStepper

using ISF::ActsFatrasSimTool::ChargedStepper = Acts::EigenStepper<Acts::EigenStepperDefaultExtension>

Definition at line 178 of file ActsFatrasSimTool.h.

◆ Generator

using ISF::ActsFatrasSimTool::Generator = std::ranlux48

Definition at line 174 of file ActsFatrasSimTool.h.

◆ Navigator

using ISF::ActsFatrasSimTool::Navigator = Acts::Navigator

Definition at line 176 of file ActsFatrasSimTool.h.

◆ NeutralInteractions

using ISF::ActsFatrasSimTool::NeutralInteractions = ActsFatras::InteractionList<ActsFatras::PhotonConversion>

Definition at line 195 of file ActsFatrasSimTool.h.

◆ NeutralPropagator

using ISF::ActsFatrasSimTool::NeutralPropagator = Acts::Propagator<NeutralStepper, Navigator>

Definition at line 182 of file ActsFatrasSimTool.h.

◆ NeutralSelector

using ISF::ActsFatrasSimTool::NeutralSelector = ActsFatras::NeutralSelector

Definition at line 194 of file ActsFatrasSimTool.h.

◆ NeutralSimulation

using ISF::ActsFatrasSimTool::NeutralSimulation = SingleParticleSimulation< NeutralPropagator, NeutralInteractions, ActsFatras::NoSurface, ActsFatras::NoDecay>

Definition at line 196 of file ActsFatrasSimTool.h.

◆ NeutralStepper

using ISF::ActsFatrasSimTool::NeutralStepper = Acts::StraightLineStepper

Definition at line 181 of file ActsFatrasSimTool.h.

◆ Simulation

using ISF::ActsFatrasSimTool::Simulation = ActsFatras::Simulation<ChargedSelector, ChargedSimulation, NeutralSelector, NeutralSimulation>

Definition at line 200 of file ActsFatrasSimTool.h.

Constructor & Destructor Documentation

◆ ActsFatrasSimTool()

ISF::ActsFatrasSimTool::ActsFatrasSimTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 18 of file ActsFatrasSimTool.cxx.

21 : BaseSimulatorTool(type, name, parent) {}

◆ ~ActsFatrasSimTool()

ISF::ActsFatrasSimTool::~ActsFatrasSimTool ( )

virtual

Definition at line 23 of file ActsFatrasSimTool.cxx.

23 {}

Member Function Documentation

◆ chronoStart()

const ChronoEntity* ISF::BaseSimulatorTool::chronoStart ( const IChronoSvc::ChronoTag & tag )

inlineinherited

wrapper call to start chrono with given tag

Definition at line 98 of file BaseSimulatorTool.h.

                                                                      {
       if (m_chrono) return m_chrono->chronoStart( tag);
       return nullptr;
     }

◆ chronoStop()

const ChronoEntity* ISF::BaseSimulatorTool::chronoStop ( const IChronoSvc::ChronoTag & tag )

inlineinherited

wrapper call to stop chrono with given tag

Definition at line 104 of file BaseSimulatorTool.h.

                                                                     {
       if (m_chrono) return m_chrono->chronoStop( tag);
       return nullptr;
     }

◆ getATLASProcessCode()

int ISF::ActsFatrasSimTool::getATLASProcessCode ( ActsFatras::ProcessType actspt )

inlineprivate

Definition at line 304 of file ActsFatrasSimTool.h.

304 {return m_processTypeMap[static_cast<uint32_t>(actspt)];};

◆ getMagneticFieldContext()

Acts::MagneticFieldContext ISF::ActsFatrasSimTool::getMagneticFieldContext ( const EventContext & ctx ) const

virtual

Definition at line 232 of file ActsFatrasSimTool.cxx.

                                                                                                     {
   SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCacheCondObjInputKey, ctx};
   if (!readHandle.isValid()) {
     ATH_MSG_ERROR(name() + ": Failed to retrieve magnetic field condition data " + m_fieldCacheCondObjInputKey.key() + ".");
   }
   else ATH_MSG_DEBUG(name() << "retrieved magnetic field condition data "<< m_fieldCacheCondObjInputKey.key());
   const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
  
   return Acts::MagneticFieldContext(fieldCondObj);
 }

◆ initialize()

StatusCode ISF::ActsFatrasSimTool::initialize ( )

overridevirtual

Reimplemented from ISF::BaseSimulatorTool.

Definition at line 25 of file ActsFatrasSimTool.cxx.

                                             {
   ATH_CHECK(BaseSimulatorTool::initialize());
   ATH_MSG_INFO("ISF::ActsFatrasSimTool update with ACTS version: v"
     << Acts::VersionMajor << "." << Acts::VersionMinor << "."
     << Acts::VersionPatch << " [" << Acts::CommitHash.value_or("unknown hash") << "]");
   // Retrieve particle filter
   if (!m_particleFilter.empty()) ATH_CHECK(m_particleFilter.retrieve());
  
   // setup logger
   m_logger = makeActsAthenaLogger(this, std::string("ActsFatras"),std::string("ActsFatrasSimTool"));
  
   // retrive tracking geo tool
   ATH_CHECK(m_trackingGeometryTool.retrieve());
   m_trackingGeometry = m_trackingGeometryTool->trackingGeometry();
   
   //retrive Magnetfield tool
   ATH_MSG_VERBOSE("Using ATLAS magnetic field service");
   ATH_CHECK( m_fieldCacheCondObjInputKey.initialize());
  
   // Random number service
   if (m_rngSvc.retrieve().isFailure()) {
     ATH_MSG_FATAL("Could not retrieve " << m_rngSvc);
     return StatusCode::FAILURE;
   }
   // Get own engine with own seeds
   m_randomEngine = m_rngSvc->getEngine(this, m_randomEngineName.value());
   if (!m_randomEngine) {
     ATH_MSG_FATAL("Could not get random engine '" << m_randomEngineName.value() << "'");
     return StatusCode::FAILURE;
   }
  
   // ISF truth service
   ATH_CHECK (m_truthRecordSvc.retrieve());
   ATH_MSG_DEBUG( "- Using ISF TruthRecordSvc : " << m_truthRecordSvc.typeAndName() );
   return StatusCode::SUCCESS;
 }

◆ releaseEvent()

virtual StatusCode ISF::ActsFatrasSimTool::releaseEvent ( const EventContext & )

inlineoverridevirtual

Release Event chain - in case of an end-of event action is needed.

Reimplemented from ISF::BaseSimulatorTool.

Definition at line 222 of file ActsFatrasSimTool.h.

                                                                     {
     std::vector<SiHitCollection> hitcolls;
     hitcolls.push_back(m_pixelSiHits);
     hitcolls.push_back(m_sctSiHits);
     ATH_CHECK(m_ActsFatrasWriteHandler->WriteHits(hitcolls,ctx));
     ATH_CHECK(m_truthRecordSvc->releaseEvent());
     return StatusCode::SUCCESS; };

◆ releaseEventST()

virtual StatusCode ISF::BaseSimulatorTool::releaseEventST ( )

inlineoverridevirtualinherited

Release Event chain - in case of an end-of event action is needed (called by ISimulationSvc)

Reimplemented in ISF::FastCaloSimV2Tool, and ISF::BaseSimulatorG4Tool.

Definition at line 94 of file BaseSimulatorTool.h.

95 { return releaseEvent(Gaudi::Hive::currentContext()); }

◆ retrieveTool()

template<class T >

StatusCode ISF::ActsFatrasSimTool::retrieveTool ( ToolHandle< T > & thandle )

inlineprivate

Definition at line 241 of file ActsFatrasSimTool.h.

                                                   {
     if (!thandle.empty() && thandle.retrieve().isFailure()) {
       ATH_MSG_FATAL("Cannot retrieve " << thandle << ". Abort.");
       return StatusCode::FAILURE;
     } else ATH_MSG_DEBUG("Successfully retrieved " << thandle);
     return StatusCode::SUCCESS;
   }

◆ setupEvent()

virtual StatusCode ISF::ActsFatrasSimTool::setupEvent ( const EventContext & )

inlineoverridevirtual

Setup Event chain - in case of a begin-of event action is needed.

Reimplemented from ISF::BaseSimulatorTool.

Definition at line 217 of file ActsFatrasSimTool.h.

                                                               {
     ATH_CHECK(m_truthRecordSvc->initializeTruthCollection());
     m_pixelSiHits.Clear();
     m_sctSiHits.Clear();
     return StatusCode::SUCCESS; };

◆ setupEventST()

virtual StatusCode ISF::BaseSimulatorTool::setupEventST ( )

inlineoverridevirtualinherited

Setup Event chain - in case of a begin-of event action is needed (called by ISimulationSvc)

Reimplemented in ISF::FastCaloSimV2Tool, and ISF::BaseSimulatorG4Tool.

Definition at line 67 of file BaseSimulatorTool.h.

68 { return setupEvent(Gaudi::Hive::currentContext()); }

◆ simFlavor()

virtual ISF::SimulationFlavor ISF::ActsFatrasSimTool::simFlavor ( ) const

inlineoverridevirtual

Definition at line 229 of file ActsFatrasSimTool.h.

229 {

230 return ISF::Fatras; };

◆ simulate()

StatusCode ISF::ActsFatrasSimTool::simulate	(	const EventContext &	ctx,
		ISFParticle &	isp,
		ISFParticleContainer &	secondaries,
		McEventCollection *	mcEventCollection
	)

overridevirtual

Reimplemented from ISF::BaseSimulatorTool.

Definition at line 62 of file ActsFatrasSimTool.cxx.

                                         {
   ATH_MSG_VERBOSE("Particle " << isp << " received for simulation.");
   // Check if particle passes filter, if there is one
   if (!m_particleFilter.empty() && !m_particleFilter->passFilter(isp)) {
     ATH_MSG_VERBOSE("ISFParticle " << isp << " does not pass selection. Ignoring.");
     return StatusCode::SUCCESS;
   }
   // Process ParticleState from particle stack
   // Wrap the input ISFParticle in an STL vector with size of 1
   const ISF::ISFParticleVector ispVector(1, &isp);
   ATH_CHECK(this->simulateVector(ctx, ispVector, secondaries, mcEventCollection));
   ATH_MSG_VERBOSE("Simulation done");
   return StatusCode::SUCCESS;
 }

◆ simulateVector()

StatusCode ISF::ActsFatrasSimTool::simulateVector	(	const EventContext &	ctx,
		const ISFParticleVector &	particles,
		ISFParticleContainer &	secondaries,
		McEventCollection *	mcEventCollection,
		McEventCollection *	= `nullptr`
	)

overridevirtual

Simulation call for vectors of particles.

Reimplemented from ISF::BaseSimulatorTool.

Definition at line 79 of file ActsFatrasSimTool.cxx.

                                                                                              {
  
   m_randomEngine->setSeed(m_randomEngineName, ctx);
   CLHEP::HepRandomEngine* randomEngine = m_randomEngine->getEngine(ctx);
   Generator generator(CLHEP::RandFlat::shoot(randomEngine->flat()));
   ATH_MSG_VERBOSE(name() << " RNG seed " << CLHEP::RandFlat::shoot(randomEngine->flat()));
   ATH_MSG_VERBOSE(name() << " received vector of size "
                << particles.size() << " particles for simulation.");
  
   // construct the ACTS simulator
   Acts::Navigator navigator( Acts::Navigator::Config{ m_trackingGeometry }, m_logger);
   auto bField = std::make_shared<ATLASMagneticFieldWrapper>();
   auto chargedStepper = ChargedStepper(std::move(bField));
   auto neutralStepper = NeutralStepper();
   auto chargedPropagator = ChargedPropagator(chargedStepper, navigator, m_logger);
   auto neutralPropagator = NeutralPropagator(neutralStepper, navigator, m_logger);
   ChargedSimulation simulatorCharged(std::move(chargedPropagator), m_logger);
   NeutralSimulation simulatorNeutral(std::move(neutralPropagator), m_logger);
   Simulation simulator=Simulation(std::move(simulatorCharged),std::move(simulatorNeutral));
   ATH_MSG_VERBOSE(name() << " Min pT for interaction " << m_interact_minPt * Acts::UnitConstants::MeV << " GeV");
   // Acts propagater options
   simulator.charged.maxStepSize = m_maxStepSize;
   simulator.charged.maxStep = m_maxStep;
   simulator.charged.pathLimit = m_pathLimit;
   simulator.charged.maxRungeKuttaStepTrials = m_maxRungeKuttaStepTrials;
   simulator.charged.loopProtection = m_loopProtection;
   simulator.charged.loopFraction = m_loopFraction;
   simulator.charged.targetTolerance = m_tolerance;
   simulator.charged.stepSizeCutOff = m_stepSizeCutOff;
   // Create interaction list
   simulator.charged.interactions = ActsFatras::makeStandardChargedElectroMagneticInteractions(m_interact_minPt * Acts::UnitConstants::MeV);
   // get Geo and Mag map
   ATH_MSG_VERBOSE(name() << " Getting per event Geo and Mag map");
   Acts::MagneticFieldContext mctx = getMagneticFieldContext(ctx);
   const ActsTrk::GeometryContext& gctx = m_trackingGeometryTool->getNominalGeometryContext();
   auto anygctx = gctx.context();
   // Loop over ISFParticleVector and process each separately
   ATH_MSG_VERBOSE(name() << " Processing particles in ISFParticleVector.");
   for (const auto isfp : particles) {
     // ====ACTSFatras Simulation====
     // //  
     // input/output particle and hits containers
     // Convert to ActsFatras::Particle
     // ISF: Energy, mass, and momentum are in MeV, position in mm
     // Acts: Energy, mass, and momentum are in GeV, position in mm
     ATH_MSG_DEBUG(name() << " Convert ISF::Particle(mass) " << isfp->id()<<"|" << isfp<<"(" << isfp->mass() << ")");
     std::vector<ActsFatras::Particle> input = std::vector<ActsFatras::Particle>{
       ActsFatras::Particle(ActsFatras::Barcode().withVertexPrimary(0).withParticle(isfp->id()), static_cast<Acts::PdgParticle>(isfp->pdgCode()),
                            isfp->charge(),isfp->mass() * Acts::UnitConstants::MeV)
         .setDirection(Acts::makeDirectionFromPhiEta(isfp->momentum().phi(), isfp->momentum().eta()))
         .setAbsoluteMomentum(isfp->momentum().mag() * Acts::UnitConstants::MeV)
         .setPosition4(ActsTrk::convertPosToActs(isfp->position(), isfp->timeStamp()))};
     ATH_MSG_DEBUG(name() << " Propagating ActsFatras::Particle  vertex|particle|generation|subparticle, " << input[0]);
     std::vector<ActsFatras::Particle> simulatedInitial;
     std::vector<ActsFatras::Particle> simulatedFinal;
     std::vector<ActsFatras::Hit> hits;
     // simulate
     auto result=simulator.simulate(anygctx, mctx, generator, input, simulatedInitial, simulatedFinal, hits);
     auto simulatedFailure=result.value();
     if (simulatedFailure.size()>0){
       for (const auto& simfail : simulatedFailure){
         auto errCode = Acts::make_error_code(Acts::PropagatorError(simfail.error.value()));
         ATH_MSG_WARNING(name() << " Particle id " <<simfail.particle.particleId()<< ": fail to be simulated during Propagation: " << errCode.message());
         ATH_MSG_WARNING(name() << " Particle vertex|particle|generation|subparticle"<<simfail.particle << " starts from position" << Acts::toString(simfail.particle.position()) << " and direction " << Acts::toString(simfail.particle.direction()));
         return StatusCode::SUCCESS;
       }
     }
  
     ATH_MSG_DEBUG(name() << " initial particle " << simulatedInitial[0]);
     ATH_MSG_DEBUG(name() << " ActsFatras simulator hits: " << hits.size());
     int i = 0;
     for (const auto& hit : hits) {
       ATH_MSG_DEBUG(name() << " hit pos: " << hit.position() );
       ++i;
       if (i>5) break;
     }
     ATH_MSG_DEBUG(name() << " No. of particles after ActsFatras simulator: " << simulatedFinal.size());
     if (!simulatedFinal.empty()){
       ATH_MSG_DEBUG(name() << " start procesing secondaries");
       auto itr = simulatedFinal.begin();
       // Save hits of isfp
       std::vector<ActsFatras::Hit> particle_hits;
       if (itr->numberOfHits() > 0) {
         std::copy(hits.begin(), hits.begin()+itr->numberOfHits(), std::back_inserter(particle_hits));
         m_ActsFatrasWriteHandler->createHits(*isfp, m_trackingGeometry,particle_hits,m_pixelSiHits,m_sctSiHits);
       }
       // Process secondaries
       auto isKilled = !itr->isAlive();
       int maxGeneration = simulatedFinal.back().particleId().generation();
       ATH_MSG_DEBUG(name() << " maxGeneration: "<< maxGeneration);
       for (int gen = 0; gen <= maxGeneration; ++gen){
         ATH_MSG_DEBUG(name() << " start with generation "<< gen << "|" << maxGeneration << ": "<< *itr);
         auto vecsecisfp = std::make_unique<ISF::ISFParticleVector>();
         while (itr != simulatedFinal.end() && static_cast<int>(itr->particleId().generation()) == gen) {
           ATH_MSG_DEBUG(name() << " genration "<< gen << "|" << maxGeneration << ": "<< *itr);
           if(itr->isSecondary()){
             // convert final particles to ISF::particle
             const auto pos = ActsTrk::convertPosFromActs(itr->fourPosition()).first;
             const auto mom = ActsTrk::convertMomFromActs(itr->fourMomentum()).first;
             double mass = itr->mass() / Acts::UnitConstants::MeV;
             double charge = itr->charge();
             int pdgid = itr->pdg();
             auto properTime = ActsTrk::timeToAthena(itr->time());
             const int status = 1 + HepMC::SIM_STATUS_THRESHOLD;
             const int id = HepMC::UNDEFINED_ID;
             auto secisfp = new ISF::ISFParticle (pos,mom,mass,charge,pdgid,status,properTime,*isfp,id);
             ATH_MSG_DEBUG(name() <<" secondaries particle (ACTS): "<<*itr<< "("<<itr->momentum()<<")|time "<<itr->time()<<"|process "<< getATLASProcessCode(itr->process()));
             ATH_MSG_DEBUG(name() <<" secondaries particle (ISF): " << *secisfp << " time "<<secisfp->timeStamp());
             vecsecisfp->push_back(secisfp);
           }
           else{
             ATH_MSG_DEBUG(name() <<" primary particle found with generation "<< gen);
           }
           ++itr;
         }
         if (!vecsecisfp->empty()) {
           ISF::ISFTruthIncident truth(*isfp,
                                       *vecsecisfp,
                                       getATLASProcessCode((itr-1)->process()),
                                       isfp->nextGeoID(),
                                       isKilled&&gen==maxGeneration ? ISF::fKillsPrimary : ISF::fPrimarySurvives
                                     );
           ATH_MSG_DEBUG(name() << " Truth incident parentPt2(MinPt2) " << truth.parentPt2() <<" (100 MeV)");
           ATH_MSG_DEBUG(name() << " Truth incident ChildPt2(MinPt2) " << truth.childrenPt2Pass(300) <<" (300 MeV)");
           m_truthRecordSvc->registerTruthIncident(truth,  true);
           truth.updateParentAfterIncidentProperties();
           truth.updateChildParticleProperties();
           for (auto *secisfp : *vecsecisfp){
             if (secisfp->getTruthBinding()) {
                 secondaries.push_back(secisfp);
             }
             else {
                 ATH_MSG_WARNING("Secondary particle not written out to truth.\n Parent (" << isfp << ")\n Secondary (" << *secisfp <<")");
             }
           } // end of truth binding 
         }// end of store truth bind secondaries
       } 
     }// end of secondaries
     ATH_MSG_VERBOSE(name() << " No. of secondaries: " << secondaries.size());
     ATH_MSG_DEBUG(name() << " End of particle " << isfp->id());
  
     std::vector<ActsFatras::Particle>().swap(input);
     std::vector<ActsFatras::Particle>().swap(simulatedInitial);
     std::vector<ActsFatras::Particle>().swap(simulatedFinal);
     std::vector<ActsFatras::Hit>().swap(hits);
   } // end of isfp loop
   return StatusCode::SUCCESS;
 }

◆ sysInitialize()

virtual StatusCode ISF::BaseSimulatorTool::sysInitialize ( )

inlineoverridevirtualinherited

Gaudi sysInitialize() methods.

Definition at line 49 of file BaseSimulatorTool.h.

     {
       ATH_CHECK( AthAlgTool::sysInitialize() );
       ATH_CHECK( m_chrono.retrieve() );
       return StatusCode::SUCCESS;
     }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

ATHRNG::RNGWrapper* m_randomEngine ISF::ActsFatrasSimTool::ATLAS_THREAD_SAFE {}

private

Definition at line 251 of file ActsFatrasSimTool.h.

◆ m_ActsFatrasWriteHandler

ToolHandle<ActsFatrasWriteHandler> ISF::ActsFatrasSimTool::m_ActsFatrasWriteHandler

private

Initial value:

{

this, "ActsFatrasWriteHandler", "ActsFatrasWriteHandler"}

Definition at line 273 of file ActsFatrasSimTool.h.

◆ m_chrono

ServiceHandle<IChronoStatSvc> ISF::BaseSimulatorTool::m_chrono {this, "ChronoStatService", "ChronoStatSvc"}

protectedinherited

The timing service for general usage.

Definition at line 115 of file BaseSimulatorTool.h.

◆ m_fieldCacheCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> ISF::ActsFatrasSimTool::m_fieldCacheCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}

private

Definition at line 261 of file ActsFatrasSimTool.h.

◆ m_includeGgradient

Gaudi::Property<bool> ISF::ActsFatrasSimTool::m_includeGgradient {this, "IncludeGgradient", true, "Boolean flag for inclusion of d(dEds)d(q/p) into energy loss"}

private

Definition at line 281 of file ActsFatrasSimTool.h.

◆ m_interact_minPt

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_interact_minPt

private

Initial value:

{this, "Interact_MinPt", 50.0,

"Min pT of the interactions (MeV)"}

Definition at line 276 of file ActsFatrasSimTool.h.

◆ m_logger

std::shared_ptr<const Acts::Logger> ISF::ActsFatrasSimTool::m_logger {nullptr}

private

Definition at line 264 of file ActsFatrasSimTool.h.

◆ m_loopFraction

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_loopFraction

private

Initial value:

{this, "LoopFraction", 0.5,

"Allowed loop fraction, 1 is a full loop"}

Definition at line 294 of file ActsFatrasSimTool.h.

◆ m_loopProtection

Gaudi::Property<bool> ISF::ActsFatrasSimTool::m_loopProtection

private

Initial value:

{this, "LoopProtection", true,

"Loop protection, it adapts the pathLimit"}

Definition at line 292 of file ActsFatrasSimTool.h.

◆ m_maxRungeKuttaStepTrials

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_maxRungeKuttaStepTrials

private

Initial value:

{this, "MaxRungeKuttaStepTrials", 10000,

"Maximum number of Runge-Kutta steps for the stepper step call"}

Definition at line 286 of file ActsFatrasSimTool.h.

◆ m_maxStep

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_maxStep

private

Initial value:

{this, "MaxSteps", 1000,

"Max number of steps"}

Definition at line 284 of file ActsFatrasSimTool.h.

◆ m_maxStepSize

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_maxStepSize

private

Initial value:

{this, "MaxStepSize", 3.0,

"Max step size (converted to Acts::UnitConstants::m)"}

Definition at line 288 of file ActsFatrasSimTool.h.

◆ m_meanEnergyLoss

Gaudi::Property<bool> ISF::ActsFatrasSimTool::m_meanEnergyLoss {this, "MeanEnergyLoss", true, "Toggle between mean and mode evaluation of energy loss"}

private

Definition at line 280 of file ActsFatrasSimTool.h.

◆ m_momentumCutOff

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_momentumCutOff {this, "MomentumCutOff", 0., "Cut-off value for the momentum in SI units"}

private

Definition at line 282 of file ActsFatrasSimTool.h.

◆ m_particleFilter

PublicToolHandle<ISF::IParticleFilter> ISF::ActsFatrasSimTool::m_particleFilter

private

Initial value:

{

this, "ParticleFilter", "", "Particle filter kinematic cuts, etc."}

Definition at line 267 of file ActsFatrasSimTool.h.

◆ m_pathLimit

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_pathLimit

private

Initial value:

{this, "PathLimit", 100.0,

"Track path limit (converted to Acts::UnitConstants::cm)"}

Definition at line 290 of file ActsFatrasSimTool.h.

◆ m_pixelSiHits

SiHitCollection ISF::ActsFatrasSimTool::m_pixelSiHits

private

Definition at line 237 of file ActsFatrasSimTool.h.

◆ m_processTypeMap

Gaudi::Property<std::map<int,int> > ISF::ActsFatrasSimTool::m_processTypeMap

private

Initial value:

{this, "ProcessTypeMap",

{{0,0}, {1,201}, {2,14}, {3,3}, {4,121}}, "proessType map <ActsFatras,G4>"}

Definition at line 301 of file ActsFatrasSimTool.h.

◆ m_randomEngineName

Gaudi::Property<std::string> ISF::ActsFatrasSimTool::m_randomEngineName

private

Initial value:

{this, "RandomEngineName",

"RandomEngineName", "Name of random number stream"}

Definition at line 252 of file ActsFatrasSimTool.h.

◆ m_rngSvc

ServiceHandle<IAthRNGSvc> ISF::ActsFatrasSimTool::m_rngSvc {this, "RNGService", "AthRNGSvc"}

private

Definition at line 250 of file ActsFatrasSimTool.h.

◆ m_sctSiHits

SiHitCollection ISF::ActsFatrasSimTool::m_sctSiHits

private

Definition at line 238 of file ActsFatrasSimTool.h.

◆ m_stepSizeCutOff

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_stepSizeCutOff

private

Initial value:

{this, "StepSizeCutOff", 0.,

"Cut-off value for the step size"}

Definition at line 298 of file ActsFatrasSimTool.h.

◆ m_tolerance

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_tolerance

private

Initial value:

{this, "Tolerance", 0.0001,

"Tolerance for the error of the integration"}

Definition at line 296 of file ActsFatrasSimTool.h.

◆ m_trackingGeometry

std::shared_ptr<const Acts::TrackingGeometry> ISF::ActsFatrasSimTool::m_trackingGeometry

private

Definition at line 258 of file ActsFatrasSimTool.h.

◆ m_trackingGeometryTool

PublicToolHandle<ActsTrk::ITrackingGeometryTool> ISF::ActsFatrasSimTool::m_trackingGeometryTool

private

Initial value:

{

this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}

Definition at line 256 of file ActsFatrasSimTool.h.

◆ m_truthRecordSvc

ServiceHandle<ISF::ITruthSvc> ISF::ActsFatrasSimTool::m_truthRecordSvc {this, "TruthRecordService", "ISF_TruthRecordSvc", ""}

private

Definition at line 270 of file ActsFatrasSimTool.h.

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

Classes

Public Types

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ ChargedInteractions

◆ ChargedPropagator

◆ ChargedSelector

◆ ChargedSimulation

◆ ChargedStepper

◆ Generator

◆ Navigator

◆ NeutralInteractions

◆ NeutralPropagator

◆ NeutralSelector

◆ NeutralSimulation

◆ NeutralStepper

◆ Simulation

Constructor & Destructor Documentation

◆ ActsFatrasSimTool()

◆ ~ActsFatrasSimTool()

Member Function Documentation

◆ chronoStart()

◆ chronoStop()

◆ getATLASProcessCode()

◆ getMagneticFieldContext()

◆ initialize()

◆ releaseEvent()

◆ releaseEventST()

◆ retrieveTool()

◆ setupEvent()

◆ setupEventST()

◆ simFlavor()

◆ simulate()

◆ simulateVector()

◆ sysInitialize()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_ActsFatrasWriteHandler

◆ m_chrono

◆ m_fieldCacheCondObjInputKey

◆ m_includeGgradient

◆ m_interact_minPt

◆ m_logger

◆ m_loopFraction

◆ m_loopProtection

◆ m_maxRungeKuttaStepTrials

◆ m_maxStep

◆ m_maxStepSize

◆ m_meanEnergyLoss

◆ m_momentumCutOff

◆ m_particleFilter

◆ m_pathLimit

◆ m_pixelSiHits

◆ m_processTypeMap

◆ m_randomEngineName

◆ m_rngSvc

◆ m_sctSiHits

◆ m_stepSizeCutOff

◆ m_tolerance

◆ m_trackingGeometry

◆ m_trackingGeometryTool

◆ m_truthRecordSvc