test
CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HGCalHitValidation.cc
Go to the documentation of this file.
1 //
3 // Package: HGCalHitValidation
4 // Class: HGCalHitValidation
5 //
26 
34 
47 
53 
54 #include <cmath>
55 #include <memory>
56 #include <iostream>
57 #include <string>
58 #include <vector>
59 
60 //#define EDM_ML_DEBUG
61 
63 
64 public:
65 
66  explicit HGCalHitValidation( const edm::ParameterSet& );
68  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
69 
70 protected:
71  void dqmBeginRun(edm::Run const&, edm::EventSetup const&) override;
72  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
73  void analyze(const edm::Event&, const edm::EventSetup&) override;
74 
75 private:
76  typedef std::tuple<float,float,float,float> HGCHitTuple;
77 
78  //HGC Geometry
79  std::vector<const HGCalDDDConstants*> hgcCons_;
80  std::vector<const HGCalGeometry*> hgcGeometry_;
84  std::vector<std::string> geometrySource_;
85  std::vector<int> ietaExcludeBH_;
86 
94 
95  //histogram related stuff
99 
104 
108 
109 };
110 
111 
113 
114  geometrySource_ = cfg.getUntrackedParameter< std::vector<std::string> >("geometrySource");
115  eeSimHitToken_ = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("eeSimHitSource"));
116  fhSimHitToken_ = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("fhSimHitSource"));
117  bhSimHitToken_ = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("bhSimHitSource"));
118  eeRecHitToken_ = consumes<HGCeeRecHitCollection>(cfg.getParameter<edm::InputTag>("eeRecHitSource"));
119  fhRecHitToken_ = consumes<HGChefRecHitCollection>(cfg.getParameter<edm::InputTag>("fhRecHitSource"));
120  bhRecHitToken_ = consumes<HBHERecHitCollection>(cfg.getParameter<edm::InputTag>("bhRecHitSource"));
121  ietaExcludeBH_ = cfg.getParameter<std::vector<int> >("ietaExcludeBH");
122 
123 #ifdef EDM_ML_DEBUG
124  edm::LogInfo("HGCalValid") << "Exclude the following "
125  << ietaExcludeBH_.size()
126  << " ieta values from BH plots";
127  for (unsigned int k=0; k<ietaExcludeBH_.size(); ++k)
128  edm::LogInfo("HGCalValid") << " [" << k << "] " << ietaExcludeBH_[k];
129 #endif
130 }
131 
133 
135  //The following says we do not know what parameters are allowed so do no validation
136  // Please change this to state exactly what you do use, even if it is no parameters
138  desc.setUnknown();
139  descriptions.addDefault(desc);
140 }
141 
143  edm::Run const&,
144  edm::EventSetup const&) {
145 
146  iB.setCurrentFolder("HGCalSimHitsV/HitValidation");
147 
148  //initiating histograms
149  heedzVsZ = iB.book2D("heedzVsZ","",7200,-360,360,100,-0.1,0.1);
150  heedyVsY = iB.book2D("heedyVsY","",400,-200,200,100,-0.02,0.02);
151  heedxVsX = iB.book2D("heedxVsX","",400,-200,200,100,-0.02,0.02);
152  heeRecVsSimZ = iB.book2D("heeRecVsSimZ","",7200,-360,360,7200,-360,360);
153  heeRecVsSimY = iB.book2D("heeRecVsSimY","",400,-200,200,400,-200,200);
154  heeRecVsSimX = iB.book2D("heeRecVsSimX","",400,-200,200,400,-200,200);
155 
156  hefdzVsZ = iB.book2D("hefdzVsZ","",8200,-410,410,100,-0.1,0.1);
157  hefdyVsY = iB.book2D("hefdyVsY","",400,-200,200,100,-0.02,0.02);
158  hefdxVsX = iB.book2D("hefdxVsX","",400,-200,200,100,-0.02,0.02);
159  hefRecVsSimZ = iB.book2D("hefRecVsSimZ","",8200,-410,410,8200,-410,410);
160  hefRecVsSimY = iB.book2D("hefRecVsSimY","",400,-200,200,400,-200,200);
161  hefRecVsSimX = iB.book2D("hefRecVsSimX","",400,-200,200,400,-200,200);
162 
163  hebdzVsZ = iB.book2D("hebdzVsZ","",1080,-540,540,100,-1.0,1.0);
164  hebdPhiVsPhi = iB.book2D("hebdPhiVsPhi","",M_PI*100,-0.5,M_PI+0.5,200,-0.2,0.2);
165  hebdEtaVsEta = iB.book2D("hebdEtaVsEta","",1000,-5,5,200,-0.1,0.1);
166  hebRecVsSimZ = iB.book2D("hebRecVsSimZ","",1080,-540,540,1080,-540,540);
167  hebRecVsSimY = iB.book2D("hebRecVsSimY","",400,-200,200,400,-200,200);
168  hebRecVsSimX = iB.book2D("hebRecVsSimX","",400,-200,200,400,-200,200);
169 
170  heeEnRec = iB.book1D("heeEnRec","",1000,0,10);
171  heeEnSim = iB.book1D("heeEnSim","",1000,0,0.01);
172  heeEnSimRec = iB.book2D("heeEnSimRec","",200,0,0.002,200,0,0.2);
173 
174  hefEnRec = iB.book1D("hefEnRec","",1000,0,10);
175  hefEnSim = iB.book1D("hefEnSim","",1000,0,0.01);
176  hefEnSimRec = iB.book2D("hefEnSimRec","",200,0,0.001,200,0,0.5);
177 
178  hebEnRec = iB.book1D("hebEnRec","",1000,0,15);
179  hebEnSim = iB.book1D("hebEnSim","",1000,0,0.01);
180  hebEnSimRec = iB.book2D("hebEnSimRec","",200,0,0.02,200,0,4);
181 
182 }
183 
185  edm::EventSetup const& iSetup) {
186  //initiating hgc Geometry
187  for (size_t i=0; i<geometrySource_.size(); i++) {
188  if (geometrySource_[i].find("Hcal") != std::string::npos) {
190  iSetup.get<HcalSimNumberingRecord>().get(pHSNDC);
191  if (pHSNDC.isValid()) {
192  hcCons_ = pHSNDC.product();
193  hgcCons_.push_back(0);
194  } else {
195  edm::LogWarning("HGCalValid") << "Cannot initiate HcalDDDSimConstants: "
196  << geometrySource_[i] << std::endl;
197  }
199  iSetup.get<HcalRecNumberingRecord>().get(pHRNDC);
200  if (pHRNDC.isValid()) {
201  hcConr_ = pHRNDC.product();
202  } else {
203  edm::LogWarning("HGCalValid") << "Cannot initiate HcalDDDRecConstants: "
204  << geometrySource_[i] << std::endl;
205  }
207  iSetup.get<CaloGeometryRecord>().get(caloG);
208  if (caloG.isValid()) {
209  const CaloGeometry* geo = caloG.product();
211  hgcGeometry_.push_back(0);
212  } else {
213  edm::LogWarning("HGCalValid") << "Cannot initiate HcalGeometry for "
214  << geometrySource_[i] << std::endl;
215  }
216  } else {
218  iSetup.get<IdealGeometryRecord>().get(geometrySource_[i],hgcCons);
219  if (hgcCons.isValid()) {
220  hgcCons_.push_back(hgcCons.product());
221  } else {
222  edm::LogWarning("HGCalValid") << "Cannot initiate HGCalDDDConstants for "
223  << geometrySource_[i] << std::endl;
224  }
226  iSetup.get<IdealGeometryRecord>().get(geometrySource_[i],hgcGeom);
227  if(hgcGeom.isValid()) {
228  hgcGeometry_.push_back(hgcGeom.product());
229  } else {
230  edm::LogWarning("HGCalValid") << "Cannot initiate HGCalGeometry for "
231  << geometrySource_[i] << std::endl;
232  }
233  }
234  }
235 }
236 
238  std::map<unsigned int, HGCHitTuple> eeHitRefs, fhHitRefs, bhHitRefs;
239 
240  //declare topology and DDD constants
241  const HGCalTopology &heeTopo=hgcGeometry_[0]->topology();
242  const HGCalTopology &hefTopo=hgcGeometry_[1]->topology();
243 
244  //Accesing ee simhits
246  iEvent.getByToken(eeSimHitToken_, eeSimHits);
247 
248  if (eeSimHits.isValid()) {
249  for (std::vector<PCaloHit>::const_iterator simHit = eeSimHits->begin(); simHit != eeSimHits->end(); ++simHit) {
250  int subdet, zside, layer, wafer, celltype, cell;
251  HGCalTestNumbering::unpackHexagonIndex(simHit->id(), subdet, zside, layer, wafer, celltype, cell);
252  std::pair<float, float> xy = hgcCons_[0]->locateCell(cell,layer,wafer,false);
253  float zp = hgcCons_[0]->waferZ(layer,false);
254  if (zside < 0) zp = -zp;
255  float xp = (zp<0) ? -xy.first/10 : xy.first/10;
256  float yp = xy.second/10.0;
257 
258  //skip this hit if after ganging it is not valid
259  std::pair<int,int> recoLayerCell=hgcCons_[0]->simToReco(cell,layer,wafer,heeTopo.detectorType());
260  cell = recoLayerCell.first;
261  layer = recoLayerCell.second;
262 
263  //skip this hit if after ganging it is not valid
264  if (layer<0 || cell<0) {
265  } else {
266 
267  //assign the RECO DetId
268  HGCalDetId id = HGCalDetId((ForwardSubdetector)(subdet),zside,layer,celltype,wafer,cell);
269  float energy = simHit->energy();
270 
271  float energySum(0);
272  if (eeHitRefs.count(id.rawId()) != 0) energySum = std::get<0>(eeHitRefs[id.rawId()]);
273  energySum += energy;
274  eeHitRefs[id.rawId()] = std::make_tuple(energySum,xp,yp,zp);
275  heeEnSim->Fill(energy);
276  }
277  }
278 #ifdef EDM_ML_DEBUG
279  for (std::map<unsigned int,HGCHitTuple>::iterator itr=eeHitRefs.begin();
280  itr != eeHitRefs.end(); ++itr) {
281  int idx = std::distance(eeHitRefs.begin(),itr);
282  edm::LogInfo("HGCalValid") << "EEHit[" << idx << "] " << std::hex
283  << itr->first << std::dec << "; Energy "
284  << std::get<0>(itr->second)
285  << "; Position (" << std::get<1>(itr->second)
286  << ", " << std::get<2>(itr->second) <<", "
287  << std::get<3>(itr->second) << ")" <<std::endl;
288  }
289 #endif
290  } else {
291  edm::LogWarning("HGCalValid") << "No EE SimHit Found " << std::endl;
292  }
293 
294  //Accesing fh simhits
296  iEvent.getByToken(fhSimHitToken_, fhSimHits);
297  if (fhSimHits.isValid()) {
298  for (std::vector<PCaloHit>::const_iterator simHit = fhSimHits->begin();
299  simHit != fhSimHits->end();++simHit) {
300  int subdet, zside, layer, wafer, celltype, cell;
301  HGCalTestNumbering::unpackHexagonIndex(simHit->id(), subdet, zside, layer, wafer, celltype, cell);
302  std::pair<float, float> xy = hgcCons_[1]->locateCell(cell,layer,wafer,false);
303  float zp = hgcCons_[1]->waferZ(layer,false);
304  if (zside < 0) zp = -zp;
305  float xp = (zp<0) ? -xy.first/10 : xy.first/10;
306  float yp = xy.second/10.0;
307 
308  //skip this hit if after ganging it is not valid
309  std::pair<int,int> recoLayerCell = hgcCons_[1]->simToReco(cell,layer,wafer,hefTopo.detectorType());
310  cell = recoLayerCell.first;
311  layer = recoLayerCell.second;
312  //skip this hit if after ganging it is not valid
313  if(layer<0 || cell<0) {
314  } else {
315  //assign the RECO DetId
316  HGCalDetId id = HGCalDetId((ForwardSubdetector)(subdet),zside,layer,celltype,wafer,cell);
317 
318  float energy = simHit->energy();
319  float energySum(0);
320  if (fhHitRefs.count(id.rawId()) != 0) energySum = std::get<0>(fhHitRefs[id.rawId()]);
321  energySum += energy;
322  fhHitRefs[id.rawId()] = std::make_tuple(energySum,xp,yp,zp);
323  hefEnSim->Fill(energy);
324  }
325  }
326 #ifdef EDM_ML_DEBUG
327  for (std::map<unsigned int,HGCHitTuple>::iterator itr=fhHitRefs.begin();
328  itr != fhHitRefs.end(); ++itr) {
329  int idx = std::distance(fhHitRefs.begin(),itr);
330  edm::LogInfo("HGCalValid") << "FHHit[" << idx << "] " << std::hex
331  << itr->first << std::dec << "; Energy "
332  << std::get<0>(itr->second) << "; Position ("
333  << std::get<1>(itr->second) << ", "
334  << std::get<2>(itr->second) <<", "
335  << std::get<3>(itr->second) << ")" <<std::endl;
336  }
337 #endif
338  } else {
339  edm::LogWarning("HGCalValid") << "No FH SimHit Found " << std::endl;
340  }
341 
342  //Accessing bh simhits
344  iEvent.getByToken(bhSimHitToken_, bhSimHits);
345  if (bhSimHits.isValid()) {
346  for (std::vector<PCaloHit>::const_iterator simHit = bhSimHits->begin();
347  simHit != bhSimHits->end(); ++simHit) {
348  int subdet, z, depth, eta, phi, lay;
349  HcalTestNumbering::unpackHcalIndex(simHit->id(), subdet, z, depth, eta, phi, lay);
350  HcalCellType::HcalCell cell = hcCons_->cell(subdet, z, lay, eta, phi);
351 
352  double zp = cell.rz/10;
353 
354  if (subdet == static_cast<int>(HcalEndcap)) {
355  HcalDDDRecConstants::HcalID idx = hcConr_->getHCID(subdet,eta,phi,lay,depth);
356  int sign = (z==0)?(-1):(1);
357  zp *= sign;
358  HcalDetId id = HcalDetId(HcalEndcap,sign*idx.eta,idx.phi,idx.depth);
359 
360  float energy = simHit->energy();
361  float energySum(0);
362  if (bhHitRefs.count(id.rawId()) != 0) energySum = std::get<0>(bhHitRefs[id.rawId()]);
363  energySum += energy;
364  hebEnSim->Fill(energy);
365  if (std::find(ietaExcludeBH_.begin(),ietaExcludeBH_.end(),idx.eta) ==
366  ietaExcludeBH_.end()) {
367  bhHitRefs[id.rawId()] = std::make_tuple(energySum,cell.eta,cell.phi,zp);
368 #ifdef EDM_ML_DEBUG
369  edm::LogInfo("HGCalValid") << "Accept " << id << std::endl;
370  } else {
371  edm::LogInfo("HGCalValid") << "Reject " << id << std::endl;
372 #endif
373  }
374  }
375  }
376 #ifdef EDM_ML_DEBUG
377  for (std::map<unsigned int,HGCHitTuple>::iterator itr=bhHitRefs.begin();
378  itr != bhHitRefs.end(); ++itr) {
379  int idx = std::distance(bhHitRefs.begin(),itr);
380  edm::LogInfo("HGCalValid") << "BHHit[" << idx << "] " << std::hex
381  << itr->first << std::dec << "; Energy "
382  << std::get<0>(itr->second) << "; Position ("
383  << std::get<1>(itr->second) << ", "
384  << std::get<2>(itr->second) <<", "
385  << std::get<3>(itr->second) << ")" <<std::endl;
386  }
387 #endif
388  } else {
389  edm::LogWarning("HGCalValid") << "No BH SimHit Found " << std::endl;
390  }
391 
392  //accessing EE Rechit information
394  iEvent.getByToken(eeRecHitToken_, eeRecHit);
395  if (eeRecHit.isValid()) {
396  const HGCeeRecHitCollection* theHits = (eeRecHit.product());
397  for (auto it = theHits->begin(); it != theHits->end(); ++it) {
398  double energy = it->energy();
399  heeEnRec->Fill(energy);
400  std::map<unsigned int, HGCHitTuple>::const_iterator itr = eeHitRefs.find(it->id().rawId());
401  if (itr != eeHitRefs.end()) {
402  GlobalPoint xyz = hgcGeometry_[0]->getPosition(it->id());
403  heeRecVsSimX->Fill(std::get<1>(itr->second),xyz.x());
404  heeRecVsSimY->Fill(std::get<2>(itr->second),xyz.y());
405  heeRecVsSimZ->Fill(std::get<3>(itr->second),xyz.z());
406  heedxVsX->Fill(std::get<1>(itr->second),(xyz.x()-std::get<1>(itr->second)));
407  heedyVsY->Fill(std::get<2>(itr->second),(xyz.y()-std::get<2>(itr->second)));
408  heedzVsZ->Fill(std::get<3>(itr->second),(xyz.z()-std::get<3>(itr->second)));
409  heeEnSimRec->Fill(std::get<0>(itr->second),energy);
410 #ifdef EDM_ML_DEBUG
411  edm::LogInfo("HGCalValid") << "EEHit: " << std::hex << it->id().rawId()
412  << std::dec << " Sim ("
413  << std::get<0>(itr->second) << ", "
414  << std::get<1>(itr->second) << ", "
415  << std::get<2>(itr->second) << ", "
416  << std::get<3>(itr->second) << ") Rec ("
417  << energy << ", " << xyz.x() << ", "
418  << xyz.y() << ", " << xyz.z() << ")\n";
419 #endif
420  }
421  }
422  } else {
423  edm::LogWarning("HGCalValid") << "No EE RecHit Found " << std::endl;
424  }
425 
426  //accessing FH Rechit information
428  iEvent.getByToken(fhRecHitToken_, fhRecHit);
429  if (fhRecHit.isValid()) {
430  const HGChefRecHitCollection* theHits = (fhRecHit.product());
431  for (auto it = theHits->begin(); it!=theHits->end(); ++it) {
432  double energy = it->energy();
433  hefEnRec->Fill(energy);
434  std::map<unsigned int, HGCHitTuple>::const_iterator itr = fhHitRefs.find(it->id().rawId());
435  if (itr != fhHitRefs.end()) {
436  GlobalPoint xyz = hgcGeometry_[1]->getPosition(it->id());
437 
438  hefRecVsSimX->Fill(std::get<1>(itr->second),xyz.x());
439  hefRecVsSimY->Fill(std::get<2>(itr->second),xyz.y());
440  hefRecVsSimZ->Fill(std::get<3>(itr->second),xyz.z());
441  hefdxVsX->Fill(std::get<1>(itr->second),(xyz.x()-std::get<1>(itr->second)));
442  hefdyVsY->Fill(std::get<2>(itr->second),(xyz.y()-std::get<2>(itr->second)));
443  hefdzVsZ->Fill(std::get<3>(itr->second),(xyz.z()-std::get<3>(itr->second)));
444  hefEnSimRec->Fill(std::get<0>(itr->second),energy);
445 #ifdef EDM_ML_DEBUG
446  edm::LogInfo("HGCalValid") << "FHHit: " << std::hex << it->id().rawId()
447  << std::dec << " Sim ("
448  << std::get<0>(itr->second) << ", "
449  << std::get<1>(itr->second) << ", "
450  << std::get<2>(itr->second) << ", "
451  << std::get<3>(itr->second) << ") Rec ("
452  << energy << "," << xyz.x() << ", "
453  << xyz.y() << ", " << xyz.z() << ")\n";
454 #endif
455  }
456  }
457  } else {
458  edm::LogWarning("HGCalValid") << "No FH RecHit Found " << std::endl;
459  }
460 
461 
462  //accessing BH Rechit information
464  iEvent.getByToken(bhRecHitToken_, bhRecHit);
465  if (bhRecHit.isValid()) {
466  const HBHERecHitCollection* theHits = (bhRecHit.product());
467 
468  for (auto it = theHits->begin(); it!=theHits->end(); ++it) {
469  DetId id = it->id();
470  if (id.subdetId() == (int)(HcalEndcap)) {
471  double energy = it->energy();
472  hebEnRec->Fill(energy);
474 
475  std::map<unsigned int, HGCHitTuple>::const_iterator itr = bhHitRefs.find(id.rawId());
476  if (itr != bhHitRefs.end()) {
477  float ang3 = xyz.phi().value(); // returns the phi in radians
478  double fac = sinh(std::get<1>(itr->second));
479  double pT = std::get<3>(itr->second) / fac;
480  double xp = pT * cos(std::get<2>(itr->second));
481  double yp = pT * sin(std::get<2>(itr->second));
482  hebRecVsSimX->Fill(xp,xyz.x());
483  hebRecVsSimY->Fill(yp,xyz.y());
484  hebRecVsSimZ->Fill(std::get<3>(itr->second),xyz.z());
485  hebdEtaVsEta->Fill(std::get<1>(itr->second),(xyz.eta()-std::get<1>(itr->second)));
486  hebdPhiVsPhi->Fill(std::get<2>(itr->second),(ang3-std::get<2>(itr->second)));
487  hebdzVsZ->Fill(std::get<3>(itr->second),(xyz.z()-std::get<3>(itr->second)));
488  hebEnSimRec->Fill(std::get<0>(itr->second),energy);
489 
490 #ifdef EDM_ML_DEBUG
491  edm::LogInfo("HGCalValid") << "BHHit: " << std::hex << id.rawId()
492  << std::dec << " Sim ("
493  << std::get<0>(itr->second) << ", "
494  << std::get<1>(itr->second) << ", "
495  << std::get<2>(itr->second) << ", "
496  << std::get<3>(itr->second) << ") Rec ("
497  << energy << ", " << xyz.x() << ", "
498  << xyz.y() << ", " << xyz.z() << ")\n";
499 #endif
500  }
501  }
502  }
503  } else {
504  edm::LogWarning("HGCalValid") << "No BH RecHit Found " << std::endl;
505  }
506 
507 }
508 
509 //define this as a plug-in
T getParameter(std::string const &) const
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:45
edm::EDGetTokenT< std::vector< PCaloHit > > fhSimHitToken_
T getUntrackedParameter(std::string const &, T const &) const
int i
Definition: DBlmapReader.cc:9
tuple cfg
Definition: looper.py:293
MonitorElement * hebRecVsSimX
MonitorElement * hefRecVsSimY
MonitorElement * hefRecVsSimX
edm::InputTag bhSimHitSource
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
MonitorElement * hebRecVsSimY
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
double sign(double x)
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
T y() const
Definition: PV3DBase.h:63
MonitorElement * hefEnSimRec
MonitorElement * heedzVsZ
bool detectorType() const
int zside(DetId const &)
MonitorElement * heeRecVsSimX
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
ForwardSubdetector
std::vector< const HGCalDDDConstants * > hgcCons_
MonitorElement * heeRecVsSimZ
const HcalDDDSimConstants * hcCons_
MonitorElement * heeRecVsSimY
MonitorElement * heeEnRec
virtual const CaloCellGeometry * getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
void Fill(long long x)
std::vector< const HGCalGeometry * > hgcGeometry_
MonitorElement * hebEnRec
edm::EDGetTokenT< HGCeeRecHitCollection > eeRecHitToken_
std::tuple< float, float, float, float > HGCHitTuple
edm::EDGetTokenT< std::vector< PCaloHit > > eeSimHitToken_
edm::EDGetTokenT< HBHERecHitCollection > bhRecHitToken_
int iEvent
Definition: GenABIO.cc:230
void addDefault(ParameterSetDescription const &psetDescription)
MonitorElement * hebEnSim
MonitorElement * hefdyVsY
MonitorElement * heeEnSimRec
std::vector< int > ietaExcludeBH_
MonitorElement * hefEnSim
HcalID getHCID(int subdet, int ieta, int iphi, int lay, int idepth) const
T z() const
Definition: PV3DBase.h:64
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
MonitorElement * hebRecVsSimZ
std::vector< std::string > geometrySource_
MonitorElement * hefEnRec
const CaloSubdetectorGeometry * hcGeometry_
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
MonitorElement * heedxVsX
static void unpackHcalIndex(const uint32_t &idx, int &det, int &z, int &depth, int &eta, int &phi, int &lay)
void analyze(const edm::Event &, const edm::EventSetup &) override
MonitorElement * heeEnSim
edm::EDGetTokenT< std::vector< PCaloHit > > bhSimHitToken_
bool isValid() const
Definition: HandleBase.h:75
MonitorElement * hebdPhiVsPhi
edm::InputTag fhSimHitSource
#define M_PI
const_iterator end() const
MonitorElement * hebEnSimRec
T value() const
Explicit access to value in case implicit conversion not OK.
Definition: Phi.h:38
Definition: DetId.h:18
edm::InputTag eeSimHitSource
HGCalHitValidation(const edm::ParameterSet &)
MonitorElement * hefdxVsX
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
T const * product() const
Definition: Handle.h:81
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
void dqmBeginRun(edm::Run const &, edm::EventSetup const &) override
const T & get() const
Definition: EventSetup.h:56
T const * product() const
Definition: ESHandle.h:86
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
T eta() const
Definition: PV3DBase.h:76
const HcalDDDRecConstants * hcConr_
MonitorElement * hefRecVsSimZ
HcalCellType::HcalCell cell(int det, int zside, int depth, int etaR, int iphi) const
MonitorElement * hefdzVsZ
double energySum(const DataFrame &df, int fs, int ls)
bool isValid() const
Definition: ESHandle.h:47
const GlobalPoint & getPosition() const
Returns the position of reference for this cell.
MonitorElement * hebdzVsZ
T x() const
Definition: PV3DBase.h:62
MonitorElement * heedyVsY
static void unpackHexagonIndex(const uint32_t &idx, int &subdet, int &z, int &lay, int &wafer, int &celltyp, int &cell)
const_iterator begin() const
Definition: Run.h:42
MonitorElement * hebdEtaVsEta
edm::EDGetTokenT< HGChefRecHitCollection > fhRecHitToken_