ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/andis/plotsqcd.c
Revision: 1.1
Committed: Tue Oct 20 09:08:16 2009 UTC (15 years, 6 months ago) by andis
Content type: text/plain
Branch: MAIN
Log Message: 
*** empty log message ***

File Contents

# User Rev Content
1 andis 1.1 #include "TCanvas.h"
2     #include "TROOT.h"
3     #include "TBrowser.h"
4     #include "TFile.h"
5     #include "TNtuple.h"
6     #include "TStyle.h"
7     #include "TFrame.h"
8     #include "TH1.h"
9     #include "TH2.h"
10     #include "TTree.h"
11     #include "TCut.h"
12     #include "TEntryList.h"
13     #include "TMath.h"
14     #include <iostream>
15     #include <fstream>
16     #include <string>
17    
18     using namespace std;
19    
20     //string filename("crab_0_091008_103310/res/all");
21    
22     // CLASSES
23    
24     class Leptons
25     {
26     public:
27     void Init_Leptons();
28     Double_t pdgid, recomuons, recoelectrons, ptgen, etagen,
29     muon1_pt, muon1_eta, muon1_phi, muon1_ep, muon1_ep_mod, muon1_hovere, muon1_bjetangle, muon1_nearestjetangle, muon1_iso_03,
30     muon1_iso_05_emEt, muon1_ipt, muon1_ipt_sigma, muon1_ip3D, muon1_ip3D_sigma,
31     electron1_pt, electron1_eta, electron1_phi, electron1_deltaeta, electron1_deltaphi, electron1_sigmaetaeta, electron1_ep,
32     electron1_ep_mod, electron1_hovere, electron1_bjetangle, electron1_nearestjetangle, electron1_iso_03, electron1_iso_04_emEt,
33     electron1_ipt, electron1_ipt_sigma, electron1_ip3D, electron1_ip3D_sigma;
34     bool addlepton,muon1_global;
35    
36     };
37    
38     Leptons leptons;
39    
40     class Met
41     {
42     public:
43     void Init_Met();
44     Double_t tc_met, calo_met, et_neutrinos, tc_met_lepton, calo_met_lepton;
45     };
46    
47     Met met;
48    
49     class Pfjets
50     {
51     public:
52     void Init_Pfjets();
53     Double_t pt1,pt2,pt3,phi1,phi2,phi3,eta1,eta2,eta3,
54     btag1,btag2,btag3,hardb_pt,hardbjet_pt,hardbjet_deltar,softb_pt,softbjet_pt,softbjet_deltar,deltar,
55     nearestjet_hcal, nearestjet_ntracks;
56     };
57    
58     Pfjets pfjets;
59    
60     class Misc
61     {
62     public:
63     void Init_Misc();
64     Double_t pthiggs,met_higgs_angle,calomet_higgs_angle,higgs_pt,higgs_invmass,higgs_eta,higgs_phi,w_pt,w_invmass,w_eta,w_phi;
65     };
66    
67     Misc misc;
68    
69     // CUTS
70    
71     void leptons_pt_iso_cut(const Double_t pt, const Double_t isoel, const Double_t isomu, TTree *maintree);
72     void leptons_isgood_cut(const Double_t deltaeta, const Double_t deltaphi, const Double_t sigmaetaeta, const Double_t ep, const Double_t hovere, TTree *maintree);
73     void leptons_bjetangle_emEt_cut(const Double_t min_bjet_angle, const Double_t max_bjet_angle, const Double_t muon_emEt_05, TTree *maintree);
74     void leptons_jets_cut(const Double_t ntracks, const Double_t hcal, const Double_t max_delta_r, TTree *maintree);
75     void pfjets_btag_cut(const Double_t btag, TTree *maintree);
76     void pfjets_pt_cut(const Double_t pt1, const Double_t pt2, TTree *maintree);
77     void w_higgs_angle_cut(const Double_t angle, TTree *maintree);
78     void higgs_invmass_cut(const Double_t invmass, TTree *maintree);
79     void higgs_pt_cut(const Double_t pt, TTree *maintree);
80     void pfjets_angle_cut(const Double_t angle, TTree *maintree);
81     //
82    
83     Double_t deltaR(Double_t eta1, Double_t phi1, Double_t eta2, Double_t phi2);
84    
85     //
86    
87     TFile *savefile = 0;
88     TFile *file = 0;
89     TTree *mytree = 0;
90     UInt_t totalevents,nentries = 0;
91     Double_t selentries = 0;
92    
93     TObjArray Hlist(0);
94     ofstream outfile;
95    
96     void plots(string filename)
97     {
98    
99     gROOT->SetStyle("Plain");
100     gROOT->ForceStyle();
101     gStyle->SetOptStat(111111);
102     gStyle->SetPalette(1);
103    
104     cout << filename + ".txt" << endl;
105    
106     outfile.open((filename + ".txt").c_str());
107     outfile << "preselection: pt(pfjets) > 20" << endl;
108     outfile.flush();
109    
110     cout << filename + "_histograms.root" << endl;
111    
112     savefile = new TFile((filename + "_histograms.root").c_str(),"RECREATE");
113    
114     cout << filename + ".root" << endl;
115    
116     file = new TFile((filename + ".root").c_str(),"READ");
117     if(file->IsZombie()) exit(-1);
118     TTree *mytree = (TTree*)file->Get("calo/recotree");
119     totalevents = mytree->GetEntries();
120     cout << totalevents << endl;
121    
122     mytree->SetBranchAddress("misc",&misc);
123     mytree->SetBranchAddress("leptons",&leptons);
124     mytree->SetBranchAddress("met",&met);
125     mytree->SetBranchAddress("pfjets",&pfjets);
126    
127     //
128    
129     TCanvas *cmisc = new TCanvas("cmisc","Miscellaneous");
130     TCanvas *cleptons = new TCanvas("cleptons","Leptons");
131     TCanvas *celectrons = new TCanvas("celectrons","Electrons");
132     TCanvas *cmuons = new TCanvas("cmuons","Muons");
133     TCanvas *cmet = new TCanvas("cmet","MET");
134     TCanvas *cjets = new TCanvas("cjets","Jets");
135     TCanvas *cuts1 = new TCanvas("cuts1","Cuts 1");
136     TCanvas *cuts2 = new TCanvas("cuts2","Cuts 2");
137     TCanvas *cuts3 = new TCanvas("cuts3","Cuts 3");
138    
139     //
140    
141     TH1D * pfjets_deltar = new TH1D("pfjets_deltar","\\Delta R of pfjets",100,0,7);
142     TH1D * ptlepton_histo = new TH1D("ptlepton_histo","p_{T} of lepton candidate",100,0,500);
143     TH1D * lepton_bjet_angle = new TH1D("lepton_bjet_angle","angle between lepton candidate and nearest btagged jet",100,0,7);
144     TH1D * higgs_invmass = new TH1D("higgs_invmass","inv. mass of higgs candidate",100,0,400);
145     TH2D * higgs_invmass_pt = new TH2D("higgs_invmass_pt","inv. mass vs. p_{T} of higgs candidate",100,0,600,100,0,600);
146    
147     TH2D * calo_met_angle_higgs_pt = new TH2D("calo_met_angle_higgs_t","cos(angle) between calo_MET and Higgs vs. Higgs p_{T}",100,-1.25,1.25,100,0,700);
148     TH2D * tc_met_angle_higgs_pt = new TH2D("tc_met_angle_higgs_pt","cos(angle) between tk_MET and Higgs vs. Higgs p_{T}",100,-1.25,1.25,100,0,500);
149     TH2D * calo_met_angle_higgs_invmass = new TH2D("calo_met_angle_higgs_invmass","cos(angle) between calo_MET/W and Higgs vs. inv. mass of Higgs",50,-1.25,1.25,100,100,500);
150     TH2D * tc_met_angle_higgs_invmass = new TH2D("tc_met_angle_higgs_invmass","cos(angle) between tk_MET/W and Higgs vs. inv. mass of Higgs",50,-1.25,1.25,100,100,500);
151    
152     cmisc->Divide(3,1);
153     cleptons->Divide(3,2);
154     celectrons->Divide(3,2);
155     cmuons->Divide(3,2);
156     cmet->Divide(3,2);
157     cjets->Divide(3,2);
158     cuts1->Divide(3,2);
159     cuts2->Divide(3,2);
160     cuts3->Divide(3,2);
161    
162     // FILL HISTOS
163    
164     /*
165     for(UInt_t i = 0; i < totalevents; ++i)
166     {
167     mytree->GetEntry(i);
168     calo_met_angle_higgs_pt->Fill(misc.calo_met_higgs_angle,misc.higgs_pt);
169     tc_met_angle_higgs_pt->Fill(misc.tc_met_higgs_angle,misc.higgs_pt);
170     calo_met_angle_higgs_invmass->Fill(misc.calo_met_higgs_angle,misc.higgs_invmass);
171     tc_met_angle_higgs_invmass->Fill(misc.tc_met_higgs_angle,misc.higgs_invmass);
172     calo_met_included_lepton->Fill(met.calo_met,met.calo_met_lepton);
173     tc_met_included_lepton->Fill(met.tc_met,met.tc_met_lepton);
174    
175     higgs_invmass->Fill(misc.higgs_invmass);
176     higgs_invmass_pt->Fill(misc.higgs_invmass,misc.higgs_pt);
177    
178     //if(pfjets.pt1 > 60 && pfjets.pt2 > 60 && pfjets.btag1 > 0.3 && pfjets.btag2 > 0.3 && (leptons.muon1_pt < leptons.electron1_pt && leptons.electron1_pt > 30 && leptons.electron1_iso_03 < 10 || leptons.muon1_pt > leptons.electron1_pt && leptons.muon1_pt > 30 && leptons.muon1_iso_03 < 10))
179    
180     }
181     */
182    
183     // SAVE HISTOS
184    
185     /*Hlist.Add(calo_met_angle_higgs_pt);
186     Hlist.Add(calo_met_angle_higgs_invmass);
187     Hlist.Add(higgs_invmass);
188     Hlist.Add(higgs_invmass_pt);
189     */
190     // DRAW HISTOS
191    
192     cmisc->cd(1);
193     mytree->Draw("misc.pthiggs","misc.pthiggs>200&&leptons.addlepton==0&&(leptons.muon1_pt>25&&leptons.muon1_iso03<15 || leptons.electron1_pt>25&&leptons.electron1_iso03<15)");
194     cmisc->cd(2);
195     mytree->Draw("misc.higgs_pt","misc.higgs_pt>180&&leptons.addlepton==0&&(leptons.muon1_pt>25&&leptons.muon1_iso03<15 || leptons.electron1_pt>25&&leptons.electron1_iso03<15)");
196     cmisc->cd(3);
197     mytree->Draw("pfjets.deltar>>pfdr(100,0,4)");
198    
199     cleptons->cd(1);
200     mytree->Draw("leptons.muon1_global>>mg1(10,-0.5,1.5)","leptons.muon1_pt > 25");
201     cleptons->cd(2);
202     mytree->Draw("leptons.electron1_deltaeta:leptons.electron1_deltaphi>>eph1(100,0,0.25,100,0,0.04)","","colz");
203     cleptons->cd(3);
204     mytree->Draw("leptons.electron1_nearestjetangle>>enja(100,0,4)");
205     cleptons->cd(4);
206     mytree->Draw("leptons.muon1_nearestjetangle>>mnja(100,0,4)");
207    
208     cjets->cd(1);
209     mytree->Draw("pfjets.nearestjet_hcal");
210     cjets->cd(2);
211     mytree->Draw("pfjets.nearestjet_ntracks");
212    
213     cmuons->cd(1);
214     mytree->Draw("leptons.muon1_iso_05_emEt>>miee(100,0,10)","leptons.muon1_global == 1");
215     cmuons->cd(2);
216     mytree->Draw("leptons.muon1_bjetangle>>mbja(100,0,4)");
217     cmuons->cd(3);
218     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ipt>>njmi(100,-10,10,100,0,6)","","colz");
219     cmuons->cd(4);
220     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ipt_sigma>>njmis(100,-300,300,100,0,6)","","colz");
221     cmuons->cd(5);
222     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ip3D>>njmi3(100,0,15,100,0,6)","","colz");
223     cmuons->cd(6);
224     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ip3D_sigma>>njmi3s(100,0,700,100,0,6)","","colz");
225    
226     celectrons->cd(1);
227     mytree->Draw("leptons.electron1_ep>>eep(100,0,10)");
228     celectrons->cd(2);
229     mytree->Draw("leptons.electron1_bjetangle>>ebja(100,0,4)");
230     celectrons->cd(3);
231     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ipt>>njei(100,-10,10,100,0,6)","","colz");
232     celectrons->cd(4);
233     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ipt_sigma>>njeis(100,-300,300,100,0,6)","","colz");
234     celectrons->cd(5);
235     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ip3D>>njei3(100,0,15,100,0,6)","","colz");
236     celectrons->cd(6);
237     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ip3D_sigma>>njei3s(100,0,700,100,0,6)","","colz");
238    
239     /*
240     cjets->cd(1);
241     mytree->Draw("higgs_eta:higgs_pt","","colz");
242    
243     cmisc->cd(1);
244     mytree->Draw("misc.higgs_invmass","misc.higgs_eta < 2.4 && misc.higgs_eta > -2.4"); //higgs_invmass->Draw();
245     cmisc->cd(2);
246     higgs_invmass_pt->Draw("colz");
247     cmisc->cd(3);
248     calo_met_angle_higgs_pt->Draw("colz");
249     cmisc->cd(4);
250     calo_met_angle_higgs_invmass->Draw("colz");
251     cmisc->cd(5);
252     tc_met_angle_higgs_pt->Draw("colz");
253     cmisc->cd(6);
254     tc_met_angle_higgs_invmass->Draw("colz");
255    
256    
257     cmet->cd(1);
258     calo_met_included_lepton->Draw("colz");
259     cmet->cd(2);
260     tc_met_included_lepton->Draw("colz");
261     cmet->cd(3);
262     mytree->Draw("misc.calo_met_higgs_angle");
263     cmet->cd(4);
264     mytree->Draw("misc.tc_met_higgs_angle");
265     cmet->cd(5);
266     mytree->Draw("met.tc_met_lepton:met.calo_met","","colz");
267    
268    
269     cleptons->cd(1);
270     ptlepton_histo->Draw();
271     cleptons->cd(3);
272     */
273    
274     // ptlepton_histo_cut->Draw("same");
275    
276     // CUTS
277    
278     savefile->cd();
279    
280     leptons_pt_iso_cut(25,15,15,mytree);
281     higgs_pt_cut(150,mytree);
282     leptons_isgood_cut(0.01,0.1,0.03,3,0.1,mytree);
283     leptons_bjetangle_emEt_cut(2,3.5,1,mytree);
284     leptons_jets_cut(12,50,0.5,mytree);
285     pfjets_btag_cut(0.7,mytree);
286     pfjets_pt_cut(50,20,mytree);
287     pfjets_angle_cut(2,mytree);
288     w_higgs_angle_cut(0,mytree);
289    
290    
291     // HISTOS AFTER CUTS
292    
293     cout << "remaining events: " << (mytree->GetEntryList())->GetN() << endl;
294    
295     /* for(UInt_t k = 0; k < (mytree->GetEntryList())->GetN(); ++k)
296     {
297     TEntryList *el = mytree->GetEntryList();
298     mytree->GetEntry(el->GetEntry(k));
299     calo_met_angle_higgs_pt->Fill(misc.calomet_higgs_angle,misc.higgs_pt);
300     tc_met_angle_higgs_pt->Fill(misc.met_higgs_angle,misc.higgs_pt);
301     calo_met_angle_higgs_invmass->Fill(misc.calomet_higgs_angle,misc.higgs_invmass);
302     tc_met_angle_higgs_invmass->Fill(misc.met_higgs_angle,misc.higgs_invmass);
303     higgs_invmass->Fill(misc.higgs_invmass);
304     higgs_invmass_pt->Fill(misc.higgs_invmass,misc.higgs_pt);
305     //pfjets_deltar->Fill(deltaR(pfjets.eta1,pfjets.phi1,pfjets.eta2,pfjets.phi2));
306    
307     if(leptons.addlepton == 0)
308     {
309     if(leptons.electron1_pt > 25)
310     {
311     ptlepton_histo->Fill(leptons.electron1_pt,misc.higgs_pt);
312     Double_t delta1 = deltaR(leptons.electron1_eta,leptons.electron1_phi,pfjets.eta1,pfjets.phi1);
313     Double_t delta2 = deltaR(leptons.electron1_eta,leptons.electron1_phi,pfjets.eta2,pfjets.phi2);
314     if(delta1 > delta2)
315     lepton_bjet_angle->Fill(delta1);
316     else
317     lepton_bjet_angle->Fill(delta2);
318     }
319     else if(leptons.muon1_pt > 25)
320     {
321     ptlepton_histo->Fill(leptons.muon1_pt,misc.higgs_pt);
322     Double_t delta1 = deltaR(leptons.muon1_eta,leptons.muon1_phi,pfjets.eta1,pfjets.phi1);
323     Double_t delta2 = deltaR(leptons.muon1_eta,leptons.muon1_phi,pfjets.eta2,pfjets.phi2);
324     if(delta1 > delta2)
325     lepton_bjet_angle->Fill(delta1);
326     else
327     lepton_bjet_angle->Fill(delta2);
328     }
329     else
330     lepton_bjet_angle->Fill(-10);
331     }
332     }
333     */
334    
335     cjets->cd(3);
336     mytree->Draw("pfjets.deltar");
337     cjets->cd(4);
338     mytree->Draw("pfjets.nearestjet_hcal");
339     cjets->cd(5);
340     mytree->Draw("pfjets.nearestjet_ntracks");
341    
342     cmet->cd(1);
343     mytree->Draw("misc.met_higgs_angle>>mha(100,-1.5,1.5)");
344     cmet->cd(2);
345     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ipt>>njmi2(100,-10,10,100,0,6)","","colz");
346     cmet->cd(3);
347     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ipt_sigma>>njmis2(100,-300,300,100,0,6)","","colz");
348     cmet->cd(4);
349     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ip3D>>njmi3_2(100,0,10,100,0,6)","","colz");
350     cmet->cd(5);
351     mytree->Draw("leptons.muon1_nearestjetangle:leptons.muon1_ip3D_sigma>>njmi3s_2(100,0,700,100,0,6)","","colz");
352    
353     /*cuts2->cd(2);
354     mytree->Draw("pfjets.hardb_pt:pfjets.hardbjet_pt","misc.higgs_pt>200","colz");
355     cuts2->cd(3);
356     mytree->Draw("pfjets.softb_pt:pfjets.softbjet_pt","misc.higgs_pt>200","colz");
357     cuts2->cd(4);
358     mytree->Draw("pfjets.hardb_pt:pfjets.hardbjet_deltar","misc.higgs_pt>200","colz");
359     cuts2->cd(5);
360     mytree->Draw("pfjets.softb_pt:pfjets.softbjet_deltar","misc.higgs_pt>200","colz");*/
361    
362     cuts1->cd(1);
363     mytree->Draw("leptons.electron1_hovere>>ehoe2(100,0,1)");
364     cuts1->cd(2);
365     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ipt>>njei2(100,-10,10,100,0,6)","","colz");
366     cuts1->cd(3);
367     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ipt_sigma>>njeis2(100,-300,300,100,0,6)","","colz");
368     cuts1->cd(4);
369     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ip3D>>njei3_2(100,0,10,100,0,6)","","colz");
370     cuts1->cd(5);
371     mytree->Draw("leptons.electron1_nearestjetangle:leptons.electron1_ip3D_sigma>>njei3s_2(100,0,700,100,0,6)","","colz");
372    
373     cuts2->cd(1);
374     mytree->Draw("leptons.muon1_iso_05_emEt>>miee(100,0,10)","leptons.muon1_global == 1");
375     cuts2->cd(2);
376     mytree->Draw("leptons.muon1_bjetangle>>mbja2(100,0,4)");
377     cuts2->cd(3);
378     mytree->Draw("leptons.electron1_bjetangle>>ebja2(100,0,4)");
379     cuts2->cd(4);
380     mytree->Draw("leptons.electron1_ep>>epm2(100,0,5)");
381     cuts2->cd(5);
382     mytree->Draw("leptons.muon1_nearestjetangle>>mnja2(100,0,4)");
383     cuts2->cd(6);
384     mytree->Draw("leptons.electron1_nearestjetangle>>enja2(100,0,4)");
385    
386    
387     /*
388     cuts1->cd(1);
389     mytree->Draw("leptons.muon1_iso_03","misc.pthiggs > 200 && leptons.pdgid != 0 && ((leptons.muon1_pt > leptons.electron1_pt && ((leptons.muon2_pt > leptons.electron1_pt && leptons.muon2_iso_03 < 5) || (leptons.muon2_pt < leptons.electron1_pt && leptons.electron1_iso_03 < 5))) || (leptons.muon1_pt < leptons.electron1_pt && ((leptons.muon1_pt > leptons.electron2_pt && leptons.muon1_iso_03 < 5) || (leptons.muon1_pt < leptons.electron2_pt && leptons.electron2_iso_03 < 5))))");
390     cuts1->cd(2);
391     mytree->Draw("leptons.muon1_iso_03","misc.pthiggs > 200 && leptons.pdgid != 0 && ((leptons.muon1_pt > leptons.electron1_pt && ((leptons.muon2_pt > leptons.electron1_pt && leptons.muon2_iso_03 < 10) || (leptons.muon2_pt < leptons.electron1_pt && leptons.electron1_iso_03 < 10))) || (leptons.muon1_pt < leptons.electron1_pt && ((leptons.muon1_pt > leptons.electron2_pt && leptons.muon1_iso_03 < 10) || (leptons.muon1_pt < leptons.electron2_pt && leptons.electron2_iso_03 < 10))))");
392    
393     cuts1->cd(1);
394     mytree->Draw("ptjets.pt1","ptjets.btag1 > 0.1 && ptjets.btag2 > 0.1 && jets.pt3 < 20");
395     cuts1->cd(2);
396     mytree->Draw("jets.pt3","jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 30 && misc.pthiggs > 200 && leptons.pdgid != 0");
397     cuts1->cd(3);
398     mytree->Draw("jets.pt3","jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 40 && misc.pthiggs > 200 && leptons.pdgid != 0");
399     cuts1->cd(4);
400     mytree->Draw("jets.pt3","(leptons.muon1_pt > 25 || leptons.electron1_pt > 25) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 20 && misc.pthiggs > 200 && leptons.pdgid != 0");
401     cuts1->cd(5);
402     mytree->Draw("jets.pt3","(leptons.muon1_pt > 25 || leptons.electron1_pt > 25) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 30 && misc.pthiggs > 200 && leptons.pdgid != 0");
403     cuts1->cd(6);
404     mytree->Draw("jets.pt3","(leptons.muon1_pt > 25 || leptons.electron1_pt > 25) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 40 && misc.pthiggs > 200 && leptons.pdgid != 0");
405     cuts2->cd(1);
406     mytree->Draw("jets.pt3","(leptons.muon1_pt > 30 || leptons.electron1_pt > 30) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 20 && misc.pthiggs > 200 && leptons.pdgid != 0");
407     cuts2->cd(2);
408     smytree->Draw("jets.pt3","(leptons.muon1_pt > 30 || leptons.electron1_pt > 30) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 30 && misc.pthiggs > 200 && leptons.pdgid != 0");
409     cuts2->cd(3);
410     mytree->Draw("jets.pt3","(leptons.muon1_pt > 30 || leptons.electron1_pt > 30) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 40 && misc.pthiggs > 200 && leptons.pdgid != 0");
411     cuts2->cd(4);
412     mytree->Draw("jets.pt3","(leptons.muon1_pt > 35 || leptons.electron1_pt > 35) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 20 && misc.pthiggs > 200 && leptons.pdgid != 0");
413     cuts2->cd(5);
414     mytree->Draw("jets.pt3","(leptons.muon1_pt > 35 || leptons.electron1_pt > 35) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 30 && misc.pthiggs > 200 && leptons.pdgid != 0");
415     cuts2->cd(6);
416     mytree->Draw("jets.pt3","(leptons.muon1_pt > 35 || leptons.electron1_pt > 35) && jets.btag1 > 0.1 && jets.btag2 > 0.1 && jets.pt3 < 40 && misc.pthiggs > 200 && leptons.pdgid != 0");
417    
418     cuts1->cd(1);
419     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 25 && (leptons.muon2_pt > 10 || leptons.electron1_pt > 10)) || (leptons.electron1_pt > 25 && (leptons.electron2_pt > 10 || leptons.muon1_pt > 10))) && misc.pthiggs > 200 && leptons.pdgid != 0");
420     cuts1->cd(2);
421     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 25 && (leptons.muon2_pt > 15 || leptons.electron1_pt > 15)) || (leptons.electron1_pt > 25 && (leptons.electron2_pt > 15 || leptons.muon1_pt > 15))) && misc.pthiggs > 200 && leptons.pdgid != 0");
422     cuts1->cd(3);
423     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 25 && (leptons.muon2_pt > 20 || leptons.electron1_pt > 20)) || (leptons.electron1_pt > 25 && (leptons.electron2_pt > 20 || leptons.muon1_pt > 20))) && misc.pthiggs > 200 && leptons.pdgid != 0");
424     cuts1->cd(4);
425     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 30 && (leptons.muon2_pt > 10 || leptons.electron1_pt > 10)) || (leptons.electron1_pt > 30 && (leptons.electron2_pt > 10 || leptons.muon1_pt > 10))) && misc.pthiggs > 200 && leptons.pdgid != 0");
426     cuts1->cd(5);
427     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 30 && (leptons.muon2_pt > 15 || leptons.electron1_pt > 15)) || (leptons.electron1_pt > 30 && (leptons.electron2_pt > 15 || leptons.muon1_pt > 15))) && misc.pthiggs > 200 && leptons.pdgid != 0");
428     cuts1->cd(6);
429     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 30 && (leptons.muon2_pt > 20 || leptons.electron1_pt > 20)) || (leptons.electron1_pt > 30 && (leptons.electron2_pt > 20 || leptons.muon1_pt > 20))) && misc.pthiggs > 200 && leptons.pdgid != 0");
430     cuts2->cd(1);
431     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 30 && (leptons.muon2_pt > 10 || leptons.electron1_pt > 10)) || (leptons.electron1_pt > 35 && (leptons.electron2_pt > 10 || leptons.muon1_pt > 10))) && misc.pthiggs > 200 && leptons.pdgid != 0");
432     cuts2->cd(2);
433     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 35 && (leptons.muon2_pt > 15 || leptons.electron1_pt > 15)) || (leptons.electron1_pt > 35 && (leptons.electron2_pt > 15 || leptons.muon1_pt > 15))) && misc.pthiggs > 200 && leptons.pdgid != 0");
434     cuts2->cd(3);
435     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 35 && (leptons.muon2_pt > 20 || leptons.electron1_pt > 20)) || (leptons.electron1_pt > 35 && (leptons.electron2_pt > 20 || leptons.muon1_pt > 20))) && misc.pthiggs > 200 && leptons.pdgid != 0");
436     cuts2->cd(4);
437     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 35 && (leptons.muon2_pt > 10 || leptons.electron1_pt > 10)) || (leptons.electron1_pt > 20 && (leptons.electron2_pt > 10 || leptons.muon1_pt > 10))) && misc.pthiggs > 200 && leptons.pdgid != 0");
438     cuts2->cd(5);
439     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 20 && (leptons.muon2_pt > 15 || leptons.electron1_pt > 15)) || (leptons.electron1_pt > 20 && (leptons.electron2_pt > 15 || leptons.muon1_pt > 15))) && misc.pthiggs > 200 && leptons.pdgid != 0");
440     cuts2->cd(6);
441     mytree->Draw("leptons.muon1_pt","((leptons.muon1_pt > 20 && (leptons.muon2_pt > 20 || leptons.electron1_pt > 20)) || (leptons.electron1_pt > 20 && (leptons.electron2_pt > 20 || leptons.muon1_pt > 20))) && misc.pthiggs > 200 && leptons.pdgid != 0");
442    
443     cuts2->cd(1);
444     mytree->Draw("leptons.muon1_pt","leptons.muon1_pt > 35 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.muon1_iso_03 < 5");
445     cuts2->cd(2);
446     mytree->Draw("leptons.muon1_pt","leptons.muon1_pt > 35 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.muon1_iso_03 < 10");
447     cuts2->cd(3);
448     mytree->Draw("leptons.muon1_pt","leptons.muon1_pt > 35 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.muon1_iso_03 < 15");
449     cuts2->cd(4);
450     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 25 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 5");
451     cuts2->cd(5);
452     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 25 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 10");
453     cuts2->cd(6);
454     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 25 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 15");
455    
456     cuts3->cd(1);
457     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 30 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 5");
458     cuts3->cd(2);
459     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 30 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 10");
460     cuts3->cd(3);
461     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 30 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 15");
462     cuts3->cd(4);
463     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 35 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 5");
464     cuts3->cd(5);
465     mytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 35 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 10");
466     cuts3->cd(6);
467     imytree->Draw("leptons.electron1_pt","leptons.electron1_pt > 35 && misc.pthiggs > 200 && leptons.pdgid != 0 && leptons.electron1_iso_03 < 15");
468     */
469     //cuts->cd();
470     //mytree->Draw("leptons.muon1_pt + leptons.electron1_pt","leptons.etagen < 2.4" && "leptons.etagen > 2.4");
471    
472     savefile->Write();
473     Hlist.Write();
474     mytree->Write();
475     //if(mytree->GetEntryList() != 0)
476     // (mytree->GetEntryList())->Write();
477    
478     outfile.close();
479    
480     }
481    
482     void Leptons::Init_Leptons()
483     {
484     pdgid = recomuons = recoelectrons = ptgen = etagen = muon1_pt = muon1_eta = muon1_phi = muon1_ep = muon1_ep_mod = muon1_hovere = muon1_bjetangle =
485     muon1_nearestjetangle = muon1_iso_03 = muon1_iso_05_emEt = muon1_ipt = muon1_ipt_sigma = muon1_ip3D = muon1_ip3D_sigma =
486     electron1_pt = electron1_eta = electron1_phi = electron1_deltaeta = electron1_deltaphi = electron1_sigmaetaeta = electron1_ep = electron1_ep_mod =
487     electron1_hovere = electron1_bjetangle = electron1_nearestjetangle = electron1_iso_03 = electron1_iso_04_emEt =
488     electron1_ipt = electron1_ipt_sigma = electron1_ip3D = electron1_ip3D_sigma = -10;
489     addlepton = muon1_global = 0;
490     }
491    
492     void Met::Init_Met()
493     {
494     tc_met = calo_met = et_neutrinos = tc_met_lepton = calo_met_lepton = -10;
495     }
496    
497     void Pfjets::Init_Pfjets()
498     {
499     pt1 = pt2 = pt3 = phi1 = phi2 = phi3 = eta1 = eta2 = eta3 = btag1 = btag2 = btag3 = hardb_pt = hardbjet_pt = hardbjet_deltar =
500     softb_pt = softbjet_pt = softbjet_deltar = deltar = nearestjet_hcal = nearestjet_ntracks = -10;
501     }
502    
503     void Misc::Init_Misc()
504     {
505     pthiggs = met_higgs_angle = calomet_higgs_angle = higgs_pt = higgs_invmass = higgs_eta = higgs_phi = w_pt = w_invmass = w_eta = w_phi = -10;
506     }
507    
508     void leptons_pt_iso_cut(const Double_t cut1, const Double_t cut2, const Double_t cut3, TTree *maintree)
509     {
510     char string1[22];
511     char string2[22];
512     char string3[22];
513     sprintf(string1,"%2.1f",cut1);
514     sprintf(string2,"%2.1f",cut2);
515     sprintf(string3,"%2.1f",cut3);
516     cout << string1 << " -- " << string2 << " -- " << string3 << endl;
517     string cut = "leptons.addlepton==0 && (leptons.electron1_pt>" + string(string1) + " && leptons.electron1_iso_03<" + string(string3) + " || leptons.muon1_pt>" + string(string1) + " && leptons.muon1_iso_03<" + string(string2) + ")";
518     if(maintree->GetEntryList() != 0 )
519     nentries = (maintree->GetEntryList())->GetN();
520     else
521     nentries = maintree->GetEntries();
522     cout << "nentries = " << nentries << endl;
523    
524     maintree->Draw(">>leptons_pt_iso_list",cut.c_str(),"entrylist");
525     TEntryList *leptons_pt_iso_list = (TEntryList*)gDirectory->Get("leptons_pt_iso_list");
526     leptons_pt_iso_list->SetReapplyCut(kTRUE);
527     maintree->SetEntryList(leptons_pt_iso_list);
528    
529     selentries = (maintree->GetEntryList())->GetN();
530     if(nentries != 0 && totalevents != 0)
531     outfile << "Lept(no add. leptons). pt > " << string1 << " iso < " << string2 << "/" << string3 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
532     else
533     outfile << "Lept. pt > " << string1 << " iso < " << string2 << "/" << string3 << "either before or after cut, #events = 0" << endl;
534     outfile.flush();
535     }
536    
537     void leptons_isgood_cut(const Double_t cut1, const Double_t cut2, const Double_t cut3, const Double_t cut4, const Double_t cut5, TTree *maintree)
538     {
539     char string1[22];
540     char string2[22];
541     char string3[22];
542     char string4[22];
543     char string5[22];
544    
545     sprintf(string1,"%2.3f",cut1);
546     sprintf(string2,"%2.3f",cut2);
547     sprintf(string3,"%2.3f",cut3);
548     sprintf(string4,"%2.3f",cut4);
549     sprintf(string5,"%2.3f",cut5);
550    
551     cout << string1 << " -- " << string2 << " -- " << string3 << " -- " << string4 << " -- " << string5 << endl;
552     string cut = "leptons.addlepton == 0 && (leptons.muon1_pt > 25 && leptons.muon1_global == 1 || leptons.electron1_pt > 25 && leptons.electron1_deltaeta<" + string(string1) + " && leptons.electron1_deltaphi<" + string(string2) + " && leptons.electron1_sigmaetaeta<" + string(string3) + " && leptons.electron1_ep<" + string(string4) + " && leptons.electron1_hovere<" + string(string5) + ")";
553     if(maintree->GetEntryList() != 0 )
554     nentries = (maintree->GetEntryList())->GetN();
555     else
556     nentries = maintree->GetEntries();
557     cout << "nentries = " << nentries << endl;
558    
559     maintree->Draw(">>leptons_isgood_list",cut.c_str(),"entrylist");
560     TEntryList *leptons_isgood_list = (TEntryList*)gDirectory->Get("leptons_isgood_list");
561     leptons_isgood_list->SetReapplyCut(kTRUE);
562     maintree->SetEntryList(leptons_isgood_list);
563    
564     selentries = (maintree->GetEntryList())->GetN();
565     if(nentries != 0 && totalevents != 0)
566     outfile << "Fake leptons. delta_eta < " << string1 << " delta_phi < " << string2 << " sigma_Ieta_Ieta < " << string3 << " E/p < " << string4 << " H/E > " << string5 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
567     else
568     outfile << "Fake leptons. delta_eta < " << string1 << " delta_phi < " << string2 << " sigma_Ieta_Ieta < " << string3 << " E/p < " << string4 << " H/E > " << string5 << " either before or after cut, #events = 0" << endl;
569     outfile.flush();
570     }
571    
572     void leptons_bjetangle_emEt_cut(const Double_t cut1, const Double_t cut2, const Double_t cut3, TTree *maintree)
573     {
574     char string1[22];
575     char string2[22];
576     char string3[22];
577    
578     sprintf(string1,"%2.3f",cut1);
579     sprintf(string2,"%2.3f",cut2);
580     sprintf(string3,"%2.3f",cut3);
581    
582     cout << string1 << " -- " << string2 << " -- " << string3 << endl;
583     string cut = "(leptons.muon1_pt>25 && leptons.muon1_bjetangle>" + string(string1) + " && leptons.muon1_bjetangle<" + string(string2) +
584     " && leptons.muon1_iso_05_emEt<" + string(string3) +
585     ") || (leptons.electron1_pt>25 && leptons.electron1_bjetangle>" + string(string1) + " && leptons.electron1_bjetangle<" +
586     string(string2) + ")";
587     if(maintree->GetEntryList() != 0 )
588     nentries = (maintree->GetEntryList())->GetN();
589     else
590     nentries = maintree->GetEntries();
591     cout << "nentries = " << nentries << endl;
592    
593     maintree->Draw(">>leptons_bjetangle_emEt_list",cut.c_str(),"entrylist");
594     TEntryList *leptons_bjetangle_emEt_list = (TEntryList*)gDirectory->Get("leptons_bjetangle_emEt_list");
595     leptons_bjetangle_emEt_list->SetReapplyCut(kTRUE);
596     maintree->SetEntryList(leptons_bjetangle_emEt_list);
597    
598     selentries = (maintree->GetEntryList())->GetN();
599     if(nentries != 0 && totalevents != 0)
600     outfile << "Fake leptons. angle(lepton,bjet) > " << string1 << " angle(lepton,bjet) < " << string2 << " muon_emEt < " << string3 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
601     else
602     outfile << "Fake leptons. angle(lepton,bjet) > " << string1 << " angle(lepton,bjet) < " << string2 << " muon_emEt < " << string3 << " either before or after cut, #events = 0" << endl;
603     outfile.flush();
604     }
605    
606     void leptons_jets_cut(const Double_t cut1, const Double_t cut2, const Double_t cut3, TTree *maintree)
607     {
608     char string1[22];
609     char string2[22];
610     char string3[22];
611    
612     sprintf(string1,"%2.3f",cut1);
613     sprintf(string2,"%2.3f",cut2);
614     sprintf(string3,"%2.3f",cut3);
615    
616     cout << string1 << " -- " << string2 << " -- " << string3 << endl;
617     string cut = "leptons.electron1_nearestjetangle>" + string(string3) + " && leptons.muon1_nearestjetangle>" + string(string3) + " || pfjets.nearestjet_ntracks<" + string(string1) + " && pfjets.nearestjet_hcal<" + string(string2);
618     if(maintree->GetEntryList() != 0 )
619     nentries = (maintree->GetEntryList())->GetN();
620     else
621     nentries = maintree->GetEntries();
622     cout << "nentries = " << nentries << endl;
623    
624     maintree->Draw(">>leptons_jets_list",cut.c_str(),"entrylist");
625     TEntryList *leptons_jets_list = (TEntryList*)gDirectory->Get("leptons_jets_list");
626     leptons_jets_list->SetReapplyCut(kTRUE);
627     maintree->SetEntryList(leptons_jets_list);
628    
629     selentries = (maintree->GetEntryList())->GetN();
630     if(nentries != 0 && totalevents != 0)
631     outfile << "Nearest jets (deltaR < " << string3 << "). Hcal < " << string1 << " ntracks < " << string2 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
632     else
633     outfile << "Nearest jets (deltaR < " << string3 << "). Hcal < " << string1 << " ntracks < " << string2 << " either before or after cut, #events = 0" << endl;
634     outfile.flush();
635     }
636    
637     void pfjets_btag_cut(const Double_t cut1, TTree *maintree)
638     {
639     char string1[22];
640     sprintf(string1,"%2.1f",cut1);
641     cout << string1 << endl;
642    
643     if(maintree->GetEntryList() != 0 )
644     nentries = (maintree->GetEntryList())->GetN();
645     else
646     nentries = maintree->GetEntries();
647     cout << "nentries = " << nentries << endl;
648    
649     string cut = "pfjets.btag1>" + string(string1) + " && pfjets.btag2>" + string(string1);
650     maintree->Draw(">>pfjets_btag_list",cut.c_str(),"entrylist");
651     TEntryList *pfjets_btag_list = (TEntryList*)gDirectory->Get("pfjets_btag_list");
652     pfjets_btag_list->SetReapplyCut(kTRUE);
653     maintree->SetEntryList(pfjets_btag_list);
654    
655     selentries = (maintree->GetEntryList())->GetN();
656     if(nentries != 0 && totalevents != 0)
657     outfile << "Pfjets. btag > " << string1 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
658     else
659     outfile << "Pfjets. btag > " << string1 << " either before or after cut, #events = 0" << endl;
660     outfile.flush();
661     }
662    
663     void pfjets_pt_cut(const Double_t cut1, const Double_t cut2, TTree *maintree)
664     {
665     char string1[22];
666     sprintf(string1,"%2.1f",cut1);
667     char string2[22];
668     sprintf(string2,"%2.1f",cut2);
669     cout << string1 << " -- " << string2 << endl;
670    
671     if(maintree->GetEntryList() != 0 )
672     nentries = (maintree->GetEntryList())->GetN();
673     else
674     nentries = maintree->GetEntries();
675     cout << "nentries = " << nentries << endl;
676    
677     string cut = "pfjets.pt1>" + string(string1) + " && pfjets.pt2>" + string(string2);
678     maintree->Draw(">>pfjets_pt_list",cut.c_str(),"entrylist");
679     TEntryList *pfjets_pt_list = (TEntryList*)gDirectory->Get("pfjets_pt_list");
680     pfjets_pt_list->SetReapplyCut(kTRUE);
681     maintree->SetEntryList(pfjets_pt_list);
682    
683     selentries = (maintree->GetEntryList())->GetN();
684     if(nentries != 0 && totalevents != 0)
685     outfile << "Pfjets. pt1 > " << string1 << " pt2 > " << string2 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
686     else
687     outfile << "Pfjets. pt1 > " << string1 << " pt2 > " << string2 << " either before or after cut, #events = 0" << endl;
688     outfile.flush();
689     }
690    
691     void w_higgs_angle_cut(const Double_t cut1, TTree *maintree)
692     {
693     char string1[22];
694     sprintf(string1,"%2.1f",cut1);
695     cout << string1 << endl;
696    
697     if(maintree->GetEntryList() != 0 )
698     nentries = (maintree->GetEntryList())->GetN();
699     else
700     nentries = maintree->GetEntries();
701     cout << "nentries = " << nentries << endl;
702    
703     string cut = "misc.met_higgs_angle<" + string(string1);
704     maintree->Draw(">>w_higgs_angle_list",cut.c_str(),"entrylist");
705     TEntryList *w_higgs_angle_list = (TEntryList*)gDirectory->Get("w_higgs_angle_list");
706     w_higgs_angle_list->SetReapplyCut(kTRUE);
707     maintree->SetEntryList(w_higgs_angle_list);
708    
709     selentries = (maintree->GetEntryList())->GetN();
710     if(nentries != 0 && totalevents != 0)
711     outfile << "WH angle. angle < " << string1 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
712     else
713     outfile << "WH angle. angle < " << string1 << " either before or after cut, #events = 0" << endl;
714     outfile.flush();
715     }
716    
717     void higgs_invmass_cut(const Double_t cut1, TTree *maintree)
718     {
719     char string1[22];
720     sprintf(string1,"%2.1f",cut1);
721     cout << string1 << endl;
722    
723     if(maintree->GetEntryList() != 0 )
724     nentries = (maintree->GetEntryList())->GetN();
725     else
726     nentries = maintree->GetEntries();
727     cout << "nentries = " << nentries << endl;
728    
729     string cut = "higgs_invmass<" + string(string1);
730     maintree->Draw(">>higgs_invmass_list",cut.c_str(),"entrylist");
731     TEntryList *higgs_invmass_list = (TEntryList*)gDirectory->Get("higgs_invmass_list");
732     higgs_invmass_list->SetReapplyCut(kTRUE);
733     maintree->SetEntryList(higgs_invmass_list);
734    
735     selentries = (maintree->GetEntryList())->GetN();
736    
737     outfile << "Higgs inv. mass. inv. mass < " << string1;
738     if(nentries != 0 || totalevents != 0)
739     outfile << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
740     else
741     outfile << " either before or after cut, #events = 0" << endl;
742     outfile.flush();
743     }
744    
745     void higgs_pt_cut(const Double_t cut1, TTree *maintree)
746     {
747     char string1[22];
748     sprintf(string1,"%2.1f",cut1);
749     cout << string1 << endl;
750    
751     if(maintree->GetEntryList() != 0 )
752     nentries = (maintree->GetEntryList())->GetN();
753     else
754     nentries = maintree->GetEntries();
755     cout << "nentries = " << nentries << endl;
756    
757     string cut = "higgs_pt>" + string(string1);
758     maintree->Draw(">>higgs_pt_list",cut.c_str(),"entrylist");
759     TEntryList *higgs_pt_list = (TEntryList*)gDirectory->Get("higgs_pt_list");
760     higgs_pt_list->SetReapplyCut(kTRUE);
761     maintree->SetEntryList(higgs_pt_list);
762    
763     selentries = (maintree->GetEntryList())->GetN();
764    
765     outfile << "Higgs pt. pt > " << string1;
766     if(nentries != 0 || totalevents != 0)
767     outfile << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
768     else
769     outfile << " either before or after cut, #events = 0" << endl;
770     outfile.flush();
771     }
772    
773     void pfjets_angle_cut(const Double_t cut1, TTree *maintree)
774     {
775     char string1[22];
776     sprintf(string1,"%2.1f",cut1);
777     cout << string1 << endl;
778    
779     if(maintree->GetEntryList() != 0 )
780     nentries = (maintree->GetEntryList())->GetN();
781     else
782     nentries = maintree->GetEntries();
783     cout << "nentries = " << nentries << endl;
784    
785     string cut = "pfjets.deltar<" + string(string1);
786     maintree->Draw(">>pfjets_angle_list",cut.c_str(),"entrylist");
787     TEntryList *pfjets_angle_list = (TEntryList*)gDirectory->Get("pfjets_angle_list");
788     pfjets_angle_list->SetReapplyCut(kTRUE);
789     maintree->SetEntryList(pfjets_angle_list);
790    
791     selentries = (maintree->GetEntryList())->GetN();
792     if(nentries != 0 && totalevents != 0)
793     outfile << "Pfjets. angle < " << string1 << " - Rel eff: " << selentries/nentries << " - Abs eff: " << selentries/totalevents << endl;
794     else
795     outfile << "Pfjets. angle < " << string1 << " either before or after cut, #events = 0" << endl;
796     outfile.flush();
797     }
798    
799    
800     Double_t deltaR(Double_t eta1, Double_t phi1, Double_t eta2, Double_t phi2)
801     {
802     Double_t deta = eta1 - eta2;
803     Double_t dphi = phi1 - phi2;
804     while (dphi > TMath::Pi()) dphi -= 2*TMath::Pi();
805     while (dphi <= -TMath::Pi()) dphi += 2*TMath::Pi();
806     return sqrt(deta*deta + dphi*dphi);
807     }
808