The core class of the new SimplifiedGeometryPropagator. More...

Inheritance diagram for FastSimProducer:

Public Member Functions
	FastSimProducer (const edm::ParameterSet &)

	~FastSimProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	beginStream (edm::StreamID id) override

virtual FSimTrack	createFSimTrack (fastsim::Particle particle, fastsim::ParticleManager particleManager, HepPDT::ParticleDataTable const &particleTable)

void	endStream () override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::unique_ptr < RandomEngineAndDistribution >	_randomEngine
	The random engine. More...

double	beamPipeRadius_
	The radius of the beampipe. More...

fastsim::Geometry	caloGeometry_
	Hack to interface "old" calo to "new" tracking. More...

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryESToken_

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	caloTopologyESToken_

fastsim::Decayer	decayer_
	Handles decays of non-stable particles using pythia. More...

double	deltaRchargedMother_
	Cut on deltaR for ClosestChargedDaughter algorithm (FastSim tracking) More...

edm::EDGetTokenT < edm::HepMCProduct >	genParticlesToken_
	Token to get the genParticles. More...

fastsim::Geometry	geometry_
	The definition of the tracker according to python config. More...

std::map< std::string, fastsim::InteractionModel * >	interactionModelMap_
	Each interaction model has a unique name. More...

std::vector< std::unique_ptr < fastsim::InteractionModel > >	interactionModels_
	All defined interaction models. More...

std::unique_ptr < CalorimetryManager >	myCalorimetry

const edm::ESGetToken < HepPDT::ParticleDataTable, edm::DefaultRecord >	particleDataTableESToken_

fastsim::ParticleFilter	particleFilter_
	Decides which particles have to be propagated. More...

bool	simulateCalorimetry

bool	simulateMuons

edm::ESWatcher < CaloGeometryRecord >	watchCaloGeometry_

edm::ESWatcher < CaloTopologyRecord >	watchCaloTopology_

Static Private Attributes
static const std::string	MESSAGECATEGORY = "FastSimulation"
	Category of debugging messages ("FastSimulation") More...

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

The core class of the new SimplifiedGeometryPropagator.

Coordinates the propagation of all particles, this means it does the following loop: 1) Get particle from ParticleManager 2) Call LayerNavigator to move particle to next intersection with layer 3) Loop over all the interactions and add secondaries to the event 4) Repeat steps 2), 3) until particle left the tracker, lost all its energy or is about to decay 5) If particle is about to decay: do decay and add secondaries to the event 6) Restart from 1) with the next particle 7) If last particle was propagated add SimTracks, SimVertices, SimHits,... to the event

Definition at line 64 of file FastSimProducer.cc.

Constructor & Destructor Documentation

FastSimProducer::FastSimProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 102 of file FastSimProducer.cc.

