CMS 3D CMS Logo
AnalyzerMinbias.cc
Go to the documentation of this file.
1 // system include files
2 #include <memory>
3 #include <string>
4 #include <iostream>
5 #include <fstream>
6 #include <sstream>
7 #include <vector>
8 #include <map>
9 
10 // user include files
11 #include "CondFormats/DataRecord/interface/HcalRespCorrsRcd.h"
12 #include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"
13 #include "CondFormats/HcalObjects/interface/HcalRespCorrs.h"
14 #include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
15 
16 #include "DataFormats/Common/interface/Ref.h"
17 #include "DataFormats/Common/interface/TriggerResults.h"
18 #include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
19 #include "DataFormats/CaloTowers/interface/CaloTowerDetId.h"
20 #include "DataFormats/Candidate/interface/Candidate.h"
21 #include "DataFormats/DetId/interface/DetId.h"
22 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
23 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
24 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
25 #include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
26 #include "DataFormats/HcalDigi/interface/HcalCalibrationEventTypes.h"
27 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
28 #include "DataFormats/JetReco/interface/Jet.h"
29 #include "DataFormats/JetReco/interface/CaloJet.h"
30 #include "DataFormats/JetReco/interface/CaloJetCollection.h"
31 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"
32 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
33 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
34 #include "DataFormats/L1GlobalTrigger/interface/L1GtfeWord.h"
35 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
36 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerRecord.h"
37 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
38 #include "DataFormats/Provenance/interface/Provenance.h"
39 
40 #include "EventFilter/HcalRawToDigi/interface/HcalDCCHeader.h"
41 
42 #include "FWCore/Common/interface/TriggerNames.h"
43 #include "FWCore/Framework/interface/Frameworkfwd.h"
44 #include "FWCore/Framework/interface/EDAnalyzer.h"
45 #include "FWCore/Framework/interface/ESHandle.h"
46 #include "FWCore/Framework/interface/Event.h"
47 #include "FWCore/Framework/interface/EventSetup.h"
48 #include "FWCore/Framework/interface/Run.h"
49 #include "FWCore/MessageLogger/interface/MessageLogger.h"
50 #include "FWCore/ParameterSet/interface/ParameterSet.h"
51 
52 #include "Geometry/Records/interface/CaloGeometryRecord.h"
53 #include "Geometry/Records/interface/IdealGeometryRecord.h"
54 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
55 #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
56 
57 #include "TFile.h"
58 #include "TH1.h"
59 #include "TH2.h"
60 #include "TTree.h"
61 
62 namespace HcalMinbias {}
63 
64 // constructors and destructor
65 class AnalyzerMinbias : public edm::EDAnalyzer {
66 public:
67  explicit AnalyzerMinbias(const edm::ParameterSet&);
68  ~AnalyzerMinbias() override;
69 
70  void analyze(const edm::Event&, const edm::EventSetup&) override;
71  void beginJob() override;
72  void endJob() override;
73 
74 private:
75  void analyzeHcal(const HcalRespCorrs* myRecalib,
76  const HBHERecHitCollection& HithbheNS,
77  const HBHERecHitCollection& HithbheMB,
78  const HFRecHitCollection& HithfNS,
79  const HFRecHitCollection& HithfMB,
80  int algoBit,
81  bool fill);
82 
83  struct myInfo {
84  double theMB0, theMB1, theMB2, theMB3, theMB4;
85  double theNS0, theNS1, theNS2, theNS3, theNS4;
86  double theDif0, theDif1, theDif2, runcheck;
87  myInfo() {
88  theMB0 = theMB1 = theMB2 = theMB3 = theMB4 = 0;
89  theNS0 = theNS1 = theNS2 = theNS3 = theNS4 = 0;
90  theDif0 = theDif1 = theDif2 = runcheck = 0;
91  }
92  };
93 
94  // ----------member data ---------------------------
95  std::string fOutputFileName, hcalfile_;
96  std::ofstream* myout_hcal;
97  TFile* hOutputFile;
98  TTree* myTree;
99  TH1D *h_Noise[4], *h_Signal[4];
100  bool runNZS_, theRecalib_, ignoreL1_;
101  double rnnum;
102 
103  // Root tree members
104  double rnnumber;
105  int mydet, mysubd, depth, iphi, ieta, cells, trigbit;
106  float phi, eta;
107  float mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, occup;
108  float mom0_Noise, mom1_Noise, mom2_Noise, mom3_Noise, mom4_Noise;
109  float mom0_Diff, mom1_Diff, mom2_Diff, mom3_Diff, mom4_Diff;
110 
111  std::map<std::pair<int, HcalDetId>, myInfo> myMap_;
112  edm::EDGetTokenT<HBHERecHitCollection> tok_hbherecoMB_, tok_hbherecoNoise_;
113  edm::EDGetTokenT<HFRecHitCollection> tok_hfrecoMB_, tok_hfrecoNoise_;
114  edm::EDGetTokenT<HORecHitCollection> tok_horecoMB_, tok_horecoNoise_;
115  edm::EDGetTokenT<HBHERecHitCollection> tok_hbheNormal_;
116  edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> tok_hltL1GtMap_;
117 };
118 
119 AnalyzerMinbias::AnalyzerMinbias(const edm::ParameterSet& iConfig) {
120  // get name of output file with histogramms
121  fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile");
122 
123  // get token names of modules, producing object collections
124  tok_hbherecoMB_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"));
125  tok_horecoMB_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"));
126  tok_hfrecoMB_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"));
127 
128  tok_hbherecoNoise_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"));
129  tok_horecoNoise_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"));
130  tok_hfrecoNoise_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"));
131 
132  theRecalib_ = iConfig.getParameter<bool>("Recalib");
133  ignoreL1_ = iConfig.getUntrackedParameter<bool>("IgnoreL1", true);
134  runNZS_ = iConfig.getUntrackedParameter<bool>("RunNZS", true);
135 
136  tok_hbheNormal_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
137  tok_hltL1GtMap_ = consumes<L1GlobalTriggerObjectMapRecord>(edm::InputTag("hltL1GtObjectMap"));
138 }
139 
140 AnalyzerMinbias::~AnalyzerMinbias() {}
141 
142 void AnalyzerMinbias::beginJob() {
143  std::string det[4] = {"HB", "HE", "HO", "HF"};
144  char name[80], title[80];
145  for (int subd = 0; subd < 4; ++subd) {
146  sprintf(name, "Noise_%s", det[subd].c_str());
147  sprintf(title, "Energy Distribution for Noise in %s", det[subd].c_str());
148  h_Noise[subd] = new TH1D(name, title, 100, -10., 10.);
149  sprintf(name, "Signal_%s", det[subd].c_str());
150  sprintf(title, "Energy Distribution for Signal in %s", det[subd].c_str());
151  h_Signal[subd] = new TH1D(name, title, 100, -10., 10.);
152  }
153 
154  hOutputFile = new TFile(fOutputFileName.c_str(), "RECREATE");
155  myTree = new TTree("RecJet", "RecJet Tree");
156  myTree->Branch("mydet", &mydet, "mydet/I");
157  myTree->Branch("mysubd", &mysubd, "mysubd/I");
158  myTree->Branch("cells", &cells, "cells");
159  myTree->Branch("depth", &depth, "depth/I");
160  myTree->Branch("ieta", &ieta, "ieta/I");
161  myTree->Branch("iphi", &iphi, "iphi/I");
162  myTree->Branch("eta", &eta, "eta/F");
163  myTree->Branch("phi", &phi, "phi/F");
164  myTree->Branch("mom0_MB", &mom0_MB, "mom0_MB/F");
165  myTree->Branch("mom1_MB", &mom1_MB, "mom1_MB/F");
166  myTree->Branch("mom2_MB", &mom2_MB, "mom2_MB/F");
167  myTree->Branch("mom3_MB", &mom3_MB, "mom3_MB/F");
168  myTree->Branch("mom4_MB", &mom4_MB, "mom4_MB/F");
169  myTree->Branch("mom0_Noise", &mom0_Noise, "mom0_Noise/F");
170  myTree->Branch("mom1_Noise", &mom1_Noise, "mom1_Noise/F");
171  myTree->Branch("mom2_Noise", &mom2_Noise, "mom2_Noise/F");
172  myTree->Branch("mom3_Noise", &mom3_Noise, "mom3_Noise/F");
173  myTree->Branch("mom4_Noise", &mom4_Noise, "mom4_Noise/F");
174  myTree->Branch("mom0_Diff", &mom0_Diff, "mom0_Diff/F");
175  myTree->Branch("mom1_Diff", &mom1_Diff, "mom1_Diff/F");
176  myTree->Branch("mom2_Diff", &mom2_Diff, "mom2_Diff/F");
177  myTree->Branch("occup", &occup, "occup/F");
178  myTree->Branch("trigbit", &trigbit, "trigbit/I");
179  myTree->Branch("rnnumber", &rnnumber, "rnnumber/D");
180 
181  myMap_.clear();
182 }
183 
184 // EndJob
185 //
186 void AnalyzerMinbias::endJob() {
187  int ii = 0;
188  for (std::map<std::pair<int, HcalDetId>, myInfo>::const_iterator itr = myMap_.begin(); itr != myMap_.end(); ++itr) {
189  LogDebug("AnalyzerMB") << "Fired trigger bit number " << itr->first.first;
190  myInfo info = itr->second;
191  if (info.theMB0 > 0) {
192  mom0_MB = info.theMB0;
193  mom1_MB = info.theMB1;
194  mom2_MB = info.theMB2;
195  mom3_MB = info.theMB3;
196  mom4_MB = info.theMB4;
197  mom0_Noise = info.theNS0;
198  mom1_Noise = info.theNS1;
199  mom2_Noise = info.theNS2;
200  mom3_Noise = info.theNS3;
201  mom4_Noise = info.theNS4;
202  mom0_Diff = info.theDif0;
203  mom1_Diff = info.theDif1;
204  mom2_Diff = info.theDif2;
205  rnnumber = info.runcheck;
206  trigbit = itr->first.first;
207  mysubd = itr->first.second.subdet();
208  depth = itr->first.second.depth();
209  ieta = itr->first.second.ieta();
210  iphi = itr->first.second.iphi();
211 
212  LogDebug("AnalyzerMB") << " Result= " << trigbit << " " << mysubd << " " << ieta << " " << iphi << " mom0 "
213  << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB << " mom3 " << mom3_MB << " mom4 "
214  << mom4_MB << " mom0_Noise " << mom0_Noise << " mom1_Noise " << mom1_Noise
215  << " mom2_Noise " << mom2_Noise << " mom3_Noise " << mom3_Noise << " mom4_Noise "
216  << mom4_Noise << " mom0_Diff " << mom0_Diff << " mom1_Diff " << mom1_Diff << " mom2_Diff "
217  << mom2_Diff;
218  myTree->Fill();
219  ii++;
220  }
221  }
222  cells = ii;
223  LogDebug("AnalyzerMB") << "cells"
224  << " " << cells;
225  hOutputFile->Write();
226  hOutputFile->cd();
227  myTree->Write();
228  for (int i = 0; i < 4; i++) {
229  h_Noise[i]->Write();
230  h_Signal[i]->Write();
231  }
232  hOutputFile->Close();
233 }
234 
235 //
236 // member functions
237 //
238 
239 // ------------ method called to produce the data ------------
240 
241 void AnalyzerMinbias::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
242  rnnum = (float)iEvent.run();
243  const HcalRespCorrs* myRecalib = nullptr;
244  if (theRecalib_) {
245  edm::ESHandle<HcalRespCorrs> recalibCorrs;
246  iSetup.get<HcalRespCorrsRcd>().get("recalibrate", recalibCorrs);
247  myRecalib = recalibCorrs.product();
248  } // theRecalib
249 
250  edm::Handle<HBHERecHitCollection> hbheNormal;
251  iEvent.getByToken(tok_hbheNormal_, hbheNormal);
252  if (!hbheNormal.isValid()) {
253  edm::LogInfo("AnalyzerMB") << " hbheNormal failed";
254  } else {
255  edm::LogInfo("AnalyzerMB") << " The size of the normal collection " << hbheNormal->size();
256  }
257 
258  edm::Handle<HBHERecHitCollection> hbheNS;
259  iEvent.getByToken(tok_hbherecoNoise_, hbheNS);
260  if (!hbheNS.isValid()) {
261  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbheNoise product!";
262  return;
263  }
264  const HBHERecHitCollection HithbheNS = *(hbheNS.product());
265  edm::LogInfo("AnalyzerMB") << "HBHE NS size of collection " << HithbheNS.size();
266  if (runNZS_ && HithbheNS.size() != 5184) {
267  edm::LogWarning("AnalyzerMB") << "HBHE NS problem " << rnnum << " size " << HithbheNS.size();
268  return;
269  }
270 
271  edm::Handle<HBHERecHitCollection> hbheMB;
272  iEvent.getByToken(tok_hbherecoMB_, hbheMB);
273  if (!hbheMB.isValid()) {
274  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe product!";
275  return;
276  }
277  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
278  edm::LogInfo("AnalyzerMB") << "HBHE MB size of collection " << HithbheMB.size();
279  if (runNZS_ && HithbheMB.size() != 5184) {
280  edm::LogWarning("AnalyzerMB") << "HBHE problem " << rnnum << " size " << HithbheMB.size();
281  return;
282  }
283 
284  edm::Handle<HFRecHitCollection> hfNS;
285  iEvent.getByToken(tok_hfrecoNoise_, hfNS);
286  if (!hfNS.isValid()) {
287  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hfNoise product!";
288  return;
289  }
290  const HFRecHitCollection HithfNS = *(hfNS.product());
291  edm::LogInfo("AnalyzerMB") << "HF NS size of collection " << HithfNS.size();
292  if (runNZS_ && HithfNS.size() != 1728) {
293  edm::LogWarning("AnalyzerMB") << "HF NS problem " << rnnum << " size " << HithfNS.size();
294  return;
295  }
296 
297  edm::Handle<HFRecHitCollection> hfMB;
298  iEvent.getByToken(tok_hfrecoMB_, hfMB);
299  if (!hfMB.isValid()) {
300  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf product!";
301  return;
302  }
303  const HFRecHitCollection HithfMB = *(hfMB.product());
304  edm::LogInfo("AnalyzerMB") << "HF MB size of collection " << HithfMB.size();
305  if (runNZS_ && HithfMB.size() != 1728) {
306  edm::LogWarning("AnalyzerMB") << "HF problem " << rnnum << " size " << HithfMB.size();
307  return;
308  }
309 
310  if (ignoreL1_) {
311  analyzeHcal(myRecalib, HithbheNS, HithbheMB, HithfNS, HithfMB, 1, true);
312  } else {
313  edm::Handle<L1GlobalTriggerObjectMapRecord> gtObjectMapRecord;
314  iEvent.getByToken(tok_hltL1GtMap_, gtObjectMapRecord);
315  if (gtObjectMapRecord.isValid()) {
316  const std::vector<L1GlobalTriggerObjectMap>& objMapVec = gtObjectMapRecord->gtObjectMap();
317  int ii(0);
318  bool ok(false), fill(true);
319  for (std::vector<L1GlobalTriggerObjectMap>::const_iterator itMap = objMapVec.begin(); itMap != objMapVec.end();
320  ++itMap, ++ii) {
321  bool resultGt = (*itMap).algoGtlResult();
322  if (resultGt == 1) {
323  ok = true;
324  int algoBit = (*itMap).algoBitNumber();
325  analyzeHcal(myRecalib, HithbheNS, HithbheMB, HithfNS, HithfMB, algoBit, fill);
326  fill = false;
327  std::string algoNameStr = (*itMap).algoName();
328  LogDebug("AnalyzerMB") << "Trigger[" << ii << "] " << algoNameStr << " bit " << algoBit << " entered";
329  }
330  }
331  if (!ok)
332  edm::LogInfo("AnalyzerMB") << "No passed L1 Triggers";
333  }
334  }
335 }
336 
337 void AnalyzerMinbias::analyzeHcal(const HcalRespCorrs* myRecalib,
338  const HBHERecHitCollection& HithbheNS,
339  const HBHERecHitCollection& HithbheMB,
340  const HFRecHitCollection& HithfNS,
341  const HFRecHitCollection& HithfMB,
342  int algoBit,
343  bool fill) {
344  // Noise part for HB HE
345  std::map<std::pair<int, HcalDetId>, myInfo> tmpMap;
346  tmpMap.clear();
347 
348  for (HBHERecHitCollection::const_iterator hbheItr = HithbheNS.begin(); hbheItr != HithbheNS.end(); hbheItr++) {
349  // Recalibration of energy
350  float icalconst = 1.;
351  DetId mydetid = hbheItr->id().rawId();
352  if (theRecalib_)
353  icalconst = myRecalib->getValues(mydetid)->getValue();
354 
355  HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
356  double energyhit = aHit.energy();
357 
358  DetId id = (*hbheItr).detid();
359  HcalDetId hid = HcalDetId(id);
360  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
361  if (itr1 == myMap_.end()) {
362  myInfo info;
363  myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
364  itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
365  }
366  itr1->second.theNS0++;
367  itr1->second.theNS1 += energyhit;
368  itr1->second.theNS2 += (energyhit * energyhit);
369  itr1->second.theNS3 += (energyhit * energyhit * energyhit);
370  itr1->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
371  itr1->second.runcheck = rnnum;
372  if (fill)
373  h_Noise[hid.subdet() - 1]->Fill(energyhit);
374 
375  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
376  if (itr2 == tmpMap.end()) {
377  myInfo info;
378  tmpMap[std::pair<int, HcalDetId>(algoBit, hid)] = info;
379  itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
380  }
381  itr2->second.theNS0++;
382  itr2->second.theNS1 += energyhit;
383  itr2->second.theNS2 += (energyhit * energyhit);
384  itr2->second.theNS3 += (energyhit * energyhit * energyhit);
385  itr2->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
386  itr2->second.runcheck = rnnum;
387 
388  } // HBHE_NS
389 
390  // Signal part for HB HE
391 
392  for (HBHERecHitCollection::const_iterator hbheItr = HithbheMB.begin(); hbheItr != HithbheMB.end(); hbheItr++) {
393  // Recalibration of energy
394  float icalconst = 1.;
395  DetId mydetid = hbheItr->id().rawId();
396  if (theRecalib_)
397  icalconst = myRecalib->getValues(mydetid)->getValue();
398 
399  HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
400  double energyhit = aHit.energy();
401 
402  DetId id = (*hbheItr).detid();
403  HcalDetId hid = HcalDetId(id);
404 
405  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
406  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
407 
408  if (itr1 == myMap_.end()) {
409  myInfo info;
410  myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
411  itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
412  }
413  itr1->second.theMB0++;
414  itr1->second.theDif0 = 0;
415  itr1->second.theMB1 += energyhit;
416  itr1->second.theMB2 += (energyhit * energyhit);
417  itr1->second.theMB3 += (energyhit * energyhit * energyhit);
418  itr1->second.theMB4 += (energyhit * energyhit * energyhit * energyhit);
419  itr1->second.runcheck = rnnum;
420  float mydiff = 0.0;
421  if (itr2 != tmpMap.end()) {
422  mydiff = energyhit - (itr2->second.theNS1);
423  itr1->second.theDif0++;
424  itr1->second.theDif1 += mydiff;
425  itr1->second.theDif2 += (mydiff * mydiff);
426  if (fill)
427  h_Signal[hid.subdet() - 1]->Fill(mydiff);
428  }
429  } // HBHE_MB
430 
431  // HF
432 
433  for (HFRecHitCollection::const_iterator hbheItr = HithfNS.begin(); hbheItr != HithfNS.end(); hbheItr++) {
434  // Recalibration of energy
435  float icalconst = 1.;
436  DetId mydetid = hbheItr->id().rawId();
437  if (theRecalib_)
438  icalconst = myRecalib->getValues(mydetid)->getValue();
439 
440  HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
441  double energyhit = aHit.energy();
442  // Remove PMT hits
443  if (fabs(energyhit) > 40.)
444  continue;
445  DetId id = (*hbheItr).detid();
446  HcalDetId hid = HcalDetId(id);
447 
448  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
449 
450  if (itr1 == myMap_.end()) {
451  myInfo info;
452  myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
453  itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
454  }
455  itr1->second.theNS0++;
456  itr1->second.theNS1 += energyhit;
457  itr1->second.theNS2 += (energyhit * energyhit);
458  itr1->second.theNS3 += (energyhit * energyhit * energyhit);
459  itr1->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
460  itr1->second.runcheck = rnnum;
461  if (fill)
462  h_Noise[hid.subdet() - 1]->Fill(energyhit);
463 
464  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
465  if (itr2 == tmpMap.end()) {
466  myInfo info;
467  tmpMap[std::pair<int, HcalDetId>(algoBit, hid)] = info;
468  itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
469  }
470  itr2->second.theNS0++;
471  itr2->second.theNS1 += energyhit;
472  itr2->second.theNS2 += (energyhit * energyhit);
473  itr2->second.theNS3 += (energyhit * energyhit * energyhit);
474  itr2->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
475  itr2->second.runcheck = rnnum;
476 
477  } // HF_NS
478 
479  // Signal part for HF
480 
481  for (HFRecHitCollection::const_iterator hbheItr = HithfMB.begin(); hbheItr != HithfMB.end(); hbheItr++) {
482  // Recalibration of energy
483  float icalconst = 1.;
484  DetId mydetid = hbheItr->id().rawId();
485  if (theRecalib_)
486  icalconst = myRecalib->getValues(mydetid)->getValue();
487  HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
488 
489  double energyhit = aHit.energy();
490  // Remove PMT hits
491  if (fabs(energyhit) > 40.)
492  continue;
493 
494  DetId id = (*hbheItr).detid();
495  HcalDetId hid = HcalDetId(id);
496 
497  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
498  std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
499 
500  if (itr1 == myMap_.end()) {
501  myInfo info;
502  myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
503  itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
504  }
505  itr1->second.theMB0++;
506  itr1->second.theDif0 = 0;
507  itr1->second.theMB1 += energyhit;
508  itr1->second.theMB2 += (energyhit * energyhit);
509  itr1->second.theMB3 += (energyhit * energyhit * energyhit);
510  itr1->second.theMB4 += (energyhit * energyhit * energyhit * energyhit);
511  itr1->second.runcheck = rnnum;
512  float mydiff = 0.0;
513  if (itr2 != tmpMap.end()) {
514  mydiff = energyhit - (itr2->second.theNS1);
515  itr1->second.theDif0++;
516  itr1->second.theDif1 += mydiff;
517  itr1->second.theDif2 += (mydiff * mydiff);
518  if (fill)
519  h_Signal[hid.subdet() - 1]->Fill(mydiff);
520  }
521  }
522 }
523 
524 //define this as a plug-in
525 #include "FWCore/Framework/interface/MakerMacros.h"
526 DEFINE_FWK_MODULE(AnalyzerMinbias);
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
L1GtTriggerMenu.h
CaloJetCollection.h
AnalyzerMinbias::mom2_Noise
float mom2_Noise
Definition: AnalyzerMinbias.cc:108
CaloJet.h
HcalRespCorrsRcd.h
HcalRespCorrs
Definition: HcalRespCorrs.h:17
mps_fire.i
i
Definition: mps_fire.py:355
AnalyzerMinbias::mom4_Noise
float mom4_Noise
Definition: AnalyzerMinbias.cc:108
AnalyzerMinbias::eta
float eta
Definition: AnalyzerMinbias.cc:106
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
MessageLogger.h
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
AnalyzerMinbias::myInfo::theNS2
double theNS2
Definition: AnalyzerMinbias.cc:85
CaloRecHit::energy
constexpr float energy() const
Definition: CaloRecHit.h:29
edm::Handle::product
T const * product() const
Definition: Handle.h:70
AnalyzerMinbias::myInfo::theNS4
double theNS4
Definition: AnalyzerMinbias.cc:85
AnalyzerMinbias::myInfo::theDif2
double theDif2
Definition: AnalyzerMinbias.cc:86
ESHandle.h
TriggerResults.h
HcalCalibrationEventTypes.h
AnalyzerMinbias::mom0_Noise
float mom0_Noise
Definition: AnalyzerMinbias.cc:108
edm::EDGetTokenT
Definition: EDGetToken.h:33
AnalyzerMinbias::mom3_Diff
float mom3_Diff
Definition: AnalyzerMinbias.cc:109
L1GtTriggerMenuRcd.h
HBHERecHit
Definition: HBHERecHit.h:13
AnalyzerMinbias::tok_hbherecoMB_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbherecoMB_
Definition: AnalyzerMinbias.cc:112
AnalyzerMinbias::myInfo::theNS0
double theNS0
Definition: AnalyzerMinbias.cc:85
AnalyzerMinbias::tok_hbheNormal_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNormal_
Definition: AnalyzerMinbias.cc:115
edm::LogInfo
Definition: MessageLogger.h:254
AnalyzerMinbias
Definition: AnalyzerMinbias.cc:65
AnalyzerMinbias::theRecalib_
bool theRecalib_
Definition: AnalyzerMinbias.cc:100
edm::SortedCollection
Definition: SortedCollection.h:49
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:215
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:152
Jet.h
AnalyzerMinbias::analyzeHcal
void analyzeHcal(const HcalRespCorrs *myRecalib, const HBHERecHitCollection &HithbheNS, const HBHERecHitCollection &HithbheMB, const HFRecHitCollection &HithfNS, const HFRecHitCollection &HithfMB, int algoBit, bool fill)
Definition: AnalyzerMinbias.cc:337
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EDAnalyzer.h
convertSQLiteXML.ok
bool ok
Definition: convertSQLiteXML.py:98
edm::Handle
Definition: AssociativeIterator.h:50
AnalyzerMinbias::tok_hbherecoNoise_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbherecoNoise_
Definition: AnalyzerMinbias.cc:112
AnalyzerMinbias::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: AnalyzerMinbias.cc:241
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
EcalRecHitCollections.h
AnalyzerMinbias::tok_hfrecoMB_
edm::EDGetTokenT< HFRecHitCollection > tok_hfrecoMB_
Definition: AnalyzerMinbias.cc:113
edm::EDAnalyzer
Definition: EDAnalyzer.h:29
AnalyzerMinbias::mom3_MB
float mom3_MB
Definition: AnalyzerMinbias.cc:107
AnalyzerMinbias::trigbit
int trigbit
Definition: AnalyzerMinbias.cc:105
DetId
Definition: DetId.h:17
MakerMacros.h
AnalyzerMinbias::rnnum
double rnnum
Definition: AnalyzerMinbias.cc:101
HFRecHit
Definition: HFRecHit.h:11
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
L1GlobalTriggerRecord.h
AnalyzerMinbias::myInfo::runcheck
double runcheck
Definition: AnalyzerMinbias.cc:86
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:262
AnalyzerMinbias::myInfo::theMB2
double theMB2
Definition: AnalyzerMinbias.cc:84
AnalyzerMinbias::AnalyzerMinbias
AnalyzerMinbias(const edm::ParameterSet &)
Definition: AnalyzerMinbias.cc:119
AnalyzerMinbias::tok_horecoNoise_
edm::EDGetTokenT< HORecHitCollection > tok_horecoNoise_
Definition: AnalyzerMinbias.cc:114
Provenance.h
Run.h
edm::ESHandle
Definition: DTSurvey.h:22
AnalyzerMinbias::h_Noise
TH1D * h_Noise[4]
Definition: AnalyzerMinbias.cc:99
AnalyzerMinbias::~AnalyzerMinbias
~AnalyzerMinbias() override
Definition: AnalyzerMinbias.cc:140
AnalyzerMinbias::h_Signal
TH1D * h_Signal[4]
Definition: AnalyzerMinbias.cc:99
AnalyzerMinbias::ieta
int ieta
Definition: AnalyzerMinbias.cc:105
AnalyzerMinbias::mom0_Diff
float mom0_Diff
Definition: AnalyzerMinbias.cc:109
AnalyzerMinbias::myInfo::theMB3
double theMB3
Definition: AnalyzerMinbias.cc:84
AnalyzerMinbias::mydet
int mydet
Definition: AnalyzerMinbias.cc:105
AnalyzerMinbias::mom4_Diff
float mom4_Diff
Definition: AnalyzerMinbias.cc:109
CaloGeometryRecord.h
AnalyzerMinbias::myInfo::theDif0
double theDif0
Definition: AnalyzerMinbias.cc:86
AnalyzerMinbias::tok_horecoMB_
edm::EDGetTokenT< HORecHitCollection > tok_horecoMB_
Definition: AnalyzerMinbias.cc:114
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::LogWarning
Definition: MessageLogger.h:141
ntuplemaker.fill
fill
Definition: ntuplemaker.py:304
AnalyzerMinbias::occup
float occup
Definition: AnalyzerMinbias.cc:107
HcalRespCorr::getValue
float getValue() const
Definition: HcalRespCorr.h:19
AnalyzerMinbias::mysubd
int mysubd
Definition: AnalyzerMinbias.cc:105
CaloSubdetectorGeometry.h
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:670
edm::ParameterSet
Definition: ParameterSet.h:36
AnalyzerMinbias::runNZS_
bool runNZS_
Definition: AnalyzerMinbias.cc:100
AnalyzerMinbias::myInfo::theMB1
double theMB1
Definition: AnalyzerMinbias.cc:84
Event.h
AnalyzerMinbias::hOutputFile
TFile * hOutputFile
Definition: AnalyzerMinbias.cc:97
AnalyzerMinbias::myInfo::theMB0
double theMB0
Definition: AnalyzerMinbias.cc:84
AnalyzerMinbias::iphi
int iphi
Definition: AnalyzerMinbias.cc:105
AnalyzerMinbias::myMap_
std::map< std::pair< int, HcalDetId >, myInfo > myMap_
Definition: AnalyzerMinbias.cc:111
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:267
HcalDetId.h
AnalyzerMinbias::mom2_MB
float mom2_MB
Definition: AnalyzerMinbias.cc:107
AnalyzerMinbias::endJob
void endJob() override
Definition: AnalyzerMinbias.cc:186
HcalDetId::subdet
constexpr HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:138
AnalyzerMinbias::tok_hltL1GtMap_
edm::EDGetTokenT< L1GlobalTriggerObjectMapRecord > tok_hltL1GtMap_
Definition: AnalyzerMinbias.cc:116
HcalDetId
Definition: HcalDetId.h:12
TriggerNames.h
AnalyzerMinbias::beginJob
void beginJob() override
Definition: AnalyzerMinbias.cc:142
iEvent
int iEvent
Definition: GenABIO.cc:224
AnalyzerMinbias::myInfo::myInfo
myInfo()
Definition: AnalyzerMinbias.cc:87
AnalyzerMinbias::phi
float phi
Definition: AnalyzerMinbias.cc:106
IdealGeometryRecord.h
CaloTowerDetId.h
edm::EventSetup
Definition: EventSetup.h:57
HcalSubdetector.h
HcalRespCorrs.h
HcalRespCorrsRcd
Definition: HcalRespCorrsRcd.h:27
AnalyzerMinbias::fOutputFileName
std::string fOutputFileName
Definition: AnalyzerMinbias.cc:95
AnalyzerMinbias::mom1_Diff
float mom1_Diff
Definition: AnalyzerMinbias.cc:109
itr
std::vector< std::pair< float, float > >::iterator itr
Definition: HGCDigitizer.cc:28
get
#define get
AnalyzerMinbias::myInfo::theMB4
double theMB4
Definition: AnalyzerMinbias.cc:84
AnalyzerMinbias::cells
int cells
Definition: AnalyzerMinbias.cc:105
L1GtfeWord.h
CaloTowerCollection.h
overlapproblemtsosanalyzer_cfi.title
title
Definition: overlapproblemtsosanalyzer_cfi.py:7
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AnalyzerMinbias::myInfo::theNS3
double theNS3
Definition: AnalyzerMinbias.cc:85
CaloCellGeometry.h
HcalMinbias
Definition: AnalyzerMinbias.cc:62
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
Ref.h
AnalyzerMinbias::mom1_MB
float mom1_MB
Definition: AnalyzerMinbias.cc:107
DetId.h
Frameworkfwd.h
AnalyzerMinbias::mom2_Diff
float mom2_Diff
Definition: AnalyzerMinbias.cc:109
AnalyzerMinbias::myTree
TTree * myTree
Definition: AnalyzerMinbias.cc:98
AnalyzerMinbias::mom1_Noise
float mom1_Noise
Definition: AnalyzerMinbias.cc:108
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EventSetup.h
L1GlobalTriggerObjectMapRecord::gtObjectMap
const std::vector< L1GlobalTriggerObjectMap > & gtObjectMap() const
get / set the vector of object maps
Definition: L1GlobalTriggerObjectMapRecord.h:64
AnalyzerMinbias::myInfo::theNS1
double theNS1
Definition: AnalyzerMinbias.cc:85
HcalRecHitCollections.h
AnalyzerMinbias::depth
int depth
Definition: AnalyzerMinbias.cc:105
AnalyzerMinbias::rnnumber
double rnnumber
Definition: AnalyzerMinbias.cc:104
AnalyzerMinbias::hcalfile_
std::string hcalfile_
Definition: AnalyzerMinbias.cc:95
AnalyzerMinbias::myInfo::theDif1
double theDif1
Definition: AnalyzerMinbias.cc:86
HcalDCCHeader.h
Candidate.h
AnalyzerMinbias::mom3_Noise
float mom3_Noise
Definition: AnalyzerMinbias.cc:108
genParticles_cff.map
map
Definition: genParticles_cff.py:11
L1GlobalTriggerObjectMapRecord.h
AnalyzerMinbias::mom4_MB
float mom4_MB
Definition: AnalyzerMinbias.cc:107
ParameterSet.h
AnalyzerMinbias::myout_hcal
std::ofstream * myout_hcal
Definition: AnalyzerMinbias.cc:96
AnalyzerMinbias::ignoreL1_
bool ignoreL1_
Definition: AnalyzerMinbias.cc:100
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
edm::Event
Definition: Event.h:73
L1GlobalTriggerObjectMap.h
AnalyzerMinbias::myInfo
Definition: AnalyzerMinbias.cc:83
AnalyzerMinbias::mom0_MB
float mom0_MB
Definition: AnalyzerMinbias.cc:107
GlobalPoint.h
cuy.ii
ii
Definition: cuy.py:590
edm::InputTag
Definition: InputTag.h:15
L1GlobalTriggerReadoutRecord.h
AnalyzerMinbias::tok_hfrecoNoise_
edm::EDGetTokenT< HFRecHitCollection > tok_hfrecoNoise_
Definition: AnalyzerMinbias.cc:113
Generated for CMSSW Reference Manual by doxygen 1.8.16