|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
#ifndef CentralityDQM_H
#define CentralityDQM_H
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/HeavyIonEvent/interface/Centrality.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DataFormats/HeavyIonEvent/interface/EvtPlane.h"
#include "RecoHI/HiEvtPlaneAlgos/interface/HiEvtPlaneList.h"
class CentralityDQM : public DQMEDAnalyzer {
public:
explicit CentralityDQM(const edm::ParameterSet& ps);
~CentralityDQM() override;
protected:
void analyze(edm::Event const& e, edm::EventSetup const& eSetup) override;
private:
void bookHistograms(DQMStore::IBooker& bei, edm::Run const&, edm::EventSetup const&) override;
// void bookHistos(DQMStore * bei );
// DQMStore* bei_;
// edm::InputTag centrality_; //CMSS_5_3x
// edm::InputTag vertex_; //CMSS_5_3x
edm::InputTag centralityTag_;
edm::EDGetTokenT<reco::Centrality> centralityToken;
edm::Handle<reco::Centrality> centrality_;
edm::InputTag vertexTag_;
edm::EDGetTokenT<std::vector<reco::Vertex> > vertexToken;
edm::Handle<std::vector<reco::Vertex> > vertex_;
edm::InputTag eventplaneTag_;
edm::EDGetTokenT<reco::EvtPlaneCollection> eventplaneToken;
edm::InputTag centralityBinTag_;
edm::EDGetTokenT<int> centralityBinToken;
edm::Handle<int> centralityBin_;
///////////////////////////
// Histograms
///////////////////////////
// Histograms - Centrality
MonitorElement* h_hiNpix;
MonitorElement* h_hiNpixelTracks;
MonitorElement* h_hiNtracks;
MonitorElement* h_hiNtracksPtCut;
MonitorElement* h_hiNtracksEtaCut;
MonitorElement* h_hiNtracksEtaPtCut;
MonitorElement* h_hiHF;
MonitorElement* h_hiHFplus;
MonitorElement* h_hiHFminus;
MonitorElement* h_hiHFplusEta4;
MonitorElement* h_hiHFminusEta4;
MonitorElement* h_hiHFhit;
MonitorElement* h_hiHFhitPlus;
MonitorElement* h_hiHFhitMinus;
MonitorElement* h_hiEB;
MonitorElement* h_hiET;
MonitorElement* h_hiEE;
MonitorElement* h_hiEEplus;
MonitorElement* h_hiEEminus;
MonitorElement* h_hiZDC;
MonitorElement* h_hiZDCplus;
MonitorElement* h_hiZDCminus;
MonitorElement* h_hiPF;
MonitorElement* h_hiPFplus;
MonitorElement* h_hiPFminus;
MonitorElement* h_vertex_x;
MonitorElement* h_vertex_y;
MonitorElement* h_vertex_z;
MonitorElement* h_cent_bin;
MonitorElement* h_ep_HFm2;
MonitorElement* h_ep_HFp2;
MonitorElement* h_ep_trackmid2;
MonitorElement* h_ep_trackm2;
MonitorElement* h_ep_trackp2;
MonitorElement* h_ep_HFm3;
MonitorElement* h_ep_HFp3;
MonitorElement* h_ep_trackmid3;
};
#endif
|