CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
HLTTauDQMSummaryPlotter Class Reference

#include <HLTTauDQMSummaryPlotter.h>

Inheritance diagram for HLTTauDQMSummaryPlotter:
HLTTauDQMPlotter

Public Member Functions

 HLTTauDQMSummaryPlotter (const edm::ParameterSet &, std::string)
 
const std::string name ()
 
void plot ()
 
 ~HLTTauDQMSummaryPlotter ()
 
- Public Member Functions inherited from HLTTauDQMPlotter
 HLTTauDQMPlotter ()
 
bool isValid () const
 
virtual ~HLTTauDQMPlotter ()
 

Private Member Functions

void bookEfficiencyHisto (std::string folder, std::string name, std::string hist1)
 
void bookTriggerBitEfficiencyHistos (std::string folder, std::string histo)
 
std::pair< double, double > calcEfficiency (float num, float denom)
 
void plotEfficiencyHisto (std::string folder, std::string name, std::string hist1, std::string hist2)
 
void plotIntegratedEffHisto (std::string folder, std::string name, std::string refHisto, std::string evCount, int bin)
 
void plotTriggerBitEfficiencyHistos (std::string folder, std::string histo)
 

Private Attributes

std::string type_
 

Additional Inherited Members

- Protected Member Functions inherited from HLTTauDQMPlotter
std::pair< bool, LV > match (const LV &, const LVColl &, double)
 
std::string triggerTag ()
 
- Protected Attributes inherited from HLTTauDQMPlotter
std::string dqmBaseFolder_
 
std::string name_
 
DQMStore * store_
 
std::string triggerTag_
 
std::string triggerTagAlias_
 
bool validity_
 

Detailed Description

Definition at line 6 of file HLTTauDQMSummaryPlotter.h.

Constructor & Destructor Documentation

HLTTauDQMSummaryPlotter::HLTTauDQMSummaryPlotter ( const edm::ParameterSet &  ps,
std::string  dqmBaseFolder 
)

Definition at line 3 of file HLTTauDQMSummaryPlotter.cc.

References bookEfficiencyHisto(), bookTriggerBitEfficiencyHistos(), HLTTauDQMPlotter::dqmBaseFolder_, alignCSCRings::e, edm::ParameterSet::getUntrackedParameter(), HLTTauDQMPlotter::name_, HLTTauDQMPlotter::store_, HLTTauDQMPlotter::triggerTag(), HLTTauDQMPlotter::triggerTag_, HLTTauDQMPlotter::triggerTagAlias_, type_, HLTTauDQMPlotter::validity_, and cms::Exception::what().