References caloGeometryESToken_, caloTopologyESToken_, beamerCreator::create(), DeDxTools::esConsumes(), Exception, get, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNames(), interactionModelMap_, interactionModels_, eostools::move(), myCalorimetry, simulateCalorimetry, and AlCaHLTBitMon_QueryRunRegistry::string.

     : genParticlesToken_(consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("src"))),
       geometry_(iConfig.getParameter<edm::ParameterSet>("trackerDefinition"), consumesCollector()),
       caloGeometry_(iConfig.getParameter<edm::ParameterSet>("caloDefinition"), consumesCollector()),
       beamPipeRadius_(iConfig.getParameter<double>("beamPipeRadius")),
       deltaRchargedMother_(iConfig.getParameter<double>("deltaRchargedMother")),
       particleFilter_(iConfig.getParameter<edm::ParameterSet>("particleFilter")),
       _randomEngine(nullptr),
       simulateCalorimetry(iConfig.getParameter<bool>("simulateCalorimetry")),
       simulateMuons(iConfig.getParameter<bool>("simulateMuons")),
       particleDataTableESToken_(esConsumes()) {
   if (simulateCalorimetry) {
     caloGeometryESToken_ = esConsumes();
     caloTopologyESToken_ = esConsumes();
   }
 
   //----------------
   // define interaction models
   //---------------
 
   const edm::ParameterSet& modelCfgs = iConfig.getParameter<edm::ParameterSet>("interactionModels");
   for (const std::string& modelName : modelCfgs.getParameterNames()) {
     const edm::ParameterSet& modelCfg = modelCfgs.getParameter<edm::ParameterSet>(modelName);
     std::string modelClassName(modelCfg.getParameter<std::string>("className"));
     // Use plugin-factory to create model
     std::unique_ptr<fastsim::InteractionModel> interactionModel(
         fastsim::InteractionModelFactory::get()->create(modelClassName, modelName, modelCfg));
     if (!interactionModel.get()) {
       throw cms::Exception("FastSimProducer") << "InteractionModel " << modelName << " could not be created";
     }
     // Add model to list
     interactionModels_.push_back(std::move(interactionModel));
     // and create the map
     interactionModelMap_[modelName] = interactionModels_.back().get();
   }
 
   //----------------
   // calorimetry
   //---------------
 
   if (simulateCalorimetry) {
     myCalorimetry =
         std::make_unique<CalorimetryManager>(nullptr,
                                              iConfig.getParameter<edm::ParameterSet>("Calorimetry"),
                                              iConfig.getParameter<edm::ParameterSet>("MaterialEffectsForMuonsInECAL"),
                                              iConfig.getParameter<edm::ParameterSet>("MaterialEffectsForMuonsInHCAL"),
                                              iConfig.getParameter<edm::ParameterSet>("GFlash"),
                                              consumesCollector());
   }
 
   //----------------
   // register products
   //----------------
 
   // SimTracks and SimVertices
   produces<edm::SimTrackContainer>();
   produces<edm::SimVertexContainer>();
   // products of interaction models, i.e. simHits
   for (auto& interactionModel : interactionModels_) {
     interactionModel->registerProducts(producesCollector());
   }
   produces<edm::PCaloHitContainer>("EcalHitsEB");
   produces<edm::PCaloHitContainer>("EcalHitsEE");
   produces<edm::PCaloHitContainer>("EcalHitsES");
   produces<edm::PCaloHitContainer>("HcalHits");
   produces<edm::SimTrackContainer>("MuonSimTracks");
 }

FastSimProducer::~FastSimProducer ( )

inlineoverride

Definition at line 67 of file FastSimProducer.cc.

67 { ; }

Member Function Documentation

void FastSimProducer::beginStream ( edm::StreamID id )

overrideprivate

Definition at line 170 of file FastSimProducer.cc.

References _randomEngine, and gpuClustering::id.

                                                       {
   _randomEngine = std::make_unique<RandomEngineAndDistribution>(id);
 }

FSimTrack FastSimProducer::createFSimTrack	(	fastsim::Particle *	particle,
		fastsim::ParticleManager *	particleManager,
		HepPDT::ParticleDataTable const &	particleTable
	)

privatevirtual

Definition at line 379 of file FastSimProducer.cc.

References _randomEngine, fastsim::ParticleManager::addSecondaries(), caloGeometry_, fastsim::Particle::charge(), fastsim::Decayer::decay(), decayer_, fastsim::SimplifiedGeometry::ECAL, fastsim::Particle::genParticleIndex(), fastsim::SimplifiedGeometry::getCaloType(), fastsim::ParticleManager::getSimTrack(), fastsim::ParticleManager::getSimVertex(), fastsim::SimplifiedGeometry::getThickness(), fastsim::SimplifiedGeometry::HCAL, fastsim::SimplifiedGeometry::isForward(), fastsim::Particle::isStable(), LogDebug, makeParticle(), MESSAGECATEGORY, fastsim::Particle::momentum(), fastsim::LayerNavigator::moveParticleToNextLayer(), FSimTrack::onEcal(), FSimTrack::onHcal(), FSimTrack::onLayer1(), FSimTrack::onLayer2(), FSimTrack::onVFcal(), fastsim::Particle::pdgId(), fastsim::Particle::position(), fastsim::SimplifiedGeometry::PRESHOWER1, fastsim::SimplifiedGeometry::PRESHOWER2, fastsim::Particle::remainingProperLifeTimeC(), FSimTrack::setEcal(), FSimTrack::setGlobal(), FSimTrack::setHcal(), FSimTrack::setLayer1(), FSimTrack::setLayer2(), FSimTrack::setVFcal(), fastsim::Particle::simTrackIndex(), fastsim::Particle::simVertexIndex(), summarizeEdmComparisonLogfiles::success, and fastsim::SimplifiedGeometry::VFCAL.

