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 = ActsFatras::StandardChargedElectroMagneticInteractions
using	ChargedSimulation = SingleParticleSimulation< ChargedPropagator, ChargedInteractions, HitSurfaceSelector, ActsFatras::NoDecay >
using	NeutralSelector = ActsFatras::NeutralSelector
using	NeutralInteractions = ActsFatras::InteractionList<>
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 &) 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
ServiceHandle< IAthRNGSvc >	m_rngSvc {this, "RNGService", "AthRNGSvc"}
ATHRNG::RNGWrapper *m_randomEngine	ATLAS_THREAD_SAFE {}
Gaudi::Property< std::string >	m_randomEngineName
ToolHandle< IActsTrackingGeometryTool >	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 68 of file ActsFatrasSimTool.h.

Member Typedef Documentation

ChargedInteractions

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

Definition at line 186 of file ActsFatrasSimTool.h.

ChargedPropagator

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

Definition at line 177 of file ActsFatrasSimTool.h.

ChargedSelector

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

Definition at line 185 of file ActsFatrasSimTool.h.

ChargedSimulation

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

Definition at line 188 of file ActsFatrasSimTool.h.

ChargedStepper

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

Definition at line 176 of file ActsFatrasSimTool.h.

Generator

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

Definition at line 172 of file ActsFatrasSimTool.h.

Navigator

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

Definition at line 174 of file ActsFatrasSimTool.h.

NeutralInteractions

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

Definition at line 193 of file ActsFatrasSimTool.h.

NeutralPropagator

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

Definition at line 180 of file ActsFatrasSimTool.h.

NeutralSelector

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

Definition at line 192 of file ActsFatrasSimTool.h.

NeutralSimulation

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

Definition at line 194 of file ActsFatrasSimTool.h.

NeutralStepper

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

Definition at line 179 of file ActsFatrasSimTool.h.

Simulation

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

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

    : BaseSimulatorTool(type, name, parent) {}

~ActsFatrasSimTool()

ISF::ActsFatrasSimTool::~ActsFatrasSimTool ( )

virtual

Definition at line 22 of file ActsFatrasSimTool.cxx.

{}

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

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

getMagneticFieldContext()

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

