#include <HiBasicGenTest.h>

Inheritance diagram for HiBasicGenTest:

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

virtual void	beginJob ()

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

virtual void	endJob ()

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

	HiBasicGenTest (const edm::ParameterSet &)

virtual	~HiBasicGenTest ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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
MonitorElement *	b [3]

DQMStore *	dbe

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 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(), and cppFunctionSkipper::operator.

 {
   dbe = 0;
   dbe = edm::Service<DQMStore>().operator->();
   generatorToken_ = consumes<edm::HepMCProduct> (
       iPSet.getParameter<edm::InputTag>("generatorLabel"));
 }

HiBasicGenTest::~HiBasicGenTest ( )

virtual

Definition at line 25 of file HiBasicGenTest.cc.

25 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 67 of file HiBasicGenTest.cc.

References b, begin, DeDxDiscriminatorTools::charge(), end, eta(), edm::Event::getByToken(), AlCaHLTBitMon_ParallelJobs::p, phi, Pi, pi, and RecoTauCleanerPlugins::pt.

 {
 
   Handle<HepMCProduct> mc;
   iEvent.getByToken(generatorToken_, mc);
   const HepMC::GenEvent *evt = mc->GetEvent();
   const HepMC::HeavyIon *hi = evt->heavy_ion();
 
   double ip = hi->impact_parameter();
   double 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::beginJob ( void )

virtual

Setting the DQM top directories

Booking the ME's

Reimplemented from edm::EDAnalyzer.

Definition at line 27 of file HiBasicGenTest.cc.

References b.

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

void HiBasicGenTest::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 59 of file HiBasicGenTest.cc.

References edm::EventSetup::getData().

 {
   iSetup.getData(pdt);
   return;
 }

void HiBasicGenTest::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 54 of file HiBasicGenTest.cc.

                            {
   // normalization of histograms can be done here (or in post-processor)
   return;
 }

void HiBasicGenTest::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 65 of file HiBasicGenTest.cc.

65 {return;}

Member Data Documentation

MonitorElement* HiBasicGenTest::b[3]

private

Definition at line 39 of file HiBasicGenTest.h.

DQMStore* HiBasicGenTest::dbe

private

Definition at line 34 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdeta[3]

private

Definition at line 37 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdphi[3]

private

Definition at line 40 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdpt[3]

private

Definition at line 38 of file HiBasicGenTest.h.

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

private

Definition at line 36 of file HiBasicGenTest.h.

edm::ESHandle< ParticleDataTable > HiBasicGenTest::pdt

private

Definition at line 43 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::rp

private

Definition at line 41 of file HiBasicGenTest.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation