CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalRecHitsValidation.h
Go to the documentation of this file.
1 #ifndef HcalRecHitsValidation_H
2 #define HcalRecHitsValidation_H
3 
14 
16 
19 
23 
30 
33 
34 #include <vector>
35 #include <utility>
36 #include <ostream>
37 #include <string>
38 #include <algorithm>
39 #include <cmath>
41 
43 // channel status
46 
49 
50 // severity level assignment for HCAL
53 
54 // severity level assignment for ECAL
56 
57 
59  public:
62  virtual void analyze(edm::Event const& ev, edm::EventSetup const& c);
63  virtual void beginJob() ;
64  virtual void endJob() ;
65  private:
66 
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 
80  bool famos_;
82 
83  //RecHit Collection input tags
91 
92  // choice of subdetector in config : noise/HB/HE/HO/HF/ALL (0/1/2/3/4/5)
93  int subdet_;
94 
95  // single/multi-particle sample (1/2)
96  int etype_;
97  int iz;
98  int imc;
99 
100  // for checking the status of ECAL and HCAL channels stored in the DB
102  // calculator of severety level for HCAL
104  int hcalSevLvl(const CaloRecHit* hit);
105 
107 
112 
113  // RecHits counters
118 
119  // ZS-specific
120 
125 
134 
143 
152 
153  // In ALL other cases : 2D ieta-iphi maps
154  // without and with cuts (a la "Scheme B") on energy
155  // - only in the cone around particle for single-part samples (mc = "yes")
156  // - for all calls in milti-particle samples (mc = "no")
157 
159 
164 
173 
182 
191 
200 
209 
218 
227 
228 
229  // also - energy in the cone around MC particle
234 
235  // for single monoenergetic particles - cone collection profile vs ieta.
243  // Single particles - deviation of cluster from MC truth
246  MonitorElement* meDeltaPhiS; // simcluster
247  MonitorElement* meDeltaEtaS; // simculster
248 
249  //----------- NOISE case
255 
256  // number of rechits above threshold 1GEV
260 
261  // number of rechits in the cone
266 
267  // time?
272 
273  // energy of rechits
278 
287 
294 
299 
306 
307 
312 
313 
320 
321 
324 
329 
330  // Histo (2D plot) for sum of RecHits vs SimHits (hcal only)
337  // profile histo (2D plot) for sum of RecHits vs SimHits (hcal only)
344 
345  // 2D plot of sum of RecHits in HCAL as function of ECAL's one
348 
349  // number of ECAL's rechits in cone 0.3
352 
354 
355  //Status word histos
361 
362  //Aux Status word histos
367 
368  // Filling vectors with essential RecHits data
369  std::vector<int> csub;
370  std::vector<int> cieta;
371  std::vector<int> ciphi;
372  std::vector<int> cdepth;
373  std::vector<double> cen;
374  std::vector<double> ceta;
375  std::vector<double> cphi;
376  std::vector<double> ctime;
377  std::vector<double> cz;
378  std::vector<uint32_t> cstwd;
379  std::vector<uint32_t> cauxstwd;
380 
381  // array or min. e-values ieta x iphi x depth x subdet
382  double emap_min[82][72][4][4];
383 
384  // counter
385  int nevtot;
386 
387 };
388 
389 #endif
MonitorElement * emean_vs_ieta_HB2
MonitorElement * RMS_vs_ieta_HB1
MonitorElement * occupancy_vs_ieta_HO
const HcalSeverityLevelComputer * theHcalSevLvlComputer
MonitorElement * meSumRecHitsEnergyConeHE
MonitorElement * occupancy_seqHB1
const HcalChannelQuality * theHcalChStatus
MonitorElement * occupancy_map_HB1
MonitorElement * occupancy_vs_ieta_HE3
edm::EDGetTokenT< edm::PCaloHitContainer > tok_hh_
MonitorElement * meRecHitsEnergyHF
MonitorElement * occupancy_vs_ieta_HF1
MonitorElement * RMS_vs_ieta_HF1
double phi12(double phi1, double en1, double phi2, double en2)
std::vector< int > hcalHBSevLvlVec
MonitorElement * meNumRecHitsConeHE
MonitorElement * meSumRecHitsEnergyConeHB
MonitorElement * meTEprofileHE_Low
std::vector< uint32_t > cstwd
MonitorElement * RMS_vs_ieta_HE3
MonitorElement * RecHit_Aux_StatusWord_HB
MonitorElement * RMS_vs_ieta_HE1
MonitorElement * emean_vs_ieta_HB1
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
MonitorElement * occupancy_seqHO
MonitorElement * occupancy_seqHE3
MonitorElement * meEnConeEtaProfile
bool ev
MonitorElement * occupancy_map_HE3
double dPhiWsign(double phi1, double phi2)
MonitorElement * meRecHitSimHitProfileHF
MonitorElement * meEcalHcalEnergyHB
int hcalSevLvl(const CaloRecHit *hit)
MonitorElement * meRecHitSimHitProfileHFS
MonitorElement * meEnergyHcalVsEcalHE
MonitorElement * occupancy_map_HF1
MonitorElement * meEcalHcalEnergyConeHE
MonitorElement * meRecHitsEnergyHB
MonitorElement * meSumRecHitsEnergyHF
double emap_min[82][72][4][4]
MonitorElement * meNumRecHitsConeHO
std::vector< double > ceta
MonitorElement * occupancy_seqHF2
MonitorElement * meNumRecHitsConeHB
MonitorElement * RMS_vs_ieta_HB2
MonitorElement * meRecHitSimHitHF
MonitorElement * meRecHitsEnergyHE
MonitorElement * occupancy_map_HE1
MonitorElement * meRecHitSimHitProfileHE
MonitorElement * meSumRecHitsEnergyHB
edm::EDGetTokenT< HORecHitCollection > tok_ho_
MonitorElement * meEnConeEtaProfile_E
MonitorElement * meRecHitSimHitHB
MonitorElement * meRecHitSimHitHFS
MonitorElement * RMS_vs_ieta_HO
MonitorElement * meSumRecHitsEnergyConeHFS
MonitorElement * RecHit_Aux_StatusWord_HO
std::vector< double > cphi
MonitorElement * meNumEcalRecHitsConeHB
MonitorElement * meNumRecHitsThreshHB
MonitorElement * meSumRecHitsEnergyHO
std::vector< int > hcalHESevLvlVec
MonitorElement * meEnergyHcalVsEcalHB
MonitorElement * RecHit_StatusWord_HF
MonitorElement * RecHit_Aux_StatusWord_HE
edm::ESHandle< CaloGeometry > geometry
MonitorElement * meSumRecHitsEnergyConeHO
MonitorElement * meEcalHcalEnergyConeHO
MonitorElement * meTEprofileHB
MonitorElement * map_econe_depth1
virtual void fillRecHitsTmp(int subdet_, edm::Event const &ev)
MonitorElement * map_econe_depth4
double dR(double eta1, double phi1, double eta2, double phi2)
MonitorElement * emean_vs_ieta_HE3
MonitorElement * meNumRecHitsThreshHO
MonitorElement * meEcalHcalEnergyConeHB
MonitorElement * meNumRecHitsThreshHE
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
std::vector< int > hcalHFSevLvlVec
MonitorElement * occupancy_map_HE2
MonitorElement * occupancy_vs_ieta_HB2
MonitorElement * meSumRecHitsEnergyConeHF
MonitorElement * RecHit_Aux_StatusWord_HF
MonitorElement * meRecHitSimHitHO
MonitorElement * meRecHitSimHitProfileHB
MonitorElement * RecHit_StatusWord_HO
MonitorElement * meEnConeEtaProfile_depth1
tuple conf
Definition: dbtoconf.py:185
MonitorElement * emean_vs_ieta_HF2
virtual void analyze(edm::Event const &ev, edm::EventSetup const &c)
MonitorElement * occupancy_map_HO
MonitorElement * occupancy_vs_ieta_HE2
MonitorElement * meTEprofileHF_Low
edm::EDGetTokenT< EBRecHitCollection > tok_EB_
MonitorElement * occupancy_map_HF2
MonitorElement * RecHit_StatusWord_HF67
MonitorElement * meTEprofileHO_High
MonitorElement * meTEprofileHE
MonitorElement * occupancy_seqHB2
MonitorElement * occupancy_map_HB2
MonitorElement * RMS_vs_ieta_HE2
MonitorElement * meNumRecHitsConeHF
MonitorElement * meTEprofileHF
MonitorElement * meEnConeEtaProfile_depth3
MonitorElement * emean_vs_ieta_HE1
MonitorElement * meRecHitSimHitProfileHO
std::vector< int > hcalHOSevLvlVec
MonitorElement * emean_vs_ieta_HE2
MonitorElement * occupancy_vs_ieta_HF2
edm::EDGetTokenT< EERecHitCollection > tok_EE_
MonitorElement * meEcalHcalEnergyConeHF
std::vector< double > ctime
MonitorElement * map_econe_depth2
MonitorElement * RecHit_StatusWord_HE
MonitorElement * meEnConeEtaProfile_EH
MonitorElement * meRecHitSimHitProfileHFL
std::vector< double > cz
MonitorElement * meEnConeEtaProfile_depth2
MonitorElement * RMS_vs_ieta_HF2
MonitorElement * meRecHitsEnergyHO
std::vector< uint32_t > cauxstwd
MonitorElement * emean_vs_ieta_HO
MonitorElement * occupancy_seqHF1
MonitorElement * meSumRecHitsEnergyHE
MonitorElement * meRecHitSimHitHFL
MonitorElement * meSumRecHitsEnergyConeHFL
MonitorElement * meNumEcalRecHitsConeHE
MonitorElement * meEcalHcalEnergyHE
MonitorElement * emean_vs_ieta_HF1
MonitorElement * meTEprofileHO
edm::EDGetTokenT< edm::HepMCProduct > tok_evt_
MonitorElement * meTEprofileHB_Low
MonitorElement * occupancy_vs_ieta_HB1
MonitorElement * occupancy_seqHE2
MonitorElement * meRecHitSimHitHE
HcalRecHitsValidation(edm::ParameterSet const &conf)
MonitorElement * map_econe_depth3
MonitorElement * occupancy_seqHE1
MonitorElement * meTEprofileHB_High
std::vector< double > cen
MonitorElement * meEnConeEtaProfile_depth4
MonitorElement * occupancy_vs_ieta_HE1
MonitorElement * RecHit_StatusWord_HB