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 
33 {
34 
35  inputMETLabel_ =iConfig.getParameter<edm::InputTag>("InputMETLabel");
36  METType_ =iConfig.getUntrackedParameter<std::string>("METType");
37 
38  std::string inputMETCollectionLabel(inputMETLabel_.label());
39 
40  isCaloMET = (std::string("calo")==METType_);
41  isPFMET = (std::string("pf") ==METType_);
42  isGenMET = (std::string("gen") ==METType_);
43  isMiniAODMET = (std::string("miniaod") ==METType_);
44 
45  pvToken_ = consumes<std::vector<reco::Vertex> >(iConfig.getParameter<edm::InputTag>("PrimaryVertices"));
46  if (isCaloMET) caloMETsToken_ = consumes<reco::CaloMETCollection> (inputMETLabel_);
47  if (isPFMET) pfMETsToken_ = consumes<reco::PFMETCollection> (inputMETLabel_);
48  if(isMiniAODMET) patMETToken_ = consumes<pat::METCollection> (inputMETLabel_);
49  if (isGenMET) genMETsToken_ = consumes<reco::GenMETCollection> (inputMETLabel_);
50  if(!isMiniAODMET){
51  genMETsTrueToken_ = consumes<reco::GenMETCollection> (edm::InputTag("genMetTrue"));
52  genMETsCaloToken_ = consumes<reco::GenMETCollection> (edm::InputTag("genMetCalo"));
53  }
54 
55  //Events variables
56  mNvertex = 0;
57 
58  //Common variables
59  mMEx = 0;
60  mMEy = 0;
61  mMETSig = 0;
62  mMET = 0;
63  mMETFine = 0;
64  mMET_Nvtx = 0;
65  mMETPhi = 0;
66  mSumET = 0;
67  mMETDifference_GenMETTrue = 0;
68  mMETDeltaPhi_GenMETTrue = 0;
69  mMETDifference_GenMETCalo = 0;
70  mMETDeltaPhi_GenMETCalo = 0;
71 
72  //CaloMET variables
73  mCaloMaxEtInEmTowers = 0;
74  mCaloMaxEtInHadTowers = 0;
75  mCaloEtFractionHadronic = 0;
76  mCaloEmEtFraction = 0;
77  mCaloHadEtInHB = 0;
78  mCaloHadEtInHO = 0;
79  mCaloHadEtInHE = 0;
80  mCaloHadEtInHF = 0;
81  mCaloEmEtInHF = 0;
82  mCaloSETInpHF = 0;
83  mCaloSETInmHF = 0;
84  mCaloEmEtInEE = 0;
85  mCaloEmEtInEB = 0;
86 
87 
88  //GenMET variables
89  mNeutralEMEtFraction=0;
90  mNeutralHadEtFraction=0;
91  mChargedEMEtFraction=0;
92  mChargedHadEtFraction=0;
93  mMuonEtFraction=0;
94  mInvisibleEtFraction=0;
95 
96  //MET variables
97 
98  //PFMET variables
99  mMETDifference_GenMETTrue_MET0to20=0;
100  mMETDifference_GenMETTrue_MET20to40=0;
101  mMETDifference_GenMETTrue_MET40to60=0;
102  mMETDifference_GenMETTrue_MET60to80=0;
103  mMETDifference_GenMETTrue_MET80to100=0;
104  mMETDifference_GenMETTrue_MET100to150=0;
105  mMETDifference_GenMETTrue_MET150to200=0;
106  mMETDifference_GenMETTrue_MET200to300=0;
107  mMETDifference_GenMETTrue_MET300to400=0;
108  mMETDifference_GenMETTrue_MET400to500=0;
109  mMETDifference_GenMETTrue_MET500=0;
110 
111 
112 }
114  edm::Run const & iRun,
115  edm::EventSetup const & /* iSetup */)
116 {
117  ibooker.setCurrentFolder("JetMET/METValidation/"+inputMETLabel_.label()) ;
118 
119  mNvertex = ibooker.book1D("Nvertex","Nvertex",80,0,80);
120  mMEx = ibooker.book1D("MEx","MEx",160,-800,800);
121  mMEy = ibooker.book1D("MEy","MEy",160,-800,800);
122  mMETSig = ibooker.book1D("METSig","METSig",25,0,24.5);
123  mMET = ibooker.book1D("MET", "MET (20 GeV binning)" , 100,0,2000);
124  mMETFine = ibooker.book1D("METFine", "MET (2 GeV binning)" , 1000,0,2000);
125  mMET_Nvtx = ibooker.bookProfile("MET_Nvtx", "MET vs. nvtx", 60, 0., 60., 0., 2000., " ");
126  mMETPhi = ibooker.book1D("METPhi","METPhi",80,-4,4);
127  mSumET = ibooker.book1D("SumET" , "SumET" , 200,0,4000); //10GeV
128  mMETDifference_GenMETTrue = ibooker.book1D("METDifference_GenMETTrue","METDifference_GenMETTrue", 500,-500,500);
129  mMETDeltaPhi_GenMETTrue = ibooker.book1D("METDeltaPhi_GenMETTrue","METDeltaPhi_GenMETTrue", 80,0,4);
130  if(!isMiniAODMET){
131  mMETDifference_GenMETCalo = ibooker.book1D("METDifference_GenMETCalo","METDifference_GenMETCalo", 500,-500,500);
132  mMETDeltaPhi_GenMETCalo = ibooker.book1D("METDeltaPhi_GenMETCalo","METDeltaPhi_GenMETCalo", 80,0,4);
133  }
134  if(!isGenMET){
135  mMETDifference_GenMETTrue_MET0to20 = ibooker.book1D("METResolution_GenMETTrue_MET0to20" , "METResolution_GenMETTrue_MET0to20" , 500,-500,500);
136  mMETDifference_GenMETTrue_MET20to40 = ibooker.book1D("METResolution_GenMETTrue_MET20to40" , "METResolution_GenMETTrue_MET20to40" , 500,-500,500);
137  mMETDifference_GenMETTrue_MET40to60 = ibooker.book1D("METResolution_GenMETTrue_MET40to60" , "METResolution_GenMETTrue_MET40to60" , 500,-500,500);
138  mMETDifference_GenMETTrue_MET60to80 = ibooker.book1D("METResolution_GenMETTrue_MET60to80" , "METResolution_GenMETTrue_MET60to80" , 500,-500,500);
139  mMETDifference_GenMETTrue_MET80to100 = ibooker.book1D("METResolution_GenMETTrue_MET80to100" , "METResolution_GenMETTrue_MET80to100" , 500,-500,500);
140  mMETDifference_GenMETTrue_MET100to150 = ibooker.book1D("METResolution_GenMETTrue_MET100to150", "METResolution_GenMETTrue_MET100to150", 500,-500,500);
141  mMETDifference_GenMETTrue_MET150to200 = ibooker.book1D("METResolution_GenMETTrue_MET150to200", "METResolution_GenMETTrue_MET150to200", 500,-500,500);
142  mMETDifference_GenMETTrue_MET200to300 = ibooker.book1D("METResolution_GenMETTrue_MET200to300", "METResolution_GenMETTrue_MET200to300", 500,-500,500);
143  mMETDifference_GenMETTrue_MET300to400 = ibooker.book1D("METResolution_GenMETTrue_MET300to400", "METResolution_GenMETTrue_MET300to400", 500,-500,500);
144  mMETDifference_GenMETTrue_MET400to500 = ibooker.book1D("METResolution_GenMETTrue_MET400to500", "METResolution_GenMETTrue_MET400to500", 500,-500,500);
145  mMETDifference_GenMETTrue_MET500 = ibooker.book1D("METResolution_GenMETTrue_MET500", "METResolution_GenMETTrue_MET500", 500,-500,500);
146  }
147  if ( isCaloMET) {
148  mCaloMaxEtInEmTowers = ibooker.book1D("CaloMaxEtInEmTowers","CaloMaxEtInEmTowers",300,0,1500); //5GeV
149  mCaloMaxEtInHadTowers = ibooker.book1D("CaloMaxEtInHadTowers","CaloMaxEtInHadTowers",300,0,1500); //5GeV
150  mCaloEtFractionHadronic = ibooker.book1D("CaloEtFractionHadronic","CaloEtFractionHadronic",100,0,1);
151  mCaloEmEtFraction = ibooker.book1D("CaloEmEtFraction","CaloEmEtFraction",100,0,1);
152  mCaloHadEtInHB = ibooker.book1D("CaloHadEtInHB","CaloHadEtInHB", 200, 0, 2000); //5GeV
153  mCaloHadEtInHE = ibooker.book1D("CaloHadEtInHE","CaloHadEtInHE", 100, 0, 500); //5GeV
154  mCaloHadEtInHO = ibooker.book1D("CaloHadEtInHO","CaloHadEtInHO", 100, 0, 200); //5GeV
155  mCaloHadEtInHF = ibooker.book1D("CaloHadEtInHF","CaloHadEtInHF", 100, 0, 200); //5GeV
156  mCaloSETInpHF = ibooker.book1D("CaloSETInpHF","CaloSETInpHF",100, 0, 500);
157  mCaloSETInmHF = ibooker.book1D("CaloSETInmHF","CaloSETInmHF",100, 0, 500);
158  mCaloEmEtInEE = ibooker.book1D("CaloEmEtInEE","CaloEmEtInEE",100, 0, 500); //5GeV
159  mCaloEmEtInEB = ibooker.book1D("CaloEmEtInEB","CaloEmEtInEB",100, 0, 500); //5GeV
160  mCaloEmEtInHF = ibooker.book1D("CaloEmEtInHF","CaloEmEtInHF",100, 0, 500); //5GeV
161  }
162 
163  if(isGenMET){
164  mNeutralEMEtFraction = ibooker.book1D("GenNeutralEMEtFraction", "GenNeutralEMEtFraction", 120, 0.0, 1.2 );
165  mNeutralHadEtFraction = ibooker.book1D("GenNeutralHadEtFraction", "GenNeutralHadEtFraction", 120, 0.0, 1.2 );
166  mChargedEMEtFraction = ibooker.book1D("GenChargedEMEtFraction", "GenChargedEMEtFraction", 120, 0.0, 1.2);
167  mChargedHadEtFraction = ibooker.book1D("GenChargedHadEtFraction", "GenChargedHadEtFraction", 120, 0.0,1.2);
168  mMuonEtFraction = ibooker.book1D("GenMuonEtFraction", "GenMuonEtFraction", 120, 0.0, 1.2 );
169  mInvisibleEtFraction = ibooker.book1D("GenInvisibleEtFraction", "GenInvisibleEtFraction", 120, 0.0, 1.2 );
170  }
171 
172  if (isPFMET || isMiniAODMET){
173  mPFphotonEtFraction = ibooker.book1D("photonEtFraction", "photonEtFraction", 100, 0, 1);
174  mPFneutralHadronEtFraction = ibooker.book1D("neutralHadronEtFraction", "neutralHadronEtFraction", 100, 0, 1);
175  mPFelectronEtFraction = ibooker.book1D("electronEtFraction", "electronEtFraction", 100, 0, 1);
176  mPFchargedHadronEtFraction = ibooker.book1D("chargedHadronEtFraction", "chargedHadronEtFraction", 100, 0, 1);
177  mPFHFHadronEtFraction = ibooker.book1D("HFHadronEtFraction", "HFHadronEtFraction", 100, 0, 1);
178  mPFmuonEtFraction = ibooker.book1D("muonEtFraction", "muonEtFraction", 100, 0, 1);
179  mPFHFEMEtFraction = ibooker.book1D("HFEMEtFraction", "HFEMEtFraction", 100, 0, 1);
180 
181  if(!isMiniAODMET){
182  mPFphotonEt = ibooker.book1D("photonEt", "photonEt", 100, 0, 1000);
183  mPFneutralHadronEt = ibooker.book1D("neutralHadronEt", "neutralHadronEt", 100, 0, 1000);
184  mPFelectronEt = ibooker.book1D("electronEt", "electronEt", 100, 0, 1000);
185  mPFchargedHadronEt = ibooker.book1D("chargedHadronEt", "chargedHadronEt", 100, 0, 1000);
186  mPFmuonEt = ibooker.book1D("muonEt", "muonEt", 100, 0, 1000);
187  mPFHFHadronEt = ibooker.book1D("HFHadronEt", "HFHadronEt", 100, 0, 500);
188  mPFHFEMEt = ibooker.book1D("HFEMEt", "HFEMEt", 100, 0, 300);
189  }
190 
191  }
192 }
193 
194 
196 { //int counter(0);
197 
199  iEvent.getByToken(pvToken_, pvHandle);
200  if (! pvHandle.isValid())
201  {
202  std::cout << __FUNCTION__ << ":" << __LINE__ << ":pvHandle handle not found!" << std::endl;
203  assert(false);
204  }
205  const int nvtx = pvHandle->size();
206  mNvertex->Fill(nvtx);
207  //Collections for all MET collections
208 
213 
214  if (isCaloMET) iEvent.getByToken(caloMETsToken_, caloMETs);
215  if (isPFMET) iEvent.getByToken(pfMETsToken_, pfMETs);
216  if (isGenMET) iEvent.getByToken(genMETsToken_, genMETs);
217  if (isMiniAODMET) iEvent.getByToken(patMETToken_, patMET);
218  if ((isCaloMET) and !caloMETs.isValid()) return;
219  if ((isPFMET) and !pfMETs.isValid()) return;
220  if ((isGenMET) and !genMETs.isValid()) return;
221  if ((isMiniAODMET) and !patMET.isValid()) return;
222 
223  reco::MET met;
224  if (isCaloMET) { met = caloMETs->front();}
225  if (isPFMET) { met = pfMETs->front() ;}
226  if (isGenMET) { met = genMETs->front() ;}
227  if (isMiniAODMET) { met = patMET->front() ;}
228 
229  const double SumET = met.sumEt();
230  const double METSig = met.mEtSig();
231  const double MET = met.pt();
232  const double MEx = met.px();
233  const double MEy = met.py();
234  const double METPhi = met.phi();
235  mMEx->Fill(MEx);
236  mMEy->Fill(MEy);
237  mMET->Fill(MET);
238  mMETFine->Fill(MET);
239  mMET_Nvtx->Fill((double)nvtx, MET);
240  mMETPhi->Fill(METPhi);
241  mSumET->Fill(SumET);
242  mMETSig->Fill(METSig);
243 
244  // Get Generated MET for Resolution plots
246  bool isvalidgenmet=false;
247 
248  if(!isMiniAODMET){
250  iEvent.getByToken(genMETsTrueToken_, genTrue);
251  if (genTrue.isValid()) {
252  isvalidgenmet=true;
253  const GenMETCollection *genmetcol = genTrue.product();
254  genMetTrue = &(genmetcol->front());
255  }
256  }else{
257  genMetTrue=patMET->front().genMET();
258  isvalidgenmet=true;
259  }
260 
261  if(isvalidgenmet){
262  double genMET = genMetTrue->pt();
263  double genMETPhi = genMetTrue->phi();
264 
265  mMETDifference_GenMETTrue->Fill( MET - genMET );
266  mMETDeltaPhi_GenMETTrue->Fill( TMath::ACos( TMath::Cos( METPhi - genMETPhi ) ) );
267 
268  if(!isGenMET){
269  //pfMET resolution in pfMET bins : Sam, Feb, 2012
270  if (MET > 0 && MET < 20) mMETDifference_GenMETTrue_MET0to20->Fill( MET - genMET );
271  else if (MET > 20 && MET < 40) mMETDifference_GenMETTrue_MET20to40->Fill( MET - genMET );
272  else if (MET > 40 && MET < 60) mMETDifference_GenMETTrue_MET40to60->Fill( MET - genMET );
273  else if (MET > 60 && MET < 80) mMETDifference_GenMETTrue_MET60to80->Fill( MET - genMET );
274  else if (MET > 80 && MET <100) mMETDifference_GenMETTrue_MET80to100->Fill( MET - genMET );
275  else if (MET >100 && MET <150) mMETDifference_GenMETTrue_MET100to150->Fill( MET - genMET );
276  else if (MET >150 && MET <200) mMETDifference_GenMETTrue_MET150to200->Fill( MET - genMET );
277  else if (MET >200 && MET <300) mMETDifference_GenMETTrue_MET200to300->Fill( MET - genMET );
278  else if (MET >300 && MET <400) mMETDifference_GenMETTrue_MET300to400->Fill( MET - genMET );
279  else if (MET >400 && MET <500) mMETDifference_GenMETTrue_MET400to500->Fill( MET - genMET );
280  else if (MET >500) mMETDifference_GenMETTrue_MET500->Fill( MET - genMET );
281 
282  } else {
283  edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task: genMetTrue";
284  }
285  }
286  if(!isMiniAODMET){
288  iEvent.getByToken(genMETsCaloToken_, genCalo);
289  if (genCalo.isValid()) {
290  const GenMETCollection *genmetcol = genCalo.product();
291  const GenMET *genMetCalo = &(genmetcol->front());
292  const double genMET = genMetCalo->pt();
293  const double genMETPhi = genMetCalo->phi();
294 
295  mMETDifference_GenMETCalo->Fill( MET - genMET );
296  mMETDeltaPhi_GenMETCalo->Fill( TMath::ACos( TMath::Cos( METPhi - genMETPhi ) ) );
297  } else {
298  edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task: genMetCalo";
299  }
300  }
301  if ( isCaloMET) {
302  const reco::CaloMET * calomet = &(caloMETs->front());
303  // ==========================================================
304  // Reconstructed MET Information
305  const double caloMaxEtInEMTowers = calomet->maxEtInEmTowers();
306  const double caloMaxEtInHadTowers = calomet->maxEtInHadTowers();
307  const double caloEtFractionHadronic = calomet->etFractionHadronic();
308  const double caloEmEtFraction = calomet->emEtFraction();
309  const double caloHadEtInHB = calomet->hadEtInHB();
310  const double caloHadEtInHO = calomet->hadEtInHO();
311  const double caloHadEtInHE = calomet->hadEtInHE();
312  const double caloHadEtInHF = calomet->hadEtInHF();
313  const double caloEmEtInEB = calomet->emEtInEB();
314  const double caloEmEtInEE = calomet->emEtInEE();
315  const double caloEmEtInHF = calomet->emEtInHF();
316  const double caloSETInpHF = calomet->CaloSETInpHF();
317  const double caloSETInmHF = calomet->CaloSETInmHF();
318 
319  mCaloMaxEtInEmTowers->Fill(caloMaxEtInEMTowers);
320  mCaloMaxEtInHadTowers->Fill(caloMaxEtInHadTowers);
321  mCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
322  mCaloEmEtFraction->Fill(caloEmEtFraction);
323  mCaloHadEtInHB->Fill(caloHadEtInHB);
324  mCaloHadEtInHO->Fill(caloHadEtInHO);
325  mCaloHadEtInHE->Fill(caloHadEtInHE);
326  mCaloHadEtInHF->Fill(caloHadEtInHF);
327  mCaloEmEtInEB->Fill(caloEmEtInEB);
328  mCaloEmEtInEE->Fill(caloEmEtInEE);
329  mCaloEmEtInHF->Fill(caloEmEtInHF);
330  mCaloSETInpHF->Fill(caloSETInpHF);
331  mCaloSETInmHF->Fill(caloSETInmHF);
332  }
333  if (isGenMET)
334  {
335  const GenMET *genmet;
336  // Get Generated MET
337  genmet = &(genMETs->front());
338 
339  const double NeutralEMEtFraction = genmet->NeutralEMEtFraction() ;
340  const double NeutralHadEtFraction = genmet->NeutralHadEtFraction() ;
341  const double ChargedEMEtFraction = genmet->ChargedEMEtFraction () ;
342  const double ChargedHadEtFraction = genmet->ChargedHadEtFraction();
343  const double MuonEtFraction = genmet->MuonEtFraction() ;
344  const double InvisibleEtFraction = genmet->InvisibleEtFraction() ;
345 
346  mNeutralEMEtFraction->Fill( NeutralEMEtFraction );
347  mNeutralHadEtFraction->Fill( NeutralHadEtFraction );
348  mChargedEMEtFraction->Fill( ChargedEMEtFraction );
349  mChargedHadEtFraction->Fill( ChargedHadEtFraction );
350  mMuonEtFraction->Fill( MuonEtFraction );
351  mInvisibleEtFraction->Fill( InvisibleEtFraction );
352  }
353  if(isPFMET)
354  {
355  const reco::PFMET * pfmet = &(pfMETs->front());
356  mPFphotonEtFraction->Fill(pfmet->photonEtFraction());
357  mPFphotonEt->Fill(pfmet->photonEt());
358  mPFneutralHadronEtFraction->Fill(pfmet->neutralHadronEtFraction());
359  mPFneutralHadronEt->Fill(pfmet->neutralHadronEt());
360  mPFelectronEtFraction->Fill(pfmet->electronEtFraction());
361  mPFelectronEt->Fill(pfmet->electronEt());
362  mPFchargedHadronEtFraction->Fill(pfmet->chargedHadronEtFraction());
363  mPFchargedHadronEt->Fill(pfmet->chargedHadronEt());
364  mPFmuonEtFraction->Fill(pfmet->muonEtFraction());
365  mPFmuonEt->Fill(pfmet->muonEt());
366  mPFHFHadronEtFraction->Fill(pfmet->HFHadronEtFraction());
367  mPFHFHadronEt->Fill(pfmet->HFHadronEt());
368  mPFHFEMEtFraction->Fill(pfmet->HFEMEtFraction());
369  mPFHFEMEt->Fill(pfmet->HFEMEt());
370  // Reconstructed MET Information
371 
372  }
373  if(isMiniAODMET){
374  const pat::MET *patmet = &(patMET->front());
375  if(patmet->isPFMET()){
376  mPFphotonEtFraction->Fill(patmet->NeutralEMFraction());
377  mPFneutralHadronEtFraction->Fill(patmet->NeutralHadEtFraction());
378  mPFelectronEtFraction->Fill(patmet->ChargedEMEtFraction());
379  mPFchargedHadronEtFraction->Fill(patmet->ChargedHadEtFraction());
380  mPFmuonEtFraction->Fill(patmet->MuonEtFraction());
381  mPFHFHadronEtFraction->Fill(patmet->Type6EtFraction());//HFHadrons
382  mPFHFEMEtFraction->Fill(patmet->Type7EtFraction());//HFEMEt
383  }
384  }
385 }
386 
387 
Analysis-level MET class.
Definition: MET.h:43
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
double hadEtInHE() const
Definition: CaloMET.h:51
Collection of Gen MET.
double HFEMEtFraction() const
Definition: PFMET.h:47
double hadEtInHF() const
Definition: CaloMET.h:53
double ChargedHadEtFraction() const
Definition: MET.h:141
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: METTester.cc:195
double CaloSETInmHF() const
Definition: CaloMET.h:65
double MuonEtFraction() const
Definition: MET.h:142
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
double maxEtInHadTowers() const
Definition: CaloMET.h:40
double neutralHadronEtFraction() const
Definition: PFMET.h:32
assert(m_qm.get())
double muonEt() const
Definition: PFMET.h:42
tuple genMetCalo
____________________________________________________________________________||
virtual double phi() const final
momentum azimuthal angle
double MuonEtFraction() const
Definition: GenMET.h:57
#define NULL
Definition: scimark2.h:8
double NeutralEMFraction() const
Definition: MET.h:138
double ChargedEMEtFraction() const
Definition: MET.h:140
double maxEtInEmTowers() const
Definition: CaloMET.h:38
double Type6EtFraction() const
Definition: MET.h:143
double Type7EtFraction() const
Definition: MET.h:144
double NeutralHadEtFraction() const
Definition: MET.h:139
bool isPFMET() const
True if this pat::MET was made from a reco::pfMET.
Definition: MET.h:90
double mEtSig() const
Definition: MET.h:58
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: METTester.cc:113
double HFHadronEtFraction() const
Definition: PFMET.h:44
int iEvent
Definition: GenABIO.cc:230
double sumEt() const
Definition: MET.h:56
double muonEtFraction() const
Definition: PFMET.h:41
double HFHadronEt() const
Definition: PFMET.h:45
Definition: MET.h:42
double photonEtFraction() const
Definition: PFMET.h:29
double ChargedEMEtFraction() const
Definition: GenMET.h:39
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
double emEtInEB() const
Definition: CaloMET.h:55
virtual double py() const final
y coordinate of momentum vector
bool isValid() const
Definition: HandleBase.h:75
double CaloSETInpHF() const
Definition: CaloMET.h:63
double InvisibleEtFraction() const
Definition: GenMET.h:63
double HFEMEt() const
Definition: PFMET.h:48
double hadEtInHO() const
Definition: CaloMET.h:49
double electronEt() const
Definition: PFMET.h:36
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
double etFractionHadronic() const
Definition: CaloMET.h:42
T const * product() const
Definition: Handle.h:81
double photonEt() const
Definition: PFMET.h:30
double NeutralEMEtFraction() const
Definition: GenMET.h:33
double emEtInEE() const
Definition: CaloMET.h:57
double electronEtFraction() const
Definition: PFMET.h:35
double chargedHadronEtFraction() const
Definition: PFMET.h:38
METTester(const edm::ParameterSet &)
Definition: METTester.cc:32
double NeutralHadEtFraction() const
Definition: GenMET.h:45
virtual double px() const final
x coordinate of momentum vector
double emEtInHF() const
Definition: CaloMET.h:59
tuple cout
Definition: gather_cfg.py:145
double neutralHadronEt() const
Definition: PFMET.h:33
double ChargedHadEtFraction() const
Definition: GenMET.h:51
double hadEtInHB() const
Definition: CaloMET.h:47
double emEtFraction() const
Definition: CaloMET.h:45
double chargedHadronEt() const
Definition: PFMET.h:39
Definition: Run.h:43
virtual double pt() const final
transverse momentum