ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitPhysics/Validation/src/JetValidationMod.cc

Comparing UserCode/MitPhysics/Validation/src/JetValidationMod.cc (file contents):
Revision 1.1 by loizides, Wed Oct 15 06:05:02 2008 UTC vs.
Revision 1.5 by loizides, Fri Jul 10 14:04:44 2009 UTC

#	Line 1 | Line 1
1		// $Id$
2
3		#include "MitPhysics/Validation/interface/JetValidationMod.h"
4	<	#include <TH1D.h>
5	<	#include <TH2D.h>
6	<	#include <TH3D.h>
4	>	#include "MitCommon/MathTools/interface/MathUtils.h"
5	>	#include "MitAna/DataTree/interface/GenJetCol.h"
6	>	#include "MitAna/DataTree/interface/JetCol.h"
7		#include "MitAna/DataTree/interface/Names.h"
8		#include "MitAna/DataCont/interface/ObjArray.h"
9	–	#include "MitCommon/MathTools/interface/MathUtils.h"
9		#include <map>
10	+	#include <TH1D.h>
11	+	#include <TH2D.h>
12	+	#include <TH3D.h>
13
14		using namespace mithep;
15
#	Line 24 | Line 26 | ClassImp(mithep::JetValidationMod)
26		fIC5Jets(0),
27		fSC5Jets(0),
28		fIC5GenJets(0),
29	<	fSC5GenJets(0),
28	<	fNEventsProcessed(0)
29	>	fSC5GenJets(0)
30		{
31		// Constructor.
32		}
33
34		//--------------------------------------------------------------------------------------------------
34	–	void JetValidationMod::Begin()
35	–	{
36	–	// Run startup code on the client machine. For this module, we dont do
37	–	// anything here.
38	–	}
39	–
40	–	//--------------------------------------------------------------------------------------------------
35		void JetValidationMod::Process()
36		{
37		// Process entries of the tree.
38
45	–	fNEventsProcessed++;
46	–	if (fNEventsProcessed % 1000 == 0 || fPrintDebug)
47	–	cerr << endl << "Process Event " << fNEventsProcessed << endl;
48	–
39		//************************************************************************************************
40		//
41		// For IterativeCone 0.5 Jets
#	Line 73 | Line 63 | void JetValidationMod::Process()
63
64		std::map <GenJet*, Jet*> IC5GenJetToIC5JetMap;
65		for (UInt_t i=0; i<GoodIC5GenJets->GetEntries(); ++i) {
66	<	double minDR = 5000;
66	>	Double_t  minDR = 5000;
67		int closestMatchIndex = -1;
68		for (UInt_t j=0; j<GoodIC5Jets->GetEntries(); j++) {
69	<	double tempDR = MathUtils::DeltaR(GoodIC5GenJets->At(i)->Mom(), GoodIC5Jets->At(j)->Mom());
69	>	Double_t  tempDR = MathUtils::DeltaR(GoodIC5GenJets->At(i)->Mom(), GoodIC5Jets->At(j)->Mom());
70		if (tempDR < minDR) {
71		closestMatchIndex = j;
72		minDR = tempDR;
73		}
74		}
75		if (closestMatchIndex > -1)
76	<	if (fPrintDebug) cerr << "GoodIC5GenJets->At(i) " << GoodIC5GenJets->At(i)->Pt() << " " << GoodIC5GenJets->At(i)->Eta() << " " << GoodIC5GenJets->At(i)->Phi() << "   GoodIC5Jets->At(" << closestMatchIndex << ") " << GoodIC5Jets->At(closestMatchIndex)->Pt() << " " << GoodIC5Jets->At(closestMatchIndex)->Eta() << " " << GoodIC5Jets->At(closestMatchIndex)->Phi() << "    DR = " << minDR << endl;
76	>	if (fPrintDebug) cerr << "GoodIC5GenJets->At(i) " << GoodIC5GenJets->At(i)->Pt() << " "
77	>	<< GoodIC5GenJets->At(i)->Eta() << " " << GoodIC5GenJets->At(i)->Phi()
78	>	<< "   GoodIC5Jets->At(" << closestMatchIndex << ") "
79	>	<< GoodIC5Jets->At(closestMatchIndex)->Pt() << " "
80	>	<< GoodIC5Jets->At(closestMatchIndex)->Eta() << " "
81	>	<< GoodIC5Jets->At(closestMatchIndex)->Phi() << "    DR = "
82	>	<< minDR << endl;
83		if (minDR < 0.5)
84		{
85		IC5GenJetToIC5JetMap[GoodIC5GenJets->At(i)] = GoodIC5Jets->At(closestMatchIndex);
#	Line 92 | Line 88 | void JetValidationMod::Process()
88
89		std::map <Jet*, GenJet*> IC5JetToIC5GenJetMap;
90		for (UInt_t i=0; i<GoodIC5Jets->GetEntries(); ++i) {
91	<	double minDR = 5000;
91	>	Double_t  minDR = 5000;
92		int closestMatchIndex = -1;
93		for (UInt_t j=0; j<GoodIC5GenJets->GetEntries(); j++) {
94	<	double tempDR = MathUtils::DeltaR(GoodIC5Jets->At(i)->Mom(), GoodIC5GenJets->At(j)->Mom());
94	>	Double_t  tempDR = MathUtils::DeltaR(GoodIC5Jets->At(i)->Mom(), GoodIC5GenJets->At(j)->Mom());
95		if (tempDR < minDR) {
96		closestMatchIndex = j;
97		minDR = tempDR;
#	Line 119 | Line 115 | void JetValidationMod::Process()
115		map<GenJet*, Jet*>::iterator iter =
116		IC5GenJetToIC5JetMap.find(GoodIC5GenJets->At(i));
117		if (iter != IC5GenJetToIC5JetMap.end()) {
118	<	double dR = MathUtils::DeltaR(GoodIC5GenJets->At(i)->Mom(), iter->second->Mom());
118	>	Double_t  dR = MathUtils::DeltaR(GoodIC5GenJets->At(i)->Mom(), iter->second->Mom());
119
120		fIC5GenJetRecoJetDeltaR->Fill(dR);
121		fIC5GenJetRecoJetDeltaEta->Fill(fabs(GoodIC5GenJets->At(i)->Eta() - iter->second->Eta()));
122	<	fIC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodIC5GenJets->At(i)->Phi(),iter->second->Phi()));
122	>	fIC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodIC5GenJets->At(i)->Phi(),
123	>	iter->second->Phi()));
124
125		if (fabs(iter->second->Eta()) <= 1.4)
126	<	fIC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),GoodIC5GenJets->At(i)->Pt());
126	>	fIC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
127	>	GoodIC5GenJets->At(i)->Pt());
128		else if (fabs(iter->second->Eta()) <= 3.0)
129	<	fIC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),GoodIC5GenJets->At(i)->Pt());
129	>	fIC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
130	>	GoodIC5GenJets->At(i)->Pt());
131		else if (fabs(iter->second->Eta()) <= 5.0)
132	<	fIC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),GoodIC5GenJets->At(i)->Pt());
132	>	fIC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
133	>	GoodIC5GenJets->At(i)->Pt());
134
135		//avoid the low pt region
136		if (iter->second->Pt() > 30)
137	<	fIC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),iter->second->Eta());
137	>	fIC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
138	>	iter->second->Eta());
139
140		if (fabs(iter->second->Eta()) <= 3.0) {
141		fIC5CentralGenJetRecoJetDeltaR->Fill(dR);
#	Line 152 | Line 153 | void JetValidationMod::Process()
153		if (GoodIC5GenJets->At(i)->Pt() > 60.0 && GoodIC5GenJets->At(i)->Pt() < 80.0)
154		fIC5NMatchedCaloJetsVsGenJetEta_Pt60To80->Fill(GoodIC5GenJets->At(i)->Eta());
155
156	<	fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(iter->second->Pt() * iter->second->L2RelativeCorrectionScale() * iter->second->L3AbsoluteCorrectionScale() / GoodIC5GenJets->At(i)->Pt(), GoodIC5GenJets->At(i)->Pt());
156	>	fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(iter->second->Pt() *
157	>	iter->second->L2RelativeCorrectionScale() *
158	>	iter->second->L3AbsoluteCorrectionScale() /
159	>	GoodIC5GenJets->At(i)->Pt(),
160	>	GoodIC5GenJets->At(i)->Pt());
161		}
162		}
163
#	Line 161 | Line 166 | void JetValidationMod::Process()
166		IC5JetToIC5GenJetMap.find(GoodIC5Jets->At(i));
167		//unmatched ones
168		if (iter == IC5JetToIC5GenJetMap.end()) {
169	<	fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(GoodIC5Jets->At(i)->Pt() * GoodIC5Jets->At(i)->L2RelativeCorrectionScale() * GoodIC5Jets->At(i)->L3AbsoluteCorrectionScale());
169	>	fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(
170	>	GoodIC5Jets->At(i)->Pt() * GoodIC5Jets->At(i)->L2RelativeCorrectionScale() *
171	>	GoodIC5Jets->At(i)->L3AbsoluteCorrectionScale());
172		fIC5NUnmatchedCalojetsVsCorrectedCaloJetEta->Fill(GoodIC5Jets->At(i)->Eta());
173		}
174		}
#	Line 194 | Line 201 | void JetValidationMod::Process()
201
202		std::map <GenJet*, Jet*> SC5GenJetToSC5JetMap;
203		for (UInt_t i=0; i<GoodSC5GenJets->GetEntries(); ++i) {
204	<	double minDR = 5000;
204	>	Double_t  minDR = 5000;
205		int closestMatchIndex = -1;
206		for (UInt_t j=0; j<GoodSC5Jets->GetEntries(); j++) {
207	<	double tempDR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), GoodSC5Jets->At(j)->Mom());
207	>	Double_t  tempDR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), GoodSC5Jets->At(j)->Mom());
208		if (tempDR < minDR) {
209		closestMatchIndex = j;
210		minDR = tempDR;
211		}
212		}
213		if (closestMatchIndex > -1)
214	<	if (fPrintDebug) cerr << "GoodSC5GenJets->At(i) " << GoodSC5GenJets->At(i)->Pt() << " " << GoodSC5GenJets->At(i)->Eta() << " " << GoodSC5GenJets->At(i)->Phi() << "   GoodSC5Jets->At(" << closestMatchIndex << ") " << GoodSC5Jets->At(closestMatchIndex)->Pt() << " " << GoodSC5Jets->At(closestMatchIndex)->Eta() << " " << GoodSC5Jets->At(closestMatchIndex)->Phi() << "    DR = " << minDR << endl;
214	>	if (fPrintDebug) cerr << "GoodSC5GenJets->At(i) " << GoodSC5GenJets->At(i)->Pt() << " "
215	>	<< GoodSC5GenJets->At(i)->Eta() << " " << GoodSC5GenJets->At(i)->Phi()
216	>	<< "   GoodSC5Jets->At(" << closestMatchIndex << ") "
217	>	<< GoodSC5Jets->At(closestMatchIndex)->Pt() << " "
218	>	<< GoodSC5Jets->At(closestMatchIndex)->Eta() << " "
219	>	<< GoodSC5Jets->At(closestMatchIndex)->Phi() << "    DR = "
220	>	<< minDR << endl;
221		if (minDR < 0.5)
222		{
223		SC5GenJetToSC5JetMap[GoodSC5GenJets->At(i)] = GoodSC5Jets->At(closestMatchIndex);
#	Line 212 | Line 225 | void JetValidationMod::Process()
225		}
226		std::map <Jet*, GenJet*> SC5JetToSC5GenJetMap;
227		for (UInt_t i=0; i<GoodSC5Jets->GetEntries(); ++i) {
228	<	double minDR = 5000;
228	>	Double_t  minDR = 5000;
229		int closestMatchIndex = -1;
230		for (UInt_t j=0; j<GoodSC5GenJets->GetEntries(); j++) {
231	<	double tempDR = MathUtils::DeltaR(GoodSC5Jets->At(i)->Mom(), GoodSC5GenJets->At(j)->Mom());
231	>	Double_t  tempDR = MathUtils::DeltaR(GoodSC5Jets->At(i)->Mom(), GoodSC5GenJets->At(j)->Mom());
232		if (tempDR < minDR) {
233		closestMatchIndex = j;
234		minDR = tempDR;
#	Line 238 | Line 251 | void JetValidationMod::Process()
251		map<GenJet*, Jet*>::iterator iter =
252		SC5GenJetToSC5JetMap.find(GoodSC5GenJets->At(i));
253		if (iter != SC5GenJetToSC5JetMap.end()) {
254	<	double dR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), iter->second->Mom());
254	>	Double_t  dR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), iter->second->Mom());
255
256		fSC5GenJetRecoJetDeltaR->Fill(dR);
257		fSC5GenJetRecoJetDeltaEta->Fill(fabs(GoodSC5GenJets->At(i)->Eta() - iter->second->Eta()));
258	<	fSC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodSC5GenJets->At(i)->Phi(),iter->second->Phi()));
258	>	fSC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodSC5GenJets->At(i)->Phi(),
259	>	iter->second->Phi()));
260
261		if (fabs(iter->second->Eta()) <= 1.4)
262	<	fSC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),GoodSC5GenJets->At(i)->Pt());
262	>	fSC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
263	>	GoodSC5GenJets->At(i)->Pt());
264		else if (fabs(iter->second->Eta()) <= 3.0)
265	<	fSC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),GoodSC5GenJets->At(i)->Pt());
265	>	fSC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
266	>	GoodSC5GenJets->At(i)->Pt());
267		else if (fabs(iter->second->Eta()) <= 5.0)
268	<	fSC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),GoodSC5GenJets->At(i)->Pt());
268	>	fSC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
269	>	GoodSC5GenJets->At(i)->Pt());
270		//avoid the low pt region
271		if (iter->second->Pt() > 30)
272	<	fSC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),iter->second->Eta());
272	>	fSC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
273	>	iter->second->Eta());
274
275		if (fabs(iter->second->Eta()) <= 3.0) {
276		fSC5CentralGenJetRecoJetDeltaR->Fill(dR);
#	Line 269 | Line 287 | void JetValidationMod::Process()
287		fSC5NMatchedCaloJetsVsGenJetEta_Pt30To40->Fill(GoodSC5GenJets->At(i)->Eta());
288		if (GoodSC5GenJets->At(i)->Pt() > 60.0 && GoodSC5GenJets->At(i)->Pt() < 80.0)
289		fSC5NMatchedCaloJetsVsGenJetEta_Pt60To80->Fill(GoodSC5GenJets->At(i)->Eta());
290	<	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(iter->second->Pt() * iter->second->L2RelativeCorrectionScale() * iter->second->L3AbsoluteCorrectionScale() / GoodSC5GenJets->At(i)->Pt(), GoodSC5GenJets->At(i)->Pt());
290	>	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(
291	>	iter->second->Pt() * iter->second->L2RelativeCorrectionScale() *
292	>	iter->second->L3AbsoluteCorrectionScale() / GoodSC5GenJets->At(i)->Pt(),
293	>	GoodSC5GenJets->At(i)->Pt());
294		}
295		}
296
#	Line 278 | Line 299 | void JetValidationMod::Process()
299		SC5JetToSC5GenJetMap.find(GoodSC5Jets->At(i));
300		//unmatched ones
301		if (iter == SC5JetToSC5GenJetMap.end()) {
302	<	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(GoodSC5Jets->At(i)->Pt() * GoodSC5Jets->At(i)->L2RelativeCorrectionScale() * GoodSC5Jets->At(i)->L3AbsoluteCorrectionScale());
302	>	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(
303	>	GoodSC5Jets->At(i)->Pt() * GoodSC5Jets->At(i)->L2RelativeCorrectionScale() *
304	>	GoodSC5Jets->At(i)->L3AbsoluteCorrectionScale());
305		fSC5NUnmatchedCalojetsVsCorrectedCaloJetEta->Fill(GoodSC5Jets->At(i)->Eta());
306		}
307		}
308
309		}
310
288	–
311		//--------------------------------------------------------------------------------------------------
312		void JetValidationMod::SlaveBegin()
313		{
#	Line 298 | Line 320 | void JetValidationMod::SlaveBegin()
320		ReqBranch(fIC5JetName,                fIC5Jets);
321		ReqBranch(fSC5JetName,                fSC5Jets);
322
323	<	fIC5GenJetRecoJetDeltaR = new TH1D("hIC5GenJetRecoJetDeltaR",";IC5GenJetRecoJetDeltaR;#",100,0.,1.0);
324	<	fIC5GenJetRecoJetDeltaEta = new TH1D("hIC5GenJetRecoJetDeltaEta",";IC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
325	<	fIC5GenJetRecoJetDeltaPhi = new TH1D("hIC5GenJetRecoJetDeltaPhi",";IC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
326	<
327	<	fIC5JetResponseVsGenJetPtInBarrel = new TH2D("hIC5JetResponseVsGenJetPtInBarrel",";IC5JetResponseVsGenJetPtInBarrel;#",500,0.,5.0,5000,0,5000);
328	<	fIC5JetResponseVsGenJetPtInEndcap = new TH2D("hIC5JetResponseVsGenJetPtInEndcap",";IC5JetResponseVsGenJetPtInEndcap;#",500,0.,5.0,5000,0,5000);
329	<	fIC5JetResponseVsGenJetPtForward = new TH2D("hIC5JetResponseVsGenJetPtForward",";IC5JetResponseVsGenJetPtForward;#",500,0.,5.0,5000,0,5000);
330	<	fIC5JetResponseVsCaloJetEta = new TH2D("hIC5JetResponseVsCaloJetEta",";IC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
331	<
332	<
333	<
334	<	fIC5CentralGenJetRecoJetDeltaR = new TH1D("hIC5CentralGenJetRecoJetDeltaR",";IC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
335	<	fIC5ForwardGenJetRecoJetDeltaR = new TH1D("hIC5ForwardGenJetRecoJetDeltaR",";IC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
336	<
337	<	fIC5NMatchedCaloJetsVsGenJetPt = new TH1D("hIC5NMatchedCaloJetsVsGenJetPt",";IC5GenJetPt;#",5000,0.,5000);
338	<	fIC5NMatchedCaloJetsVsGenJetEta = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta",";IC5GenJetEta;#",100,-5.0,5.0);
339	<	fIC5NMatchedCaloJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt20To30",";IC5GenJetEta;#",100,-5.0,5.0);
340	<	fIC5NMatchedCaloJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt30To40",";IC5GenJetEta;#",100,-5.0,5.0);
341	<	fIC5NMatchedCaloJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt60To80",";IC5GenJetEta;#",100,-5.0,5.0);
323	>	fIC5GenJetRecoJetDeltaR =
324	>	new TH1D("hIC5GenJetRecoJetDeltaR",";IC5GenJetRecoJetDeltaR;#",100,0.,1.0);
325	>	fIC5GenJetRecoJetDeltaEta =
326	>	new TH1D("hIC5GenJetRecoJetDeltaEta",";IC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
327	>	fIC5GenJetRecoJetDeltaPhi =
328	>	new TH1D("hIC5GenJetRecoJetDeltaPhi",";IC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
329	>
330	>	fIC5JetResponseVsGenJetPtInBarrel = new TH2D("hIC5JetResponseVsGenJetPtInBarrel",
331	>	";IC5JetResponseVsGenJetPtInBarrel;#",
332	>	500,0.,5.0,5000,0,5000);
333	>	fIC5JetResponseVsGenJetPtInEndcap = new TH2D("hIC5JetResponseVsGenJetPtInEndcap",
334	>	";IC5JetResponseVsGenJetPtInEndcap;#",
335	>	500,0.,5.0,5000,0,5000);
336	>	fIC5JetResponseVsGenJetPtForward = new TH2D("hIC5JetResponseVsGenJetPtForward",
337	>	";IC5JetResponseVsGenJetPtForward;#",
338	>	500,0.,5.0,5000,0,5000);
339	>	fIC5JetResponseVsCaloJetEta = new TH2D("hIC5JetResponseVsCaloJetEta",
340	>	";IC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
341	>	fIC5CentralGenJetRecoJetDeltaR =
342	>	new TH1D("hIC5CentralGenJetRecoJetDeltaR",";IC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
343	>	fIC5ForwardGenJetRecoJetDeltaR = new TH1D("hIC5ForwardGenJetRecoJetDeltaR",
344	>	";IC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
345	>	fIC5NMatchedCaloJetsVsGenJetPt = new TH1D("hIC5NMatchedCaloJetsVsGenJetPt",
346	>	";IC5GenJetPt;#",5000,0.,5000);
347	>	fIC5NMatchedCaloJetsVsGenJetEta = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta",
348	>	";IC5GenJetEta;#",100,-5.0,5.0);
349	>	fIC5NMatchedCaloJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt20To30",
350	>	";IC5GenJetEta;#",100,-5.0,5.0);
351	>	fIC5NMatchedCaloJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt30To40",
352	>	";IC5GenJetEta;#",100,-5.0,5.0);
353	>	fIC5NMatchedCaloJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt60To80",
354	>	";IC5GenJetEta;#",100,-5.0,5.0);
355		fIC5NGenJetsVsGenJetPt = new TH1D("hIC5NGenJetsVsGenJetPt",";IC5GenJetPt;#",5000,0.,5000);
356		fIC5NGenJetsVsGenJetEta = new TH1D("hIC5NGenJetsVsGenJetEta",";IC5GenJetEta;#",100,-5.,5.0);
357	<	fIC5NGenJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt20To30",";IC5GenJetEta;#",100,-5.,5.0);
358	<	fIC5NGenJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt30To40",";IC5GenJetEta;#",100,-5.,5.0);
359	<	fIC5NGenJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt60To80",";IC5GenJetEta;#",100,-5.,5.0);
360	<
357	>	fIC5NGenJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt20To30",
358	>	";IC5GenJetEta;#",100,-5.,5.0);
359	>	fIC5NGenJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt30To40",
360	>	";IC5GenJetEta;#",100,-5.,5.0);
361	>	fIC5NGenJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt60To80",
362	>	";IC5GenJetEta;#",100,-5.,5.0);
363		fIC5CaloJetsPt = new TH1D("hIC5CaloJetsPt",";IC5CaloJetPt;#",5000,0.,5000);
364		fIC5CaloJetsEta = new TH1D("hIC5CaloJetsEta",";IC5CaloJetEta;#",100,-5.0,5.0);
365	<	fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt = new TH1D("hIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",";IC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
366	<	fIC5NUnmatchedCalojetsVsCorrectedCaloJetEta = new TH1D("hIC5NUnmatchedCalojetsVsCorrectedCaloJetEta",";IC5CorrectedCaloJetEta;#",100,-5.0,5.0);
367	<
368	<	fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt = new TH2D("hIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",";IC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
369	<
365	>	fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt =
366	>	new TH1D("hIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",
367	>	";IC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
368	>	fIC5NUnmatchedCalojetsVsCorrectedCaloJetEta =
369	>	new TH1D("hIC5NUnmatchedCalojetsVsCorrectedCaloJetEta",
370	>	";IC5CorrectedCaloJetEta;#",100,-5.0,5.0);
371	>	fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt =
372	>	new TH2D("hIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",
373	>	";IC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
374		fIC5NCaloJets = new TH1D("hIC5NCaloJets",";IC5NCaloJets;#",100,0.,15);
375	<
376	<
377	<	fSC5GenJetRecoJetDeltaR = new TH1D("hSC5GenJetRecoJetDeltaR",";SC5GenJetRecoJetDeltaR;#",100,0.,1.0);
378	<	fSC5GenJetRecoJetDeltaEta = new TH1D("hSC5GenJetRecoJetDeltaEta",";SC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
379	<	fSC5GenJetRecoJetDeltaPhi = new TH1D("hSC5GenJetRecoJetDeltaPhi",";SC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
380	<
381	<	fSC5JetResponseVsGenJetPtInBarrel = new TH2D("hSC5JetResponseVsGenJetPtInBarrel",";SC5JetResponseVsGenJetPtInBarrel;#",500,0.,5.0,5000,0,5000);
382	<	fSC5JetResponseVsGenJetPtInEndcap = new TH2D("hSC5JetResponseVsGenJetPtInEndcap",";SC5JetResponseVsGenJetPtInEndcap;#",500,0.,5.0,5000,0,5000);
383	<	fSC5JetResponseVsGenJetPtForward = new TH2D("hSC5JetResponseVsGenJetPtForward",";SC5JetResponseVsGenJetPtForward;#",500,0.,5.0,5000,0,5000);
384	<
385	<	fSC5JetResponseVsCaloJetEta = new TH2D("hSC5JetResponseVsCaloJetEta",";SC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
386	<
387	<
388	<
389	<	fSC5CentralGenJetRecoJetDeltaR = new TH1D("hSC5CentralGenJetRecoJetDeltaR",";SC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
390	<	fSC5ForwardGenJetRecoJetDeltaR = new TH1D("hSC5ForwardGenJetRecoJetDeltaR",";SC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
391	<
392	<	fSC5NMatchedCaloJetsVsGenJetPt = new TH1D("hSC5NMatchedCaloJetsVsGenJetPt",";SC5GenJetPt;#",5000,0.,5000);
393	<	fSC5NMatchedCaloJetsVsGenJetEta = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta",";SC5GenJetEta;#",100,-5.0,5.0);
394	<	fSC5NMatchedCaloJetsVsGenJetEta_Pt20To30 = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.0,5.0);
395	<	fSC5NMatchedCaloJetsVsGenJetEta_Pt30To40 = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.0,5.0);
396	<	fSC5NMatchedCaloJetsVsGenJetEta_Pt60To80 = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.0,5.0);
375	>	fSC5GenJetRecoJetDeltaR = new TH1D("hSC5GenJetRecoJetDeltaR",
376	>	";SC5GenJetRecoJetDeltaR;#",100,0.,1.0);
377	>	fSC5GenJetRecoJetDeltaEta = new TH1D("hSC5GenJetRecoJetDeltaEta",
378	>	";SC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
379	>	fSC5GenJetRecoJetDeltaPhi = new TH1D("hSC5GenJetRecoJetDeltaPhi",
380	>	";SC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
381	>	fSC5JetResponseVsGenJetPtInBarrel =
382	>	new TH2D("hSC5JetResponseVsGenJetPtInBarrel",
383	>	";SC5JetResponseVsGenJetPtInBarrel;#",500,0.,5.0,5000,0,5000);
384	>	fSC5JetResponseVsGenJetPtInEndcap =
385	>	new TH2D("hSC5JetResponseVsGenJetPtInEndcap",
386	>	";SC5JetResponseVsGenJetPtInEndcap;#",500,0.,5.0,5000,0,5000);
387	>	fSC5JetResponseVsGenJetPtForward =
388	>	new TH2D("hSC5JetResponseVsGenJetPtForward",
389	>	";SC5JetResponseVsGenJetPtForward;#",500,0.,5.0,5000,0,5000);
390	>
391	>	fSC5JetResponseVsCaloJetEta = new TH2D("hSC5JetResponseVsCaloJetEta",
392	>	";SC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
393	>	fSC5CentralGenJetRecoJetDeltaR =
394	>	new TH1D("hSC5CentralGenJetRecoJetDeltaR",
395	>	";SC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
396	>	fSC5ForwardGenJetRecoJetDeltaR =
397	>	new TH1D("hSC5ForwardGenJetRecoJetDeltaR",
398	>	";SC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
399	>	fSC5NMatchedCaloJetsVsGenJetPt =
400	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetPt",";SC5GenJetPt;#",5000,0.,5000);
401	>	fSC5NMatchedCaloJetsVsGenJetEta =
402	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta",";SC5GenJetEta;#",100,-5.0,5.0);
403	>	fSC5NMatchedCaloJetsVsGenJetEta_Pt20To30 =
404	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.0,5.0);
405	>	fSC5NMatchedCaloJetsVsGenJetEta_Pt30To40 =
406	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.0,5.0);
407	>	fSC5NMatchedCaloJetsVsGenJetEta_Pt60To80 =
408	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.0,5.0);
409		fSC5NGenJetsVsGenJetPt = new TH1D("hSC5NGenJetsVsGenJetPt",";SC5GenJetPt;#",5000,0.,5000);
410		fSC5NGenJetsVsGenJetEta = new TH1D("hSC5NGenJetsVsGenJetEta",";SC5GenJetEta;#",100,-5.,5.0);
411	<	fSC5NGenJetsVsGenJetEta_Pt20To30 = new TH1D("hSC5NGenJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.,5.0);
412	<	fSC5NGenJetsVsGenJetEta_Pt30To40 = new TH1D("hSC5NGenJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.,5.0);
413	<	fSC5NGenJetsVsGenJetEta_Pt60To80 = new TH1D("hSC5NGenJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.,5.0);
414	<
411	>	fSC5NGenJetsVsGenJetEta_Pt20To30 =
412	>	new TH1D("hSC5NGenJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.,5.0);
413	>	fSC5NGenJetsVsGenJetEta_Pt30To40 =
414	>	new TH1D("hSC5NGenJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.,5.0);
415	>	fSC5NGenJetsVsGenJetEta_Pt60To80 =
416	>	new TH1D("hSC5NGenJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.,5.0);
417		fSC5CaloJetsPt = new TH1D("hSC5CaloJetsPt",";SC5CaloJetPt;#",5000,0.,5000);
418		fSC5CaloJetsEta = new TH1D("hSC5CaloJetsEta",";SC5CaloJetEta;#",100,-5.0,5.0);
419	<	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt = new TH1D("hSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",";SC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
420	<	fSC5NUnmatchedCalojetsVsCorrectedCaloJetEta = new TH1D("hSC5NUnmatchedCalojetsVsCorrectedCaloJetEta",";SC5CorrectedCaloJetEta;#",100,-5.0,5.0);
421	<
422	<	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt = new TH2D("hSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",";SC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
423	<
419	>	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt =
420	>	new TH1D("hSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",
421	>	";SC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
422	>	fSC5NUnmatchedCalojetsVsCorrectedCaloJetEta =
423	>	new TH1D("hSC5NUnmatchedCalojetsVsCorrectedCaloJetEta",
424	>	";SC5CorrectedCaloJetEta;#",100,-5.0,5.0);
425	>	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt =
426	>	new TH2D("hSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",
427	>	";SC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
428		fSC5NCaloJets = new TH1D("hSC5NCaloJets",";SC5NCaloJets;#",100,0.,15);
429
430		AddOutput(fIC5GenJetRecoJetDeltaR);
#	Line 419 | Line 478 | void JetValidationMod::SlaveBegin()
478		AddOutput(fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt);
479		AddOutput(fSC5NCaloJets);
480		}
422	–
423	–	//--------------------------------------------------------------------------------------------------
424	–	void JetValidationMod::SlaveTerminate()
425	–	{
426	–	// Run finishing code on the computer (slave) that did the analysis. For this
427	–	// module, we dont do anything here.
428	–	cout << "Events Processed: " << fNEventsProcessed << endl;
429	–
430	–	}
431	–
432	–	//--------------------------------------------------------------------------------------------------
433	–	void JetValidationMod::Terminate()
434	–	{
435	–	// Run finishing code on the client computer. For this module, we dont do
436	–	// anything here.
437	–	}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines