ATLAS Offline Software
Loading...
Searching...
No Matches
ISF::ActsFatrasSimTool Class Reference

#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
using ChargedSimulation
using NeutralSelector = ActsFatras::NeutralSelector
using NeutralInteractions = ActsFatras::InteractionList<ActsFatras::PhotonConversion>
using NeutralSimulation
using Simulation

Public Member Functions

Amg::Vector3D convertMom3FromActs (const Acts::Vector3 &actsMom)
 Convert ACTS momentum to Athena momentum.
Amg::Vector3D convertPos3FromActs (const Acts::Vector3 &actsPos)
 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.
virtual StatusCode setupEvent (const EventContext &) override
 Setup Event chain - in case of a begin-of event action is needed.
virtual StatusCode releaseEvent (const EventContext &ctx) override
 Release Event chain - in case of an end-of event action is needed.
virtual ISF::SimulationFlavor simFlavor () const override
virtual Acts::MagneticFieldContext getMagneticFieldContext (const EventContext &) const
virtual StatusCode sysInitialize () override
 Gaudi sysInitialize() methods.
virtual StatusCode setupEventST () override
 Setup Event chain - in case of a begin-of event action is needed (called by ISimulationSvc)
virtual StatusCode releaseEventST () override
 Release Event chain - in case of an end-of event action is needed (called by ISimulationSvc)
const ChronoEntity * chronoStart (const IChronoSvc::ChronoTag &tag)
 wrapper call to start chrono with given tag
const ChronoEntity * chronoStop (const IChronoSvc::ChronoTag &tag)
 wrapper call to stop chrono with given tag

Protected Attributes

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

Private Member Functions

template<class T>
StatusCode retrieveTool (ToolHandle< T > &thandle)
bool checkStartSurface (const Acts::MagneticFieldContext &mctx, const Acts::GeometryContext &anygctx, const ChargedPropagator &chargedPropagator, const Acts::BoundTrackParameters &startParameters, Acts::Direction navDir=Acts::Direction::Forward(), double pathLimit=std::numeric_limits< double >::max()) const
int getATLASProcessCode (ActsFatras::ProcessType actspt)

Private Attributes

SiHitCollection m_pixelSiHits
SiHitCollection m_sctSiHits
ServiceHandle< IAthRNGSvcm_rngSvc {this, "RNGService", "AthRNGSvc"}
ATHRNG::RNGWrapper *m_randomEngine ATLAS_THREAD_SAFE {}
Gaudi::Property< std::string > m_randomEngineName
ServiceHandle< ISF::IGeoIDSvcm_geoIDSvc {this, "GeoIDSvc", "ISF::GeoIDSvc"}
PublicToolHandle< ActsTrk::IExtrapolationToolm_extrapolationTool {this, "ExtrapolationTool", "ActsExtrapolationTool"}
PublicToolHandle< ActsTrk::ITrackingGeometryToolm_trackingGeometryTool
std::shared_ptr< const Acts::TrackingGeometry > m_trackingGeometry
SG::ReadCondHandleKey< AtlasFieldCacheCondObjm_fieldCacheCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}
std::shared_ptr< const Acts::Logger > m_logger {nullptr}
PublicToolHandle< ISF::IParticleFilterm_particleFilter
ServiceHandle< ISF::ITruthSvcm_truthRecordSvc {this, "TruthRecordService", "ISF_TruthRecordSvc", ""}
ToolHandle< ActsFatrasWriteHandlerm_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 73 of file ActsFatrasSimTool.h.

Member Typedef Documentation

◆ ChargedInteractions

Initial value:
ActsFatras::StandardChargedElectroMagneticInteractions

Definition at line 191 of file ActsFatrasSimTool.h.

◆ ChargedPropagator

Definition at line 182 of file ActsFatrasSimTool.h.

◆ ChargedSelector

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

Definition at line 190 of file ActsFatrasSimTool.h.

◆ ChargedSimulation

Initial value:
ActsFatras::NoDecay>
ActsFatras::StandardChargedElectroMagneticInteractions ChargedInteractions
Acts::Propagator< ChargedStepper, Navigator > ChargedPropagator
Simple struct to select surfaces where hits should be generated.
Single particle simulation with fixed propagator, interactions, and decay.

Definition at line 193 of file ActsFatrasSimTool.h.

◆ ChargedStepper

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

Definition at line 181 of file ActsFatrasSimTool.h.

◆ Generator

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

Definition at line 177 of file ActsFatrasSimTool.h.

◆ Navigator

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

Definition at line 179 of file ActsFatrasSimTool.h.

◆ NeutralInteractions

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

Definition at line 198 of file ActsFatrasSimTool.h.

◆ NeutralPropagator

Definition at line 185 of file ActsFatrasSimTool.h.

◆ NeutralSelector

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

Definition at line 197 of file ActsFatrasSimTool.h.

◆ NeutralSimulation

Initial value:
NeutralPropagator, NeutralInteractions, ActsFatras::NoSurface,
ActsFatras::NoDecay>
Acts::Propagator< NeutralStepper, Navigator > NeutralPropagator
ActsFatras::InteractionList< ActsFatras::PhotonConversion > NeutralInteractions

Definition at line 199 of file ActsFatrasSimTool.h.

◆ NeutralStepper

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

Definition at line 184 of file ActsFatrasSimTool.h.

◆ Simulation

Initial value:
ActsFatras::Simulation<ChargedSelector, ChargedSimulation,
ActsFatras::NeutralSelector NeutralSelector
ActsFatras::ChargedSelector ChargedSelector
SingleParticleSimulation< NeutralPropagator, NeutralInteractions, ActsFatras::NoSurface, ActsFatras::NoDecay > NeutralSimulation
SingleParticleSimulation< ChargedPropagator, ChargedInteractions, HitSurfaceSelector, ActsFatras::NoDecay > ChargedSimulation

