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