CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
HcalRecHitsAnalyzer.h
Go to the documentation of this file.
1 #ifndef _DQMOFFLINE_HCAL_HCALRECHITSANALYZER_H_
2 #define _DQMOFFLINE_HCAL_HCALRECHITSANALYZER_H_
3 
11 
13 
16 
20 
32 
34 #include <algorithm>
35 #include <cmath>
36 #include <ostream>
37 #include <string>
38 #include <utility>
39 #include <vector>
40 
42 // channel status
45 
48 
49 // severity level assignment for HCAL
52 
53 // severity level assignment for ECAL
55 
57 public:
59 
60  void analyze(edm::Event const &ev, edm::EventSetup const &) override;
61  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
62  // virtual void beginRun(edm::Run const& run, edm::EventSetup const&)
63  // override;
64  void dqmBeginRun(const edm::Run &run, const edm::EventSetup &) override;
65 
66 private:
67  virtual void fillRecHitsTmp(int subdet_, edm::Event const &ev);
68  double dR(double eta1, double phi1, double eta2, double phi2);
69  double phi12(double phi1, double en1, double phi2, double en2);
70  double dPhiWsign(double phi1, double phi2);
71 
73 
79  bool hep17_;
81  bool famos_;
82 
84 
85  int nChannels_[5]; // 0:any, 1:HB, 2:HE
86 
89 
92 
93  // RecHit Collection input tags
99 
106 
107  // choice of subdetector in config : noise/HB/HE/HO/HF/ALL (0/1/2/3/4/5)
108  int subdet_;
109 
110  // single/multi-particle sample (1/2)
111  int etype_;
112  int iz;
113  int imc;
114 
115  // Hcal topology
116  const HcalTopology *theHcalTopology = nullptr;
117  // for checking the status of ECAL and HCAL channels stored in the DB
119  // calculator of severety level for HCAL
121  int hcalSevLvl(const CaloRecHit *hit);
122 
124 
129 
130  // RecHits counters
131  std::vector<MonitorElement *> Nhb;
132  std::vector<MonitorElement *> Nhe;
133  std::vector<MonitorElement *> Nho;
134  std::vector<MonitorElement *> Nhf;
135 
136  // In ALL other cases : 2D ieta-iphi maps
137  // without and with cuts (a la "Scheme B") on energy
138  // - only in the cone around particle for single-part samples (mc = "yes")
139  // - for all calls in milti-particle samples (mc = "no")
140 
142 
143  std::vector<MonitorElement *> emap;
145 
146  std::vector<MonitorElement *> emean_vs_ieta_HB;
147  std::vector<MonitorElement *> emean_vs_ieta_HBM0;
148  std::vector<MonitorElement *> emean_vs_ieta_HBM3;
149  std::vector<MonitorElement *> emean_vs_ieta_HE;
150  std::vector<MonitorElement *> emean_vs_ieta_HEM0;
151  std::vector<MonitorElement *> emean_vs_ieta_HEM3;
152  std::vector<MonitorElement *> emean_vs_ieta_HEP17;
153  std::vector<MonitorElement *> emean_vs_ieta_HEP17M0;
154  std::vector<MonitorElement *> emean_vs_ieta_HEP17M3;
155  std::vector<MonitorElement *> emean_vs_ieta_HF;
157 
158  std::vector<MonitorElement *> occupancy_map_HB;
159  std::vector<MonitorElement *> occupancy_map_HE;
160  std::vector<MonitorElement *> occupancy_map_HF;
162 
163  std::vector<MonitorElement *> occupancy_vs_ieta_HB;
164  std::vector<MonitorElement *> occupancy_vs_ieta_HE;
165  std::vector<MonitorElement *> occupancy_vs_ieta_HF;
167 
168  std::vector<MonitorElement *> nrechits_vs_iphi_HBP, nrechits_vs_iphi_HBM;
169  std::vector<MonitorElement *> nrechits_vs_iphi_HEP, nrechits_vs_iphi_HEM;
170  std::vector<MonitorElement *> nrechits_vs_iphi_HFP, nrechits_vs_iphi_HFM;
172 
173  // for single monoenergetic particles - cone collection profile vs ieta.
177  // Single particles - deviation of cluster from MC truth
180 
181  // time?
186 
187  // energy of rechits
196 
201  std::vector<MonitorElement *> meRecHitsEnergyHEP17;
202  std::vector<MonitorElement *> meRecHitsEnergyHEP17M0;
203  std::vector<MonitorElement *> meRecHitsEnergyHEP17M3;
208 
211 
214 
222 
228 
233 
238 
243 
250 
253 
258 
259  // 2D plot of sum of RecHits in HCAL as function of ECAL's one
262 
263  // number of ECAL's rechits in cone 0.3
266 
267  CaloGeometry const *geometry = nullptr;
268 
269  // Status word histos
275 
276  // Status word correlation
279 
280  // Aux Status word histos
285 
286  // Filling vectors with essential RecHits data
287  std::vector<int> csub;
288  std::vector<int> cieta;
289  std::vector<int> ciphi;
290  std::vector<int> cdepth;
291  std::vector<double> cen;
292  std::vector<double> cenM0;
293  std::vector<double> cenM3;
294  std::vector<double> cchi2;
295  std::vector<double> ceta;
296  std::vector<double> cphi;
297  std::vector<double> ctime;
298  std::vector<double> cz;
299  std::vector<uint32_t> cstwd;
300  std::vector<uint32_t> cauxstwd;
301  std::vector<int> csevlev;
302 
303  // counter
304  int nevtot;
305 };
306 
307 #endif
MonitorElement * nrechits_vs_iphi_HOP
MonitorElement * sevLvl_HF
std::vector< MonitorElement * > emean_vs_ieta_HEP17M3
MonitorElement * meTEprofileHB_High
MonitorElement * RecHit_StatusWord_HE
std::vector< int > csub
std::vector< int > csevlev
MonitorElement * meRecHitsM2Chi2HB
MonitorElement * meTimeHF
std::vector< double > cphi
std::vector< int > hcalHOSevLvlVec
double dR(double eta1, double phi1, double eta2, double phi2)
MonitorElement * RecHit_StatusWord_HF
MonitorElement * meSumRecHitsEnergyConeHO
MonitorElement * meRecHitsEnergyM3vM2HE
std::vector< MonitorElement * > Nho
std::vector< MonitorElement * > emean_vs_ieta_HBM3
MonitorElement * sevLvl_HE
std::vector< MonitorElement * > occupancy_vs_ieta_HB
MonitorElement * RecHit_Aux_StatusWord_HO
double dPhiWsign(double phi1, double phi2)
MonitorElement * meNumEcalRecHitsConeHE
MonitorElement * meTE_Low_HE
MonitorElement * meRecHitsEnergyHEM3
std::vector< int > cdepth
MonitorElement * meSumRecHitsEnergyHB
std::vector< MonitorElement * > occupancy_map_HE
const HcalTopology * theHcalTopology
int hcalSevLvl(const CaloRecHit *hit)
MonitorElement * meSumRecHitsEnergyHF
edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd > hcalChannelQualityToken_
MonitorElement * meSumRecHitsEnergyConeHF
std::vector< int > hcalHBSevLvlVec
MonitorElement * occupancy_map_HO
MonitorElement * meEcalHcalEnergyConeHO
std::vector< MonitorElement * > occupancy_map_HF
std::vector< double > ceta
edm::EDGetTokenT< EBRecHitCollection > tok_EB_
MonitorElement * meEnConeEtaProfile_EH
void analyze(edm::Event const &ev, edm::EventSetup const &) override
std::vector< int > hcalHESevLvlVec
double phi12(double phi1, double en1, double phi2, double en2)
MonitorElement * meRecHitsEnergyM3vM0HE
std::vector< MonitorElement * > Nhe
MonitorElement * meTE_Low_HB
MonitorElement * meTEprofileHF_Low
const HcalSeverityLevelComputer * theHcalSevLvlComputer
std::vector< MonitorElement * > Nhf
MonitorElement * meTE_HB
MonitorElement * meSumRecHitsEnergyConeHB
MonitorElement * meRecHitsEnergyM2vM0HB
std::vector< MonitorElement * > nrechits_vs_iphi_HEP
MonitorElement * sevLvl_HO
virtual void fillRecHitsTmp(int subdet_, edm::Event const &ev)
std::vector< MonitorElement * > emean_vs_ieta_HEP17
MonitorElement * meEcalHcalEnergyHB
std::vector< MonitorElement * > emean_vs_ieta_HBM0
MonitorElement * meRecHitsEnergyM2vM0HE
MonitorElement * meEcalHcalEnergyHE
std::vector< MonitorElement * > occupancy_vs_ieta_HF
MonitorElement * emap_HO
std::vector< MonitorElement * > occupancy_vs_ieta_HE
MonitorElement * meSumRecHitsEnergyConeHE
std::vector< MonitorElement * > emean_vs_ieta_HB
MonitorElement * meTE_HO
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > caloGeometryRunToken_
const HcalChannelQuality * theHcalChStatus
std::vector< MonitorElement * > occupancy_map_HB
MonitorElement * meEnConeEtaProfile_E
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > caloGeometryEventToken_
std::vector< MonitorElement * > emap
std::vector< MonitorElement * > nrechits_vs_iphi_HFP
MonitorElement * meTEprofileHO
MonitorElement * meTEprofileHB
edm::EDGetTokenT< HORecHitCollection > tok_ho_
std::vector< double > cchi2
MonitorElement * meLog10Chi2profileHE
std::vector< MonitorElement * > nrechits_vs_iphi_HFM
MonitorElement * meRecHitsEnergyHBM0
MonitorElement * meTE_High_HO
MonitorElement * meRecHitsEnergyHO
MonitorElement * meRecHitsM2Chi2HE
MonitorElement * meDeltaEta
MonitorElement * meTE_HF
std::vector< uint32_t > cauxstwd
MonitorElement * meRecHitsEnergyHEM0
MonitorElement * occupancy_vs_ieta_HO
MonitorElement * meTEprofileHE_Low
MonitorElement * RecHit_StatusWord_HO
MonitorElement * meEcalHcalEnergyConeHF
MonitorElement * meRecHitsEnergyHB
std::vector< MonitorElement * > emean_vs_ieta_HE
std::vector< int > ciphi
MonitorElement * nrechits_vs_iphi_HOM
MonitorElement * emean_vs_ieta_HO
MonitorElement * sevLvl_HB
std::vector< MonitorElement * > emean_vs_ieta_HEP17M0
MonitorElement * meDeltaPhi
std::vector< double > cz
void dqmBeginRun(const edm::Run &run, const edm::EventSetup &) override
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
MonitorElement * map_ecal
std::vector< MonitorElement * > nrechits_vs_iphi_HBP
MonitorElement * meTimeHE
std::vector< int > cieta
MonitorElement * meEcalHcalEnergyConeHE
MonitorElement * meSumRecHitsEnergyHE
std::vector< MonitorElement * > nrechits_vs_iphi_HEM
MonitorElement * meTimeHO
MonitorElement * meRecHitsEnergyHF
std::vector< double > cenM3
HcalRecHitsAnalyzer(edm::ParameterSet const &conf)
MonitorElement * meEnConeEtaProfile
std::vector< MonitorElement * > emean_vs_ieta_HEM3
MonitorElement * meRecHitsCleanedEnergyHF
MonitorElement * meTEprofileHE
std::vector< uint32_t > cstwd
std::vector< MonitorElement * > meRecHitsEnergyHEP17
MonitorElement * RecHit_Aux_StatusWord_HB
MonitorElement * meEcalHcalEnergyConeHB
std::vector< MonitorElement * > emean_vs_ieta_HF
MonitorElement * meRecHitsEnergyM3vM2HB
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
std::vector< MonitorElement * > meRecHitsEnergyHEP17M0
MonitorElement * meLog10Chi2profileHB
std::vector< MonitorElement * > nrechits_vs_iphi_HBM
MonitorElement * RecHit_StatusWord_HF67
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > hcalDDDRecConstantsToken_
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
MonitorElement * meRecHitsCleanedEnergyHB
MonitorElement * meTE_High_HB
std::vector< int > hcalHFSevLvlVec
MonitorElement * RecHit_StatusWord_HB
MonitorElement * meTEprofileHO_High
MonitorElement * meEnergyHcalVsEcalHE
MonitorElement * meTEprofileHB_Low
edm::EDGetTokenT< EERecHitCollection > tok_EE_
std::vector< MonitorElement * > Nhb
MonitorElement * meRecHitsCleanedEnergyHE
std::vector< double > ctime
std::vector< double > cen
MonitorElement * meRecHitsEnergyM3vM0HB
MonitorElement * meRecHitsEnergyHBM3
std::vector< double > cenM0
MonitorElement * RecHit_StatusWordCorr_HE
MonitorElement * meTimeHB
MonitorElement * meTEprofileHF
MonitorElement * RecHit_Aux_StatusWord_HF
MonitorElement * RecHit_Aux_StatusWord_HE
edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd > hcalSeverityLevelComputerToken_
MonitorElement * RecHit_StatusWordCorr_HB
MonitorElement * meEnergyHcalVsEcalHB
std::vector< MonitorElement * > meRecHitsEnergyHEP17M3
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > hcalTopologyToken_
Definition: Run.h:45
MonitorElement * meSumRecHitsEnergyConeHFL
MonitorElement * meRecHitsCleanedEnergyHO
MonitorElement * meTE_HE
MonitorElement * meNumEcalRecHitsConeHB
std::vector< MonitorElement * > emean_vs_ieta_HEM0
MonitorElement * meSumRecHitsEnergyHO
MonitorElement * meSumRecHitsEnergyConeHFS
MonitorElement * meRecHitsEnergyHE
MonitorElement * meTE_Low_HF