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