CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
GlobalRecHitsHistogrammer.cc
Go to the documentation of this file.
1 
10 using namespace std;
11 #include "Validation/GlobalRecHits/interface/GlobalRecHitsHistogrammer.h"
12 #include "DQMServices/Core/interface/DQMStore.h"
13 #include "DQMServices/Core/interface/MonitorElement.h"
14 
15 GlobalRecHitsHistogrammer::GlobalRecHitsHistogrammer(const edm::ParameterSet& iPSet) :
16  fName(""), verbosity(0), frequency(0), label(""), getAllProvenances(false),
17  printProvenanceInfo(false), count(0)
18 {
19  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_GlobalRecHitsAnalyzer";
20 
21  // get information from parameter set
22  fName = iPSet.getUntrackedParameter<std::string>("Name");
23  verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
24  frequency = iPSet.getUntrackedParameter<int>("Frequency");
25  outputfile = iPSet.getParameter<std::string>("outputFile");
26  doOutput = iPSet.getParameter<bool>("DoOutput");
27  edm::ParameterSet m_Prov =
28  iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
29  getAllProvenances =
30  m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
31  printProvenanceInfo =
32  m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");
33 
34  //get Labels to use to extract information
35  GlobalRecHitSrc_ = iPSet.getParameter<edm::InputTag>("GlobalRecHitSrc");
36  // ECalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalEBSrc");
37  //ECalUncalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalUncalEBSrc");
38  //ECalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalEESrc");
39  //ECalUncalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalUncalEESrc");
40  //ECalESSrc_ = iPSet.getParameter<edm::InputTag>("ECalESSrc");
41  //HCalSrc_ = iPSet.getParameter<edm::InputTag>("HCalSrc");
42  //SiStripSrc_ = iPSet.getParameter<edm::InputTag>("SiStripSrc");
43  //SiPxlSrc_ = iPSet.getParameter<edm::InputTag>("SiPxlSrc");
44  //MuDTSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSrc");
45  //MuDTSimSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSimSrc");
46  //MuCSCSrc_ = iPSet.getParameter<edm::InputTag>("MuCSCSrc");
47  //MuRPCSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSrc");
48  //MuRPCSimSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSimSrc");
49 
50  //conf_ = iPSet;
51 
52  // use value of first digit to determine default output level (inclusive)
53  // 0 is none, 1 is basic, 2 is fill output, 3 is gather output
54  verbosity %= 10;
55 
56  // create persistent object
57  // produces<PGlobalRecHit>(label);
58 
59  // print out Parameter Set information being used
60  if (verbosity >= 0) {
61  edm::LogInfo(MsgLoggerCat)
62  << "\n===============================\n"
63  << "Initialized as EDProducer with parameter values:\n"
64  << " Name = " << fName << "\n"
65  << " Verbosity = " << verbosity << "\n"
66  << " Frequency = " << frequency << "\n"
67  << " OutputFile = " << outputfile << "\n"
68  << " DoOutput = " << doOutput << "\n"
69  << " GetProv = " << getAllProvenances << "\n"
70  << " PrintProv = " << printProvenanceInfo << "\n"
71  << " Global Src = " << GlobalRecHitSrc_ << "\n"
72  << "===============================\n";
73 
74  }
75  //Put in analyzer stuff here....
76 
77  dbe = 0;
78 dbe = edm::Service<DQMStore>().operator->();
79 if (dbe) {
80  if (verbosity > 0 ) {
81  dbe->setVerbose(1);
82  } else {
83  dbe->setVerbose(0);
84  }
85 }
86 if (dbe) {
87  if (verbosity > 0 ) dbe->showDirStructure();
88  }
89 
90  Char_t hname[100];
91  Char_t htitle[100];
92 
93 //monitor elements
94 
95 //Si Strip
96  if(dbe)
97  {
98 string SiStripString[19] = {"TECW1", "TECW2", "TECW3", "TECW4", "TECW5", "TECW6", "TECW7", "TECW8", "TIBL1", "TIBL2", "TIBL3", "TIBL4", "TIDW1", "TIDW2", "TIDW3", "TOBL1", "TOBL2", "TOBL3", "TOBL4"};
99 for(int i = 0; i<19; ++i)
100 {
101  mehSiStripn[i]=0;
102  mehSiStripResX[i]=0;
103  mehSiStripResY[i]=0;
104 }
105  string hcharname, hchartitle;
106 dbe->setCurrentFolder("GlobalRecHitsV/SiStrips");
107 for(int amend = 0; amend < 19; ++amend)
108 {
109  hcharname = "hSiStripn_"+SiStripString[amend];
110  hchartitle= SiStripString[amend]+" rechits";
111  sprintf(hname, hcharname.c_str());
112  sprintf(htitle, hchartitle.c_str());
113  mehSiStripn[amend] = dbe->book1D(hname,htitle,20,0.,20.);
114  mehSiStripn[amend]->setAxisTitle("Number of hits in "+SiStripString[amend],1);
115  mehSiStripn[amend]->setAxisTitle("Count",2);
116  hcharname = "hSiStripResX_"+SiStripString[amend];
117  hchartitle= SiStripString[amend]+" rechit x resolution";
118  sprintf(hname, hcharname.c_str());
119  sprintf(htitle, hchartitle.c_str());
120  mehSiStripResX[amend] = dbe->book1D(hname,htitle,200,-0.02,.02);
121  mehSiStripResX[amend]->setAxisTitle("X-resolution in "+SiStripString[amend],1);
122  mehSiStripResX[amend]->setAxisTitle("Count",2);
123  hcharname = "hSiStripResY_"+SiStripString[amend];
124  hchartitle= SiStripString[amend]+" rechit y resolution";
125  sprintf(hname, hcharname.c_str());
126  sprintf(htitle, hchartitle.c_str());
127  mehSiStripResY[amend] = dbe->book1D(hname,htitle,200,-0.02,.02);
128  mehSiStripResY[amend]->setAxisTitle("Y-resolution in "+SiStripString[amend],1);
129  mehSiStripResY[amend]->setAxisTitle("Count",2);
130 }
131 
132 
133 //HCal
134 //string hcharname, hchartitle;
135 string HCalString[4]={"HB", "HE", "HF", "HO"};
136 float HCalnUpper[4]={3000.,3000.,3000.,2000.}; float HCalnLower[4]={2000.,2000.,2000.,1000.};
137 for(int j =0; j <4; ++j)
138 {
139  mehHcaln[j]=0;
140  mehHcalRes[j]=0;
141 }
142 
143 dbe->setCurrentFolder("GlobalRecHitsV/HCals");
144 for(int amend = 0; amend < 4; ++amend)
145 {
146  hcharname = "hHcaln_"+HCalString[amend];
147  hchartitle= HCalString[amend]+" rechits";
148  sprintf(hname, hcharname.c_str());
149  sprintf(htitle, hchartitle.c_str());
150  mehHcaln[amend] = dbe->book1D(hname,htitle, 500, HCalnLower[amend], HCalnUpper[amend]);
151  mehHcaln[amend]->setAxisTitle("Number of RecHits",1);
152  mehHcaln[amend]->setAxisTitle("Count",2);
153  hcharname = "hHcalRes_"+HCalString[amend];
154  hchartitle= HCalString[amend]+" rechit resolution";
155  sprintf(hname, hcharname.c_str());
156  sprintf(htitle, hchartitle.c_str());
157  mehHcalRes[amend] = dbe->book1D(hname,htitle, 25, -2., 2.);
158  mehHcalRes[amend]->setAxisTitle("RecHit E - SimHit E",1);
159  mehHcalRes[amend]->setAxisTitle("Count",2);
160 }
161 
162 
163 //Ecal
164 string ECalString[3] = {"EB","EE", "ES"};
165 int ECalnBins[3] = {700,100,50};
166 float ECalnUpper[3] = {20000., 62000., 300.};
167 float ECalnLower[3] = {6000., 60000., 100.};
168 int ECalResBins[3] = {200,200,200};
169 float ECalResUpper[3] = {1., 0.3, .0002};
170 float ECalResLower[3] = {-1., -0.3, -.0002};
171 for(int i =0; i<3; ++i)
172 {
173  mehEcaln[i]=0;
174  mehEcalRes[i]=0;
175 }
176 dbe->setCurrentFolder("GlobalRecHitsV/ECals");
177 
178 for(int amend = 0; amend < 3; ++amend)
179 {
180  hcharname = "hEcaln_"+ECalString[amend];
181  hchartitle= ECalString[amend]+" rechits";
182  sprintf(hname, hcharname.c_str());
183  sprintf(htitle, hchartitle.c_str());
184  mehEcaln[amend] = dbe->book1D(hname,htitle, ECalnBins[amend], ECalnLower[amend], ECalnUpper[amend]);
185  mehEcaln[amend]->setAxisTitle("Number of RecHits",1);
186  mehEcaln[amend]->setAxisTitle("Count",2);
187  hcharname = "hEcalRes_"+ECalString[amend];
188  hchartitle= ECalString[amend]+" rechit resolution";
189  sprintf(hname, hcharname.c_str());
190  sprintf(htitle, hchartitle.c_str());
191  mehEcalRes[amend] = dbe->book1D(hname,htitle,ECalResBins[amend], ECalResLower[amend], ECalResUpper[amend]);
192  mehEcalRes[amend]->setAxisTitle("RecHit E - SimHit E",1);
193  mehEcalRes[amend]->setAxisTitle("Count",2);
194 }
195 
196 
197 //Si Pixels
198 string SiPixelString[7] = {"BRL1", "BRL2", "BRL3", "FWD1n", "FWD1p", "FWD2n", "FWD2p"};
199 for(int j =0; j<7; ++j)
200 {
201  mehSiPixeln[j]=0;
202  mehSiPixelResX[j]=0;
203  mehSiPixelResY[j]=0;
204 }
205 
206 dbe->setCurrentFolder("GlobalRecHitsV/SiPixels");
207 for(int amend = 0; amend < 7; ++amend)
208 {
209  hcharname = "hSiPixeln_"+SiPixelString[amend];
210  hchartitle= SiPixelString[amend]+" rechits";
211  sprintf(hname, hcharname.c_str());
212  sprintf(htitle, hchartitle.c_str());
213  mehSiPixeln[amend] = dbe->book1D(hname,htitle,20,0.,20.);
214  mehSiPixeln[amend]->setAxisTitle("Number of hits in "+SiPixelString[amend],1);
215  mehSiPixeln[amend]->setAxisTitle("Count",2);
216  hcharname = "hSiPixelResX_"+SiPixelString[amend];
217  hchartitle= SiPixelString[amend]+" rechit x resolution";
218  sprintf(hname, hcharname.c_str());
219  sprintf(htitle, hchartitle.c_str());
220  mehSiPixelResX[amend] = dbe->book1D(hname,htitle,200,-0.02,.02);
221  mehSiPixelResX[amend]->setAxisTitle("X-resolution in "+SiPixelString[amend],1);
222  mehSiPixelResX[amend]->setAxisTitle("Count",2);
223  hcharname = "hSiPixelResY_"+SiPixelString[amend];
224  hchartitle= SiPixelString[amend]+" rechit y resolution";
225 
226  sprintf(hname, hcharname.c_str());
227  sprintf(htitle, hchartitle.c_str());
228  mehSiPixelResY[amend] = dbe->book1D(hname,htitle,200,-0.02,.02);
229  mehSiPixelResY[amend]->setAxisTitle("Y-resolution in "+SiPixelString[amend],1);
230  mehSiPixelResY[amend]->setAxisTitle("Count",2);
231 }
232 //Muons
233 dbe->setCurrentFolder("GlobalRecHitsV/Muons");
234 
235 mehDtMuonn = 0;
236 mehCSCn = 0;
237 mehRPCn = 0;
238 
239 //std::vector<MonitorElement *> me_List = {mehDtMuonn, mehCSCn, mehRPCn};
240 string n_List[3] = {"hDtMuonn", "hCSCn", "hRPCn"};
241 //float hist_prop[3] = [25., 0., 50.];
242 string hist_string[3] = {"Dt", "CSC", "RPC"};
243 
244 for(int amend=0; amend<3; ++amend)
245 {
246  hchartitle = hist_string[amend]+" rechits";
247  sprintf(hname, n_List[amend].c_str());
248  sprintf(htitle, hchartitle.c_str());
249  if(amend==0)
250  {
251  mehDtMuonn=dbe->book1D(hname,htitle,25, 0., 50.);
252  mehDtMuonn->setAxisTitle("Number of Rechits",1);
253  mehDtMuonn->setAxisTitle("Count",2);
254  }
255 if(amend==1)
256  {
257  mehCSCn=dbe->book1D(hname,htitle,25, 0., 50.);
258  mehCSCn->setAxisTitle("Number of Rechits",1);
259  mehCSCn->setAxisTitle("Count",2);
260  }
261 if(amend==2)
262  {
263  mehRPCn=dbe->book1D(hname,htitle,25, 0., 50.);
264  mehRPCn->setAxisTitle("Number of Rechits",1);
265  mehRPCn->setAxisTitle("Count",2);
266  }
267 }
268 
269 mehDtMuonRes=0;
270 mehCSCResRDPhi=0;
271 mehRPCResX=0;
272 
273 sprintf(hname, "hDtMuonRes");
274 sprintf(htitle, "DT wire distance resolution");
275 mehDtMuonRes = dbe->book1D(hname, htitle, 200, -0.2, 0.2);
276 sprintf(hname, "CSCResRDPhi");
277 sprintf(htitle, "CSC perp*dphi resolution");
278 mehCSCResRDPhi = dbe->book1D(hname, htitle, 200, -0.2, 0.2);
279 sprintf(hname,"hRPCResX");
280 sprintf(htitle, "RPC rechits x resolution");
281 mehRPCResX = dbe->book1D(hname, htitle, 50, -5., 5.);
282 }
283 
284 }
285 
286 
287 GlobalRecHitsHistogrammer::~GlobalRecHitsHistogrammer()
288 {
289  if (doOutput)
290  if (outputfile.size() != 0 && dbe) dbe->save(outputfile);
291 }
292 
293 void GlobalRecHitsHistogrammer::beginJob()
294 {
295  return;
296 }
297 
298 void GlobalRecHitsHistogrammer::endJob()
299 {
300  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_endJob";
301  if (verbosity >= 0)
302  edm::LogInfo(MsgLoggerCat)
303  << "Terminating having processed " << count << " events.";
304  return;
305 }
306 
307 void GlobalRecHitsHistogrammer::analyze(const edm::Event& iEvent,
308  const edm::EventSetup& iSetup)
309 {
310  std::string MsgLoggerCat = "GlobalRecHitsHistogrammer_analyze";
311 
312  // keep track of number of events processed
313  ++count;
314 
315  // get event id information
316  int nrun = iEvent.id().run();
317  int nevt = iEvent.id().event();
318 
319  if (verbosity > 0) {
320  edm::LogInfo(MsgLoggerCat)
321  << "Processing run " << nrun << ", event " << nevt
322  << " (" << count << " events total)";
323  } else if (verbosity == 0) {
324  if (nevt%frequency == 0 || nevt == 1) {
325  edm::LogInfo(MsgLoggerCat)
326  << "Processing run " << nrun << ", event " << nevt
327  << " (" << count << " events total)";
328  }
329  }
330 
331  // clear event holders
332  //clear(); Not in example I'm using, thus I comment it out.
333 
334  // look at information available in the event
335  if (getAllProvenances) {
336 
337  std::vector<const edm::Provenance*> AllProv;
338  iEvent.getAllProvenance(AllProv);
339 
340  if (verbosity >= 0)
341  edm::LogInfo(MsgLoggerCat)
342  << "Number of Provenances = " << AllProv.size();
343 
344  if (printProvenanceInfo && (verbosity >= 0)) {
345  TString eventout("\nProvenance info:\n");
346 
347  for (unsigned int i = 0; i < AllProv.size(); ++i) {
348  eventout += "\n ******************************";
349  eventout += "\n Module : ";
350  //eventout += (AllProv[i]->product).moduleLabel();
351  eventout += AllProv[i]->moduleLabel();
352  eventout += "\n ProductID : ";
353  //eventout += (AllProv[i]->product).productID_.id_;
354  eventout += AllProv[i]->productID().id();
355  eventout += "\n ClassName : ";
356  //eventout += (AllProv[i]->product).fullClassName_;
357  eventout += AllProv[i]->className();
358  eventout += "\n InstanceName : ";
359  //eventout += (AllProv[i]->product).productInstanceName_;
360  eventout += AllProv[i]->productInstanceName();
361  eventout += "\n BranchName : ";
362  //eventout += (AllProv[i]->product).branchName_;
363  eventout += AllProv[i]->branchName();
364  }
365  eventout += "\n ******************************\n";
366  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
367  printProvenanceInfo = false;
368  }
369  getAllProvenances = false;
370  }
371 
372 edm::Handle<PGlobalRecHit> srcGlobalRecHits;
373  iEvent.getByLabel(GlobalRecHitSrc_,srcGlobalRecHits);
374  if (!srcGlobalRecHits.isValid()) {
375  edm::LogWarning(MsgLoggerCat)
376  << "Unable to find PGlobalRecHit in event!";
377  return;
378  }
379 
380  int nEBCalRecHits = srcGlobalRecHits->getnEBCalRecHits();
381  int nEECalRecHits = srcGlobalRecHits->getnEECalRecHits();
382  int nESCalRecHits = srcGlobalRecHits->getnESCalRecHits();
383 
384  int nHBCalRecHits = srcGlobalRecHits->getnHBCalRecHits();
385  int nHECalRecHits = srcGlobalRecHits->getnHECalRecHits();
386  int nHOCalRecHits = srcGlobalRecHits->getnHOCalRecHits();
387  int nHFCalRecHits = srcGlobalRecHits->getnHFCalRecHits();
388 
389  int nTIBL1RecHits = srcGlobalRecHits->getnTIBL1RecHits();
390  int nTIBL2RecHits = srcGlobalRecHits->getnTIBL2RecHits();
391  int nTIBL3RecHits = srcGlobalRecHits->getnTIBL3RecHits();
392  int nTIBL4RecHits = srcGlobalRecHits->getnTIBL4RecHits();
393  int nTOBL1RecHits = srcGlobalRecHits->getnTOBL1RecHits();
394  int nTOBL2RecHits = srcGlobalRecHits->getnTOBL2RecHits();
395  int nTOBL3RecHits = srcGlobalRecHits->getnTOBL3RecHits();
396  int nTOBL4RecHits = srcGlobalRecHits->getnTOBL4RecHits();
397  int nTIDW1RecHits = srcGlobalRecHits->getnTIDW1RecHits();
398  int nTIDW2RecHits = srcGlobalRecHits->getnTIDW2RecHits();
399  int nTIDW3RecHits = srcGlobalRecHits->getnTIDW3RecHits();
400  int nTECW1RecHits = srcGlobalRecHits->getnTECW1RecHits();
401  int nTECW2RecHits = srcGlobalRecHits->getnTECW2RecHits();
402  int nTECW3RecHits = srcGlobalRecHits->getnTECW3RecHits();
403  int nTECW4RecHits = srcGlobalRecHits->getnTECW4RecHits();
404  int nTECW5RecHits = srcGlobalRecHits->getnTECW5RecHits();
405  int nTECW6RecHits = srcGlobalRecHits->getnTECW6RecHits();
406  int nTECW7RecHits = srcGlobalRecHits->getnTECW7RecHits();
407  int nTECW8RecHits = srcGlobalRecHits->getnTECW8RecHits();
408 
409  int nBRL1RecHits = srcGlobalRecHits->getnBRL1RecHits();
410  int nBRL2RecHits = srcGlobalRecHits->getnBRL2RecHits();
411  int nBRL3RecHits = srcGlobalRecHits->getnBRL3RecHits();
412  int nFWD1nRecHits = srcGlobalRecHits->getnFWD1nRecHits();
413  int nFWD1pRecHits = srcGlobalRecHits->getnFWD1pRecHits();
414  int nFWD2nRecHits = srcGlobalRecHits->getnFWD2nRecHits();
415  int nFWD2pRecHits = srcGlobalRecHits->getnFWD2pRecHits();
416 
417  int nDTRecHits = srcGlobalRecHits->getnDTRecHits();
418 
419  int nCSCRecHits = srcGlobalRecHits->getnCSCRecHits();
420 
421  int nRPCRecHits = srcGlobalRecHits->getnRPCRecHits();
422 
423  // get Ecal info
424  std::vector<PGlobalRecHit::ECalRecHit> EECalRecHits =
425  srcGlobalRecHits->getEECalRecHits();
426  mehEcaln[0]->Fill((float)nEECalRecHits);
427  for (unsigned int i = 0; i < EECalRecHits.size(); ++i) {
428  mehEcalRes[0]->Fill(EECalRecHits[i].RE - EECalRecHits[i].SHE);
429  }
430 
431  std::vector<PGlobalRecHit::ECalRecHit> EBCalRecHits =
432  srcGlobalRecHits->getEBCalRecHits();
433  mehEcaln[1]->Fill((float)nEBCalRecHits);
434  for (unsigned int i = 0; i < EBCalRecHits.size(); ++i) {
435  mehEcalRes[1]->Fill(EBCalRecHits[i].RE - EBCalRecHits[i].SHE);
436  }
437 
438  std::vector<PGlobalRecHit::ECalRecHit> ESCalRecHits =
439  srcGlobalRecHits->getESCalRecHits();
440  mehEcaln[2]->Fill((float)nESCalRecHits);
441  for (unsigned int i = 0; i < ESCalRecHits.size(); ++i) {
442  mehEcalRes[2]->Fill(ESCalRecHits[i].RE - ESCalRecHits[i].SHE);
443  }
444 
445  // Get HCal info
446  std::vector<PGlobalRecHit::HCalRecHit> HBCalRecHits =
447  srcGlobalRecHits->getHBCalRecHits();
448  mehHcaln[0]->Fill((float)nHBCalRecHits);
449  for (unsigned int i = 0; i < HBCalRecHits.size(); ++i) {
450  mehHcalRes[0]->Fill(HBCalRecHits[i].REC - HBCalRecHits[i].SHE);
451  }
452 
453  std::vector<PGlobalRecHit::HCalRecHit> HECalRecHits =
454  srcGlobalRecHits->getHECalRecHits();
455  mehHcaln[1]->Fill((float)nHECalRecHits);
456  for (unsigned int i = 0; i < HECalRecHits.size(); ++i) {
457  mehHcalRes[1]->Fill(HECalRecHits[i].REC - HECalRecHits[i].SHE);
458  }
459 
460  std::vector<PGlobalRecHit::HCalRecHit> HOCalRecHits =
461  srcGlobalRecHits->getHOCalRecHits();
462  mehHcaln[2]->Fill((float)nHOCalRecHits);
463  for (unsigned int i = 0; i < HOCalRecHits.size(); ++i) {
464  mehHcalRes[2]->Fill(HOCalRecHits[i].REC - HOCalRecHits[i].SHE);
465  }
466 
467  std::vector<PGlobalRecHit::HCalRecHit> HFCalRecHits =
468  srcGlobalRecHits->getHFCalRecHits();
469  mehHcaln[3]->Fill((float)nHFCalRecHits);
470  for (unsigned int i = 0; i < HFCalRecHits.size(); ++i) {
471  mehHcalRes[3]->Fill(HFCalRecHits[i].REC - HFCalRecHits[i].SHE);
472  }
473 
474  // get SiStrip info
475  std::vector<PGlobalRecHit::SiStripRecHit> TIBL1RecHits =
476  srcGlobalRecHits->getTIBL1RecHits();
477  mehSiStripn[0]->Fill((float)nTIBL1RecHits);
478  for (unsigned int i = 0; i < TIBL1RecHits.size(); ++i) {
479  mehSiStripResX[0]->Fill(TIBL1RecHits[i].RX - TIBL1RecHits[i].SX);
480  mehSiStripResY[0]->Fill(TIBL1RecHits[i].RY - TIBL1RecHits[i].SY);
481  }
482 
483  std::vector<PGlobalRecHit::SiStripRecHit> TIBL2RecHits =
484  srcGlobalRecHits->getTIBL2RecHits();
485  mehSiStripn[1]->Fill((float)nTIBL2RecHits);
486  for (unsigned int i = 0; i < TIBL2RecHits.size(); ++i) {
487  mehSiStripResX[1]->Fill(TIBL2RecHits[i].RX - TIBL2RecHits[i].SX);
488  mehSiStripResY[1]->Fill(TIBL2RecHits[i].RY - TIBL2RecHits[i].SY);
489  }
490 
491  std::vector<PGlobalRecHit::SiStripRecHit> TIBL3RecHits =
492  srcGlobalRecHits->getTIBL3RecHits();
493  mehSiStripn[2]->Fill((float)nTIBL3RecHits);
494  for (unsigned int i = 0; i < TIBL3RecHits.size(); ++i) {
495  mehSiStripResX[2]->Fill(TIBL3RecHits[i].RX - TIBL3RecHits[i].SX);
496  mehSiStripResY[2]->Fill(TIBL3RecHits[i].RY - TIBL3RecHits[i].SY);
497  }
498 
499  std::vector<PGlobalRecHit::SiStripRecHit> TIBL4RecHits =
500  srcGlobalRecHits->getTIBL4RecHits();
501  mehSiStripn[3]->Fill((float)nTIBL4RecHits);
502  for (unsigned int i = 0; i < TIBL4RecHits.size(); ++i) {
503  mehSiStripResX[3]->Fill(TIBL4RecHits[i].RX - TIBL4RecHits[i].SX);
504  mehSiStripResY[3]->Fill(TIBL4RecHits[i].RY - TIBL4RecHits[i].SY);
505  }
506 
507  std::vector<PGlobalRecHit::SiStripRecHit> TOBL1RecHits =
508  srcGlobalRecHits->getTOBL1RecHits();
509  mehSiStripn[4]->Fill((float)nTOBL1RecHits);
510  for (unsigned int i = 0; i < TOBL1RecHits.size(); ++i) {
511  mehSiStripResX[4]->Fill(TOBL1RecHits[i].RX - TOBL1RecHits[i].SX);
512  mehSiStripResY[4]->Fill(TOBL1RecHits[i].RY - TOBL1RecHits[i].SY);
513  }
514 
515  std::vector<PGlobalRecHit::SiStripRecHit> TOBL2RecHits =
516  srcGlobalRecHits->getTOBL2RecHits();
517  mehSiStripn[5]->Fill((float)nTOBL2RecHits);
518  for (unsigned int i = 0; i < TOBL2RecHits.size(); ++i) {
519  mehSiStripResX[5]->Fill(TOBL2RecHits[i].RX - TOBL2RecHits[i].SX);
520  mehSiStripResY[5]->Fill(TOBL2RecHits[i].RY - TOBL2RecHits[i].SY);
521  }
522 
523  std::vector<PGlobalRecHit::SiStripRecHit> TOBL3RecHits =
524  srcGlobalRecHits->getTOBL3RecHits();
525  mehSiStripn[6]->Fill((float)nTOBL3RecHits);
526  for (unsigned int i = 0; i < TOBL3RecHits.size(); ++i) {
527  mehSiStripResX[6]->Fill(TOBL3RecHits[i].RX - TOBL3RecHits[i].SX);
528  mehSiStripResY[6]->Fill(TOBL3RecHits[i].RY - TOBL3RecHits[i].SY);
529  }
530 
531  std::vector<PGlobalRecHit::SiStripRecHit> TOBL4RecHits =
532  srcGlobalRecHits->getTOBL4RecHits();
533  mehSiStripn[7]->Fill((float)nTOBL4RecHits);
534  for (unsigned int i = 0; i < TOBL4RecHits.size(); ++i) {
535  mehSiStripResX[7]->Fill(TOBL4RecHits[i].RX - TOBL4RecHits[i].SX);
536  mehSiStripResY[7]->Fill(TOBL4RecHits[i].RY - TOBL4RecHits[i].SY);
537  }
538 
539  std::vector<PGlobalRecHit::SiStripRecHit> TIDW1RecHits =
540  srcGlobalRecHits->getTIDW1RecHits();
541  mehSiStripn[8]->Fill((float)nTIDW1RecHits);
542  for (unsigned int i = 0; i < TIDW1RecHits.size(); ++i) {
543  mehSiStripResX[8]->Fill(TIDW1RecHits[i].RX - TIDW1RecHits[i].SX);
544  mehSiStripResY[8]->Fill(TIDW1RecHits[i].RY - TIDW1RecHits[i].SY);
545  }
546 
547  std::vector<PGlobalRecHit::SiStripRecHit> TIDW2RecHits =
548  srcGlobalRecHits->getTIDW2RecHits();
549  mehSiStripn[9]->Fill((float)nTIDW2RecHits);
550  for (unsigned int i = 0; i < TIDW2RecHits.size(); ++i) {
551  mehSiStripResX[9]->Fill(TIDW2RecHits[i].RX - TIDW2RecHits[i].SX);
552  mehSiStripResY[9]->Fill(TIDW2RecHits[i].RY - TIDW2RecHits[i].SY);
553  }
554 
555  std::vector<PGlobalRecHit::SiStripRecHit> TIDW3RecHits =
556  srcGlobalRecHits->getTIDW3RecHits();
557  mehSiStripn[10]->Fill((float)nTIDW3RecHits);
558  for (unsigned int i = 0; i < TIDW3RecHits.size(); ++i) {
559  mehSiStripResX[10]->Fill(TIDW3RecHits[i].RX - TIDW3RecHits[i].SX);
560  mehSiStripResY[10]->Fill(TIDW3RecHits[i].RY - TIDW3RecHits[i].SY);
561  }
562 
563  std::vector<PGlobalRecHit::SiStripRecHit> TECW1RecHits =
564  srcGlobalRecHits->getTECW1RecHits();
565  mehSiStripn[11]->Fill((float)nTECW1RecHits);
566  for (unsigned int i = 0; i < TECW1RecHits.size(); ++i) {
567  mehSiStripResX[11]->Fill(TECW1RecHits[i].RX - TECW1RecHits[i].SX);
568  mehSiStripResY[11]->Fill(TECW1RecHits[i].RY - TECW1RecHits[i].SY);
569  }
570 
571  std::vector<PGlobalRecHit::SiStripRecHit> TECW2RecHits =
572  srcGlobalRecHits->getTECW2RecHits();
573  mehSiStripn[12]->Fill((float)nTECW2RecHits);
574  for (unsigned int i = 0; i < TECW2RecHits.size(); ++i) {
575  mehSiStripResX[12]->Fill(TECW2RecHits[i].RX - TECW2RecHits[i].SX);
576  mehSiStripResY[12]->Fill(TECW2RecHits[i].RY - TECW2RecHits[i].SY);
577  }
578 
579  std::vector<PGlobalRecHit::SiStripRecHit> TECW3RecHits =
580  srcGlobalRecHits->getTECW3RecHits();
581  mehSiStripn[13]->Fill((float)nTECW3RecHits);
582  for (unsigned int i = 0; i < TECW3RecHits.size(); ++i) {
583  mehSiStripResX[13]->Fill(TECW3RecHits[i].RX - TECW3RecHits[i].SX);
584  mehSiStripResY[13]->Fill(TECW3RecHits[i].RY - TECW3RecHits[i].SY);
585  }
586 
587  std::vector<PGlobalRecHit::SiStripRecHit> TECW4RecHits =
588  srcGlobalRecHits->getTECW4RecHits();
589  mehSiStripn[14]->Fill((float)nTECW4RecHits);
590  for (unsigned int i = 0; i < TECW4RecHits.size(); ++i) {
591  mehSiStripResX[14]->Fill(TECW4RecHits[i].RX - TECW4RecHits[i].SX);
592  mehSiStripResY[14]->Fill(TECW4RecHits[i].RY - TECW4RecHits[i].SY);
593  }
594 
595  std::vector<PGlobalRecHit::SiStripRecHit> TECW5RecHits =
596  srcGlobalRecHits->getTECW5RecHits();
597  mehSiStripn[15]->Fill((float)nTECW5RecHits);
598  for (unsigned int i = 0; i < TECW5RecHits.size(); ++i) {
599  mehSiStripResX[15]->Fill(TECW5RecHits[i].RX - TECW5RecHits[i].SX);
600  mehSiStripResY[15]->Fill(TECW5RecHits[i].RY - TECW5RecHits[i].SY);
601  }
602 
603  std::vector<PGlobalRecHit::SiStripRecHit> TECW6RecHits =
604  srcGlobalRecHits->getTECW6RecHits();
605  mehSiStripn[16]->Fill((float)nTECW6RecHits);
606  for (unsigned int i = 0; i < TECW6RecHits.size(); ++i) {
607  mehSiStripResX[16]->Fill(TECW6RecHits[i].RX - TECW6RecHits[i].SX);
608  mehSiStripResY[16]->Fill(TECW6RecHits[i].RY - TECW6RecHits[i].SY);
609  }
610 
611  std::vector<PGlobalRecHit::SiStripRecHit> TECW7RecHits =
612  srcGlobalRecHits->getTECW7RecHits();
613  mehSiStripn[17]->Fill((float)nTECW7RecHits);
614  for (unsigned int i = 0; i < TECW7RecHits.size(); ++i) {
615  mehSiStripResX[17]->Fill(TECW7RecHits[i].RX - TECW7RecHits[i].SX);
616  mehSiStripResY[17]->Fill(TECW7RecHits[i].RY - TECW7RecHits[i].SY);
617  }
618 
619  std::vector<PGlobalRecHit::SiStripRecHit> TECW8RecHits =
620  srcGlobalRecHits->getTECW8RecHits();
621  mehSiStripn[18]->Fill((float)nTECW8RecHits);
622  for (unsigned int i = 0; i < TECW8RecHits.size(); ++i) {
623  mehSiStripResX[18]->Fill(TECW8RecHits[i].RX - TECW8RecHits[i].SX);
624  mehSiStripResY[18]->Fill(TECW8RecHits[i].RY - TECW8RecHits[i].SY);
625  }
626 
627  // get SiPixel info
628  std::vector<PGlobalRecHit::SiPixelRecHit> BRL1RecHits =
629  srcGlobalRecHits->getBRL1RecHits();
630  mehSiPixeln[0]->Fill((float)nBRL1RecHits);
631  for (unsigned int i = 0; i < BRL1RecHits.size(); ++i) {
632  mehSiPixelResX[0]->Fill(BRL1RecHits[i].RX - BRL1RecHits[i].SX);
633  mehSiPixelResY[0]->Fill(BRL1RecHits[i].RY - BRL1RecHits[i].SY);
634  }
635 
636  std::vector<PGlobalRecHit::SiPixelRecHit> BRL2RecHits =
637  srcGlobalRecHits->getBRL2RecHits();
638  mehSiPixeln[1]->Fill((float)nBRL2RecHits);
639  for (unsigned int i = 0; i < BRL2RecHits.size(); ++i) {
640  mehSiPixelResX[1]->Fill(BRL2RecHits[i].RX - BRL2RecHits[i].SX);
641  mehSiPixelResY[1]->Fill(BRL2RecHits[i].RY - BRL2RecHits[i].SY);
642  }
643 
644  std::vector<PGlobalRecHit::SiPixelRecHit> BRL3RecHits =
645  srcGlobalRecHits->getBRL3RecHits();
646  mehSiPixeln[2]->Fill((float)nBRL3RecHits);
647  for (unsigned int i = 0; i < BRL3RecHits.size(); ++i) {
648  mehSiPixelResX[2]->Fill(BRL3RecHits[i].RX - BRL3RecHits[i].SX);
649  mehSiPixelResY[2]->Fill(BRL3RecHits[i].RY - BRL3RecHits[i].SY);
650  }
651 
652  std::vector<PGlobalRecHit::SiPixelRecHit> FWD1pRecHits =
653  srcGlobalRecHits->getFWD1pRecHits();
654  mehSiPixeln[3]->Fill((float)nFWD1pRecHits);
655  for (unsigned int i = 0; i < FWD1pRecHits.size(); ++i) {
656  mehSiPixelResX[3]->Fill(FWD1pRecHits[i].RX - FWD1pRecHits[i].SX);
657  mehSiPixelResY[3]->Fill(FWD1pRecHits[i].RY - FWD1pRecHits[i].SY);
658  }
659 
660  std::vector<PGlobalRecHit::SiPixelRecHit> FWD1nRecHits =
661  srcGlobalRecHits->getFWD1nRecHits();
662  mehSiPixeln[4]->Fill((float)nFWD1nRecHits);
663  for (unsigned int i = 0; i < FWD1nRecHits.size(); ++i) {
664  mehSiPixelResX[4]->Fill(FWD1nRecHits[i].RX - FWD1nRecHits[i].SX);
665  mehSiPixelResY[4]->Fill(FWD1nRecHits[i].RY - FWD1nRecHits[i].SY);
666  }
667 
668  std::vector<PGlobalRecHit::SiPixelRecHit> FWD2pRecHits =
669  srcGlobalRecHits->getFWD2pRecHits();
670  mehSiPixeln[5]->Fill((float)nFWD2pRecHits);
671  for (unsigned int i = 0; i < FWD2pRecHits.size(); ++i) {
672  mehSiPixelResX[5]->Fill(FWD2pRecHits[i].RX - FWD2pRecHits[i].SX);
673  mehSiPixelResY[5]->Fill(FWD2pRecHits[i].RY - FWD2pRecHits[i].SY);
674  }
675 
676  std::vector<PGlobalRecHit::SiPixelRecHit> FWD2nRecHits =
677  srcGlobalRecHits->getFWD2nRecHits();
678  mehSiPixeln[6]->Fill((float)nFWD2nRecHits);
679  for (unsigned int i = 0; i < FWD2nRecHits.size(); ++i) {
680  mehSiPixelResX[6]->Fill(FWD2nRecHits[i].RX - FWD2nRecHits[i].SX);
681  mehSiPixelResY[6]->Fill(FWD2nRecHits[i].RY - FWD2nRecHits[i].SY);
682  }
683 
684  // get DtMuon info
685  std::vector<PGlobalRecHit::DTRecHit> DTRecHits =
686  srcGlobalRecHits->getDTRecHits();
687  mehDtMuonn->Fill((float)nDTRecHits);
688  for (unsigned int i = 0; i < DTRecHits.size(); ++i) {
689  mehDtMuonRes->Fill(DTRecHits[i].RHD - DTRecHits[i].SHD);
690  }
691 
692  // get CSC info
693  std::vector<PGlobalRecHit::CSCRecHit> CSCRecHits =
694  srcGlobalRecHits->getCSCRecHits();
695  mehCSCn->Fill((float)nCSCRecHits);
696  for (unsigned int i = 0; i < CSCRecHits.size(); ++i) {
697  mehCSCResRDPhi->Fill(CSCRecHits[i].RHPERP *
698  (CSCRecHits[i].RHPHI - CSCRecHits[i].SHPHI));
699  }
700 
701  // get RPC info
702  std::vector<PGlobalRecHit::RPCRecHit> RPCRecHits =
703  srcGlobalRecHits->getRPCRecHits();
704  mehRPCn->Fill((float)nRPCRecHits);
705  for (unsigned int i = 0; i < RPCRecHits.size(); ++i) {
706  mehRPCResX->Fill(RPCRecHits[i].RHX - RPCRecHits[i].SHX);
707  }
708 
709  if (verbosity > 0)
710  edm::LogInfo (MsgLoggerCat)
711  << "Done gathering data from event.";
712 
713 }
714 
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:44
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
edm::Event::getAllProvenance
void getAllProvenance(std::vector< Provenance const * > &provenances) const
Definition: Event.cc:71
i
int i
Definition: DBlmapReader.cc:9
label
const std::string & label
Definition: MVAComputer.cc:186
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:519
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:1883
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
iEvent
int iEvent
Definition: GenABIO.cc:243
j
int j
Definition: DBlmapReader.cc:9
GlobalRecHitsHistogrammer::GlobalRecHitsHistogrammer
GlobalRecHitsHistogrammer(const edm::ParameterSet &)
Definition: GlobalRecHitsHistogrammer.cc:15
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:201
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
verbosity
const int verbosity
Definition: DQMFileLoader.cc:25
GlobalRecHitsHistogrammer::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: GlobalRecHitsHistogrammer.cc:307
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
DQMStore::showDirStructure
void showDirStructure(void) const
Definition: DQMStore.cc:2539
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:903
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:237
Generated for CMSSW Reference Manual by doxygen 1.8.5