CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PhotonOfflineClient.cc
Go to the documentation of this file.
1 #include <iostream>
2 //
3 
4 #include "DQMOffline/EGamma/plugins/PhotonOfflineClient.h"
5 
6 
7 //#define TWOPI 6.283185308
8 //
9 
23 using namespace std;
24 using std::cout;
25 
26 PhotonOfflineClient::PhotonOfflineClient(const edm::ParameterSet& pset)
27 {
28 
29  dbe_ = 0;
30  dbe_ = edm::Service<DQMStore>().operator->();
31  dbe_->setVerbose(0);
32  parameters_ = pset;
33 
34  cutStep_ = pset.getParameter<double>("cutStep");
35  numberOfSteps_ = pset.getParameter<int>("numberOfSteps");
36 
37  etMin = pset.getParameter<double>("etMin");
38  etMax = pset.getParameter<double>("etMax");
39  etBin = pset.getParameter<int>("etBin");
40  etaMin = pset.getParameter<double>("etaMin");
41  etaMax = pset.getParameter<double>("etaMax");
42  etaBin = pset.getParameter<int>("etaBin");
43  phiMin = pset.getParameter<double>("phiMin");
44  phiMax = pset.getParameter<double>("phiMax");
45  phiBin = pset.getParameter<int>("phiBin");
46 
47  standAlone_ = pset.getParameter<bool>("standAlone");
48  batch_ = pset.getParameter<bool>("batch");
49  excludeBkgHistos_ = pset.getParameter<bool>("excludeBkgHistos");
50 
51  outputFileName_ = pset.getParameter<string>("OutputFileName");
52  inputFileName_ = pset.getUntrackedParameter<string>("InputFileName");
53 
54  histo_index_photons_ = 0;
55  histo_index_conversions_ = 0;
56  histo_index_efficiency_ = 0;
57  histo_index_invMass_ = 0;
58 
59  types_.push_back("All");
60  types_.push_back("GoodCandidate");
61  if (!excludeBkgHistos_) types_.push_back("Background");
62 
63  parts_.push_back("AllEcal");
64  parts_.push_back("Barrel");
65  parts_.push_back("Endcaps");
66 
67 }
68 
69 PhotonOfflineClient::~PhotonOfflineClient()
70 {}
71 void PhotonOfflineClient::beginJob()
72 {}
73 void PhotonOfflineClient::analyze(const edm::Event& e, const edm::EventSetup& esup)
74 {}
75 void PhotonOfflineClient::endLuminosityBlock( const edm::LuminosityBlock& , const edm::EventSetup& setup)
76 {}
77 
78 
79 
80 void PhotonOfflineClient::endJob()
81 {
82  if(standAlone_) runClient();
83 }
84 void PhotonOfflineClient::endRun(const edm::Run& run, const edm::EventSetup& setup)
85 {
86  if(!standAlone_) runClient();
87 }
88 
89 
90 
91 
92 void PhotonOfflineClient::runClient()
93 {
94  if(!dbe_) return;
95 
96  if(batch_) dbe_->open(inputFileName_);
97 
98  if(!dbe_->dirExists("Egamma/PhotonAnalyzer")){
99  //cout << "egamma directory doesn't exist..." << std::endl;
100  return;
101  }
102 
103  //find out how many histograms are in the various folders
104  histo_index_photons_ = dbe_->get("Egamma/PhotonAnalyzer/numberOfHistogramsInPhotonsFolder")->getIntValue();
105  histo_index_conversions_ = dbe_->get("Egamma/PhotonAnalyzer/numberOfHistogramsInConversionsFolder")->getIntValue();
106  histo_index_efficiency_ = dbe_->get("Egamma/PhotonAnalyzer/numberOfHistogramsInEfficiencyFolder")->getIntValue();
107  histo_index_invMass_ = dbe_->get("Egamma/PhotonAnalyzer/numberOfHistogramsInInvMassFolder")->getIntValue();
108 
109  dbe_->setCurrentFolder("Egamma/PhotonAnalyzer/");
110  dbe_->removeElement("numberOfHistogramsInPhotonsFolder");
111  dbe_->removeElement("numberOfHistogramsInConversionsFolder");
112  dbe_->removeElement("numberOfHistogramsInEfficiencyFolder");
113  dbe_->removeElement("numberOfHistogramsInInvMassFolder");
114 
115 
116  string AllPath = "Egamma/PhotonAnalyzer/AllPhotons/";
117  string IsoPath = "Egamma/PhotonAnalyzer/GoodCandidatePhotons/";
118  string NonisoPath = "Egamma/PhotonAnalyzer/BackgroundPhotons/";
119  string EffPath = "Egamma/PhotonAnalyzer/Efficiencies/";
120 
121 
122  //booking efficiency histograms
123  dbe_->setCurrentFolder(EffPath);
124 
125  p_efficiencyVsEtaLoose_ = bookHisto("EfficiencyVsEtaLoose","Fraction of Photons passing Loose Isolation vs #eta;#eta",etaBin,etaMin, etaMax);
126  p_efficiencyVsEtLoose_ = bookHisto("EfficiencyVsEtLoose", "Fraction of Photons passing Loose Isolation vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
127  p_efficiencyVsEtaTight_ = bookHisto("EfficiencyVsEtaTight","Fraction of Photons passing Tight Isolation vs #eta;#eta",etaBin,etaMin, etaMax);
128  p_efficiencyVsEtTight_ = bookHisto("EfficiencyVsEtTight", "Fraction of Photons passing Tight Isolation vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
129 
130 
131  p_efficiencyVsEtaHLT_ = bookHisto("EfficiencyVsEtaHLT","Fraction of Photons firing HLT vs #eta;#eta",etaBin,etaMin, etaMax);
132  p_efficiencyVsEtHLT_ = bookHisto("EfficiencyVsEtHLT", "Fraction of Photons firing HLT vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
133 
134 
135  p_convFractionVsEtaLoose_ = bookHisto("ConvFractionVsEtaLoose","Fraction of Loosely Isolated Photons which are matched to two tracks vs #eta;#eta",etaBin,etaMin, etaMax);
136  p_convFractionVsEtLoose_ = bookHisto("ConvFractionVsEtLoose", "Fraction of Loosely Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
137  p_convFractionVsEtaTight_ = bookHisto("ConvFractionVsEtaTight","Fraction of Tightly Isolated Photons which are matched to two tracks vs #eta;#eta",etaBin,etaMin, etaMax);
138  p_convFractionVsEtTight_ = bookHisto("ConvFractionVsEtTight", "Fraction of Tightly Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
139 
140  p_vertexReconstructionEfficiencyVsEta_ = bookHisto("VertexReconstructionEfficiencyVsEta","Fraction of Converted Photons which have a valid vertex vs #eta;#eta",etaBin,etaMin, etaMax);
141 
142 
143  //booking conversion fraction histograms
144  dbe_->setCurrentFolder(AllPath+"Et above 20 GeV/Conversions");
145  p_convFractionVsEt_ = book2DHistoVector("1D","convFractionVsEt", "Fraction of Converted Photons vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
146  p_convFractionVsPhi_ = book3DHistoVector("1D","convFractionVsPhi","Fraction of Converted Photons vs #phi;#phi",phiBin,phiMin,phiMax);
147  p_convFractionVsEta_ = book2DHistoVector("1D","convFractionVsEta","Fraction of Converted Photons vs #eta;#eta",etaBin,etaMin,etaMax);
148 
149 
150  //booking bad channel fraction histograms
151  dbe_->setCurrentFolder(AllPath+"Et above 20 GeV/");
152  p_badChannelsFractionVsPhi_ = book2DHistoVector("1D","badChannelsFractionVsPhi","Fraction of Photons which have at least one bad channel vs #phi;#phi",phiBin,phiMin,phiMax);
153  p_badChannelsFractionVsEta_ = book2DHistoVector("1D","badChannelsFractionVsEta","Fraction of Photons which have at least one bad channel vs #eta;#eta",etaBin,etaMin, etaMax);
154  p_badChannelsFractionVsEt_ = book2DHistoVector("1D","badChannelsFractionVsEt", "Fraction of Photons which have at least one bad channel vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
155 
156 
157 
158  //making efficiency plots
159 
160  MonitorElement * dividend;
161  MonitorElement * numerator;
162  MonitorElement * denominator;
163 
164  currentFolder_.str("");
165  currentFolder_ << AllPath << "Et above 20 GeV/";
166 
167  //HLT efficiency plots
168  dividend = retrieveHisto(EffPath,"EfficiencyVsEtaHLT");
169  numerator = retrieveHisto(EffPath,"phoEtaPostHLT");
170  denominator = retrieveHisto(EffPath,"phoEtaPreHLT");
171  dividePlots(dividend,numerator,denominator);
172 
173  dividend = retrieveHisto(EffPath,"EfficiencyVsEtHLT");
174  numerator = retrieveHisto(EffPath,"phoEtPostHLT");
175  denominator = retrieveHisto(EffPath,"phoEtPreHLT");
176  dividePlots(dividend,numerator,denominator);
177 
178 
179  //efficiencies vs Eta
180  denominator = retrieveHisto(currentFolder_.str(),"phoEta");
181 
182  dividend = retrieveHisto(EffPath,"EfficiencyVsEtaLoose");
183  numerator = retrieveHisto(EffPath,"phoEtaLoose");
184  dividePlots(dividend,numerator,denominator);
185 
186  dividend = retrieveHisto(EffPath,"EfficiencyVsEtaTight");
187  numerator = retrieveHisto(EffPath,"phoEtaTight");
188  dividePlots(dividend,numerator,denominator);
189 
190 
191  //efficiencies vs Et
192  denominator = retrieveHisto(currentFolder_.str(),"phoEtAllEcal");
193 
194  dividend = retrieveHisto(EffPath,"EfficiencyVsEtLoose");
195  numerator = retrieveHisto(EffPath,"phoEtLoose");
196  dividePlots(dividend,numerator,denominator);
197 
198  dividend = retrieveHisto(EffPath,"EfficiencyVsEtTight");
199  numerator = retrieveHisto(EffPath,"phoEtTight");
200  dividePlots(dividend,numerator,denominator);
201 
202 
203  //conversion fractions vs Eta
204  dividend = retrieveHisto(EffPath,"ConvFractionVsEtaLoose");
205  numerator = retrieveHisto(EffPath,"convEtaLoose");
206  denominator = retrieveHisto(EffPath,"phoEtaLoose");
207  dividePlots(dividend,numerator,denominator);
208 
209  dividend = retrieveHisto(EffPath,"ConvFractionVsEtaTight");
210  numerator = retrieveHisto(EffPath,"convEtaTight");
211  denominator = retrieveHisto(EffPath,"phoEtaTight");
212  dividePlots(dividend,numerator,denominator);
213 
214 
215  //conversion fractions vs Et
216  dividend = retrieveHisto(EffPath,"ConvFractionVsEtLoose");
217  numerator = retrieveHisto(EffPath,"convEtLoose");
218  denominator = retrieveHisto(EffPath,"phoEtLoose");
219  dividePlots(dividend,numerator,denominator);
220 
221  dividend = retrieveHisto(EffPath,"ConvFractionVsEtTight");
222  numerator = retrieveHisto(EffPath,"convEtTight");
223  denominator = retrieveHisto(EffPath,"phoEtTight");
224  dividePlots(dividend,numerator,denominator);
225 
226 
227  //conversion vertex recontruction efficiency
228  dividend = retrieveHisto(EffPath,"VertexReconstructionEfficiencyVsEta");
229  numerator = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvEta");
230  denominator = retrieveHisto(EffPath,"phoEtaVertex");
231  dividePlots(dividend,numerator,denominator);
232 
233 
234  dbe_->setCurrentFolder(EffPath);
235 
236 
237  dbe_->removeElement("phoEtaPreHLT");
238  dbe_->removeElement("phoEtPreHLT");
239  dbe_->removeElement("phoEtaPostHLT");
240  dbe_->removeElement("phoEtPostHLT");
241  dbe_->removeElement("phoEtaLoose");
242  dbe_->removeElement("phoEtaTight");
243  dbe_->removeElement("phoEtLoose");
244  dbe_->removeElement("phoEtTight");
245  dbe_->removeElement("phoEtaVertex");
246  dbe_->removeElement("convEtaLoose");
247  dbe_->removeElement("convEtaTight");
248  dbe_->removeElement("convEtLoose");
249  dbe_->removeElement("convEtTight");
250 
251 
252  for(uint type=0;type!=types_.size();++type){
253 
254  for (int cut = 0; cut !=numberOfSteps_; ++cut) {
255 
256  currentFolder_.str("");
257  currentFolder_ << "Egamma/PhotonAnalyzer/" << types_[type] << "Photons/Et above " << (cut+1)*cutStep_ << " GeV/";
258 
259  //making bad channel histograms
260 
261  //vs Et
262  dividend = retrieveHisto(currentFolder_.str(),"badChannelsFractionVsEt");
263  numerator = retrieveHisto(currentFolder_.str(),"phoEtBadChannels");
264  denominator = retrieveHisto(currentFolder_.str(),"phoEtAllEcal");
265  dividePlots(dividend,numerator,denominator);
266 
267  //vs eta
268  dividend = retrieveHisto(currentFolder_.str(),"badChannelsFractionVsEta");
269  numerator = retrieveHisto(currentFolder_.str(),"phoEtaBadChannels");
270  denominator = retrieveHisto(currentFolder_.str(),"phoEta");
271  dividePlots(dividend,numerator,denominator);
272 
273  //vs phi
274  dividend = retrieveHisto(currentFolder_.str(),"badChannelsFractionVsPhi");
275  numerator = retrieveHisto(currentFolder_.str(),"phoPhiBadChannels");
276  denominator = retrieveHisto(currentFolder_.str(),"phoPhiAllEcal");
277  dividePlots(dividend,numerator,denominator);
278 
279  //making conversion fraction histograms
280 
281  //vs Et
282  dividend = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsEt");
283  numerator = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvEtAllEcal");
284  denominator = retrieveHisto(currentFolder_.str(),"phoEtAllEcal");
285  dividePlots(dividend,numerator,denominator);
286 
287  //vs eta
288  dividend = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsEta");
289  numerator = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvEtaForEfficiency");
290  denominator = retrieveHisto(currentFolder_.str(),"phoEta");
291  dividePlots(dividend,numerator,denominator);
292 
293  //vs phi
294  dividend = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsPhiAllEcal");
295  numerator = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyAllEcal");
296  denominator = retrieveHisto(currentFolder_.str(),"phoPhiAllEcal");
297  dividePlots(dividend,numerator,denominator);
298  dividend = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsPhiBarrel");
299  numerator = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyBarrel");
300  denominator = retrieveHisto(currentFolder_.str(),"phoPhiBarrel");
301  dividePlots(dividend,numerator,denominator);
302  dividend = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsPhiEndcaps");
303  numerator = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyEndcaps");
304  denominator = retrieveHisto(currentFolder_.str(),"phoPhiEndcaps");
305  dividePlots(dividend,numerator,denominator);
306 
307 
308  dbe_->setCurrentFolder(currentFolder_.str()+"Conversions/");
309  dbe_->removeElement("phoConvEtaForEfficiency");
310  dbe_->removeElement("phoConvPhiForEfficiencyAllEcal");
311  dbe_->removeElement("phoConvPhiForEfficiencyBarrel");
312  dbe_->removeElement("phoConvPhiForEfficiencyEndcaps");
313 
314  }
315 
316 
317  }
318 
319 
320  if(standAlone_) dbe_->save(outputFileName_);
321  else if(batch_) dbe_->save(inputFileName_);
322 
323 
324 }
325 
326 
327 
328 
329 void PhotonOfflineClient::dividePlots(MonitorElement* dividend, MonitorElement* numerator, MonitorElement* denominator){
330  double value,err;
331 
332  dividend->setEfficiencyFlag();
333  if(denominator->getEntries()==0) return;
334 
335  for (int j=1; j<=numerator->getNbinsX(); j++){
336  if (denominator->getBinContent(j)!=0){
337  value = ((double) numerator->getBinContent(j))/((double) denominator->getBinContent(j));
338  err = sqrt( value*(1-value) / ((double) denominator->getBinContent(j)) );
339  dividend->setBinContent(j, value);
340  dividend->setBinError(j,err);
341  }
342  else {
343  dividend->setBinContent(j, 0);
344  dividend->setBinError(j,0);
345  }
346  dividend->setEntries(numerator->getEntries());
347  }
348 }
349 
350 
351 void PhotonOfflineClient::dividePlots(MonitorElement* dividend, MonitorElement* numerator, double denominator){
352  double value,err;
353 
354  dividend->setEfficiencyFlag();
355  for (int j=1; j<=numerator->getNbinsX(); j++){
356  if (denominator!=0){
357  value = ((double) numerator->getBinContent(j))/denominator;
358  err = sqrt( value*(1-value) / denominator);
359  dividend->setBinContent(j, value);
360  dividend->setBinError(j,err);
361  }
362  else {
363  dividend->setBinContent(j, 0);
364  }
365  }
366 
367 }
368 
369 MonitorElement* PhotonOfflineClient::bookHisto(string histoName, string title, int bin, double min, double max)
370 {
371 
372  int histo_index = 0;
373  stringstream histo_number_stream;
374 
375 
376  //determining which folder we're in
377  if(dbe_->pwd().find( "InvMass" ) != string::npos){
378  histo_index_invMass_++;
379  histo_index = histo_index_invMass_;
380  }
381  if(dbe_->pwd().find( "Efficiencies" ) != string::npos){
382  histo_index_efficiency_++;
383  histo_index = histo_index_efficiency_;
384  }
385  histo_number_stream << "h_";
386  if(histo_index<10) histo_number_stream << "0";
387  histo_number_stream << histo_index;
388 
389  return dbe_->book1D(histo_number_stream.str()+"_"+histoName,title,bin,min,max);
390 
391 }
392 
393 vector<vector<MonitorElement*> > PhotonOfflineClient::book2DHistoVector(string histoType, string histoName, string title,
394  int xbin, double xmin,double xmax,
395  int ybin, double ymin, double ymax)
396 {
397  int histo_index = 0;
398 
399  vector<MonitorElement*> temp1DVector;
400  vector<vector<MonitorElement*> > temp2DVector;
401 
402  //determining which folder we're in
403  bool conversionPlot = false;
404  if(dbe_->pwd().find( "Conversions" ) != string::npos) conversionPlot = true;
405 
406 
407  if(conversionPlot){
408  histo_index_conversions_++;
409  histo_index = histo_index_conversions_;
410  }
411  else{
412  histo_index_photons_++;
413  histo_index = histo_index_photons_;
414  }
415 
416 
417  stringstream histo_number_stream;
418  histo_number_stream << "h_";
419  if(histo_index<10) histo_number_stream << "0";
420  histo_number_stream << histo_index << "_";
421 
422 
423  for(int cut = 0; cut != numberOfSteps_; ++cut){ //looping over Et cut values
424 
425  for(uint type=0;type!=types_.size();++type){ //looping over isolation type
426 
427  currentFolder_.str("");
428  currentFolder_ << "Egamma/PhotonAnalyzer/" << types_[type] << "Photons/Et above " << (cut+1)*cutStep_ << " GeV";
429  if(conversionPlot) currentFolder_ << "/Conversions";
430 
431  dbe_->setCurrentFolder(currentFolder_.str());
432 
433  string kind;
434  if(conversionPlot) kind = " Conversions: ";
435  else kind = " Photons: ";
436 
437  if(histoType=="1D") temp1DVector.push_back(dbe_->book1D(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax));
438  else if(histoType=="2D") temp1DVector.push_back(dbe_->book2D(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax,ybin,ymin,ymax));
439  else if(histoType=="Profile") temp1DVector.push_back(dbe_->bookProfile(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax,ybin,ymin,ymax,""));
440  //else cout << "bad histoType\n";
441  }
442 
443  temp2DVector.push_back(temp1DVector);
444  temp1DVector.clear();
445  }
446 
447  return temp2DVector;
448 
449 }
450 
451 
452 
453 vector<vector<vector<MonitorElement*> > > PhotonOfflineClient::book3DHistoVector(string histoType, string histoName, string title,
454  int xbin, double xmin,double xmax,
455  int ybin, double ymin, double ymax)
456 {
457  int histo_index = 0;
458 
459  vector<MonitorElement*> temp1DVector;
460  vector<vector<MonitorElement*> > temp2DVector;
461  vector<vector<vector<MonitorElement*> > > temp3DVector;
462 
463 
464  //determining which folder we're in
465  bool conversionPlot = false;
466  if(dbe_->pwd().find( "Conversions" ) != string::npos) conversionPlot = true;
467 
468 
469  if(conversionPlot){
470  histo_index_conversions_++;
471  histo_index = histo_index_conversions_;
472  }
473  else{
474  histo_index_photons_++;
475  histo_index = histo_index_photons_;
476  }
477 
478  stringstream histo_number_stream;
479  histo_number_stream << "h_";
480  if(histo_index<10) histo_number_stream << "0";
481  histo_number_stream << histo_index << "_";
482 
483 
484  for(int cut = 0; cut != numberOfSteps_; ++cut){ //looping over Et cut values
485 
486  for(uint type=0;type!=types_.size();++type){ //looping over isolation type
487 
488  for(uint part=0;part!=parts_.size();++part){ //looping over different parts of the ecal
489 
490  currentFolder_.str("");
491  currentFolder_ << "Egamma/PhotonAnalyzer/" << types_[type] << "Photons/Et above " << (cut+1)*cutStep_ << " GeV";
492  if(conversionPlot) currentFolder_ << "/Conversions";
493 
494  dbe_->setCurrentFolder(currentFolder_.str());
495 
496  string kind;
497  if(conversionPlot) kind = " Conversions: ";
498  else kind = " Photons: ";
499 
500  if(histoType=="1D") temp1DVector.push_back(dbe_->book1D( histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax));
501  else if(histoType=="2D") temp1DVector.push_back(dbe_->book2D( histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax,ybin,ymin,ymax));
502  else if(histoType=="Profile") temp1DVector.push_back(dbe_->bookProfile( histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax,ybin,ymin,ymax,""));
503  //else cout << "bad histoType\n";
504 
505 
506  }
507 
508  temp2DVector.push_back(temp1DVector);
509  temp1DVector.clear();
510  }
511 
512  temp3DVector.push_back(temp2DVector);
513  temp2DVector.clear();
514  }
515 
516  return temp3DVector;
517 }
518 
519 
520 MonitorElement* PhotonOfflineClient::retrieveHisto(string dir, string name){
521  //cout << "dir = " << dir << endl;
522  //cout << "name = " << name << endl;
523  vector<MonitorElement*> histoVector;
524  uint indexOfRelevantHistogram=0;
525  string fullMEName = "";
526  histoVector = dbe_->getContents(dir);
527  for(uint index=0;index!=histoVector.size();index++){
528  string MEName = histoVector[index]->getName();
529  if(MEName.find( name ) != string::npos){
530  indexOfRelevantHistogram = index;
531  break;
532  }
533  }
534  return histoVector[indexOfRelevantHistogram];
535 }
536 
type
type
Definition: HCALResponse.h:21
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:813
cuy.numerator
list numerator
Definition: cuy.py:483
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
PhotonOfflineClient::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: PhotonOfflineClient.cc:73
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:2118
PhotonOfflineClient::PhotonOfflineClient
PhotonOfflineClient(const edm::ParameterSet &pset)
Definition: PhotonOfflineClient.cc:26
jptDQMConfig_cff.etMax
int etMax
Definition: jptDQMConfig_cff.py:19
min
#define min(a, b)
Definition: mlp_lapack.h:161
jptDQMConfig_cff.phiMax
float phiMax
Definition: jptDQMConfig_cff.py:7
MonitorElement::getEntries
double getEntries(void) const
get # of entries
Definition: MonitorElement.cc:748
SiStripMonitorClusterAlca_cfi.ymax
tuple ymax
Definition: SiStripMonitorClusterAlca_cfi.py:37
SiStripMonitorClusterAlca_cfi.xmin
tuple xmin
Definition: SiStripMonitorClusterAlca_cfi.py:34
PhotonOfflineClient::endRun
virtual void endRun(const edm::Run &, const edm::EventSetup &)
Definition: PhotonOfflineClient.cc:84
PhotonOfflineClient::book3DHistoVector
std::vector< std::vector< std::vector< MonitorElement * > > > book3DHistoVector(std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)
Definition: PhotonOfflineClient.cc:453
PhotonOfflineClient::book2DHistoVector
std::vector< std::vector< MonitorElement * > > book2DHistoVector(std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)
Definition: PhotonOfflineClient.cc:393
cuy.denominator
list denominator
Definition: cuy.py:484
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2577
PhotonOfflineClient::dividePlots
void dividePlots(MonitorElement *dividend, MonitorElement *numerator, MonitorElement *denominator)
Definition: PhotonOfflineClient.cc:329
j
int j
Definition: DBlmapReader.cc:9
MonitorElement::setBinError
void setBinError(int binx, double error)
set uncertainty on content of bin (1-D)
Definition: MonitorElement.cc:838
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1036
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:398
MonitorElement::setEntries
void setEntries(double nentries)
set # of entries
Definition: MonitorElement.cc:880
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1473
SiStripMonitorClusterAlca_cfi.ymin
tuple ymin
Definition: SiStripMonitorClusterAlca_cfi.py:36
DQMStore::getContents
std::vector< MonitorElement * > getContents(const std::string &path) const
Definition: DQMStore.cc:1502
dbe_
DQMStore * dbe_
Definition: PFJetBenchmarkAnalyzer.cc:78
reco::tau::qcuts::etMin
bool etMin(const PFCandidate &cand, double cut)
Definition: RecoTauQualityCuts.cc:36
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:498
PhotonOfflineClient::retrieveHisto
MonitorElement * retrieveHisto(std::string dir, std::string name)
Definition: PhotonOfflineClient.cc:520
jptDQMConfig_cff.etaMin
int etaMin
Definition: jptDQMConfig_cff.py:52
PhotonOfflineClient::bookHisto
MonitorElement * bookHisto(std::string histoName, std::string title, int bin, double min, double max)
Definition: PhotonOfflineClient.cc:369
part
part
Definition: HCALResponse.h:20
MonitorElement::getIntValue
int64_t getIntValue(void) const
Definition: MonitorElement.h:333
SiStripMonitorClusterAlca_cfi.xmax
tuple xmax
Definition: SiStripMonitorClusterAlca_cfi.py:35
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:712
jptDQMConfig_cff.etaMax
float etaMax
Definition: jptDQMConfig_cff.py:9
PhotonOfflineClient::endLuminosityBlock
virtual void endLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &)
Definition: PhotonOfflineClient.cc:75
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:694
jptDQMConfig_cff.phiMin
int phiMin
Definition: jptDQMConfig_cff.py:50
DQMStore::open
bool open(const std::string &filename, bool overwrite=false, const std::string &path="", const std::string &prepend="", OpenRunDirs stripdirs=KeepRunDirs, bool fileMustExist=true)
Definition: DQMStore.cc:2437
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
MonitorElement::setEfficiencyFlag
void setEfficiencyFlag(void)
Definition: MonitorElement.h:141
dir
dbl *** dir
Definition: mlp_gen.cc:35
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:850
HcalObjRepresent::setup
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
Definition: HcalObjRepresent.h:628
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
edm::Run
Definition: Run.h:36
GOODCOLL_filter_cfg.cut
tuple cut
Definition: GOODCOLL_filter_cfg.py:65
DQMStore::pwd
const std::string & pwd(void) const
Definition: DQMStore.cc:406
Generated for CMSSW Reference Manual by doxygen 1.8.5