ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/TreeReader.C

Comparing UserCode/LJMet/MultivariateAnalysis/root/TreeReader.C (file contents):
Revision 1.1 by jindal, Wed Jan 28 04:14:16 2009 UTC vs.
Revision 1.2 by jindal, Mon Feb 9 22:09:28 2009 UTC

#	Line 14 | Line 14
14		int TreeReader(char *infile = "temp.root")
15		{
16
17	<	TFile *outf = new TFile("Ntuple_goodglobalmuon.root","RECREATE");
17	>	TFile *outf = new TFile("eidRobustTight.root","RECREATE");
18	>
19	>
20		// TH1D * hEt3;
21	<	TH1D * hlepton_calo_iso = new TH1D( "hlepton_calo_iso", "lepton calo iso", 120, 0., 120.);
22	<	TH1D * hlepton_track_iso = new TH1D( "hlepton_track_iso", "lepton track iso", 120, 0., 120.);
21	>	TH1D * hlepton_calo_iso = new TH1D("hlepton_calo_iso", "lepton calo iso", 120, 0., 120.);
22	>	TH1D * hlepton_track_iso = new TH1D("hlepton_track_iso", "lepton track iso", 120, 0., 120.);
23		// TH1 * hmuon_isGood_AllGlobalMuons = new TH1( "hmuon_isGood_AllGlobalMuons", "muon isGood AllGlobalMuons", 120, 0., 120.);
24		//TH1 * hmuon_isGood_AllStandAloneMuons = new TH1( "hmuon_isGood_AllStandAloneMuons", "muon isGood AllStandAloneMuons", 120, 0., 120.);
25
26
27	<	TH1D * hMt = new TH1D( "hMt", "Mt 3 (all) ", 100, 0., 90.);
28	<	TH1D * hMt_rel_iso = new TH1D( "hMt_rel_iso", "Mt 3 (relative iso)", 100, 0., 90.);
29	<	TH1D * hMt_trk = new TH1D( "hMt_trk", "Mt 3 (track & calo iso)", 100, 0., 90.);
30	<	TH1D * hMt_rel_iso_rev = new TH1D( "hMt_rel_iso_rev", "Mt 3 (relative iso reversed)", 100, 0., 90.);
31	<	TH1D * hMt_trk_rev = new TH1D( "hMt_trk_rev", "Mt 3 (track & calo iso reversed)", 100, 0., 90.);
32	<	TH1D * hMt_goodmuon = new TH1D( "hMt_goodmuon", "Mt 3 (good Global Muon)", 100, 0., 90.);
33	<	TH1D * hMt_goodmuon_rev = new TH1D( "hMt_goodmuon_rev", "Mt 3 (good Global Muon reversed)", 100, 0., 90.);
27	>	TH1D * hMt = new TH1D( "hMt", "Mt (all) ", 100, 0., 90.);
28	>	TH1D * hMt_nominal  = new TH1D( "hMt_nominal", "Mt (nominal selection)", 100, 0., 90.);
29	>	TH1D * hMt_reversed = new TH1D( "hMt_reversed", "Mt (reversed selection)", 100, 0., 90.);
30	>
31	>
32	>	TH1D * hht          = new TH1D( "hht", "ht  (all) ", 500, 0., 2000.);
33	>	TH1D * hht_nominal  = new TH1D( "hht_nominal", "ht  (nominal) ", 500, 0., 2000.);
34	>	TH1D * hht_reversed = new TH1D( "hht_reversed", "ht  (reversed) ", 500, 0., 2000.);
35
36	+	TH1D * hM123inv = new TH1D( "hM123inv", "M123inv  (all) ", 500, 0., 4500.);
37	+	TH1D * hM123inv_nominal  = new TH1D( "hM123inv_nominal", "M123inv (nominal selection)", 500, 0., 4500.);
38	+	TH1D * hM123inv_reversed = new TH1D( "hM123inv_reversed", "M123inv (reversed selection)", 500, 0., 4500.);
39	+
40	+	TH1D * hgetHt3          = new TH1D( "hgetHt3", "Ht3 (all) ", 100, 0., 900.);
41	+	TH1D * hgetHt3_nominal  = new TH1D( "hgetHt3_nominal", "Ht3 (nominal selection)", 100, 0., 900.);
42	+	TH1D * hgetHt3_reversed = new TH1D( "hgetHt3_reversed", "Ht3 (reversed selection)", 100, 0., 900.);
43
44	<	TH1D * hgetHt3 = new TH1D( "hgetHt3", "Ht3 (all) ", 100, 0., 900.);
45	<	TH1D * hgetHt3_rel_iso = new TH1D( "hgetHt3_rel_iso", "Ht3 (relative iso)", 100, 0., 900.);
46	<	TH1D * hgetHt3_trk = new TH1D( "hgetHt3_trk", "Ht3 (track & calo iso)", 100, 0., 900.);
47	<	TH1D * hgetHt3_rel_iso_rev = new TH1D( "hgetHt3_rel_iso_rev", "Ht3 (relative iso reversed)", 100, 0., 900.);
48	<	TH1D * hgetHt3_trk_rev = new TH1D( "hgetHt3_trk_rev", "Ht3 (track & calo iso reversed)", 100, 0., 900.);
39	<	TH1D * hgetHt3_goodmuon = new TH1D( "hgetHt3_goodmuon", "Ht3 (good Global muon )", 100, 0., 900.);
40	<	TH1D * hgetHt3_goodmuon_rev = new TH1D( "hgetHt3_goodmuon_rev", "Ht3 (good Global muon reversed)", 100, 0., 900.);
41	<
42	<	TH1D * hktMinPrime = new TH1D( "hktMinPrime", "ktMinPrime (all) ", 100, 0., 25.);
43	<	TH1D * hktMinPrime_rel_iso = new TH1D( "hktMinPrime_rel_iso", "ktMinPrime (relative iso)", 100, 0., 25.);
44	<	TH1D * hktMinPrime_trk = new TH1D( "hktMinPrime_trk", "ktMinPrime (track & calo iso)", 100, 0., 25.);
45	<	TH1D * hktMinPrime_rel_iso_rev= new TH1D( "hktMinPrime_rel_iso_rev", "ktMinPrime (relative iso reversed)", 100, 0., 25.);
46	<	TH1D * hktMinPrime_trk_rev = new TH1D( "hktMinPrime_trk_rev", "ktMinPrime (track & calo iso reversed)", 100, 0., 25.);
47	<	TH1D * hktMinPrime_goodmuon = new TH1D( "hktMinPrime_goodmuon", "ktMinPrime (Good Global muon)", 100, 0., 25.);
48	<	TH1D * hktMinPrime_goodmuon_rev = new TH1D( "hktMinPrime_goodmuon_rev", "ktMinPrime (Good Global muon reversed)", 100, 0., 25.);
44	>
45	>	TH1D * hktMinPrime          = new TH1D( "hktMinPrime", "ktMinPrime (all) ", 100, 0., 25.);
46	>	TH1D * hktMinPrime_nominal  = new TH1D( "hktMinPrime_nominal", "ktMinPrime (nominal selection)", 100, 0., 25.);
47	>	TH1D * hktMinPrime_reversed = new TH1D( "hktMinPrime_reversed", "ktMinPrime (reversed selection)", 100, 0., 25.);
48	>
49
50	<	TH1D * hht2p = new TH1D( "hht2p", "ht2p (all) ", 100, 0., 4.);
51	<	TH1D * hht2p_rel_iso = new TH1D( "hht2p_rel_iso", "ht2p (relative iso)", 100, 0., 4.);
52	<	TH1D * hht2p_trk = new TH1D( "hht2p_trk", "ht2p (track & calo iso)", 100, 0., 4.);
53	<	TH1D * hht2p_rel_iso_rev = new TH1D( "hht2p_rel_iso_rev", "ht2p (relative iso reversed)", 100, 0., 4.);
54	<	TH1D * hht2p_trk_rev = new TH1D( "hht2p_trk_rev", "ht2p (track & calo iso reversed)", 100, 0., 4.);
55	<	TH1D * hht2p_goodmuon = new TH1D( "hht2p_goodmuon", "ht2p (good global muon )", 100, 0., 4.);
56	<	TH1D * hht2p_goodmuon_rev = new TH1D( "hht2p_goodmuon_rev", "ht2p (good global muon reversed)", 100, 0., 4.);
57	<
58	<	TH1D * hDphiJMET = new TH1D( "hDphiJMET", "DphiJMET (all) ", 100, 0., 3.5);
59	<	TH1D * hDphiJMET_rel_iso = new TH1D( "hDphiJMET_rel_iso", "DphiJMET (relative iso)", 100, 0., 3.5);
60	<	TH1D * hDphiJMET_trk = new TH1D( "hDphiJMET_trk", "DphiJMET (track & calo iso)", 100, 0., 3.5);
61	<	TH1D * hDphiJMET_rel_iso_rev = new TH1D( "hDphiJMET_rel_iso_rev", "DphiJMET (relative iso reveresed)", 100, 0., 3.5);
62	<	TH1D * hDphiJMET_trk_rev = new TH1D( "hDphiJMET_trk_rev", "DphiJMET (track & calo iso reversed)", 100, 0., 3.5);
63	<	TH1D * hDphiJMET_goodmuon = new TH1D( "hDphiJMET_goodmuon", "DphiJMET (good global muon )", 100, 0., 3.5);
64	<	TH1D * hDphiJMET_goodmuon_rev = new TH1D( "hDphiJMET_goodmuon_rev", "DphiJMET (good global muon reveresed)", 100, 0., 3.5);
65	<
66	<	TH1D * haplanarity = new TH1D( "haplanarity", "aplanarity (all) ", 100, 0., 0.5);
67	<	TH1D * haplanarity_rel_iso = new TH1D( "haplanarity_rel_iso", "aplanarity (relative iso)", 100, 0., 0.5);
68	<	TH1D * haplanarity_trk = new TH1D( "haplanarity_trk", "aplanarity (track & calo iso)", 100, 0., 0.5);
69	<	TH1D * haplanarity_rel_iso_rev = new TH1D( "haplanarity_rel_iso_rev", "aplanarity (relative iso reversed)", 100, 0., 0.5);
70	<	TH1D * haplanarity_trk_rev = new TH1D( "haplanarity_trk_rev", "aplanarity (track & calo iso reversed)", 100, 0., 0.5);
71	<	TH1D * haplanarity_goodmuon = new TH1D( "haplanarity_goodmuon", "aplanarity (good global muon )", 100, 0., 0.5);
72	<	TH1D * haplanarity_goodmuon_rev = new TH1D( "haplanarity_goodmuon_rev", "aplanarity (good global muon reversed)", 100, 0., 0.5);
73	<
74	<	TH1D * hhtPlusLepton = new TH1D( "hhtPlusLepton", "htPlusLepton (all) ", 100, 0., 20000.);
75	<	TH1D * hhtPlusLepton_rel_iso = new TH1D( "hhtPlusLepton_rel_iso", "htPlusLepton (relative iso)", 100, 0., 20000.);
76	<	TH1D * hhtPlusLepton_trk = new TH1D( "hhtPlusLepton_trk", "htPlusLepton (track & calo iso)", 100, 0., 20000.);
77	<	TH1D * hhtPlusLepton_rel_iso_rev = new TH1D( "hhtPlusLepton_rel_iso_rev", "htPlusLepton (relative iso reversed)", 100, 0., 20000.);
78	<	TH1D * hhtPlusLepton_trk_rev = new TH1D( "hhtPlusLepton_trk_rev", "htPlusLepton (track & calo iso reversed)", 100, 0., 20000.);
79	<	TH1D * hhtPlusLepton_goodmuon = new TH1D( "hhtPlusLepton_goodmuon", "htPlusLepton (good global muon)", 100, 0., 20000.);
80	<	TH1D * hhtPlusLepton_goodmuon_rev = new TH1D( "hhtPlusLepton_goodmuon_rev", "htPlusLepton (good global muon reversed)", 100, 0., 20000.);
81	<
82	<	TH1D * hdPhiLMet = new TH1D( "hdPhiLMet", "dPhiLMet (all) ", 100, 0., 3.5);
83	<	TH1D * hdPhiLMet_rel_iso = new TH1D( "hdPhiLMet_rel_iso", "dPhiLMet (relative iso)", 100, 0., 3.5);
84	<	TH1D * hdPhiLMet_trk = new TH1D( "hdPhiLMet_trk", "dPhiLMet (track & calo iso)", 100, 0., 3.5);
85	<	TH1D * hdPhiLMet_rel_iso_rev = new TH1D( "hdPhiLMet_rel_iso_rev", "dPhiLMet (relative iso reversed)", 100, 0., 3.5);
86	<	TH1D * hdPhiLMet_trk_rev = new TH1D( "hdPhiLMet_trk_rev", "dPhiLMet (track & calo iso reversed)", 100, 0., 3.5);
87	<	TH1D * hdPhiLMet_goodmuon = new TH1D( "hdPhiLMet_goodmuon", "dPhiLMet (good global muon )", 100, 0., 3.5);
88	<	TH1D * hdPhiLMet_goodmuon_rev = new TH1D( "hdPhiLMet_goodmuon_rev", "dPhiLMet (good global muon reversed)", 100, 0., 3.5);
50	>	TH1D * hht2p          = new TH1D( "hht2p", "ht2p (all) ", 100, 0., 4.);
51	>	TH1D * hht2p_nominal  = new TH1D( "hht2p_nominal", "ht2p (nominal selection)", 100, 0., 4.);
52	>	TH1D * hht2p_reversed = new TH1D( "hht2p_reversed", "ht2p (reversed selection)", 100, 0., 4.);
53	>
54	>	TH1D * hcentrality          = new TH1D( "hcentrality", "centrality (all) ", 100, 0., 1.);
55	>	TH1D * hcentrality_nominal  = new TH1D( "hcentrality_nominal", "centrality (nominal selection)", 100, 0., 1.);
56	>	TH1D * hcentrality_reversed = new TH1D( "hcentrality_reversed", "centrality (reversed selection)", 100, 0., 1.);
57	>
58	>	TH1D * hsphericity          = new TH1D( "hsphericity", "sphericity (all) ", 100, 0., 1.);
59	>	TH1D * hsphericity_nominal  = new TH1D( "hsphericity_nominal", "sphericity (nominal selection)", 100, 0., 1.);
60	>	TH1D * hsphericity_reversed = new TH1D( "hsphericity_reversed", "sphericity (reversed selection)", 100, 0., 1.);
61	>
62	>	TH1D * hhtPlusLepton = new TH1D( "hhtPlusLepton", "htPlusLepton (all) ", 500, 0., 3000.);
63	>	TH1D * hhtPlusLepton_nominal  = new TH1D( "hhtPlusLepton_nominal", "htPlusLepton (nominal selection)", 500, 0., 3000.);
64	>	TH1D * hhtPlusLepton_reversed = new TH1D( "hhtPlusLepton_reversed", "htPlusLepton (reversed selection)", 500, 0., 3000.);
65	>
66	>	TH1D * hdPhiLMet          = new TH1D( "hdPhiLMet", "dPhiLMet (all) ", 100, 0., 3.5);
67	>	TH1D * hdPhiLMet_nominal  = new TH1D( "hdPhiLMet_nominal", "dPhiLMet (nominal selection)", 100, 0., 3.5);
68	>	TH1D * hdPhiLMet_reversed = new TH1D( "hdPhiLMet_reversed", "dPhiLMet (reversed selection)", 100, 0., 3.5);
69	>
70	>	TH1D * hDphiJMET          = new TH1D( "hDphiJMET", "DphiJMET (all) ", 100, 0., 3.5);
71	>	TH1D * hDphiJMET_nominal  = new TH1D( "hDphiJMET_nominal", "DphiJMET (nominal selection)", 100, 0., 3.5);
72	>	TH1D * hDphiJMET_reversed = new TH1D( "hDphiJMET_reversed", "DphiJMET (reversed selection)", 100, 0., 3.5);
73	>
74	>
75	>	TH1D * haplanarity          = new TH1D( "haplanarity", "aplanarity (all) ", 100, 0., 0.5);
76	>	TH1D * haplanarity_nominal  = new TH1D( "haplanarity_nominal", "aplanarity (nominal selection)", 100, 0., 0.5);
77	>	TH1D * haplanarity_reversed = new TH1D( "haplanarity_reversed", "aplanarity (reversed selection)", 100, 0., 0.5);
78	>
79	>	TH1D * hW_MT          = new TH1D( "hW_MT", "W_MT (all) ", 500, 0., 5000);
80	>	TH1D * hW_MT_nominal  = new TH1D( "hW_MT_nominal", "W_MT (nominal selection)", 500, 0., 5000);
81	>	TH1D * hW_MT_reversed = new TH1D( "hW_MT_reversed", "W_MT (reversed selection)", 500, 0., 5000);
82	>
83	>	TH1D * hmetHtPlusLepton = new TH1D( "hmetHtPlusLepton", "metHtPlusLepton (all) ", 500, 0., 4000.);
84	>	TH1D * hmetHtPlusLepton_nominal  = new TH1D( "hmetHtPlusLepton_nominal", "metHtPlusLepton (nominal selection)", 500, 0., 4000.);
85	>	TH1D * hmetHtPlusLepton_reversed = new TH1D( "hmetHtPlusLepton_reversed", "metHtPlusLepton (reversed selection)", 500, 0., 4000.);
86	>
87	>
88
89		TFile *fin = TFile::Open(infile);
90	<	TTree * t1 = (TTree*)fin->Get("qcd");
90	>	TTree * t1 = (TTree*)fin->Get("ttljets");
91
92	<	Double_t Mt, lepton_track_iso, lepton_calo_iso;
93	<	Double_t getHt3, ktMinPrime, ht2p,aplanarity, htPlusLepton,dPhiLMet,DphiJMET,lepton_pt;
94	<	Bool_t muon_isGood_AllGlobalMuons,muon_isGood_AllStandAloneMuons;
95	<
96	<	t1->SetBranchAddress("getMt",&Mt);
92	>	Double_t lepton_track_iso, lepton_calo_iso;
93	>	Double_t ht, Mt, M123inv,getHt3, ktMinPrime, ht2p,centrality,aplanarity, sphericity,htPlusLepton;
94	>	Double_t dPhiLMet,DphiJMET,W_MT,metHtPlusLepton,leptonJetDeltaR,muon_track_Chi2,muon_track_nChi2,muon_track_ndof,lepton_pt;
95	>	Bool_t muon_isGood_AllGlobalMuons,muon_isGood_TMLastStationTight,muon_isGood_TMLastStationLoose;
96	>	Float_t eidRobustLoose,eidRobustTight;
97	>	Bool_t muon_isGood_TM2DCompatibilityTight,muon_isGood_TM2DCompatibilityLoose;
98	>
99		t1->SetBranchAddress("lepton_track_iso",&lepton_track_iso);
100		t1->SetBranchAddress("lepton_calo_iso",&lepton_calo_iso);
101	+	t1->SetBranchAddress("ht",&ht);
102	+	t1->SetBranchAddress("getMt",&Mt);
103	+	t1->SetBranchAddress("getM123inv",&M123inv);
104		t1->SetBranchAddress("getHt3",&getHt3);
105		t1->SetBranchAddress("ktMinPrime",&ktMinPrime);
106		t1->SetBranchAddress("ht2p",&ht2p);
107	<	t1->SetBranchAddress("aplanarity",&aplanarity);
107	>	t1->SetBranchAddress("centrality",&centrality);
108	>	t1->SetBranchAddress("aplanarity",&aplanarity);
109	>	t1->SetBranchAddress("sphericity",&sphericity);
110		t1->SetBranchAddress("htPlusLepton",&htPlusLepton);
111		t1->SetBranchAddress("dPhiLMet",&dPhiLMet);
112		t1->SetBranchAddress("DphiJMET",&DphiJMET);
113	+	t1->SetBranchAddress("W_MT",&W_MT);
114	+	t1->SetBranchAddress("metHtPlusLepton",&metHtPlusLepton);
115		t1->SetBranchAddress("muon_isGood_AllGlobalMuons",&muon_isGood_AllGlobalMuons);
116	<	t1->SetBranchAddress("muon_isGood_AllStandAloneMuons",&muon_isGood_AllStandAloneMuons);
116	>	t1->SetBranchAddress("leptonJetDeltaR",&leptonJetDeltaR);
117	>	t1->SetBranchAddress("muon_track_nChi2",&muon_track_nChi2);
118	>	t1->SetBranchAddress("muon_track_ndof",&muon_track_ndof);
119	>	t1->SetBranchAddress("muon_track_chi2",&muon_track_Chi2);
120	>	t1->SetBranchAddress("muon_isGood_TMLastStationTight",&muon_isGood_TMLastStationTight);
121	>	t1->SetBranchAddress("muon_isGood_TMLastStationLoose",&muon_isGood_TMLastStationLoose);
122	>	t1->SetBranchAddress("muon_isGood_TM2DCompatibilityTight",&muon_isGood_TM2DCompatibilityTight);
123	>	t1->SetBranchAddress("muon_isGood_TM2DCompatibilityLoose",&muon_isGood_TM2DCompatibilityLoose);
124		t1->SetBranchAddress("lepton_pt",&lepton_pt);
125	<
125	>	t1->SetBranchAddress("eidRobustLoose",&eidRobustLoose);
126	>	t1->SetBranchAddress("eidRobustTight",&eidRobustTight);
127
128
129		Int_t n_gen = (Int_t)t1->GetEntries();
130		cout<<"No of entries from trees "<<n_gen<<endl;
131	+	bool leptonid_cut;
132	+
133
134		for(Int_t i=0; i<n_gen; i++){
135		t1->GetEntry(i);
136
137	+	hht->Fill(ht);
138		hMt->Fill(Mt);
139	+	hM123inv->Fill(M123inv);
140		hgetHt3->Fill(getHt3);
141	<	hktMinPrime->Fill(ktMinPrime);
142	<	hht2p->Fill(ht2p);
143	<	hDphiJMET->Fill(DphiJMET);
144	<	haplanarity->Fill(aplanarity);
145	<	hhtPlusLepton->Fill(htPlusLepton);
146	<	hdPhiLMet->Fill(dPhiLMet);
147	<
148	<	if(1.0/(1.0+(lepton_calo_iso+lepton_track_iso)/lepton_pt) > 0.67 ){
149	<	hMt_rel_iso->Fill(Mt);
150	<	hgetHt3_rel_iso->Fill(getHt3);
151	<	hktMinPrime_rel_iso->Fill(ktMinPrime);
152	<	hht2p_rel_iso->Fill(ht2p);
153	<	hDphiJMET_rel_iso->Fill(DphiJMET);
154	<	haplanarity_rel_iso->Fill(aplanarity);
155	<	hhtPlusLepton_rel_iso->Fill(htPlusLepton);
156	<	hdPhiLMet_rel_iso->Fill(dPhiLMet);
157	<	}
141	>	hktMinPrime->Fill(ktMinPrime);
142	>	hht2p->Fill(ht2p);
143	>	haplanarity->Fill(aplanarity);
144	>	hcentrality->Fill(centrality);
145	>	hsphericity->Fill(sphericity);
146	>	hdPhiLMet->Fill(dPhiLMet);
147	>	hDphiJMET->Fill(DphiJMET);
148	>	hW_MT->Fill(W_MT);
149	>	hhtPlusLepton->Fill(htPlusLepton);
150	>	hmetHtPlusLepton->Fill(metHtPlusLepton);
151	>	//leptonid_cut = (muon_isGood_TM2DCompatibilityLoose && (1.0/(1.0+(lepton_calo_iso+lepton_track_iso)/lepton_pt)) > 0.67 );
152	>	//leptonid_cut = (muon_isGood_TMLastStationLoose && muon_track_nChi2 < 5 );
153	>	leptonid_cut = (eidRobustTight == 1);
154	>	if(leptonid_cut){
155	>	hht_nominal->Fill(ht);
156	>	hMt_nominal->Fill(Mt);
157	>	hM123inv_nominal->Fill(M123inv);
158	>	hgetHt3_nominal->Fill(getHt3);
159	>	hktMinPrime_nominal->Fill(ktMinPrime);
160	>	hht2p_nominal->Fill(ht2p);
161	>	haplanarity_nominal->Fill(aplanarity);
162	>	hcentrality_nominal->Fill(centrality);
163	>	hsphericity_nominal->Fill(sphericity);
164	>	hdPhiLMet_nominal->Fill(dPhiLMet);
165	>	hDphiJMET_nominal->Fill(DphiJMET);
166	>	hW_MT_nominal->Fill(W_MT);
167	>	hhtPlusLepton_nominal->Fill(htPlusLepton);
168	>	hmetHtPlusLepton_nominal->Fill(metHtPlusLepton);
169	>	}
170		//     hlepton_iso->Fill(lepton_calo_iso);
171	<	if(lepton_calo_iso<7 && lepton_track_iso < 3 && muon_isGood_AllStandAloneMuons){
172	<	hMt_trk->Fill(Mt);
173	<	hgetHt3_trk->Fill(getHt3);
174	<	hktMinPrime_trk->Fill(ktMinPrime);
175	<	hht2p_trk->Fill(ht2p);
176	<	hDphiJMET_trk->Fill(DphiJMET);
177	<	haplanarity_trk->Fill(aplanarity);
178	<	hhtPlusLepton_trk->Fill(htPlusLepton);
179	<	hdPhiLMet_trk->Fill(dPhiLMet); }
180	<	if(1.0/(1.0+(lepton_calo_iso+lepton_track_iso)/lepton_pt) < 0.67 ){
181	<	hMt_rel_iso_rev->Fill(Mt);
182	<	hgetHt3_rel_iso_rev->Fill(getHt3);
183	<	hktMinPrime_rel_iso_rev->Fill(ktMinPrime);
184	<	hht2p_rel_iso_rev->Fill(ht2p);
185	<	hDphiJMET_rel_iso_rev->Fill(DphiJMET);
186	<	haplanarity_rel_iso_rev->Fill(aplanarity);
155	<	hhtPlusLepton_rel_iso_rev->Fill(htPlusLepton);
156	<	hdPhiLMet_rel_iso_rev->Fill(dPhiLMet);
157	<	}
158	<	//     hlepton_iso->Fill(lepton_calo_iso);
159	<	if(lepton_calo_iso>7 && lepton_track_iso > 3 && muon_isGood_AllStandAloneMuons){
160	<	hMt_trk_rev->Fill(Mt);
161	<	hgetHt3_trk_rev->Fill(getHt3);
162	<	hktMinPrime_trk_rev->Fill(ktMinPrime);
163	<	hht2p_trk_rev->Fill(ht2p);
164	<	hDphiJMET_trk_rev->Fill(DphiJMET);
165	<	haplanarity_trk_rev->Fill(aplanarity);
166	<	hhtPlusLepton_trk_rev->Fill(htPlusLepton);
167	<	hdPhiLMet_trk_rev->Fill(dPhiLMet); }
168	<	if(muon_isGood_AllGlobalMuons){
169	<	hMt_goodmuon->Fill(Mt);
170	<	hgetHt3_goodmuon->Fill(getHt3);
171	<	hktMinPrime_goodmuon->Fill(ktMinPrime);
172	<	hht2p_goodmuon->Fill(ht2p);
173	<	hDphiJMET_goodmuon->Fill(DphiJMET);
174	<	haplanarity_goodmuon->Fill(aplanarity);
175	<	hhtPlusLepton_goodmuon->Fill(htPlusLepton);
176	<	hdPhiLMet_goodmuon->Fill(dPhiLMet); }
177	<	if(!muon_isGood_AllGlobalMuons){
178	<	hMt_goodmuon_rev->Fill(Mt);
179	<	hgetHt3_goodmuon_rev->Fill(getHt3);
180	<	hktMinPrime_goodmuon_rev->Fill(ktMinPrime);
181	<	hht2p_goodmuon_rev->Fill(ht2p);
182	<	hDphiJMET_goodmuon_rev->Fill(DphiJMET);
183	<	haplanarity_goodmuon_rev->Fill(aplanarity);
184	<	hhtPlusLepton_goodmuon_rev->Fill(htPlusLepton);
185	<	hdPhiLMet_goodmuon_rev->Fill(dPhiLMet); }
171	>	if(!leptonid_cut){
172	>	hht_reversed->Fill(ht);
173	>	hMt_reversed->Fill(Mt);
174	>	hM123inv_reversed->Fill(M123inv);
175	>	hgetHt3_reversed->Fill(getHt3);
176	>	hktMinPrime_reversed->Fill(ktMinPrime);
177	>	hht2p_reversed->Fill(ht2p);
178	>	haplanarity_reversed->Fill(aplanarity);
179	>	hcentrality_reversed->Fill(centrality);
180	>	hsphericity_reversed->Fill(sphericity);
181	>	hdPhiLMet_reversed->Fill(dPhiLMet);
182	>	hDphiJMET_reversed->Fill(DphiJMET);
183	>	hW_MT_reversed->Fill(W_MT);
184	>	hhtPlusLepton_reversed->Fill(htPlusLepton);
185	>	hmetHtPlusLepton_reversed->Fill(metHtPlusLepton);
186	>	}
187
188
189		}
189	–
190	–
191	–
190
191		outf->Write();
192		outf->Close();

Diff Legend

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