virtual

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

                                            {
  ATH_CHECK(BaseSimulatorTool::initialize());
  ATH_MSG_INFO("ISF::ActsFatrasSimTool update with ACTS 37.0.0");
  // 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 217 of file ActsFatrasSimTool.h.

                                                                {
    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.

Definition at line 94 of file BaseSimulatorTool.h.

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

retrieveTool()

template<class T >

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

inlineprivate

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

                                                              {
    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.

Definition at line 67 of file BaseSimulatorTool.h.

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

simFlavor()

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

inlineoverridevirtual

Definition at line 219 of file ActsFatrasSimTool.h.

                                                        {
    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 59 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(m_truthRecordSvc->initializeTruthCollection());
  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 77 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->clone(Acts::Logging::Level::FATAL));
  auto neutralPropagator = NeutralPropagator(neutralStepper, navigator, m_logger->clone(Acts::Logging::Level::FATAL));
  ChargedSimulation simulatorCharged(std::move(chargedPropagator), m_logger->clone());
  NeutralSimulation simulatorNeutral(std::move(neutralPropagator), m_logger->clone());
  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 ActsGeometryContext& gctx = m_trackingGeometryTool->getNominalGeometryContext();
  auto anygctx = gctx.context();
  // Loop over ISFParticleVector and process each separately
  ATH_MSG_VERBOSE(name() << " Processing particles in ISFParticleVector.");
  // For sihit creation
  SiHitCollection pixelSiHits;
  SiHitCollection sctSiHits;
  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
    std::vector<ActsFatras::Particle> input = std::vector<ActsFatras::Particle>{
      ActsFatras::Particle(ActsFatras::Barcode().setVertexPrimary(0).setParticle(
                           HepMC::barcode(isfp)), 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(isfp->position().x(), isfp->position().y(),
                    isfp->position().z(), isfp->timeStamp())};
    std::vector<ActsFatras::Particle> simulatedInitial;
    std::vector<ActsFatras::Particle> simulatedFinal;
    std::vector<ActsFatras::Hit> hits;
    ATH_MSG_DEBUG(name() << " Convert ISF::Particle " << isfp->barcode() << " to ActsFatras::Particle " << input[0].particleId().value());
    ATH_MSG_DEBUG(name() << " Propagating ActsFatras::Particle  vertex|particle|generation|subparticle, " << input[0]);
    // 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.size()>1){
      ATH_MSG_DEBUG(name() << " start procesing secondaries");
    // convert final particles to ISF::particle
      auto isAlive = ISF::fPrimarySurvives;
      // int generation = simulatedFinal[-1].particleId().generation();
      int n = 1;
      for (std::vector<ActsFatras::Particle>::iterator itr=simulatedFinal.begin()+1;itr!=simulatedFinal.end();++itr){
        ATH_MSG_DEBUG(name() << " secondaries particle " <<n<<"/"<< simulatedFinal.size()-1<<": "<< *itr);
        ++n;
        const auto pos = Amg::Vector3D(itr->position()[Acts::ePos0],itr->position()[Acts::ePos1],itr->position()[Acts::ePos2]);
        const auto mom = Amg::Vector3D(itr->fourMomentum()[Acts::eMom0] / Acts::UnitConstants::MeV,itr->fourMomentum()[Acts::eMom1] / Acts::UnitConstants::MeV,itr->fourMomentum()[Acts::eMom2] / Acts::UnitConstants::MeV);
        double mass = itr->mass() / Acts::UnitConstants::MeV;
        double charge = itr->charge();
        int pdgid = itr->pdg();
        double properTime = itr->properTime();
        const int status = 1 + HepMC::SIM_STATUS_THRESHOLD;
        const int id = HepMC::UNDEFINED_ID;
        ATH_MSG_DEBUG(name() << " secondaries particle process code " << itr->process());
        auto secisfp = new ISF::ISFParticle (pos,mom,mass,charge,pdgid,status,properTime,*isfp,id);
        auto vecsecisfp = std::make_unique<ISF::ISFParticleVector>();
        vecsecisfp->push_back(secisfp);
        ATH_MSG_DEBUG(name() << " vecsecisfp barcode|p: " << (*vecsecisfp)[0]->barcode() <<"|"<< (*vecsecisfp)[0]->momentum().mag());
        if(!itr->isAlive()) isAlive = ISF::fKillsPrimary;
        ISF::ISFTruthIncident truth(*isfp,
                                    *vecsecisfp,
                                    getATLASProcessCode(itr->process()),
                                    isfp->nextGeoID(),  // inherits from the parent
                                    isAlive
                                    );
        ATH_MSG_DEBUG(name() << " Truth incident physicsProcessCode()" << truth.physicsProcessCode());
        m_truthRecordSvc->registerTruthIncident(truth, true);
        truth.updateChildParticleProperties();
        ATH_MSG_DEBUG(name() << " Create secondariy ISF::Particle " << secisfp->barcode());
        if (secisfp->getTruthBinding()) {
            ATH_MSG_DEBUG(name() << " Save secondariy " << *secisfp);
            secondaries.push_back(secisfp);
        }
        else {
            ATH_MSG_WARNING("Secondary particle not written out to truth.\n Parent (" << isfp << ")\n Secondary (" << *secisfp <<")");
        }
      }
    }//end of secondaries
    ATH_MSG_VERBOSE(name() << " No. of secondaries: " << secondaries.size());
    ATH_MSG_DEBUG(name() << " End of particle " << isfp->barcode());
    m_ActsFatrasWriteHandler->createHits(*isfp, m_trackingGeometry,hits,pixelSiHits,sctSiHits);
 
    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
  std::vector<SiHitCollection> hitcolls;
  hitcolls.push_back(pixelSiHits);
  hitcolls.push_back(sctSiHits);
  ATH_CHECK(m_ActsFatrasWriteHandler->WriteHits(hitcolls,ctx));
  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 238 of file ActsFatrasSimTool.h.

m_ActsFatrasWriteHandler

ToolHandle<ActsFatrasWriteHandler> ISF::ActsFatrasSimTool::m_ActsFatrasWriteHandler

private

Initial value:

{
      this, "ActsFatrasWriteHandler", "ActsFatrasWriteHandler"}

Definition at line 260 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 248 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 268 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 263 of file ActsFatrasSimTool.h.

m_logger

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

private

Definition at line 251 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 281 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 279 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 273 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 271 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 275 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 267 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 269 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 254 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 277 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 288 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 239 of file ActsFatrasSimTool.h.

m_rngSvc

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

private

Definition at line 237 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 285 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 283 of file ActsFatrasSimTool.h.

m_trackingGeometry

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

private

Definition at line 245 of file ActsFatrasSimTool.h.

m_trackingGeometryTool

ToolHandle<IActsTrackingGeometryTool> ISF::ActsFatrasSimTool::m_trackingGeometryTool

private

Initial value:

{
      this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}

Definition at line 243 of file ActsFatrasSimTool.h.

m_truthRecordSvc

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

private

Definition at line 257 of file ActsFatrasSimTool.h.

---

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

• ActsFatrasSimTool.h
• ActsFatrasSimTool.cxx