CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
METTester.cc
Go to the documentation of this file.
1 // author: Mike Schmitt, University of Florida
2 // first version 8/24/2006
3 // modification: Bobby Scurlock
4 // date: 03.11.2006
5 // note: added RMS(METx) vs SumET capability
6 // modification: Rick Cavanaugh
7 // date: 05.11.2006
8 // note: cleaned up constructor and beginJob, removed int conv. warning
9 // added configuration params
10 // modification: Mike Schmitt
11 // date: 02.28.2007
12 // note: code rewrite. Now uses STL map for monitoring element container.
13 // modification: Bobby Scurlock
14 // date: 04.03.2007
15 // note: Eliminated automated resolution fitting. This is now done in a ROOT script.
16 
17 // date: 02.04.2009
18 // note: Added option to use fine binning or course binning for histos
19 //
20 // modification: Samantha Hewamanage, Florida International University
21 // date: 01.30.2012
22 // note: Added few hists for various nvtx ranges to study PU effects.
23 // Cleaned up the code by making it readable and const'ing the
24 // variables that should be changed.
25 // Changed the number of bins from odd to even. Odd number of bins
26 // makes it impossible to rebin a hist.
27 #include "METTester.h"
28 using namespace reco;
29 using namespace std;
30 using namespace edm;
31 
32 METTester::METTester(const edm::ParameterSet& iConfig)
33 {
34 
35  inputMETLabel_ =iConfig.getParameter<edm::InputTag>("InputMETLabel");
36  METType_ =iConfig.getUntrackedParameter<std::string>("METType");
37  mOutputFile =iConfig.getUntrackedParameter<std::string>("OutputFile","");
38 
39  std::string inputMETCollectionLabel(inputMETLabel_.label());
40 
41  isCaloMET = (std::string("calo")==METType_);
42  isPFMET = (std::string("pf") ==METType_);
43  isGenMET = (std::string("gen") ==METType_);
44 
45  pvToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag("offlinePrimaryVertices"));
46  if (isCaloMET) caloMETsToken_ = consumes<reco::CaloMETCollection> (inputMETLabel_);
47  if (isPFMET) pfMETsToken_ = consumes<reco::PFMETCollection> (inputMETLabel_);
48  if (isGenMET) genMETsToken_ = consumes<reco::GenMETCollection> (inputMETLabel_);
49  genMETsTrueToken_ = consumes<reco::GenMETCollection> (edm::InputTag("genMetTrue"));
50  genMETsCaloToken_ = consumes<reco::GenMETCollection> (edm::InputTag("genMetCalo"));
51  //Events variables
52  mNvertex = 0;
53 
54  //Common variables
55  mMEx = 0;
56  mMEy = 0;
57  mMETSig = 0;
58  mMET = 0;
59  mMETFine = 0;
60  mMET_Nvtx = 0;
61  mMETPhi = 0;
62  mSumET = 0;
63  mMETDifference_GenMETTrue = 0;
64  mMETDeltaPhi_GenMETTrue = 0;
65  mMETDifference_GenMETCalo = 0;
66  mMETDeltaPhi_GenMETCalo = 0;
67 
68  //CaloMET variables
69  mCaloMaxEtInEmTowers = 0;
70  mCaloMaxEtInHadTowers = 0;
71  mCaloEtFractionHadronic = 0;
72  mCaloEmEtFraction = 0;
73  mCaloHadEtInHB = 0;
74  mCaloHadEtInHO = 0;
75  mCaloHadEtInHE = 0;
76  mCaloHadEtInHF = 0;
77  mCaloEmEtInHF = 0;
78  mCaloSETInpHF = 0;
79  mCaloSETInmHF = 0;
80  mCaloEmEtInEE = 0;
81  mCaloEmEtInEB = 0;
82 
83 
84  //GenMET variables
85  mNeutralEMEtFraction=0;
86  mNeutralHadEtFraction=0;
87  mChargedEMEtFraction=0;
88  mChargedHadEtFraction=0;
89  mMuonEtFraction=0;
90  mInvisibleEtFraction=0;
91 
92  //MET variables
93 
94  //PFMET variables
95  mMETDifference_GenMETTrue_MET0to20=0;
96  mMETDifference_GenMETTrue_MET20to40=0;
97  mMETDifference_GenMETTrue_MET40to60=0;
98  mMETDifference_GenMETTrue_MET60to80=0;
99  mMETDifference_GenMETTrue_MET80to100=0;
100  mMETDifference_GenMETTrue_MET100to150=0;
101  mMETDifference_GenMETTrue_MET150to200=0;
102  mMETDifference_GenMETTrue_MET200to300=0;
103  mMETDifference_GenMETTrue_MET300to400=0;
104  mMETDifference_GenMETTrue_MET400to500=0;
105  mMETDifference_GenMETTrue_METResolution=0;
106 
107  //TCMET specific variables
108  mMExCorrection=0;
109  mMEyCorrection=0;
110  mMuonCorrectionFlag=0;
111  mtrkPt=0;
112  mtrkEta=0;
113  mtrkNhits=0;
114  mtrkChi2=0;
115  mtrkD0=0;
116  mtrkQuality=0;
117  mtrkAlgo=0;
118  mtrkPtErr=0;
119  melePt=0;
120  meleEta=0;
121  meleHoE=0;
122 
123  mmuPt=0;
124  mmuEta=0;
125  mmuNhits=0;
126  mmuChi2=0;
127  mmuD0=0;
128  mnMus=0;
129  mnMusPis=0;
130  mnEls=0;
131  mfracTrks=0;
132  mdMET=0;
133  mdMETx=0;
134  mdMETy=0;
135  mdMEy=0;
136  mdMUx=0;
137  mdMUy=0;
138 
139 }
140 void METTester::bookHistograms(DQMStore::IBooker & ibooker,
141  edm::Run const & iRun,
142  edm::EventSetup const & /* iSetup */)
143 {
144  ibooker.setCurrentFolder("JetMET/METValidation/"+inputMETLabel_.label()) ;
145 
146  mNvertex = ibooker.book1D("Nvertex","Nvertex",80,0,80);
147  mMEx = ibooker.book1D("MEx","MEx",160,-800,800);
148  mMEy = ibooker.book1D("MEy","MEy",160,-800,800);
149  mMETSig = ibooker.book1D("METSig","METSig",25,0,24.5);
150  mMET = ibooker.book1D("MET", "MET (20 GeV binning)" , 100,0,2000);
151  mMETFine = ibooker.book1D("METFine", "MET (2 GeV binning)" , 1000,0,2000);
152  mMET_Nvtx = ibooker.bookProfile("MET_Nvtx", "MET vs. nvtx", 60, 0., 60., 0., 2000., " ");
153  mMETPhi = ibooker.book1D("METPhi","METPhi",80,-4,4);
154  mSumET = ibooker.book1D("SumET" , "SumET" , 200,0,4000); //10GeV
155  mMETDifference_GenMETTrue = ibooker.book1D("METDifference_GenMETTrue","METDifference_GenMETTrue", 500,-500,500);
156  mMETDeltaPhi_GenMETTrue = ibooker.book1D("METDeltaPhi_GenMETTrue","METDeltaPhi_GenMETTrue", 80,0,4);
157  mMETDifference_GenMETCalo = ibooker.book1D("METDifference_GenMETCalo","METDifference_GenMETCalo", 500,-500,500);
158  mMETDeltaPhi_GenMETCalo = ibooker.book1D("METDeltaPhi_GenMETCalo","METDeltaPhi_GenMETCalo", 80,0,4);
159 
160  mMETDifference_GenMETTrue_MET0to20 = ibooker.book1D("METResolution_GenMETTrue_MET0to20" , "METResolution_GenMETTrue_MET0to20" , 500,-500,500);
161  mMETDifference_GenMETTrue_MET20to40 = ibooker.book1D("METResolution_GenMETTrue_MET20to40" , "METResolution_GenMETTrue_MET20to40" , 500,-500,500);
162  mMETDifference_GenMETTrue_MET40to60 = ibooker.book1D("METResolution_GenMETTrue_MET40to60" , "METResolution_GenMETTrue_MET40to60" , 500,-500,500);
163  mMETDifference_GenMETTrue_MET60to80 = ibooker.book1D("METResolution_GenMETTrue_MET60to80" , "METResolution_GenMETTrue_MET60to80" , 500,-500,500);
164  mMETDifference_GenMETTrue_MET80to100 = ibooker.book1D("METResolution_GenMETTrue_MET80to100" , "METResolution_GenMETTrue_MET80to100" , 500,-500,500);
165  mMETDifference_GenMETTrue_MET100to150 = ibooker.book1D("METResolution_GenMETTrue_MET100to150", "METResolution_GenMETTrue_MET100to150", 500,-500,500);
166  mMETDifference_GenMETTrue_MET150to200 = ibooker.book1D("METResolution_GenMETTrue_MET150to200", "METResolution_GenMETTrue_MET150to200", 500,-500,500);
167  mMETDifference_GenMETTrue_MET200to300 = ibooker.book1D("METResolution_GenMETTrue_MET200to300", "METResolution_GenMETTrue_MET200to300", 500,-500,500);
168  mMETDifference_GenMETTrue_MET300to400 = ibooker.book1D("METResolution_GenMETTrue_MET300to400", "METResolution_GenMETTrue_MET300to400", 500,-500,500);
169  mMETDifference_GenMETTrue_MET400to500 = ibooker.book1D("METResolution_GenMETTrue_MET400to500", "METResolution_GenMETTrue_MET400to500", 500,-500,500);
170  //this will be filled at the end of the job using info from above hists
171  int nBins = 10;
172  float bins[] = {0.,20.,40.,60.,80.,100.,150.,200.,300.,400.,500.};
173  mMETDifference_GenMETTrue_METResolution = ibooker.book1D("METResolution_GenMETTrue_InMETBins","METResolution_GenMETTrue_InMETBins",nBins, bins);
174 
175  if ( isCaloMET) {
176  mCaloMaxEtInEmTowers = ibooker.book1D("CaloMaxEtInEmTowers","CaloMaxEtInEmTowers",300,0,1500); //5GeV
177  mCaloMaxEtInHadTowers = ibooker.book1D("CaloMaxEtInHadTowers","CaloMaxEtInHadTowers",300,0,1500); //5GeV
178  mCaloEtFractionHadronic = ibooker.book1D("CaloEtFractionHadronic","CaloEtFractionHadronic",100,0,1);
179  mCaloEmEtFraction = ibooker.book1D("CaloEmEtFraction","CaloEmEtFraction",100,0,1);
180  mCaloHadEtInHB = ibooker.book1D("CaloHadEtInHB","CaloHadEtInHB", 200, 0, 2000); //5GeV
181  mCaloHadEtInHE = ibooker.book1D("CaloHadEtInHE","CaloHadEtInHE", 100, 0, 500); //5GeV
182  mCaloHadEtInHO = ibooker.book1D("CaloHadEtInHO","CaloHadEtInHO", 100, 0, 200); //5GeV
183  mCaloHadEtInHF = ibooker.book1D("CaloHadEtInHF","CaloHadEtInHF", 100, 0, 200); //5GeV
184  mCaloSETInpHF = ibooker.book1D("CaloSETInpHF","CaloSETInpHF",100, 0, 500);
185  mCaloSETInmHF = ibooker.book1D("CaloSETInmHF","CaloSETInmHF",100, 0, 500);
186  mCaloEmEtInEE = ibooker.book1D("CaloEmEtInEE","CaloEmEtInEE",100, 0, 500); //5GeV
187  mCaloEmEtInEB = ibooker.book1D("CaloEmEtInEB","CaloEmEtInEB",100, 0, 500); //5GeV
188  mCaloEmEtInHF = ibooker.book1D("CaloEmEtInHF","CaloEmEtInHF",100, 0, 500); //5GeV
189  }
190 
191  if(isGenMET){
192  mNeutralEMEtFraction = ibooker.book1D("GenNeutralEMEtFraction", "GenNeutralEMEtFraction", 120, 0.0, 1.2 );
193  mNeutralHadEtFraction = ibooker.book1D("GenNeutralHadEtFraction", "GenNeutralHadEtFraction", 120, 0.0, 1.2 );
194  mChargedEMEtFraction = ibooker.book1D("GenChargedEMEtFraction", "GenChargedEMEtFraction", 120, 0.0, 1.2);
195  mChargedHadEtFraction = ibooker.book1D("GenChargedHadEtFraction", "GenChargedHadEtFraction", 120, 0.0,1.2);
196  mMuonEtFraction = ibooker.book1D("GenMuonEtFraction", "GenMuonEtFraction", 120, 0.0, 1.2 );
197  mInvisibleEtFraction = ibooker.book1D("GenInvisibleEtFraction", "GenInvisibleEtFraction", 120, 0.0, 1.2 );
198  }
199 
200  if (isPFMET){
201  mPFphotonEtFraction = ibooker.book1D("photonEtFraction", "photonEtFraction", 100, 0, 1);
202  mPFphotonEt = ibooker.book1D("photonEt", "photonEt", 100, 0, 1000);
203  mPFneutralHadronEtFraction = ibooker.book1D("neutralHadronEtFraction", "neutralHadronEtFraction", 100, 0, 1);
204  mPFneutralHadronEt = ibooker.book1D("neutralHadronEt", "neutralHadronEt", 100, 0, 1000);
205  mPFelectronEtFraction = ibooker.book1D("electronEtFraction", "electronEtFraction", 100, 0, 1);
206  mPFelectronEt = ibooker.book1D("electronEt", "electronEt", 100, 0, 1000);
207  mPFchargedHadronEtFraction = ibooker.book1D("chargedHadronEtFraction", "chargedHadronEtFraction", 100, 0, 1);
208  mPFchargedHadronEt = ibooker.book1D("chargedHadronEt", "chargedHadronEt", 100, 0, 1000);
209  mPFmuonEtFraction = ibooker.book1D("muonEtFraction", "muonEtFraction", 100, 0, 1);
210  mPFmuonEt = ibooker.book1D("muonEt", "muonEt", 100, 0, 1000);
211  mPFHFHadronEtFraction = ibooker.book1D("HFHadronEtFraction", "HFHadronEtFraction", 100, 0, 1);
212  mPFHFHadronEt = ibooker.book1D("HFHadronEt", "HFHadronEt", 100, 0, 500);
213  mPFHFEMEtFraction = ibooker.book1D("HFEMEtFraction", "HFEMEtFraction", 100, 0, 1);
214  mPFHFEMEt = ibooker.book1D("HFEMEt", "HFEMEt", 100, 0, 300);
215 
216  }
217  }
218 
219 
220 void METTester::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
221 { //int counter(0);
222  edm::Handle<reco::VertexCollection> pvHandle;
223  iEvent.getByToken(pvToken_, pvHandle);
224  if (! pvHandle.isValid())
225  {
226  std::cout << __FUNCTION__ << ":" << __LINE__ << ":pvHandle handle not found!" << std::endl;
227  assert(false);
228  }
229  const int nvtx = pvHandle->size();
230  mNvertex->Fill(nvtx);
231  //Collections for all MET collections
232 
233  edm::Handle<CaloMETCollection> caloMETs;
234  edm::Handle<PFMETCollection> pfMETs;
235  edm::Handle<GenMETCollection> genMETs;
236 
237  if (isCaloMET) iEvent.getByToken(caloMETsToken_, caloMETs);
238  if (isPFMET) iEvent.getByToken(pfMETsToken_, pfMETs);
239  if (isGenMET) iEvent.getByToken(genMETsToken_, genMETs);
240  if ((isCaloMET) and !caloMETs.isValid()) return;
241  if ((isPFMET) and !pfMETs.isValid()) return;
242  if ((isGenMET) and !genMETs.isValid()) return;
243 
244  reco::MET met;
245  if (isCaloMET) { met = caloMETs->front();}
246  if (isPFMET) { met = pfMETs->front() ;}
247  if (isGenMET) { met = genMETs->front() ;}
248 
249  const double SumET = met.sumEt();
250  const double METSig = met.mEtSig();
251  const double MET = met.pt();
252  const double MEx = met.px();
253  const double MEy = met.py();
254  const double METPhi = met.phi();
255  mMEx->Fill(MEx);
256  mMEy->Fill(MEy);
257  mMET->Fill(MET);
258  mMETFine->Fill(MET);
259  mMET_Nvtx->Fill((double)nvtx, MET);
260  mMETPhi->Fill(METPhi);
261  mSumET->Fill(SumET);
262  mMETSig->Fill(METSig);
263 
264  // Get Generated MET for Resolution plots
265  edm::Handle<GenMETCollection> genTrue;
266  iEvent.getByToken(genMETsTrueToken_, genTrue);
267  if (genTrue.isValid()) {
268  const GenMETCollection *genmetcol = genTrue.product();
269  const GenMET *genMetTrue = &(genmetcol->front());
270  double genMET = genMetTrue->pt();
271  double genMETPhi = genMetTrue->phi();
272 
273  mMETDifference_GenMETTrue->Fill( MET - genMET );
274  mMETDeltaPhi_GenMETTrue->Fill( TMath::ACos( TMath::Cos( METPhi - genMETPhi ) ) );
275 
276  //pfMET resolution in pfMET bins : Sam, Feb, 2012
277  if (MET > 0 && MET < 20) mMETDifference_GenMETTrue_MET0to20->Fill( MET - genMET );
278  else if (MET > 20 && MET < 40) mMETDifference_GenMETTrue_MET20to40->Fill( MET - genMET );
279  else if (MET > 40 && MET < 60) mMETDifference_GenMETTrue_MET40to60->Fill( MET - genMET );
280  else if (MET > 60 && MET < 80) mMETDifference_GenMETTrue_MET60to80->Fill( MET - genMET );
281  else if (MET > 80 && MET <100) mMETDifference_GenMETTrue_MET80to100->Fill( MET - genMET );
282  else if (MET >100 && MET <150) mMETDifference_GenMETTrue_MET100to150->Fill( MET - genMET );
283  else if (MET >150 && MET <200) mMETDifference_GenMETTrue_MET150to200->Fill( MET - genMET );
284  else if (MET >200 && MET <300) mMETDifference_GenMETTrue_MET200to300->Fill( MET - genMET );
285  else if (MET >300 && MET <400) mMETDifference_GenMETTrue_MET300to400->Fill( MET - genMET );
286  else if (MET >400 && MET <500) mMETDifference_GenMETTrue_MET400to500->Fill( MET - genMET );
287 
288  } else {
289  edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task: genMetTrue";
290  }
291  edm::Handle<GenMETCollection> genCalo;
292  iEvent.getByToken(genMETsCaloToken_, genCalo);
293  if (genCalo.isValid()) {
294  const GenMETCollection *genmetcol = genCalo.product();
295  const GenMET *genMetCalo = &(genmetcol->front());
296  const double genMET = genMetCalo->pt();
297  const double genMETPhi = genMetCalo->phi();
298 
299  mMETDifference_GenMETCalo->Fill( MET - genMET );
300  mMETDeltaPhi_GenMETCalo->Fill( TMath::ACos( TMath::Cos( METPhi - genMETPhi ) ) );
301  } else {
302  edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task: genMetCalo";
303  }
304 
305  if ( isCaloMET) {
306  const reco::CaloMET * calomet = &(caloMETs->front());
307  // ==========================================================
308  // Reconstructed MET Information
309  const double caloMaxEtInEMTowers = calomet->maxEtInEmTowers();
310  const double caloMaxEtInHadTowers = calomet->maxEtInHadTowers();
311  const double caloEtFractionHadronic = calomet->etFractionHadronic();
312  const double caloEmEtFraction = calomet->emEtFraction();
313  const double caloHadEtInHB = calomet->hadEtInHB();
314  const double caloHadEtInHO = calomet->hadEtInHO();
315  const double caloHadEtInHE = calomet->hadEtInHE();
316  const double caloHadEtInHF = calomet->hadEtInHF();
317  const double caloEmEtInEB = calomet->emEtInEB();
318  const double caloEmEtInEE = calomet->emEtInEE();
319  const double caloEmEtInHF = calomet->emEtInHF();
320  const double caloSETInpHF = calomet->CaloSETInpHF();
321  const double caloSETInmHF = calomet->CaloSETInmHF();
322 
323  mCaloMaxEtInEmTowers->Fill(caloMaxEtInEMTowers);
324  mCaloMaxEtInHadTowers->Fill(caloMaxEtInHadTowers);
325  mCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
326  mCaloEmEtFraction->Fill(caloEmEtFraction);
327  mCaloHadEtInHB->Fill(caloHadEtInHB);
328  mCaloHadEtInHO->Fill(caloHadEtInHO);
329  mCaloHadEtInHE->Fill(caloHadEtInHE);
330  mCaloHadEtInHF->Fill(caloHadEtInHF);
331  mCaloEmEtInEB->Fill(caloEmEtInEB);
332  mCaloEmEtInEE->Fill(caloEmEtInEE);
333  mCaloEmEtInHF->Fill(caloEmEtInHF);
334  mCaloSETInpHF->Fill(caloSETInpHF);
335  mCaloSETInmHF->Fill(caloSETInmHF);
336  }
337  if (isGenMET)
338  {
339  const GenMET *genmet;
340  // Get Generated MET
341  genmet = &(genMETs->front());
342 
343  const double NeutralEMEtFraction = genmet->NeutralEMEtFraction() ;
344  const double NeutralHadEtFraction = genmet->NeutralHadEtFraction() ;
345  const double ChargedEMEtFraction = genmet->ChargedEMEtFraction () ;
346  const double ChargedHadEtFraction = genmet->ChargedHadEtFraction();
347  const double MuonEtFraction = genmet->MuonEtFraction() ;
348  const double InvisibleEtFraction = genmet->InvisibleEtFraction() ;
349 
350  mNeutralEMEtFraction->Fill( NeutralEMEtFraction );
351  mNeutralHadEtFraction->Fill( NeutralHadEtFraction );
352  mChargedEMEtFraction->Fill( ChargedEMEtFraction );
353  mChargedHadEtFraction->Fill( ChargedHadEtFraction );
354  mMuonEtFraction->Fill( MuonEtFraction );
355  mInvisibleEtFraction->Fill( InvisibleEtFraction );
356  }
357  if(isPFMET)
358  {
359  const reco::PFMET * pfmet = &(pfMETs->front());
360  mPFphotonEtFraction->Fill(pfmet->photonEtFraction());
361  mPFphotonEt->Fill(pfmet->photonEt());
362  mPFneutralHadronEtFraction->Fill(pfmet->neutralHadronEtFraction());
363  mPFneutralHadronEt->Fill(pfmet->neutralHadronEt());
364  mPFelectronEtFraction->Fill(pfmet->electronEtFraction());
365  mPFelectronEt->Fill(pfmet->electronEt());
366  mPFchargedHadronEtFraction->Fill(pfmet->chargedHadronEtFraction());
367  mPFchargedHadronEt->Fill(pfmet->chargedHadronEt());
368  mPFmuonEtFraction->Fill(pfmet->muonEtFraction());
369  mPFmuonEt->Fill(pfmet->muonEt());
370  mPFHFHadronEtFraction->Fill(pfmet->HFHadronEtFraction());
371  mPFHFHadronEt->Fill(pfmet->HFHadronEt());
372  mPFHFEMEtFraction->Fill(pfmet->HFEMEtFraction());
373  mPFHFEMEt->Fill(pfmet->HFEMEt());
374  // Reconstructed MET Information
375 
376  }
377 
378  //This is so dirty I could cry. It should be called only ONCE in endJob. But the MonitorElements don't exist then any more.
379  FillMETRes();
380 }
381 
382 //void METTester::endRun(const edm::Run& iRun, const edm::EventSetup& iSetup)
383 void METTester::FillMETRes()
384 {
385  mMETDifference_GenMETTrue_METResolution->setBinContent(1, mMETDifference_GenMETTrue_MET0to20->getMean());
386  mMETDifference_GenMETTrue_METResolution->setBinContent(2, mMETDifference_GenMETTrue_MET20to40->getMean());
387  mMETDifference_GenMETTrue_METResolution->setBinContent(3, mMETDifference_GenMETTrue_MET40to60->getMean());
388  mMETDifference_GenMETTrue_METResolution->setBinContent(4, mMETDifference_GenMETTrue_MET60to80->getMean());
389  mMETDifference_GenMETTrue_METResolution->setBinContent(5, mMETDifference_GenMETTrue_MET80to100->getMean());
390  mMETDifference_GenMETTrue_METResolution->setBinContent(6, mMETDifference_GenMETTrue_MET100to150->getMean());
391  mMETDifference_GenMETTrue_METResolution->setBinContent(7, mMETDifference_GenMETTrue_MET150to200->getMean());
392  mMETDifference_GenMETTrue_METResolution->setBinContent(8, mMETDifference_GenMETTrue_MET200to300->getMean());
393  mMETDifference_GenMETTrue_METResolution->setBinContent(9, mMETDifference_GenMETTrue_MET300to400->getMean());
394  mMETDifference_GenMETTrue_METResolution->setBinContent(10, mMETDifference_GenMETTrue_MET400to500->getMean());
395 
396  //the error computation should be done in a postProcessor in the harvesting step otherwise the histograms will be just summed
397  mMETDifference_GenMETTrue_METResolution->setBinError(1, mMETDifference_GenMETTrue_MET0to20->getRMS());
398  mMETDifference_GenMETTrue_METResolution->setBinError(2, mMETDifference_GenMETTrue_MET20to40->getRMS());
399  mMETDifference_GenMETTrue_METResolution->setBinError(3, mMETDifference_GenMETTrue_MET40to60->getRMS());
400  mMETDifference_GenMETTrue_METResolution->setBinError(4, mMETDifference_GenMETTrue_MET60to80->getRMS());
401  mMETDifference_GenMETTrue_METResolution->setBinError(5, mMETDifference_GenMETTrue_MET80to100->getRMS());
402  mMETDifference_GenMETTrue_METResolution->setBinError(6, mMETDifference_GenMETTrue_MET100to150->getRMS());
403  mMETDifference_GenMETTrue_METResolution->setBinError(7, mMETDifference_GenMETTrue_MET150to200->getRMS());
404  mMETDifference_GenMETTrue_METResolution->setBinError(8, mMETDifference_GenMETTrue_MET200to300->getRMS());
405  mMETDifference_GenMETTrue_METResolution->setBinError(9, mMETDifference_GenMETTrue_MET300to400->getRMS());
406  mMETDifference_GenMETTrue_METResolution->setBinError(10, mMETDifference_GenMETTrue_MET400to500->getRMS());
407 
408 }
409 
410 //determines if track is "good" - i.e. passes quality and kinematic cuts
411 bool METTester::isGoodTrack( const reco::TrackRef track, float d0corr ) {
412 
413  if( fabs( d0corr ) > maxd0_ ) return false;
414  if( track->numberOfValidHits() < minhits_ ) return false;
415  if( track->normalizedChi2() > maxchi2_ ) return false;
416  if( fabs( track->eta() ) > maxeta_ ) return false;
417  if( track->pt() > maxpt_ ) return false;
418  if( (track->ptError() / track->pt()) > maxPtErr_ ) return false;
419 
420  int cut = 0;
421  for( unsigned int i = 0; i < trkQuality_.size(); i++ ) {
422 
423  cut |= (1 << trkQuality_.at(i));
424  }
425 
426  if( !( ( track->qualityMask() & cut ) == cut ) ) return false;
427 
428  bool isGoodAlgo = false;
429  if( trkAlgos_.size() == 0 ) isGoodAlgo = true;
430  for( unsigned int i = 0; i < trkAlgos_.size(); i++ ) {
431 
432  if( track->algo() == trkAlgos_.at(i) ) isGoodAlgo = true;
433  }
434 
435  if( !isGoodAlgo ) return false;
436 
437  return true;
438  }
439 
440 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
reco::CaloMET::hadEtInHE
double hadEtInHE() const
Definition: CaloMET.h:51
GenMETCollection
Collection of Gen MET.
reco::PFMET::HFEMEtFraction
double HFEMEtFraction() const
Definition: PFMET.h:47
reco::CaloMET::hadEtInHF
double hadEtInHF() const
Definition: CaloMET.h:53
METTester::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: METTester.cc:220
reco::CaloMET::CaloSETInmHF
double CaloSETInmHF() const
Definition: CaloMET.h:65
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:446
METTester::FillMETRes
void FillMETRes()
Definition: METTester.cc:383
reco::CaloMET::maxEtInHadTowers
double maxEtInHadTowers() const
Definition: CaloMET.h:40
reco::PFMET::neutralHadronEtFraction
double neutralHadronEtFraction() const
Definition: PFMET.h:32
reco::PFMET::muonEt
double muonEt() const
Definition: PFMET.h:42
genMetCalo_cfi.genMetCalo
tuple genMetCalo
____________________________________________________________________________||
Definition: genMetCalo_cfi.py:4
reco::GenMET::MuonEtFraction
double MuonEtFraction() const
Definition: GenMET.h:57
reco::CaloMET::maxEtInEmTowers
double maxEtInEmTowers() const
Definition: CaloMET.h:38
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
reco::LeafCandidate::py
virtual double py() const GCC11_FINAL
y coordinate of momentum vector
Definition: LeafCandidate.h:174
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
reco::LeafCandidate::phi
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
Definition: LeafCandidate.h:180
reco::MET::mEtSig
double mEtSig() const
Definition: MET.h:58
METTester::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: METTester.cc:140
reco::PFMET::HFHadronEtFraction
double HFHadronEtFraction() const
Definition: PFMET.h:44
iEvent
int iEvent
Definition: GenABIO.cc:230
reco::MET::sumEt
double sumEt() const
Definition: MET.h:56
reco::PFMET::muonEtFraction
double muonEtFraction() const
Definition: PFMET.h:41
reco::PFMET::HFHadronEt
double HFHadronEt() const
Definition: PFMET.h:45
reco::MET
Definition: MET.h:42
reco::PFMET::photonEtFraction
double photonEtFraction() const
Definition: PFMET.h:29
reco::GenMET::ChargedEMEtFraction
double ChargedEMEtFraction() const
Definition: GenMET.h:39
genMetTrue_cff.genMetTrue
tuple genMetTrue
Definition: genMetTrue_cff.py:21
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
reco::LeafCandidate::px
virtual double px() const GCC11_FINAL
x coordinate of momentum vector
Definition: LeafCandidate.h:172
reco::CaloMET::emEtInEB
double emEtInEB() const
Definition: CaloMET.h:55
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
reco::CaloMET::CaloSETInpHF
double CaloSETInpHF() const
Definition: CaloMET.h:63
reco::GenMET::InvisibleEtFraction
double InvisibleEtFraction() const
Definition: GenMET.h:63
reco::PFMET::HFEMEt
double HFEMEt() const
Definition: PFMET.h:48
reco::CaloMET::hadEtInHO
double hadEtInHO() const
Definition: CaloMET.h:49
reco::PFMET::electronEt
double electronEt() const
Definition: PFMET.h:36
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
reco::CaloMET::etFractionHadronic
double etFractionHadronic() const
Definition: CaloMET.h:42
edm::Handle::product
T const * product() const
Definition: Handle.h:81
reco::PFMET::photonEt
double photonEt() const
Definition: PFMET.h:30
reco::GenMET::NeutralEMEtFraction
double NeutralEMEtFraction() const
Definition: GenMET.h:33
reco::CaloMET::emEtInEE
double emEtInEE() const
Definition: CaloMET.h:57
reco::PFMET::electronEtFraction
double electronEtFraction() const
Definition: PFMET.h:35
reco::PFMET::chargedHadronEtFraction
double chargedHadronEtFraction() const
Definition: PFMET.h:38
METTester::METTester
METTester(const edm::ParameterSet &)
Definition: METTester.cc:32
reco::GenMET::NeutralHadEtFraction
double NeutralHadEtFraction() const
Definition: GenMET.h:45
METTester::isGoodTrack
bool isGoodTrack(const reco::TrackRef, float d0corr)
Definition: METTester.cc:411
reco::CaloMET::emEtInHF
double emEtInHF() const
Definition: CaloMET.h:59
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
reco::PFMET::neutralHadronEt
double neutralHadronEt() const
Definition: PFMET.h:33
caloMETBenchmarkGeneric_cfi.met
string met
Definition: caloMETBenchmarkGeneric_cfi.py:3
reco::GenMET::ChargedHadEtFraction
double ChargedHadEtFraction() const
Definition: GenMET.h:51
reco::LeafCandidate::pt
virtual float pt() const GCC11_FINAL
transverse momentum
Definition: LeafCandidate.h:178
reco::CaloMET::hadEtInHB
double hadEtInHB() const
Definition: CaloMET.h:47
reco::CaloMET::emEtFraction
double emEtFraction() const
Definition: CaloMET.h:45
reco::PFMET::chargedHadronEt
double chargedHadronEt() const
Definition: PFMET.h:39
edm::Run
Definition: Run.h:41
GOODCOLL_filter_cfg.cut
tuple cut
Definition: GOODCOLL_filter_cfg.py:65
Generated for CMSSW Reference Manual by doxygen 1.8.5