FamosProducer Class Reference

#include <FamosProducer.h>

Inheritance diagram for FamosProducer:

Public Member Functions
virtual void	beginRun (edm::Run const &run, const edm::EventSetup &es) override
virtual void	endJob () override
	FamosProducer (edm::ParameterSet const &p)
virtual void	produce (edm::Event &e, const edm::EventSetup &c) override
virtual	~FamosProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
edm::InputTag	beamSpotLabel
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken
HepMC::GenEvent *	evt_
FamosManager *	famosManager_
edm::InputTag	genParticleLabel
edm::EDGetTokenT < reco::GenParticleCollection >	genParticleToken
edm::EDGetTokenT < reco::GenParticleCollection >	mixGenParticleToken
edm::EDGetTokenT < CrossingFrame < edm::HepMCProduct > >	mixSourceToken
edm::EDGetTokenT < edm::HepMCProduct >	puToken
bool	simulateMuons
edm::InputTag	sourceLabel
edm::EDGetTokenT < edm::HepMCProduct >	sourceToken

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 20 of file FamosProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 32 of file FamosProducer.cc.

References beamSpotLabel, beamSpotToken, famosManager_, genParticleLabel, genParticleToken, edm::ParameterSet::getParameter(), mixGenParticleToken, mixSourceToken, puToken, simulateMuons, sourceLabel, 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");

    // The generator input label
    sourceLabel = p.getParameter<edm::InputTag>("SourceLabel");
    genParticleLabel = p.getParameter<edm::InputTag>("GenParticleLabel");
    beamSpotLabel = p.getParameter<edm::InputTag>("BeamSpotLabel");

    // consume declarations
    beamSpotToken = consumes<reco::BeamSpot>(beamSpotLabel);
    genParticleToken = consumes<reco::GenParticleCollection>(genParticleLabel);
    // FUTURE OBSOLETE CODE
    sourceToken = consumes<edm::HepMCProduct>(sourceLabel);
    edm::InputTag _label = edm::InputTag("famosPileUp","PileUpEvents");
    puToken = consumes<edm::HepMCProduct>(_label);
    // OBSOLETE CODE
    _label = edm::InputTag("mixGenPU","generator");
    mixSourceToken = consumes<CrossingFrame<edm::HepMCProduct> >(_label);
    _label = edm::InputTag("genParticlesFromMixingModule");
    mixGenParticleToken = consumes<reco::GenParticleCollection>(_label);
    
    // famos manager
    famosManager_ = new FamosManager(p);
}

FamosProducer::~FamosProducer ( )

virtual

Definition at line 69 of file FamosProducer.cc.

References famosManager_.

{ if ( famosManager_ ) delete famosManager_; }

Member Function Documentation

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

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 73 of file FamosProducer.cc.

References famosManager_, and FamosManager::setupGeometryAndField().

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

void FamosProducer::endJob ( void  )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 77 of file FamosProducer.cc.

{ 
}

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

overridevirtual

Implements edm::EDProducer.

Definition at line 81 of file FamosProducer.cc.

