CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
CaloJetTester.cc
Go to the documentation of this file.
1 // Producer for validation histograms for CaloJet objects
2 // F. Ratnikov, Sept. 7, 2006
3 // Modified by J F Novak July 10, 2008
4 // $Id: CaloJetTester.cc,v 1.23 2010/06/16 19:06:12 srappocc Exp $
5 
6 #include "FWCore/MessageLogger/interface/MessageLogger.h"
7 #include "FWCore/Framework/interface/Event.h"
8 #include "DataFormats/Common/interface/Handle.h"
9 #include "FWCore/ParameterSet/interface/ParameterSet.h"
10 #include "FWCore/ServiceRegistry/interface/Service.h"
11 
12 #include "DataFormats/Math/interface/deltaR.h"
13 
14 #include "DQMServices/Core/interface/DQMStore.h"
15 #include "DQMServices/Core/interface/MonitorElement.h"
16 
17 #include "DataFormats/JetReco/interface/CaloJet.h"
18 #include "DataFormats/JetReco/interface/CaloJetCollection.h"
19 #include "DataFormats/JetReco/interface/GenJet.h"
20 #include "DataFormats/JetReco/interface/GenJetCollection.h"
21 
22 //I don't know which of these I actually need yet
23 #include "DataFormats/METReco/interface/CaloMET.h"
24 #include "DataFormats/METReco/interface/CaloMETCollection.h"
25 #include "DataFormats/METReco/interface/GenMET.h"
26 #include "DataFormats/METReco/interface/GenMETCollection.h"
27 #include "DataFormats/METReco/interface/MET.h"
28 #include "DataFormats/METReco/interface/METCollection.h"
29 
30 #include "RecoJets/JetProducers/interface/JetMatchingTools.h"
31 
32 #include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
33 
34 #include "CaloJetTester.h"
35 
36 #include <cmath>
37 
38 using namespace edm;
39 using namespace reco;
40 using namespace std;
41 
42 namespace {
43  bool is_B (const reco::Jet& fJet) {return fabs (fJet.eta()) < 1.3;}
44  bool is_E (const reco::Jet& fJet) {return fabs (fJet.eta()) >= 1.3 && fabs (fJet.eta()) < 3.;}
45  bool is_F (const reco::Jet& fJet) {return fabs (fJet.eta()) >= 3.;}
46 }
47 
48 CaloJetTester::CaloJetTester(const edm::ParameterSet& iConfig)
49  : mInputCollection (iConfig.getParameter<edm::InputTag>( "src" )),
50  mInputGenCollection (iConfig.getParameter<edm::InputTag>( "srcGen" )),
51  mOutputFile (iConfig.getUntrackedParameter<std::string>("outputFile", "")),
52  mMatchGenPtThreshold (iConfig.getParameter<double>("genPtThreshold")),
53  mGenEnergyFractionThreshold (iConfig.getParameter<double>("genEnergyFractionThreshold")),
54  mReverseEnergyFractionThreshold (iConfig.getParameter<double>("reverseEnergyFractionThreshold")),
55  mRThreshold (iConfig.getParameter<double>("RThreshold")),
56  mTurnOnEverything (iConfig.getUntrackedParameter<std::string>("TurnOnEverything",""))
57 {
58  numberofevents
59  = mEta = mEtaFineBin = mPhi = mPhiFineBin = mE = mE_80 = mE_3000
60  = mP = mP_80 = mP_3000 = mPt = mPt_80 = mPt_3000
61  = mMass = mMass_80 = mMass_3000 = mConstituents = mConstituents_80
62  = mEtaFirst = mPhiFirst = mEFirst = mEFirst_80 = mEFirst_3000 = mPtFirst = mPtFirst_80 = mPtFirst_3000
63  = mMjj = mMjj_3000 = mDelEta = mDelPhi = mDelPt
64  = mMaxEInEmTowers = mMaxEInHadTowers
65  = mHadEnergyInHO = mHadEnergyInHB = mHadEnergyInHF = mHadEnergyInHE
66  = mHadEnergyInHO_80 = mHadEnergyInHB_80 = mHadEnergyInHE_80
67  = mHadEnergyInHO_3000 = mHadEnergyInHB_3000 = mHadEnergyInHE_3000
68  = mEmEnergyInEB = mEmEnergyInEE = mEmEnergyInHF
69  = mEmEnergyInEB_80 = mEmEnergyInEE_80
70  = mEmEnergyInEB_3000 = mEmEnergyInEE_3000
71  = mEnergyFractionHadronic = mEnergyFractionEm
72  = mHFLong = mHFTotal = mHFLong_80 = mHFLong_3000 = mHFShort = mHFShort_80 = mHFShort_3000
73  = mN90
74  = mCaloMEx = mCaloMEx_3000 = mCaloMEy = mCaloMEy_3000 = mCaloMETSig = mCaloMETSig_3000
75  = mCaloMET = mCaloMET_3000 = mCaloMETPhi = mCaloSumET = mCaloSumET_3000
76  = mHadTiming = mEmTiming
77  = mNJetsEtaC = mNJetsEtaF = mNJets1 = mNJets2
78  = mAllGenJetsPt = mMatchedGenJetsPt = mAllGenJetsEta = mMatchedGenJetsEta
79  = mGenJetMatchEnergyFraction = mReverseMatchEnergyFraction = mRMatch
80  = mDeltaEta = mDeltaPhi = mEScale = mlinEScale = mDeltaE
81  = mHadEnergyProfile = mEmEnergyProfile = mJetEnergyProfile = mHadJetEnergyProfile = mEMJetEnergyProfile
82  = mEScale_pt10 = mEScaleFineBin
83  = mpTScaleB_s = mpTScaleE_s = mpTScaleF_s
84  = mpTScaleB_d = mpTScaleE_d = mpTScaleF_d
85  = mpTScale_60_120_s = mpTScale_200_300_s = mpTScale_600_900_s = mpTScale_2700_3500_s
86  = mpTScale_60_120_d = mpTScale_200_300_d = mpTScale_600_900_d = mpTScale_2700_3500_d
87  = mpTScale1DB_60_120 = mpTScale1DE_60_120 = mpTScale1DF_60_120
88  = mpTScale1DB_200_300 = mpTScale1DE_200_300 = mpTScale1DF_200_300
89  = mpTScale1DB_600_900 = mpTScale1DE_600_900 = mpTScale1DF_600_900
90  = mpTScale1DB_2700_3500 = mpTScale1DE_2700_3500 = mpTScale1DF_2700_3500
91  = mpTScale1D_60_120 = mpTScale1D_200_300 = mpTScale1D_600_900 = mpTScale1D_2700_3500
92  = mHBEne = mHBTime = mHEEne = mHETime = mHFEne = mHFTime = mHOEne = mHOTime
93  = mEBEne = mEBTime = mEEEne = mEETime
94  = mPthat_80 = mPthat_3000
95  = 0;
96 
97  DQMStore* dbe = &*edm::Service<DQMStore>();
98  if (dbe) {
99  dbe->setCurrentFolder("JetMET/RecoJetsV/CaloJetTask_" + mInputCollection.label());
100  //
101  numberofevents = dbe->book1D("numberofevents","numberofevents", 3, 0 , 2);
102  //
103  mEta = dbe->book1D("Eta", "Eta", 100, -5, 5);
104  mEtaFineBin = dbe->book1D("EtaFineBin_Pt10", "EtaFineBin_Pt10", 500, -5, 5);
105  mEtaFineBin1p = dbe->book1D("EtaFineBin1p_Pt10", "EtaFineBin1p_Pt10", 100, 0, 1.3);
106  mEtaFineBin2p = dbe->book1D("EtaFineBin2p_Pt10", "EtaFineBin2p_Pt10", 100, 1.3, 3);
107  mEtaFineBin3p = dbe->book1D("EtaFineBin3p_Pt10", "EtaFineBin3p_Pt10", 100, 3, 5);
108  mEtaFineBin1m = dbe->book1D("EtaFineBin1m_Pt10", "EtaFineBin1m_Pt10", 100, -1.3, 0);
109  mEtaFineBin2m = dbe->book1D("EtaFineBin2m_Pt10", "EtaFineBin2m_Pt10", 100, -3, -1.3);
110  mEtaFineBin3m = dbe->book1D("EtaFineBin3m_Pt10", "EtaFineBin3m_Pt10", 100, -5, -3);
111  //
112  mPhi = dbe->book1D("Phi", "Phi", 70, -3.5, 3.5);
113  mPhiFineBin = dbe->book1D("PhiFineBin_Pt10", "PhiFineBin_Pt10", 350, -3.5, 3.5);
114  //
115  mE = dbe->book1D("E", "E", 100, 0, 500);
116  mE_80 = dbe->book1D("E_80", "E_80", 100, 0, 4500);
117  mE_3000 = dbe->book1D("E_3000", "E_3000", 100, 0, 6000);
118  //
119  mP = dbe->book1D("P", "P", 100, 0, 500);
120  mP_80 = dbe->book1D("P_80", "P_80", 100, 0, 4500);
121  mP_3000 = dbe->book1D("P_3000", "P_3000", 100, 0, 6000);
122  //
123  mPt = dbe->book1D("Pt", "Pt", 100, 0, 50);
124  mPt_80 = dbe->book1D("Pt_80", "Pt_80", 100, 0, 140);
125  mPt_3000 = dbe->book1D("Pt_3000", "Pt_3000", 100, 0, 4000);
126  //
127  mMass = dbe->book1D("Mass", "Mass", 100, 0, 25);
128  mMass_80 = dbe->book1D("Mass_80", "Mass_80", 100, 0, 120);
129  mMass_3000 = dbe->book1D("Mass_3000", "Mass_3000", 100, 0, 500);
130  //
131  mConstituents = dbe->book1D("Constituents", "# of Constituents", 100, 0, 100);
132  mConstituents_80 = dbe->book1D("Constituents_80", "# of Constituents_80", 40, 0, 40);
133  //
134  mEtaFirst = dbe->book1D("EtaFirst", "EtaFirst", 100, -5, 5);
135  mPhiFirst = dbe->book1D("PhiFirst", "PhiFirst", 70, -3.5, 3.5);
136  mEFirst = dbe->book1D("EFirst", "EFirst", 100, 0, 1000);
137  mEFirst_80 = dbe->book1D("EFirst_80", "EFirst_80", 100, 0, 180);
138  mEFirst_3000 = dbe->book1D("EFirst_3000", "EFirst_3000", 100, 0, 4000);
139  mPtFirst = dbe->book1D("PtFirst", "PtFirst", 100, 0, 50);
140  mPtFirst_80 = dbe->book1D("PtFirst_80", "PtFirst_80", 100, 0, 140);
141  mPtFirst_3000 = dbe->book1D("PtFirst_3000", "PtFirst_3000", 100, 0, 4000);
142  //
143  mMjj = dbe->book1D("Mjj", "Mjj", 100, 0, 2000);
144  mMjj_3000 = dbe->book1D("Mjj_3000", "Mjj_3000", 100, 0, 10000);
145  mDelEta = dbe->book1D("DelEta", "DelEta", 100, -.5, .5);
146  mDelPhi = dbe->book1D("DelPhi", "DelPhi", 100, -.5, .5);
147  mDelPt = dbe->book1D("DelPt", "DelPt", 100, -1, 1);
148  //
149  mMaxEInEmTowers = dbe->book1D("MaxEInEmTowers", "MaxEInEmTowers", 100, 0, 100);
150  mMaxEInHadTowers = dbe->book1D("MaxEInHadTowers", "MaxEInHadTowers", 100, 0, 100);
151  mHadEnergyInHO = dbe->book1D("HadEnergyInHO", "HadEnergyInHO", 100, 0, 10);
152  mHadEnergyInHB = dbe->book1D("HadEnergyInHB", "HadEnergyInHB", 100, 0, 50);
153  mHadEnergyInHF = dbe->book1D("HadEnergyInHF", "HadEnergyInHF", 100, 0, 50);
154  mHadEnergyInHE = dbe->book1D("HadEnergyInHE", "HadEnergyInHE", 100, 0, 100);
155  //
156  mHadEnergyInHO_80 = dbe->book1D("HadEnergyInHO_80", "HadEnergyInHO_80", 100, 0, 50);
157  mHadEnergyInHB_80 = dbe->book1D("HadEnergyInHB_80", "HadEnergyInHB_80", 100, 0, 200);
158  mHadEnergyInHE_80 = dbe->book1D("HadEnergyInHE_80", "HadEnergyInHE_80", 100, 0, 1000);
159  mHadEnergyInHO_3000 = dbe->book1D("HadEnergyInHO_3000", "HadEnergyInHO_3000", 100, 0, 500);
160  mHadEnergyInHB_3000 = dbe->book1D("HadEnergyInHB_3000", "HadEnergyInHB_3000", 100, 0, 3000);
161  mHadEnergyInHE_3000 = dbe->book1D("HadEnergyInHE_3000", "HadEnergyInHE_3000", 100, 0, 2000);
162  //
163  mEmEnergyInEB = dbe->book1D("EmEnergyInEB", "EmEnergyInEB", 100, 0, 50);
164  mEmEnergyInEE = dbe->book1D("EmEnergyInEE", "EmEnergyInEE", 100, 0, 50);
165  mEmEnergyInHF = dbe->book1D("EmEnergyInHF", "EmEnergyInHF", 120, -20, 100);
166  mEmEnergyInEB_80 = dbe->book1D("EmEnergyInEB_80", "EmEnergyInEB_80", 100, 0, 200);
167  mEmEnergyInEE_80 = dbe->book1D("EmEnergyInEE_80", "EmEnergyInEE_80", 100, 0, 1000);
168  mEmEnergyInEB_3000= dbe->book1D("EmEnergyInEB_3000", "EmEnergyInEB_3000", 100, 0, 3000);
169  mEmEnergyInEE_3000= dbe->book1D("EmEnergyInEE_3000", "EmEnergyInEE_3000", 100, 0, 2000);
170  mEnergyFractionHadronic = dbe->book1D("EnergyFractionHadronic", "EnergyFractionHadronic", 120, -0.1, 1.1);
171  mEnergyFractionEm = dbe->book1D("EnergyFractionEm", "EnergyFractionEm", 120, -0.1, 1.1);
172  //
173  mHFTotal = dbe->book1D("HFTotal", "HFTotal", 100, 0, 500);
174  mHFTotal_80 = dbe->book1D("HFTotal_80", "HFTotal_80", 100, 0, 3000);
175  mHFTotal_3000 = dbe->book1D("HFTotal_3000", "HFTotal_3000", 100, 0, 6000);
176  mHFLong = dbe->book1D("HFLong", "HFLong", 100, 0, 500);
177  mHFLong_80 = dbe->book1D("HFLong_80", "HFLong_80", 100, 0, 200);
178  mHFLong_3000 = dbe->book1D("HFLong_3000", "HFLong_3000", 100, 0, 1500);
179  mHFShort = dbe->book1D("HFShort", "HFShort", 100, 0, 500);
180  mHFShort_80 = dbe->book1D("HFShort_80", "HFShort_80", 100, 0, 200);
181  mHFShort_3000 = dbe->book1D("HFShort_3000", "HFShort_3000", 100, 0, 1500);
182  //
183  mN90 = dbe->book1D("N90", "N90", 50, 0, 50);
184  //
185  mGenEta = dbe->book1D("GenEta", "GenEta", 100, -5, 5);
186  mGenPhi = dbe->book1D("GenPhi", "GenPhi", 70, -3.5, 3.5);
187  mGenPt = dbe->book1D("GenPt", "GenPt", 100, 0, 50);
188  mGenPt_80 = dbe->book1D("GenPt_80", "GenPt_80", 100, 0, 140);
189  mGenPt_3000 = dbe->book1D("GenPt_3000", "GenPt_3000", 100, 0, 1500);
190  //
191  mGenEtaFirst = dbe->book1D("GenEtaFirst", "GenEtaFirst", 100, -5, 5);
192  mGenPhiFirst = dbe->book1D("GenPhiFirst", "GenPhiFirst", 70, -3.5, 3.5);
193  //
194  mCaloMEx = dbe->book1D("CaloMEx","CaloMEx",200,-150,150);
195  mCaloMEx_3000 = dbe->book1D("CaloMEx_3000","CaloMEx_3000",100,-500,500);
196  mCaloMEy = dbe->book1D("CaloMEy","CaloMEy",200,-150,150);
197  mCaloMEy_3000 = dbe->book1D("CaloMEy_3000","CaloMEy_3000",100,-500,500);
198  mCaloMETSig = dbe->book1D("CaloMETSig","CaloMETSig",100,0,15);
199  mCaloMETSig_3000 = dbe->book1D("CaloMETSig_3000","CaloMETSig_3000",100,0,50);
200  mCaloMET = dbe->book1D("CaloMET","CaloMET",100,0,150);
201  mCaloMET_3000 = dbe->book1D("CaloMET_3000","CaloMET_3000",100,0,1000);
202  mCaloMETPhi = dbe->book1D("CaloMETPhi","CaloMETPhi",70, -3.5, 3.5);
203  mCaloSumET = dbe->book1D("CaloSumET","CaloSumET",100,0,500);
204  mCaloSumET_3000 = dbe->book1D("CaloSumET_3000","CaloSumET_3000",100,3000,8000);
205  //
206  mHadTiming = dbe->book1D("HadTiming", "HadTiming", 75, -50, 100);
207  mEmTiming = dbe->book1D("EMTiming", "EMTiming", 75, -50, 100);
208  //
209  mNJetsEtaC = dbe->book1D("NJetsEtaC_Pt10", "NJetsEtaC_Pt10", 15, 0, 15);
210  mNJetsEtaF = dbe->book1D("NJetsEtaF_Pt10", "NJetsEtaF_Pt10", 15, 0, 15);
211  //
212  mNJets1 = dbe->bookProfile("NJets1", "NJets1", 100, 0, 200, 100, 0, 50, "s");
213  mNJets2 = dbe->bookProfile("NJets2", "NJets2", 100, 0, 4000, 100, 0, 50, "s");
214  //
215  mHBEne = dbe->book1D( "HBEne", "HBEne", 1000, -20, 100 );
216  mHBTime = dbe->book1D( "HBTime", "HBTime", 200, -200, 200 );
217  mHEEne = dbe->book1D( "HEEne", "HEEne", 1000, -20, 100 );
218  mHETime = dbe->book1D( "HETime", "HETime", 200, -200, 200 );
219  mHOEne = dbe->book1D( "HOEne", "HOEne", 1000, -20, 100 );
220  mHOTime = dbe->book1D( "HOTime", "HOTime", 200, -200, 200 );
221  mHFEne = dbe->book1D( "HFEne", "HFEne", 1000, -20, 100 );
222  mHFTime = dbe->book1D( "HFTime", "HFTime", 200, -200, 200 );
223  mEBEne = dbe->book1D( "EBEne", "EBEne", 1000, -20, 100 );
224  mEBTime = dbe->book1D( "EBTime", "EBTime", 200, -200, 200 );
225  mEEEne = dbe->book1D( "EEEne", "EEEne", 1000, -20, 100 );
226  mEETime = dbe->book1D( "EETime", "EETime", 200, -200, 200 );
227  //
228  mPthat_80 = dbe->book1D("Pthat_80", "Pthat_80", 100, 40.0, 160.0);
229  mPthat_3000 = dbe->book1D("Pthat_3000", "Pthat_3000", 100, 2500.0, 4000.0);
230  //
231  double log10PtMin = 0.5; //=3.1622766
232  double log10PtMax = 3.75; //=5623.41325
233  int log10PtBins = 26;
234  double etaMin = -5.;
235  double etaMax = 5.;
236  int etaBins = 50;
237 
238  double linPtMin = 5;
239  double linPtMax = 155;
240  int linPtBins = 15;
241 
242  int log10PtFineBins = 50;
243 
244  mAllGenJetsPt = dbe->book1D("GenJetLOGpT", "GenJet LOG(pT_gen)",
245  log10PtBins, log10PtMin, log10PtMax);
246  mMatchedGenJetsPt = dbe->book1D("MatchedGenJetLOGpT", "MatchedGenJet LOG(pT_gen)",
247  log10PtBins, log10PtMin, log10PtMax);
248  mAllGenJetsEta = dbe->book2D("GenJetEta", "GenJet Eta vs LOG(pT_gen)",
249  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax);
250  mMatchedGenJetsEta = dbe->book2D("MatchedGenJetEta", "MatchedGenJet Eta vs LOG(pT_gen)",
251  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax);
252  //
253  if (mTurnOnEverything.compare("yes")==0) {
254  mHadEnergyProfile = dbe->bookProfile2D("HadEnergyProfile", "HadEnergyProfile", 82, -41, 41, 73, 0, 73, 100, 0, 10000, "s");
255  mEmEnergyProfile = dbe->bookProfile2D("EmEnergyProfile", "EmEnergyProfile", 82, -41, 41, 73, 0, 73, 100, 0, 10000, "s");
256  }
257  mJetEnergyProfile = dbe->bookProfile2D("JetEnergyProfile", "JetEnergyProfile", 50, -5, 5, 36, -3.1415987, 3.1415987, 100, 0, 10000, "s");
258  mHadJetEnergyProfile = dbe->bookProfile2D("HadJetEnergyProfile", "HadJetEnergyProfile", 50, -5, 5, 36, -3.1415987, 3.1415987, 100, 0, 10000, "s");
259  mEMJetEnergyProfile = dbe->bookProfile2D("EMJetEnergyProfile", "EMJetEnergyProfile", 50, -5, 5, 36, -3.1415987, 3.1415987, 100, 0, 10000, "s");
260  //
261  if (mTurnOnEverything.compare("yes")==0) {
262  mGenJetMatchEnergyFraction = dbe->book3D("GenJetMatchEnergyFraction", "GenJetMatchEnergyFraction vs LOG(pT_gen) vs eta",
263  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 101, 0, 1.01);
264  mReverseMatchEnergyFraction = dbe->book3D("ReverseMatchEnergyFraction", "ReverseMatchEnergyFraction vs LOG(pT_gen) vs eta",
265  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 101, 0, 1.01);
266  mRMatch = dbe->book3D("RMatch", "delta(R)(Gen-Calo) vs LOG(pT_gen) vs eta",
267  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 60, 0, 3);
268  mDeltaEta = dbe->book3D("DeltaEta", "DeltaEta vs LOG(pT_gen) vs eta",
269  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, -1, 1);
270  mDeltaPhi = dbe->book3D("DeltaPhi", "DeltaPhi vs LOG(pT_gen) vs eta",
271  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, -1, 1);
272  mEScale = dbe->book3D("EScale", "EnergyScale vs LOG(pT_gen) vs eta",
273  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, 0, 2);
274  mlinEScale = dbe->book3D("linEScale", "EnergyScale vs LOG(pT_gen) vs eta",
275  linPtBins, linPtMin, linPtMax, etaBins, etaMin, etaMax, 100, 0, 2);
276  mDeltaE = dbe->book3D("DeltaE", "DeltaE vs LOG(pT_gen) vs eta",
277  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 2000, -200, 200);
278  //
279  mEScale_pt10 = dbe->book3D("EScale_pt10", "EnergyScale vs LOG(pT_gen) vs eta",
280  log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, 0, 2);
281  mEScaleFineBin = dbe->book3D("EScaleFineBins", "EnergyScale vs LOG(pT_gen) vs eta",
282  log10PtFineBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, 0, 2);
283  }
284 
285  mpTScaleB_s = dbe->bookProfile("pTScaleB_s", "pTScale_s_0<|eta|<1.3",
286  log10PtBins, log10PtMin, log10PtMax, 0, 2, "s");
287  mpTScaleE_s = dbe->bookProfile("pTScaleE_s", "pTScale_s_1.3<|eta|<3.0",
288  log10PtBins, log10PtMin, log10PtMax, 0, 2, "s");
289  mpTScaleF_s = dbe->bookProfile("pTScaleF_s", "pTScale_s_3.0<|eta|<5.0",
290  log10PtBins, log10PtMin, log10PtMax, 0, 2, "s");
291  mpTScaleB_d = dbe->bookProfile("pTScaleB_d", "pTScale_d_0<|eta|<1.3",
292  log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
293  mpTScaleE_d = dbe->bookProfile("pTScaleE_d", "pTScale_d_1.3<|eta|<3.0",
294  log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
295  mpTScaleF_d = dbe->bookProfile("pTScaleF_d", "pTScale_d_3.0<|eta|<5.0",
296  log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
297 
298  mpTScale_60_120_s = dbe->bookProfile("pTScale_60_120_s", "pTScale_s_60<pT<120",
299  etaBins, etaMin, etaMax, 0., 2., "s");
300  mpTScale_200_300_s = dbe->bookProfile("pTScale_200_300_s", "pTScale_s_200<pT<300",
301  etaBins, etaMin, etaMax, 0., 2., "s");
302  mpTScale_600_900_s = dbe->bookProfile("pTScale_600_900_s", "pTScale_s_600<pT<900",
303  etaBins, etaMin, etaMax, 0., 2., "s");
304  mpTScale_2700_3500_s = dbe->bookProfile("pTScale_2700_3500_s", "pTScale_s_2700<pt<3500",
305  etaBins, etaMin, etaMax, 0., 2., "s");
306  mpTScale_60_120_d = dbe->bookProfile("pTScale_60_120_d", "pTScale_d_60<pT<120",
307  etaBins, etaMin, etaMax, 0., 2., " ");
308  mpTScale_200_300_d = dbe->bookProfile("pTScale_200_300_d", "pTScale_d_200<pT<300",
309  etaBins, etaMin, etaMax, 0., 2., " ");
310  mpTScale_600_900_d = dbe->bookProfile("pTScale_600_900_d", "pTScale_d_600<pT<900",
311  etaBins, etaMin, etaMax, 0., 2., " ");
312  mpTScale_2700_3500_d = dbe->bookProfile("pTScale_2700_3500_d", "pTScale_d_2700<pt<3500",
313  etaBins, etaMin, etaMax, 0., 2., " ");
314 
315  mpTScale1DB_60_120 = dbe->book1D("pTScale1DB_60_120", "pTScale_distribution_for_0<|eta|<1.3_60_120",
316  100, 0, 2);
317  mpTScale1DE_60_120 = dbe->book1D("pTScale1DE_60_120", "pTScale_distribution_for_1.3<|eta|<3.0_60_120",
318  100, 0, 2);
319  mpTScale1DF_60_120 = dbe->book1D("pTScale1DF_60_120", "pTScale_distribution_for_3.0<|eta|<5.0_60_120",
320  100, 0, 2);
321 
322  mpTScale1DB_200_300 = dbe->book1D("pTScale1DB_200_300", "pTScale_distribution_for_0<|eta|<1.3_200_300",
323  100, 0, 2);
324  mpTScale1DE_200_300 = dbe->book1D("pTScale1DE_200_300", "pTScale_distribution_for_1.3<|eta|<3.0_200_300",
325  100, 0, 2);
326  mpTScale1DF_200_300 = dbe->book1D("pTScale1DF_200_300", "pTScale_distribution_for_3.0<|eta|<5.0_200_300",
327  100, 0, 2);
328 
329  mpTScale1DB_600_900 = dbe->book1D("pTScale1DB_600_900", "pTScale_distribution_for_0<|eta|<1.3_600_900",
330  100, 0, 2);
331  mpTScale1DE_600_900 = dbe->book1D("pTScale1DE_600_900", "pTScale_distribution_for_1.3<|eta|<3.0_600_900",
332  100, 0, 2);
333  mpTScale1DF_600_900 = dbe->book1D("pTScale1DF_600_900", "pTScale_distribution_for_3.0<|eta|<5.0_600_900",
334  100, 0, 2);
335 
336  mpTScale1DB_2700_3500 = dbe->book1D("pTScale1DB_2700_3500", "pTScale_distribution_for_0<|eta|<1.3_2700_3500",
337  100, 0, 2);
338  mpTScale1DE_2700_3500 = dbe->book1D("pTScale1DE_2700_3500", "pTScale_distribution_for_1.3<|eta|<3.0_2700_3500",
339  100, 0, 2);
340  mpTScale1DF_2700_3500 = dbe->book1D("pTScale1DF_2700_3500", "pTScale_distribution_for_3.0<|eta|<5.0_2700_3500",
341  100, 0, 2);
342 
343  mpTScale1D_60_120 = dbe->book1D("pTScale1D_60_120", "pTScale_distribution_for_60<pT<120",
344  100, 0, 2);
345  mpTScale1D_200_300 = dbe->book1D("pTScale1D_200_300", "pTScale_distribution_for_200<pT<300",
346  100, 0, 2);
347  mpTScale1D_600_900 = dbe->book1D("pTScale1D_600_900", "pTScale_distribution_for_600<pT<900",
348  100, 0, 2);
349  mpTScale1D_2700_3500 = dbe->book1D("pTScale1D_2700_3500", "pTScale_distribution_for_2700<pt<3500",
350  100, 0, 2);
351 
352  }
353 
354  if (mOutputFile.empty ()) {
355  LogInfo("OutputInfo") << " CaloJet histograms will NOT be saved";
356  }
357  else {
358  LogInfo("OutputInfo") << " CaloJethistograms will be saved to file:" << mOutputFile;
359  }
360 }
361 
362 CaloJetTester::~CaloJetTester()
363 {
364 }
365 
366 void CaloJetTester::beginJob(){
367 }
368 
369 void CaloJetTester::endJob() {
370  if (!mOutputFile.empty() && &*edm::Service<DQMStore>()) edm::Service<DQMStore>()->save (mOutputFile);
371 }
372 
373 
374 void CaloJetTester::analyze(const edm::Event& mEvent, const edm::EventSetup& mSetup)
375 {
376  double countsfornumberofevents = 1;
377  numberofevents->Fill(countsfornumberofevents);
378  // *********************************
379  // *** Get pThat
380  // *********************************
381 
382  edm::Handle<HepMCProduct> evt;
383  mEvent.getByLabel("generator", evt);
384  if (evt.isValid()) {
385  HepMC::GenEvent * myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
386 
387  double pthat = myGenEvent->event_scale();
388 
389  mPthat_80->Fill(pthat);
390  mPthat_3000->Fill(pthat);
391 
392  delete myGenEvent;
393  }
394 
395  // ***********************************
396  // *** Get CaloMET
397  // ***********************************
398 
399  const CaloMET *calomet;
400  edm::Handle<CaloMETCollection> calo;
401  mEvent.getByLabel("met", calo);
402  if (!calo.isValid()) {
403  edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task";
404  edm::LogInfo("OutputInfo") << " MET Task cannot continue...!";
405  } else {
406  const CaloMETCollection *calometcol = calo.product();
407  calomet = &(calometcol->front());
408 
409  double caloSumET = calomet->sumEt();
410  double caloMETSig = calomet->mEtSig();
411  double caloMET = calomet->pt();
412  double caloMEx = calomet->px();
413  double caloMEy = calomet->py();
414  double caloMETPhi = calomet->phi();
415 
416  mCaloMEx->Fill(caloMEx);
417  mCaloMEx_3000->Fill(caloMEx);
418  mCaloMEy->Fill(caloMEy);
419  mCaloMEy_3000->Fill(caloMEy);
420  mCaloMET->Fill(caloMET);
421  mCaloMET_3000->Fill(caloMET);
422  mCaloMETPhi->Fill(caloMETPhi);
423  mCaloSumET->Fill(caloSumET);
424  mCaloSumET_3000->Fill(caloSumET);
425  mCaloMETSig->Fill(caloMETSig);
426  mCaloMETSig_3000->Fill(caloMETSig);
427  }
428 
429  // ***********************************
430  // *** Get the CaloTower collection
431  // ***********************************
432  Handle<CaloTowerCollection> caloTowers;
433  mEvent.getByLabel( "towerMaker", caloTowers );
434  if (caloTowers.isValid()) {
435  for( CaloTowerCollection::const_iterator cal = caloTowers->begin(); cal != caloTowers->end(); ++ cal ){
436 
437  //To compensate for the index
438  if (mTurnOnEverything.compare("yes")==0) {
439  if (cal->ieta() >> 0 ){mHadEnergyProfile->Fill (cal->ieta()-1, cal->iphi(), cal->hadEnergy());
440  mEmEnergyProfile->Fill (cal->ieta()-1, cal->iphi(), cal->emEnergy());}
441  mHadEnergyProfile->Fill (cal->ieta(), cal->iphi(), cal->hadEnergy());
442  mEmEnergyProfile->Fill (cal->ieta(), cal->iphi(), cal->emEnergy());
443  }
444 
445  mHadTiming->Fill (cal->hcalTime());
446  mEmTiming->Fill (cal->ecalTime());
447  }
448  }
449 
450  // ***********************************
451  // *** Get the RecHits collection
452  // ***********************************
453  try {
454  std::vector<edm::Handle<HBHERecHitCollection> > colls;
455  mEvent.getManyByType(colls);
456  std::vector<edm::Handle<HBHERecHitCollection> >::iterator i;
457  for (i=colls.begin(); i!=colls.end(); i++) {
458  for (HBHERecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
459  // std::cout << *j << std::endl;
460  if (j->id().subdet() == HcalBarrel) {
461  mHBEne->Fill(j->energy());
462  mHBTime->Fill(j->time());
463  }
464  if (j->id().subdet() == HcalEndcap) {
465  mHEEne->Fill(j->energy());
466  mHETime->Fill(j->time());
467  }
468 
469  }
470  }
471  } catch (...) {
472  edm::LogInfo("OutputInfo") << " No HB/HE RecHits.";
473  }
474 
475  try {
476  std::vector<edm::Handle<HFRecHitCollection> > colls;
477  mEvent.getManyByType(colls);
478  std::vector<edm::Handle<HFRecHitCollection> >::iterator i;
479  for (i=colls.begin(); i!=colls.end(); i++) {
480  for (HFRecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
481  // std::cout << *j << std::endl;
482  if (j->id().subdet() == HcalForward) {
483  mHFEne->Fill(j->energy());
484  mHFTime->Fill(j->time());
485  }
486  }
487  }
488  } catch (...) {
489  edm::LogInfo("OutputInfo") << " No HF RecHits.";
490  }
491 
492  try {
493  std::vector<edm::Handle<HORecHitCollection> > colls;
494  mEvent.getManyByType(colls);
495  std::vector<edm::Handle<HORecHitCollection> >::iterator i;
496  for (i=colls.begin(); i!=colls.end(); i++) {
497  for (HORecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
498  if (j->id().subdet() == HcalOuter) {
499  mHOEne->Fill(j->energy());
500  mHOTime->Fill(j->time());
501  }
502  }
503  }
504  } catch (...) {
505  edm::LogInfo("OutputInfo") << " No HO RecHits.";
506  }
507  try {
508  std::vector<edm::Handle<EBRecHitCollection> > colls;
509  mEvent.getManyByType(colls);
510  std::vector<edm::Handle<EBRecHitCollection> >::iterator i;
511  for (i=colls.begin(); i!=colls.end(); i++) {
512  for (EBRecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
513  // if (j->id() == EcalBarrel) {
514  mEBEne->Fill(j->energy());
515  mEBTime->Fill(j->time());
516  // }
517  // std::cout << *j << std::endl;
518  // std::cout << j->id() << std::endl;
519  }
520  }
521  } catch (...) {
522  edm::LogInfo("OutputInfo") << " No EB RecHits.";
523  }
524 
525  try {
526  std::vector<edm::Handle<EERecHitCollection> > colls;
527  mEvent.getManyByType(colls);
528  std::vector<edm::Handle<EERecHitCollection> >::iterator i;
529  for (i=colls.begin(); i!=colls.end(); i++) {
530  for (EERecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
531  // if (j->id().subdet() == EcalEndcap) {
532  mEEEne->Fill(j->energy());
533  mEETime->Fill(j->time());
534  // }
535  // std::cout << *j << std::endl;
536  }
537  }
538  } catch (...) {
539  edm::LogInfo("OutputInfo") << " No EE RecHits.";
540  }
541 
542 
543  //***********************************
544  //*** Get the Jet collection
545  //***********************************
546  math::XYZTLorentzVector p4tmp[2];
547  Handle<CaloJetCollection> caloJets;
548  mEvent.getByLabel(mInputCollection, caloJets);
549  if (!caloJets.isValid()) return;
550  CaloJetCollection::const_iterator jet = caloJets->begin ();
551  int jetIndex = 0;
552  int nJet = 0;
553  int nJetF = 0;
554  int nJetC = 0;
555  for (; jet != caloJets->end (); jet++, jetIndex++) {
556  if (mEta) mEta->Fill (jet->eta());
557 
558  if (jet->pt() > 10.) {
559  if (fabs(jet->eta()) > 1.3)
560  nJetF++;
561  else
562  nJetC++;
563  }
564  if (jet->pt() > 10.) {
565  if (mEtaFineBin) mEtaFineBin->Fill (jet->eta());
566  if (mEtaFineBin1p) mEtaFineBin1p->Fill (jet->eta());
567  if (mEtaFineBin2p) mEtaFineBin2p->Fill (jet->eta());
568  if (mEtaFineBin3p) mEtaFineBin3p->Fill (jet->eta());
569  if (mEtaFineBin1m) mEtaFineBin1m->Fill (jet->eta());
570  if (mEtaFineBin2m) mEtaFineBin2m->Fill (jet->eta());
571  if (mEtaFineBin3m) mEtaFineBin3m->Fill (jet->eta());
572  if (mPhiFineBin) mPhiFineBin->Fill (jet->phi());
573  }
574  if (mPhi) mPhi->Fill (jet->phi());
575  if (mE) mE->Fill (jet->energy());
576  if (mE_80) mE_80->Fill (jet->energy());
577  if (mE_3000) mE_3000->Fill (jet->energy());
578  if (mP) mP->Fill (jet->p());
579  if (mP_80) mP_80->Fill (jet->p());
580  if (mP_3000) mP_3000->Fill (jet->p());
581  if (mPt) mPt->Fill (jet->pt());
582  if (mPt_80) mPt_80->Fill (jet->pt());
583  if (mPt_3000) mPt_3000->Fill (jet->pt());
584  if (mMass) mMass->Fill (jet->mass());
585  if (mMass_80) mMass_80->Fill (jet->mass());
586  if (mMass_3000) mMass_3000->Fill (jet->mass());
587  if (mConstituents) mConstituents->Fill (jet->nConstituents());
588  if (mConstituents_80) mConstituents_80->Fill (jet->nConstituents());
589  if (jet == caloJets->begin ()) { // first jet
590  if (mEtaFirst) mEtaFirst->Fill (jet->eta());
591  if (mPhiFirst) mPhiFirst->Fill (jet->phi());
592  if (mEFirst) mEFirst->Fill (jet->energy());
593  if (mEFirst_80) mEFirst_80->Fill (jet->energy());
594  if (mEFirst_3000) mEFirst_3000->Fill (jet->energy());
595  if (mPtFirst) mPtFirst->Fill (jet->pt());
596  if (mPtFirst_80) mPtFirst_80->Fill (jet->pt());
597  if (mPtFirst_3000) mPtFirst_3000->Fill (jet->pt());
598  }
599  if (jetIndex == 0) {
600  nJet++;
601  p4tmp[0] = jet->p4();
602  }
603  if (jetIndex == 1) {
604  nJet++;
605  p4tmp[1] = jet->p4();
606  }
607 
608  if (mMaxEInEmTowers) mMaxEInEmTowers->Fill (jet->maxEInEmTowers());
609  if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->maxEInHadTowers());
610  if (mHadEnergyInHO) mHadEnergyInHO->Fill (jet->hadEnergyInHO());
611  if (mHadEnergyInHO_80) mHadEnergyInHO_80->Fill (jet->hadEnergyInHO());
612  if (mHadEnergyInHO_3000) mHadEnergyInHO_3000->Fill (jet->hadEnergyInHO());
613  if (mHadEnergyInHB) mHadEnergyInHB->Fill (jet->hadEnergyInHB());
614  if (mHadEnergyInHB_80) mHadEnergyInHB_80->Fill (jet->hadEnergyInHB());
615  if (mHadEnergyInHB_3000) mHadEnergyInHB_3000->Fill (jet->hadEnergyInHB());
616  if (mHadEnergyInHF) mHadEnergyInHF->Fill (jet->hadEnergyInHF());
617  if (mHadEnergyInHE) mHadEnergyInHE->Fill (jet->hadEnergyInHE());
618  if (mHadEnergyInHE_80) mHadEnergyInHE_80->Fill (jet->hadEnergyInHE());
619  if (mHadEnergyInHE_3000) mHadEnergyInHE_3000->Fill (jet->hadEnergyInHE());
620  if (mEmEnergyInEB) mEmEnergyInEB->Fill (jet->emEnergyInEB());
621  if (mEmEnergyInEB_80) mEmEnergyInEB_80->Fill (jet->emEnergyInEB());
622  if (mEmEnergyInEB_3000) mEmEnergyInEB_3000->Fill (jet->emEnergyInEB());
623  if (mEmEnergyInEE) mEmEnergyInEE->Fill (jet->emEnergyInEE());
624  if (mEmEnergyInEE_80) mEmEnergyInEE_80->Fill (jet->emEnergyInEE());
625  if (mEmEnergyInEE_3000) mEmEnergyInEE_3000->Fill (jet->emEnergyInEE());
626  if (mEmEnergyInHF) mEmEnergyInHF->Fill (jet->emEnergyInHF());
627  if (mEnergyFractionHadronic) mEnergyFractionHadronic->Fill (jet->energyFractionHadronic());
628  if (mEnergyFractionEm) mEnergyFractionEm->Fill (jet->emEnergyFraction());
629 
630  if (mHFTotal) mHFTotal->Fill (jet->hadEnergyInHF()+jet->emEnergyInHF());
631  if (mHFTotal_80) mHFTotal_80->Fill (jet->hadEnergyInHF()+jet->emEnergyInHF());
632  if (mHFTotal_3000) mHFTotal_3000->Fill (jet->hadEnergyInHF()+jet->emEnergyInHF());
633  if (mHFLong) mHFLong->Fill (jet->hadEnergyInHF()*0.5+jet->emEnergyInHF());
634  if (mHFLong_80) mHFLong_80->Fill (jet->hadEnergyInHF()*0.5+jet->emEnergyInHF());
635  if (mHFLong_3000) mHFLong_3000->Fill (jet->hadEnergyInHF()*0.5+jet->emEnergyInHF());
636  if (mHFShort) mHFShort->Fill (jet->hadEnergyInHF()*0.5);
637  if (mHFShort_80) mHFShort_80->Fill (jet->hadEnergyInHF()*0.5);
638  if (mHFShort_3000) mHFShort_3000->Fill (jet->hadEnergyInHF()*0.5);
639 
640 
641 
642  if (mN90) mN90->Fill (jet->n90());
643  mJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->energy());
644  mHadJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->hadEnergyInHO()+jet->hadEnergyInHB()+jet->hadEnergyInHF()+jet->hadEnergyInHE());
645  mEMJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->emEnergyInEB()+jet->emEnergyInEE()+jet->emEnergyInHF());
646  }
647 
648 
649 
650  if (mNJetsEtaC) mNJetsEtaC->Fill( nJetC );
651  if (mNJetsEtaF) mNJetsEtaF->Fill( nJetF );
652 
653  if (nJet == 2) {
654  if (mMjj) mMjj->Fill( (p4tmp[0]+p4tmp[1]).mass() );
655  if (mMjj_3000) mMjj_3000->Fill( (p4tmp[0]+p4tmp[1]).mass() );
656  }
657 
658  // Count Jets above Pt cut
659  for (int istep = 0; istep < 100; ++istep) {
660  int njet = 0;
661  float ptStep = (istep * (200./100.));
662 
663  for ( CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
664  if ( cal->pt() > ptStep ) njet++;
665  }
666  mNJets1->Fill( ptStep, njet );
667  }
668 
669  for (int istep = 0; istep < 100; ++istep) {
670  int njet = 0;
671  float ptStep = (istep * (4000./100.));
672  for ( CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
673  if ( cal->pt() > ptStep ) njet++;
674  }
675  mNJets2->Fill( ptStep, njet );
676  }
677 
678  // Gen jet analysis
679  Handle<GenJetCollection> genJets;
680  mEvent.getByLabel(mInputGenCollection, genJets);
681  if (!genJets.isValid()) return;
682  GenJetCollection::const_iterator gjet = genJets->begin ();
683  int gjetIndex = 0;
684  for (; gjet != genJets->end (); gjet++, gjetIndex++) {
685  if (mGenEta) mGenEta->Fill (gjet->eta());
686  if (mGenPhi) mGenPhi->Fill (gjet->phi());
687  if (mGenPt) mGenPt->Fill (gjet->pt());
688  if (mGenPt_80) mGenPt_80->Fill (gjet->pt());
689  if (mGenPt_3000) mGenPt_3000->Fill (gjet->pt());
690  if (gjet == genJets->begin ()) { // first jet
691  if (mGenEtaFirst) mGenEtaFirst->Fill (gjet->eta());
692  if (mGenPhiFirst) mGenPhiFirst->Fill (gjet->phi());
693  }
694  }
695 
696 
697  // now match CaloJets to GenJets
698  JetMatchingTools jetMatching (mEvent);
699  if (!(mInputGenCollection.label().empty())) {
700  // Handle<GenJetCollection> genJets;
701  // mEvent.getByLabel(mInputGenCollection, genJets);
702 
703  std::vector <std::vector <const reco::GenParticle*> > genJetConstituents (genJets->size());
704  std::vector <std::vector <const reco::GenParticle*> > caloJetConstituents (caloJets->size());
705  if (mRThreshold > 0) {
706  }
707  else {
708  for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
709  genJetConstituents [iGenJet] = jetMatching.getGenParticles ((*genJets) [iGenJet]);
710  }
711 
712  for (unsigned iCaloJet = 0; iCaloJet < caloJets->size(); ++iCaloJet) {
713  caloJetConstituents [iCaloJet] = jetMatching.getGenParticles ((*caloJets) [iCaloJet], false);
714  }
715  }
716 
717  for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) { //****************************************************************
718  //for (unsigned iGenJet = 0; iGenJet < 1; ++iGenJet) { // only FIRST Jet !!!!
719  const GenJet& genJet = (*genJets) [iGenJet];
720  double genJetPt = genJet.pt();
721 
722  //std::cout << iGenJet <<". Genjet: pT = " << genJetPt << "GeV" << std::endl; // *****************************************************
723 
724  if (fabs(genJet.eta()) > 5.) continue; // out of detector
725  if (genJetPt < mMatchGenPtThreshold) continue; // no low momentum
726  double logPtGen = log10 (genJetPt);
727  mAllGenJetsPt->Fill (logPtGen);
728  mAllGenJetsEta->Fill (logPtGen, genJet.eta());
729  if (caloJets->size() <= 0) continue; // no CaloJets - nothing to match
730  if (mRThreshold > 0) {
731  unsigned iCaloJetBest = 0;
732  double deltaRBest = 999.;
733  for (unsigned iCaloJet = 0; iCaloJet < caloJets->size(); ++iCaloJet) {
734  double dR = deltaR (genJet.eta(), genJet.phi(), (*caloJets) [iCaloJet].eta(), (*caloJets) [iCaloJet].phi());
735  if (deltaRBest < mRThreshold && dR < mRThreshold && genJet.pt() > 5.) {
736  /*
737  std::cout << "Yet another matched jet for GenJet pt=" << genJet.pt()
738  << " previous CaloJet pt/dr: " << (*caloJets) [iCaloJetBest].pt() << '/' << deltaRBest
739  << " new CaloJet pt/dr: " << (*caloJets) [iCaloJet].pt() << '/' << dR
740  << std::endl;
741  */
742  }
743  if (dR < deltaRBest) {
744  iCaloJetBest = iCaloJet;
745  deltaRBest = dR;
746  }
747  }
748  if (mTurnOnEverything.compare("yes")==0) {
749  mRMatch->Fill (logPtGen, genJet.eta(), deltaRBest);
750  }
751  if (deltaRBest < mRThreshold) { // Matched
752  fillMatchHists (genJet, (*caloJets) [iCaloJetBest]);
753  }
754  }
755  else {
756  unsigned iCaloJetBest = 0;
757  double energyFractionBest = 0.;
758  for (unsigned iCaloJet = 0; iCaloJet < caloJets->size(); ++iCaloJet) {
759  double energyFraction = jetMatching.overlapEnergyFraction (genJetConstituents [iGenJet],
760  caloJetConstituents [iCaloJet]);
761  if (energyFraction > energyFractionBest) {
762  iCaloJetBest = iCaloJet;
763  energyFractionBest = energyFraction;
764  }
765  }
766  if (mTurnOnEverything.compare("yes")==0) {
767  mGenJetMatchEnergyFraction->Fill (logPtGen, genJet.eta(), energyFractionBest);
768  }
769  if (energyFractionBest > mGenEnergyFractionThreshold) { // good enough
770  double reverseEnergyFraction = jetMatching.overlapEnergyFraction (caloJetConstituents [iCaloJetBest],
771  genJetConstituents [iGenJet]);
772  if (mTurnOnEverything.compare("yes")==0) {
773  mReverseMatchEnergyFraction->Fill (logPtGen, genJet.eta(), reverseEnergyFraction);
774  }
775  if (reverseEnergyFraction > mReverseEnergyFractionThreshold) { // Matched
776  fillMatchHists (genJet, (*caloJets) [iCaloJetBest]);
777  }
778  }
779  }
780  }
781  }
782 }
783 
784 
785 void CaloJetTester::fillMatchHists (const reco::GenJet& fGenJet, const reco::CaloJet& fCaloJet) {
786  double logPtGen = log10 (fGenJet.pt());
787  double PtGen = fGenJet.pt();
788  double PtCalo = fCaloJet.pt();
789  mMatchedGenJetsPt->Fill (logPtGen);
790  mMatchedGenJetsEta->Fill (logPtGen, fGenJet.eta());
791 
792  double PtThreshold = 10.;
793 
794  if (mTurnOnEverything.compare("yes")==0) {
795  mDeltaEta->Fill (logPtGen, fGenJet.eta(), fCaloJet.eta()-fGenJet.eta());
796  mDeltaPhi->Fill (logPtGen, fGenJet.eta(), fCaloJet.phi()-fGenJet.phi());
797  mEScale->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
798  mlinEScale->Fill (fGenJet.pt(), fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
799  mDeltaE->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()-fGenJet.energy());
800 
801  mEScaleFineBin->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
802 
803  if (fGenJet.pt()>PtThreshold) {
804  mEScale_pt10->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
805 
806  }
807 
808  }
809  if (fCaloJet.pt() > PtThreshold) {
810  mDelEta->Fill (fGenJet.eta()-fCaloJet.eta());
811  mDelPhi->Fill (fGenJet.phi()-fCaloJet.phi());
812  mDelPt->Fill ((fGenJet.pt()-fCaloJet.pt())/fGenJet.pt());
813  }
814 
815  if (fabs(fGenJet.eta())<1.3) {
816 
817  mpTScaleB_s->Fill (log10(PtGen), PtCalo/PtGen);
818  mpTScaleB_d->Fill (log10(PtGen), PtCalo/PtGen);
819  if (PtGen>60.0 && PtGen<120.0) {
820  mpTScale1DB_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
821  }
822  if (PtGen>200.0 && PtGen<300.0) {
823  mpTScale1DB_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
824  }
825  if (PtGen>600.0 && PtGen<900.0) {
826  mpTScale1DB_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
827  }
828  if (PtGen>2700.0 && PtGen<3500.0) {
829  mpTScale1DB_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
830  }
831  }
832 
833  if (fabs(fGenJet.eta())>1.3 && fabs(fGenJet.eta())<3.0) {
834 
835  mpTScaleE_s->Fill (log10(PtGen), PtCalo/PtGen);
836  mpTScaleE_d->Fill (log10(PtGen), PtCalo/PtGen);
837  if (PtGen>60.0 && PtGen<120.0) {
838  mpTScale1DE_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
839  }
840  if (PtGen>200.0 && PtGen<300.0) {
841  mpTScale1DE_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
842  }
843  if (PtGen>600.0 && PtGen<900.0) {
844  mpTScale1DE_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
845  }
846  if (PtGen>2700.0 && PtGen<3500.0) {
847  mpTScale1DE_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
848  }
849  }
850 
851  if (fabs(fGenJet.eta())>3.0 && fabs(fGenJet.eta())<5.0) {
852 
853  mpTScaleF_s->Fill (log10(PtGen), PtCalo/PtGen);
854  mpTScaleF_d->Fill (log10(PtGen), PtCalo/PtGen);
855  if (PtGen>60.0 && PtGen<120.0) {
856  mpTScale1DF_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
857  }
858  if (PtGen>200.0 && PtGen<300.0) {
859  mpTScale1DF_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
860  }
861  if (PtGen>600.0 && PtGen<900.0) {
862  mpTScale1DF_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
863  }
864  if (PtGen>2700.0 && PtGen<3500.0) {
865  mpTScale1DF_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
866  }
867  }
868 
869  if (fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
870  mpTScale_60_120_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
871  mpTScale_60_120_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
872  mpTScale1D_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
873  }
874 
875  if (fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
876  mpTScale_200_300_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
877  mpTScale_200_300_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
878  mpTScale1D_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
879  }
880 
881  if (fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
882  mpTScale_600_900_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
883  mpTScale_600_900_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
884  mpTScale1D_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
885  }
886 
887  if (fGenJet.pt()>2700.0 && fGenJet.pt()<3500.0) {
888  mpTScale_2700_3500_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
889  mpTScale_2700_3500_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
890  mpTScale1D_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
891  }
892 
893 
894 
895 }
CaloJetTester::mEtaFineBin3m
MonitorElement * mEtaFineBin3m
Definition: CaloJetTester.h:54
edm::Event::getManyByType
void getManyByType(std::vector< Handle< PROD > > &results) const
Definition: Event.h:407
CaloJetTester::mpTScaleB_s
MonitorElement * mpTScaleB_s
Definition: CaloJetTester.h:160
i
int i
Definition: DBlmapReader.cc:9
CaloJetTester::mpTScale1DF_60_120
MonitorElement * mpTScale1DF_60_120
Definition: CaloJetTester.h:179
CaloJetTester::mEmEnergyInEE_3000
MonitorElement * mEmEnergyInEE_3000
Definition: CaloJetTester.h:114
CaloJetTester::mMatchedGenJetsPt
MonitorElement * mMatchedGenJetsPt
Definition: CaloJetTester.h:145
reco::CaloJet
Jets made from CaloTowers.
Definition: CaloJet.h:30
CaloJetTester::fillMatchHists
void fillMatchHists(const reco::GenJet &fGenJet, const reco::CaloJet &fCaloJet)
Definition: CaloJetTester.cc:785
CaloJetTester::mEEEne
MonitorElement * mEEEne
Definition: CaloJetTester.h:238
CaloJetTester::mInputCollection
edm::InputTag mInputCollection
Definition: CaloJetTester.h:35
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:514
CaloJetTester::mHadEnergyInHE_80
MonitorElement * mHadEnergyInHE_80
Definition: CaloJetTester.h:104
CaloJetTester::mCaloMEy_3000
MonitorElement * mCaloMEy_3000
Definition: CaloJetTester.h:218
CaloJetTester::mHadEnergyInHB_80
MonitorElement * mHadEnergyInHB_80
Definition: CaloJetTester.h:103
CaloJetTester::mMaxEInHadTowers
MonitorElement * mMaxEInHadTowers
Definition: CaloJetTester.h:97
CaloJetTester::mReverseEnergyFractionThreshold
double mReverseEnergyFractionThreshold
Definition: CaloJetTester.h:201
CaloJetTester::mEtaFineBin2p
MonitorElement * mEtaFineBin2p
Definition: CaloJetTester.h:50
CaloJetTester::mpTScale1DF_200_300
MonitorElement * mpTScale1DF_200_300
Definition: CaloJetTester.h:182
CaloJetTester::mConstituents
MonitorElement * mConstituents
Definition: CaloJetTester.h:69
CaloJetTester::mpTScale1D_60_120
MonitorElement * mpTScale1D_60_120
Definition: CaloJetTester.h:189
CaloJetTester::mPtFirst
MonitorElement * mPtFirst
Definition: CaloJetTester.h:81
DQMStore::book3D
MonitorElement * book3D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ)
Book 3D histogram.
Definition: DQMStore.cc:776
CaloJetTester::mDelEta
MonitorElement * mDelEta
Definition: CaloJetTester.h:194
CaloJetTester::mNJets2
MonitorElement * mNJets2
Definition: CaloJetTester.h:89
CaloJetTester::mCaloMEx
MonitorElement * mCaloMEx
Definition: CaloJetTester.h:215
CaloJetTester::mGenPt_3000
MonitorElement * mGenPt_3000
Definition: CaloJetTester.h:137
CaloJetTester::mpTScale1DE_60_120
MonitorElement * mpTScale1DE_60_120
Definition: CaloJetTester.h:178
CaloJetTester::mEmEnergyInEB_3000
MonitorElement * mEmEnergyInEB_3000
Definition: CaloJetTester.h:113
CaloJetTester::mpTScale_600_900_d
MonitorElement * mpTScale_600_900_d
Definition: CaloJetTester.h:174
CaloJetTester::mHadEnergyInHE
MonitorElement * mHadEnergyInHE
Definition: CaloJetTester.h:101
reco::Jet
Base class for all types of Jets.
Definition: Jet.h:21
CaloJetTester::mRThreshold
double mRThreshold
Definition: CaloJetTester.h:202
CaloJetTester::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: CaloJetTester.cc:374
CaloJetTester::mNJetsEtaF
MonitorElement * mNJetsEtaF
Definition: CaloJetTester.h:86
CaloJetTester::mEtaFineBin
MonitorElement * mEtaFineBin
Definition: CaloJetTester.h:48
CaloJetTester::mpTScaleB_d
MonitorElement * mpTScaleB_d
Definition: CaloJetTester.h:163
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:81
CaloJetTester::mOutputFile
std::string mOutputFile
Definition: CaloJetTester.h:37
CaloJetTester::mpTScale1DE_600_900
MonitorElement * mpTScale1DE_600_900
Definition: CaloJetTester.h:184
CaloJetTester::mEtaFineBin2m
MonitorElement * mEtaFineBin2m
Definition: CaloJetTester.h:53
CaloJetTester::mpTScale_600_900_s
MonitorElement * mpTScale_600_900_s
Definition: CaloJetTester.h:169
CaloJetTester::mHFTotal_80
MonitorElement * mHFTotal_80
Definition: CaloJetTester.h:118
CaloJetTester::mDelPt
MonitorElement * mDelPt
Definition: CaloJetTester.h:196
CaloJetTester::mPtFirst_3000
MonitorElement * mPtFirst_3000
Definition: CaloJetTester.h:83
CaloJetTester::mTurnOnEverything
std::string mTurnOnEverything
Definition: CaloJetTester.h:205
CaloJetTester::mHadEnergyInHB
MonitorElement * mHadEnergyInHB
Definition: CaloJetTester.h:99
CaloJetTester::mHOEne
MonitorElement * mHOEne
Definition: CaloJetTester.h:232
CaloJetTester::mEBEne
MonitorElement * mEBEne
Definition: CaloJetTester.h:236
CaloJetTester::mDelPhi
MonitorElement * mDelPhi
Definition: CaloJetTester.h:195
CaloJetTester::mPt_80
MonitorElement * mPt_80
Definition: CaloJetTester.h:64
CaloJetTester::mHadEnergyInHE_3000
MonitorElement * mHadEnergyInHE_3000
Definition: CaloJetTester.h:107
CaloJetTester::mpTScale1DB_60_120
MonitorElement * mpTScale1DB_60_120
Definition: CaloJetTester.h:177
CaloJetTester::mpTScale_2700_3500_s
MonitorElement * mpTScale_2700_3500_s
Definition: CaloJetTester.h:170
CaloJetTester::mPt
MonitorElement * mPt
Definition: CaloJetTester.h:63
CaloJetTester::mEScale_pt10
MonitorElement * mEScale_pt10
Definition: CaloJetTester.h:157
CaloJetTester::mNJets1
MonitorElement * mNJets1
Definition: CaloJetTester.h:88
CaloJetTester::mHFShort_80
MonitorElement * mHFShort_80
Definition: CaloJetTester.h:124
CaloJetTester::mMatchGenPtThreshold
double mMatchGenPtThreshold
Definition: CaloJetTester.h:199
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:101
CaloJetTester::mEmEnergyInHF
MonitorElement * mEmEnergyInHF
Definition: CaloJetTester.h:110
CaloJetTester::mpTScale1D_200_300
MonitorElement * mpTScale1D_200_300
Definition: CaloJetTester.h:190
CaloJetTester::mHadEnergyInHO_3000
MonitorElement * mHadEnergyInHO_3000
Definition: CaloJetTester.h:105
CaloJetTester::mpTScaleF_s
MonitorElement * mpTScaleF_s
Definition: CaloJetTester.h:162
njet
int njet
Definition: HydjetWrapper.h:91
CaloJetTester::mpTScale_2700_3500_d
MonitorElement * mpTScale_2700_3500_d
Definition: CaloJetTester.h:175
CaloJetTester::mpTScale1DE_2700_3500
MonitorElement * mpTScale1DE_2700_3500
Definition: CaloJetTester.h:187
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:143
CaloJetTester::mHadEnergyInHB_3000
MonitorElement * mHadEnergyInHB_3000
Definition: CaloJetTester.h:106
CaloJetTester::mpTScale1DF_2700_3500
MonitorElement * mpTScale1DF_2700_3500
Definition: CaloJetTester.h:188
CaloJetTester::mEnergyFractionEm
MonitorElement * mEnergyFractionEm
Definition: CaloJetTester.h:116
CaloJetTester::mCaloMEy
MonitorElement * mCaloMEy
Definition: CaloJetTester.h:217
CaloJetTester::mCaloMETSig
MonitorElement * mCaloMETSig
Definition: CaloJetTester.h:224
CaloJetTester::mEmTiming
MonitorElement * mEmTiming
Definition: CaloJetTester.h:73
CaloJetTester::mRMatch
MonitorElement * mRMatch
Definition: CaloJetTester.h:150
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
CaloJetTester::mAllGenJetsPt
MonitorElement * mAllGenJetsPt
Definition: CaloJetTester.h:144
CaloJetTester::mHadTiming
MonitorElement * mHadTiming
Definition: CaloJetTester.h:72
CaloJetTester::mHFTotal_3000
MonitorElement * mHFTotal_3000
Definition: CaloJetTester.h:119
CaloJetTester::mHFShort_3000
MonitorElement * mHFShort_3000
Definition: CaloJetTester.h:125
CaloMETCollection
Collection of Calo MET.
CaloJetTester::mEFirst
MonitorElement * mEFirst
Definition: CaloJetTester.h:78
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:30
CaloJetTester::mCaloMETSig_3000
MonitorElement * mCaloMETSig_3000
Definition: CaloJetTester.h:225
CaloJetTester::mCaloMET
MonitorElement * mCaloMET
Definition: CaloJetTester.h:219
CaloJetTester::mMjj
MonitorElement * mMjj
Definition: CaloJetTester.h:92
reco::LeafCandidate::energy
virtual double energy() const
energy
Definition: LeafCandidate.h:119
reco::MET::mEtSig
double mEtSig() const
Definition: MET.h:50
CaloJetTester::endJob
virtual void endJob()
Definition: CaloJetTester.cc:369
CaloJetTester::mPhiFirst
MonitorElement * mPhiFirst
Definition: CaloJetTester.h:77
CaloJetTester::mpTScale_200_300_s
MonitorElement * mpTScale_200_300_s
Definition: CaloJetTester.h:168
CaloJetTester::mHadJetEnergyProfile
MonitorElement * mHadJetEnergyProfile
Definition: CaloJetTester.h:211
CaloJetTester::mEmEnergyProfile
MonitorElement * mEmEnergyProfile
Definition: CaloJetTester.h:209
reco::MET::sumEt
double sumEt() const
Definition: MET.h:48
CaloJetTester::mDeltaE
MonitorElement * mDeltaE
Definition: CaloJetTester.h:155
CaloJetTester::mHFTime
MonitorElement * mHFTime
Definition: CaloJetTester.h:235
CaloJetTester::mEnergyFractionHadronic
MonitorElement * mEnergyFractionHadronic
Definition: CaloJetTester.h:115
CaloJetTester::mHFLong_3000
MonitorElement * mHFLong_3000
Definition: CaloJetTester.h:122
CaloJetTester::mE_80
MonitorElement * mE_80
Definition: CaloJetTester.h:58
CaloJetTester::mP_80
MonitorElement * mP_80
Definition: CaloJetTester.h:61
CaloJetTester::mGenEtaFirst
MonitorElement * mGenEtaFirst
Definition: CaloJetTester.h:140
CaloJetTester::mGenJetMatchEnergyFraction
MonitorElement * mGenJetMatchEnergyFraction
Definition: CaloJetTester.h:148
CaloJetTester::mMass
MonitorElement * mMass
Definition: CaloJetTester.h:66
CaloJetTester::mE_3000
MonitorElement * mE_3000
Definition: CaloJetTester.h:59
CaloJetTester::mConstituents_80
MonitorElement * mConstituents_80
Definition: CaloJetTester.h:70
CaloJetTester::mCaloMEx_3000
MonitorElement * mCaloMEx_3000
Definition: CaloJetTester.h:216
CaloJetTester::mEmEnergyInEB_80
MonitorElement * mEmEnergyInEB_80
Definition: CaloJetTester.h:111
CaloJetTester::mPhiFineBin
MonitorElement * mPhiFineBin
Definition: CaloJetTester.h:56
CaloJetTester::mpTScale1DB_600_900
MonitorElement * mpTScale1DB_600_900
Definition: CaloJetTester.h:183
CaloJetTester::mEmEnergyInEE_80
MonitorElement * mEmEnergyInEE_80
Definition: CaloJetTester.h:112
j
int j
Definition: DBlmapReader.cc:9
CaloJetTester::mHadEnergyInHO
MonitorElement * mHadEnergyInHO
Definition: CaloJetTester.h:98
reco::GenJet
Jets made from MC generator particles.
Definition: GenJet.h:25
CaloJetTester::mHBTime
MonitorElement * mHBTime
Definition: CaloJetTester.h:229
CaloJetTester::mpTScale1DF_600_900
MonitorElement * mpTScale1DF_600_900
Definition: CaloJetTester.h:185
CaloJetTester::mEMJetEnergyProfile
MonitorElement * mEMJetEnergyProfile
Definition: CaloJetTester.h:212
CaloJetTester::mGenEta
MonitorElement * mGenEta
Definition: CaloJetTester.h:133
CaloJetTester::mHadEnergyInHO_80
MonitorElement * mHadEnergyInHO_80
Definition: CaloJetTester.h:102
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:828
CaloJetTester::mReverseMatchEnergyFraction
MonitorElement * mReverseMatchEnergyFraction
Definition: CaloJetTester.h:149
CaloJetTester::mCaloMETPhi
MonitorElement * mCaloMETPhi
Definition: CaloJetTester.h:221
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
CaloJetTester::mHETime
MonitorElement * mHETime
Definition: CaloJetTester.h:231
CaloJetTester::mHadEnergyInHF
MonitorElement * mHadEnergyInHF
Definition: CaloJetTester.h:100
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:355
CaloJetTester::mMjj_3000
MonitorElement * mMjj_3000
Definition: CaloJetTester.h:93
CaloJetTester::mEETime
MonitorElement * mEETime
Definition: CaloJetTester.h:239
CaloJetTester::mGenPt
MonitorElement * mGenPt
Definition: CaloJetTester.h:135
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
CaloJetTester::mPthat_3000
MonitorElement * mPthat_3000
Definition: CaloJetTester.h:130
CaloJetTester::mPhi
MonitorElement * mPhi
Definition: CaloJetTester.h:55
CaloJetTester::mHFShort
MonitorElement * mHFShort
Definition: CaloJetTester.h:123
CaloJetTester::mEScaleFineBin
MonitorElement * mEScaleFineBin
Definition: CaloJetTester.h:158
CaloJetTester::mEtaFineBin1p
MonitorElement * mEtaFineBin1p
Definition: CaloJetTester.h:49
CaloJetTester::numberofevents
MonitorElement * numberofevents
Definition: CaloJetTester.h:44
jptDQMConfig_cff.etaMin
int etaMin
Definition: jptDQMConfig_cff.py:52
CaloJetTester::mEBTime
MonitorElement * mEBTime
Definition: CaloJetTester.h:237
CaloJetTester::mHadEnergyProfile
MonitorElement * mHadEnergyProfile
Definition: CaloJetTester.h:208
CaloJetTester::mPt_3000
MonitorElement * mPt_3000
Definition: CaloJetTester.h:65
CaloJetTester::mpTScale_200_300_d
MonitorElement * mpTScale_200_300_d
Definition: CaloJetTester.h:173
CaloJetTester::mMass_3000
MonitorElement * mMass_3000
Definition: CaloJetTester.h:68
reco::LeafCandidate::px
virtual double px() const
x coordinate of momentum vector
Definition: LeafCandidate.h:131
CaloJetTester::mpTScaleE_d
MonitorElement * mpTScaleE_d
Definition: CaloJetTester.h:164
CaloJetTester::mPthat_80
MonitorElement * mPthat_80
Definition: CaloJetTester.h:129
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.cc:127
CaloJetTester::mpTScale_60_120_d
MonitorElement * mpTScale_60_120_d
Definition: CaloJetTester.h:172
CaloJetTester::mGenPt_80
MonitorElement * mGenPt_80
Definition: CaloJetTester.h:136
CaloJetTester::mHOTime
MonitorElement * mHOTime
Definition: CaloJetTester.h:233
CaloJetTester::mlinEScale
MonitorElement * mlinEScale
Definition: CaloJetTester.h:154
CaloJetTester::mDeltaEta
MonitorElement * mDeltaEta
Definition: CaloJetTester.h:151
CaloJetTester::mpTScale1DB_2700_3500
MonitorElement * mpTScale1DB_2700_3500
Definition: CaloJetTester.h:186
edm::Handle::product
T const * product() const
Definition: Handle.h:74
CaloJetTester::mHFEne
MonitorElement * mHFEne
Definition: CaloJetTester.h:234
CaloJetTester::beginJob
virtual void beginJob()
Definition: CaloJetTester.cc:366
CaloJetTester::mpTScale1DE_200_300
MonitorElement * mpTScale1DE_200_300
Definition: CaloJetTester.h:181
JetMatchingTools::getGenParticles
std::vector< const reco::GenParticle * > getGenParticles(const reco::CaloJet &fJet, bool fVerbose=true)
GenParticles for CaloJet.
Definition: JetMatchingTools.cc:249
CaloJetTester::mHFLong_80
MonitorElement * mHFLong_80
Definition: CaloJetTester.h:121
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:25
CaloJetTester::mCaloSumET_3000
MonitorElement * mCaloSumET_3000
Definition: CaloJetTester.h:223
CaloJetTester::mpTScale1DB_200_300
MonitorElement * mpTScale1DB_200_300
Definition: CaloJetTester.h:180
CaloJetTester::mpTScaleE_s
MonitorElement * mpTScaleE_s
Definition: CaloJetTester.h:161
CaloJetTester::mHFTotal
MonitorElement * mHFTotal
Definition: CaloJetTester.h:117
CaloJetTester::mEmEnergyInEE
MonitorElement * mEmEnergyInEE
Definition: CaloJetTester.h:109
CaloJetTester::mMass_80
MonitorElement * mMass_80
Definition: CaloJetTester.h:67
CaloJetTester::mNJetsEtaC
MonitorElement * mNJetsEtaC
Definition: CaloJetTester.h:85
jptDQMConfig_cff.etaMax
float etaMax
Definition: jptDQMConfig_cff.py:9
CaloJetTester::mHEEne
MonitorElement * mHEEne
Definition: CaloJetTester.h:230
CaloJetTester::mEFirst_80
MonitorElement * mEFirst_80
Definition: CaloJetTester.h:79
CaloJetTester::mE
MonitorElement * mE
Definition: CaloJetTester.h:57
CaloJetTester::mEta
MonitorElement * mEta
Definition: CaloJetTester.h:47
CaloJetTester::mN90
MonitorElement * mN90
Definition: CaloJetTester.h:126
CaloJetTester::mEtaFineBin1m
MonitorElement * mEtaFineBin1m
Definition: CaloJetTester.h:52
CaloJetTester::mJetEnergyProfile
MonitorElement * mJetEnergyProfile
Definition: CaloJetTester.h:210
CaloJetTester::mCaloMET_3000
MonitorElement * mCaloMET_3000
Definition: CaloJetTester.h:220
CaloJetTester::mP_3000
MonitorElement * mP_3000
Definition: CaloJetTester.h:62
CaloJetTester::mPtFirst_80
MonitorElement * mPtFirst_80
Definition: CaloJetTester.h:82
CaloJetTester::mP
MonitorElement * mP
Definition: CaloJetTester.h:60
CaloJetTester::mGenPhiFirst
MonitorElement * mGenPhiFirst
Definition: CaloJetTester.h:141
CaloJetTester::mMaxEInEmTowers
MonitorElement * mMaxEInEmTowers
Definition: CaloJetTester.h:96
jptDQMConfig_cff.etaBins
tuple etaBins
Definition: jptDQMConfig_cff.py:53
HiRecoJets_cff.caloTowers
tuple caloTowers
Definition: HiRecoJets_cff.py:14
CaloJetTester::mEtaFirst
MonitorElement * mEtaFirst
Definition: CaloJetTester.h:76
CaloJetTester::mpTScaleF_d
MonitorElement * mpTScaleF_d
Definition: CaloJetTester.h:165
CaloJetTester::mGenEnergyFractionThreshold
double mGenEnergyFractionThreshold
Definition: CaloJetTester.h:200
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:642
CaloJetTester::mEtaFineBin3p
MonitorElement * mEtaFineBin3p
Definition: CaloJetTester.h:51
CaloJetTester::mpTScale1D_600_900
MonitorElement * mpTScale1D_600_900
Definition: CaloJetTester.h:191
CaloJetTester::mEFirst_3000
MonitorElement * mEFirst_3000
Definition: CaloJetTester.h:80
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:139
CaloJetTester::mpTScale_60_120_s
MonitorElement * mpTScale_60_120_s
Definition: CaloJetTester.h:167
JetMatchingTools::overlapEnergyFraction
double overlapEnergyFraction(const std::vector< const reco::GenParticle * > &fObject, const std::vector< const reco::GenParticle * > &fReference) const
energy overlap
Definition: JetMatchingTools.cc:316
CaloJetTester::mHFLong
MonitorElement * mHFLong
Definition: CaloJetTester.h:120
CaloJetTester::mHBEne
MonitorElement * mHBEne
Definition: CaloJetTester.h:228
CaloJetTester::mMatchedGenJetsEta
MonitorElement * mMatchedGenJetsEta
Definition: CaloJetTester.h:147
reco::LeafCandidate::py
virtual double py() const
y coordinate of momentum vector
Definition: LeafCandidate.h:133
CaloJetTester::mpTScale1D_2700_3500
MonitorElement * mpTScale1D_2700_3500
Definition: CaloJetTester.h:192
CaloJetTester::mCaloSumET
MonitorElement * mCaloSumET
Definition: CaloJetTester.h:222
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:232
CaloJetTester::mDeltaPhi
MonitorElement * mDeltaPhi
Definition: CaloJetTester.h:152
CaloJetTester::mEmEnergyInEB
MonitorElement * mEmEnergyInEB
Definition: CaloJetTester.h:108
CaloJetTester::mAllGenJetsEta
MonitorElement * mAllGenJetsEta
Definition: CaloJetTester.h:146
CaloJetTester::mGenPhi
MonitorElement * mGenPhi
Definition: CaloJetTester.h:134
CaloJetTester::CaloJetTester
CaloJetTester(const edm::ParameterSet &)
Definition: CaloJetTester.cc:48
CaloJetTester::mEScale
MonitorElement * mEScale
Definition: CaloJetTester.h:153
DQMStore::bookProfile2D
MonitorElement * bookProfile2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ, const char *option="s")
Definition: DQMStore.cc:972
CaloJetTester::mInputGenCollection
edm::InputTag mInputGenCollection
Definition: CaloJetTester.h:36
Generated for CMSSW Reference Manual by doxygen 1.8.5