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 Attributes
JetMETDQMPostProcessor Class Reference

#include <JetMETDQMPostProcessor.h>

Inheritance diagram for JetMETDQMPostProcessor:
edm::EDAnalyzer

Public Member Functions

void analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
 
void endJob ()
 
void endRun (edm::Run const &, edm::EventSetup const &)
 
 JetMETDQMPostProcessor (const edm::ParameterSet &pset)
 
 ~JetMETDQMPostProcessor ()
 
- Public Member Functions inherited from edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 

Private Attributes

std::string subDir_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzer > WorkerType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const * currentContext () const
 

Detailed Description

Definition at line 8 of file JetMETDQMPostProcessor.h.

Constructor & Destructor Documentation

JetMETDQMPostProcessor::JetMETDQMPostProcessor ( const edm::ParameterSet &  pset)

Definition at line 16 of file JetMETDQMPostProcessor.cc.

References edm::ParameterSet::getUntrackedParameter(), and subDir_.

17 {
18  subDir_ = pset.getUntrackedParameter<std::string>("subDir");
19 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
JetMETDQMPostProcessor::~JetMETDQMPostProcessor ( )
inline

Definition at line 11 of file JetMETDQMPostProcessor.h.

11 {};

Member Function Documentation

void JetMETDQMPostProcessor::analyze ( const edm::Event &  event,
const edm::EventSetup &  eventSetup 
)
inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 13 of file JetMETDQMPostProcessor.h.

13 {};
void JetMETDQMPostProcessor::endJob ( void  )
virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 435 of file JetMETDQMPostProcessor.cc.

436 {
437  /*
439  // setup DQM stor //
441 
442  DQMStore * dqm = 0;
443  dqm = edm::Service<DQMStore>().operator->();
444 
445  if ( ! dqm ) {
446  edm::LogInfo("JetMETDQMPostProcessor") << "Cannot create DQMStore instance\n";
447  return;
448  }
449 
450 
451  //go to the directory to be processed
452  if(dqm->dirExists(subDir_)) dqm->cd(subDir_);
453  else {
454  edm::LogWarning("JetMETDQMPostProcessor") << "cannot find directory: " << subDir_ << " , skipping";
455  return;
456  }
457 
458  std::vector<std::string> subdirectories = dqm->getSubdirs();
459  for(std::vector<std::string>::iterator dir = subdirectories.begin() ;dir!= subdirectories.end(); dir++ ){
460  dqm->cd(*dir);
461 
462  MonitorElement*_meTurnOnMET = dqm->book1D("_meTurnOnMET","Missing ET Turn-On",100,0,500);
463  MonitorElement*_meTurnOnJetPt = dqm->book1D("_meTurnOnJetPt","Jet Pt Turn-On",100,0,500);
464 
465  //std::vector<std::string> mes = dqm->getMEs();
466  //for(std::vector<std::string>::iterator me = mes.begin() ;me!= mes.end(); me++ )
467  // std::cout <<*me <<std::endl;
468  //std::cout <<std::endl;
469 
470  _meTurnOnMET->getTH1F()->Add(dqm->get(dqm->pwd() + "/_meGenMETTrg")->getTH1F(),1);
471  _meTurnOnMET->getTH1F()->Sumw2();
472  dqm->get(dqm->pwd() + "/_meGenMET")->getTH1F()->Sumw2();
473  _meTurnOnMET->getTH1F()->Divide(_meTurnOnMET->getTH1F(),dqm->get(dqm->pwd() + "/_meGenMET")->getTH1F(),1,1,"B");
474 
475  _meTurnOnJetPt->getTH1F()->Add(dqm->get(dqm->pwd() + "/_meGenJetPtTrg")->getTH1F(),1);
476  _meTurnOnJetPt->getTH1F()->Sumw2();
477  dqm->get(dqm->pwd() + "/_meGenJetPt")->getTH1F()->Sumw2();
478  _meTurnOnJetPt->getTH1F()->Divide(_meTurnOnJetPt->getTH1F(),dqm->get(dqm->pwd() + "/_meGenJetPt")->getTH1F(),1,1,"B");
479 
480  dqm->goUp();
481  }
482  */
483 }
void JetMETDQMPostProcessor::endRun ( edm::Run const &  run,
edm::EventSetup const &  es 
)
virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 21 of file JetMETDQMPostProcessor.cc.

References DQMStore::bookProfile(), DQMStore::cd(), dir, DQMStore::dirExists(), DQMStore::get(), DQMStore::getSubdirs(), getTH1F(), ght, DQMStore::goUp(), cmsCodeRules.cppFunctionSkipper::operator, DQMStore::pwd(), mathSSE::sqrt(), and subDir_.

22 {
24  // setup DQM stor //
26 
27  DQMStore * dqm = 0;
28  dqm = edm::Service<DQMStore>().operator->();
29 
30  if ( ! dqm ) {
31  edm::LogInfo("JetMETDQMPostProcessor") << "Cannot create DQMStore instance\n";
32  return;
33  }
34 
35 
36  //go to the directory to be processed
37  if(dqm->dirExists(subDir_)) dqm->cd(subDir_);
38  else {
39  edm::LogWarning("JetMETDQMPostProcessor") << "cannot find directory: " << subDir_ << " , skipping";
40  return;
41  }
42 
43  std::vector<std::string> subdirectories = dqm->getSubdirs();
44  for(std::vector<std::string>::iterator dir = subdirectories.begin() ;dir!= subdirectories.end(); dir++ ){
45  dqm->cd(*dir);
46 
47  TH1F* gmetTrg = new TH1F("gmetTrg","gmetTrg",100,0,500);
48  TH1F* gmetTrg2 = new TH1F("gmetTrg2","gmetTrg2",100,0,500);
49  TH1F* gmetTrgLow = new TH1F("gmetTrgLow","gmetTrgLow",100,0,500);
50  TH1F* gmet = new TH1F("gmet","gmet",100,0,500);
51  TH1F* gjetTrg = new TH1F("gjetTrg","gjetTrg",100,0,500);
52  TH1F* gjetTrg2 = new TH1F("gjetTrg2","gjetTrg2",100,0,500);
53  TH1F* gjetTrgLow = new TH1F("gjetTrgLow","gjetTrgLow",100,0,500);
54  TH1F* gjet = new TH1F("gjet","gjet",100,0,500);
55  TH1F* gjetEtaTrg = new TH1F("gjetEtaTrg","gjetEtaTrg",100,-10,10);
56  TH1F* gjetEtaTrg2 = new TH1F("gjetEtaTrg2","gjetEtaTrg2",100,-10,10);
57  TH1F* gjetEtaTrgLow = new TH1F("gjetEtaTrgLow","gjetEtaTrgLow",100,-10,10);
58  TH1F* gjetEta = new TH1F("gjetEta","gjetEta",100,-10,10);
59  TH1F* gjetPhiTrg = new TH1F("gjetPhiTrg","gjetPhiTrg",100,-4,4);
60  TH1F* gjetPhiTrg2 = new TH1F("gjetPhiTrg2","gjetPhiTrg2",100,-4,4);
61  TH1F* gjetPhiTrgLow = new TH1F("gjetPhiTrgLow","gjetPhiTrgLow",100,-4,4);
62  TH1F* gjetPhi = new TH1F("gjetPhi","gjetPhi",100,-4,4);
63  TH1F* ghtTrg = new TH1F("ghtTrg","ghtTrg",100,0,1000);
64  TH1F* ghtTrg2 = new TH1F("ghtTrg2","ghtTrg2",100,0,1000);
65  TH1F* ghtTrgLow = new TH1F("ghtTrgLow","ghtTrgLow",100,0,1000);
66  TH1F* ght = new TH1F("ght","ght",100,0,1000);
67  TH1F* rmetTrg = new TH1F("rmetTrg","rmetTrg",100,0,500);
68  TH1F* rmetTrg2 = new TH1F("rmetTrg2","rmetTrg2",100,0,500);
69  TH1F* rmetTrgLow = new TH1F("rmetTrgLow","rmetTrgLow",100,0,500);
70  TH1F* rmet = new TH1F("rmet","rmet",100,0,500);
71  TH1F* rjetTrg = new TH1F("rjetTrg","rjetTrg",100,0,500);
72  TH1F* rjetTrg2 = new TH1F("rjetTrg2","rjetTrg2",100,0,500);
73  TH1F* rjetTrgLow = new TH1F("rjetTrgLow","rjetTrgLow",100,0,500);
74  TH1F* rjet = new TH1F("rjet","rjet",100,0,500);
75  TH1F* rjetEtaTrg = new TH1F("rjetEtaTrg","rjetEtaTrg",100,-10,10);
76  TH1F* rjetEtaTrg2 = new TH1F("rjetEtaTrg2","rjetEtaTrg2",100,-10,10);
77  TH1F* rjetEtaTrgLow = new TH1F("rjetEtaTrgLow","rjetEtaTrgLow",100,-10,10);
78  TH1F* rjetEta = new TH1F("rjetEta","rjetEta",100,-10,10);
79  TH1F* rjetPhiTrg = new TH1F("rjetPhiTrg","rjetPhiTrg",100,-4,4);
80  TH1F* rjetPhiTrg2 = new TH1F("rjetPhiTrg2","rjetPhiTrg2",100,-4,4);
81  TH1F* rjetPhiTrgLow = new TH1F("rjetPhiTrgLow","rjetPhiTrgLow",100,-4,4);
82  TH1F* rjetPhi = new TH1F("rjetPhi","rjetPhi",100,-4,4);
83  TH1F* rhtTrg = new TH1F("rhtTrg","rhtTrg",100,0,1000);
84  TH1F* rhtTrg2 = new TH1F("rhtTrg2","rhtTrg2",100,0,1000);
85  TH1F* rhtTrgLow = new TH1F("rhtTrgLow","rhtTrgLow",100,0,1000);
86  TH1F* rht = new TH1F("rht","rht",100,0,1000);
87 
88  TProfile* gmto = new TProfile("_meTurnOngMET","Gen Missing ET Turn-On RelVal",100,0,500,0,1);
89  TProfile* gmtol = new TProfile("_meTurnOngMETLow","Gen Missing ET Turn-On Data",100,0,500,0,1);
90  TProfile* gjto = new TProfile("_meTurnOngJetPt","Gen Jet Pt Turn-On RelVal",100,0,500,0,1);
91  TProfile* gjtol = new TProfile("_meTurnOngJetPtLow","Gen Jet Pt Turn-On Data",100,0,500,0,1);
92  TProfile* gjeto = new TProfile("_meTurnOngJetEta","Gen Jet Eta Turn-On RelVal",100,-10,10,0,1);
93  TProfile* gjetol = new TProfile("_meTurnOngJetEtaLow","Gen Jet Eta Turn-On Data",100,-10,10,0,1);
94  TProfile* gjpto = new TProfile("_meTurnOngJetPhi","Gen Jet Phi Turn-On RelVal",100,-4,4,0,1);
95  TProfile* gjptol = new TProfile("_meTurnOngJetPhiLow","Gen Jet Phi Turn-On Data",100,-4,4,0,1);
96  TProfile* ghto = new TProfile("_meTurnOngHT","Gen HT Turn-On RelVal",100,0,1000,0,1);
97  TProfile* ghtol = new TProfile("_meTurnOngHTLow","Gen HT Turn-On Data",100,0,1000,0,1);
98  TProfile* rmto = new TProfile("_meTurnOnrMET","Reco Missing ET Turn-On RelVal",100,0,500,0,1);
99  TProfile* rmtol = new TProfile("_meTurnOnrMETLow","Reco Missing ET Turn-On Data",100,0,500,0,1);
100  TProfile* rjto = new TProfile("_meTurnOnrJetPt","Reco Jet Pt Turn-On RelVal",100,0,500,0,1);
101  TProfile* rjtol = new TProfile("_meTurnOnrJetPtLow","Reco Jet Pt Turn-On Data",100,0,500,0,1);
102  TProfile* rjeto = new TProfile("_meTurnOnrJetEta","Reco Jet Eta Turn-On RelVal",100,-10,10,0,1);
103  TProfile* rjetol = new TProfile("_meTurnOnrJetEtaLow","Reco Jet Eta Turn-On Data",100,-10,10,0,1);
104  TProfile* rjpto = new TProfile("_meTurnOnrJetPhi","Reco Jet Phi Turn-On RelVal",100,-4,4,0,1);
105  TProfile* rjptol = new TProfile("_meTurnOnrJetPhiLow","Reco Jet Phi Turn-On Data",100,-4,4,0,1);
106  TProfile* rhto = new TProfile("_meTurnOnrHT","Reco HT Turn-On RelVal",100,0,1000,0,1);
107  TProfile* rhtol = new TProfile("_meTurnOnrHTLow","Reco HT Turn-On Data",100,0,1000,0,1);
108 
109  //std::vector<std::string> mes = dqm->getMEs();
110  //for(std::vector<std::string>::iterator me = mes.begin() ;me!= mes.end(); me++ )
111  // std::cout <<*me <<std::endl;
112  //std::cout <<std::endl;
113 
114  gmetTrg->Add(dqm->get(dqm->pwd() + "/_meGenMETTrg")->getTH1F(),1);
115  gmetTrg->Sumw2();
116  gmetTrg2->Add(dqm->get(dqm->pwd() + "/_meGenMETTrg")->getTH1F(),1);
117  gmetTrg2->Sumw2();
118  gmetTrgLow->Add(dqm->get(dqm->pwd() + "/_meGenMETTrgLow")->getTH1F(),1);
119  gmetTrgLow->Sumw2();
120  gmet->Add(dqm->get(dqm->pwd() + "/_meGenMET")->getTH1F(),1);
121  gmet->Sumw2();
122  gmetTrg->Divide(gmetTrg,gmet,1,1,"B");
123  gmetTrg2->Divide(gmetTrg2,gmetTrgLow,1,1,"B");
124 
125  rmetTrg->Add(dqm->get(dqm->pwd() + "/_meRecoMETTrg")->getTH1F(),1);
126  rmetTrg->Sumw2();
127  rmetTrg2->Add(dqm->get(dqm->pwd() + "/_meRecoMETTrg")->getTH1F(),1);
128  rmetTrg2->Sumw2();
129  rmetTrgLow->Add(dqm->get(dqm->pwd() + "/_meRecoMETTrgLow")->getTH1F(),1);
130  rmetTrgLow->Sumw2();
131  rmet->Add(dqm->get(dqm->pwd() + "/_meRecoMET")->getTH1F(),1);
132  rmet->Sumw2();
133  rmetTrg->Divide(rmetTrg,rmet,1,1,"B");
134  rmetTrg2->Divide(rmetTrg2,rmetTrgLow,1,1,"B");
135 
136  gjetTrg->Add(dqm->get(dqm->pwd() + "/_meGenJetPtTrg")->getTH1F(),1);
137  gjetTrg->Sumw2();
138  gjetTrg2->Add(dqm->get(dqm->pwd() + "/_meGenJetPtTrg")->getTH1F(),1);
139  gjetTrg2->Sumw2();
140  gjetTrgLow->Add(dqm->get(dqm->pwd() + "/_meGenJetPtTrgLow")->getTH1F(),1);
141  gjetTrgLow->Sumw2();
142  gjet->Add(dqm->get(dqm->pwd() + "/_meGenJetPt")->getTH1F(),1);
143  gjet->Sumw2();
144  gjetTrg->Divide(gjetTrg,gjet,1,1,"B");
145  gjetTrg2->Divide(gjetTrg2,gjetTrgLow,1,1,"B");
146 
147  rjetTrg->Add(dqm->get(dqm->pwd() + "/_meRecoJetPtTrg")->getTH1F(),1);
148  rjetTrg->Sumw2();
149  rjetTrg2->Add(dqm->get(dqm->pwd() + "/_meRecoJetPtTrg")->getTH1F(),1);
150  rjetTrg2->Sumw2();
151  rjetTrgLow->Add(dqm->get(dqm->pwd() + "/_meRecoJetPtTrgLow")->getTH1F(),1);
152  rjetTrgLow->Sumw2();
153  rjet->Add(dqm->get(dqm->pwd() + "/_meRecoJetPt")->getTH1F(),1);
154  rjet->Sumw2();
155  rjetTrg->Divide(rjetTrg,rjet,1,1,"B");
156  rjetTrg2->Divide(rjetTrg2,rjetTrgLow,1,1,"B");
157 
158  gjetEtaTrg->Add(dqm->get(dqm->pwd() + "/_meGenJetEtaTrg")->getTH1F(),1);
159  gjetEtaTrg->Sumw2();
160  gjetEtaTrg2->Add(dqm->get(dqm->pwd() + "/_meGenJetEtaTrg")->getTH1F(),1);
161  gjetEtaTrg2->Sumw2();
162  gjetEtaTrgLow->Add(dqm->get(dqm->pwd() + "/_meGenJetEtaTrgLow")->getTH1F(),1);
163  gjetEtaTrgLow->Sumw2();
164  gjetEta->Add(dqm->get(dqm->pwd() + "/_meGenJetEta")->getTH1F(),1);
165  gjetEta->Sumw2();
166  gjetEtaTrg->Divide(gjetEtaTrg,gjetEta,1,1,"B");
167  gjetEtaTrg2->Divide(gjetEtaTrg2,gjetEtaTrgLow,1,1,"B");
168 
169  rjetEtaTrg->Add(dqm->get(dqm->pwd() + "/_meRecoJetEtaTrg")->getTH1F(),1);
170  rjetEtaTrg->Sumw2();
171  rjetEtaTrg2->Add(dqm->get(dqm->pwd() + "/_meRecoJetEtaTrg")->getTH1F(),1);
172  rjetEtaTrg2->Sumw2();
173  rjetEtaTrgLow->Add(dqm->get(dqm->pwd() + "/_meRecoJetEtaTrgLow")->getTH1F(),1);
174  rjetEtaTrgLow->Sumw2();
175  rjetEta->Add(dqm->get(dqm->pwd() + "/_meRecoJetEta")->getTH1F(),1);
176  rjetEta->Sumw2();
177  rjetEtaTrg->Divide(rjetEtaTrg,rjetEta,1,1,"B");
178  rjetEtaTrg2->Divide(rjetEtaTrg2,rjetEtaTrgLow,1,1,"B");
179 
180  gjetPhiTrg->Add(dqm->get(dqm->pwd() + "/_meGenJetPhiTrg")->getTH1F(),1);
181  gjetPhiTrg->Sumw2();
182  gjetPhiTrg2->Add(dqm->get(dqm->pwd() + "/_meGenJetPhiTrg")->getTH1F(),1);
183  gjetPhiTrg2->Sumw2();
184  gjetPhiTrgLow->Add(dqm->get(dqm->pwd() + "/_meGenJetPhiTrgLow")->getTH1F(),1);
185  gjetPhiTrgLow->Sumw2();
186  gjetPhi->Add(dqm->get(dqm->pwd() + "/_meGenJetPhi")->getTH1F(),1);
187  gjetPhi->Sumw2();
188  gjetPhiTrg->Divide(gjetPhiTrg,gjetPhi,1,1,"B");
189  gjetPhiTrg2->Divide(gjetPhiTrg2,gjetPhiTrgLow,1,1,"B");
190 
191  rjetPhiTrg->Add(dqm->get(dqm->pwd() + "/_meRecoJetPhiTrg")->getTH1F(),1);
192  rjetPhiTrg->Sumw2();
193  rjetPhiTrg2->Add(dqm->get(dqm->pwd() + "/_meRecoJetPhiTrg")->getTH1F(),1);
194  rjetPhiTrg2->Sumw2();
195  rjetPhiTrgLow->Add(dqm->get(dqm->pwd() + "/_meRecoJetPhiTrgLow")->getTH1F(),1);
196  rjetPhiTrgLow->Sumw2();
197  rjetPhi->Add(dqm->get(dqm->pwd() + "/_meRecoJetPhi")->getTH1F(),1);
198  rjetPhi->Sumw2();
199  rjetPhiTrg->Divide(rjetPhiTrg,rjetPhi,1,1,"B");
200  rjetPhiTrg2->Divide(rjetPhiTrg2,rjetPhiTrgLow,1,1,"B");
201 
202  ghtTrg->Add(dqm->get(dqm->pwd() + "/_meGenHTTrg")->getTH1F(),1);
203  ghtTrg->Sumw2();
204  ghtTrg2->Add(dqm->get(dqm->pwd() + "/_meGenHTTrg")->getTH1F(),1);
205  ghtTrg2->Sumw2();
206  ghtTrgLow->Add(dqm->get(dqm->pwd() + "/_meGenHTTrgLow")->getTH1F(),1);
207  ghtTrgLow->Sumw2();
208  ght->Add(dqm->get(dqm->pwd() + "/_meGenHT")->getTH1F(),1);
209  ght->Sumw2();
210  ghtTrg->Divide(ghtTrg,ght,1,1,"B");
211  ghtTrg2->Divide(ghtTrg2,ghtTrgLow,1,1,"B");
212 
213  rhtTrg->Add(dqm->get(dqm->pwd() + "/_meRecoHTTrg")->getTH1F(),1);
214  rhtTrg->Sumw2();
215  rhtTrg2->Add(dqm->get(dqm->pwd() + "/_meRecoHTTrg")->getTH1F(),1);
216  rhtTrg2->Sumw2();
217  rhtTrgLow->Add(dqm->get(dqm->pwd() + "/_meRecoHTTrgLow")->getTH1F(),1);
218  rhtTrgLow->Sumw2();
219  rht->Add(dqm->get(dqm->pwd() + "/_meRecoHT")->getTH1F(),1);
220  rht->Sumw2();
221  rhtTrg->Divide(rhtTrg,rht,1,1,"B");
222  rhtTrg2->Divide(rhtTrg2,rhtTrgLow,1,1,"B");
223 
224  double val,err;
225  for (int ib=0;ib<100;ib++) {
226  //genmet relval
227  val = gmetTrg->GetBinContent(ib+1);
228  gmto->SetBinContent(ib+1,val);
229  gmto->SetBinEntries(ib+1,1);
230  err = gmetTrg->GetBinError(ib+1);
231  gmto->SetBinError(ib+1,sqrt(err*err+val*val));
232  //genmet data
233  val = gmetTrg2->GetBinContent(ib+1);
234  gmtol->SetBinContent(ib+1,val);
235  gmtol->SetBinEntries(ib+1,1);
236  err = gmetTrg2->GetBinError(ib+1);
237  gmtol->SetBinError(ib+1,sqrt(err*err+val*val));
238  //recmet relval
239  val = rmetTrg->GetBinContent(ib+1);
240  rmto->SetBinContent(ib+1,val);
241  rmto->SetBinEntries(ib+1,1);
242  err = rmetTrg->GetBinError(ib+1);
243  rmto->SetBinError(ib+1,sqrt(err*err+val*val));
244  //recmet data
245  val = rmetTrg2->GetBinContent(ib+1);
246  rmtol->SetBinContent(ib+1,val);
247  rmtol->SetBinEntries(ib+1,1);
248  err = rmetTrg2->GetBinError(ib+1);
249  rmtol->SetBinError(ib+1,sqrt(err*err+val*val));
250  //genjet relval
251  val = gjetTrg->GetBinContent(ib+1);
252  gjto->SetBinContent(ib+1,val);
253  gjto->SetBinEntries(ib+1,1);
254  err = gjetTrg->GetBinError(ib+1);
255  gjto->SetBinError(ib+1,sqrt(err*err+val*val));
256  //genjet data
257  val = gjetTrg2->GetBinContent(ib+1);
258  gjtol->SetBinContent(ib+1,val);
259  gjtol->SetBinEntries(ib+1,1);
260  err = gjetTrg2->GetBinError(ib+1);
261  gjtol->SetBinError(ib+1,sqrt(err*err+val*val));
262  //recjet relval
263  val = rjetTrg->GetBinContent(ib+1);
264  rjto->SetBinContent(ib+1,val);
265  rjto->SetBinEntries(ib+1,1);
266  err = rjetTrg->GetBinError(ib+1);
267  rjto->SetBinError(ib+1,sqrt(err*err+val*val));
268  //recjet data
269  val = rjetTrg2->GetBinContent(ib+1);
270  rjtol->SetBinContent(ib+1,val);
271  rjtol->SetBinEntries(ib+1,1);
272  err = rjetTrg2->GetBinError(ib+1);
273  rjtol->SetBinError(ib+1,sqrt(err*err+val*val));
274  //genjeteta relval
275  val = gjetEtaTrg->GetBinContent(ib+1);
276  gjeto->SetBinContent(ib+1,val);
277  gjeto->SetBinEntries(ib+1,1);
278  err = gjetEtaTrg->GetBinError(ib+1);
279  gjeto->SetBinError(ib+1,sqrt(err*err+val*val));
280  //genjeteta data
281  val = gjetEtaTrg2->GetBinContent(ib+1);
282  gjetol->SetBinContent(ib+1,val);
283  gjetol->SetBinEntries(ib+1,1);
284  err = gjetEtaTrg2->GetBinError(ib+1);
285  gjetol->SetBinError(ib+1,sqrt(err*err+val*val));
286  //recjeteta relval
287  val = rjetEtaTrg->GetBinContent(ib+1);
288  rjeto->SetBinContent(ib+1,val);
289  rjeto->SetBinEntries(ib+1,1);
290  err = rjetEtaTrg->GetBinError(ib+1);
291  rjeto->SetBinError(ib+1,sqrt(err*err+val*val));
292  //recjeteta data
293  val = rjetEtaTrg2->GetBinContent(ib+1);
294  rjetol->SetBinContent(ib+1,val);
295  rjetol->SetBinEntries(ib+1,1);
296  err = rjetEtaTrg2->GetBinError(ib+1);
297  rjetol->SetBinError(ib+1,sqrt(err*err+val*val));
298  //genjetphi relval
299  val = gjetPhiTrg->GetBinContent(ib+1);
300  gjpto->SetBinContent(ib+1,val);
301  gjpto->SetBinEntries(ib+1,1);
302  err = gjetPhiTrg->GetBinError(ib+1);
303  gjpto->SetBinError(ib+1,sqrt(err*err+val*val));
304  //genjetphi data
305  val = gjetPhiTrg2->GetBinContent(ib+1);
306  gjptol->SetBinContent(ib+1,val);
307  gjptol->SetBinEntries(ib+1,1);
308  err = gjetPhiTrg2->GetBinError(ib+1);
309  gjptol->SetBinError(ib+1,sqrt(err*err+val*val));
310  //recjetphi relval
311  val = rjetPhiTrg->GetBinContent(ib+1);
312  rjpto->SetBinContent(ib+1,val);
313  rjpto->SetBinEntries(ib+1,1);
314  err = rjetPhiTrg->GetBinError(ib+1);
315  rjpto->SetBinError(ib+1,sqrt(err*err+val*val));
316  //recjetphi data
317  val = rjetPhiTrg2->GetBinContent(ib+1);
318  rjptol->SetBinContent(ib+1,val);
319  rjptol->SetBinEntries(ib+1,1);
320  err = rjetPhiTrg2->GetBinError(ib+1);
321  rjptol->SetBinError(ib+1,sqrt(err*err+val*val));
322  //genht relval
323  val = ghtTrg->GetBinContent(ib+1);
324  ghto->SetBinContent(ib+1,val);
325  ghto->SetBinEntries(ib+1,1);
326  err = ghtTrg->GetBinError(ib+1);
327  ghto->SetBinError(ib+1,sqrt(err*err+val*val));
328  //genht data
329  val = ghtTrg2->GetBinContent(ib+1);
330  ghtol->SetBinContent(ib+1,val);
331  ghtol->SetBinEntries(ib+1,1);
332  err = ghtTrg2->GetBinError(ib+1);
333  ghtol->SetBinError(ib+1,sqrt(err*err+val*val));
334  //recht relval
335  val = rhtTrg->GetBinContent(ib+1);
336  rhto->SetBinContent(ib+1,val);
337  rhto->SetBinEntries(ib+1,1);
338  err = rhtTrg->GetBinError(ib+1);
339  rhto->SetBinError(ib+1,sqrt(err*err+val*val));
340  //recht data
341  val = rhtTrg2->GetBinContent(ib+1);
342  rhtol->SetBinContent(ib+1,val);
343  rhtol->SetBinEntries(ib+1,1);
344  err = rhtTrg2->GetBinError(ib+1);
345  rhtol->SetBinError(ib+1,sqrt(err*err+val*val));
346 
347 
348  //std::cout <<"MET:"<<_meTurnOnMET->getBinContent(ib+1)<<" "<<gmetTrg->GetBinContent(ib+1)<<" "<<_meTurnOnMET->getBinError(ib+1)<<" "<<gmetTrg->GetBinError(ib+1)<<std::endl;
349  //std::cout <<"JET:"<<_meTurnOnJetPt->getBinContent(ib+1)<<" "<<gjetTrg->GetBinContent(ib+1)<<" "<<_meTurnOnJetPt->getBinError(ib+1)<<" "<<gjetTrg->GetBinError(ib+1)<<std::endl;
350  }
351  dqm->bookProfile("Gen Missing ET Turn-On RelVal",gmto);
352  dqm->bookProfile("Gen Missing ET Turn-On Data",gmtol);
353  dqm->bookProfile("Reco Missing ET Turn-On RelVal",rmto);
354  dqm->bookProfile("Reco Missing ET Turn-On Data",rmtol);
355  dqm->bookProfile("Gen Jet Pt Turn-On RelVal",gjto);
356  dqm->bookProfile("Gen Jet Pt Turn-On Data",gjtol);
357  dqm->bookProfile("Reco Jet Pt Turn-On RelVal",rjto);
358  dqm->bookProfile("Reco Jet Pt Turn-On Data",rjtol);
359  dqm->bookProfile("Gen Jet Eta Turn-On RelVal",gjeto);
360  dqm->bookProfile("Gen Jet Eta Turn-On Data",gjetol);
361  dqm->bookProfile("Reco Jet Eta Turn-On RelVal",rjeto);
362  dqm->bookProfile("Reco Jet Eta Turn-On Data",rjetol);
363  dqm->bookProfile("Gen Jet Phi Turn-On RelVal",gjpto);
364  dqm->bookProfile("Gen Jet Phi Turn-On Data",gjptol);
365  dqm->bookProfile("Reco Jet Phi Turn-On RelVal",rjpto);
366  dqm->bookProfile("Reco Jet Phi Turn-On Data",rjptol);
367  dqm->bookProfile("Gen HT Turn-On RelVal",ghto);
368  dqm->bookProfile("Gen HT Turn-On Data",ghtol);
369  dqm->bookProfile("Reco HT Turn-On RelVal",rhto);
370  dqm->bookProfile("Reco HT Turn-On Data",rhtol);
371  delete gjto;
372  delete gjtol;
373  delete gjeto;
374  delete gjetol;
375  delete gjpto;
376  delete gjptol;
377  delete rjto;
378  delete rjtol;
379  delete rjeto;
380  delete rjetol;
381  delete rjpto;
382  delete rjptol;
383  delete gmto;
384  delete gmtol;
385  delete rmto;
386  delete rmtol;
387  delete ghto;
388  delete ghtol;
389  delete rhto;
390  delete rhtol;
391  delete gmetTrg;
392  delete gmetTrg2;
393  delete gmetTrgLow;
394  delete gmet;
395  delete rmetTrg;
396  delete rmetTrg2;
397  delete rmetTrgLow;
398  delete rmet;
399  delete gjetTrg;
400  delete gjetTrg2;
401  delete gjetTrgLow;
402  delete gjet;
403  delete rjetTrg;
404  delete rjetTrg2;
405  delete rjetTrgLow;
406  delete rjet;
407  delete gjetEtaTrg;
408  delete gjetEtaTrg2;
409  delete gjetEtaTrgLow;
410  delete gjetEta;
411  delete rjetEtaTrg;
412  delete rjetEtaTrg2;
413  delete rjetEtaTrgLow;
414  delete rjetEta;
415  delete gjetPhiTrg;
416  delete gjetPhiTrg2;
417  delete gjetPhiTrgLow;
418  delete gjetPhi;
419  delete rjetPhiTrg;
420  delete rjetPhiTrg2;
421  delete rjetPhiTrgLow;
422  delete rjetPhi;
423  delete ghtTrg;
424  delete ghtTrg2;
425  delete ghtTrgLow;
426  delete ght;
427  delete rhtTrg;
428  delete rhtTrg2;
429  delete rhtTrgLow;
430  delete rht;
431 
432  dqm->goUp();
433  }
434 }
ght
double ght
Definition: hdecay.h:106
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
getTH1F
TH1F * getTH1F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
Definition: HcalHistoUtils.h:57
dir
dbl *** dir
Definition: mlp_gen.cc:35

Member Data Documentation

std::string JetMETDQMPostProcessor::subDir_
private

Definition at line 18 of file JetMETDQMPostProcessor.h.

Referenced by endRun(), and JetMETDQMPostProcessor().

Generated for CMSSW Reference Manual by doxygen 1.8.5