ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/TreeReader_New.C
Revision: 1.1
Committed: Fri Feb 13 00:53:03 2009 UTC (16 years, 2 months ago) by kukartse
Content type: text/plain
Branch: MAIN
CVS Tags: V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15, V00-01-14, V00-01-13, V00-01-12, V00-01-11, V00-01-10, gak031009, gak030509, gak022309, gak021209, HEAD
Branch point for: ZMorph-V00-03-01, CMSSW_22X_branch
Log Message: 
*** empty log message ***

File Contents

# Content
1 #include <fstream>
2 #include <iomanip>
3 #include <iostream>
4 #include <string.h>
5 #include <fstream>
6
7 #include <TFile.h>
8 #include <TH1.h>
9 #include <TTree.h>
10 #include <TCanvas.h>
11 #include <TF1.h>
12 #include <TStyle.h>
13 #include <TCut.h>
14 #include "cuts.C"
15
16 int TreeReader(char *infile = "/uscms/home/kukarzev/work/disk_ljmet/kukarzev/tmp/ljmet_tmva_InclusiveMuPt15-12feb2009.root",char *outfile = "test.root")
17 {
18 TFile *fin = TFile::Open(infile);
19 TFile *outf = new TFile(outfile,"RECREATE");
20
21
22
23 TH1D * hMt_nominal = new TH1D( "hMt_nominal", "Mt (nominal selection)", 100, 0., 90.);
24 TH1D * hMt_reversed = new TH1D( "hMt_reversed", "Mt (reversed selection)", 100, 0., 90.);
25
26
27
28 TH1D * hht_nominal = new TH1D( "hht_nominal", "ht (nominal) ", 500, 0., 2000.);
29 TH1D * hht_reversed = new TH1D( "hht_reversed", "ht (reversed) ", 500, 0., 2000.);
30
31
32 TH1D * hM123inv_nominal = new TH1D( "hM123inv_nominal", "M123inv (nominal selection)", 500, 0., 4500.);
33 TH1D * hM123inv_reversed = new TH1D( "hM123inv_reversed", "M123inv (reversed selection)", 500, 0., 4500.);
34
35
36 TH1D * hgetHt3_nominal = new TH1D( "hgetHt3_nominal", "Ht3 (nominal selection)", 100, 0., 900.);
37 TH1D * hgetHt3_reversed = new TH1D( "hgetHt3_reversed", "Ht3 (reversed selection)", 100, 0., 900.);
38
39
40
41 TH1D * hktMinPrime_nominal = new TH1D( "hktMinPrime_nominal", "ktMinPrime (nominal selection)", 100, 0., 25.);
42 TH1D * hktMinPrime_reversed = new TH1D( "hktMinPrime_reversed", "ktMinPrime (reversed selection)", 100, 0., 25.);
43
44
45
46 TH1D * hht2p_nominal = new TH1D( "hht2p_nominal", "ht2p (nominal selection)", 100, 0., 4.);
47 TH1D * hht2p_reversed = new TH1D( "hht2p_reversed", "ht2p (reversed selection)", 100, 0., 4.);
48
49
50 TH1D * hcentrality_nominal = new TH1D( "hcentrality_nominal", "centrality (nominal selection)", 100, 0., 1.);
51 TH1D * hcentrality_reversed = new TH1D( "hcentrality_reversed", "centrality (reversed selection)", 100, 0., 1.);
52
53
54 TH1D * hsphericity_nominal = new TH1D( "hsphericity_nominal", "sphericity (nominal selection)", 100, 0., 1.);
55 TH1D * hsphericity_reversed = new TH1D( "hsphericity_reversed", "sphericity (reversed selection)", 100, 0., 1.);
56
57
58 TH1D * hhtPlusLepton_nominal = new TH1D( "hhtPlusLepton_nominal", "htPlusLepton (nominal selection)", 500, 0., 3000.);
59 TH1D * hhtPlusLepton_reversed = new TH1D( "hhtPlusLepton_reversed", "htPlusLepton (reversed selection)", 500, 0., 3000.);
60
61
62 TH1D * hdPhiLMet_nominal = new TH1D( "hdPhiLMet_nominal", "dPhiLMet (nominal selection)", 100, 0., 3.5);
63 TH1D * hdPhiLMet_reversed = new TH1D( "hdPhiLMet_reversed", "dPhiLMet (reversed selection)", 100, 0., 3.5);
64
65
66 TH1D * hDphiJMET_nominal = new TH1D( "hDphiJMET_nominal", "DphiJMET (nominal selection)", 100, 0., 3.5);
67 TH1D * hDphiJMET_reversed = new TH1D( "hDphiJMET_reversed", "DphiJMET (reversed selection)", 100, 0., 3.5);
68
69
70
71 TH1D * haplanarity_nominal = new TH1D( "haplanarity_nominal", "aplanarity (nominal selection)", 100, 0., 0.5);
72 TH1D * haplanarity_reversed = new TH1D( "haplanarity_reversed", "aplanarity (reversed selection)", 100, 0., 0.5);
73
74
75 TH1D * hW_MT_nominal = new TH1D( "hW_MT_nominal", "W_MT (nominal selection)", 500, 0., 5000);
76 TH1D * hW_MT_reversed = new TH1D( "hW_MT_reversed", "W_MT (reversed selection)", 500, 0., 5000);
77
78
79 TH1D * hmetHtPlusLepton_nominal = new TH1D( "hmetHtPlusLepton_nominal", "metHtPlusLepton (nominal selection)", 500, 0., 4000.);
80 TH1D * hmetHtPlusLepton_reversed = new TH1D( "hmetHtPlusLepton_reversed", "metHtPlusLepton (reversed selection)", 500, 0., 4000.);
81
82
83
84
85 TTree * t1 = (TTree*)fin->Get("ttljets");
86
87 t1->Draw("ht>> hht_nominal",semilep_cut && reversible_cut);
88 t1->Draw("getMt>> hMt_nominal",semilep_cut && reversible_cut);
89 t1->Draw("getM123inv>> hM123inv_nominal",semilep_cut && reversible_cut);
90 t1->Draw("getHt3>> hgetHt3_nominal",semilep_cut && reversible_cut);
91 t1->Draw("ktMinPrime>> hktMinPrime_nominal",semilep_cut && reversible_cut);
92 t1->Draw("ht2p>> hht2p_nominal",semilep_cut && reversible_cut);
93 t1->Draw("aplanarity>> haplanarity_nominal",semilep_cut && reversible_cut);
94 t1->Draw("centrality>> hcentrality_nominal",semilep_cut && reversible_cut);
95 t1->Draw("sphericity>> hsphericity_nominal",semilep_cut && reversible_cut);
96 t1->Draw("dPhiLMet>> hdPhiLMet_nominal",semilep_cut && reversible_cut);
97 t1->Draw("DphiJMET>> hDphiJMET_nominal",semilep_cut && reversible_cut);
98 t1->Draw("W_MT>> hW_MT_nominal",semilep_cut && reversible_cut);
99 t1->Draw("htPlusLepton>> hhtPlusLepton_nominal",semilep_cut && reversible_cut);
100 t1->Draw("metHtPlusLepton>> hmetHtPlusLepton_nominal",semilep_cut && reversible_cut);
101
102
103
104 t1->Draw("ht>> hht_reversed",semilep_cut && !reversible_cut);
105 t1->Draw("getMt>> hMt_reversed",semilep_cut && !reversible_cut);
106 t1->Draw("getM123inv>> hM123inv_reversed",semilep_cut && !reversible_cut);
107 t1->Draw("getHt3>> hgetHt3_reversed",semilep_cut && !reversible_cut);
108 t1->Draw("ktMinPrime>> hktMinPrime_reversed",semilep_cut && !reversible_cut);
109 t1->Draw("ht2p>> hht2p_reversed",semilep_cut && !reversible_cut);
110 t1->Draw("aplanarity>> haplanarity_reversed",semilep_cut && !reversible_cut);
111 t1->Draw("centrality>> hcentrality_reversed",semilep_cut && !reversible_cut);
112 t1->Draw("sphericity>> hsphericity_reversed",semilep_cut && !reversible_cut);
113 t1->Draw("dPhiLMet>> hdPhiLMet_reversed",semilep_cut && !reversible_cut);
114 t1->Draw("DphiJMET>> hDphiJMET_reversed",semilep_cut && !reversible_cut);
115 t1->Draw("W_MT>> hW_MT_reversed",semilep_cut && !reversible_cut);
116 t1->Draw("htPlusLepton>> hhtPlusLepton_reversed",semilep_cut && !reversible_cut);
117 t1->Draw("metHtPlusLepton>> hmetHtPlusLepton_reversed",semilep_cut && !reversible_cut);
118
119
120 outf->Write();
121 outf->Close();
122
123
124 }