Referenced by produce().

                                                                                          {
   FSimTrack myFSimTrack(particle->pdgId(),
                         particleManager->getSimTrack(particle->simTrackIndex()).momentum(),
                         particle->simVertexIndex(),
                         particle->genParticleIndex(),
                         particle->simTrackIndex(),
                         particle->charge(),
                         particle->position(),
                         particle->momentum(),
                         particleManager->getSimVertex(particle->simVertexIndex()));
 
   // move the particle through the caloLayers
   fastsim::LayerNavigator caloLayerNavigator(caloGeometry_);
   const fastsim::SimplifiedGeometry* caloLayer = nullptr;
 
   // moveParticleToNextLayer(..) returns 0 in case that particle decays
   // in this case particle is propagated up to its decay vertex
   while (caloLayerNavigator.moveParticleToNextLayer(*particle, caloLayer)) {
     LogDebug(MESSAGECATEGORY) << "   moved to next caloLayer: " << *caloLayer;
     LogDebug(MESSAGECATEGORY) << "   new state: " << *particle;
 
     // break after 25 ns: only happens for particles stuck in loops
     if (particle->position().T() > 50) {
       caloLayer = nullptr;
       break;
     }
 
     // Define ParticlePropagators (RawParticle) needed for CalorimetryManager and save them
 
     RawParticle PP = makeParticle(&particleTable, particle->pdgId(), particle->momentum(), particle->position());
 
     // no material
     if (caloLayer->getThickness(particle->position(), particle->momentum()) < 1E-10) {
       // unfortunately needed for CalorimetryManager
       if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::ECAL) {
         if (!myFSimTrack.onEcal()) {
           myFSimTrack.setEcal(PP, 0);
         }
       } else if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::HCAL) {
         if (!myFSimTrack.onHcal()) {
           myFSimTrack.setHcal(PP, 0);
         }
       } else if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::VFCAL) {
         if (!myFSimTrack.onVFcal()) {
           myFSimTrack.setVFcal(PP, 0);
         }
       }
 
       // not necessary to continue propagation
       if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::VFCAL) {
         myFSimTrack.setGlobal();
         caloLayer = nullptr;
         break;
       }
 
       continue;
     }
 
     // Stupid variable used by the old propagator
     // For details check BaseParticlePropagator.h
     int success = 0;
     if (caloLayer->isForward()) {
       success = 2;
       // particle moves inwards
       if (particle->position().Z() * particle->momentum().Z() < 0) {
         success *= -1;
       }
     } else {
       success = 1;
       // particle moves inwards
       if (particle->momentum().X() * particle->position().X() + particle->momentum().Y() * particle->position().Y() <
           0) {
         success *= -1;
       }
     }
 
     // Save the hit
     if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::PRESHOWER1) {
       if (!myFSimTrack.onLayer1()) {
         myFSimTrack.setLayer1(PP, success);
       }
     }
 
     if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::PRESHOWER2) {
       if (!myFSimTrack.onLayer2()) {
         myFSimTrack.setLayer2(PP, success);
       }
     }
 
     if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::ECAL) {
       if (!myFSimTrack.onEcal()) {
         myFSimTrack.setEcal(PP, success);
       }
     }
 
     if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::HCAL) {
       if (!myFSimTrack.onHcal()) {
         myFSimTrack.setHcal(PP, success);
       }
     }
 
     if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::VFCAL) {
       if (!myFSimTrack.onVFcal()) {
         myFSimTrack.setVFcal(PP, success);
       }
     }
 
     // Particle reached end of detector
     if (caloLayer->getCaloType() == fastsim::SimplifiedGeometry::VFCAL) {
       myFSimTrack.setGlobal();
       caloLayer = nullptr;
       break;
     }
 
     LogDebug(MESSAGECATEGORY) << "--------------------------------"
                               << "\n-------------------------------";
   }
 
   // do decays
   // don't have to worry about daughters if particle already within the calorimetry
   // since they will be rejected by the vertex cut of the ParticleFilter
   if (!particle->isStable() && particle->remainingProperLifeTimeC() < 1E-10) {
     LogDebug(MESSAGECATEGORY) << "Decaying particle...";
     std::vector<std::unique_ptr<fastsim::Particle> > secondaries;
     decayer_.decay(*particle, secondaries, _randomEngine->theEngine());
     LogDebug(MESSAGECATEGORY) << "   decay has " << secondaries.size() << " products";
     particleManager->addSecondaries(particle->position(), particle->simTrackIndex(), secondaries);
   }
 
   return myFSimTrack;
 }