3  {
4  //Initialize Plotter
5  name_ = "HLTTauDQMSummaryPlotter";
6 
7  //Process PSet
8  try {
9  triggerTag_ = ps.getUntrackedParameter<std::string>("DQMFolder");
10  triggerTagAlias_ = ps.getUntrackedParameter<std::string>("Alias","");
11  type_ = ps.getUntrackedParameter<std::string>("ConfigType");
12  dqmBaseFolder_ = dqmBaseFolder;
13  validity_ = true;
14  } catch ( cms::Exception &e ) {
15  edm::LogInfo("HLTTauDQMSummaryPlotter::HLTTauDQMSummaryPlotter") << e.what() << std::endl;
16  validity_ = false;
17  return;
18  }
19 
20  if (store_) {
21  //Path Summary
22  if ( type_ == "Path" ) {
23  bookTriggerBitEfficiencyHistos(triggerTag(),"MatchedTriggerBits");
24  }
25 
26  //Lite Path Summary
27  if ( type_ == "LitePath" ) {
28  bookEfficiencyHisto(triggerTag(),"PathEfficiency","MatchedPathTriggerBits");
29  bookEfficiencyHisto(triggerTag(),"TrigTauEtEff","EfficiencyHelpers/TrigTauEtEffNum");
30  bookEfficiencyHisto(triggerTag(),"TrigTauEtaEff","EfficiencyHelpers/TrigTauEtaEffNum");
31  bookEfficiencyHisto(triggerTag(),"TrigTauPhiEff","EfficiencyHelpers/TrigTauPhiEffNum");
32  }
33 
34  //L1 Summary
35  if ( type_ == "L1" ) {
36  bookEfficiencyHisto(triggerTag(),"L1TauEtEff","EfficiencyHelpers/L1TauEtEffNum");
37  bookEfficiencyHisto(triggerTag(),"L1TauEtaEff","EfficiencyHelpers/L1TauEtaEffNum");
38  bookEfficiencyHisto(triggerTag(),"L1TauPhiEff","EfficiencyHelpers/L1TauPhiEffNum");
39 
40  bookEfficiencyHisto(triggerTag(),"L1JetEtEff","EfficiencyHelpers/L1JetEtEffNum");
41  bookEfficiencyHisto(triggerTag(),"L1JetEtaEff","EfficiencyHelpers/L1JetEtaEffNum");
42  bookEfficiencyHisto(triggerTag(),"L1JetPhiEff","EfficiencyHelpers/L1JetPhiEffNum");
43 
44  bookEfficiencyHisto(triggerTag(),"L1SingleTauEff","L1LeadTauEt");
45  bookEfficiencyHisto(triggerTag(),"L1DoubleTauEff","L1SecondTauEt");
46  }
47 
48  //L2 Summary
49  if ( type_ == "Calo" ) {
50  bookEfficiencyHisto(triggerTag(),"L2RecoTauEtEff","EfficiencyHelpers/L2RecoTauEtEffNum");
51  bookEfficiencyHisto(triggerTag(),"L2RecoTauEtaEff","EfficiencyHelpers/L2RecoTauEtaEffNum");
52  bookEfficiencyHisto(triggerTag(),"L2RecoTauPhiEff","EfficiencyHelpers/L2RecoTauPhiEffNum");
53 
54  bookEfficiencyHisto(triggerTag(),"L2IsoTauEtEff","EfficiencyHelpers/L2IsoTauEtEffNum");
55  bookEfficiencyHisto(triggerTag(),"L2IsoTauEtaEff","EfficiencyHelpers/L2IsoTauEtaEffNum");
56  bookEfficiencyHisto(triggerTag(),"L2IsoTauPhiEff","EfficiencyHelpers/L2IsoTauPhiEffNum");
57  }
58 
59  //L25/3 Summary
60  if ( type_ == "Track" ) {
61  bookEfficiencyHisto(triggerTag(),"L25TauEtEff","L25TauEtEffNum");
62  bookEfficiencyHisto(triggerTag(),"L25TauEtaEff","L25TauEtaEffNum");
63  bookEfficiencyHisto(triggerTag(),"L25TauPhiEff","L25TauPhiEffNum");
64  bookEfficiencyHisto(triggerTag(),"L3TauEtEff","L3TauEtEffNum");
65  bookEfficiencyHisto(triggerTag(),"L3TauEtaEff","L3TauEtaEffNum");
66  bookEfficiencyHisto(triggerTag(),"L3TauPhiEff","L3TauPhiEffNum");
67  }
68  }
69 }
cms::Exception::what
virtual char const * what() const
Definition: Exception.cc:141
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HLTTauDQMSummaryPlotter::bookTriggerBitEfficiencyHistos
void bookTriggerBitEfficiencyHistos(std::string folder, std::string histo)
Definition: HLTTauDQMSummaryPlotter.cc:202
HLTTauDQMPlotter::triggerTagAlias_
std::string triggerTagAlias_
Definition: HLTTauDQMPlotter.h:73
HLTTauDQMPlotter::triggerTag
std::string triggerTag()
Definition: HLTTauDQMPlotter.cc:26
HLTTauDQMSummaryPlotter::bookEfficiencyHisto
void bookEfficiencyHisto(std::string folder, std::string name, std::string hist1)
Definition: HLTTauDQMSummaryPlotter.cc:134
HLTTauDQMPlotter::triggerTag_
std::string triggerTag_
Definition: HLTTauDQMPlotter.h:72
HLTTauDQMPlotter::dqmBaseFolder_
std::string dqmBaseFolder_
Definition: HLTTauDQMPlotter.h:71
HLTTauDQMSummaryPlotter::~HLTTauDQMSummaryPlotter ( )

