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