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