Definition at line 203 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 20 of file ActsFatrasSimTool.cxx.

23 : BaseSimulatorTool(type, name, parent) {}
BaseSimulatorTool()
Default constructor.

◆ ~ActsFatrasSimTool()

ISF::ActsFatrasSimTool::~ActsFatrasSimTool ( )
virtual

Definition at line 25 of file ActsFatrasSimTool.cxx.

25{}

Member Function Documentation

◆ checkStartSurface()

bool ISF::ActsFatrasSimTool::checkStartSurface ( const Acts::MagneticFieldContext & mctx,
const Acts::GeometryContext & anygctx,
const ChargedPropagator & chargedPropagator,
const Acts::BoundTrackParameters & startParameters,
Acts::Direction navDir = Acts::Direction::Forward(),
double pathLimit = std::numeric_limits<double>::max() ) const
private

Definition at line 413 of file ActsFatrasSimTool.cxx.

419{
420
421 ATH_MSG_VERBOSE(name() << "::" << __FUNCTION__ << " begin");
422 using ActorList =
423 Acts::ActorList<Acts::detail::SteppingLogger, Acts::MaterialInteractor>;
424 using PropagatorOptions = typename ChargedPropagator::template Options<ActorList>;
425
426 ATH_MSG_VERBOSE(name() << "::" << __FUNCTION__ << " Setting up propagator options for start surface check.");
427 PropagatorOptions options(anygctx, mctx);
428 options.loopProtection = (Acts::VectorHelpers::perp(startParameters.momentum()) < 300 * 1_MeV);
429 options.direction = navDir;
430 options.pathLimit = pathLimit;
431
432 // The state creation triggers the initial volume/surface lookup
433 ATH_MSG_VERBOSE(name() << "::" << __FUNCTION__ << " Initializing propagator state with start parameters: position "
434 << startParameters.position(anygctx).transpose() << ", momentum "
435 << startParameters.momentum().transpose());
436 auto state = chargedPropagator.makeState(options);
437 ATH_MSG_VERBOSE(name() << "::" << __FUNCTION__ << " Created propagator state. Now initializing with start parameters.");
438 auto initResult = chargedPropagator.initialize(state, startParameters);
439 if (!initResult.ok()) {
440 ATH_MSG_WARNING(name() << "::" << __FUNCTION__ << " Failed to initialize propagator state: "
441 << initResult.error().message());
442 return false;
443 }
444 return true;
445}
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)

◆ chronoStart()

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

wrapper call to start chrono with given tag

Definition at line 97 of file BaseSimulatorTool.h.

97 {
98 if (m_chrono) return m_chrono->chronoStart( tag);
99 return nullptr;
100 }
ServiceHandle< IChronoStatSvc > m_chrono
The timing service for general usage.

◆ chronoStop()

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

wrapper call to stop chrono with given tag

Definition at line 103 of file BaseSimulatorTool.h.

103 {
104 if (m_chrono) return m_chrono->chronoStop( tag);
105 return nullptr;
106 }

◆ convertMom3FromActs()

Amg::Vector3D ISF::ActsFatrasSimTool::convertMom3FromActs ( const Acts::Vector3 & actsMom)
inline

Convert ACTS momentum to Athena momentum.

Definition at line 207 of file ActsFatrasSimTool.h.

207 {
208 Amg::Vector3D threeMom{Amg::Vector3D::Zero()};
209 threeMom[Amg::x] = ActsTrk::energyToAthena(actsMom[Acts::eMom0]);
210 threeMom[Amg::y] = ActsTrk::energyToAthena(actsMom[Acts::eMom1]);
211 threeMom[Amg::z] = ActsTrk::energyToAthena(actsMom[Acts::eMom2]);
212 return threeMom;
213 }
constexpr double energyToAthena(const double actsE)
Converts an energy scalar from Acts to Athena units.
Eigen::Matrix< double, 3, 1 > Vector3D

◆ convertPos3FromActs()

Amg::Vector3D ISF::ActsFatrasSimTool::convertPos3FromActs ( const Acts::Vector3 & actsPos)
inline

Definition at line 215 of file ActsFatrasSimTool.h.

215 {
216 Amg::Vector3D pos{Amg::Vector3D::Zero()};
217 pos[Amg::x] = ActsTrk::lengthToAthena(actsPos[Acts::ePos0]);
218 pos[Amg::y] = ActsTrk::lengthToAthena(actsPos[Acts::ePos1]);
219 pos[Amg::z] = ActsTrk::lengthToAthena(actsPos[Acts::ePos2]);
220 return pos;
221 }
constexpr double lengthToAthena(const double actsL)
Converts a length scalar from Acts to Athena units.

◆ getATLASProcessCode()

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

Definition at line 337 of file ActsFatrasSimTool.h.

337{return m_processTypeMap[static_cast<uint32_t>(actspt)];};
Gaudi::Property< std::map< int, int > > m_processTypeMap
setEventNumber uint32_t

◆ getMagneticFieldContext()

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

Definition at line 402 of file ActsFatrasSimTool.cxx.

402 {
403 SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCacheCondObjInputKey, ctx};
404 if (!readHandle.isValid()) {
405 ATH_MSG_ERROR(name() + ": Failed to retrieve magnetic field condition data " + m_fieldCacheCondObjInputKey.key() + ".");
406 }
407 else ATH_MSG_DEBUG(name() << "retrieved magnetic field condition data "<< m_fieldCacheCondObjInputKey.key());
408 const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
409
410 return Acts::MagneticFieldContext(fieldCondObj);
411}
#define ATH_MSG_ERROR(x)
#define ATH_MSG_DEBUG(x)
SG::ReadCondHandleKey< AtlasFieldCacheCondObj > m_fieldCacheCondObjInputKey

◆ initialize()

StatusCode ISF::ActsFatrasSimTool::initialize ( )
overridevirtual