Definition at line 71 of file HLTTauDQMSummaryPlotter.cc.

71  {
72 }

Member Function Documentation

void HLTTauDQMSummaryPlotter::bookEfficiencyHisto ( std::string  folder,
std::string  name,
std::string  hist1 
)
private

Definition at line 134 of file HLTTauDQMSummaryPlotter.cc.

References DQMStore::bookProfile(), DQMStore::dirExists(), DQMStore::get(), MonitorElement::getTH1F(), DQMStore::setCurrentFolder(), MonitorElement::setTitle(), HLTTauDQMPlotter::store_, and tmp.

Referenced by HLTTauDQMSummaryPlotter().

134  {
135  if ( store_->dirExists(folder) ) {
136  store_->setCurrentFolder(folder);
137 
138  MonitorElement * effnum = store_->get(folder+"/"+hist1);
139 
140  if ( effnum ) {
141  MonitorElement *tmp = store_->bookProfile(name,name,effnum->getTH1F()->GetNbinsX(),effnum->getTH1F()->GetXaxis()->GetXmin(),effnum->getTH1F()->GetXaxis()->GetXmax(),105,0,1.05);
142 
143  tmp->setTitle(name);
144  }
145  }
146 }
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:1031
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:1468
MonitorElement::setTitle
void setTitle(const std::string &title)
set (ie. change) histogram/profile title
Definition: MonitorElement.cc:952
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:493
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1486
tmp
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
void HLTTauDQMSummaryPlotter::bookTriggerBitEfficiencyHistos ( std::string  folder,
std::string  histo 
)
private

Definition at line 202 of file HLTTauDQMSummaryPlotter.cc.

References DQMStore::bookProfile(), DQMStore::dirExists(), interpolateCardsSimple::eff, DQMStore::get(), MonitorElement::getNbinsX(), DQMStore::setCurrentFolder(), and HLTTauDQMPlotter::store_.

Referenced by HLTTauDQMSummaryPlotter().

202  {
203  if ( store_->dirExists(folder) ) {
204  store_->setCurrentFolder(folder);
205 
206  MonitorElement * eff = store_->get(folder+"/"+histo);
207 
208  if ( eff ) {
209  store_->bookProfile("EfficiencyRefInput","Efficiency with Matching",eff->getNbinsX()-1,0,eff->getNbinsX()-1,100,0,1);
210  store_->bookProfile("EfficiencyRefL1","Efficiency with Matching Ref to L1",eff->getNbinsX()-2,0,eff->getNbinsX()-2,100,0,1);
211  store_->bookProfile("EfficiencyRefPrevious","Efficiency with Matching Ref to Previous",eff->getNbinsX()-1,0,eff->getNbinsX()-1,100,0,1);
212  }
213  }
214 }
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:1031
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:1468
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:493
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:694
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
std::pair< double, double > HLTTauDQMSummaryPlotter::calcEfficiency ( float  num,
float  denom 
)
private

Definition at line 260 of file HLTTauDQMSummaryPlotter.cc.

References mathSSE::sqrt().

Referenced by plotEfficiencyHisto(), and plotTriggerBitEfficiencyHistos().

260  {
261  if ( denom != 0.0 ) {
262  return std::pair<double,double>(num/denom,sqrt(num/denom*(1.0-num/denom)/denom));
263  }
264  return std::pair<double,double>(0.0,0.0);
265 }
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
evf::utils::num
long long int num
Definition: procUtils.cc:71
const std::string HLTTauDQMSummaryPlotter::name ( void  )
inlinevirtual

Implements HLTTauDQMPlotter.

Definition at line 10 of file HLTTauDQMSummaryPlotter.h.

References HLTTauDQMPlotter::name_.

Referenced by BeautifulSoup.Tag::_invert(), and Vispa.Views.PropertyView.Property::valueChanged().

10 { return name_; }
void HLTTauDQMSummaryPlotter::plot ( )

Definition at line 74 of file HLTTauDQMSummaryPlotter.cc.

References plotEfficiencyHisto(), plotIntegratedEffHisto(), plotTriggerBitEfficiencyHistos(), HLTTauDQMPlotter::store_, HLTTauDQMPlotter::triggerTag(), and type_.

