#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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
HepMC::GenEvent *	evt_

FamosManager *	famosManager_

bool	simulateMuons

edm::InputTag	theBeamSpotLabel

edm::InputTag	theGenParticleLabel

edm::InputTag	theSourceLabel

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 15 of file FamosProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file FamosProducer.cc.

References famosManager_, edm::ParameterSet::getParameter(), simulateMuons, theBeamSpotLabel, theGenParticleLabel, and theSourceLabel.

 {    
 
     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
     theSourceLabel = p.getParameter<edm::InputTag>("SourceLabel");
     theGenParticleLabel = p.getParameter<edm::InputTag>("GenParticleLabel");
     theBeamSpotLabel = p.getParameter<edm::InputTag>("BeamSpotLabel");
 
     famosManager_ = new FamosManager(p);
   
 }

FamosProducer::~FamosProducer ( )

virtual

Definition at line 60 of file FamosProducer.cc.

References famosManager_.

61 { if ( famosManager_ ) delete famosManager_; }

FamosProducer::famosManager_

FamosManager * famosManager_

Definition: FamosProducer.h:28

Member Function Documentation

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

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 64 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 68 of file FamosProducer.cc.

69 {

70 }

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

overridevirtual

Implements edm::EDProducer.

Definition at line 72 of file FamosProducer.cc.

References FamosManager::calorimetryManager(), gather_cfg::cout, famosManager_, TtGenEvtProducer_cfi::genEvt, edm::EventSetup::get(), edm::Event::getByLabel(), p1, p2, p3, p4, edm::Event::put(), random, FamosManager::reconstruct(), FamosManager::simEvent(), simulateMuons, source, edm::Event::streamID(), theBeamSpotLabel, theGenParticleLabel, theSourceLabel, patCandidatesForDimuonsSequences_cff::tracker, and FamosManager::trackerManager().

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

 {
    using namespace edm;
 
    RandomEngineAndDistribution random(iEvent.streamID());
 
    //  // The beam spot position
    edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
    iEvent.getByLabel(theBeamSpotLabel,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;
    
    PrimaryVertexGenerator* theVertexGenerator = fevt->thePrimaryVertexGenerator();
    
    
    const reco::GenParticleCollection* myGenParticlesXF = 0;
    const reco::GenParticleCollection* myGenParticles = 0;
    const HepMC::GenEvent* thePUEvents = 0;
 
    Handle<CrossingFrame<HepMCProduct> > theHepMCProductCrossingFrame;
    bool isPileUpXF = iEvent.getByLabel("mixGenPU","generator",theHepMCProductCrossingFrame);
 
    if (isPileUpXF){// take the GenParticle from crossingframe event collection, if it exists 
      Handle<reco::GenParticleCollection> genEvtXF;
      bool genPartXF = iEvent.getByLabel("genParticlesFromMixingModule",genEvtXF);
      if(genPartXF) myGenParticlesXF = &(*genEvtXF);
    }
    else{// otherwise, use the old famos PU     
      // Get the generated signal event
      bool source = iEvent.getByLabel(theSourceLabel,theHepMCProduct);
      if ( source ) { 
        myGenEvent = theHepMCProduct->GetEvent();
        // First rotate in case of beam crossing angle (except if done already)
        if ( theVertexGenerator ) { 
      TMatrixD* boost = theVertexGenerator->boost();
      if ( boost ) theHepMCProduct->boostToLab(boost,"momentum");
        }          
        myGenEvent = theHepMCProduct->GetEvent();
      } 
 
      fevt->setBeamSpot(BSPosition_);
        
      //     In case there is no HepMCProduct, seek a genParticle Candidate Collection
      bool genPart = false;
      if ( !myGenEvent ) { 
        // Look for the particle CandidateCollection
        Handle<reco::GenParticleCollection> genEvt;
        genPart = iEvent.getByLabel(theGenParticleLabel,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;
        
      // 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.getByLabel("famosPileUp","PileUpEvents",thePileUpEvents);
      thePUEvents = isPileUp ? thePileUpEvents->GetEvent() : 0;
       
    }//end else
 
    // pass the event to the Famos Manager for propagation and simulation
    if (myGenParticlesXF) {
      famosManager_->reconstruct(myGenParticlesXF,tTopo, &random);
    } else {
      famosManager_->reconstruct(myGenEvent,myGenParticles,thePUEvents,tTopo, &random);
    }
 
    // Set the vertex back to the HepMCProduct (except if it was smeared already)
    if ( myGenEvent ) { 
      if ( theVertexGenerator ) { 
        HepMC::FourVector theVertex(
                    (theVertexGenerator->X()-theVertexGenerator->beamSpot().X()+BSPosition_.X())*10.,
                    (theVertexGenerator->Y()-theVertexGenerator->beamSpot().Y()+BSPosition_.Y())*10.,
                    (theVertexGenerator->Z()-theVertexGenerator->beamSpot().Z()+BSPosition_.Z())*10.,
                    0.);
        if ( fabs(theVertexGenerator->Z()) > 1E-10 ) theHepMCProduct->applyVtxGen( &theVertex );
      }
    }
    
    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");
 
 }