Reimplemented from ISF::BaseSimulatorTool.

Definition at line 27 of file ActsFatrasSimTool.cxx.

27 {
29 ATH_MSG_INFO("ISF::ActsFatrasSimTool update with ACTS version: v"
30 << Acts::VersionMajor << "." << Acts::VersionMinor << "."
31 << Acts::VersionPatch << " [" << Acts::CommitHash.value_or("unknown hash") << "]");
32 // Retrieve particle filter
33 ATH_CHECK(m_particleFilter.retrieve());
34 ATH_MSG_INFO("Using particle filter: " << m_particleFilter.typeAndName());
35 // setup logger
36 m_logger = makeActsAthenaLogger(this, std::string("ActsFatras"),std::string("ActsFatrasSimTool"));
37
38 // Geometry identifier service
39 if ( !m_geoIDSvc.empty() && m_geoIDSvc.retrieve().isFailure()){
40 ATH_MSG_FATAL ("Could not retrieve " << m_geoIDSvc);
41 return StatusCode::FAILURE;
42 }
43
44 // Acts Extrapolator
46 ATH_MSG_INFO( "- ActsExtrapolationTool : " << m_extrapolationTool.typeAndName() );
47
48 // retrive tracking geo tool
50 m_trackingGeometry = m_trackingGeometryTool->trackingGeometry();
51
52 //retrive Magnetfield tool
53 ATH_MSG_VERBOSE("Using ATLAS magnetic field service");
55
56 // Random number service
57 if (m_rngSvc.retrieve().isFailure()) {
58 ATH_MSG_FATAL("Could not retrieve " << m_rngSvc);
59 return StatusCode::FAILURE;
60 }
61 // Get own engine with own seeds
62 m_randomEngine = m_rngSvc->getEngine(this, m_randomEngineName.value());
63 if (!m_randomEngine) {
64 ATH_MSG_FATAL("Could not get random engine '" << m_randomEngineName.value() << "'");
65 return StatusCode::FAILURE;
66 }
67
68 // ISF truth service
69 ATH_CHECK (m_truthRecordSvc.retrieve());
70 ATH_MSG_DEBUG( "- Using ISF TruthRecordSvc : " << m_truthRecordSvc.typeAndName() );
71 return StatusCode::SUCCESS;
72}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_FATAL(x)
#define ATH_MSG_INFO(x)
std::unique_ptr< const Acts::Logger > makeActsAthenaLogger(IMessageSvc *svc, const std::string &name, int level, std::optional< std::string > parent_name)
PublicToolHandle< ActsTrk::ITrackingGeometryTool > m_trackingGeometryTool
ServiceHandle< IAthRNGSvc > m_rngSvc
PublicToolHandle< ActsTrk::IExtrapolationTool > m_extrapolationTool
std::shared_ptr< const Acts::TrackingGeometry > m_trackingGeometry
std::shared_ptr< const Acts::Logger > m_logger
ServiceHandle< ISF::ITruthSvc > m_truthRecordSvc
PublicToolHandle< ISF::IParticleFilter > m_particleFilter
ServiceHandle< ISF::IGeoIDSvc > m_geoIDSvc
Gaudi::Property< std::string > m_randomEngineName
virtual StatusCode initialize() override

◆ 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 241 of file ActsFatrasSimTool.h.

241 {
242 std::vector<SiHitCollection> hitcolls;
243 hitcolls.push_back(m_pixelSiHits);
244 hitcolls.push_back(m_sctSiHits);
245 ATH_CHECK(m_ActsFatrasWriteHandler->WriteHits(hitcolls,ctx));
246 ATH_CHECK(m_truthRecordSvc->releaseEvent());
247 return StatusCode::SUCCESS; };
SiHitCollection m_pixelSiHits
ToolHandle< ActsFatrasWriteHandler > m_ActsFatrasWriteHandler

◆ 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::BaseSimulatorG4Tool, and ISF::FastCaloSimV2Tool.

Definition at line 93 of file BaseSimulatorTool.h.

94 { return releaseEvent(Gaudi::Hive::currentContext()); }
virtual StatusCode releaseEvent(const EventContext &) override
Release Event chain - in case of an end-of event action is needed.

◆ retrieveTool()

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

Definition at line 260 of file ActsFatrasSimTool.h.

260 {
261 if (!thandle.empty() && thandle.retrieve().isFailure()) {
262 ATH_MSG_FATAL("Cannot retrieve " << thandle << ". Abort.");
263 return StatusCode::FAILURE;
264 } else ATH_MSG_DEBUG("Successfully retrieved " << thandle);
265 return StatusCode::SUCCESS;
266 }

◆ 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 236 of file ActsFatrasSimTool.h.