74  {
75  if (store_) {
76  //Path Summary
77  if ( type_ == "Path" ) {
78  plotTriggerBitEfficiencyHistos(triggerTag(),"MatchedTriggerBits");
79  }
80 
81  //Lite Path Summary
82  if ( type_ == "LitePath" ) {
83  plotEfficiencyHisto(triggerTag(),"PathEfficiency","MatchedPathTriggerBits","RefEvents");
84  plotEfficiencyHisto(triggerTag(),"TrigTauEtEff","EfficiencyHelpers/TrigTauEtEffNum","EfficiencyHelpers/TrigTauEtEffDenom");
85  plotEfficiencyHisto(triggerTag(),"TrigTauEtaEff","EfficiencyHelpers/TrigTauEtaEffNum","EfficiencyHelpers/TrigTauEtaEffDenom");
86  plotEfficiencyHisto(triggerTag(),"TrigTauPhiEff","EfficiencyHelpers/TrigTauPhiEffNum","EfficiencyHelpers/TrigTauPhiEffDenom");
87  }
88 
89  //L1 Summary
90  if ( type_ == "L1" ) {
91  plotEfficiencyHisto(triggerTag(),"L1TauEtEff","EfficiencyHelpers/L1TauEtEffNum","EfficiencyHelpers/L1TauEtEffDenom");
92  plotEfficiencyHisto(triggerTag(),"L1TauEtaEff","EfficiencyHelpers/L1TauEtaEffNum","EfficiencyHelpers/L1TauEtaEffDenom");
93  plotEfficiencyHisto(triggerTag(),"L1TauPhiEff","EfficiencyHelpers/L1TauPhiEffNum","EfficiencyHelpers/L1TauPhiEffDenom");
94 
95  plotEfficiencyHisto(triggerTag(),"L1JetEtEff","EfficiencyHelpers/L1JetEtEffNum","EfficiencyHelpers/L1JetEtEffDenom");
96  plotEfficiencyHisto(triggerTag(),"L1JetEtaEff","EfficiencyHelpers/L1JetEtaEffNum","EfficiencyHelpers/L1JetEtaEffDenom");
97  plotEfficiencyHisto(triggerTag(),"L1JetPhiEff","EfficiencyHelpers/L1JetPhiEffNum","EfficiencyHelpers/L1JetPhiEffDenom");
98 
99  plotEfficiencyHisto(triggerTag(),"L1ElectronEtEff","EfficiencyHelpers/L1ElectronEtEffNum","EfficiencyHelpers/L1ElectronEtEffDenom");
100  plotEfficiencyHisto(triggerTag(),"L1ElectronEtaEff","EfficiencyHelpers/L1ElectronEtaEffNum","EfficiencyHelpers/L1ElectronEtaEffDenom");
101  plotEfficiencyHisto(triggerTag(),"L1ElectronPhiEff","EfficiencyHelpers/L1ElectronPhiEffNum","EfficiencyHelpers/L1ElectronPhiEffDenom");
102 
103  plotEfficiencyHisto(triggerTag(),"L1MuonEtEff","EfficiencyHelpers/L1MuonEtEffNum","EfficiencyHelpers/L1MuonEtEffDenom");
104  plotEfficiencyHisto(triggerTag(),"L1MuonEtaEff","EfficiencyHelpers/L1MuonEtaEffNum","EfficiencyHelpers/L1MuonEtaEffDenom");
105  plotEfficiencyHisto(triggerTag(),"L1MuonPhiEff","EfficiencyHelpers/L1MuonPhiEffNum","EfficiencyHelpers/L1MuonPhiEffDenom");
106 
107  plotIntegratedEffHisto(triggerTag(),"L1SingleTauEff","L1LeadTauEt","InputEvents",1);
108  plotIntegratedEffHisto(triggerTag(),"L1DoubleTauEff","L1SecondTauEt","InputEvents",2);
109  }
110 
111  //L2 Summary
112  if ( type_ == "Calo" ) {
113  plotEfficiencyHisto(triggerTag(),"L2RecoTauEtEff","EfficiencyHelpers/L2RecoTauEtEffNum","EfficiencyHelpers/L2RecoTauEtEffDenom");
114  plotEfficiencyHisto(triggerTag(),"L2RecoTauEtaEff","EfficiencyHelpers/L2RecoTauEtaEffNum","EfficiencyHelpers/L2RecoTauEtaEffDenom");
115  plotEfficiencyHisto(triggerTag(),"L2RecoTauPhiEff","EfficiencyHelpers/L2RecoTauPhiEffNum","EfficiencyHelpers/L2RecoTauPhiEffDenom");
116 
117  plotEfficiencyHisto(triggerTag(),"L2IsoTauEtEff","EfficiencyHelpers/L2IsoTauEtEffNum","EfficiencyHelpers/L2IsoTauEtEffDenom");
118  plotEfficiencyHisto(triggerTag(),"L2IsoTauEtaEff","EfficiencyHelpers/L2IsoTauEtaEffNum","EfficiencyHelpers/L2IsoTauEtaEffDenom");
119  plotEfficiencyHisto(triggerTag(),"L2IsoTauPhiEff","EfficiencyHelpers/L2IsoTauPhiEffNum","EfficiencyHelpers/L2IsoTauPhiEffDenom");
120  }
121 
122  //L25/3 Summary
123  if ( type_ == "Track" ) {
124  plotEfficiencyHisto(triggerTag(),"L25TauEtEff","L25TauEtEffNum","L25TauEtEffDenom");
125  plotEfficiencyHisto(triggerTag(),"L25TauEtaEff","L25TauEtaEffNum","L25TauEtaEffDenom");
126  plotEfficiencyHisto(triggerTag(),"L25TauPhiEff","L25TauPhiEffNum","L25TauPhiEffDenom");
127  plotEfficiencyHisto(triggerTag(),"L3TauEtEff","L3TauEtEffNum","L3TauEtEffDenom");
128  plotEfficiencyHisto(triggerTag(),"L3TauEtaEff","L3TauEtaEffNum","L3TauEtaEffDenom");
129  plotEfficiencyHisto(triggerTag(),"L3TauPhiEff","L3TauPhiEffNum","L3TauPhiEffDenom");
130  }
131  }
132 }
HLTTauDQMSummaryPlotter::plotTriggerBitEfficiencyHistos
void plotTriggerBitEfficiencyHistos(std::string folder, std::string histo)
Definition: HLTTauDQMSummaryPlotter.cc:216
HLTTauDQMSummaryPlotter::plotEfficiencyHisto
void plotEfficiencyHisto(std::string folder, std::string name, std::string hist1, std::string hist2)
Definition: HLTTauDQMSummaryPlotter.cc:148
HLTTauDQMPlotter::triggerTag
std::string triggerTag()
Definition: HLTTauDQMPlotter.cc:26
HLTTauDQMSummaryPlotter::plotIntegratedEffHisto
void plotIntegratedEffHisto(std::string folder, std::string name, std::string refHisto, std::string evCount, int bin)
Definition: HLTTauDQMSummaryPlotter.cc:168
void HLTTauDQMSummaryPlotter::plotEfficiencyHisto ( std::string  folder,
std::string  name,
std::string  hist1,
std::string  hist2 
)
private

