#include <yetkin/GenHIEventProducer/src/GenHIEventProducer.cc>

Inheritance diagram for GenHIEventProducer:

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

	~GenHIEventProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	doParticleInfo_

edm::EDGetTokenT< CrossingFrame< edm::HepMCProduct > >	hepmcSrc_

edm::ESHandle< ParticleDataTable >	pdt

double	ptCut_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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::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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 49 of file GenHIEventProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 75 of file GenHIEventProducer.cc.

References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

 {
     produces<edm::GenHIEvent>();
     hepmcSrc_ = consumes<CrossingFrame<edm::HepMCProduct> >(iConfig.getParameter<edm::InputTag>("src"));
     doParticleInfo_ = iConfig.getUntrackedParameter<bool>("doParticleInfo",false);
     if(doParticleInfo_){
       ptCut_ = iConfig.getParameter<double> ("ptCut");
     }
 }

GenHIEventProducer::~GenHIEventProducer ( )

Definition at line 86 of file GenHIEventProducer.cc.

 {
 
     // do anything here that needs to be done at desctruction time
     // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

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

overrideprivatevirtual

Definition at line 101 of file GenHIEventProducer.cc.

References b, begin, ALCARECOTkAlJpsiMuMu_cff::charge, DEFINE_FWK_MODULE, end, stringResolutionProvider_cfi::eta, Exception, edm::Event::getByToken(), edm::EventSetup::getData(), edm::HepMCProduct::GetEvent(), ecalTB2006H4_GenSimDigiReco_cfg::mix, eostools::move(), hpstanc_transforms::nCharged, np, HiggsValidation_cfi::pdg_id, edm::Handle< T >::product(), EnergyCorrector::pt, and edm::Event::put().

 {
     using namespace edm;
 
     if(!(pdt.isValid())) iSetup.getData(pdt);
 
     double b = -1;
     int npart = -1;
     int ncoll = 0;
     int nhard = 0;
     double phi = 0;
     double ecc = -1;
 
     int nCharged = 0;
     int nChargedMR = 0;
     int nChargedPtCut = 0; // NchargedPtCut bym
     int nChargedPtCutMR = 0; // NchargedPtCutMR bym
 
     double meanPt = 0;
     double meanPtMR = 0;
     double EtMR = 0; // Normalized of total energy bym
     double TotEnergy = 0; // Total energy bym
 
     Handle<CrossingFrame<edm::HepMCProduct> > hepmc;
     iEvent.getByToken(hepmcSrc_,hepmc);
     MixCollection<HepMCProduct> mix(hepmc.product());
     
     if(mix.size() < 1){
       throw cms::Exception("MatchVtx")
     <<"Mixing has "<<mix.size()<<" sub-events, should have been at least 1"
     <<endl;
     }
 
     const HepMCProduct& hievt = mix.getObject(mix.size()-1);
     const HepMC::GenEvent* evt = hievt.GetEvent();
     if(doParticleInfo_){
       HepMC::GenEvent::particle_const_iterator begin = evt->particles_begin();
       HepMC::GenEvent::particle_const_iterator end = evt->particles_end();
       for(HepMC::GenEvent::particle_const_iterator it = begin; it != end; ++it){
         if((*it)->status() != 1) continue;
         int pdg_id = (*it)->pdg_id();
             const ParticleData * part = pdt->particle(pdg_id );
             int charge = static_cast<int>(part->charge());
         
             if(charge == 0) continue;
             float pt = (*it)->momentum().perp();
             float eta = (*it)->momentum().eta();
             float energy = (*it)->momentum().e(); // energy bym
             //float energy = (*it)->momentum().energy(); // energy bym
             nCharged++;
             meanPt += pt;
             // Get the total energy bym
             if(fabs(eta)<1.0){
           TotEnergy += energy;
             }
             if(pt>ptCut_){
           nChargedPtCut++;
           if(fabs(eta)<0.5){
           nChargedPtCutMR++;
                 }
             }
             // end bym
         
             if(fabs(eta) > 0.5) continue;
             nChargedMR++;
             meanPtMR += pt;
       }
     }
         const HepMC::HeavyIon* hi = evt->heavy_ion();
 
         if(hi){
             ncoll = ncoll + hi->Ncoll();
             nhard = nhard + hi->Ncoll_hard();
             int np = hi->Npart_proj() + hi->Npart_targ();
             if(np >= 0){
                 npart = np;
                 b = hi->impact_parameter();
                 phi = hi->event_plane_angle();
         ecc = hi->eccentricity();
             }
         }
 
     // Get the normalized total energy bym
     if(TotEnergy != 0){
         EtMR = TotEnergy/2;
     }
 
     if(nChargedMR != 0){
         meanPtMR /= nChargedMR;
     }
     if(nCharged != 0){
         meanPt /= nCharged;
     }
 
     std::unique_ptr<edm::GenHIEvent> pGenHI(new edm::GenHIEvent(b,
                                   npart,
                                   ncoll,
                                   nhard,
                                   phi, 
                                   ecc,
                                   nCharged,
                                   nChargedMR,
                                   meanPt,
                                   meanPtMR,   
                                   EtMR, 
                                   nChargedPtCut,
                                   nChargedPtCutMR
                                   ));
 
     iEvent.put(std::move(pGenHI));
 
 }