236 {
237 ATH_CHECK(m_truthRecordSvc->initializeTruthCollection());
238 m_pixelSiHits.Clear();
239 m_sctSiHits.Clear();
240 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::BaseSimulatorG4Tool, and ISF::FastCaloSimV2Tool.

Definition at line 66 of file BaseSimulatorTool.h.

67 { return setupEvent(Gaudi::Hive::currentContext()); }
virtual StatusCode setupEvent(const EventContext &) override
Setup Event chain - in case of a begin-of event action is needed.

◆ simFlavor()

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

Definition at line 248 of file ActsFatrasSimTool.h.

248 {
249 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 74 of file ActsFatrasSimTool.cxx.

76 {
77 ATH_MSG_VERBOSE("Particle " << isp << " received for simulation.");
78 // Check if particle passes filter, if there is one
79 if (!m_particleFilter.empty() && !m_particleFilter->passFilter(isp)) {
80 ATH_MSG_VERBOSE("ISFParticle " << isp << " does not pass selection. Ignoring.");
81 return StatusCode::SUCCESS;
82 }
83 // Process ParticleState from particle stack
84 // Wrap the input ISFParticle in an STL vector with size of 1
85 const ISF::ISFParticleVector ispVector(1, &isp);
86 ATH_CHECK(this->simulateVector(ctx, ispVector, secondaries, mcEventCollection));
87 ATH_MSG_VERBOSE("Simulation done");
88 return StatusCode::SUCCESS;
89}
virtual StatusCode simulateVector(const EventContext &ctx, const ISFParticleVector &particles, ISFParticleContainer &secondaries, McEventCollection *mcEventCollection, McEventCollection *shadowTruth=nullptr) override
Simulation call for vectors of particles.
std::vector< ISF::ISFParticle * > ISFParticleVector
ISFParticle vector.

◆ 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 91 of file ActsFatrasSimTool.cxx.

95 {
96 // filter particles
97 std::vector<ISFParticle*> selectedParticles;
98 for (const auto isfp : particles) {
99 if (!m_particleFilter.empty() && !m_particleFilter->passFilter(*isfp)) {
100 ATH_MSG_VERBOSE("ISFParticle " << *isfp << " does not pass selection. Ignoring.");
101 continue;
102 }
103 selectedParticles.push_back(isfp);
104 }
105 if (selectedParticles.empty()) {
106 ATH_MSG_VERBOSE("No particles passed selection. Ignoring.");
107 return StatusCode::SUCCESS;
108 }
109 // Set random seed for current event
110 m_randomEngine->setSeed(m_randomEngineName, ctx);
111 CLHEP::HepRandomEngine* randomEngine = m_randomEngine->getEngine(ctx);
112 Generator generator(CLHEP::RandFlat::shoot(randomEngine->flat()));
113 ATH_MSG_VERBOSE(name() << " RNG seed " << CLHEP::RandFlat::shoot(randomEngine->flat()));
114 ATH_MSG_VERBOSE(name() << " received vector of size "
115 << selectedParticles.size() << " particles for simulation.");
116
117 // construct the ACTS simulator
118 Acts::Navigator navigator( Acts::Navigator::Config{ m_trackingGeometry }, m_logger);
119 auto bField = std::make_shared<ATLASMagneticFieldWrapper>();
120 auto chargedStepper = ChargedStepper(std::move(bField));
121 auto neutralStepper = NeutralStepper();
122 auto chargedPropagator = ChargedPropagator(chargedStepper, navigator, m_logger);
123 auto neutralPropagator = NeutralPropagator(neutralStepper, navigator, m_logger);
124 ChargedSimulation simulatorCharged(std::move(chargedPropagator), m_logger);
125 NeutralSimulation simulatorNeutral(std::move(neutralPropagator), m_logger);
126 Simulation simulator=Simulation(std::move(simulatorCharged),std::move(simulatorNeutral));
127 ATH_MSG_VERBOSE(name() << " Min pT for interaction " << m_interact_minPt * Acts::UnitConstants::MeV << " GeV");
128 // Acts propagater options
129 simulator.charged.maxStepSize = m_maxStepSize;
130 simulator.charged.maxStep = m_maxStep;
131 simulator.charged.pathLimit = m_pathLimit;
132 simulator.charged.maxRungeKuttaStepTrials = m_maxRungeKuttaStepTrials;
133 simulator.charged.loopProtection = m_loopProtection;
134 simulator.charged.loopFraction = m_loopFraction;
135 simulator.charged.targetTolerance = m_tolerance;
136 simulator.charged.stepSizeCutOff = m_stepSizeCutOff;
137 // Create interaction list
138 simulator.charged.interactions = ActsFatras::makeStandardChargedElectroMagneticInteractions(m_interact_minPt * Acts::UnitConstants::MeV);
139
140 // Construct the ACTS propagator for starting surface check
141 auto surfaceCheckPropagator = ChargedPropagator(chargedStepper, navigator);
142
143 // get Geo and Mag map
144 ATH_MSG_VERBOSE(name() << " Getting per event Geo and Mag map");
145 Acts::MagneticFieldContext mctx = getMagneticFieldContext(ctx);
146 const auto& gctx = m_trackingGeometryTool->getGeometryContext(ctx);
147 auto anygctx = gctx.context();
148 // Loop over ISFParticleVector and process each separately
149 ATH_MSG_VERBOSE(name() << " Processing particles in ISFParticleVector.");
150 for (const auto isfp : selectedParticles) {
151 // ====ACTSFatras Simulation====
152 // //
153 // input/output particle and hits containers
154 // Convert to ActsFatras::Particle
155 // ISF: Energy, mass, and momentum are in MeV, position in mm
156 // Acts: Energy, mass, and momentum are in GeV, position in mm
157 ATH_MSG_DEBUG(name() << " Convert ISF::Particle(mass) " << isfp->id()<<"|" << *isfp<<"(" << isfp->mass() << ")");
158 std::vector<ActsFatras::Particle> input = std::vector<ActsFatras::Particle>{
159 ActsFatras::Particle(ActsFatras::Barcode().withVertexPrimary(0).withParticle(isfp->id()), static_cast<Acts::PdgParticle>(isfp->pdgCode()),
160 isfp->charge(),isfp->mass() * Acts::UnitConstants::MeV)
161 .setDirection(Acts::makeDirectionFromPhiEta(isfp->momentum().phi(), isfp->momentum().eta()))
162 .setAbsoluteMomentum(isfp->momentum().mag() * Acts::UnitConstants::MeV)
163 .setPosition4(ActsTrk::convertPosToActs(isfp->position(), isfp->timeStamp()))};
164 ATH_MSG_DEBUG(name() << " Propagating ActsFatras::Particle vertex|particle|generation|subparticle, " << input[0]);
165 std::vector<ActsFatras::Particle> simulatedInitial;
166 std::vector<ActsFatras::Particle> simulatedFinal;
167 std::vector<ActsFatras::Hit> hits;
168 // simulate
169 auto result=simulator.simulate(anygctx, mctx, generator, input, simulatedInitial, simulatedFinal, hits);
170 auto simulatedFailure=result.value();
171 if (simulatedFailure.size()>0){
172 for (const auto& simfail : simulatedFailure){
173 auto errCode = Acts::make_error_code(Acts::PropagatorError(simfail.error.value()));
174 ATH_MSG_WARNING(name() << " Particle id " <<simfail.particle.particleId()<< ": fail to be simulated during Propagation: " << errCode.message());
175 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()));
176 return StatusCode::SUCCESS;
177 }
178 }
179
180 ATH_MSG_DEBUG(name() << " initial particle " << simulatedInitial[0]);
181 ATH_MSG_DEBUG(name() << " ActsFatras simulator hits: " << hits.size());
182 int i = 0;
183 for (const auto& hit : hits) {
184 ATH_MSG_DEBUG(name() << " hit pos: " << hit.position() );
185 ++i;
186 if (i>5) break;
187 }
188 ATH_MSG_DEBUG(name() << " No. of particles after ActsFatras simulator: " << simulatedFinal.size());
189 if (!simulatedFinal.empty()){
190 ATH_MSG_DEBUG(name() << " start procesing secondaries");
191 auto itr = simulatedFinal.begin();
192 // Save hits of isfp
193 std::vector<ActsFatras::Hit> particle_hits;
194 if (itr->numberOfHits() > 0) {
195 std::copy(hits.begin(), hits.begin()+itr->numberOfHits(), std::back_inserter(particle_hits));
197 }
198 // Process secondaries
199 auto isKilled = !itr->isAlive();
200 int maxGeneration = simulatedFinal.back().particleId().generation();
201 ATH_MSG_DEBUG(name() << " maxGeneration: "<< maxGeneration);
202 for (int gen = 0; gen <= maxGeneration; ++gen){
203 ATH_MSG_DEBUG(name() << " start with generation "<< gen << "|" << maxGeneration << ": "<< *itr);
204 auto vecsecisfp = std::make_unique<ISF::ISFParticleVector>();
205 std::unique_ptr<ISF::ISFParticle> newisfp = nullptr; // Boundary crossing particle
206
207 while (itr != simulatedFinal.end() && static_cast<int>(itr->particleId().generation()) == gen) {
208 ATH_MSG_DEBUG(name() << " genration "<< gen << "|" << maxGeneration << ": "<< *itr);
209 if(itr->isSecondary()){
210 // convert final particles to ISF::particle
211 const auto pos = ActsTrk::convertPosFromActs(itr->fourPosition()).first;
212 const auto mom = ActsTrk::convertMomFromActs(itr->fourMomentum()).first;
213 double mass = itr->mass() / Acts::UnitConstants::MeV;
214 double charge = itr->charge();
215 int pdgid = itr->pdg();
216 auto properTime = ActsTrk::timeToAthena(itr->time());
217 const int status = 1 + HepMC::SIM_STATUS_THRESHOLD;
218 const int id = HepMC::UNDEFINED_ID;
219 auto secisfp = std::make_unique<ISF::ISFParticle>(pos,mom,mass,charge,pdgid,status,properTime,*isfp,id);
220 secisfp->setNextGeoID(m_geoIDSvc->identifyNextGeoID(*secisfp));
221 ATH_MSG_DEBUG(name() <<" secondaries particle (ACTS): "<<*itr<< "("<<itr->momentum()<<")|time "<<itr->time()<<"|process "<< getATLASProcessCode(itr->process()));
222 ATH_MSG_DEBUG(name() <<" secondaries particle (ISF): pdg=" << secisfp->pdgCode()
223 << " pos=" << secisfp->position() << " mom=" << secisfp->momentum()
224 << " GeoID=" << m_geoIDSvc->identifyNextGeoID(*secisfp));
225 vecsecisfp->push_back(secisfp.release());
226 }
227 else{
228 // Primary particle handling
229 ATH_MSG_DEBUG(name() <<" primary particle found with generation ("<< gen <<")");
230 // After simulation, check particle's final state
231 if (!isKilled) {
232 auto fisfp = std::make_unique<ISF::ISFParticle>(*isfp);
233 fisfp->updateMomentum(ActsTrk::convertMomFromActs(itr->fourMomentum()).first);
234 fisfp->updatePosition(ActsTrk::convertPosFromActs(itr->fourPosition()).first);
235 ATH_MSG_DEBUG(name() << " After simulation, primary particle state: " << *fisfp);
236 if (!m_particleFilter.empty() && !m_particleFilter->passFilter(*fisfp)) {
237 ATH_MSG_VERBOSE("ISFParticle" << fisfp << " after simulation does not pass selection. Ignoring for boundary check.");
238 continue;
239 }
240 ATH_MSG_DEBUG(name() << " [ISF] original GeoID: " << m_geoIDSvc->identifyGeoID(*isfp)
241 << " new particle GeoID: " << m_geoIDSvc->identifyGeoID(*fisfp)
242 << ", nextGeoID: " << m_geoIDSvc->identifyNextGeoID(*fisfp));
243
244 // Use ActsExtrapolationTool
245 ATH_MSG_DEBUG(name() << " Extrapolating using ActsExtrapolationTool");
246
247 // Convert to ACTS BoundTrackParameters for extrapolation
248 Acts::BoundTrackParameters startParams = Acts::BoundTrackParameters::createCurvilinear(
249 itr->fourPosition(), itr->direction(), itr->qOverP(), std::nullopt, itr->hypothesis());
250
251 //=============== try if a starting surface exist before passing to Extrapolator ==
252 auto do_exit_startsurface = checkStartSurface(mctx, anygctx, surfaceCheckPropagator, startParams);
253 ATH_MSG_DEBUG(name() << " checkStartSurface returned: " << do_exit_startsurface);
254 if (do_exit_startsurface) {
255 ATH_MSG_DEBUG(name() << " Particle starts at a valid surface, doing extrapolation...");
256
257 // Extrapolate and get propagation steps
259 Amg::Vector3D entryPos{Amg::Vector3D::Zero()};
260 try {
261 auto stepsResult = m_extrapolationTool->propagationSteps(ctx, startParams, Acts::Direction::Forward());
262 auto steps = stepsResult.value().first;
263 ATH_MSG_DEBUG(name() << " Number of propagation steps: " << steps.size());
264 if (steps.size() != 0) {
265 for (const auto& step : steps) {
266 ATH_MSG_DEBUG(name() << " [Acts] Step at position " << step.position
267 << " (eta " << Acts::VectorHelpers::eta(step.position)
268 << ") with GeoID " << step.geoID);
269 entryPos = convertPos3FromActs(step.position);
270 nextGeoID = m_geoIDSvc->identifyGeoID(entryPos);
271 ATH_MSG_DEBUG(name() << " [Acts] GeoID from service: " << nextGeoID);
272 if (nextGeoID > AtlasDetDescr::fAtlasID) { // Valid boundary crossing
273 ATH_MSG_DEBUG(name() << " Boundary crossing detected at GeoID " << nextGeoID);
274 break;
275 }
276 }
277 } else {
278 ATH_MSG_WARNING(name() << " No propagation steps returned by ActsExtrapolationTool");
279 }
280 }
281 catch (const std::exception& e) {
282 ATH_MSG_WARNING(name() << " extrapolation [" << m_extrapolationTool.name() << "] failed: " << e.what() << "\nSkip boundary check for " << *fisfp);
283 break; // Skip boundary check for this particle and continue with next one
284 }
285
286 if (fisfp && nextGeoID > AtlasDetDescr::fAtlasID){
287 const auto mom = ActsTrk::convertMomFromActs(itr->fourMomentum()).first;
288 double mass = itr->mass() / Acts::UnitConstants::MeV;
289 double charge = itr->charge();
290 int pdgid = itr->pdg();
291 auto properTime = ActsTrk::timeToAthena(itr->time());
292
293 // Create boundary crossing particle
294 newisfp = std::make_unique<ISF::ISFParticle>(entryPos, mom, mass, charge, pdgid, isfp->status(), properTime, *isfp, isfp->id(), isfp->barcode());
295 newisfp->setNextGeoID(nextGeoID);
296 ATH_MSG_DEBUG(name() << " Truthbinding of parent ISFParticle: " << (isfp->getTruthBinding() ? "exists" : "null"));
297 if (isfp->getTruthBinding()) {
298 ATH_MSG_DEBUG(name() << " Current GenParticle: " << isfp->getTruthBinding()->getCurrentGenParticle());
299 }
300 ATH_MSG_DEBUG(name() << " Created new ISFParticle at boundary with nextGeoID: "
302 << "(" << nextGeoID << ")");
303 }
304
305 // Handle boundary crossing particle separately - DON'T add to vecsecisfp yet
306 if (newisfp && nextGeoID > AtlasDetDescr::fAtlasID) {
307 ATH_MSG_DEBUG(name() << " [ISF] Processing boundary particle with nextGeoID: "
308 << AtlasDetDescr::AtlasRegionHelper::getName(newisfp->nextGeoID())
309 << "(" << newisfp->nextGeoID() << ")");
310
311 // Identify Entrylayer
313
314 switch(nextGeoID) {
316 entryLayer = ISF::fAtlasCaloEntry;
317 ATH_MSG_DEBUG("Particle crossing to Calorimeter");
318 break;
320 entryLayer = ISF::fAtlasMuonEntry;
321 ATH_MSG_DEBUG("Particle crossing to Muon System");
322 break;
323 default:
324 ATH_MSG_DEBUG("Particle at unspecified boundary");
325 break;
326 }
327
328 if (entryLayer != ISF::fUnsetEntryLayer) {
329 vecsecisfp->push_back(newisfp.release()); // Add boundary crossing particle to secondaries vector
330 } else {
331 ATH_MSG_WARNING("Invalid entry layer for boundary particle");
332 }
333 }
334 }
335 else {
336 ATH_MSG_DEBUG(name() << " No starting surface found, skipping boundary check and extrapolation.");
337 }
338 } // end of !isKilled
339 } // end of primary vs secondary
340 ++itr;
341 } // end of while loop over particles in generation
342
343 // Process truth for this generation
344 if (!vecsecisfp->empty()) {
345 // Determine process code and geoID based on whether we have boundary crossing
346 int processCode = 0;
348 auto isParentKilled = ISF::fPrimarySurvives;
349
350 if (newisfp) {
351 // Boundary crossing - use boundary info
352 processCode = 91; // Boundary crossing has no process
353 geoID = newisfp->nextGeoID();
354 isParentKilled = ISF::fKillsPrimary;
355 } else {
356 // Regular secondaries - use process from last particle
357 processCode = getATLASProcessCode((itr-1)->process());
358 geoID = (isfp->nextGeoID() <= AtlasDetDescr::fUndefinedAtlasRegion) ?
359 AtlasDetDescr::fAtlasID : isfp->nextGeoID();
360 isParentKilled = isKilled && gen==maxGeneration ? ISF::fKillsPrimary : ISF::fPrimarySurvives;
361 }
362
363 ISF::ISFTruthIncident truth(*isfp,
364 *vecsecisfp,
365 processCode,
366 geoID,
367 isParentKilled);
368
369 ATH_MSG_DEBUG(name() << " Truth incident parentPt2(MinPt2) " << truth.parentPt2() <<" (100 MeV)");
370 ATH_MSG_DEBUG(name() << " Truth incident ChildPt2(MinPt2) " << truth.childrenPt2Pass(300) <<" (300 MeV)");
371 m_truthRecordSvc->registerTruthIncident(truth, true);
372 truth.updateParentAfterIncidentProperties();
373 truth.updateChildParticleProperties();
374 for (auto *secisfp : *vecsecisfp){
375 if (secisfp->getTruthBinding()) {
376 secondaries.push_back(secisfp);
377 ATH_MSG_DEBUG(name() << " Secondary particle written out to truth.\n Parent ("
378 << *isfp << ")\n Secondary (" << *secisfp <<")");
379
380 ATH_MSG_DEBUG("Secondary particle push back to ISF, TruthBinding: " << secisfp->getTruthBinding()->getCurrentGenParticle() << " (current) | " << secisfp->getTruthBinding()->getPrimaryGenParticle() << " (primary) | " << secisfp->getTruthBinding()->getGenerationZeroGenParticle() << " (zero)");
381 if (secisfp->getTruthBinding()->getCurrentGenParticle() != nullptr) ATH_MSG_DEBUG("Secondary particle GenParticle EndVertex: " << (secisfp->getTruthBinding()->getCurrentGenParticle()->end_vertex() ? HepMC::barcode(secisfp->getTruthBinding()->getCurrentGenParticle()->end_vertex()) : 1));
382 } else {
383 ATH_MSG_WARNING("Secondary particle not written out to truth.\n Parent ("
384 << *isfp << ")\n Secondary (" << *secisfp <<")");
385 delete secisfp; // Clean up particles without truth binding
386 }
387 }
388 }
389 } // end of generation loop
390 } // end of !simulatedFinal.empty()
391 ATH_MSG_VERBOSE(name() << " No. of secondaries: " << secondaries.size());
392 ATH_MSG_DEBUG(name() << " End of particle " << isfp->id());
393
394 std::vector<ActsFatras::Particle>().swap(input);
395 std::vector<ActsFatras::Particle>().swap(simulatedInitial);
396 std::vector<ActsFatras::Particle>().swap(simulatedFinal);
397 std::vector<ActsFatras::Hit>().swap(hits);
398 } // end of isfp loop
399 return StatusCode::SUCCESS;
400}
double charge(const T &p)
Definition AtlasPID.h:997
static const char * getName(int region)
Gaudi::Property< double > m_maxStepSize
Acts::EigenStepper< Acts::EigenStepperDefaultExtension > ChargedStepper
Gaudi::Property< double > m_tolerance
Gaudi::Property< double > m_interact_minPt
ActsFatras::Simulation< ChargedSelector, ChargedSimulation, NeutralSelector, NeutralSimulation > Simulation
bool checkStartSurface(const Acts::MagneticFieldContext &mctx, const Acts::GeometryContext &anygctx, const ChargedPropagator &chargedPropagator, const Acts::BoundTrackParameters &startParameters, Acts::Direction navDir=Acts::Direction::Forward(), double pathLimit=std::numeric_limits< double >::max()) const
Gaudi::Property< double > m_maxStep
int getATLASProcessCode(ActsFatras::ProcessType actspt)
Gaudi::Property< double > m_pathLimit
Gaudi::Property< double > m_stepSizeCutOff
Gaudi::Property< bool > m_loopProtection
Gaudi::Property< double > m_maxRungeKuttaStepTrials
virtual Acts::MagneticFieldContext getMagneticFieldContext(const EventContext &) const
Gaudi::Property< double > m_loopFraction
Amg::Vector3D convertPos3FromActs(const Acts::Vector3 &actsPos)
Acts::StraightLineStepper NeutralStepper
const std::string process
constexpr double timeToAthena(const double actsT)
Converts a time unit from Acts to Athena units.
std::pair< Amg::Vector3D, double > convertMomFromActs(const Acts::Vector4 &actsMom)
Converts an Acts four-momentum vector into an pair of an Athena three-momentum and the paritcle's ene...
std::pair< Amg::Vector3D, double > convertPosFromActs(const Acts::Vector4 &actsPos)
Converts an Acts 4-vector into a pair of an Athena spatial vector and the passed time.
Acts::Vector4 convertPosToActs(const Amg::Vector3D &athenaPos, const double athenaTime=0.)
Converts a position vector & time from Athena units into Acts units.
AtlasRegion
A simple enum of ATLAS regions and sub-detectors.
Definition AtlasRegion.h:21
int barcode(const T *p)
Definition Barcode.h:16
constexpr int UNDEFINED_ID
constexpr int SIM_STATUS_THRESHOLD
Constant definiting the status threshold for simulated particles, eg. can be used to separate generat...
EntryLayer
Identifiers for the TrackRecordCollections on the boundaries between CaloEntry: Inner Detector - Calo...
Definition EntryLayer.h:31
@ fUnsetEntryLayer
Definition EntryLayer.h:33
@ fAtlasCaloEntry
Definition EntryLayer.h:37
@ fAtlasMuonEntry
Definition EntryLayer.h:38
@ fKillsPrimary
@ fPrimarySurvives
generator
Configure Herwig7 These are the commands corresponding to what would go into the regular Herwig infil...
status
Definition merge.py:16