Definition at line 148 of file HLTTauDQMSummaryPlotter.cc.

References calcEfficiency(), DQMStore::dirExists(), interpolateCardsSimple::eff, postValidation_cfi::efficiency, first, DQMStore::get(), MonitorElement::getTH1F(), MonitorElement::getTProfile(), i, edm::second(), DQMStore::setCurrentFolder(), mathSSE::sqrt(), and HLTTauDQMPlotter::store_.

Referenced by plot().

148  {
149  if ( store_->dirExists(folder) ) {
150  store_->setCurrentFolder(folder);
151 
152  MonitorElement * effnum = store_->get(folder+"/"+hist1);
153  MonitorElement * effdenom = store_->get(folder+"/"+hist2);
154  MonitorElement * eff = store_->get(folder+"/"+name);
155 
156  if ( effnum && effdenom && eff ) {
157  for ( int i = 1; i <= effnum->getTH1F()->GetNbinsX(); ++i ) {
158  double efficiency = calcEfficiency(effnum->getTH1F()->GetBinContent(i),effdenom->getTH1F()->GetBinContent(i)).first;
159  double err = calcEfficiency(effnum->getTH1F()->GetBinContent(i),effdenom->getTH1F()->GetBinContent(i)).second;
160  eff->getTProfile()->SetBinContent(i,efficiency);
161  eff->getTProfile()->SetBinEntries(i,1);
162  eff->getTProfile()->SetBinError(i,sqrt(efficiency*efficiency+err*err));
163  }
164  }
165  }
166 }
i
int i
Definition: DBlmapReader.cc:9
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:234
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
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:1468
first
bool first
Definition: L1TdeRCT.cc:94
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:493
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1486
MonitorElement::getTProfile
TProfile * getTProfile(void) const
Definition: MonitorElement.cc:1542
HLTTauDQMSummaryPlotter::calcEfficiency
std::pair< double, double > calcEfficiency(float num, float denom)
Definition: HLTTauDQMSummaryPlotter.cc:260
postValidation_cfi.efficiency
tuple efficiency
Definition: postValidation_cfi.py:73
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
void HLTTauDQMSummaryPlotter::plotIntegratedEffHisto ( std::string  folder,
std::string  name,
std::string  refHisto,
std::string  evCount,
int  bin 
)
private