void FastSimProducer::endStream ( )

overrideprivate

Definition at line 377 of file FastSimProducer.cc.

References _randomEngine.

377 { _randomEngine.reset(); }

FastSimProducer::_randomEngine

std::unique_ptr< RandomEngineAndDistribution > _randomEngine

The random engine.

Definition: FastSimProducer.cc:83

void FastSimProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 174 of file FastSimProducer.cc.

References _randomEngine, funct::abs(), beamPipeRadius_, caloGeometry_, caloGeometryESToken_, caloTopologyESToken_, edm::ESWatcher< T >::check(), createFSimTrack(), fastsim::Decayer::decay(), decayer_, deltaRchargedMother_, genParticleCandidates2GenParticles_cfi::genParticles, genParticlesToken_, geometry_, edm::Event::getByToken(), fastsim::SimplifiedGeometry::getCaloType(), edm::EventSetup::getData(), fastsim::SimplifiedGeometry::getInteractionModels(), fastsim::Geometry::getMagneticFieldZ(), fastsim::SimplifiedGeometry::getThickness(), interactionModelMap_, interactionModels_, phase1PixelTopology::layer, LogDebug, MESSAGECATEGORY, eostools::move(), fastsim::LayerNavigator::moveParticleToNextLayer(), myCalorimetry, particleDataTableESToken_, particleFilter_, edm::Event::put(), simulateCalorimetry, simulateMuons, fastsim::SimplifiedGeometry::TRACKERBOUNDARY, fastsim::Geometry::update(), watchCaloGeometry_, and watchCaloTopology_.

                                                                            {
   LogDebug(MESSAGECATEGORY) << "   produce";
 
   geometry_.update(iSetup, interactionModelMap_);
   caloGeometry_.update(iSetup, interactionModelMap_);
 
   // Define containers for SimTracks, SimVertices
   std::unique_ptr<edm::SimTrackContainer> simTracks_(new edm::SimTrackContainer);
   std::unique_ptr<edm::SimVertexContainer> simVertices_(new edm::SimVertexContainer);
 
   // Get the particle data table (in case lifetime or charge of GenParticles not set)
   auto const& pdt = iSetup.getData(particleDataTableESToken_);
 
   // Get the GenParticle collection
   edm::Handle<edm::HepMCProduct> genParticles;
   iEvent.getByToken(genParticlesToken_, genParticles);
 
   // Load the ParticleManager which returns the particles that have to be propagated
   // Creates a fastsim::Particle out of a GenParticle/secondary
   fastsim::ParticleManager particleManager(*genParticles->GetEvent(),
                                            pdt,
                                            beamPipeRadius_,
                                            deltaRchargedMother_,
                                            particleFilter_,
                                            *simTracks_,
                                            *simVertices_);
 
   //  Initialize the calorimeter geometry
   if (simulateCalorimetry) {
     if (watchCaloGeometry_.check(iSetup) || watchCaloTopology_.check(iSetup)) {
       auto const& pG = iSetup.getData(caloGeometryESToken_);
       myCalorimetry->getCalorimeter()->setupGeometry(pG);
 
       auto const& theCaloTopology = iSetup.getData(caloTopologyESToken_);
       myCalorimetry->getCalorimeter()->setupTopology(theCaloTopology);
       myCalorimetry->getCalorimeter()->initialize(geometry_.getMagneticFieldZ(math::XYZTLorentzVector(0., 0., 0., 0.)));
 
       myCalorimetry->getHFShowerLibrary()->initHFShowerLibrary(iSetup);
     }
 
     // Important: this also cleans the calorimetry information from the last event
     myCalorimetry->initialize(_randomEngine.get());
   }
 
   // The vector of SimTracks needed for the CalorimetryManager
   std::vector<FSimTrack> myFSimTracks;
 
   LogDebug(MESSAGECATEGORY) << "################################"
                             << "\n###############################";
 
   // loop over particles
   for (std::unique_ptr<fastsim::Particle> particle = particleManager.nextParticle(*_randomEngine); particle != nullptr;
        particle = particleManager.nextParticle(*_randomEngine)) {
     LogDebug(MESSAGECATEGORY) << "\n   moving NEXT particle: " << *particle;
 
     // -----------------------------
     // This condition is necessary because of hack for calorimetry
     // -> The CalorimetryManager should also be implemented based on this new FastSim classes (Particle.h) in a future project.
     // A second loop (below) loops over all parts of the calorimetry in order to create a track of the old FastSim class FSimTrack.
     // The condition below (R<128, z<302) makes sure that the particle geometrically is outside the tracker boundaries
     // -----------------------------
 
     if (particle->position().Perp2() < 128. * 128. && std::abs(particle->position().Z()) < 302.) {
       // move the particle through the layers
       fastsim::LayerNavigator layerNavigator(geometry_);
       const fastsim::SimplifiedGeometry* layer = nullptr;
 
       // moveParticleToNextLayer(..) returns 0 in case that particle decays
       // in this case particle is propagated up to its decay vertex
       while (layerNavigator.moveParticleToNextLayer(*particle, layer)) {
         LogDebug(MESSAGECATEGORY) << "   moved to next layer: " << *layer;
         LogDebug(MESSAGECATEGORY) << "   new state: " << *particle;
 
         // Hack to interface "old" calo to "new" tracking
         // Particle reached calorimetry so stop further propagation
         if (layer->getCaloType() == fastsim::SimplifiedGeometry::TRACKERBOUNDARY) {
           layer = nullptr;
           // particle no longer is on a layer
           particle->resetOnLayer();
           break;
         }
 
         // break after 25 ns: only happens for particles stuck in loops
         if (particle->position().T() > 25) {
           layer = nullptr;
           // particle no longer is on a layer
           particle->resetOnLayer();
           break;
         }
 
         // perform interaction between layer and particle
         // do only if there is actual material
         if (layer->getThickness(particle->position(), particle->momentum()) > 1E-10) {
           int nSecondaries = 0;
           // loop on interaction models
           for (fastsim::InteractionModel* interactionModel : layer->getInteractionModels()) {
             LogDebug(MESSAGECATEGORY) << "   interact with " << *interactionModel;
             std::vector<std::unique_ptr<fastsim::Particle> > secondaries;
             interactionModel->interact(*particle, *layer, secondaries, *_randomEngine);
             nSecondaries += secondaries.size();
             particleManager.addSecondaries(particle->position(), particle->simTrackIndex(), secondaries, layer);
           }
 
           // kinematic cuts: particle might e.g. lost all its energy
           if (!particleFilter_.acceptsEn(*particle)) {
             // Add endvertex if particle did not create any secondaries
             if (nSecondaries == 0)
               particleManager.addEndVertex(particle.get());
             layer = nullptr;
             break;
           }
         }
 
         LogDebug(MESSAGECATEGORY) << "--------------------------------"
                                   << "\n-------------------------------";
       }
 
       // do decays
       if (!particle->isStable() && particle->remainingProperLifeTimeC() < 1E-10) {
         LogDebug(MESSAGECATEGORY) << "Decaying particle...";
         std::vector<std::unique_ptr<fastsim::Particle> > secondaries;
         decayer_.decay(*particle, secondaries, _randomEngine->theEngine());
         LogDebug(MESSAGECATEGORY) << "   decay has " << secondaries.size() << " products";
         particleManager.addSecondaries(particle->position(), particle->simTrackIndex(), secondaries);
         continue;
       }
 
       LogDebug(MESSAGECATEGORY) << "################################"
                                 << "\n###############################";
     }
 
     // -----------------------------
     // Hack to interface "old" calorimetry with "new" propagation in tracker
     // The CalorimetryManager has to know which particle could in principle hit which parts of the calorimeter
     // I think it's a bit strange to propagate the particle even further (and even decay it) if it already hits
     // some part of the calorimetry but this is how the code works...
     // -----------------------------
 
     if (particle->position().Perp2() >= 128. * 128. || std::abs(particle->position().Z()) >= 302.) {
       LogDebug(MESSAGECATEGORY) << "\n   moving particle to calorimetry: " << *particle;
 
       // create FSimTrack (this is the object the old propagation uses)
       myFSimTracks.push_back(createFSimTrack(particle.get(), &particleManager, pdt));
       // particle was decayed
       if (!particle->isStable() && particle->remainingProperLifeTimeC() < 1E-10) {
         continue;
       }
 
       LogDebug(MESSAGECATEGORY) << "################################"
                                 << "\n###############################";
     }
 
     // -----------------------------
     // End Hack
     // -----------------------------
 
     LogDebug(MESSAGECATEGORY) << "################################"
                               << "\n###############################";
   }
 
   // store simTracks and simVertices
   iEvent.put(std::move(simTracks_));
   iEvent.put(std::move(simVertices_));
   // store products of interaction models, i.e. simHits
   for (auto& interactionModel : interactionModels_) {
     interactionModel->storeProducts(iEvent);
   }
 
   // -----------------------------
   // Calorimetry Manager
   // -----------------------------
   if (simulateCalorimetry) {
     for (auto myFSimTrack : myFSimTracks) {
       myCalorimetry->reconstructTrack(myFSimTrack, _randomEngine.get());
     }
   }
 
   // -----------------------------
   // Store Hits
   // -----------------------------
   std::unique_ptr<edm::PCaloHitContainer> p4(new edm::PCaloHitContainer);
   std::unique_ptr<edm::PCaloHitContainer> p5(new edm::PCaloHitContainer);
   std::unique_ptr<edm::PCaloHitContainer> p6(new edm::PCaloHitContainer);
   std::unique_ptr<edm::PCaloHitContainer> p7(new edm::PCaloHitContainer);
 
   std::unique_ptr<edm::SimTrackContainer> m1(new edm::SimTrackContainer);
 
   if (simulateCalorimetry) {
     myCalorimetry->loadFromEcalBarrel(*p4);
     myCalorimetry->loadFromEcalEndcap(*p5);
     myCalorimetry->loadFromPreshower(*p6);
     myCalorimetry->loadFromHcal(*p7);
     if (simulateMuons) {
       myCalorimetry->harvestMuonSimTracks(*m1);
     }
   }
   iEvent.put(std::move(p4), "EcalHitsEB");
   iEvent.put(std::move(p5), "EcalHitsEE");
   iEvent.put(std::move(p6), "EcalHitsES");
   iEvent.put(std::move(p7), "HcalHits");
   iEvent.put(std::move(m1), "MuonSimTracks");
 }

