CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Attributes
HiBasicGenTest Class Reference

#include <HiBasicGenTest.h>

Inheritance diagram for HiBasicGenTest:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

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 &event, edm::EventSetup const &setup) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
void beginStream (edm::StreamID id) final
 
 DQMEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProduceroperator= (const EDProducer &)=delete
 

Private Attributes

MonitorElementb [3]
 
MonitorElementdnchdeta [3]
 
MonitorElementdnchdphi [3]
 
MonitorElementdnchdpt [3]
 
edm::EDGetTokenT
< edm::HepMCProduct
generatorToken_
 
edm::ESHandle< ParticleDataTablepdt
 
edm::ESGetToken
< ParticleDataTable,
edm::DefaultRecord
pdtToken_
 
MonitorElementrp
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using CacheTypes = CacheContexts< T...>
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T...>
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static void globalEndJob (DQMEDAnalyzerGlobalCache const *)
 
static void globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
 
static void globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
 
static std::unique_ptr
< DQMEDAnalyzerGlobalCache
initializeGlobalCache (edm::ParameterSet const &)
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMTokenlumiToken_
 
edm::EDPutTokenT< DQMTokenrunToken_
 
unsigned int streamId_
 

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

17  {
18  generatorToken_ = consumes<edm::HepMCProduct>(iPSet.getParameter<edm::InputTag>("generatorLabel"));
19  pdtToken_ = esConsumes<edm::Transition::BeginRun>();
20 }
edm::EDGetTokenT< edm::HepMCProduct > generatorToken_
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::ESGetToken< ParticleDataTable, edm::DefaultRecord > pdtToken_
HiBasicGenTest::~HiBasicGenTest ( )
override

Definition at line 22 of file HiBasicGenTest.cc.

22 {}

Member Function Documentation

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

Reimplemented from DQMEDAnalyzer.

Definition at line 45 of file HiBasicGenTest.cc.

References b, SplitLinear::begin, RecoTauCleanerPlugins::charge, dataset::end, PVValHelper::eta, edm::Event::getByToken(), AlCaHLTBitMon_ParallelJobs::p, Pi, pi, and DiDispStaMuonMonitor_cfi::pt.

45  {
47  iEvent.getByToken(generatorToken_, mc);
48  const HepMC::GenEvent* evt = mc->GetEvent();
49  const HepMC::HeavyIon* hi = evt->heavy_ion();
50 
51  int cbin = 0;
52  double phi0 = 0.;
53 
54  if (hi) {
55  double ip = hi->impact_parameter();
56  phi0 = hi->event_plane_angle();
57 
58  // fill reaction plane distribution
59  rp->Fill(phi0);
60 
61  // if the event is in one of the centrality bins of interest fill hists
62  int cbin = -1;
63  if (ip < 5.045)
64  cbin = 0;
65  else if (ip < 7.145 && ip > 5.045)
66  cbin = 1;
67  else if (ip < 15.202 && ip > 14.283)
68  cbin = 2;
69  if (cbin < 0)
70  return;
71 
72  // fill impact parameter distributions
73  b[cbin]->Fill(ip);
74  }
75 
76  // loop over particles
77  HepMC::GenEvent::particle_const_iterator begin = evt->particles_begin();
78  HepMC::GenEvent::particle_const_iterator end = evt->particles_end();
79  for (HepMC::GenEvent::particle_const_iterator it = begin; it != end; ++it) {
80  // only fill hists for status=1 particles
81  if ((*it)->status() != 1)
82  continue;
83 
84  // only fill hists for charged particles
85  int pdg_id = (*it)->pdg_id();
86  const ParticleData* part = pdt->particle(pdg_id);
87  int charge = static_cast<int>(part->charge());
88  if (charge == 0)
89  continue;
90 
91  float eta = (*it)->momentum().eta();
92  float phi = (*it)->momentum().phi();
93  float pt = (*it)->momentum().perp();
94 
95  dnchdeta[cbin]->Fill(eta);
96  dnchdpt[cbin]->Fill(pt);
97 
98  double pi = TMath::Pi();
99  double p = phi - phi0;
100  if (p > pi)
101  p = p - 2 * pi;
102  if (p < -1 * pi)
103  p = p + 2 * pi;
104  dnchdphi[cbin]->Fill(p);
105  }
106 
107  return;
108 }
const double Pi
MonitorElement * rp
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
edm::EDGetTokenT< edm::HepMCProduct > generatorToken_
MonitorElement * dnchdpt[3]
const Double_t pi
void Fill(long long x)
MonitorElement * dnchdeta[3]
MonitorElement * dnchdphi[3]
MonitorElement * b[3]
HepPDT::ParticleData ParticleData
part
Definition: HCALResponse.h:20
string end
Definition: dataset.py:937
edm::ESHandle< ParticleDataTable > pdt
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 24 of file HiBasicGenTest.cc.

References b, dqm::implementation::IBooker::book1D(), and dqm::implementation::NavigatorBase::setCurrentFolder().

24  {
26  ibooker.setCurrentFolder("Generator/Particles");
27 
29  for (int ibin = 0; ibin < 3; ibin++) {
30  dnchdeta[ibin] = ibooker.book1D(Form("dnchdeta%d", ibin), ";#eta;dN^{ch}/d#eta", 100, -6.0, 6.0);
31 
32  dnchdpt[ibin] = ibooker.book1D(Form("dnchdpt%d", ibin), ";p_{T};dN^{ch}/dp_{T}", 200, 0.0, 100.0);
33 
34  b[ibin] = ibooker.book1D(Form("b%d", ibin), ";b[fm];events", 100, 0.0, 20.0);
35  dnchdphi[ibin] = ibooker.book1D(Form("dnchdphi%d", ibin), ";#phi;dN^{ch}/d#phi", 100, -3.2, 3.2);
36  }
37 
38  rp = ibooker.book1D("phi0", ";#phi_{RP};events", 100, -3.2, 3.2);
39 }
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
MonitorElement * rp
MonitorElement * dnchdpt[3]
MonitorElement * dnchdeta[3]
MonitorElement * dnchdphi[3]
MonitorElement * b[3]
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
void HiBasicGenTest::dqmBeginRun ( const edm::Run r,
const edm::EventSetup c 
)
overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 41 of file HiBasicGenTest.cc.

References edm::EventSetup::getHandle().

41  {
42  pdt = iSetup.getHandle(pdtToken_);
43 }
edm::ESGetToken< ParticleDataTable, edm::DefaultRecord > pdtToken_
edm::ESHandle< ParticleDataTable > pdt

Member Data Documentation

MonitorElement* HiBasicGenTest::b[3]
private

Definition at line 33 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdeta[3]
private

Definition at line 31 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdphi[3]
private

Definition at line 34 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::dnchdpt[3]
private

Definition at line 32 of file HiBasicGenTest.h.

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

Definition at line 30 of file HiBasicGenTest.h.

edm::ESHandle<ParticleDataTable> HiBasicGenTest::pdt
private

Definition at line 37 of file HiBasicGenTest.h.

edm::ESGetToken<ParticleDataTable, edm::DefaultRecord> HiBasicGenTest::pdtToken_
private

Definition at line 38 of file HiBasicGenTest.h.

MonitorElement* HiBasicGenTest::rp
private

Definition at line 35 of file HiBasicGenTest.h.