Definition at line 168 of file HLTTauDQMSummaryPlotter.cc.

References DQMStore::dirExists(), interpolateCardsSimple::eff, DQMStore::get(), MonitorElement::getTH1F(), MonitorElement::getTProfile(), interpolateCardsSimple::histo, i, funct::integral(), pileupCalc::nbins, DQMStore::setCurrentFolder(), mathSSE::sqrt(), and HLTTauDQMPlotter::store_.

Referenced by plot().

168  {
169  if ( store_->dirExists(folder) ) {
170  store_->setCurrentFolder(folder);
171 
172  MonitorElement * refH = store_->get(folder+"/"+refHisto);
173  MonitorElement * evC = store_->get(folder+"/"+evCount);
174  MonitorElement * eff = store_->get(folder+"/"+name);
175 
176  if ( refH && evC && eff ) {
177  TH1F *histo = refH->getTH1F();
178  float nGenerated = evC->getTH1F()->GetBinContent(bin);
179  // Assuming that the histogram is incremented with weight=1 for each event
180  // this function integrates the histogram contents above every bin and stores it
181  // in that bin. The result is plot of integral rate versus threshold plot.
182  int nbins = histo->GetNbinsX();
183  double integral = histo->GetBinContent(nbins+1); // Initialize to overflow
184  if (nGenerated<=0.0) nGenerated=1.0;
185  for ( int i = nbins; i >= 1; i-- ) {
186  double thisBin = histo->GetBinContent(i);
187  integral += thisBin;
188  double integralEff;
189  double integralError;
190  integralEff = (integral / nGenerated);
191  eff->getTProfile()->SetBinContent(i, integralEff);
192  eff->getTProfile()->SetBinEntries(i, 1);
193  // error
194  integralError = (sqrt(integral) / nGenerated);
195 
196  eff->getTProfile()->SetBinError(i, sqrt(integralEff*integralEff+integralError*integralError));
197  }
198  }
199  }
200 }
i
int i
Definition: DBlmapReader.cc:9
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
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:1468
funct::integral
Integral< F, X >::type integral(const F &f)
Definition: Integral.h:69
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:493
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1486
MonitorElement::getTProfile
TProfile * getTProfile(void) const
Definition: MonitorElement.cc:1542
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
void HLTTauDQMSummaryPlotter::plotTriggerBitEfficiencyHistos ( std::string  folder,
std::string  histo 
)
private

Definition at line 216 of file HLTTauDQMSummaryPlotter.cc.

References calcEfficiency(), DQMStore::dirExists(), interpolateCardsSimple::eff, postValidation_cfi::efficiency, first, DQMStore::get(), MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::getTH1F(), MonitorElement::getTProfile(), i, edm::second(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), mathSSE::sqrt(), and HLTTauDQMPlotter::store_.

Referenced by plot().

