ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitPhysics/Utils/src/JetIDMVA.cc

Comparing UserCode/MitPhysics/Utils/src/JetIDMVA.cc (file contents):
Revision 1.3 by pharris, Wed Apr 4 09:51:31 2012 UTC vs.
Revision 1.10 by ceballos, Sat May 12 07:26:56 2012 UTC

#	Line 1 | Line 1
1	+	#include "FWCore/PythonParameterSet/interface/MakeParameterSets.h"
2	+	#include "FWCore/ParameterSet/interface/ParameterSet.h"
3		#include "MitPhysics/Utils/interface/JetIDMVA.h"
4		#include "MitPhysics/Utils/interface/JetTools.h"
5		#include "MitAna/DataTree/interface/StableData.h"
#	Line 13 | Line 15 | using namespace mithep;
15
16		//--------------------------------------------------------------------------------------------------
17		JetIDMVA::JetIDMVA() :
18	<	fJetPtMin(10), //We need to lower this
19	<	fMethodName("JetIDMA"),
18	>	fJetPtMin(0.)  , //We need to lower this
19	>	fDZCut   (0.2),
20	>	fLowPtMethodName ("JetIDMVALowPt" ),
21	>	fHighPtMethodName("JetIDMVAHighPt"),
22		fIsInitialized(kFALSE),
23	<	fNPV    (0),
24	<	fJPt1   (0),
25	<	fJEta1  (0),
26	<	fJPhi1  (0),
27	<	fJD01   (0),
28	<	fJDZ1   (0),
29	<	fJM1    (0),
30	<	fNPart1 (0),
31	<	fLPt1   (0),
32	<	fLEta1  (0),
33	<	fLPhi1  (0),
34	<	fSPt1   (0),
35	<	fSEta1  (0),
36	<	fSPhi1  (0),
37	<	fNEPt1  (0),
38	<	fNEEta1 (0),
35	<	fNEPhi1 (0),
36	<	fEMPt1  (0),
37	<	fEMEta1 (0),
38	<	fEMPhi1 (0),
39	<	fChPt1  (0),
40	<	fChPhi1 (0),
41	<	fLFr1   (0),
42	<	fDRLC1  (0),
43	<	fDRLS1  (0),
44	<	fDRM1   (0),
45	<	fDRNE1  (0),
46	<	fDREM1  (0),
47	<	fDRCH1  (0)
23	>	fNVtx     (0),
24	>	fJPt1     (0),
25	>	fJEta1    (0),
26	>	fJPhi1    (0),
27	>	fJD01     (0),
28	>	fJDZ1     (0),
29	>	fBeta     (0),
30	>	fBetaStar (0),
31	>	fNCharged (0),
32	>	fNNeutrals(0),
33	>	fDRMean   (0),
34	>	fFrac01   (0),
35	>	fFrac02   (0),
36	>	fFrac03   (0),
37	>	fFrac04   (0),
38	>	fFrac05   (0)
39		{
40		fReader = 0;
41		}
42		//--------------------------------------------------------------------------------------------------
43	<	JetIDMVA::~JetIDMVA()
44	<	{
43	>	JetIDMVA::~JetIDMVA() {
44	>
45		fReader = 0;
46		}
47
48		//--------------------------------------------------------------------------------------------------
49	<	void JetIDMVA::Initialize( TString iMethodName,
50	<	TString iWeights,
51	<	JetIDMVA::MVAType iType)
49	>	void JetIDMVA::Initialize( JetIDMVA::CutType iCutType,
50	>	TString iLowPtWeights,
51	>	TString iHighPtWeights,
52	>	JetIDMVA::MVAType iType,
53	>	TString iCutFileName)
54		{
55
56		fIsInitialized = kTRUE;
64	–	fMethodName    = iMethodName;
57		fType          = iType;
58	+	fCutType       = iCutType;
59		fReader        = 0;
60		fReader        = new TMVA::Reader( "!Color:!Silent:Error" );
61		if (fType == kBaseline) {
62	<	fReader->AddSpectator( "npv"     , &fNPV    );
63	<	fReader->AddVariable( "jspt_1"   , &fJPt1   );
64	<	fReader->AddVariable( "jseta_1"  , &fJEta1  );
65	<	fReader->AddVariable( "jsphi_1"  , &fJPhi1  );
66	<	fReader->AddVariable( "jd0_1"    , &fJD01   );
67	<	fReader->AddVariable( "jdZ_1"    , &fJDZ1   );
68	<	fReader->AddVariable( "jm_1"     , &fJM1    );
69	<	fReader->AddVariable( "npart_1"  , &fNPart1 );
70	<	fReader->AddVariable( "lpt_1"    , &fLPt1   );
71	<	fReader->AddVariable( "leta_1"   , &fLEta1  );
72	<	fReader->AddVariable( "lphi_1"   , &fLPhi1  );
73	<	fReader->AddVariable( "spt_1"    , &fSPt1   );
74	<	fReader->AddVariable( "seta_1"   , &fSEta1  );
75	<	fReader->AddVariable( "sphi_1"   , &fSPhi1  );
76	<	fReader->AddVariable( "lnept_1"  , &fNEPt1  );
77	<	fReader->AddVariable( "lneeta_1" , &fNEEta1 );
85	<	fReader->AddVariable( "lnephi_1" , &fNEPhi1 );
86	<	fReader->AddVariable( "lempt_1"  , &fEMPt1  );
87	<	fReader->AddVariable( "lemeta_1" , &fEMEta1 );
88	<	fReader->AddVariable( "lemphi_1" , &fEMPhi1 );
89	<	fReader->AddVariable( "lchpt_1"  , &fChPt1  );
90	<	fReader->AddVariable( "lchphi_1" , &fChPhi1 );
91	<	fReader->AddVariable( "lLfr_1"   , &fLFr1   );
92	<	fReader->AddVariable( "drlc_1"   , &fDRLC1  );
93	<	fReader->AddVariable( "drls_1"   , &fDRLS1  );
94	<	fReader->AddVariable( "drm_1"    , &fDRM1   );
95	<	fReader->AddVariable( "drmne_1"  , &fDRNE1  );
96	<	fReader->AddVariable( "drem_1"   , &fDREM1  );
97	<	fReader->AddVariable( "drch_1"   , &fDRCH1  );
62	>	fReader->AddVariable( "nvtx"     , &fNVtx      );
63	>	fReader->AddVariable( "jetPt"    , &fJPt1      );
64	>	fReader->AddVariable( "jetEta"   , &fJEta1     );
65	>	fReader->AddVariable( "jetPhi"   , &fJPhi1     );
66	>	fReader->AddVariable( "dZ"       , &fJDZ1      );
67	>	fReader->AddVariable( "d0"       , &fJD01      );
68	>	fReader->AddVariable( "beta"     , &fBeta      );
69	>	fReader->AddVariable( "betaStar" , &fBetaStar  );
70	>	fReader->AddVariable( "nCharged" , &fNCharged  );
71	>	fReader->AddVariable( "nNeutrals", &fNNeutrals );
72	>	fReader->AddVariable( "dRMean"   , &fDRMean    );
73	>	fReader->AddVariable( "frac01"   , &fFrac01    );
74	>	fReader->AddVariable( "frac02"   , &fFrac02    );
75	>	fReader->AddVariable( "frac03"   , &fFrac03    );
76	>	fReader->AddVariable( "frac04"   , &fFrac04    );
77	>	fReader->AddVariable( "frac05"   , &fFrac05    );
78		}
79	<	fReader->BookMVA(fMethodName , iWeights );
79	>	fReader->BookMVA(fLowPtMethodName  , iLowPtWeights );
80	>	fReader->BookMVA(fHighPtMethodName , iHighPtWeights );
81		std::cout << "Jet ID MVA Initialization\n";
82	<	std::cout << "MethodName : " << fMethodName << " , type == " << fType << std::endl;
82	>	std::cout << "MethodName : " << fLowPtMethodName << " , type == " << fType << std::endl;
83	>
84	>	//Load Cut Matrix
85	>	edm::ParameterSet lConfig = edm::readPSetsFrom(iCutFileName.Data())->getParameter<edm::ParameterSet>("JetIdParams");
86	>	std::string lCutType = "Tight";
87	>	if(fCutType == kMedium) lCutType = "Medium";
88	>	if(fCutType == kLoose ) lCutType = "Loose";
89	>	if(fCutType == kMET   ) lCutType = "MET";
90	>	std::vector<double> lPt010  = lConfig.getParameter<std::vector<double> >(("Pt010_" +lCutType).c_str());
91	>	std::vector<double> lPt1020 = lConfig.getParameter<std::vector<double> >(("Pt1020_"+lCutType).c_str());
92	>	std::vector<double> lPt2030 = lConfig.getParameter<std::vector<double> >(("Pt2030_"+lCutType).c_str());
93	>	std::vector<double> lPt3050 = lConfig.getParameter<std::vector<double> >(("Pt3050_"+lCutType).c_str());
94	>	for(int i0 = 0; i0 < 4; i0++) fMVACut[0][i0] = lPt010 [i0];
95	>	for(int i0 = 0; i0 < 4; i0++) fMVACut[1][i0] = lPt1020[i0];
96	>	for(int i0 = 0; i0 < 4; i0++) fMVACut[2][i0] = lPt2030[i0];
97	>	for(int i0 = 0; i0 < 4; i0++) fMVACut[3][i0] = lPt3050[i0];
98	>	//std::cout << " Working Points : << " << std::endl;
99	>	//for(int i0 = 0; i0 < 4; i0++) for(int i1 = 0; i1 < 4; i1++)
100	>	//  std::cout << " ==> " << i0 << " -- " << i1 << " -- " << fMVACut[i0][i1] << std::endl;
101		}
102
103		//--------------------------------------------------------------------------------------------------
#	Line 109 | Line 108 | Double_t JetIDMVA::MVAValue(
108		Float_t iJPhi1  ,
109		Float_t iJD01   ,
110		Float_t iJDZ1   ,
111	<	Float_t iJM1    ,
112	<	Float_t iNPart1 ,
113	<	Float_t iLPt1   ,
114	<	Float_t iLEta1  ,
115	<	Float_t iLPhi1  ,
116	<	Float_t iSPt1   ,
117	<	Float_t iSEta1  ,
118	<	Float_t iSPhi1  ,
119	<	Float_t iNEPt1  ,
120	<	Float_t iNEEta1 ,
122	<	Float_t iNEPhi1 ,
123	<	Float_t iEMPt1  ,
124	<	Float_t iEMEta1 ,
125	<	Float_t iEMPhi1 ,
126	<	Float_t iChPt1  ,
127	<	Float_t iChPhi1 ,
128	<	Float_t iLFr1   ,
129	<	Float_t iDRLC1  ,
130	<	Float_t iDRLS1  ,
131	<	Float_t iDRM1   ,
132	<	Float_t iDRNE1 ,
133	<	Float_t iDREM1  ,
134	<	Float_t iDRCH1
111	>	Float_t iBeta   ,
112	>	Float_t iBetaStar,
113	>	Float_t iNCharged,
114	>	Float_t iNNeutrals,
115	>	Float_t iDRMean  ,
116	>	Float_t iFrac01  ,
117	>	Float_t iFrac02  ,
118	>	Float_t iFrac03  ,
119	>	Float_t iFrac04  ,
120	>	Float_t iFrac05
121		){
122
123	<	if (!fIsInitialized) {
123	>	if(!fIsInitialized) {
124		std::cout << "Error: JetIDMVA not properly initialized.\n";
125		return -9999;
126		}
127
128	<	fNPV    = iNPV;
129	<	fJPt1   = iJPt1;
130	<	fJEta1  = iJEta1;
131	<	fJPhi1  = fJPhi1;
132	<	fJD01   = iJD01;
133	<	fJDZ1   = iJDZ1;
134	<	fJM1    = iJM1 ;
135	<	fNPart1 = iNPart1;
136	<	fLPt1   = iLPt1;
137	<	fLEta1  = iLEta1;
138	<	fLPhi1  = iLPhi1;
139	<	fSPt1   = iSPt1;
140	<	fSEta1  = iSEta1;
141	<	fSPhi1  = iSPhi1;
142	<	fNEPt1  = iNEPt1;
143	<	fNEEta1 = iNEEta1;
158	<	fNEPhi1 = iNEPhi1;
159	<	fEMPt1  = iEMPt1;
160	<	fEMEta1 = iEMEta1;
161	<	fEMPhi1 = iEMPhi1;
162	<	fChPt1  = iChPt1;
163	<	fChPhi1 = iChPhi1;
164	<	fLFr1   = iLFr1;
165	<	fDRLC1  = iDRLC1;
166	<	fDRLS1  = iDRLS1;
167	<	fDRM1   = iDRM1;
168	<	fDRNE1  = iDRNE1;
169	<	fDREM1  = iDREM1;
170	<	fDRCH1  = iDRCH1;
128	>	fNVtx      = iNPV;
129	>	fJPt1      = iJPt1;
130	>	fJEta1     = iJEta1;
131	>	fJPhi1     = fJPhi1;
132	>	fJD01      = iJD01;
133	>	fJDZ1      = iJDZ1;
134	>	fBeta      = iBeta;
135	>	fBetaStar  = iBetaStar;
136	>	fNCharged  = iNCharged;
137	>	fNNeutrals = iNNeutrals;
138	>	fDRMean    = iDRMean;
139	>	fFrac01    = iFrac01;
140	>	fFrac02    = iFrac02;
141	>	fFrac03    = iFrac03;
142	>	fFrac04    = iFrac04;
143	>	fFrac05    = iFrac05;
144
145		Double_t lMVA = -9999;
146	<	lMVA = fReader->EvaluateMVA( fMethodName );
146	>	if(iJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
147	>	if(iJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
148
149		return lMVA;
150		}
151		//--------------------------------------------------------------------------------------------------
152	<	Bool_t JetIDMVA::pass(const PFJet *iJet,const Vertex *iVertex,
152	>	Bool_t JetIDMVA::pass(const PFJet *iJet,const Vertex *iVertex,const VertexCol *iVertices,
153		FactorizedJetCorrector *iJetCorrector,
154	<	const PileupEnergyDensityCol *iPileupEnergyDensity) {
154	>	const PileupEnergyDensityCol *iPileupEnergyDensity,
155	>	RhoUtilities::RhoType type) {
156	>
157		if(!JetTools::passPFLooseId(iJet))                 return false;
158	<	if(correctedPt(iJet,iJetCorrector,iPileupEnergyDensity) < fJetPtMin && iJet->TrackCountingHighEffBJetTagsDisc() == -100) return false; //This line is a bug in the Met training
159	<	if( fabs(JetTools::impactParameter(iJet,iVertex,true)) < 0.2) return true;
160	<	double lMVA = MVAValue(iJet,iVertex,iJetCorrector,iPileupEnergyDensity);
161	<	if(lMVA < -0.8)                            return false;
162	<	if(lMVA < -0.5 && fabs(iJet->Eta()) > 3.0) return false;
158	>	if(fabs(iJet->Eta()) > 4.99)                       return false;
159	>
160	>	double lMVA = MVAValue   (iJet,iVertex,iVertices,iJetCorrector,iPileupEnergyDensity);
161	>	double lPt  = correctedPt(iJet,                  iJetCorrector,iPileupEnergyDensity,type);
162	>	if(lPt < fJetPtMin)                         return false;
163	>
164	>	int lPtId = 0;
165	>	if(lPt > 10 && lPt < 20) lPtId = 1;
166	>	if(lPt > 20 && lPt < 30) lPtId = 2;
167	>	if(lPt > 30                   ) lPtId = 3;
168	>
169	>	int lEtaId = 0;
170	>	if(fabs(iJet->Eta()) > 2.5  && fabs(iJet->Eta()) < 2.75) lEtaId = 1;
171	>	if(fabs(iJet->Eta()) > 2.75 && fabs(iJet->Eta()) < 3.0 ) lEtaId = 2;
172	>	if(fabs(iJet->Eta()) > 3.0  && fabs(iJet->Eta()) < 5.0 ) lEtaId = 3;
173	>
174	>	double lMVACut = fMVACut[lPtId][lEtaId];
175	>	if(lMVA < lMVACut) return false;
176		return true;
177		}
178		//--------------------------------------------------------------------------------------------------
179	<	Bool_t JetIDMVA::pass(const PFJet *iJet,const Vertex *iVertex) {
179	>	Bool_t JetIDMVA::pass(const PFJet *iJet,const Vertex *iVertex,const VertexCol *iVertices) {
180		if(!JetTools::passPFLooseId(iJet))                 return false;
181	<	if(iJet->Pt() < fJetPtMin && iJet->TrackCountingHighEffBJetTagsDisc() == -100) return false; //This line is a bug in the Met training
182	<	if( fabs(JetTools::impactParameter(iJet,iVertex,true)) < 0.2) return true;
183	<	double lMVA = MVAValue(iJet,iVertex);
184	<	if(lMVA < -0.8)                            return false;
185	<	if(lMVA < -0.5 && fabs(iJet->Eta()) > 3.0) return false;
181	>	if(iJet->Pt()        < fJetPtMin) return false;
182	>	if(fabs(iJet->Eta()) > 4.99)      return false;
183	>	double lMVA = MVAValue(iJet,iVertex,iVertices);
184	>
185	>	int lPtId = 0;
186	>	if(iJet->Pt() > 10 && iJet->Pt() < 20) lPtId = 1;
187	>	if(iJet->Pt() > 20 && iJet->Pt() < 30) lPtId = 2;
188	>	if(iJet->Pt() > 30                   ) lPtId = 3;
189	>
190	>	int lEtaId = 0;
191	>	if(fabs(iJet->Eta()) > 2.5  && fabs(iJet->Eta()) < 2.75) lEtaId = 1;
192	>	if(fabs(iJet->Eta()) > 2.75 && fabs(iJet->Eta()) < 3.0 ) lEtaId = 2;
193	>	if(fabs(iJet->Eta()) > 3.0  && fabs(iJet->Eta()) < 5.0 ) lEtaId = 3;
194	>
195	>	double lMVACut = fMVACut[lPtId][lEtaId];
196	>	if(lMVA < lMVACut) return false;
197		return true;
198	+	//if(lMVA < -0.8)                            return false;
199	+	//if(lMVA < -0.5 && fabs(iJet->Eta()) > 3.0) return false;
200	+	//if(iJet->Pt() < fJetPtMin && iJet->TrackCountingHighEffBJetTagsDisc() == -100) return false; //This line is a bug in the Met training
201	+	//if( fabs(JetTools::impactParameter(iJet,iVertex,true)) < 0.2) return true;
202		}
203		//--------------------------------------------------------------------------------------------------
204	<	Double_t JetIDMVA::MVAValue(const PFJet *iJet,const Vertex *iVertex, //Vertex here is the PV
204	>	Double_t JetIDMVA::MVAValue(const PFJet *iJet,const Vertex *iVertex,const VertexCol *iVertices, //Vertex here is the PV
205		FactorizedJetCorrector *iJetCorrector,
206		const PileupEnergyDensityCol *iPileupEnergyDensity,
207		Bool_t printDebug) {
#	Line 209 | Line 213 | Double_t JetIDMVA::MVAValue(const PFJet
213		if(!JetTools::passPFLooseId(iJet)) return -2.;
214
215		//set all input variables
216	+	fNVtx      = iVertices->GetEntries();
217		fJPt1      = correctedPt(iJet,iJetCorrector,iPileupEnergyDensity);
218		fJEta1     = iJet->RawMom().Eta();
219		fJPhi1     = iJet->RawMom().Phi();
220	<	fJM1       = iJet->Mass();
221	<
222	<	const mithep::PFCandidate *lLead     = JetTools::leadCand(iJet,-1);
223	<	const mithep::PFCandidate *lSecond   = JetTools::leadCand(iJet,-1,true);
224	<	const mithep::PFCandidate *lLeadNeut = JetTools::leadCand(iJet ,5);
225	<	const mithep::PFCandidate *lLeadEm   = JetTools::leadCand(iJet ,4);
226	<	const mithep::PFCandidate *lLeadCh   = JetTools::leadCand(iJet ,1);
227	<
228	<	fJD01         = JetTools::impactParameter(iJet,iVertex);
229	<	fJDZ1         = JetTools::impactParameter(iJet,iVertex,true);
230	<	fNPart1       = iJet->NPFCands();
231	<	fLPt1         = lLead    ->Pt();
232	<	fLEta1        = lLead    ->Eta();
233	<	fLPhi1        = lLead    ->Phi();
234	<	fSPt1         = lSecond  ->Pt();
235	<	fSEta1        = lSecond  ->Eta();
236	<	fSPhi1        = lSecond  ->Phi();
232	<	fNEPt1        = lLeadNeut->Pt();
233	<	fNEEta1       = lLeadNeut->Eta();
234	<	fNEPhi1       = lLeadNeut->Phi();
235	<	fEMPt1        = lLeadEm  ->Pt();
236	<	fEMEta1       = lLeadEm  ->Eta();
237	<	fEMPhi1       = lLeadEm  ->Phi();
238	<	fChPt1        = lLeadCh  ->Pt();
239	<	//fChEta1       = lLeadCh  ->Eta();
240	<	fChPhi1       = lLeadCh  ->Phi();
241	<	fLFr1         = lLead->Pt()/iJet->Pt();
242	<
243	<	fDRLC1        = MathUtils::DeltaR(iJet->Mom(),lLead  ->Mom());
244	<	fDRLS1        = MathUtils::DeltaR(iJet->Mom(),lSecond->Mom());
245	<	fDRM1         = JetTools::dRMean (iJet,-1);
246	<	fDRNE1        = JetTools::dRMean (iJet, 5);
247	<	fDREM1        = JetTools::dRMean (iJet, 4);
248	<	fDRCH1        = JetTools::dRMean (iJet, 1);
249	<
250	<	double lMVA = fReader->EvaluateMVA( fMethodName );
251	<
220	>	fJD01      = JetTools::impactParameter(iJet,iVertex);
221	>	fJDZ1      = JetTools::impactParameter(iJet,iVertex,true);
222	>	fBeta      = JetTools::Beta(iJet,iVertex,fDZCut);
223	>	fBetaStar  = JetTools::betaStar(iJet,iVertex,iVertices,fDZCut);
224	>	fNCharged  = iJet->ChargedMultiplicity();
225	>	fNNeutrals = iJet->NeutralMultiplicity();
226	>
227	>	fDRMean    = JetTools::dRMean(iJet,-1);
228	>	fFrac01    = JetTools::frac  (iJet,0.1,0. ,-1);
229	>	fFrac02    = JetTools::frac  (iJet,0.2,0.1,-1);
230	>	fFrac03    = JetTools::frac  (iJet,0.3,0.2,-1);
231	>	fFrac04    = JetTools::frac  (iJet,0.4,0.3,-1);
232	>	fFrac05    = JetTools::frac  (iJet,0.5,0.4,-1);
233	>
234	>	double lMVA = 0;
235	>	if(fJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
236	>	if(fJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
237		if (printDebug == kTRUE) {
238		std::cout << "Debug Jet MVA: "
239	<	<< fNPV    << " "
240	<	<< fJPt1   << " "
241	<	<< fJEta1  << " "
242	<	<< fJPhi1  << " "
243	<	<< fJD01   << " "
244	<	<< fJDZ1   << " "
245	<	<< fJM1    << " "
246	<	<< fNPart1 << " "
247	<	<< fLPt1   << " "
248	<	<< fLEta1  << " "
249	<	<< fLPhi1  << " "
250	<	<< fSPt1   << " "
251	<	<< fSEta1  << " "
252	<	<< fSPhi1  << " "
253	<	<< fNEPt1  << " "
254	<	<< fNEEta1 << " "
270	<	<< fNEPhi1 << " "
271	<	<< fEMPt1  << " "
272	<	<< fEMEta1 << " "
273	<	<< fEMPhi1 << " "
274	<	<< fChPt1  << " "
275	<	<< fChPhi1 << " "
276	<	<< fLFr1   << " "
277	<	<< fDRLC1  << " "
278	<	<< fDRLS1  << " "
279	<	<< fDRM1   << " "
280	<	<< fDRNE1 << " "
281	<	<< fDREM1  << " "
282	<	<< fDRCH1  << " "
239	>	<< fNVtx      << " "
240	>	<< fJPt1      << " "
241	>	<< fJEta1     << " "
242	>	<< fJPhi1     << " "
243	>	<< fJD01      << " "
244	>	<< fJDZ1      << " "
245	>	<< fBeta      << " "
246	>	<< fBetaStar  << " "
247	>	<< fNCharged  << " "
248	>	<< fNNeutrals << " "
249	>	<< fDRMean    << " "
250	>	<< fFrac01    << " "
251	>	<< fFrac02    << " "
252	>	<< fFrac03    << " "
253	>	<< fFrac04    << " "
254	>	<< fFrac05
255		<< " === : === "
256		<< lMVA << " "
257		<< std::endl;
#	Line 287 | Line 259 | Double_t JetIDMVA::MVAValue(const PFJet
259
260		return lMVA;
261		}
262	<	Double_t JetIDMVA::MVAValue(const PFJet *iJet,const Vertex *iVertex, //Vertex here is the PV
262	>	Double_t JetIDMVA::MVAValue(const PFJet *iJet,const Vertex *iVertex, const VertexCol *iVertices,//Vertex here is the PV
263		Bool_t printDebug) {
264
265		if (!fIsInitialized) {
#	Line 297 | Line 269 | Double_t JetIDMVA::MVAValue(const PFJet
269		if(!JetTools::passPFLooseId(iJet)) return -2.;
270
271		//set all input variables
272	+	fNVtx      = iVertices->GetEntries();
273		fJPt1      = iJet->Pt();
274		fJEta1     = iJet->RawMom().Eta();
275		fJPhi1     = iJet->RawMom().Phi();
276	<	fJM1       = iJet->Mass();
277	<
278	<	const mithep::PFCandidate *lLead     = JetTools::leadCand(iJet,-1);
279	<	const mithep::PFCandidate *lSecond   = JetTools::leadCand(iJet,-1,true);
280	<	const mithep::PFCandidate *lLeadNeut = JetTools::leadCand(iJet ,5);
281	<	const mithep::PFCandidate *lLeadEm   = JetTools::leadCand(iJet ,4);
282	<	const mithep::PFCandidate *lLeadCh   = JetTools::leadCand(iJet ,1);
283	<
284	<	fJD01         = JetTools::impactParameter(iJet,iVertex);
285	<	fJDZ1         = JetTools::impactParameter(iJet,iVertex,true);
286	<	fNPart1       = iJet->NPFCands();
287	<	fLPt1         = lLead    ->Pt();
288	<	fLEta1        = lLead    ->Eta();
289	<	fLPhi1        = lLead    ->Phi();
290	<	fSPt1         = lSecond  ->Pt();
291	<	fSEta1        = lSecond  ->Eta();
292	<	fSPhi1        = lSecond  ->Phi();
293	<	fNEPt1        = lLeadNeut->Pt();
321	<	fNEEta1       = lLeadNeut->Eta();
322	<	fNEPhi1       = lLeadNeut->Phi();
323	<	fEMPt1        = lLeadEm  ->Pt();
324	<	fEMEta1       = lLeadEm  ->Eta();
325	<	fEMPhi1       = lLeadEm  ->Phi();
326	<	fChPt1        = lLeadCh  ->Pt();
327	<	//fChEta1       = lLeadCh  ->Eta();
328	<	fChPhi1       = lLeadCh  ->Phi();
329	<	fLFr1         = lLead->Pt()/iJet->RawMom().Pt();
330	<
331	<	fDRLC1        = MathUtils::DeltaR(iJet->Mom(),lLead  ->Mom());
332	<	fDRLS1        = MathUtils::DeltaR(iJet->Mom(),lSecond->Mom());
333	<	fDRM1         = JetTools::dRMean (iJet,-1);
334	<	fDRNE1        = JetTools::dRMean (iJet, 5);
335	<	fDREM1        = JetTools::dRMean (iJet, 4);
336	<	fDRCH1        = JetTools::dRMean (iJet, 1);
337	<
338	<	double lMVA = fReader->EvaluateMVA( fMethodName );
339	<
276	>	fJD01      = JetTools::impactParameter(iJet,iVertex);
277	>	fJDZ1      = JetTools::impactParameter(iJet,iVertex,true);
278	>	fBeta      = JetTools::Beta(iJet,iVertex,fDZCut);
279	>	fBetaStar  = JetTools::betaStar(iJet,iVertex,iVertices,fDZCut);
280	>	fNCharged  = iJet->ChargedMultiplicity();
281	>	fNNeutrals = iJet->NeutralMultiplicity();
282	>
283	>	fDRMean    = JetTools::dRMean(iJet,-1);
284	>	fFrac01    = JetTools::frac  (iJet,0.1,0. ,-1);
285	>	fFrac02    = JetTools::frac  (iJet,0.2,0.1,-1);
286	>	fFrac03    = JetTools::frac  (iJet,0.3,0.2,-1);
287	>	fFrac04    = JetTools::frac  (iJet,0.4,0.3,-1);
288	>	fFrac05    = JetTools::frac  (iJet,0.5,0.4,-1);
289	>
290	>	double lMVA = 0;
291	>	if(fJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
292	>	if(fJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
293	>
294		if (printDebug == kTRUE) {
295		std::cout << "Debug Jet MVA: "
296	<	<< fNPV    << " "
297	<	<< fJPt1   << " "
298	<	<< fJEta1  << " "
299	<	<< fJPhi1  << " "
300	<	<< fJD01   << " "
301	<	<< fJDZ1   << " "
302	<	<< fJM1    << " "
303	<	<< fNPart1 << " "
304	<	<< fLPt1   << " "
305	<	<< fLEta1  << " "
306	<	<< fLPhi1  << " "
307	<	<< fSPt1   << " "
308	<	<< fSEta1  << " "
309	<	<< fSPhi1  << " "
310	<	<< fNEPt1  << " "
311	<	<< fNEEta1 << " "
358	<	<< fNEPhi1 << " "
359	<	<< fEMPt1  << " "
360	<	<< fEMEta1 << " "
361	<	<< fEMPhi1 << " "
362	<	<< fChPt1  << " "
363	<	<< fChPhi1 << " "
364	<	<< fLFr1   << " "
365	<	<< fDRLC1  << " "
366	<	<< fDRLS1  << " "
367	<	<< fDRM1   << " "
368	<	<< fDRNE1 << " "
369	<	<< fDREM1  << " "
370	<	<< fDRCH1  << " "
296	>	<< fNVtx      << " "
297	>	<< fJPt1      << " "
298	>	<< fJEta1     << " "
299	>	<< fJPhi1     << " "
300	>	<< fJD01      << " "
301	>	<< fJDZ1      << " "
302	>	<< fBeta      << " "
303	>	<< fBetaStar  << " "
304	>	<< fNCharged  << " "
305	>	<< fNNeutrals << " "
306	>	<< fDRMean    << " "
307	>	<< fFrac01    << " "
308	>	<< fFrac02    << " "
309	>	<< fFrac03    << " "
310	>	<< fFrac04    << " "
311	>	<< fFrac05
312		<< " === : === "
313		<< lMVA << " "
314		<< std::endl;
#	Line 375 | Line 316 | Double_t JetIDMVA::MVAValue(const PFJet
316
317		return lMVA;
318		}
319	<	Double_t JetIDMVA::correctedPt(const PFJet *iJet, FactorizedJetCorrector *iJetCorrector,const PileupEnergyDensityCol *iPUEnergyDensity) {
319	>	Double_t JetIDMVA::correctedPt(const PFJet *iJet, FactorizedJetCorrector *iJetCorrector,
320	>	const PileupEnergyDensityCol *iPUEnergyDensity,
321	>	RhoUtilities::RhoType type) {
322	>	Double_t Rho = 0.0;
323	>	switch(type) {
324	>	case RhoUtilities::MIT_RHO_VORONOI_HIGH_ETA:
325	>	Rho = iPUEnergyDensity->At(0)->Rho();
326	>	break;
327	>	case RhoUtilities::MIT_RHO_VORONOI_LOW_ETA:
328	>	Rho = iPUEnergyDensity->At(0)->RhoLowEta();
329	>	break;
330	>	case RhoUtilities::MIT_RHO_RANDOM_HIGH_ETA:
331	>	Rho = iPUEnergyDensity->At(0)->RhoRandom();
332	>	break;
333	>	case RhoUtilities::MIT_RHO_RANDOM_LOW_ETA:
334	>	Rho = iPUEnergyDensity->At(0)->RhoRandomLowEta();
335	>	break;
336	>	case RhoUtilities::CMS_RHO_RHOKT6PFJETS:
337	>	Rho = iPUEnergyDensity->At(0)->RhoKt6PFJets();
338	>	break;
339	>	default:
340	>	// use the old default
341	>	Rho = iPUEnergyDensity->At(0)->Rho();
342	>	break;
343	>	}
344	>
345		const FourVectorM rawMom = iJet->RawMom();
346		iJetCorrector->setJetEta(rawMom.Eta());
347		iJetCorrector->setJetPt (rawMom.Pt());
348		iJetCorrector->setJetPhi(rawMom.Phi());
349		iJetCorrector->setJetE  (rawMom.E());
350	<	iJetCorrector->setRho   (iPUEnergyDensity->At(0)->RhoHighEta());
350	>	iJetCorrector->setRho   (Rho);
351		iJetCorrector->setJetA  (iJet->JetArea());
352		iJetCorrector->setJetEMF(-99.0);
353		Double_t correction = iJetCorrector->getCorrection();

Diff Legend

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