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.4 by loizides, Mon Jun 15 15:00:23 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()));
#	Line 194 | Line 190 | void JetValidationMod::Process()
190
191		std::map <GenJet*, Jet*> SC5GenJetToSC5JetMap;
192		for (UInt_t i=0; i<GoodSC5GenJets->GetEntries(); ++i) {
193	<	double minDR = 5000;
193	>	Double_t  minDR = 5000;
194		int closestMatchIndex = -1;
195		for (UInt_t j=0; j<GoodSC5Jets->GetEntries(); j++) {
196	<	double tempDR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), GoodSC5Jets->At(j)->Mom());
196	>	Double_t  tempDR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), GoodSC5Jets->At(j)->Mom());
197		if (tempDR < minDR) {
198		closestMatchIndex = j;
199		minDR = tempDR;
200		}
201		}
202		if (closestMatchIndex > -1)
203	<	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;
203	>	if (fPrintDebug) cerr << "GoodSC5GenJets->At(i) " << GoodSC5GenJets->At(i)->Pt() << " "
204	>	<< GoodSC5GenJets->At(i)->Eta() << " " << GoodSC5GenJets->At(i)->Phi()
205	>	<< "   GoodSC5Jets->At(" << closestMatchIndex << ") "
206	>	<< GoodSC5Jets->At(closestMatchIndex)->Pt() << " "
207	>	<< GoodSC5Jets->At(closestMatchIndex)->Eta() << " "
208	>	<< GoodSC5Jets->At(closestMatchIndex)->Phi() << "    DR = "
209	>	<< minDR << endl;
210		if (minDR < 0.5)
211		{
212		SC5GenJetToSC5JetMap[GoodSC5GenJets->At(i)] = GoodSC5Jets->At(closestMatchIndex);
#	Line 212 | Line 214 | void JetValidationMod::Process()
214		}
215		std::map <Jet*, GenJet*> SC5JetToSC5GenJetMap;
216		for (UInt_t i=0; i<GoodSC5Jets->GetEntries(); ++i) {
217	<	double minDR = 5000;
217	>	Double_t  minDR = 5000;
218		int closestMatchIndex = -1;
219		for (UInt_t j=0; j<GoodSC5GenJets->GetEntries(); j++) {
220	<	double tempDR = MathUtils::DeltaR(GoodSC5Jets->At(i)->Mom(), GoodSC5GenJets->At(j)->Mom());
220	>	Double_t  tempDR = MathUtils::DeltaR(GoodSC5Jets->At(i)->Mom(), GoodSC5GenJets->At(j)->Mom());
221		if (tempDR < minDR) {
222		closestMatchIndex = j;
223		minDR = tempDR;
#	Line 238 | Line 240 | void JetValidationMod::Process()
240		map<GenJet*, Jet*>::iterator iter =
241		SC5GenJetToSC5JetMap.find(GoodSC5GenJets->At(i));
242		if (iter != SC5GenJetToSC5JetMap.end()) {
243	<	double dR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), iter->second->Mom());
243	>	Double_t  dR = MathUtils::DeltaR(GoodSC5GenJets->At(i)->Mom(), iter->second->Mom());
244
245		fSC5GenJetRecoJetDeltaR->Fill(dR);
246		fSC5GenJetRecoJetDeltaEta->Fill(fabs(GoodSC5GenJets->At(i)->Eta() - iter->second->Eta()));
247	<	fSC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodSC5GenJets->At(i)->Phi(),iter->second->Phi()));
247	>	fSC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodSC5GenJets->At(i)->Phi(),
248	>	iter->second->Phi()));
249
250		if (fabs(iter->second->Eta()) <= 1.4)
251	<	fSC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),GoodSC5GenJets->At(i)->Pt());
251	>	fSC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
252	>	GoodSC5GenJets->At(i)->Pt());
253		else if (fabs(iter->second->Eta()) <= 3.0)
254	<	fSC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),GoodSC5GenJets->At(i)->Pt());
254	>	fSC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
255	>	GoodSC5GenJets->At(i)->Pt());
256		else if (fabs(iter->second->Eta()) <= 5.0)
257	<	fSC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),GoodSC5GenJets->At(i)->Pt());
257	>	fSC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
258	>	GoodSC5GenJets->At(i)->Pt());
259		//avoid the low pt region
260		if (iter->second->Pt() > 30)
261	<	fSC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),iter->second->Eta());
261	>	fSC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodSC5GenJets->At(i)->Et(),
262	>	iter->second->Eta());
263
264		if (fabs(iter->second->Eta()) <= 3.0) {
265		fSC5CentralGenJetRecoJetDeltaR->Fill(dR);
#	Line 269 | Line 276 | void JetValidationMod::Process()
276		fSC5NMatchedCaloJetsVsGenJetEta_Pt30To40->Fill(GoodSC5GenJets->At(i)->Eta());
277		if (GoodSC5GenJets->At(i)->Pt() > 60.0 && GoodSC5GenJets->At(i)->Pt() < 80.0)
278		fSC5NMatchedCaloJetsVsGenJetEta_Pt60To80->Fill(GoodSC5GenJets->At(i)->Eta());
279	<	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(iter->second->Pt() * iter->second->L2RelativeCorrectionScale() * iter->second->L3AbsoluteCorrectionScale() / GoodSC5GenJets->At(i)->Pt(), GoodSC5GenJets->At(i)->Pt());
279	>	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(
280	>	iter->second->Pt() * iter->second->L2RelativeCorrectionScale() *
281	>	iter->second->L3AbsoluteCorrectionScale() / GoodSC5GenJets->At(i)->Pt(),
282	>	GoodSC5GenJets->At(i)->Pt());
283		}
284		}
285
#	Line 278 | Line 288 | void JetValidationMod::Process()
288		SC5JetToSC5GenJetMap.find(GoodSC5Jets->At(i));
289		//unmatched ones
290		if (iter == SC5JetToSC5GenJetMap.end()) {
291	<	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(GoodSC5Jets->At(i)->Pt() * GoodSC5Jets->At(i)->L2RelativeCorrectionScale() * GoodSC5Jets->At(i)->L3AbsoluteCorrectionScale());
291	>	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(
292	>	GoodSC5Jets->At(i)->Pt() * GoodSC5Jets->At(i)->L2RelativeCorrectionScale() *
293	>	GoodSC5Jets->At(i)->L3AbsoluteCorrectionScale());
294		fSC5NUnmatchedCalojetsVsCorrectedCaloJetEta->Fill(GoodSC5Jets->At(i)->Eta());
295		}
296		}
297
298		}
299
288	–
300		//--------------------------------------------------------------------------------------------------
301		void JetValidationMod::SlaveBegin()
302		{
#	Line 298 | Line 309 | void JetValidationMod::SlaveBegin()
309		ReqBranch(fIC5JetName,                fIC5Jets);
310		ReqBranch(fSC5JetName,                fSC5Jets);
311
312	<	fIC5GenJetRecoJetDeltaR = new TH1D("hIC5GenJetRecoJetDeltaR",";IC5GenJetRecoJetDeltaR;#",100,0.,1.0);
313	<	fIC5GenJetRecoJetDeltaEta = new TH1D("hIC5GenJetRecoJetDeltaEta",";IC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
314	<	fIC5GenJetRecoJetDeltaPhi = new TH1D("hIC5GenJetRecoJetDeltaPhi",";IC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
315	<
316	<	fIC5JetResponseVsGenJetPtInBarrel = new TH2D("hIC5JetResponseVsGenJetPtInBarrel",";IC5JetResponseVsGenJetPtInBarrel;#",500,0.,5.0,5000,0,5000);
317	<	fIC5JetResponseVsGenJetPtInEndcap = new TH2D("hIC5JetResponseVsGenJetPtInEndcap",";IC5JetResponseVsGenJetPtInEndcap;#",500,0.,5.0,5000,0,5000);
318	<	fIC5JetResponseVsGenJetPtForward = new TH2D("hIC5JetResponseVsGenJetPtForward",";IC5JetResponseVsGenJetPtForward;#",500,0.,5.0,5000,0,5000);
319	<	fIC5JetResponseVsCaloJetEta = new TH2D("hIC5JetResponseVsCaloJetEta",";IC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
320	<
321	<
322	<
323	<	fIC5CentralGenJetRecoJetDeltaR = new TH1D("hIC5CentralGenJetRecoJetDeltaR",";IC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
324	<	fIC5ForwardGenJetRecoJetDeltaR = new TH1D("hIC5ForwardGenJetRecoJetDeltaR",";IC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
325	<
326	<	fIC5NMatchedCaloJetsVsGenJetPt = new TH1D("hIC5NMatchedCaloJetsVsGenJetPt",";IC5GenJetPt;#",5000,0.,5000);
327	<	fIC5NMatchedCaloJetsVsGenJetEta = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta",";IC5GenJetEta;#",100,-5.0,5.0);
328	<	fIC5NMatchedCaloJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt20To30",";IC5GenJetEta;#",100,-5.0,5.0);
329	<	fIC5NMatchedCaloJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt30To40",";IC5GenJetEta;#",100,-5.0,5.0);
330	<	fIC5NMatchedCaloJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt60To80",";IC5GenJetEta;#",100,-5.0,5.0);
312	>	fIC5GenJetRecoJetDeltaR =
313	>	new TH1D("hIC5GenJetRecoJetDeltaR",";IC5GenJetRecoJetDeltaR;#",100,0.,1.0);
314	>	fIC5GenJetRecoJetDeltaEta =
315	>	new TH1D("hIC5GenJetRecoJetDeltaEta",";IC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
316	>	fIC5GenJetRecoJetDeltaPhi =
317	>	new TH1D("hIC5GenJetRecoJetDeltaPhi",";IC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
318	>
319	>	fIC5JetResponseVsGenJetPtInBarrel = new TH2D("hIC5JetResponseVsGenJetPtInBarrel",
320	>	";IC5JetResponseVsGenJetPtInBarrel;#",
321	>	500,0.,5.0,5000,0,5000);
322	>	fIC5JetResponseVsGenJetPtInEndcap = new TH2D("hIC5JetResponseVsGenJetPtInEndcap",
323	>	";IC5JetResponseVsGenJetPtInEndcap;#",
324	>	500,0.,5.0,5000,0,5000);
325	>	fIC5JetResponseVsGenJetPtForward = new TH2D("hIC5JetResponseVsGenJetPtForward",
326	>	";IC5JetResponseVsGenJetPtForward;#",
327	>	500,0.,5.0,5000,0,5000);
328	>	fIC5JetResponseVsCaloJetEta = new TH2D("hIC5JetResponseVsCaloJetEta",
329	>	";IC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
330	>	fIC5CentralGenJetRecoJetDeltaR =
331	>	new TH1D("hIC5CentralGenJetRecoJetDeltaR",";IC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
332	>	fIC5ForwardGenJetRecoJetDeltaR = new TH1D("hIC5ForwardGenJetRecoJetDeltaR",
333	>	";IC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
334	>	fIC5NMatchedCaloJetsVsGenJetPt = new TH1D("hIC5NMatchedCaloJetsVsGenJetPt",
335	>	";IC5GenJetPt;#",5000,0.,5000);
336	>	fIC5NMatchedCaloJetsVsGenJetEta = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta",
337	>	";IC5GenJetEta;#",100,-5.0,5.0);
338	>	fIC5NMatchedCaloJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt20To30",
339	>	";IC5GenJetEta;#",100,-5.0,5.0);
340	>	fIC5NMatchedCaloJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt30To40",
341	>	";IC5GenJetEta;#",100,-5.0,5.0);
342	>	fIC5NMatchedCaloJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NMatchedCaloJetsVsGenJetEta_Pt60To80",
343	>	";IC5GenJetEta;#",100,-5.0,5.0);
344		fIC5NGenJetsVsGenJetPt = new TH1D("hIC5NGenJetsVsGenJetPt",";IC5GenJetPt;#",5000,0.,5000);
345		fIC5NGenJetsVsGenJetEta = new TH1D("hIC5NGenJetsVsGenJetEta",";IC5GenJetEta;#",100,-5.,5.0);
346	<	fIC5NGenJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt20To30",";IC5GenJetEta;#",100,-5.,5.0);
347	<	fIC5NGenJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt30To40",";IC5GenJetEta;#",100,-5.,5.0);
348	<	fIC5NGenJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt60To80",";IC5GenJetEta;#",100,-5.,5.0);
349	<
346	>	fIC5NGenJetsVsGenJetEta_Pt20To30 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt20To30",
347	>	";IC5GenJetEta;#",100,-5.,5.0);
348	>	fIC5NGenJetsVsGenJetEta_Pt30To40 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt30To40",
349	>	";IC5GenJetEta;#",100,-5.,5.0);
350	>	fIC5NGenJetsVsGenJetEta_Pt60To80 = new TH1D("hIC5NGenJetsVsGenJetEta_Pt60To80",
351	>	";IC5GenJetEta;#",100,-5.,5.0);
352		fIC5CaloJetsPt = new TH1D("hIC5CaloJetsPt",";IC5CaloJetPt;#",5000,0.,5000);
353		fIC5CaloJetsEta = new TH1D("hIC5CaloJetsEta",";IC5CaloJetEta;#",100,-5.0,5.0);
354	<	fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt = new TH1D("hIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",";IC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
355	<	fIC5NUnmatchedCalojetsVsCorrectedCaloJetEta = new TH1D("hIC5NUnmatchedCalojetsVsCorrectedCaloJetEta",";IC5CorrectedCaloJetEta;#",100,-5.0,5.0);
356	<
357	<	fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt = new TH2D("hIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",";IC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
358	<
354	>	fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt =
355	>	new TH1D("hIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",
356	>	";IC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
357	>	fIC5NUnmatchedCalojetsVsCorrectedCaloJetEta =
358	>	new TH1D("hIC5NUnmatchedCalojetsVsCorrectedCaloJetEta",
359	>	";IC5CorrectedCaloJetEta;#",100,-5.0,5.0);
360	>	fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt =
361	>	new TH2D("hIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",
362	>	";IC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
363		fIC5NCaloJets = new TH1D("hIC5NCaloJets",";IC5NCaloJets;#",100,0.,15);
364	<
365	<
366	<	fSC5GenJetRecoJetDeltaR = new TH1D("hSC5GenJetRecoJetDeltaR",";SC5GenJetRecoJetDeltaR;#",100,0.,1.0);
367	<	fSC5GenJetRecoJetDeltaEta = new TH1D("hSC5GenJetRecoJetDeltaEta",";SC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
368	<	fSC5GenJetRecoJetDeltaPhi = new TH1D("hSC5GenJetRecoJetDeltaPhi",";SC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
369	<
370	<	fSC5JetResponseVsGenJetPtInBarrel = new TH2D("hSC5JetResponseVsGenJetPtInBarrel",";SC5JetResponseVsGenJetPtInBarrel;#",500,0.,5.0,5000,0,5000);
371	<	fSC5JetResponseVsGenJetPtInEndcap = new TH2D("hSC5JetResponseVsGenJetPtInEndcap",";SC5JetResponseVsGenJetPtInEndcap;#",500,0.,5.0,5000,0,5000);
372	<	fSC5JetResponseVsGenJetPtForward = new TH2D("hSC5JetResponseVsGenJetPtForward",";SC5JetResponseVsGenJetPtForward;#",500,0.,5.0,5000,0,5000);
373	<
374	<	fSC5JetResponseVsCaloJetEta = new TH2D("hSC5JetResponseVsCaloJetEta",";SC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
375	<
376	<
377	<
378	<	fSC5CentralGenJetRecoJetDeltaR = new TH1D("hSC5CentralGenJetRecoJetDeltaR",";SC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
379	<	fSC5ForwardGenJetRecoJetDeltaR = new TH1D("hSC5ForwardGenJetRecoJetDeltaR",";SC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
380	<
381	<	fSC5NMatchedCaloJetsVsGenJetPt = new TH1D("hSC5NMatchedCaloJetsVsGenJetPt",";SC5GenJetPt;#",5000,0.,5000);
382	<	fSC5NMatchedCaloJetsVsGenJetEta = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta",";SC5GenJetEta;#",100,-5.0,5.0);
383	<	fSC5NMatchedCaloJetsVsGenJetEta_Pt20To30 = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.0,5.0);
384	<	fSC5NMatchedCaloJetsVsGenJetEta_Pt30To40 = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.0,5.0);
385	<	fSC5NMatchedCaloJetsVsGenJetEta_Pt60To80 = new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.0,5.0);
364	>	fSC5GenJetRecoJetDeltaR = new TH1D("hSC5GenJetRecoJetDeltaR",
365	>	";SC5GenJetRecoJetDeltaR;#",100,0.,1.0);
366	>	fSC5GenJetRecoJetDeltaEta = new TH1D("hSC5GenJetRecoJetDeltaEta",
367	>	";SC5GenJetRecoJetDeltaEta;#",100,0.,0.5);
368	>	fSC5GenJetRecoJetDeltaPhi = new TH1D("hSC5GenJetRecoJetDeltaPhi",
369	>	";SC5GenJetRecoJetDeltaPhi;#",100,0.,0.5);
370	>	fSC5JetResponseVsGenJetPtInBarrel =
371	>	new TH2D("hSC5JetResponseVsGenJetPtInBarrel",
372	>	";SC5JetResponseVsGenJetPtInBarrel;#",500,0.,5.0,5000,0,5000);
373	>	fSC5JetResponseVsGenJetPtInEndcap =
374	>	new TH2D("hSC5JetResponseVsGenJetPtInEndcap",
375	>	";SC5JetResponseVsGenJetPtInEndcap;#",500,0.,5.0,5000,0,5000);
376	>	fSC5JetResponseVsGenJetPtForward =
377	>	new TH2D("hSC5JetResponseVsGenJetPtForward",
378	>	";SC5JetResponseVsGenJetPtForward;#",500,0.,5.0,5000,0,5000);
379	>
380	>	fSC5JetResponseVsCaloJetEta = new TH2D("hSC5JetResponseVsCaloJetEta",
381	>	";SC5JetResponseVsCaloJetEta;#",500,0.,5.0,100,-5,5);
382	>	fSC5CentralGenJetRecoJetDeltaR =
383	>	new TH1D("hSC5CentralGenJetRecoJetDeltaR",
384	>	";SC5CentralGenJetRecoJetDeltaR;#",100,0.,1.0);
385	>	fSC5ForwardGenJetRecoJetDeltaR =
386	>	new TH1D("hSC5ForwardGenJetRecoJetDeltaR",
387	>	";SC5ForwardGenJetRecoJetDeltaR;#",100,0.,1.0);
388	>	fSC5NMatchedCaloJetsVsGenJetPt =
389	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetPt",";SC5GenJetPt;#",5000,0.,5000);
390	>	fSC5NMatchedCaloJetsVsGenJetEta =
391	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta",";SC5GenJetEta;#",100,-5.0,5.0);
392	>	fSC5NMatchedCaloJetsVsGenJetEta_Pt20To30 =
393	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.0,5.0);
394	>	fSC5NMatchedCaloJetsVsGenJetEta_Pt30To40 =
395	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.0,5.0);
396	>	fSC5NMatchedCaloJetsVsGenJetEta_Pt60To80 =
397	>	new TH1D("hSC5NMatchedCaloJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.0,5.0);
398		fSC5NGenJetsVsGenJetPt = new TH1D("hSC5NGenJetsVsGenJetPt",";SC5GenJetPt;#",5000,0.,5000);
399		fSC5NGenJetsVsGenJetEta = new TH1D("hSC5NGenJetsVsGenJetEta",";SC5GenJetEta;#",100,-5.,5.0);
400	<	fSC5NGenJetsVsGenJetEta_Pt20To30 = new TH1D("hSC5NGenJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.,5.0);
401	<	fSC5NGenJetsVsGenJetEta_Pt30To40 = new TH1D("hSC5NGenJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.,5.0);
402	<	fSC5NGenJetsVsGenJetEta_Pt60To80 = new TH1D("hSC5NGenJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.,5.0);
403	<
400	>	fSC5NGenJetsVsGenJetEta_Pt20To30 =
401	>	new TH1D("hSC5NGenJetsVsGenJetEta_Pt20To30",";SC5GenJetEta;#",100,-5.,5.0);
402	>	fSC5NGenJetsVsGenJetEta_Pt30To40 =
403	>	new TH1D("hSC5NGenJetsVsGenJetEta_Pt30To40",";SC5GenJetEta;#",100,-5.,5.0);
404	>	fSC5NGenJetsVsGenJetEta_Pt60To80 =
405	>	new TH1D("hSC5NGenJetsVsGenJetEta_Pt60To80",";SC5GenJetEta;#",100,-5.,5.0);
406		fSC5CaloJetsPt = new TH1D("hSC5CaloJetsPt",";SC5CaloJetPt;#",5000,0.,5000);
407		fSC5CaloJetsEta = new TH1D("hSC5CaloJetsEta",";SC5CaloJetEta;#",100,-5.0,5.0);
408	<	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt = new TH1D("hSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",";SC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
409	<	fSC5NUnmatchedCalojetsVsCorrectedCaloJetEta = new TH1D("hSC5NUnmatchedCalojetsVsCorrectedCaloJetEta",";SC5CorrectedCaloJetEta;#",100,-5.0,5.0);
410	<
411	<	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt = new TH2D("hSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",";SC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
412	<
408	>	fSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt =
409	>	new TH1D("hSC5NUnmatchedCaloJetsVsCorrectedCaloJetPt",
410	>	";SC5CorrectedCaloJetPt;#",5000,0.0,5000.0);
411	>	fSC5NUnmatchedCalojetsVsCorrectedCaloJetEta =
412	>	new TH1D("hSC5NUnmatchedCalojetsVsCorrectedCaloJetEta",
413	>	";SC5CorrectedCaloJetEta;#",100,-5.0,5.0);
414	>	fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt =
415	>	new TH2D("hSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt",
416	>	";SC5GenJetPt;#",1000,0,1000.0, 100,0.0,5.0);
417		fSC5NCaloJets = new TH1D("hSC5NCaloJets",";SC5NCaloJets;#",100,0.,15);
418
419		AddOutput(fIC5GenJetRecoJetDeltaR);
#	Line 419 | Line 467 | void JetValidationMod::SlaveBegin()
467		AddOutput(fSC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt);
468		AddOutput(fSC5NCaloJets);
469		}
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