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.1 by bortigno, Fri Apr 27 14:25:30 2012 UTC vs.
Revision 1.6 by bortigno, Mon Jun 11 14:51:58 2012 UTC

#	Line 4 | Line 4
4		#include <TH1F.h>
5		/* #include "DataFormats/GeometryVector/interface/VectorUtil.h" */
6		/* #include "DataFormats/Math/interface/deltaR.h" */
7	<	#include "CutsAndHistos.h"
7	>	#include "../interface/CutsAndHistos.h"
8		#include "TLorentzVector.h"
9		#include <sstream>
10		#include "TKey.h"
#	Line 250 | Line 250 | public:
250		StH_nPVs = new TH1F(("nPVs"+suffix).c_str(),("Number of reconstructed PV ("+suffix+")").c_str(), bin_npvs, min_npvs, max_npvs );
251		StH_addJets = new TH1F(("addJets"+suffix).c_str(),("Number of additional AK5 ("+suffix+")").c_str(), bin_addJets, min_addJets, max_addJets );
252
253	+	StH_typeImet_et = new TH1F(("typeIMET_et"+suffix).c_str(),("Type corrected missing transverse energy ("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.5 );
254		StH_met_et = new TH1F(("MET_et"+suffix).c_str(),("Missing transverse energy ("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.5 );
255		StH_met_sig = new TH1F(("MET_sig"+suffix).c_str(),("MET significance("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.1 );
256
#	Line 278 | Line 279 | public:
279		StH_simpleJet2_vtxMass = new TH1F(("SimpleJet2_vtxMass"+suffix).c_str(),("Simple Jet2 vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass );
280
281		StH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix).c_str(),("Simple Jets bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
282	+	StH_simpleJets_bTagSum = new TH1F(("SimpleJets_bTagSum"+suffix).c_str(),("Simple Jets bTagSum ("+suffix+")").c_str(), bin_btag*2, min_btag, max_btag*2 );
283		StH_simpleJets_vtxMass = new TH1F(("SimpleJets_vtxMass"+suffix).c_str(),("Simple Jets vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass );
284
285		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 293 | Line 295 | public:
295		StH_HdPhi = new TH1F(("HdPhi"+suffix).c_str(),("Higgs deltaPhi ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
296		StH_HdEta = new TH1F(("HdEta"+suffix).c_str(),("Higgs deltaEta ("+suffix+")").c_str(), bin_deltaEta, min_deltaEta, max_deltaEta );
297		StH_HHel = new TH1F(("HiggsHel"+suffix).c_str(),("Higgs helicity angle ("+suffix+")").c_str(), bin_hel, min_hel, max_hel );
298	<	StH_HPullAngle = new TH1F(("HiggsPullAngle"+suffix).c_str(),("Higgs pull angle ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
298	>	StH_HPullAngle1 = new TH1F(("HiggsPullAngle1"+suffix).c_str(),("Higgs pull angle1 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
299	>	StH_HPullAngle2 = new TH1F(("HiggsPullAngle2"+suffix).c_str(),("Higgs pull angle2 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
300	>	StH_HPullAK7Angle1 = new TH1F(("HiggsPullAK7Angle1"+suffix).c_str(),("Higgs pullAK7 angle1 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
301	>	StH_HPullAK7Angle2 = new TH1F(("HiggsPullAK7Angle2"+suffix).c_str(),("Higgs pullAK7 angle2 ("+suffix+")").c_str(), bin_phi, min_phi, max_phi );
302		StH_HdR_addJets = new TH1F(("HdR_addJets"+suffix).c_str(),("Higgs deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
303	+	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 );
304	+	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 );
305	+	StH_HdRnorm_addJets = new TH1F(("HdRnorm_addJets"+suffix).c_str(),("Higgs normalized deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
306		StH_HMass_addJets = new TH1F(("HiggsMass_addJets"+suffix).c_str(),(" Higgs Mass with addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
307	+	StH_HMass_leading_addJets = new TH1F(("HiggsMass_leading_addJets"+suffix).c_str(),(" Higgs Mass with leading addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
308	+	StH_HMass_following_addJets = new TH1F(("HiggsMass_following_addJets"+suffix).c_str(),(" Higgs Mass with following addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
309	+	StH_HMass_both_addJets = new TH1F(("HiggsMass_both_addJets"+suffix).c_str(),(" Higgs Mass with both addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
310
311		StH_massDrop = new TH1F(("MassDrop"+suffix).c_str(),(" Higgs Mass Drop ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
312		StH_y12 = new TH1F(("Y12"+suffix).c_str(),(" Higgs Y12 ("+suffix+")").c_str(), bin_btag, min_btag, max_btag*1.2 );
#	Line 315 | Line 326 | public:
326		StH_nPVs->Fill(iEvent.nPVs, w);
327		StH_addJets->Fill(iEvent.CountAddJets(), w);
328
329	<	StH_met_et->Fill( TMath::Min(iEvent.METnoPUCh_et, iEvent.METnoPU_et ), w);
329	>	StH_typeImet_et->Fill(iEvent.typeIcorrMET(), w);
330	>	//    StH_met_et->Fill( TMath::Min(iEvent.METnoPUCh_et, iEvent.METnoPU_et ), w);
331	>	StH_met_et->Fill( iEvent.MET_et, w);
332		StH_met_sig->Fill(iEvent.MET_sig, w);
333
334		StH_vectorLepton1_pt->Fill(iEvent.vLepton_pt[0], w);
#	Line 349 | Line 362 | public:
362		StH_simpleJet2_vtxMass->Fill(iEvent.hJet_vtxMass[1], w);
363		}
364		//here I fill both jets in one histo
365	+	StH_simpleJets_bTagSum->Fill(iEvent.hJet_csv[0]+iEvent.hJet_csv[1], w);
366		StH_simpleJets_bTag->Fill(iEvent.hJet_csv[0], w);
367		StH_simpleJets_bTag->Fill(iEvent.hJet_csv[1], w);
368		StH_simpleJets_vtxMass->Fill(iEvent.hJet_vtxMass[0], w);
#	Line 361 | Line 375 | public:
375		StH_HdEta->Fill(iEvent.H_dEta, w);
376
377		StH_HHel->Fill(iEvent.hJet_cosTheta[0], w);
378	<	StH_HPullAngle->Fill(iEvent.deltaPullAngle, w);
378	>
379	>	StH_HPullAngle1->Fill(iEvent.deltaPullAngle, w);
380	>	StH_HPullAngle2->Fill(iEvent.deltaPullAngle2, w);
381	>	StH_HPullAK7Angle1->Fill(iEvent.deltaPullAngleAK7, w);
382	>	StH_HPullAK7Angle2->Fill(iEvent.deltaPullAngle2AK7, w);
383
384		bool debug_ = false;
385		TLorentzVector higgs, hJet1, hJet2, addJet;
#	Line 394 | Line 412 | public:
412		}
413		if(iEvent.CountAddJets() > 0 ){
414		StH_HdR_addJets->Fill(higgs.DeltaR(addJets.at(dRmin_idx)), w);
415	<	if( iEvent.aJet_csv[dRmin_idx] > 0.7 )
415	>	StH_HdR_leading_addJets->Fill(hJet1.DeltaR(addJets.at(dRmin_idx)), w);
416	>	StH_HdR_following_addJets->Fill(hJet2.DeltaR(addJets.at(dRmin_idx)), w);
417	>	StH_HdRnorm_addJets->Fill(higgs.DeltaR(addJets.at(dRmin_idx)) / iEvent.H_dR, w);
418	>	if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet1.DeltaR(addJets.at(dRmin_idx)) < 2 )
419	>	StH_HMass_leading_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
420	>	else StH_HMass_leading_addJets->Fill(higgs.M(), w);
421	>	if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet2.DeltaR(addJets.at(dRmin_idx)) < 2 )
422	>	StH_HMass_following_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
423	>	else StH_HMass_following_addJets->Fill(higgs.M(), w);
424	>	if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet2.DeltaR(addJets.at(dRmin_idx)) < 2  ||  hJet1.DeltaR(addJets.at(dRmin_idx)) < 2 )
425	>	StH_HMass_both_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
426	>	else StH_HMass_both_addJets->Fill(higgs.M(), w);
427	>	if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  higgs.DeltaR(addJets.at(dRmin_idx)) < 2 )
428		StH_HMass_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
429		else StH_HMass_addJets->Fill(higgs.M(), w);
430		StH_addJet1_pt->Fill(iEvent.aJet_pt[dRmin_idx], w);
#	Line 427 | Line 457 | public:
457
458		TH1F * StH_nPVs;
459		TH1F * StH_addJets;
460	+	TH1F * StH_typeImet_et;
461		TH1F * StH_met_et;
462		TH1F * StH_met_sig;
463
#	Line 455 | Line 486 | public:
486		TH1F * StH_simpleJet2_vtxMass;
487
488		TH1F * StH_simpleJets_bTag;
489	+	TH1F * StH_simpleJets_bTagSum;
490		TH1F * StH_simpleJets_vtxMass;
491
492		TH1F * StH_addJet1_pt;
#	Line 470 | Line 502 | public:
502		TH1F * StH_HdPhi;
503		TH1F * StH_HdEta;
504		TH1F * StH_HHel;
505	<	TH1F * StH_HPullAngle;
505	>	TH1F * StH_HPullAngle1;
506	>	TH1F * StH_HPullAngle2;
507	>	TH1F * StH_HPullAK7Angle1;
508	>	TH1F * StH_HPullAK7Angle2;
509		TH1F * StH_HdR_addJets;
510	+	TH1F * StH_HdR_leading_addJets;
511	+	TH1F * StH_HdR_following_addJets;
512	+
513	+	TH1F * StH_HdRnorm_addJets;
514		TH1F * StH_HMass_addJets;
515	+	TH1F * StH_HMass_leading_addJets;
516	+	TH1F * StH_HMass_following_addJets;
517	+	TH1F * StH_HMass_both_addJets;
518
519		TH1F * StH_massDrop;
520		TH1F * StH_y12;
#	Line 614 | Line 656 | public:
656		SystUPH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystUP").c_str(),(" ZH delta Phi ("+suffix+"SystUP"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
657		SystUPH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystUP").c_str(),("Higgs deltaR addJets ("+suffix+"SystUP"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
658
659	+	SystJERUPH_simpleJet1_pt = new TH1F(("SimpleJet1_pt"+suffix+"SystJERUP").c_str(),("Simple Jet1 pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
660	+	SystJERUPH_simpleJet2_pt = new TH1F(("SimpleJet2_pt"+suffix+"SystJERUP").c_str(),("Simple Jet2 pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
661	+	SystJERUPH_HMass = new TH1F(("HiggsMass"+suffix+"SystJERUP").c_str(),(" Higgs Mass ("+suffix+"SystJERUP"+")").c_str(), bin_mass, min_mass, max_mass );
662	+	SystJERUPH_HPt = new TH1F(("HiggsPt"+suffix+"SystJERUP").c_str(),(" Higgs Pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt );
663	+	SystJERUPH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystJERUP").c_str(),(" ZH delta Phi ("+suffix+"SystJERUP"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
664	+
665		SystUPH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix+"SystUP").c_str(),("Simple Jet1 bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
666		SystUPH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystUP").c_str(),("Simple Jet2 bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
667		SystUPH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix+"SystUP").c_str(),("Simple Jets bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
#	Line 627 | Line 675 | public:
675		SystDOWNH_HMass = new TH1F(("HiggsMass"+suffix+"SystDOWN").c_str(),(" Higgs Mass ("+suffix+"SystDOWN"+")").c_str(), bin_mass, min_mass, max_mass );
676		SystDOWNH_HPt = new TH1F(("HiggsPt"+suffix+"SystDOWN").c_str(),(" Higgs Pt ("+suffix+"SystDOWN"+")").c_str(), bin_pt, min_pt, max_pt );
677		SystDOWNH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystDOWN").c_str(),(" ZH delta Phi ("+suffix+"SystDOWN"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
678	<	SystDOWNH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystDOWN").c_str(),("Higgs deltaR addJets ("+suffix+"SystDOWN"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
678	>	SystDOWNH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystDWON").c_str(),("Higgs deltaR addJets ("+suffix+"SystDOWN"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
679	>
680	>	SystJERDOWNH_simpleJet1_pt = new TH1F(("SimpleJet1_pt"+suffix+"SystJERDOWN").c_str(),("Simple Jet1 pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
681	>	SystJERDOWNH_simpleJet2_pt = new TH1F(("SimpleJet2_pt"+suffix+"SystJERDOWN").c_str(),("Simple Jet2 pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
682	>	SystJERDOWNH_HMass = new TH1F(("HiggsMass"+suffix+"SystJERDOWN").c_str(),(" Higgs Mass ("+suffix+"SystJERDOWN"+")").c_str(), bin_mass, min_mass, max_mass );
683	>	SystJERDOWNH_HPt = new TH1F(("HiggsPt"+suffix+"SystJERDOWN").c_str(),(" Higgs Pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt );
684	>	SystJERDOWNH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystJERDOWN").c_str(),(" ZH delta Phi ("+suffix+"SystJERDOWN"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
685
686		SystDOWNH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix+"SystDOWN").c_str(),("Simple Jet1 bTag ("+suffix+"SystDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
687		SystDOWNH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystDOWN").c_str(),("Simple Jet2 bTag ("+suffix+"SystDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
#	Line 636 | Line 690 | public:
690		SystFDOWNH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystFDOWN").c_str(),("Simple Jet2 bTag ("+suffix+"SystFDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
691		SystFDOWNH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix+"SystFDOWN").c_str(),("Simple Jets bTag ("+suffix+"SystFDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
692
639	–
693		}
694
695		virtual void fill(ntupleReader &iEvent,float w) {
696
697		TLorentzVector higgs, addJet;
698
646	–
699		//UPSCALE
700		SystUPH_addJets->Fill(iEvent.CountAddJets_jec(+1), w);
701		SystUPH_simpleJet1_pt->Fill(iEvent.hJet_pt_jecUP(0), w);
702		SystUPH_simpleJet2_pt->Fill(iEvent.hJet_pt_jecUP(1), w);
703	+	SystJERUPH_simpleJet1_pt->Fill(iEvent.hJet_jer(0,+1).Pt(), w);
704	+	SystJERUPH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,+1).Pt(), w);
705	+
706		//btag sorted by btag
707		if( iEvent.hJet_csv_cUP(0) > iEvent.hJet_csv_cUP(1) ){
708		SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cUP(0), w);
#	Line 675 | Line 730 | public:
730		SystUPH_HMass->Fill(iEvent.H_jecUP().M(), w);
731		SystUPH_HPt->Fill(iEvent.H_jecUP().Pt(), w);
732		SystUPH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jecUP().DeltaPhi( iEvent.VectorBoson() ) ), w);
733	+
734	+	SystJERUPH_HMass->Fill(iEvent.H_jer(+1).M(), w);
735	+	SystJERUPH_HPt->Fill(iEvent.H_jer(+1).Pt(), w);
736	+	SystJERUPH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jer(+1).DeltaPhi( iEvent.VectorBoson() ) ), w);
737	+
738		higgs=iEvent.H_jecUP();
739		addJet.SetPtEtaPhiE(iEvent.aJet_pt_jecUP(0),iEvent.aJet_eta[0],iEvent.aJet_phi[0],iEvent.aJet_e[0]);
740		SystUPH_HdR_addJets->Fill(higgs.DeltaR(addJet), w);
#	Line 683 | Line 743 | public:
743		SystDOWNH_addJets->Fill(iEvent.CountAddJets_jec(-1), w);
744		SystDOWNH_simpleJet1_pt->Fill(iEvent.hJet_pt_jecDOWN(0), w);
745		SystDOWNH_simpleJet2_pt->Fill(iEvent.hJet_pt_jecDOWN(1), w);
746	+	SystJERDOWNH_simpleJet1_pt->Fill(iEvent.hJet_jer(0,-1).Pt(), w);
747	+	SystJERDOWNH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,-1).Pt(), w);
748	+
749		//btag sorted by btag
750		if( iEvent.hJet_csv_cDOWN(0) > iEvent.hJet_csv_cDOWN(1) ){
751		SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cDOWN(0), w);
#	Line 709 | Line 772 | public:
772		SystDOWNH_HMass->Fill(iEvent.H_jecDOWN().M(), w);
773		SystDOWNH_HPt->Fill(iEvent.H_jecDOWN().Pt(), w);
774		SystDOWNH_ZH_dPhi->Fill( TMath::Abs( iEvent.H_jecDOWN().DeltaPhi( iEvent.VectorBoson() ) ), w);
775	+
776	+	SystJERDOWNH_HMass->Fill(iEvent.H_jer(-1).M(), w);
777	+	SystJERDOWNH_HPt->Fill(iEvent.H_jer(-1).Pt(), w);
778	+	SystJERDOWNH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jer(-1).DeltaPhi( iEvent.VectorBoson() ) ), w);
779	+
780		higgs=iEvent.H_jecDOWN();
781		addJet.SetPtEtaPhiE(iEvent.aJet_pt_jecDOWN(0),iEvent.aJet_eta[0],iEvent.aJet_phi[0],iEvent.aJet_e[0]);
782		SystDOWNH_HdR_addJets->Fill(higgs.DeltaR(addJet), w);
783	<
783	>
784		}
785
786		TH1F * SystUPH_addJets;
#	Line 728 | Line 796 | public:
796		TH1F * SystFUPH_simpleJet1_bTag;
797		TH1F * SystFUPH_simpleJet2_bTag;
798		TH1F * SystFUPH_simpleJets_bTag;
799	+	TH1F * SystJERUPH_simpleJet1_pt;
800	+	TH1F * SystJERUPH_simpleJet2_pt;
801	+	TH1F * SystJERUPH_HMass;
802	+	TH1F * SystJERUPH_HPt;
803	+	TH1F * SystJERUPH_ZH_dPhi;
804
805		TH1F * SystDOWNH_addJets;
806		TH1F * SystDOWNH_simpleJet1_pt;
#	Line 742 | Line 815 | public:
815		TH1F * SystFDOWNH_simpleJet1_bTag;
816		TH1F * SystFDOWNH_simpleJet2_bTag;
817		TH1F * SystFDOWNH_simpleJets_bTag;
818	+	TH1F * SystJERDOWNH_simpleJet1_pt;
819	+	TH1F * SystJERDOWNH_simpleJet2_pt;
820	+	TH1F * SystJERDOWNH_HMass;
821	+	TH1F * SystJERDOWNH_HPt;
822	+	TH1F * SystJERDOWNH_ZH_dPhi;
823
824
825		private:

Diff Legend

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