 |
ATLAS Offline Software
|
Go to the documentation of this file.
5 #ifndef ISF_ACTSTOOLS_ACTSFATRASSIMTOOL_H
6 #define ISF_ACTSTOOLS_ACTSFATRASSIMTOOL_H
9 #include "GaudiKernel/ServiceHandle.h"
10 #include "GaudiKernel/ToolHandle.h"
25 #include "Acts/Utilities/UnitVectors.hpp"
28 #include "Acts/Geometry/GeometryContext.hpp"
29 #include "Acts/MagneticField/MagneticFieldContext.hpp"
30 #include "Acts/EventData/TrackParameters.hpp"
31 #include "Acts/Propagator/Navigator.hpp"
32 #include "Acts/Propagator/EigenStepper.hpp"
33 #include "Acts/Propagator/EigenStepperDefaultExtension.hpp"
34 #include "Acts/Propagator/StraightLineStepper.hpp"
35 #include "Acts/Propagator/detail/SteppingLogger.hpp"
36 #include "Acts/Propagator/ActorList.hpp"
37 #include "Acts/Propagator/Propagator.hpp"
38 #include "Acts/Definitions/ParticleData.hpp"
39 #include "ActsFatras/EventData/ProcessType.hpp"
40 #include "ActsFatras/Kernel/InteractionList.hpp"
41 #include "ActsFatras/Kernel/Simulation.hpp"
42 #include "ActsFatras/Kernel/SimulationResult.hpp"
43 #include "ActsFatras/Physics/Decay/NoDecay.hpp"
44 #include "ActsFatras/Physics/StandardInteractions.hpp"
45 #include "ActsFatras/Physics/ElectroMagnetic/PhotonConversion.hpp"
46 #include "ActsFatras/Selectors/SurfaceSelectors.hpp"
62 class IParticleHelper;
77 bool isSensitive = surface.associatedDetectorElement() !=
nullptr;
88 template <
typename propagator_t,
typename interactions_t,
89 typename hit_surface_selector_t,
typename decay_t>
118 std::shared_ptr<const Acts::Logger> localLogger_)
133 template <
typename generator_t>
134 Acts::Result<ActsFatras::SimulationResult>
simulate(
135 const Acts::GeometryContext &geoCtx,
136 const Acts::MagneticFieldContext &magCtx, generator_t &
generator,
139 ACTS_VERBOSE(
"Using ActsFatrasSimTool simulate()");
141 using SteppingLogger = Acts::detail::SteppingLogger;
142 using SimulationActor = ActsFatras::detail::SimulationActor<generator_t, decay_t, interactions_t, hit_surface_selector_t>;
143 using Result =
typename SimulationActor::result_type;
144 using Actions = Acts::ActorList<SteppingLogger, SimulationActor, Acts::EndOfWorldReached>;
148 PropagatorOptions
options(geoCtx, magCtx);
150 auto &actor =
options.actorList.template get<SimulationActor>();
157 Acts::BoundTrackParameters startPoint = Acts::BoundTrackParameters::createCurvilinear(
169 return result.value().template get<Result>();
189 ActsFatras::StandardChargedElectroMagneticInteractions;
192 ActsFatras::NoDecay>;
198 ActsFatras::NoDecay>;
205 const IInterface*
parent);
213 const EventContext& ctx,
221 return StatusCode::SUCCESS; };
223 std::vector<SiHitCollection> hitcolls;
228 return StatusCode::SUCCESS; };
233 const EventContext&)
const;
242 if (!thandle.empty() && thandle.retrieve().isFailure()) {
244 return StatusCode::FAILURE;
246 return StatusCode::SUCCESS;
253 "RandomEngineName",
"Name of random number stream"};
257 this,
"TrackingGeometryTool",
"ActsTrackingGeometryTool"};
264 std::shared_ptr<const Acts::Logger>
m_logger{
nullptr};
268 this,
"ParticleFilter",
"",
"Particle filter kinematic cuts, etc."};
274 this,
"ActsFatrasWriteHandler",
"ActsFatrasWriteHandler"};
277 "Min pT of the interactions (MeV)"};
280 Gaudi::Property<bool>
m_meanEnergyLoss{
this,
"MeanEnergyLoss",
true,
"Toggle between mean and mode evaluation of energy loss"};
281 Gaudi::Property<bool>
m_includeGgradient{
this,
"IncludeGgradient",
true,
"Boolean flag for inclusion of d(dEds)d(q/p) into energy loss"};
282 Gaudi::Property<double>
m_momentumCutOff{
this,
"MomentumCutOff", 0.,
"Cut-off value for the momentum in SI units"};
284 Gaudi::Property<double>
m_maxStep{
this,
"MaxSteps", 1000,
285 "Max number of steps"};
287 "Maximum number of Runge-Kutta steps for the stepper step call"};
289 "Max step size (converted to Acts::UnitConstants::m)"};
291 "Track path limit (converted to Acts::UnitConstants::cm)"};
293 "Loop protection, it adapts the pathLimit"};
295 "Allowed loop fraction, 1 is a full loop"};
297 "Tolerance for the error of the integration"};
299 "Cut-off value for the step size"};
302 {{0,0}, {1,201}, {2,14}, {3,3}, {4,121}},
"proessType map <ActsFatras,G4>"};
std::list< ISF::ISFParticle * > ISFParticleContainer
generic ISFParticle container (not necessarily a std::list!)
std::shared_ptr< const Acts::Logger > localLogger
Local logger for debug output.
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
const Acts::Logger & logger() const
Provide access to the local logger instance, e.g. for logging macros.
interactions_t interactions
Interaction list containing the simulated interactions.
Acts::Navigator Navigator
SingleParticleSimulation(propagator_t &&propagator_, std::shared_ptr< const Acts::Logger > localLogger_)
Alternatively construct the simulator with an external logger.
Particle_v1 Particle
Define the latest version of the particle class.
propagator_t propagator
How and within which geometry to propagate the particle.
double maxRungeKuttaStepTrials
::StatusCode StatusCode
StatusCode definition for legacy code.
Acts::Result< ActsFatras::SimulationResult > simulate(const Acts::GeometryContext &geoCtx, const Acts::MagneticFieldContext &magCtx, generator_t &generator, const ActsFatras::Particle &particle) const
Simulate a single particle without secondaries.
std::vector< ISF::ISFParticle * > ISFParticleVector
ISFParticle vector.
Single particle simulation with fixed propagator, interactions, and decay.
This defines the McEventCollection, which is really just an ObjectVector of McEvent objects.
A wrapper class for event-slot-local random engines.
hit_surface_selector_t selectHitSurface
Selector for surfaces that should generate hits.
generator
Configure Herwig7 These are the commands corresponding to what would go into the regular Herwig infil...
decay_t decay
Decay module.
double maxStepSize
parameters for propagator options
Define macros for attributes used to control the static checker.
int SimulationFlavor
Identifier type for simulation flavor.
This class provides an interface to generate or decode an identifier for the upper levels of the dete...