Member Data Documentation

std::unique_ptr<RandomEngineAndDistribution> FastSimProducer::_randomEngine

private

The random engine.

Definition at line 83 of file FastSimProducer.cc.

Referenced by beginStream(), createFSimTrack(), endStream(), and produce().

double FastSimProducer::beamPipeRadius_

private

The radius of the beampipe.

Definition at line 80 of file FastSimProducer.cc.

Referenced by produce().

fastsim::Geometry FastSimProducer::caloGeometry_

private

Hack to interface "old" calo to "new" tracking.

Definition at line 79 of file FastSimProducer.cc.

Referenced by createFSimTrack(), and produce().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> FastSimProducer::caloGeometryESToken_

private

Definition at line 96 of file FastSimProducer.cc.

Referenced by FastSimProducer(), and produce().

edm::ESGetToken<CaloTopology, CaloTopologyRecord> FastSimProducer::caloTopologyESToken_

private

Definition at line 97 of file FastSimProducer.cc.

Referenced by FastSimProducer(), and produce().

fastsim::Decayer FastSimProducer::decayer_

private

Handles decays of non-stable particles using pythia.

Definition at line 91 of file FastSimProducer.cc.

Referenced by createFSimTrack(), and produce().

double FastSimProducer::deltaRchargedMother_