216  {
217  if ( store_->dirExists(folder) ) {
218  store_->setCurrentFolder(folder);
219  MonitorElement * eff = store_->get(folder+"/"+histo);
220  MonitorElement * effRefTruth = store_->get(folder+"/EfficiencyRefInput");
221  MonitorElement * effRefL1 = store_->get(folder+"/EfficiencyRefL1");
222  MonitorElement * effRefPrevious = store_->get(folder+"/EfficiencyRefPrevious");
223 
224  if ( eff && effRefTruth && effRefL1 && effRefPrevious ) {
225  //Calculate efficiencies with ref to matched objects
226  for ( int i = 2; i <= eff->getNbinsX(); ++i ) {
227  double efficiency = calcEfficiency(eff->getBinContent(i),eff->getBinContent(1)).first;
228  double err = calcEfficiency(eff->getBinContent(i),eff->getBinContent(1)).second;
229 
230  effRefTruth->getTProfile()->SetBinContent(i-1,efficiency);
231  effRefTruth->getTProfile()->SetBinEntries(i-1,1);
232  effRefTruth->getTProfile()->SetBinError(i-1,sqrt(efficiency*efficiency+err*err));
233  effRefTruth->setBinLabel(i-1,eff->getTH1F()->GetXaxis()->GetBinLabel(i));
234 
235  }
236  //Calculate efficiencies with ref to L1
237  for ( int i = 3; i <= eff->getNbinsX(); ++i ) {
238  double efficiency = calcEfficiency(eff->getBinContent(i),eff->getBinContent(2)).first;
239  double err = calcEfficiency(eff->getBinContent(i),eff->getBinContent(2)).second;
240 
241  effRefL1->getTProfile()->SetBinContent(i-2,efficiency);
242  effRefL1->getTProfile()->SetBinEntries(i-2,1);
243  effRefL1->getTProfile()->SetBinError(i-2,sqrt(efficiency*efficiency+err*err));
244  effRefL1->setBinLabel(i-2,eff->getTH1F()->GetXaxis()->GetBinLabel(i));
245  }
246  //Calculate efficiencies with ref to previous
247  for ( int i = 2; i <= eff->getNbinsX(); ++i ) {
248  double efficiency = calcEfficiency(eff->getBinContent(i),eff->getBinContent(i-1)).first;
249  double err = calcEfficiency(eff->getBinContent(i),eff->getBinContent(i-1)).second;
250 
251  effRefPrevious->getTProfile()->SetBinContent(i-1,efficiency);
252  effRefPrevious->getTProfile()->SetBinEntries(i-1,1);
253  effRefPrevious->getTProfile()->SetBinError(i-1,sqrt(efficiency*efficiency+err*err));
254  effRefPrevious->setBinLabel(i-1,eff->getTH1F()->GetXaxis()->GetBinLabel(i));
255  }
256  }
257  }
258 }
i
int i
Definition: DBlmapReader.cc:9
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:889
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:234
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
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:1468
first
bool first
Definition: L1TdeRCT.cc:94
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:493
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1486
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:712
MonitorElement::getTProfile
TProfile * getTProfile(void) const
Definition: MonitorElement.cc:1542
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:694
HLTTauDQMSummaryPlotter::calcEfficiency
std::pair< double, double > calcEfficiency(float num, float denom)
Definition: HLTTauDQMSummaryPlotter.cc:260
postValidation_cfi.efficiency
tuple efficiency
Definition: postValidation_cfi.py:73
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429

Member Data Documentation

std::string HLTTauDQMSummaryPlotter::type_
private

Definition at line 21 of file HLTTauDQMSummaryPlotter.h.

Referenced by Modules.Service::_placeImpl(), Modules.ESSource::_placeImpl(), Modules.ESProducer::_placeImpl(), Modules.ESPrefer::copy(), Mixins._TypedParameterizable::dumpPython(), Modules.ESPrefer::dumpPythonAs(), HLTTauDQMSummaryPlotter(), Modules.Service::insertInto(), RandomRunSource.RandomRunSource::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::moduleLabel_(), Modules.ESProducer::moduleLabel_(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), Modules.ESPrefer::nameInProcessDesc_(), and plot().

Generated for CMSSW Reference Manual by doxygen 1.8.5