◆ sysInitialize()

virtual StatusCode ISF::BaseSimulatorTool::sysInitialize ( )
inlineoverridevirtualinherited

Gaudi sysInitialize() methods.

Definition at line 48 of file BaseSimulatorTool.h.

49 {
51 ATH_CHECK( m_chrono.retrieve() );
52 return StatusCode::SUCCESS;
53 }
virtual StatusCode sysInitialize() override

Member Data Documentation

◆ ATLAS_THREAD_SAFE

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

Definition at line 277 of file ActsFatrasSimTool.h.

277{};

◆ m_ActsFatrasWriteHandler

ToolHandle<ActsFatrasWriteHandler> ISF::ActsFatrasSimTool::m_ActsFatrasWriteHandler
private
Initial value:
{
this, "ActsFatrasWriteHandler", "ActsFatrasWriteHandler"}

Definition at line 305 of file ActsFatrasSimTool.h.

305 {
306 this, "ActsFatrasWriteHandler", "ActsFatrasWriteHandler"};

◆ m_chrono

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

The timing service for general usage.

Definition at line 114 of file BaseSimulatorTool.h.

114{this, "ChronoStatService", "ChronoStatSvc"};

◆ m_extrapolationTool

PublicToolHandle<ActsTrk::IExtrapolationTool> ISF::ActsFatrasSimTool::m_extrapolationTool {this, "ExtrapolationTool", "ActsExtrapolationTool"}
private

