#include <FamosProducer.h>

Inheritance diagram for FamosProducer:

Public Member Functions
void	beginRun (edm::Run const &run, const edm::EventSetup &es) override

	FamosProducer (edm::ParameterSet const &p)

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

	~FamosProducer () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
HepMC::GenEvent *	evt_

FamosManager *	famosManager_

bool	simulateMuons

edm::EDGetTokenT< edm::HepMCProduct >	sourceToken

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

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 17 of file FamosProducer.h.

Constructor & Destructor Documentation

FamosProducer::FamosProducer ( edm::ParameterSet const & p )

explicit

Definition at line 32 of file FamosProducer.cc.

References famosManager_, edm::ParameterSet::getParameter(), simulateMuons, and sourceToken.

 {    
 
     produces<edm::SimTrackContainer>();
     produces<edm::SimVertexContainer>();
     produces<FSimVertexTypeCollection>("VertexTypes");
     produces<edm::PSimHitContainer>("TrackerHits");
     produces<edm::PCaloHitContainer>("EcalHitsEB");
     produces<edm::PCaloHitContainer>("EcalHitsEE");
     produces<edm::PCaloHitContainer>("EcalHitsES");
     produces<edm::PCaloHitContainer>("HcalHits");
     // Temporary facility to allow for the crossing frame to work...
     simulateMuons = p.getParameter<bool>("SimulateMuons");
     if ( simulateMuons ) produces<edm::SimTrackContainer>("MuonSimTracks");
 
     // hepmc event from signal event
     edm::InputTag sourceLabel = p.getParameter<edm::InputTag>("SourceLabel");
     sourceToken = consumes<edm::HepMCProduct>(sourceLabel);
     
     // famos manager
     famosManager_ = new FamosManager(p);
 }

FamosProducer::~FamosProducer ( )

override

Definition at line 55 of file FamosProducer.cc.

References famosManager_.

56 { delete famosManager_; }

FamosProducer::famosManager_

FamosManager * famosManager_

Definition: FamosProducer.h:29

Member Function Documentation

void FamosProducer::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	es
	)

override

Definition at line 59 of file FamosProducer.cc.

References famosManager_, and FamosManager::setupGeometryAndField().

                                                                    {
   famosManager_->setupGeometryAndField(run,es);
 }

void FamosProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

override

Definition at line 63 of file FamosProducer.cc.

References FamosManager::calorimetryManager(), DEFINE_FWK_MODULE, famosManager_, edm::EventSetup::get(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), FSimEvent::load(), CalorimetryManager::loadFromEcalBarrel(), CalorimetryManager::loadFromEcalEndcap(), CalorimetryManager::loadFromHcal(), CalorimetryManager::loadFromPreshower(), CalorimetryManager::loadMuonSimTracks(), TrajectoryManager::loadSimHits(), eostools::move(), p1, p2, p3, p4, edm::Event::put(), random, FamosManager::reconstruct(), FamosManager::simEvent(), simulateMuons, sourceToken, edm::Event::streamID(), FBaseSimEvent::theTable(), trackingTruthProducer_cfi::tracker, and FamosManager::trackerManager().

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

 {
   ParticleTable::Sentry ptable(famosManager_->simEvent()->theTable());
    using namespace edm;
 
    RandomEngineAndDistribution random(iEvent.streamID());
 
    //Retrieve tracker topology from geometry
    edm::ESHandle<TrackerTopology> tTopoHand;
    es.get<TrackerTopologyRcd>().get(tTopoHand);
    const TrackerTopology *tTopo=tTopoHand.product();
 
    // get the signal event
    Handle<HepMCProduct> theHepMCProduct;
    iEvent.getByToken(sourceToken,theHepMCProduct);
    const HepMC::GenEvent * myGenEvent = theHepMCProduct->GetEvent();
    
    // do the simulation
    famosManager_->reconstruct(myGenEvent,tTopo, &random);
 
    // get the hits, simtracks and simvertices and put in the event
    CalorimetryManager * calo = famosManager_->calorimetryManager();
    TrajectoryManager * tracker = famosManager_->trackerManager();
    
    std::unique_ptr<edm::SimTrackContainer> p1(new edm::SimTrackContainer);
    std::unique_ptr<edm::SimTrackContainer> m1(new edm::SimTrackContainer);
    std::unique_ptr<edm::SimVertexContainer> p2(new edm::SimVertexContainer);
    std::unique_ptr<FSimVertexTypeCollection> v1(new FSimVertexTypeCollection);
    std::unique_ptr<edm::PSimHitContainer> p3(new edm::PSimHitContainer);
    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);
    
    FSimEvent* fevt = famosManager_->simEvent();
    fevt->load(*p1,*m1);
    fevt->load(*p2);
    fevt->load(*v1);
    tracker->loadSimHits(*p3);
 
 
    if ( calo ) {  
      calo->loadFromEcalBarrel(*p4);
      calo->loadFromEcalEndcap(*p5);
      calo->loadFromPreshower(*p6);
      calo->loadFromHcal(*p7);
      calo->loadMuonSimTracks(*m1);
    }
 
    if ( simulateMuons ) iEvent.put(std::move(m1),"MuonSimTracks");
    iEvent.put(std::move(p1));
    iEvent.put(std::move(p2));
    iEvent.put(std::move(p3),"TrackerHits");
    iEvent.put(std::move(v1),"VertexTypes");
    iEvent.put(std::move(p4),"EcalHitsEB");
    iEvent.put(std::move(p5),"EcalHitsEE");
    iEvent.put(std::move(p6),"EcalHitsES");
    iEvent.put(std::move(p7),"HcalHits");
 }