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