private

Cut on deltaR for ClosestChargedDaughter algorithm (FastSim tracking)

Definition at line 81 of file FastSimProducer.cc.

Referenced by produce().

edm::EDGetTokenT<edm::HepMCProduct> FastSimProducer::genParticlesToken_

private

Token to get the genParticles.

Definition at line 77 of file FastSimProducer.cc.

Referenced by produce().

fastsim::Geometry FastSimProducer::geometry_

private

The definition of the tracker according to python config.

Definition at line 78 of file FastSimProducer.cc.

Referenced by produce().

std::map<std::string, fastsim::InteractionModel*> FastSimProducer::interactionModelMap_

private

Each interaction model has a unique name.

Definition at line 93 of file FastSimProducer.cc.

Referenced by FastSimProducer(), and produce().

std::vector<std::unique_ptr<fastsim::InteractionModel> > FastSimProducer::interactionModels_

private

All defined interaction models.

Definition at line 92 of file FastSimProducer.cc.

Referenced by FastSimProducer(), and produce().

const std::string FastSimProducer::MESSAGECATEGORY = "FastSimulation"

staticprivate

Category of debugging messages ("FastSimulation")

Definition at line 94 of file FastSimProducer.cc.

Referenced by createFSimTrack(), and produce().

std::unique_ptr<CalorimetryManager> FastSimProducer::myCalorimetry

private

Definition at line 88 of file FastSimProducer.cc.

Referenced by FastSimProducer(), and produce().

const edm::ESGetToken<HepPDT::ParticleDataTable, edm::DefaultRecord> FastSimProducer::particleDataTableESToken_

private

Definition at line 95 of file FastSimProducer.cc.

Referenced by produce().

fastsim::ParticleFilter FastSimProducer::particleFilter_

private

Decides which particles have to be propagated.

Definition at line 82 of file FastSimProducer.cc.

Referenced by produce().

bool FastSimProducer::simulateCalorimetry

private

Definition at line 85 of file FastSimProducer.cc.

Referenced by FastSimProducer(), and produce().

bool FastSimProducer::simulateMuons

private

Definition at line 89 of file FastSimProducer.cc.

Referenced by produce().

edm::ESWatcher<CaloGeometryRecord> FastSimProducer::watchCaloGeometry_

private

Definition at line 86 of file FastSimProducer.cc.

Referenced by produce().

edm::ESWatcher<CaloTopologyRecord> FastSimProducer::watchCaloTopology_

private

Definition at line 87 of file FastSimProducer.cc.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation