HiBasicGenTest Class Reference

#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
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
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.

{}

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, end, eta(), edm::Event::getByToken(), AlCaHLTBitMon_ParallelJobs::p, phi, 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::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.

{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.