#include <HiBasicGenTest.h>

Inheritance diagram for HiBasicGenTest:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override

	HiBasicGenTest (const edm::ParameterSet &)

	~HiBasicGenTest () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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 ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	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 &	itemsToGetFrom (BranchType iType) 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
MonitorElement *	b [3]

MonitorElement *	dnchdeta [3]

MonitorElement *	dnchdphi [3]

MonitorElement *	dnchdpt [3]

edm::EDGetTokenT< edm::HepMCProduct >	generatorToken_

edm::ESHandle< ParticleDataTable >	pdt

MonitorElement *	rp

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 21 of file HiBasicGenTest.h.

Constructor & Destructor Documentation

HiBasicGenTest::HiBasicGenTest ( const edm::ParameterSet & iPSet )

explicit

Definition at line 17 of file HiBasicGenTest.cc.

References edm::ParameterSet::getParameter().

 {
   generatorToken_ = consumes<edm::HepMCProduct> (
       iPSet.getParameter<edm::InputTag>("generatorLabel"));
 }

HiBasicGenTest::~HiBasicGenTest ( )

override

Definition at line 23 of file HiBasicGenTest.cc.

23 {}

Member Function Documentation

void HiBasicGenTest::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 59 of file HiBasicGenTest.cc.

References b, begin, ALCARECOTkAlJpsiMuMu_cff::charge, DEFINE_FWK_MODULE, end, PVValHelper::eta, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), CaloTowersParam_cfi::mc, AlCaHLTBitMon_ParallelJobs::p, HiggsValidation_cfi::pdg_id, Pi, pi, and EnergyCorrector::pt.

 {
 
   Handle<HepMCProduct> mc;
   iEvent.getByToken(generatorToken_, mc);
   const HepMC::GenEvent *evt = mc->GetEvent();
   const HepMC::HeavyIon *hi = evt->heavy_ion();
 
   int cbin = 0;
   double phi0 =0.;
 
   if(hi){
 
     double ip = hi->impact_parameter();
     phi0 = hi->event_plane_angle();
     
     // fill reaction plane distribution
     rp->Fill(phi0);
     
     // if the event is in one of the centrality bins of interest fill hists
     int cbin=-1;
     if(ip < 5.045) cbin=0;
     else if (ip < 7.145 && ip > 5.045) cbin=1;
     else if (ip < 15.202 && ip > 14.283) cbin=2;
     if(cbin<0) return;
     
     // fill impact parameter distributions
     b[cbin]->Fill(ip);
   }
 
   // loop over particles
   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){
 
     // only fill hists for status=1 particles
     if((*it)->status() != 1) continue;
 
     // only fill hists for charged particles
     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 eta = (*it)->momentum().eta();
     float phi = (*it)->momentum().phi();
     float pt = (*it)->momentum().perp();
 
     dnchdeta[cbin]->Fill(eta);
     dnchdpt[cbin]->Fill(pt);
 
     double pi = TMath::Pi();
     double p = phi-phi0;
     if(p > pi) p = p - 2*pi;
     if(p < -1*pi) p = p + 2*pi;
     dnchdphi[cbin]->Fill(p);
 
   }
 
   return;
 
 }

void HiBasicGenTest::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Setting the DQM top directories

Booking the ME's

Implements DQMEDAnalyzer.

Definition at line 25 of file HiBasicGenTest.cc.

References b, DQMStore::IBooker::book1D(), DQMStore::IBooker::setCurrentFolder(), and DQMStore::IBooker::tag().

                                          {
 
   ibooker.setCurrentFolder("Generator/Particles");
 
   for (int ibin = 0; ibin < 3; ibin++) {
     dnchdeta[ibin] = ibooker.book1D(Form("dnchdeta%d", ibin), ";#eta;dN^{ch}/d#eta",
         100, -6.0, 6.0);
 
     dnchdpt[ibin] = ibooker.book1D(Form("dnchdpt%d", ibin), ";p_{T};dN^{ch}/dp_{T}",
         200, 0.0, 100.0);
 
     b[ibin] = ibooker.book1D(Form("b%d",ibin),";b[fm];events", 100, 0.0, 20.0);
     dnchdphi[ibin] = ibooker.book1D(Form("dnchdphi%d", ibin), ";#phi;dN^{ch}/d#phi",
         100, -3.2, 3.2);
 
     ibooker.tag(dnchdeta[ibin], 1+ibin*4);
     ibooker.tag(dnchdpt[ibin], 2+ibin*4);
     ibooker.tag(b[ibin], 3+ibin*4);
     ibooker.tag(dnchdphi[ibin], 4+ibin*4);
   }
 
   rp = ibooker.book1D("phi0", ";#phi_{RP};events", 100, -3.2, 3.2);
   ibooker.tag(rp, 13);
 }

void HiBasicGenTest::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

override

Definition at line 53 of file HiBasicGenTest.cc.

References edm::EventSetup::getData().

 {
   iSetup.getData(pdt);
 }

Member Data Documentation

MonitorElement* HiBasicGenTest::b[3]

private

Definition at line 37 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdeta[3]

private

Definition at line 35 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdphi[3]

private

Definition at line 38 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdpt[3]

private

Definition at line 36 of file HiBasicGenTest.h.

edm::EDGetTokenT<edm::HepMCProduct> HiBasicGenTest::generatorToken_

private

Definition at line 34 of file HiBasicGenTest.h.

edm::ESHandle< ParticleDataTable > HiBasicGenTest::pdt

private

Definition at line 41 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::rp

private

Definition at line 39 of file HiBasicGenTest.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation