forked from JosephineWittkowski/analysis_SUSYTools
-
Notifications
You must be signed in to change notification settings - Fork 0
/
cutflow_data.C
164 lines (140 loc) · 15.3 KB
/
cutflow_data.C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
{
// test
TFile hn("/data/etp/jwittkowski/analysis_SUSYTools_03_04/histos_cutflow_data.root");
hn.Print();
TH1F *alpha_raw = (TH1F*)hn.Get("cutflow_EE");
TH1F *beta_raw = (TH1F*)hn.Get("cutflow_EM");
TH1F *gamma_raw = (TH1F*)hn.Get("cutflow_MM");
cout << "raw events:" << endl;
cout << fixed;
cout << "events :" << beta_raw->GetBinContent( 1) << endl;
cout << "cut grl :" << beta_raw->GetBinContent( 2) << endl;
cout << "TileTripReader :" << beta_raw->GetBinContent( 3) << endl;
cout << "incomplete veto :" << beta_raw->GetBinContent( 4) << endl;
cout << "cut LArTileError :" << beta_raw->GetBinContent( 5) << endl;
cout << "cut TileCalHotSpot :" << beta_raw->GetBinContent( 6) << endl;
cout << "cut BadJets :" << beta_raw->GetBinContent( 7) << endl;
cout << "cut CaloJets :" << beta_raw->GetBinContent( 8) << endl;
cout << "cut PrimaryVertex :" << beta_raw->GetBinContent( 9) << endl;
cout << "cut BadMuons :" << beta_raw->GetBinContent( 10) << endl;
cout << "Cosmic Muons :" << beta_raw->GetBinContent(11) << endl;
cout << "hfor veto :" << beta_raw->GetBinContent(12) << endl;
cout << ">= 2 lep :" << beta_raw->GetBinContent(13) << endl;
cout << "== 2 lept :" << beta_raw->GetBinContent(14) << endl;
cout << "mll cut :" << beta_raw->GetBinContent(15) << endl;
cout << "---------------- SS --------------------------" << endl;
cout << " EE MM EM " << endl;
cout << "----------------------------------------------" << endl;
cout << "pass category :" << alpha_raw->GetBinContent(16) << " " << gamma_raw->GetBinContent(16) << " " << beta_raw->GetBinContent(16) << endl;
cout << "flavor cut :" << alpha_raw->GetBinContent(17) << " " << gamma_raw->GetBinContent(17) << " " << beta_raw->GetBinContent(17) << endl;
cout << "tau veto :" << alpha_raw->GetBinContent(18) << " " << gamma_raw->GetBinContent(18) << " " << beta_raw->GetBinContent(18) << endl;
cout << "valid pT :" << alpha_raw->GetBinContent(19) << " " << gamma_raw->GetBinContent(19) << " " << beta_raw->GetBinContent(19) << endl;
cout << "Trigger :" << alpha_raw->GetBinContent(20) << " " << gamma_raw->GetBinContent(20) << " " << beta_raw->GetBinContent(20) << endl;
cout << "RealLeptonsCheck :" << alpha_raw->GetBinContent(21) << " " << gamma_raw->GetBinContent(21) << " " << beta_raw->GetBinContent(21) << endl;
cout << "SameSign :" << alpha_raw->GetBinContent(22) << " " << gamma_raw->GetBinContent(22) << " " << beta_raw->GetBinContent(22) << endl;
cout << "Iso :" << alpha_raw->GetBinContent(23) << " " << gamma_raw->GetBinContent(23) << " " << beta_raw->GetBinContent(23) << endl;
cout << "dOS :" << alpha_raw->GetBinContent(24) << " " << gamma_raw->GetBinContent(24) << " " << beta_raw->GetBinContent(24) << endl;
cout << "FVeto :" << alpha_raw->GetBinContent(25) << " " << gamma_raw->GetBinContent(25) << " " << beta_raw->GetBinContent(25) << endl;
cout << "bVeto :" << alpha_raw->GetBinContent(26) << " " << gamma_raw->GetBinContent(26) << " " << beta_raw->GetBinContent(26) << endl;
cout << ">=1jet :" << alpha_raw->GetBinContent(27) << " " << gamma_raw->GetBinContent(27) << " " << beta_raw->GetBinContent(27) << endl;
cout << "lep pt :" << alpha_raw->GetBinContent(28) << " " << gamma_raw->GetBinContent(28) << " " << beta_raw->GetBinContent(28) << endl;
cout << "mll :" << alpha_raw->GetBinContent(29) << " " << gamma_raw->GetBinContent(29) << " " << beta_raw->GetBinContent(29) << endl;
cout << "mWWT :" << alpha_raw->GetBinContent(30) << " " << gamma_raw->GetBinContent(30) << " " << beta_raw->GetBinContent(30) << endl;
cout << "HT :" << alpha_raw->GetBinContent(31) << " " << gamma_raw->GetBinContent(31) << " " << beta_raw->GetBinContent(31) << endl;
cout << "metREL :" << alpha_raw->GetBinContent(32) << " " << gamma_raw->GetBinContent(32) << " " << beta_raw->GetBinContent(32) << endl;
cout << "MM SRSS2: " << endl;
cout << "mWWT :" << gamma_raw->GetBinContent(33) << endl;
cout << "HT :" << gamma_raw->GetBinContent(34) << endl;
cout << "metREL :" << gamma_raw->GetBinContent(35) << endl;
cout << "MM SRSS3: " << endl;
cout << "mWWT :" << gamma_raw->GetBinContent(36) << endl;
cout << "HT :" << gamma_raw->GetBinContent(37) << endl;
cout << "metREL :" << gamma_raw->GetBinContent(38) << endl;
cout << "MM SRSS4: " << endl;
cout << "mWWT :" << gamma_raw->GetBinContent(39) << endl;
cout << "HT :" << gamma_raw->GetBinContent(40) << endl;
cout << "metREL :" << gamma_raw->GetBinContent(41) << endl;
cout << "EE SRSS:" << endl;
cout << "mt2>90 :" << alpha_raw->GetBinContent(33) << endl;
cout << "---------------- OS --------------------------" << endl;
cout << " EE MM EM " << endl;
cout << "----------------------------------------------" << endl;
cout << "pass OS :" << alpha_raw->GetBinContent(51) << " " << gamma_raw->GetBinContent(51) << " " << beta_raw->GetBinContent(51) << endl;
cout << "pass Iso :" << alpha_raw->GetBinContent(52) << " " << gamma_raw->GetBinContent(52) << " " << beta_raw->GetBinContent(52) << endl;
cout << "pass FVeto :" << alpha_raw->GetBinContent(53) << " " << gamma_raw->GetBinContent(53) << " " << beta_raw->GetBinContent(53) << endl;
cout << "pass bVeto :" << alpha_raw->GetBinContent(54) << " " << gamma_raw->GetBinContent(54) << " " << beta_raw->GetBinContent(54) << endl;
cout << "pass >=2 J :" << alpha_raw->GetBinContent(55) << " " << gamma_raw->GetBinContent(55) << " " << beta_raw->GetBinContent(55) << endl;
cout << "pass Mjj :" << alpha_raw->GetBinContent(56) << " " << gamma_raw->GetBinContent(56) << " " << beta_raw->GetBinContent(56) << endl;
cout << "pass lepPt :" << alpha_raw->GetBinContent(57) << " " << gamma_raw->GetBinContent(57) << " " << beta_raw->GetBinContent(57) << endl;
cout << "pass dRll :" << alpha_raw->GetBinContent(58) << " " << gamma_raw->GetBinContent(58) << " " << beta_raw->GetBinContent(58) << endl;
cout << "pass min(mT1,mT2) :" << alpha_raw->GetBinContent(59) << " " << gamma_raw->GetBinContent(59) << " " << beta_raw->GetBinContent(59) << endl;
cout << "pass dPhi(ll,Met) :" << alpha_raw->GetBinContent(60) << " " << gamma_raw->GetBinContent(60) << " " << beta_raw->GetBinContent(60) << endl;
cout << "pass Met :" << alpha_raw->GetBinContent(61) << " " << gamma_raw->GetBinContent(61) << " " << beta_raw->GetBinContent(61) << endl;
cout << "#############################################" << endl;
cout << "#############################################" << endl;
cout << "MC event Weights: nt.evt()->w only" << endl;
TH1F *alpha_mcEventWeight = (TH1F*)hn.Get("cutflow_EE_ALL");
TH1F *beta_mcEventWeight = (TH1F*)hn.Get("cutflow_EM_ALL");
TH1F *gamma_mcEventWeight = (TH1F*)hn.Get("cutflow_MM_ALL");
//
cout << "events :" << beta_mcEventWeight->GetBinContent( 1) << endl;
cout << "cut grl :" << beta_mcEventWeight->GetBinContent( 2) << endl;
cout << "TileTripReader :" << beta_mcEventWeight->GetBinContent( 3) << endl;
cout << "incomplete veto :" << beta_mcEventWeight->GetBinContent( 4) << endl;
cout << "cut LArTileError :" << beta_mcEventWeight->GetBinContent( 5) << endl;
cout << "cut TileCalHotSpot :" << beta_mcEventWeight->GetBinContent( 6) << endl;
cout << "cut BadJets :" << beta_mcEventWeight->GetBinContent( 7) << endl;
cout << "cut CaloJets :" << beta_mcEventWeight->GetBinContent( 8) << endl;
cout << "cut PrimaryVertex :" << beta_mcEventWeight->GetBinContent( 9) << endl;
cout << "cut BadMuons :" << beta_mcEventWeight->GetBinContent( 10) << endl;
cout << "Cosmic Muons :" << beta_mcEventWeight->GetBinContent(11) << endl;
cout << "hfor veto :" << beta_mcEventWeight->GetBinContent(12) << endl;
cout << ">= 2 lep :" << beta_mcEventWeight->GetBinContent(13) << endl;
cout << "== 2 lept :" << beta_mcEventWeight->GetBinContent(14) << endl;
cout << "mll cut :" << beta_mcEventWeight->GetBinContent(15) << endl;
cout << "---------------- SS --------------------------" << endl;
cout << " EE MM EM " << endl;
cout << "----------------------------------------------" << endl;
cout << "pass category :" << alpha_mcEventWeight->GetBinContent(16) << " " << gamma_mcEventWeight->GetBinContent(16) << " " << beta_mcEventWeight->GetBinContent(16) << endl;
cout << "flavor cut :" << alpha_mcEventWeight->GetBinContent(17) << " " << gamma_mcEventWeight->GetBinContent(17) << " " << beta_mcEventWeight->GetBinContent(17) << endl;
cout << "tau veto :" << alpha_mcEventWeight->GetBinContent(18) << " " << gamma_mcEventWeight->GetBinContent(18) << " " << beta_mcEventWeight->GetBinContent(18) << endl;
cout << "valid pT :" << alpha_mcEventWeight->GetBinContent(19) << " " << gamma_mcEventWeight->GetBinContent(19) << " " << beta_mcEventWeight->GetBinContent(19) << endl;
cout << "Trigger :" << alpha_mcEventWeight->GetBinContent(20) << " " << gamma_mcEventWeight->GetBinContent(20) << " " << beta_mcEventWeight->GetBinContent(20) << endl;
cout << "RealLeptonsCheck :" << alpha_mcEventWeight->GetBinContent(21) << " " << gamma_mcEventWeight->GetBinContent(21) << " " << beta_mcEventWeight->GetBinContent(21) << endl;
cout << "SameSign :" << alpha_mcEventWeight->GetBinContent(22) << " " << gamma_mcEventWeight->GetBinContent(22) << " " << beta_mcEventWeight->GetBinContent(22) << endl;
cout << "Iso :" << alpha_mcEventWeight->GetBinContent(23) << " " << gamma_mcEventWeight->GetBinContent(23) << " " << beta_mcEventWeight->GetBinContent(23) << endl;
cout << "dOS :" << alpha_mcEventWeight->GetBinContent(24) << " " << gamma_mcEventWeight->GetBinContent(24) << " " << beta_mcEventWeight->GetBinContent(24) << endl;
cout << "FVeto :" << alpha_mcEventWeight->GetBinContent(25) << " " << gamma_mcEventWeight->GetBinContent(25) << " " << beta_mcEventWeight->GetBinContent(25) << endl;
cout << "bVeto :" << alpha_mcEventWeight->GetBinContent(26) << " " << gamma_mcEventWeight->GetBinContent(26) << " " << beta_mcEventWeight->GetBinContent(26) << endl;
cout << ">=1jet :" << alpha_mcEventWeight->GetBinContent(27) << " " << gamma_mcEventWeight->GetBinContent(27) << " " << beta_mcEventWeight->GetBinContent(27) << endl;
cout << "lep pt :" << alpha_mcEventWeight->GetBinContent(28) << " " << gamma_mcEventWeight->GetBinContent(28) << " " << beta_mcEventWeight->GetBinContent(28) << endl;
cout << "mll :" << alpha_mcEventWeight->GetBinContent(29) << " " << gamma_mcEventWeight->GetBinContent(29) << " " << beta_mcEventWeight->GetBinContent(29) << endl;
cout << "mWWT :" << alpha_mcEventWeight->GetBinContent(30) << " " << gamma_mcEventWeight->GetBinContent(30) << " " << beta_mcEventWeight->GetBinContent(30) << endl;
cout << "HT :" << alpha_mcEventWeight->GetBinContent(31) << " " << gamma_mcEventWeight->GetBinContent(31) << " " << beta_mcEventWeight->GetBinContent(31) << endl;
cout << "metREL :" << alpha_mcEventWeight->GetBinContent(32) << " " << gamma_mcEventWeight->GetBinContent(32) << " " << beta_mcEventWeight->GetBinContent(32) << endl;
cout << "MM SRSS2: " << endl;
cout << "mWWT :" << gamma_mcEventWeight->GetBinContent(33) << endl;
cout << "HT :" << gamma_mcEventWeight->GetBinContent(34) << endl;
cout << "MM SRSS3: " << endl;
cout << "mWWT :" << gamma_mcEventWeight->GetBinContent(36) << endl;
cout << "HT :" << gamma_mcEventWeight->GetBinContent(37) << endl;
cout << "MM SRSS4: " << endl;
cout << "mWWT :" << gamma_mcEventWeight->GetBinContent(39) << endl;
cout << "HT :" << gamma_mcEventWeight->GetBinContent(40) << endl;
cout << "EE SRSS:" << endl;
cout << "mt2>90 :" << alpha_mcEventWeight->GetBinContent(33) << endl;
cout << "---------------- OS --------------------------" << endl;
cout << " EE MM EM " << endl;
cout << "----------------------------------------------" << endl;
cout << "pass OS :" << alpha_mcEventWeight->GetBinContent(51) << " " << gamma_mcEventWeight->GetBinContent(51) << " " << beta_mcEventWeight->GetBinContent(51) << endl;
cout << "pass Iso :" << alpha_mcEventWeight->GetBinContent(52) << " " << gamma_mcEventWeight->GetBinContent(52) << " " << beta_mcEventWeight->GetBinContent(52) << endl;
cout << "pass FVeto :" << alpha_mcEventWeight->GetBinContent(53) << " " << gamma_mcEventWeight->GetBinContent(53) << " " << beta_mcEventWeight->GetBinContent(53) << endl;
cout << "pass bVeto :" << alpha_mcEventWeight->GetBinContent(54) << " " << gamma_mcEventWeight->GetBinContent(54) << " " << beta_mcEventWeight->GetBinContent(54) << endl;
cout << "pass >=2 J :" << alpha_mcEventWeight->GetBinContent(55) << " " << gamma_mcEventWeight->GetBinContent(55) << " " << beta_mcEventWeight->GetBinContent(55) << endl;
cout << "pass Mjj :" << alpha_mcEventWeight->GetBinContent(56) << " " << gamma_mcEventWeight->GetBinContent(56) << " " << beta_mcEventWeight->GetBinContent(56) << endl;
cout << "pass lepPt :" << alpha_mcEventWeight->GetBinContent(57) << " " << gamma_mcEventWeight->GetBinContent(57) << " " << beta_mcEventWeight->GetBinContent(57) << endl;
cout << "pass dRll :" << alpha_mcEventWeight->GetBinContent(58) << " " << gamma_mcEventWeight->GetBinContent(58) << " " << beta_mcEventWeight->GetBinContent(58) << endl;
cout << "pass min(mT1,mT2) :" << alpha_mcEventWeight->GetBinContent(59) << " " << gamma_mcEventWeight->GetBinContent(59) << " " << beta_mcEventWeight->GetBinContent(59) << endl;
cout << "pass dPhi(ll,Met) :" << alpha_mcEventWeight->GetBinContent(60) << " " << gamma_mcEventWeight->GetBinContent(60) << " " << beta_mcEventWeight->GetBinContent(60) << endl;
cout << "pass Met :" << alpha_mcEventWeight->GetBinContent(61) << " " << gamma_mcEventWeight->GetBinContent(61) << " " << beta_mcEventWeight->GetBinContent(61) << endl;
}