Definition at line 285 of file ActsFatrasSimTool.h.

285{this, "ExtrapolationTool", "ActsExtrapolationTool"};

◆ m_fieldCacheCondObjInputKey

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

Definition at line 293 of file ActsFatrasSimTool.h.

293{this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"};

◆ m_geoIDSvc

ServiceHandle<ISF::IGeoIDSvc> ISF::ActsFatrasSimTool::m_geoIDSvc {this, "GeoIDSvc", "ISF::GeoIDSvc"}
private

Definition at line 282 of file ActsFatrasSimTool.h.

282{this, "GeoIDSvc", "ISF::GeoIDSvc"};

◆ 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 313 of file ActsFatrasSimTool.h.

313{this, "IncludeGgradient", true, "Boolean flag for inclusion of d(dEds)d(q/p) into energy loss"};

◆ 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 308 of file ActsFatrasSimTool.h.

308 {this, "Interact_MinPt", 50.0,
309 "Min pT of the interactions (MeV)"};

◆ m_logger

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

Definition at line 296 of file ActsFatrasSimTool.h.

296{nullptr};

◆ 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 326 of file ActsFatrasSimTool.h.

326 {this, "LoopFraction", 0.5,
327 "Allowed loop fraction, 1 is a full loop"};

◆ m_loopProtection

Gaudi::Property<bool> ISF::ActsFatrasSimTool::m_loopProtection
private
Initial value:
{this, "LoopProtection", true,
"Loop protection, it adapts the pathLimit"}

Definition at line 324 of file ActsFatrasSimTool.h.

324 {this, "LoopProtection", true,
325 "Loop protection, it adapts the pathLimit"};

◆ 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 318 of file ActsFatrasSimTool.h.

318 {this, "MaxRungeKuttaStepTrials", 10000,
319 "Maximum number of Runge-Kutta steps for the stepper step call"};

◆ m_maxStep

Gaudi::Property<double> ISF::ActsFatrasSimTool::m_maxStep
private
Initial value:
{this, "MaxSteps", 1000,
"Max number of steps"}

Definition at line 316 of file ActsFatrasSimTool.h.

316 {this, "MaxSteps", 1000,
317 "Max number of steps"};

◆ 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 320 of file ActsFatrasSimTool.h.

320 {this, "MaxStepSize", 3.0,
321 "Max step size (converted to Acts::UnitConstants::m)"};

◆ 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 312 of file ActsFatrasSimTool.h.

312{this, "MeanEnergyLoss", true, "Toggle between mean and mode evaluation of energy loss"};

◆ 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 314 of file ActsFatrasSimTool.h.

314{this, "MomentumCutOff", 0., "Cut-off value for the momentum in SI units"};

◆ m_particleFilter

PublicToolHandle<ISF::IParticleFilter> ISF::ActsFatrasSimTool::m_particleFilter
private
Initial value:
{
this, "ParticleFilter", "", "Particle filter kinematic cuts, etc."}

Definition at line 299 of file ActsFatrasSimTool.h.

299 {
300 this, "ParticleFilter", "", "Particle filter kinematic cuts, etc."};

◆ 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 322 of file ActsFatrasSimTool.h.

322 {this, "PathLimit", 100.0,
323 "Track path limit (converted to Acts::UnitConstants::cm)"};

◆ m_pixelSiHits

SiHitCollection ISF::ActsFatrasSimTool::m_pixelSiHits
private

Definition at line 256 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 334 of file ActsFatrasSimTool.h.

334 {this, "ProcessTypeMap",
335 {{0,0}, {1,201}, {2,14}, {3,3}, {4,121}}, "proessType map <ActsFatras,G4>"};

◆ m_randomEngineName

Gaudi::Property<std::string> ISF::ActsFatrasSimTool::m_randomEngineName
private
Initial value:
{this, "RandomEngineName",
"RandomEngineName", "Name of random number stream"}

Definition at line 278 of file ActsFatrasSimTool.h.

278 {this, "RandomEngineName",
279 "RandomEngineName", "Name of random number stream"};

◆ m_rngSvc

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

Definition at line 276 of file ActsFatrasSimTool.h.

276{this, "RNGService", "AthRNGSvc"};

◆ m_sctSiHits

SiHitCollection ISF::ActsFatrasSimTool::m_sctSiHits
private

Definition at line 257 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 330 of file ActsFatrasSimTool.h.

330 {this, "StepSizeCutOff", 0.,
331 "Cut-off value for the step size"};

◆ 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 328 of file ActsFatrasSimTool.h.

328 {this, "Tolerance", 0.0001,
329 "Tolerance for the error of the integration"};

◆ m_trackingGeometry

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

Definition at line 290 of file ActsFatrasSimTool.h.

◆ m_trackingGeometryTool

PublicToolHandle<ActsTrk::ITrackingGeometryTool> ISF::ActsFatrasSimTool::m_trackingGeometryTool
private
Initial value:
{
this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}

Definition at line 288 of file ActsFatrasSimTool.h.

288 {
289 this, "TrackingGeometryTool", "ActsTrackingGeometryTool"};

◆ m_truthRecordSvc

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

Definition at line 302 of file ActsFatrasSimTool.h.

302{this, "TruthRecordService", "ISF_TruthRecordSvc", ""};

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