References beamSpotToken, FamosManager::calorimetryManager(), gather_cfg::cout, famosManager_, TtGenEvtProducer_cfi::genEvt, genParticleToken, edm::EventSetup::get(), edm::Event::getByToken(), mixGenParticleToken, mixSourceToken, p1, p2, p3, p4, edm::Event::put(), puToken, random, FamosManager::reconstruct(), FamosManager::simEvent(), simulateMuons, source, sourceToken, edm::Event::streamID(), FBaseSimEvent::theTable(), patCandidatesForDimuonsSequences_cff::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());

   //  // The beam spot position
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   iEvent.getByToken(beamSpotToken,recoBeamSpotHandle); 
   math::XYZPoint BSPosition_ = recoBeamSpotHandle->position();

   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHand;
   es.get<IdealGeometryRecord>().get(tTopoHand);
   const TrackerTopology *tTopo=tTopoHand.product();

      
   const HepMC::GenEvent* myGenEvent = 0;
   FSimEvent* fevt = famosManager_->simEvent();
   //   fevt->setBeamSpot(BSPosition_);     
   
   // Get the generated event(s) from the edm::Event
   // 1. Check if a HepMCProduct exists
   //    a. Take the VtxSmeared if it exists
   //    b. Take the source  otherwise
   // 2. Otherwise go for the CandidateCollection
   
   Handle<HepMCProduct> theHepMCProduct;
   
   const reco::GenParticleCollection* myGenParticlesXF = 0; //OBSOLETE
   const reco::GenParticleCollection* myGenParticles = 0;
   const HepMC::GenEvent* thePUEvents = 0;

   // BEGIN OBSOLETE CODE
   Handle<CrossingFrame<HepMCProduct> > theHepMCProductCrossingFrame;
   bool isPileUpXF = iEvent.getByToken(mixSourceToken,theHepMCProductCrossingFrame);
   if (isPileUpXF){// take the GenParticle from crossingframe event collection, if it exists 
     Handle<reco::GenParticleCollection> genEvtXF;
     bool genPartXF = iEvent.getByToken(mixGenParticleToken,genEvtXF);
     if(genPartXF) myGenParticlesXF = &(*genEvtXF);
   }
   else{// otherwise, use the old famos PU     
   // END OBSOLETE CODE
     // Get the generated signal event
     // BEGIN FUTURE OBSOLETE CODE
     bool source = iEvent.getByToken(sourceToken,theHepMCProduct);
     if ( source ) { 
       myGenEvent = theHepMCProduct->GetEvent();
     } 

     fevt->setBeamSpot(BSPosition_);
     // GEN LEVEL INFO NOT IN HEPMC FORMAT
     //     In case there is no HepMCProduct, seek a genParticle Candidate Collection
     bool genPart = false;
     if ( !myGenEvent ) { 
       //END FUTURE OBSOLETE CODE
       // Look for the particle CandidateCollection
       Handle<reco::GenParticleCollection> genEvt;
       genPart = iEvent.getByToken(genParticleToken,genEvt);
       if ( genPart ) myGenParticles = &(*genEvt);
     }
       
     if ( !myGenEvent && !genPart )
       std::cout << "There is no generator input for this event, under " 
         << "any form (HepMCProduct, genParticles)" << std::endl
         << "Please check SourceLabel or GenParticleLabel" << std::endl;
       
     // BEGIN FUTURE OBSOLETE CODE
     // Get the pile-up events from the pile-up producer
     // There might be no pile-up events, by the way, in that case, just continue
     Handle<HepMCProduct> thePileUpEvents;
     bool isPileUp = iEvent.getByToken(puToken,thePileUpEvents);
     thePUEvents = isPileUp ? thePileUpEvents->GetEvent() : 0;
     // END FUTURE OBSOLETE CODE
   }//end else
   

   // pass the event to the Famos Manager for propagation and simulation
   if (myGenParticlesXF) {// OBSOLETE OPTION
     famosManager_->reconstruct(myGenParticlesXF,tTopo, &random);
   } else {
     famosManager_->reconstruct(myGenEvent,myGenParticles,thePUEvents,tTopo, &random); // FUTURE OBSOLETE ARGUMENTS: myGenEvents, thePUEvents 
   }

   CalorimetryManager * calo = famosManager_->calorimetryManager();
   TrajectoryManager * tracker = famosManager_->trackerManager();
   
   // Save everything in the edm::Event
   std::auto_ptr<edm::SimTrackContainer> p1(new edm::SimTrackContainer);
   std::auto_ptr<edm::SimTrackContainer> m1(new edm::SimTrackContainer);
   std::auto_ptr<edm::SimVertexContainer> p2(new edm::SimVertexContainer);
   std::auto_ptr<FSimVertexTypeCollection> v1(new FSimVertexTypeCollection);
   std::auto_ptr<edm::PSimHitContainer> p3(new edm::PSimHitContainer);
   std::auto_ptr<edm::PCaloHitContainer> p4(new edm::PCaloHitContainer);
   std::auto_ptr<edm::PCaloHitContainer> p5(new edm::PCaloHitContainer);
   std::auto_ptr<edm::PCaloHitContainer> p6(new edm::PCaloHitContainer); 
   std::auto_ptr<edm::PCaloHitContainer> p7(new edm::PCaloHitContainer);
   
   fevt->load(*p1,*m1);
   fevt->load(*p2);
   fevt->load(*v1);
   //   fevt->print();
   tracker->loadSimHits(*p3);

   //   fevt->print();

   if ( calo ) {  
     calo->loadFromEcalBarrel(*p4);
     calo->loadFromEcalEndcap(*p5);
     calo->loadFromPreshower(*p6);
     calo->loadFromHcal(*p7);
     // update the muon SimTracks
     calo->loadMuonSimTracks(*m1);
   }

   // Write muon first, to allow tracking particles to work... (pending MixingModule fix)
   if ( simulateMuons ) iEvent.put(m1,"MuonSimTracks");
   iEvent.put(p1);
   iEvent.put(p2);
   iEvent.put(p3,"TrackerHits");
   iEvent.put(v1,"VertexTypes");
   iEvent.put(p4,"EcalHitsEB");
   iEvent.put(p5,"EcalHitsEE");
   iEvent.put(p6,"EcalHitsES");
   iEvent.put(p7,"HcalHits");

}

Member Data Documentation

edm::InputTag FamosProducer::beamSpotLabel

private

Definition at line 40 of file FamosProducer.h.

Referenced by FamosProducer().

edm::EDGetTokenT<reco::BeamSpot> FamosProducer::beamSpotToken

private

Definition at line 43 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().

HepMC::GenEvent* FamosProducer::evt_

private

Definition at line 34 of file FamosProducer.h.

FamosManager* FamosProducer::famosManager_

private

Definition at line 33 of file FamosProducer.h.

Referenced by beginRun(), FamosProducer(), produce(), and ~FamosProducer().

edm::InputTag FamosProducer::genParticleLabel

private

Definition at line 39 of file FamosProducer.h.

Referenced by FamosProducer().

edm::EDGetTokenT<reco::GenParticleCollection> FamosProducer::genParticleToken

private

Definition at line 44 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().

edm::EDGetTokenT<reco::GenParticleCollection> FamosProducer::mixGenParticleToken

private

Definition at line 50 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().

edm::EDGetTokenT<CrossingFrame<edm::HepMCProduct> > FamosProducer::mixSourceToken

private

Definition at line 49 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().

edm::EDGetTokenT<edm::HepMCProduct> FamosProducer::puToken

private

Definition at line 47 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().

bool FamosProducer::simulateMuons

private

Definition at line 35 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().

edm::InputTag FamosProducer::sourceLabel

private

Definition at line 38 of file FamosProducer.h.

Referenced by FamosProducer().

edm::EDGetTokenT<edm::HepMCProduct> FamosProducer::sourceToken

private

Definition at line 46 of file FamosProducer.h.

Referenced by FamosProducer(), and produce().