CMS 3D CMS Logo
Analyzer_minbias.cc
Go to the documentation of this file.
1 // system include files
2 #include <memory>
3 #include <string>
4 #include <iostream>
5 
6 // user include files
7 #include "Geometry/Records/interface/IdealGeometryRecord.h"
8 #include "Calibration/HcalCalibAlgos/plugins/Analyzer_minbias.h"
9 #include "DataFormats/Provenance/interface/StableProvenance.h"
10 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
11 #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
12 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
13 #include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
14 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
15 #include "DataFormats/L1GlobalTrigger/interface/L1GtfeWord.h"
16 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
17 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerRecord.h"
18 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
19 #include "DataFormats/FEDRawData/interface/FEDRawData.h"
20 #include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
21 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
22 #include "DataFormats/HcalDigi/interface/HcalCalibrationEventTypes.h"
23 #include "EventFilter/HcalRawToDigi/interface/HcalDCCHeader.h"
24 #include "FWCore/MessageLogger/interface/MessageLogger.h"
25 
26 #include "TFile.h"
27 #include "TH1.h"
28 #include "TH2.h"
29 #include <fstream>
30 #include <sstream>
31 
32 #include "CondFormats/HcalObjects/interface/HcalRespCorrs.h"
33 #include "CondFormats/DataRecord/interface/HcalRespCorrsRcd.h"
34 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"
35 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
36 
37 //
38 // constructors and destructor
39 //
40 namespace cms {
41  Analyzer_minbias::Analyzer_minbias(const edm::ParameterSet& iConfig) {
42  usesResource(TFileService::kSharedResource);
43  // get name of output file with histogramms
44  fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile");
45  // get names of modules, producing object collections
46 
47  tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"));
48  tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"));
49  tok_hf_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"));
50  tok_data_ = consumes<FEDRawDataCollection>(edm::InputTag(iConfig.getParameter<std::string>("InputLabel")));
51 
52  tok_hbheNoise_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"));
53  tok_hoNoise_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"));
54  tok_hfNoise_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"));
55 
56  // this was hardcodded..
57  tok_gtRec_ = consumes<L1GlobalTriggerReadoutRecord>(edm::InputTag("gtDigisAlCaMB"));
58  tok_hbheNorm_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
59 
60  tok_respCorr_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
61  tok_l1gt_ = esConsumes<L1GtTriggerMenu, L1GtTriggerMenuRcd>();
62 
63  theRecalib = iConfig.getParameter<bool>("Recalib");
64 
65  //
66  //
67  for (int i = 0; i < 73; i++) {
68  for (int j = 0; j < 43; j++) {
69  noise_min[i][j] = 0.;
70  noise_pl[i][j] = 0.;
71  }
72  }
73  //
74  //
75  }
76 
77  Analyzer_minbias::~Analyzer_minbias() {
78  // do anything here that needs to be done at desctruction time
79  // (e.g. close files, deallocate resources etc.)
80  }
81 
82  void Analyzer_minbias::beginRun(const edm::Run&, const edm::EventSetup&) { nevent_run = 0; }
83  void Analyzer_minbias::endRun(const edm::Run& r, const edm::EventSetup&) {
84  edm::LogVerbatim("AnalyzerMB") << " Runnumber " << r.run() << " Nevents " << nevent_run;
85  }
86 
87  void Analyzer_minbias::beginJob() {
88  // hOutputFile = new TFile( fOutputFileName.c_str(), "RECREATE" ) ;
89  myTree = fs->make<TTree>("RecJet", "RecJet Tree");
90  myTree->Branch("mydet", &mydet, "mydet/I");
91  myTree->Branch("mysubd", &mysubd, "mysubd/I");
92  myTree->Branch("depth", &depth, "depth/I");
93  myTree->Branch("ieta", &ieta, "ieta/I");
94  myTree->Branch("iphi", &iphi, "iphi/I");
95  myTree->Branch("eta", &eta, "eta/F");
96  myTree->Branch("phi", &phi, "phi/F");
97 
98  myTree->Branch("mom0_MB", &mom0_MB, "mom0_MB/F");
99  myTree->Branch("mom1_MB", &mom1_MB, "mom1_MB/F");
100  myTree->Branch("mom2_MB", &mom2_MB, "mom2_MB/F");
101  myTree->Branch("mom4_MB", &mom4_MB, "mom4_MB/F");
102 
103  myTree->Branch("mom0_Noise", &mom0_Noise, "mom0_Noise/F");
104  myTree->Branch("mom1_Noise", &mom1_Noise, "mom1_Noise/F");
105  myTree->Branch("mom2_Noise", &mom2_Noise, "mom2_Noise/F");
106  myTree->Branch("mom4_Noise", &mom4_Noise, "mom4_Noise/F");
107 
108  myTree->Branch("mom0_Diff", &mom0_Diff, "mom0_Diff/F");
109  myTree->Branch("mom1_Diff", &mom1_Diff, "mom1_Diff/F");
110  myTree->Branch("mom2_Diff", &mom2_Diff, "mom2_Diff/F");
111 
112  myTree->Branch("occup", &occup, "occup/F");
113 
114  edm::LogVerbatim("AnalyzerMB") << " Before ordering Histos ";
115 
116  char str0[32];
117  char str1[32];
118 
119  char str10[32];
120  char str11[32];
121 
122  int k = 0;
123  nevent = 0;
124  // Size of collections
125 
126  hHBHEsize_vs_run =
127  fs->make<TH2F>("hHBHEsize_vs_run", "hHBHEsize_vs_run", 500, 111500., 112000., 6101, -100.5, 6000.5);
128  hHFsize_vs_run = fs->make<TH2F>("hHFsize_vs_run", "hHFsize_vs_run", 500, 111500., 112000., 6101, -100.5, 6000.5);
129 
130  for (int i = 1; i < 73; i++) {
131  for (int j = 1; j < 43; j++) {
132  meannoise_pl[i][j] = 0.;
133  meannoise_min[i][j] = 0.;
134 
135  // for(int l=1;l<5;l++){
136  k = i * 1000 + j;
137  sprintf(str0, "mpl%d", k);
138  sprintf(str1, "mmin%d", k);
139 
140  sprintf(str10, "vpl%d", k);
141  sprintf(str11, "vmin%d", k);
142  // edm::LogVerbatim("AnalyzerMB")<<" "<<i<<" "<<j;
143  if (j < 30) {
144  // first order moment
145  hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
146  hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
147 
148  // second order moment
149  hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 20.);
150  hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 20.);
151  } else {
152  // HF
153  // first order moment
154  // edm::LogVerbatim("AnalyzerMB")<<" "<<i<<" "<<j<<" "<<k;
155  if (j < 40) {
156  hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
157  hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
158  //
159  // second order moment
160  hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 40.);
161  hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 40.);
162  } else {
163  hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
164  hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
165 
166  // second order moment
167  hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 120.);
168  hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 120.);
169  }
170  } // HE/HF boundary
171  // } // l
172  } // j
173  } // i
174 
175  hbheNoiseE = fs->make<TH1F>("hbheNoiseE", "hbheNoiseE", 320, -10., 10.);
176  hfNoiseE = fs->make<TH1F>("hfNoiseE", "hfNoiseE", 320, -10., 10.);
177  hbheSignalE = fs->make<TH1F>("hbheSignalE", "hbheSignalE", 320, -10., 10.);
178  hfSignalE = fs->make<TH1F>("hfSignalE", "hfSignalE", 320, -10., 10.);
179 
180  edm::LogVerbatim("AnalyzerMB") << " After ordering Histos ";
181 
182  std::string ccc = "noise_0.dat";
183 
184  myout_hcal = new std::ofstream(ccc.c_str());
185  if (!myout_hcal)
186  edm::LogVerbatim("AnalyzerMB") << " Output file not open!!! ";
187 
188  //
189  for (int i = 0; i < 5; i++) {
190  for (int j = 0; j < 5; j++) {
191  for (int k = 0; k < 73; k++) {
192  for (int l = 0; l < 43; l++) {
193  theMBFillDetMapPl0[i][j][k][l] = 0.;
194  theMBFillDetMapPl1[i][j][k][l] = 0.;
195  theMBFillDetMapPl2[i][j][k][l] = 0.;
196  theMBFillDetMapPl4[i][j][k][l] = 0.;
197 
198  theMBFillDetMapMin0[i][j][k][l] = 0.;
199  theMBFillDetMapMin1[i][j][k][l] = 0.;
200  theMBFillDetMapMin2[i][j][k][l] = 0.;
201  theMBFillDetMapMin4[i][j][k][l] = 0.;
202 
203  theNSFillDetMapPl0[i][j][k][l] = 0.;
204  theNSFillDetMapPl1[i][j][k][l] = 0.;
205  theNSFillDetMapPl2[i][j][k][l] = 0.;
206  theNSFillDetMapPl4[i][j][k][l] = 0.;
207 
208  theNSFillDetMapMin0[i][j][k][l] = 0.;
209  theNSFillDetMapMin1[i][j][k][l] = 0.;
210  theNSFillDetMapMin2[i][j][k][l] = 0.;
211  theNSFillDetMapMin4[i][j][k][l] = 0.;
212 
213  theDFFillDetMapPl0[i][j][k][l] = 0.;
214  theDFFillDetMapPl1[i][j][k][l] = 0.;
215  theDFFillDetMapPl2[i][j][k][l] = 0.;
216  theDFFillDetMapMin0[i][j][k][l] = 0.;
217  theDFFillDetMapMin1[i][j][k][l] = 0.;
218  theDFFillDetMapMin2[i][j][k][l] = 0.;
219  }
220  }
221  }
222  }
223 
224  return;
225  }
226  //
227  // EndJob
228  //
229  void Analyzer_minbias::endJob() {
230  int ii = 0;
231 
232  for (int i = 1; i < 5; i++) {
233  for (int j = 1; j < 5; j++) {
234  for (int k = 1; k < 73; k++) {
235  for (int l = 1; l < 43; l++) {
236  if (theMBFillDetMapPl0[i][j][k][l] > 0) {
237  mom0_MB = theMBFillDetMapPl0[i][j][k][l];
238  mom1_MB = theMBFillDetMapPl1[i][j][k][l];
239  mom2_MB = theMBFillDetMapPl2[i][j][k][l];
240  mom4_MB = theMBFillDetMapPl4[i][j][k][l];
241  mom0_Noise = theNSFillDetMapPl0[i][j][k][l];
242  mom1_Noise = theNSFillDetMapPl1[i][j][k][l];
243  mom2_Noise = theNSFillDetMapPl2[i][j][k][l];
244  mom4_Noise = theNSFillDetMapPl4[i][j][k][l];
245  mom0_Diff = theDFFillDetMapPl0[i][j][k][l];
246  mom1_Diff = theDFFillDetMapPl1[i][j][k][l];
247  mom2_Diff = theDFFillDetMapPl2[i][j][k][l];
248 
249  mysubd = i;
250  depth = j;
251  ieta = l;
252  iphi = k;
253  edm::LogVerbatim("AnalyzerMB") << " Result Plus= " << mysubd << " " << ieta << " " << iphi << " mom0 "
254  << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB;
255  myTree->Fill();
256  ii++;
257  } // Pl > 0
258 
259  if (theMBFillDetMapMin0[i][j][k][l] > 0) {
260  mom0_MB = theMBFillDetMapMin0[i][j][k][l];
261  mom1_MB = theMBFillDetMapMin1[i][j][k][l];
262  mom2_MB = theMBFillDetMapMin2[i][j][k][l];
263  mom4_MB = theMBFillDetMapMin4[i][j][k][l];
264  mom0_Noise = theNSFillDetMapMin0[i][j][k][l];
265  mom1_Noise = theNSFillDetMapMin1[i][j][k][l];
266  mom2_Noise = theNSFillDetMapMin2[i][j][k][l];
267  mom4_Noise = theNSFillDetMapMin4[i][j][k][l];
268  mom0_Diff = theDFFillDetMapMin0[i][j][k][l];
269  mom1_Diff = theDFFillDetMapMin1[i][j][k][l];
270  mom2_Diff = theDFFillDetMapMin2[i][j][k][l];
271 
272  mysubd = i;
273  depth = j;
274  ieta = -1 * l;
275  iphi = k;
276  edm::LogVerbatim("AnalyzerMB") << " Result Minus= " << mysubd << " " << ieta << " " << iphi << " mom0 "
277  << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB;
278  myTree->Fill();
279  ii++;
280 
281  } // Min>0
282  } // ieta
283  } // iphi
284  } // depth
285  } //subd
286 
287  edm::LogVerbatim("AnalyzerMB") << " Number of cells " << ii;
288 
289  // hOutputFile->Write();
290  // hOutputFile->cd();
291  // myTree->Write();
292  // hHBHEsize_vs_run->Write() ;
293  // hHFsize_vs_run->Write() ;
294 
295  for (int i = 1; i < 73; i++) {
296  for (int j = 1; j < 43; j++) {
297  hCalo1[i][j]->Write();
298  hCalo2[i][j]->Write();
299  hCalo1mom2[i][j]->Write();
300  hCalo2mom2[i][j]->Write();
301  }
302  }
303 
304  // hbheNoiseE->Write() ;
305  // hfNoiseE->Write() ;
306  // hbheSignalE->Write() ;
307  // hfSignalE->Write() ;
308  // hOutputFile->Close() ;
309 
310  edm::LogVerbatim("AnalyzerMB") << " File is closed ";
311 
312  return;
313  }
314 
315  //
316  // member functions
317  //
318 
319  // ------------ method called to produce the data ------------
320  void Analyzer_minbias::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
321  edm::LogVerbatim("AnalyzerMB") << " Start Analyzer_minbias::analyze " << nevent;
322  nevent++;
323  nevent_run++;
324 
325  float rnnum = (float)iEvent.run();
326 
327  std::vector<edm::StableProvenance const*> theProvenance;
328  iEvent.getAllStableProvenance(theProvenance);
329 
330  for (auto const& provenance : theProvenance) {
331  edm::LogVerbatim("AnalyzerMB") << " Print all process/modulelabel/product names " << provenance->processName()
332  << " , " << provenance->moduleLabel() << " , "
333  << provenance->productInstanceName();
334  }
335  /*
336  edm::Handle<FEDRawDataCollection> rawdata;
337  iEvent.getByToken(tok_data_,rawdata);
338 
339  if (!rawdata.isValid()) {
340  edm::LogVerbatim("AnalyzerMB")<<" No valid collection ";
341  } else {
342  edm::LogVerbatim("AnalyzerMB")<<" Valid collection ";
343  int calibType = -1 ; int numEmptyFEDs = 0 ;
344  std::vector<int> calibTypeCounter(8,0) ;
345  for (int i=FEDNumbering::MINHCALFEDID;
346  i<=FEDNumbering::MAXHCALFEDID; i++) {
347  const FEDRawData& fedData = rawdata->FEDData(i) ;
348  if ( fedData.size() < 24 ) numEmptyFEDs++ ;
349  if ( fedData.size() < 24 ) continue ;
350  // int value = ((const HcalDCCHeader*)(fedData.data()))->getCalibType() ;
351  // calibTypeCounter.at(value)++ ; // increment the counter for this calib type
352  }
353  edm::LogVerbatim("AnalyzerMB")<<" NumFed "<<numEmptyFEDs<<" "<<calibType;
354  }
355  */
356  /*
357  std::vector<edm::Handle<FEDRawDataCollection> > rawdata1;
358  iEvent.getManyByType(rawdata1);
359 
360  for(std::vector<edm::Handle<FEDRawDataCollection> >::const_iterator it = rawdata1.begin();it != rawdata1.end(); it++) {
361 
362  edm::LogVerbatim("AnalyzerMB")<<" Many by Type product name "<< (*it).provenance()->processName()<<
363  " "<<(*it).provenance()->moduleLabel();
364 
365  if((*it).provenance()->moduleLabel() == "hltHcalCalibrationRaw") {
366  int calibType = -1 ; int numEmptyFEDs = 0 ;
367 
368  for (int i=FEDNumbering::MINHCALFEDID;
369  i<=FEDNumbering::MAXHCALFEDID; i++) {
370  const FEDRawData& fedData = (*it)->FEDData(i) ;
371  edm::LogVerbatim("AnalyzerMB")<<" FED size "<<fedData.size();
372  if ( fedData.size() < 24 ) numEmptyFEDs++ ;
373  if ( fedData.size() < 24 ) continue ;
374  int value = ((const HcalDCCHeader*)(fedData.data()))->getCalibType() ;
375  edm::LogVerbatim("AnalyzerMB")<<" Value "<<value;
376  }
377  edm::LogVerbatim("AnalyzerMB")<<" Many by Type NumFed "<<numEmptyFEDs<<" "<<calibType;
378  }
379  }
380 
381 
382  */
383 
384  /*
385  const L1GtTriggerMenu* menu = &iSetup.getData(tok_l1gt_);
386  const AlgorithmMap& bitMap = menu->gtAlgorithmMap();
387 
388  edm::Handle<L1GlobalTriggerReadoutRecord> gtRecord;
389  iEvent.getByToken(tok_gtRec_, gtRecord);
390 
391  if (!gtRecord.isValid()) {
392 
393  // LogDebug("L1GlobalTriggerRecordProducer")
394  // << "\n\n Error: no L1GlobalTriggerReadoutRecord found with input tag "
395  // << m_l1GtReadoutRecord
396  // << "\n Returning empty L1GlobalTriggerRecord.\n\n";
397  edm::LogVerbatim("AnalyzerMB")<<" No L1 trigger record ";
398  } else {
399 
400  const DecisionWord dWord = gtRecord->decisionWord();
401 
402  for (CItAlgo itAlgo = bitMap.begin(); itAlgo != bitMap.end(); itAlgo++) {
403  bool decision=menu->gtAlgorithmResult(itAlgo->first,dWord);
404  if(decision == 1) edm::LogVerbatim("AnalyzerMB")<<" Trigger "<<itAlgo->first<<" "<<decision;
405  }
406 
407  }
408  */
409 
410  const HcalRespCorrs* myRecalib = nullptr;
411  if (theRecalib) {
412  myRecalib = &iSetup.getData(tok_respCorr_);
413  } // theRecalib
414 
415  // Noise part for HB HE
416 
417  double tmpNSFillDetMapPl1[5][5][73][43];
418  double tmpNSFillDetMapMin1[5][5][73][43];
419 
420  for (int i = 0; i < 5; i++) {
421  for (int j = 0; j < 5; j++) {
422  for (int k = 0; k < 73; k++) {
423  for (int l = 0; l < 43; l++) {
424  tmpNSFillDetMapPl1[i][j][k][l] = 0.;
425  tmpNSFillDetMapMin1[i][j][k][l] = 0.;
426  }
427  }
428  }
429  }
430 
431  edm::Handle<HBHERecHitCollection> hbheNormal;
432  iEvent.getByToken(tok_hbheNorm_, hbheNormal);
433  if (!hbheNormal.isValid()) {
434  edm::LogWarning("AnalyzerMB") << " hbheNormal failed ";
435  } else {
436  edm::LogVerbatim("AnalyzerMB") << " The size of the normal collection " << hbheNormal->size();
437  }
438 
439  edm::Handle<HBHERecHitCollection> hbheNS;
440  iEvent.getByToken(tok_hbheNoise_, hbheNS);
441 
442  if (!hbheNS.isValid()) {
443  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe"
444  << " product! No HBHE MS ";
445  return;
446  }
447 
448  const HBHERecHitCollection HithbheNS = *(hbheNS.product());
449  edm::LogVerbatim("AnalyzerMB") << " HBHE NS size of collection " << HithbheNS.size();
450  hHBHEsize_vs_run->Fill(rnnum, (float)HithbheNS.size());
451 
452  if (HithbheNS.size() != 5184) {
453  edm::LogWarning("AnalyzerMB") << " HBHE problem " << rnnum << " " << HithbheNS.size();
454  // return;
455  }
456  edm::Handle<HBHERecHitCollection> hbheMB;
457  iEvent.getByToken(tok_hbhe_, hbheMB);
458 
459  if (!hbheMB.isValid()) {
460  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe"
461  << " product! No HBHE MB";
462  // return ;
463  }
464 
465  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
466  edm::LogVerbatim("AnalyzerMB") << " HBHE MB size of collection " << HithbheMB.size();
467  if (HithbheMB.size() != 5184) {
468  edm::LogWarning("AnalyzerMB") << " HBHE problem " << rnnum << " " << HithbheMB.size();
469  // return;
470  }
471 
472  edm::Handle<HFRecHitCollection> hfNS;
473  iEvent.getByToken(tok_hfNoise_, hfNS);
474 
475  if (!hfNS.isValid()) {
476  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf"
477  << " product! No HF NS ";
478  // return ;
479  }
480 
481  const HFRecHitCollection HithfNS = *(hfNS.product());
482  edm::LogVerbatim("AnalyzerMB") << " HFE NS size of collection " << HithfNS.size();
483  hHFsize_vs_run->Fill(rnnum, (float)HithfNS.size());
484  if (HithfNS.size() != 1728) {
485  edm::LogWarning("AnalyzerMB") << " HF problem " << rnnum << " " << HithfNS.size();
486  // return;
487  }
488 
489  edm::Handle<HFRecHitCollection> hfMB;
490  iEvent.getByToken(tok_hf_, hfMB);
491 
492  if (!hfMB.isValid()) {
493  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf"
494  << " product! No HF MB";
495  // return ;
496  }
497 
498  const HFRecHitCollection HithfMB = *(hfMB.product());
499  edm::LogVerbatim("AnalyzerMB") << " HF MB size of collection " << HithfMB.size();
500  if (HithfMB.size() != 1728) {
501  edm::LogWarning("AnalyzerMB") << " HF problem " << rnnum << " " << HithfMB.size();
502  // return;
503  }
504 
505  for (HBHERecHitCollection::const_iterator hbheItr = HithbheNS.begin(); hbheItr != HithbheNS.end(); hbheItr++) {
506  // Recalibration of energy
507  float icalconst = 1.;
508  DetId mydetid = hbheItr->id().rawId();
509  if (theRecalib)
510  icalconst = myRecalib->getValues(mydetid)->getValue();
511 
512  HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
513 
514  double energyhit = aHit.energy();
515 
516  DetId id = (*hbheItr).detid();
517  HcalDetId hid = HcalDetId(id);
518 
519  int mysu = ((hid).rawId() >> 25) & 0x7;
520  if (hid.ieta() > 0) {
521  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
522  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
523  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
524  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
525  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
526  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
527  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
528  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
529 
530  tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
531 
532  } else {
533  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
534  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
535  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
536  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
537  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
538  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
539  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
540  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
541 
542  tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
543  }
544 
545  if (hid.depth() == 1) {
546  hbheNoiseE->Fill(energyhit);
547 
548  if (energyhit < -2.)
549  edm::LogVerbatim("AnalyzerMB") << " Run " << rnnum << " ieta,iphi " << hid.ieta() << " " << hid.iphi()
550  << energyhit;
551 
552  // if( hid.ieta() > 0 ) {
553  // hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit-noise_pl[hid.iphi()][hid.ieta()]);
554  // hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
555  // } else {
556  // hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit-noise_min[hid.iphi()][abs(hid.ieta())]);
557  // hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
558  // } // eta><0
559 
560  } // depth=1
561 
562  } // HBHE_NS
563 
564  // Signal part for HB HE
565 
566  for (HBHERecHitCollection::const_iterator hbheItr = HithbheMB.begin(); hbheItr != HithbheMB.end(); hbheItr++) {
567  // Recalibration of energy
568  float icalconst = 1.;
569  DetId mydetid = hbheItr->id().rawId();
570  if (theRecalib)
571  icalconst = myRecalib->getValues(mydetid)->getValue();
572 
573  HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
574 
575  double energyhit = aHit.energy();
576 
577  DetId id = (*hbheItr).detid();
578  HcalDetId hid = HcalDetId(id);
579 
580  int mysu = ((hid).rawId() >> 25) & 0x7;
581  if (hid.ieta() > 0) {
582  theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] += 1.;
583  theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] += energyhit;
584  theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] += pow(energyhit, 2);
585  theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] += pow(energyhit, 4);
586  float mydiff = energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
587 
588  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
589  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + mydiff;
590  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
591  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(mydiff, 2);
592  } else {
593  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
594  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
595  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
596  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
597  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
598  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
599  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
600  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
601 
602  float mydiff = energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
603  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
604  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + mydiff;
605  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
606  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(mydiff, 2);
607  }
608 
609  if (hid.depth() == 1) {
610  hbheSignalE->Fill(energyhit);
611 
612  if (hid.ieta() > 0) {
613  hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
614  hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit, 2));
615  } else {
616  hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
617  hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit, 2));
618  } // eta><0
619 
620  } // depth=1
621 
622  } // HBHE_MB
623 
624  // HF
625 
626  for (HFRecHitCollection::const_iterator hbheItr = HithfNS.begin(); hbheItr != HithfNS.end(); hbheItr++) {
627  // Recalibration of energy
628  float icalconst = 1.;
629  DetId mydetid = hbheItr->id().rawId();
630  if (theRecalib)
631  icalconst = myRecalib->getValues(mydetid)->getValue();
632 
633  HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
634 
635  double energyhit = aHit.energy();
636  //
637  // Remove PMT hits
638  //
639  DetId id = (*hbheItr).detid();
640  HcalDetId hid = HcalDetId(id);
641 
642  if (fabs(energyhit) > 40.)
643  continue;
644 
645  int mysu = hid.subdetId();
646  if (hid.ieta() > 0) {
647  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
648  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
649  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
650  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
651  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
652  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
653  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
654  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
655 
656  tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
657 
658  } else {
659  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
660  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
661  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
662  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
663  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
664  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
665  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
666  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
667 
668  tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
669  }
670 
671  if (hid.depth() == 1) {
672  hfNoiseE->Fill(energyhit);
673 
674  //if( hid.ieta() > 0 ) {
675  // hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit-noise_pl[hid.iphi()][hid.ieta()]);
676  // hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
677  //} else {
678  // hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit-noise_min[hid.iphi()][abs(hid.ieta())]);
679  // hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
680  //} // eta><0
681 
682  } // depth=1
683 
684  } // HBHE_NS
685 
686  // Signal part for HB HE
687 
688  for (HFRecHitCollection::const_iterator hbheItr = HithfMB.begin(); hbheItr != HithfMB.end(); hbheItr++) {
689  // Recalibration of energy
690  float icalconst = 1.;
691  DetId mydetid = hbheItr->id().rawId();
692  if (theRecalib)
693  icalconst = myRecalib->getValues(mydetid)->getValue();
694 
695  HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
696 
697  double energyhit = aHit.energy();
698  //
699  // Remove PMT hits
700  //
701  if (fabs(energyhit) > 40.)
702  continue;
703 
704  DetId id = (*hbheItr).detid();
705  HcalDetId hid = HcalDetId(id);
706 
707  int mysu = ((hid).rawId() >> 25) & 0x7;
708  if (hid.ieta() > 0) {
709  theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
710  theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
711  theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
712  theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
713  theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
714  theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
715  theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
716  theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
717 
718  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
719  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit -
720  tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
721  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
722  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] +
723  pow((energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]), 2);
724  } else {
725  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
726  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
727  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
728  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
729  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
730  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
731  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
732  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
733 
734  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
735  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit -
736  tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
737  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
738  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] +
739  pow((energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()]), 2);
740  }
741 
742  if (hid.depth() == 1) {
743  hfSignalE->Fill(energyhit);
744 
745  if (hid.ieta() > 0) {
746  hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
747  hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit, 2));
748  } else {
749  hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
750  hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit, 2));
751  } // eta><0
752 
753  } // depth=1
754 
755  } // HF_MB
756 
757  edm::LogVerbatim("AnalyzerMB") << " Event is finished ";
758  }
759 } // namespace cms
760 
761 #include "FWCore/Framework/interface/MakerMacros.h"
762 
763 using cms::Analyzer_minbias;
764 
765 DEFINE_FWK_MODULE(Analyzer_minbias);
cms::Analyzer_minbias::noise_min
double noise_min[73][43]
Definition: Analyzer_minbias.h:108
cms::Analyzer_minbias::endJob
void endJob() override
Definition: Analyzer_minbias.cc:229
cms::Analyzer_minbias::theDFFillDetMapMin1
double theDFFillDetMapMin1[5][5][73][43]
Definition: Analyzer_minbias.h:137
FEDNumbering.h
cms::Analyzer_minbias::theMBFillDetMapMin4
double theMBFillDetMapMin4[5][5][73][43]
Definition: Analyzer_minbias.h:121
HcalRespCorrsRcd.h
HcalRespCorrs
Definition: HcalRespCorrs.h:17
mps_fire.i
i
Definition: mps_fire.py:428
cms::Analyzer_minbias::theMBFillDetMapPl1
double theMBFillDetMapPl1[5][5][73][43]
Definition: Analyzer_minbias.h:114
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
MessageLogger.h
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
cms::Analyzer_minbias::theDFFillDetMapMin0
double theDFFillDetMapMin0[5][5][73][43]
Definition: Analyzer_minbias.h:136
cms::Analyzer_minbias::hCalo1mom2
TH1F * hCalo1mom2[73][43]
Definition: Analyzer_minbias.h:88
cms::Analyzer_minbias::beginJob
void beginJob() override
Definition: Analyzer_minbias.cc:87
CaloRecHit::energy
constexpr float energy() const
Definition: CaloRecHit.h:29
edm::Handle::product
T const * product() const
Definition: Handle.h:70
cms::Analyzer_minbias::phi
float phi
Definition: Analyzer_minbias.h:100
cms::Analyzer_minbias::tok_hoNoise_
edm::EDGetTokenT< HORecHitCollection > tok_hoNoise_
Definition: Analyzer_minbias.h:145
cms::Analyzer_minbias::mom0_Diff
float mom0_Diff
Definition: Analyzer_minbias.h:103
cms::Analyzer_minbias::theMBFillDetMapMin0
double theMBFillDetMapMin0[5][5][73][43]
Definition: Analyzer_minbias.h:118
cms::Analyzer_minbias::meannoise_pl
double meannoise_pl[73][43]
Definition: Analyzer_minbias.h:107
HcalDetId::iphi
constexpr int iphi() const
get the cell iphi
Definition: HcalDetId.h:157
HcalCalibrationEventTypes.h
cms::Analyzer_minbias::tok_hbheNorm_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNorm_
Definition: Analyzer_minbias.h:150
cms::Analyzer_minbias::myout_hcal
std::ofstream * myout_hcal
Definition: Analyzer_minbias.h:77
edm::Run
Definition: Run.h:45
cms::Analyzer_minbias::mom4_MB
float mom4_MB
Definition: Analyzer_minbias.h:101
HBHERecHit
Definition: HBHERecHit.h:13
cms::Analyzer_minbias::theMBFillDetMapPl4
double theMBFillDetMapPl4[5][5][73][43]
Definition: Analyzer_minbias.h:116
cms::Analyzer_minbias::mom2_Noise
float mom2_Noise
Definition: Analyzer_minbias.h:102
edm::SortedCollection
Definition: SortedCollection.h:49
cms::Analyzer_minbias::mom4_Noise
float mom4_Noise
Definition: Analyzer_minbias.h:102
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:215
cms::Analyzer_minbias::fs
edm::Service< TFileService > fs
Definition: Analyzer_minbias.h:82
HcalDetId::depth
constexpr int depth() const
get the tower depth
Definition: HcalDetId.h:164
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
cms::Analyzer_minbias::hbheSignalE
TH1F * hbheSignalE
Definition: Analyzer_minbias.h:91
cms::Analyzer_minbias::noise_pl
double noise_pl[73][43]
Definition: Analyzer_minbias.h:108
cms::Analyzer_minbias::hHFsize_vs_run
TH2F * hHFsize_vs_run
Definition: Analyzer_minbias.h:96
FEDRawData.h
edm::Handle
Definition: AssociativeIterator.h:50
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
cms::Analyzer_minbias::tok_gtRec_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_gtRec_
Definition: Analyzer_minbias.h:149
cms::Analyzer_minbias::occup
float occup
Definition: Analyzer_minbias.h:101
cms::Analyzer_minbias::theMBFillDetMapMin1
double theMBFillDetMapMin1[5][5][73][43]
Definition: Analyzer_minbias.h:119
cms::Analyzer_minbias::analyze
void analyze(edm::Event const &, edm::EventSetup const &) override
Definition: Analyzer_minbias.cc:320
cms::Analyzer_minbias::hCalo2
TH1F * hCalo2[73][43]
Definition: Analyzer_minbias.h:87
cms::Analyzer_minbias::theNSFillDetMapPl1
double theNSFillDetMapPl1[5][5][73][43]
Definition: Analyzer_minbias.h:124
DetId
Definition: DetId.h:17
cms::Analyzer_minbias::tok_ho_
edm::EDGetTokenT< HORecHitCollection > tok_ho_
Definition: Analyzer_minbias.h:141
MakerMacros.h
cms::Analyzer_minbias::beginRun
void beginRun(edm::Run const &, edm::EventSetup const &) override
Definition: Analyzer_minbias.cc:82
cms::Analyzer_minbias::nevent
double nevent
Definition: Analyzer_minbias.h:112
cms::Analyzer_minbias::hCalo2mom2
TH1F * hCalo2mom2[73][43]
Definition: Analyzer_minbias.h:89
Analyzer_minbias.h
cms::Analyzer_minbias::Analyzer_minbias
Analyzer_minbias(const edm::ParameterSet &)
Definition: Analyzer_minbias.cc:41
cms::Analyzer_minbias::mom0_Noise
float mom0_Noise
Definition: Analyzer_minbias.h:102
cms::Analyzer_minbias::hCalo1
TH1F * hCalo1[73][43]
Definition: Analyzer_minbias.h:86
HFRecHit
Definition: HFRecHit.h:11
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
L1GlobalTriggerRecord.h
cms::Analyzer_minbias::theMBFillDetMapPl0
double theMBFillDetMapPl0[5][5][73][43]
Definition: Analyzer_minbias.h:113
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:262
cms::Analyzer_minbias::theNSFillDetMapMin1
double theNSFillDetMapMin1[5][5][73][43]
Definition: Analyzer_minbias.h:129
cms::Analyzer_minbias::nevent_run
int nevent_run
Definition: Analyzer_minbias.h:98
cms::Analyzer_minbias::eta
float eta
Definition: Analyzer_minbias.h:100
cms::Analyzer_minbias::theDFFillDetMapMin2
double theDFFillDetMapMin2[5][5][73][43]
Definition: Analyzer_minbias.h:138
cms::Analyzer_minbias::theNSFillDetMapPl2
double theNSFillDetMapPl2[5][5][73][43]
Definition: Analyzer_minbias.h:125
cms::Analyzer_minbias::~Analyzer_minbias
~Analyzer_minbias() override
Definition: Analyzer_minbias.cc:77
cms::Analyzer_minbias::myTree
TTree * myTree
Definition: Analyzer_minbias.h:85
dqmdumpme.k
k
Definition: dqmdumpme.py:60
cms::Analyzer_minbias::theNSFillDetMapMin0
double theNSFillDetMapMin0[5][5][73][43]
Definition: Analyzer_minbias.h:128
cms::Analyzer_minbias::theNSFillDetMapMin2
double theNSFillDetMapMin2[5][5][73][43]
Definition: Analyzer_minbias.h:130
cms::Analyzer_minbias::theNSFillDetMapMin4
double theNSFillDetMapMin4[5][5][73][43]
Definition: Analyzer_minbias.h:131
cms::Analyzer_minbias::mysubd
int mysubd
Definition: Analyzer_minbias.h:99
cms::Analyzer_minbias::mydet
int mydet
Definition: Analyzer_minbias.h:99
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
cms::Analyzer_minbias::theDFFillDetMapPl0
double theDFFillDetMapPl0[5][5][73][43]
Definition: Analyzer_minbias.h:133
HcalRespCorr::getValue
float getValue() const
Definition: HcalRespCorr.h:19
CaloSubdetectorGeometry.h
cms::Analyzer_minbias::meannoise_min
double meannoise_min[73][43]
Definition: Analyzer_minbias.h:107
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
edm::ParameterSet
Definition: ParameterSet.h:47
cms::Analyzer_minbias::theMBFillDetMapPl2
double theMBFillDetMapPl2[5][5][73][43]
Definition: Analyzer_minbias.h:115
cms::Analyzer_minbias::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: Analyzer_minbias.h:140
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:267
HcalDetId.h
cms::Analyzer_minbias::ieta
int ieta
Definition: Analyzer_minbias.h:99
HcalDetId
Definition: HcalDetId.h:12
iEvent
int iEvent
Definition: GenABIO.cc:224
cms::Analyzer_minbias::hbheNoiseE
TH1F * hbheNoiseE
Definition: Analyzer_minbias.h:90
FEDRawDataCollection.h
cms::Analyzer_minbias
Definition: Analyzer_minbias.h:62
IdealGeometryRecord.h
edm::EventSetup
Definition: EventSetup.h:58
HcalSubdetector.h
HcalRespCorrs.h
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:1053
L1GtfeWord.h
cms::Analyzer_minbias::depth
int depth
Definition: Analyzer_minbias.h:99
cmsLHEtoEOSManager.l
l
Definition: cmsLHEtoEOSManager.py:204
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:127
alignCSCRings.r
r
Definition: alignCSCRings.py:93
CaloCellGeometry.h
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
cms::Analyzer_minbias::hfSignalE
TH1F * hfSignalE
Definition: Analyzer_minbias.h:93
cms::Analyzer_minbias::endRun
void endRun(edm::Run const &, edm::EventSetup const &) override
Definition: Analyzer_minbias.cc:83
cms::Analyzer_minbias::mom2_MB
float mom2_MB
Definition: Analyzer_minbias.h:101
cms::Analyzer_minbias::theNSFillDetMapPl0
double theNSFillDetMapPl0[5][5][73][43]
Definition: Analyzer_minbias.h:123
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
cms::Analyzer_minbias::tok_data_
edm::EDGetTokenT< FEDRawDataCollection > tok_data_
Definition: Analyzer_minbias.h:79
CaloGeometry.h
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
cms::Analyzer_minbias::tok_respCorr_
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respCorr_
Definition: Analyzer_minbias.h:152
cms::Analyzer_minbias::mom1_MB
float mom1_MB
Definition: Analyzer_minbias.h:101
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
HcalDCCHeader.h
cms::Analyzer_minbias::tok_hfNoise_
edm::EDGetTokenT< HFRecHitCollection > tok_hfNoise_
Definition: Analyzer_minbias.h:146
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
cms::Analyzer_minbias::fOutputFileName
std::string fOutputFileName
Definition: Analyzer_minbias.h:75
cms::Analyzer_minbias::theDFFillDetMapPl1
double theDFFillDetMapPl1[5][5][73][43]
Definition: Analyzer_minbias.h:134
L1GlobalTriggerObjectMapRecord.h
cms::Analyzer_minbias::mom2_Diff
float mom2_Diff
Definition: Analyzer_minbias.h:103
cms::Analyzer_minbias::theRecalib
bool theRecalib
Definition: Analyzer_minbias.h:155
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
cms::Analyzer_minbias::hHBHEsize_vs_run
TH2F * hHBHEsize_vs_run
Definition: Analyzer_minbias.h:95
cms::Analyzer_minbias::tok_hbheNoise_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
Definition: Analyzer_minbias.h:144
cms::Analyzer_minbias::theMBFillDetMapMin2
double theMBFillDetMapMin2[5][5][73][43]
Definition: Analyzer_minbias.h:120
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
edm::Event
Definition: Event.h:73
L1GlobalTriggerObjectMap.h
cms::Analyzer_minbias::tok_hf_
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
Definition: Analyzer_minbias.h:142
cms::Analyzer_minbias::theNSFillDetMapPl4
double theNSFillDetMapPl4[5][5][73][43]
Definition: Analyzer_minbias.h:126
cuy.ii
ii
Definition: cuy.py:589
cms::Analyzer_minbias::hfNoiseE
TH1F * hfNoiseE
Definition: Analyzer_minbias.h:92
StableProvenance.h
edm::InputTag
Definition: InputTag.h:15
L1GlobalTriggerReadoutRecord.h
cms::Analyzer_minbias::mom0_MB
float mom0_MB
Definition: Analyzer_minbias.h:101
cms::Analyzer_minbias::mom1_Noise
float mom1_Noise
Definition: Analyzer_minbias.h:102
cms::Analyzer_minbias::mom1_Diff
float mom1_Diff
Definition: Analyzer_minbias.h:103
cms::Analyzer_minbias::theDFFillDetMapPl2
double theDFFillDetMapPl2[5][5][73][43]
Definition: Analyzer_minbias.h:135
cms
Namespace of DDCMS conversion namespace.
Definition: ProducerAnalyzer.cc:21
cms::Analyzer_minbias::iphi
int iphi
Definition: Analyzer_minbias.h:99
cms::Analyzer_minbias::tok_l1gt_
edm::ESGetToken< L1GtTriggerMenu, L1GtTriggerMenuRcd > tok_l1gt_
Definition: Analyzer_minbias.h:153
Generated for CMSSW Reference Manual by doxygen 1.8.16