ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/VHbb/plugins/Histos.h

Comparing UserCode/VHbb/plugins/Histos.h (file contents):
Revision 1.4 by bortigno, Fri May 11 09:46:37 2012 UTC vs.
Revision 1.7 by bortigno, Fri Jun 22 15:49:50 2012 UTC

#	Line 96 | Line 96 | public:
96		subDir = f.GetDirectory(suffix.c_str());
97		subDir->cd();
98
99	<	bin_btag_prob = 20;
100	<	min_btag_prob = 0;
101	<	max_btag_prob = 1;
102	<
103	<	bin_btag_count = 10;
104	<	min_btag_count = 0;
105	<	max_btag_count = 10;
99	>	bin_btag = 40;
100	>	min_btag = -1;
101	>	max_btag = 1;
102
103		//Candidates
104	<	BTagH_bJet1_csv = new TH1F(("BJet1_CSV"+suffix).c_str(),("BJet1 CSV ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob );
105	<	BTagH_bJet2_csv = new TH1F(("BJet2_CSV"+suffix).c_str(),("BJet2 CSV ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob );
106	<	/*     BTagH_bJet1_csvmva = new TH1F(("BJet1_CSVMVA"+suffix).c_str(),("BJet1 CSVMVA ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
107	<	/*     BTagH_bJet2_csvmva = new TH1F(("BJet2_CSVMVA"+suffix).c_str(),("BJet2 CSVMVA ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
112	<	/*     BTagH_bJet1_ssvhe = new TH1F(("BJet1_SSVHE"+suffix).c_str(),("BJet1 SSVHE ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
113	<	/*     BTagH_bJet2_ssvhe = new TH1F(("BJet2_SSVHE"+suffix).c_str(),("BJet2 SSVHE ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
114	<	/*     BTagH_bJet1_jpb = new TH1F(("BJet1_JPB"+suffix).c_str(),("BJet1 JPB ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
115	<	/*     BTagH_bJet2_jpb = new TH1F(("BJet2_JPB"+suffix).c_str(),("BJet2 JPB ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
116	<	/*     BTagH_bJet1_tche = new TH1F(("BJet1_TCHE"+suffix).c_str(),("BJet1 TCHE ("+suffix+")").c_str(), bin_btag_count, min_btag_count, max_btag_count ); */
117	<	/*     BTagH_bJet2_tche = new TH1F(("BJet2_TCHE"+suffix).c_str(),("BJet2 TCHE ("+suffix+")").c_str(), bin_btag_count, min_btag_count, max_btag_count ); */
118	<	/*     BTagH_bJet1_jp = new TH1F(("BJet1_JP"+suffix).c_str(),("BJet1 JP ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
119	<	/*     BTagH_bJet2_jp = new TH1F(("BJet2_JP"+suffix).c_str(),("BJet2 JP ("+suffix+")").c_str(), bin_btag_prob, min_btag_prob, max_btag_prob ); */
120	<	/*     BTagH_bJet1_tchp = new TH1F(("BJet1_TCHP"+suffix).c_str(),("BJet1 TCHP ("+suffix+")").c_str(), bin_btag_count, min_btag_count, max_btag_count ); */
121	<	/*     BTagH_bJet2_tchp = new TH1F(("BJet2_TCHP"+suffix).c_str(),("BJet2 TCHP ("+suffix+")").c_str(), bin_btag_count, min_btag_count, max_btag_count ); */
104	>	StH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix).c_str(),("Simple Jet1 bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
105	>	StH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix).c_str(),("Simple Jet2 bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
106	>	/*     StH_simpleJet1_bTag_reshaped = new TH1F(("SimpleJet1_bTag_reshaped"+suffix).c_str(),("Simple Jet1 bTag reshaped ("+suffix+")").c_str(), bin_btag, min_btag, max_btag ); */
107	>	/*     StH_simpleJet2_bTag_reshaped = new TH1F(("SimpleJet2_bTag_reshaped"+suffix).c_str(),("Simple Jet2 bTag reshaped ("+suffix+")").c_str(), bin_btag, min_btag, max_btag ); */
108
109		}
110
111		virtual void fill(ntupleReader &iEvent,float w) {
112
113	<	BTagH_bJet1_csv->Fill(iEvent.hJet_csv[0], w);
114	<	BTagH_bJet2_csv->Fill(iEvent.hJet_csv[1], w);
115	<	/*       BTagH_bJet1_csvmva->Fill(iEvent.H.jets.at(0).csvmva, w); */
116	<	/*       BTagH_bJet2_csvmva->Fill(iEvent.H.jets.at(1).csvmva, w); */
117	<	/*       BTagH_bJet1_ssvhe->Fill(iEvent.H.jets.at(0).ssvhe, w); */
118	<	/*       BTagH_bJet2_ssvhe->Fill(iEvent.H.jets.at(1).ssvhe, w); */
119	<	/*       BTagH_bJet1_tche->Fill(iEvent.H.jets.at(0).tche, w); */
120	<	/*       BTagH_bJet2_tche->Fill(iEvent.H.jets.at(1).tche, w); */
121	<	/*       BTagH_bJet1_tchp->Fill(iEvent.H.jets.at(0).tchp, w); */
136	<	/*       BTagH_bJet2_tchp->Fill(iEvent.H.jets.at(1).tchp, w); */
137	<	/*       BTagH_bJet1_jpb->Fill(iEvent.H.jets.at(0).jpb, w); */
138	<	/*       BTagH_bJet2_jpb->Fill(iEvent.H.jets.at(1).jpb, w); */
139	<	/*       BTagH_bJet1_jp->Fill(iEvent.H.jets.at(0).jp, w); */
140	<	/*       BTagH_bJet2_jp->Fill(iEvent.H.jets.at(1).jp, w); */
141	<
113	>	//vertex mas and btag sorted by btag
114	>	if( iEvent.hJet_CSV(0,0) > iEvent.hJet_CSV(1,0) ){
115	>	StH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(0,0), w);
116	>	StH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(1,0), w);
117	>	}
118	>	else{
119	>	StH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(1,0), w);
120	>	StH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(0,0), w);
121	>	}
122
123		}
124	<
125	<	TH1F * BTagH_bJet1_csv;
126	<	TH1F * BTagH_bJet2_csv;
127	<	/*   TH1F * BTagH_bJet1_csvmva; */
148	<	/*   TH1F * BTagH_bJet2_csvmva; */
149	<	/*   TH1F * BTagH_bJet1_ssvhe; */
150	<	/*   TH1F * BTagH_bJet2_ssvhe; */
151	<	/*   TH1F * BTagH_bJet1_jpb; */
152	<	/*   TH1F * BTagH_bJet2_jpb; */
153	<	/*   TH1F * BTagH_bJet1_tche; */
154	<	/*   TH1F * BTagH_bJet2_tche; */
155	<	/*   TH1F * BTagH_bJet1_jp; */
156	<	/*   TH1F * BTagH_bJet2_jp; */
157	<	/*   TH1F * BTagH_bJet1_tchp; */
158	<	/*   TH1F * BTagH_bJet2_tchp; */
159	<
124	>
125	>	TH1F * StH_simpleJet1_bTag;
126	>	TH1F * StH_simpleJet2_bTag;
127	>
128		private:
129
130	<	Int_t bin_btag_prob;
131	<	Double_t min_btag_prob;
132	<	Double_t max_btag_prob;
133	<
166	<	Int_t bin_btag_count;
167	<	Double_t min_btag_count;
168	<	Double_t max_btag_count;
169	<
130	>	Int_t bin_btag;
131	>	Double_t min_btag;
132	>	Double_t max_btag;
133	>
134		};
135
136
#	Line 275 | Line 239 | public:
239		StH_simpleJet2_eta = new TH1F(("SimpleJet2_eta"+suffix).c_str(),("Simple Jet2 eta ("+suffix+")").c_str(), bin_eta, min_eta, max_eta );
240		StH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix).c_str(),("Simple Jet1 bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
241		StH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix).c_str(),("Simple Jet2 bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
242	+	StH_simpleJet1_bTag_reshaped = new TH1F(("SimpleJet1_bTag_reshaped"+suffix).c_str(),("Simple Jet1 bTag reshaped ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
243	+	StH_simpleJet2_bTag_reshaped = new TH1F(("SimpleJet2_bTag_reshaped"+suffix).c_str(),("Simple Jet2 bTag reshaped ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
244		StH_simpleJet1_vtxMass = new TH1F(("SimpleJet1_vtxMass"+suffix).c_str(),("Simple Jet1 vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass );
245		StH_simpleJet2_vtxMass = new TH1F(("SimpleJet2_vtxMass"+suffix).c_str(),("Simple Jet2 vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass );
246
247		StH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix).c_str(),("Simple Jets bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
248		StH_simpleJets_bTagSum = new TH1F(("SimpleJets_bTagSum"+suffix).c_str(),("Simple Jets bTagSum ("+suffix+")").c_str(), bin_btag*2, min_btag, max_btag*2 );
249	+	StH_simpleJets_bTag_reshaped = new TH1F(("SimpleJets_bTag_reshaped"+suffix).c_str(),("Simple Jets bTag reshaped ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
250	+	StH_simpleJets_bTag_reshapedSum = new TH1F(("SimpleJets_bTag_reshapedSum"+suffix).c_str(),("Simple Jets bTag reshapedSum ("+suffix+")").c_str(), bin_btag*2, min_btag, max_btag*2 );
251		StH_simpleJets_vtxMass = new TH1F(("SimpleJets_vtxMass"+suffix).c_str(),("Simple Jets vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass );
252
253		StH_addJet1_pt = new TH1F(("AddJet1_pt"+suffix).c_str(),("Additional Jet1 pt ("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
#	Line 295 | Line 263 | public:
263		StH_HdPhi = new TH1F(("HdPhi"+suffix).c_str(),("Higgs deltaPhi ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
264		StH_HdEta = new TH1F(("HdEta"+suffix).c_str(),("Higgs deltaEta ("+suffix+")").c_str(), bin_deltaEta, min_deltaEta, max_deltaEta );
265		StH_HHel = new TH1F(("HiggsHel"+suffix).c_str(),("Higgs helicity angle ("+suffix+")").c_str(), bin_hel, min_hel, max_hel );
266	<	StH_HPullAngle = new TH1F(("HiggsPullAngle"+suffix).c_str(),("Higgs pull angle ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
266	>	StH_HPullAngle1 = new TH1F(("HiggsPullAngle1"+suffix).c_str(),("Higgs pull angle1 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
267	>	StH_HPullAngle2 = new TH1F(("HiggsPullAngle2"+suffix).c_str(),("Higgs pull angle2 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
268	>	StH_HPullAK7Angle1 = new TH1F(("HiggsPullAK7Angle1"+suffix).c_str(),("Higgs pullAK7 angle1 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
269	>	StH_HPullAK7Angle2 = new TH1F(("HiggsPullAK7Angle2"+suffix).c_str(),("Higgs pullAK7 angle2 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
270		StH_HdR_addJets = new TH1F(("HdR_addJets"+suffix).c_str(),("Higgs deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
271		StH_HdR_leading_addJets = new TH1F(("HdR_leading_addJets"+suffix).c_str(),("Higgs deltaR leading-add Jets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
272		StH_HdR_following_addJets = new TH1F(("HdR_following_addJets"+suffix).c_str(),("Higgs deltaR following-add Jets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
#	Line 308 | Line 279 | public:
279		StH_massDrop = new TH1F(("MassDrop"+suffix).c_str(),(" Higgs Mass Drop ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
280		StH_y12 = new TH1F(("Y12"+suffix).c_str(),(" Higgs Y12 ("+suffix+")").c_str(), bin_btag, min_btag, max_btag*1.2 );
281
311	–
282		StH_ZMass = new TH1F(("ZMass"+suffix).c_str(),(" Z Mass ("+suffix+")").c_str(), 100, 50, 150 );
283		StH_ZPt = new TH1F(("ZPt"+suffix).c_str(),(" Z Pt ("+suffix+")").c_str(), bin_pt, min_pt, max_pt );
284		StH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix).c_str(),(" ZH delta Phi ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
#	Line 347 | Line 317 | public:
317		StH_simpleJet2_phi->Fill(iEvent.hJet_phi[1], w);
318		StH_simpleJet1_eta->Fill(iEvent.hJet_eta[0], w);
319		StH_simpleJet2_eta->Fill(iEvent.hJet_eta[1], w);
320	+
321		//vertex mas and btag sorted by btag
322	<	if( iEvent.hJet_csv[0] > iEvent.hJet_csv[1] ){
323	<	StH_simpleJet1_bTag->Fill(iEvent.hJet_csv[0], w);
324	<	StH_simpleJet2_bTag->Fill(iEvent.hJet_csv[1], w);
322	>	if( iEvent.hJet_CSV(0,0) > iEvent.hJet_CSV(1,0) ){
323	>	StH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(0,0), w);
324	>	StH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(1,0), w);
325		StH_simpleJet1_vtxMass->Fill(iEvent.hJet_vtxMass[0], w);
326		}
327		else{
328	<	StH_simpleJet1_bTag->Fill(iEvent.hJet_csv[1], w);
329	<	StH_simpleJet2_bTag->Fill(iEvent.hJet_csv[0], w);
328	>	StH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(1,0), w);
329	>	StH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(0,0), w);
330		StH_simpleJet2_vtxMass->Fill(iEvent.hJet_vtxMass[1], w);
331	<	}
331	>	}
332	>
333		//here I fill both jets in one histo
334	<	StH_simpleJets_bTagSum->Fill(iEvent.hJet_csv[0]+iEvent.hJet_csv[1], w);
335	<	StH_simpleJets_bTag->Fill(iEvent.hJet_csv[0], w);
336	<	StH_simpleJets_bTag->Fill(iEvent.hJet_csv[1], w);
334	>	StH_simpleJets_bTagSum->Fill(iEvent.hJet_CSV(0,0)+iEvent.hJet_CSV(1,0), w);
335	>	StH_simpleJets_bTag->Fill(iEvent.hJet_CSV(0,0), w);
336	>	StH_simpleJets_bTag->Fill(iEvent.hJet_CSV(1,0), w);
337		StH_simpleJets_vtxMass->Fill(iEvent.hJet_vtxMass[0], w);
338		StH_simpleJets_vtxMass->Fill(iEvent.hJet_vtxMass[1], w);
339	+
340	+	//reshaped btag
341	+	//vertex mas and btag sorted by btag
342	+	if( iEvent.hJet_csv_nominal[0] > iEvent.hJet_csv_nominal[1] ){
343	+	StH_simpleJet1_bTag_reshaped->Fill(iEvent.hJet_csv_nominal[0], w);
344	+	StH_simpleJet2_bTag_reshaped->Fill(iEvent.hJet_csv_nominal[1], w);
345	+	}
346	+	else{
347	+	StH_simpleJet1_bTag_reshaped->Fill(iEvent.hJet_csv_nominal[1], w);
348	+	StH_simpleJet2_bTag_reshaped->Fill(iEvent.hJet_csv_nominal[0], w);
349	+	}
350	+	StH_simpleJets_bTag_reshapedSum->Fill(iEvent.hJet_csv_nominal[0]+iEvent.hJet_csv_nominal[1], w);
351	+	StH_simpleJets_bTag_reshaped->Fill(iEvent.hJet_csv_nominal[0], w);
352	+	StH_simpleJets_bTag_reshaped->Fill(iEvent.hJet_csv_nominal[1], w);
353
354		StH_HMass->Fill(iEvent.H_mass, w);
355		StH_HPt->Fill(iEvent.H_pt, w);
#	Line 372 | Line 358 | public:
358		StH_HdEta->Fill(iEvent.H_dEta, w);
359
360		StH_HHel->Fill(iEvent.hJet_cosTheta[0], w);
361	<	StH_HPullAngle->Fill(iEvent.deltaPullAngle, w);
361	>
362	>	StH_HPullAngle1->Fill(iEvent.deltaPullAngle, w);
363	>	StH_HPullAngle2->Fill(iEvent.deltaPullAngle2, w);
364	>	StH_HPullAK7Angle1->Fill(iEvent.deltaPullAngleAK7, w);
365	>	StH_HPullAK7Angle2->Fill(iEvent.deltaPullAngle2AK7, w);
366
367		bool debug_ = false;
368		TLorentzVector higgs, hJet1, hJet2, addJet;
#	Line 475 | Line 465 | public:
465		TH1F * StH_simpleJet2_eta;
466		TH1F * StH_simpleJet1_bTag;
467		TH1F * StH_simpleJet2_bTag;
468	+	TH1F * StH_simpleJet1_bTag_reshaped;
469	+	TH1F * StH_simpleJet2_bTag_reshaped;
470		TH1F * StH_simpleJet1_vtxMass;
471		TH1F * StH_simpleJet2_vtxMass;
472
473		TH1F * StH_simpleJets_bTag;
474		TH1F * StH_simpleJets_bTagSum;
475	+	TH1F * StH_simpleJets_bTag_reshaped;
476	+	TH1F * StH_simpleJets_bTag_reshapedSum;
477		TH1F * StH_simpleJets_vtxMass;
478
479		TH1F * StH_addJet1_pt;
#	Line 495 | Line 489 | public:
489		TH1F * StH_HdPhi;
490		TH1F * StH_HdEta;
491		TH1F * StH_HHel;
492	<	TH1F * StH_HPullAngle;
492	>	TH1F * StH_HPullAngle1;
493	>	TH1F * StH_HPullAngle2;
494	>	TH1F * StH_HPullAK7Angle1;
495	>	TH1F * StH_HPullAK7Angle2;
496		TH1F * StH_HdR_addJets;
497		TH1F * StH_HdR_leading_addJets;
498		TH1F * StH_HdR_following_addJets;
#	Line 694 | Line 691 | public:
691		SystJERUPH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,+1).Pt(), w);
692
693		//btag sorted by btag
694	<	if( iEvent.hJet_csv_cUP(0) > iEvent.hJet_csv_cUP(1) ){
695	<	SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cUP(0), w);
696	<	SystUPH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cUP(1), w);
694	>	if( iEvent.hJet_CSV(0,1) > iEvent.hJet_CSV(1,1) ){
695	>	SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(0,1), w);
696	>	SystUPH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(1,1), w);
697		}
698		else{
699	<	SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cUP(1), w);
700	<	SystUPH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cUP(0), w);
699	>	SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(1,1), w);
700	>	SystUPH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(0,1), w);
701		}
702	<	if( iEvent.hJet_csv_cFUP(0) > iEvent.hJet_csv_cFUP(1) ){
703	<	SystFUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cFUP(0), w);
704	<	SystFUPH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cFUP(1), w);
702	>	if( iEvent.hJet_CSV(0,2) > iEvent.hJet_CSV(1,2) ){
703	>	SystFUPH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(0,2), w);
704	>	SystFUPH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(1,2), w);
705		}
706		else{
707	<	SystFUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cFUP(1), w);
708	<	SystFUPH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cFUP(0), w);
707	>	SystFUPH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(1,2), w);
708	>	SystFUPH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(0,2), w);
709		}
710
711		//here I fill both jets in one histo
712	<	SystUPH_simpleJets_bTag->Fill(iEvent.hJet_csv_cUP(0), w);
713	<	SystUPH_simpleJets_bTag->Fill(iEvent.hJet_csv_cUP(1), w);
714	<	SystFUPH_simpleJets_bTag->Fill(iEvent.hJet_csv_cFUP(0), w);
715	<	SystFUPH_simpleJets_bTag->Fill(iEvent.hJet_csv_cFUP(1), w);
712	>	SystUPH_simpleJets_bTag->Fill(iEvent.hJet_CSV(0,1), w);
713	>	SystUPH_simpleJets_bTag->Fill(iEvent.hJet_CSV(1,1), w);
714	>	SystFUPH_simpleJets_bTag->Fill(iEvent.hJet_CSV(0,2), w);
715	>	SystFUPH_simpleJets_bTag->Fill(iEvent.hJet_CSV(1,2), w);
716
717		SystUPH_HMass->Fill(iEvent.H_jecUP().M(), w);
718		SystUPH_HPt->Fill(iEvent.H_jecUP().Pt(), w);
#	Line 737 | Line 734 | public:
734		SystJERDOWNH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,-1).Pt(), w);
735
736		//btag sorted by btag
737	<	if( iEvent.hJet_csv_cDOWN(0) > iEvent.hJet_csv_cDOWN(1) ){
738	<	SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cDOWN(0), w);
739	<	SystDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cDOWN(1), w);
737	>	if( iEvent.hJet_CSV(0,-1) > iEvent.hJet_CSV(1,-1) ){
738	>	SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(0,-1), w);
739	>	SystDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(1,-1), w);
740		}
741		else{
742	<	SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cDOWN(1), w);
743	<	SystDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cDOWN(0), w);
742	>	SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(1,-1), w);
743	>	SystDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(0,-1), w);
744		}
745	<	if( iEvent.hJet_csv_cFDOWN(0) > iEvent.hJet_csv_cFDOWN(1) ){
746	<	SystFDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cFDOWN(0), w);
747	<	SystFDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cFDOWN(1), w);
745	>	if( iEvent.hJet_CSV(0,-2) > iEvent.hJet_CSV(1,-2) ){
746	>	SystFDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(0,-2), w);
747	>	SystFDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(1,-2), w);
748		}
749		else{
750	<	SystFDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cFDOWN(1), w);
751	<	SystFDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_csv_cFDOWN(0), w);
750	>	SystFDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_CSV(1,-2), w);
751	>	SystFDOWNH_simpleJet2_bTag->Fill(iEvent.hJet_CSV(0,-2), w);
752		}
753
754		//here I fill both jets in one histo
755	<	SystDOWNH_simpleJets_bTag->Fill(iEvent.hJet_csv_cDOWN(0), w);
756	<	SystDOWNH_simpleJets_bTag->Fill(iEvent.hJet_csv_cDOWN(1), w);
757	<	SystFDOWNH_simpleJets_bTag->Fill(iEvent.hJet_csv_cFDOWN(0), w);
758	<	SystFDOWNH_simpleJets_bTag->Fill(iEvent.hJet_csv_cFDOWN(1), w);
755	>	SystDOWNH_simpleJets_bTag->Fill(iEvent.hJet_CSV(0,-1), w);
756	>	SystDOWNH_simpleJets_bTag->Fill(iEvent.hJet_CSV(1,-1), w);
757	>	SystFDOWNH_simpleJets_bTag->Fill(iEvent.hJet_CSV(0,-2), w);
758	>	SystFDOWNH_simpleJets_bTag->Fill(iEvent.hJet_CSV(1,-2), w);
759		SystDOWNH_HMass->Fill(iEvent.H_jecDOWN().M(), w);
760		SystDOWNH_HPt->Fill(iEvent.H_jecDOWN().Pt(), w);
761		SystDOWNH_ZH_dPhi->Fill( TMath::Abs( iEvent.H_jecDOWN().DeltaPhi( iEvent.VectorBoson() ) ), w);

Diff Legend

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