-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathHighMass-hgg_m200_models_Bkg_8TeV_test.rs
265 lines (163 loc) · 8.28 KB
/
HighMass-hgg_m200_models_Bkg_8TeV_test.rs
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
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
PhotonsMass[0,500];
nvtx[0,50];
evweight[0, 1000];
PhotonsMass_sig_m0[200, 197.0, 203.0];
PhotonsMass_sig_sigma0[2., 0.0, 5.0];
PhotonsMass_sig_m1[200, 197.0, 203.0];
PhotonsMass_sig_sigma1[2.5, 0.0, 5.0];
PhotonsMass_sig_frac[0.5, 0.0, 1.0];
PhotonsMassGaussSig = Gaussian(PhotonsMass, PhotonsMass_sig_m0, PhotonsMass_sig_sigma0);
PhotonsMassGaussSig_bis = Gaussian(PhotonsMass, PhotonsMass_sig_m1, PhotonsMass_sig_sigma1);
PhotonsMassSig = AddPdf(PhotonsMassGaussSig, PhotonsMassGaussSig_bis,PhotonsMass_sig_frac);
PhotonsMass_sig_m0_cat0[199.926, 198.5, 201.0];
PhotonsMass_sig_sigma0_cat0[2.36, 2., 2.8];
PhotonsMass_sig_m1_cat0[199.926, 198.5, 201.0];
PhotonsMass_sig_sigma1_cat0[7.18, 6.7, 7.8];
PhotonsMass_sig_frac_cat0[0.9, 0.7, 1.];
PhotonsMassGaussSig_cat0 = Gaussian(PhotonsMass, PhotonsMass_sig_m0_cat0, PhotonsMass_sig_sigma0_cat0);
PhotonsMassGaussSig_cat0_bis = Gaussian(PhotonsMass, PhotonsMass_sig_m1_cat0, PhotonsMass_sig_sigma1_cat0);
#PhotonsMassSig_cat0 = AddPdf(PhotonsMassGaussSig_cat0,PhotonsMassGaussSig_cat0_bis,PhotonsMass_sig_frac_cat0);
PhotonsMass_sig_m0_cat1[199.936, 198.5, 201.0];
PhotonsMass_sig_sigma0_cat1[2.514, 2., 3.0];
PhotonsMass_sig_m1_cat1[198.268, 198.5, 201.0];
PhotonsMass_sig_sigma1_cat1[7.53, 6.5, 8.5];
PhotonsMass_sig_frac_cat1[0.862, 0.7, 1.];
PhotonsMassGaussSig_cat1 = Gaussian(PhotonsMass, PhotonsMass_sig_m0_cat1, PhotonsMass_sig_sigma0_cat1);
PhotonsMassGaussSig_cat1_bis = Gaussian(PhotonsMass, PhotonsMass_sig_m1_cat1, PhotonsMass_sig_sigma1_cat1);
#PhotonsMassSig_cat1 = AddPdf(PhotonsMassGaussSig_cat1,PhotonsMassGaussSig_cat1_bis,PhotonsMass_sig_frac_cat1);
PhotonsMass_sig_m0_cat2[199.714, 198.5, 201.0];
PhotonsMass_sig_sigma0_cat2[5.99, 5.2, 6.4];
PhotonsMass_sig_m1_cat2[199.952, 197.0, 203.0];
PhotonsMass_sig_sigma1_cat2[3.6, 3.1, 4.];
PhotonsMass_sig_frac_cat2[0.52, 0.5, 0.7];
PhotonsMassGaussSig_cat2 = Gaussian(PhotonsMass, PhotonsMass_sig_m0_cat2, PhotonsMass_sig_sigma0_cat2);
PhotonsMassGaussSig_cat2_bis = Gaussian(PhotonsMass, PhotonsMass_sig_m1_cat2, PhotonsMass_sig_sigma1_cat2);
#PhotonsMassSig_cat2 = AddPdf(PhotonsMassGaussSig_cat2,PhotonsMassGaussSig_cat2_bis,PhotonsMass_sig_frac_cat2 );
PhotonsMass_sig_m0_cat3[199.953, 199.5, 200.5];
PhotonsMass_sig_sigma0_cat3[4.62, 4., 5.2];
PhotonsMass_sig_m1_cat3[199.953, 199.5, 200.5];
PhotonsMass_sig_sigma1_cat3[13.82, 12.5, 15.];
PhotonsMass_sig_frac_cat3[0.8, 0.7, 1.];
PhotonsMassGaussSig_cat3 = Gaussian(PhotonsMass, PhotonsMass_sig_m0_cat3, PhotonsMass_sig_sigma0_cat3);
PhotonsMassGaussSig_cat3_bis = Gaussian(PhotonsMass, PhotonsMass_sig_m1_cat3, PhotonsMass_sig_sigma1_cat3);
#PhotonsMassSig_cat3 = AddPdf(PhotonsMassGaussSig_cat3,PhotonsMassGaussSig_cat3_bis,PhotonsMass_sig_frac_cat3);
nsig_cat0[50, 0, 100];
PhotonsMassSigExt_cat0 = RooExtendPdf(PhotonsMassSig_cat0, nsig_cat0);
nsig_cat1[50, 0, 100];
PhotonsMassSigExt_cat1 = RooExtendPdf(PhotonsMassSig_cat1, nsig_cat1);
nsig_cat2[50, 0, 100];
PhotonsMassSigExt_cat2 = RooExtendPdf(PhotonsMassSig_cat2, nsig_cat2);
nsig_cat3[50, 0, 100];
PhotonsMassSigExt_cat3 = RooExtendPdf(PhotonsMassSig_cat3, nsig_cat3);
PhotonsMass_sig_alphaCB_cat0[1.57, 1., 2.];
PhotonsMass_sig_alphaC_cat0[0.119, 0.08, 0.16];
PhotonsMass_sig_mean_cat0[199.81, 199., 200.5];
PhotonsMass_sig_n_cat0[9.24, 7., 11];
PhotonsMass_sig_sigma_cat0[2.2, 1.4,3.];
PhotonsMass_sig_alphaCB_cat1[1.59, 1.2, 1.8];
PhotonsMass_sig_alphaC_cat1[0.136, 0.08, 0.16];
PhotonsMass_sig_mean_cat1[199.81, 199., 200.5];
PhotonsMass_sig_n_cat1[14.24, 10., 15.];
PhotonsMass_sig_sigma_cat1[2.49, 1.7,3.1];
PhotonsMass_sig_alphaCB_cat2[1.77, 1.1, 2.];
PhotonsMass_sig_alphaC_cat2[0.07, 0.03, 0.12];
PhotonsMass_sig_mean_cat2[199.81, 199., 200.5];
PhotonsMass_sig_n_cat2[47.24, 30., 50];
PhotonsMass_sig_sigma_cat2[4.42, 3.8,5.];
PhotonsMass_sig_alphaCB_cat3[1.8, 1., 2.];
PhotonsMass_sig_alphaC_cat3[0.095, 0.05, 0.13];
PhotonsMass_sig_mean_cat3[199.81, 199., 200.5];
PhotonsMass_sig_n_cat3[53.24, 40., 60.];
PhotonsMass_sig_sigma_cat3[4.42, 3.5,5.4];
ReducedMass_sig_alphaCB_cat0[2.1, 1., 3.];
ReducedMass_sig_alphaC_cat0[0.063, 0., 0.2];
ReducedMass_sig_mean_cat0[0., -0.1, 0.1];
ReducedMass_sig_n_cat0[4.59, 0., 10];
ReducedMass_sig_sigma_cat0[0.012 0.,3.];
ReducedMass_sig_alphaCB_cat1[2.23, 1., 3.5];
ReducedMass_sig_alphaC_cat1[0.0675, 0., 0.2];
ReducedMass_sig_mean_cat1[0., -0.1, 0.1];
ReducedMass_sig_n_cat1[4.26, 0., 10];
ReducedMass_sig_sigma_cat1[0.014, 0.,3.];
ReducedMass_sig_alphaCB_cat2[2.38, 1., 3.5];
ReducedMass_sig_alphaC_cat2[0.0611, 0., 0.2];
ReducedMass_sig_mean_cat2[0., -0.1, 0.1];
ReducedMass_sig_n_cat2[3.21, 0., 10];
ReducedMass_sig_sigma_cat2[0.021, 0.,3..];
ReducedMass_sig_alphaCB_cat3[2.35, 1., 3.5];
ReducedMass_sig_alphaC_cat3[0.064, 0., 0.2];
ReducedMass_sig_mean_cat3[0., -0.1, 0.1];
ReducedMass_sig_n_cat3[3.55, 0., 10];
ReducedMass_sig_sigma_cat3[0.022, 0.,3.];
nsigCBC_cat0[50, 0, 100];
nsigCBC_cat1[50, 0, 100];
nsigCBC_cat2[50, 0, 100];
nsigCBC_cat3[50, 0, 100];
wei[1,0,10];
#PhotonsMass_bkg_exp[-0.026,-0.05, 0.];
#PhotonsMassBkg = Exponential(PhotonsMass, PhotonsMass_bkg_exp);
#PhotonsMass_bkg_exp_cat0[-0.018,-0.05, 0.];
#PhotonsMassBkg_cat0_exp = Exponential(PhotonsMass, PhotonsMass_bkg_exp_cat0);
#PhotonsMass_bkg_exp_cat1[-0.016,-0.05, 0.];
#PhotonsMassBkg_cat1_exp = Exponential(PhotonsMass, PhotonsMass_bkg_exp_cat1);
#PhotonsMass_bkg_exp_cat2[-0.015,-0.05, 0.];
#PhotonsMassBkg_cat2_exp = Exponential(PhotonsMass, PhotonsMass_bkg_exp_cat2);
#PhotonsMass_bkg_exp_cat3[-0.013,-0.05, 0.];
#PhotonsMassBkg_cat3_exp = Exponential(PhotonsMass, PhotonsMass_bkg_exp_cat3);
PhotonsMass_bkg_8TeV_c0_cat0[4.75, 3., 8.];
PhotonsMass_bkg_8TeV_c1_cat0[0.48, -2.,2.];
PhotonsMass_bkg_8TeV_c2_cat0[0.93, 0.,2.];
PhotonsMass_bkg_8TeV_c3_cat0[0.32, 0.,2.];
PhotonsMass_bkg_8TeV_c0_cat1[5., 3., 8.];
PhotonsMass_bkg_8TeV_c1_cat1[0.47, -2.,2.];
PhotonsMass_bkg_8TeV_c2_cat1[0.8, -2.,2.];
PhotonsMass_bkg_8TeV_c3_cat1[0.32, 0.,2.];
PhotonsMass_bkg_8TeV_c0_cat2[4.84, 2., 10.];
PhotonsMass_bkg_8TeV_c1_cat2[0.48, 0.,3.];
PhotonsMass_bkg_8TeV_c2_cat2[0.93, 0.,3.];
PhotonsMass_bkg_8TeV_c3_cat2[0.32, 0.,3.];
PhotonsMass_bkg_8TeV_c0_cat3[4.84, 2., 10.];
PhotonsMass_bkg_8TeV_c1_cat3[0.48, 0.,3.];
PhotonsMass_bkg_8TeV_c2_cat3[0.93, 0.,3.];
PhotonsMass_bkg_8TeV_c3_cat3[0.32, 0.,3.];
PhotonsMass_bkg_8TeV_slope1_cat0[5.7,1., 20.];
PhotonsMass_bkg_8TeV_slope2_cat0[0.000002, 0., 0.00005];
PhotonsMass_bkg_8TeV_slope3_cat0[0.24, 0., 1.];
PhotonsMass_bkg_8TeV_norm_cat0[ 14739. ,1.e+04,1.e+05];
PhotonsMass_bkg_8TeV_slope1_cat1[9.3,1., 20.];
PhotonsMass_bkg_8TeV_slope2_cat1[0.000002, 0., 0.00005];
PhotonsMass_bkg_8TeV_slope3_cat1[0.6, 0., 1.];
PhotonsMass_bkg_8TeV_norm_cat1[ 25418. ,2.e+04,2.e+05];
PhotonsMass_bkg_8TeV_slope1_cat2[8.4,1., 20.];
PhotonsMass_bkg_8TeV_slope2_cat2[0.000002, 0., 0.00005];
PhotonsMass_bkg_8TeV_slope3_cat2[0.6, 0., 1.];
PhotonsMass_bkg_8TeV_norm_cat2[ 15416. ,1.e+04,1.e+05];
PhotonsMass_bkg_8TeV_slope1_cat3[9.44,3., 20.];
PhotonsMass_bkg_8TeV_slope2_cat3[0.000002, 0., 0.00005];
PhotonsMass_bkg_8TeV_slope3_cat3[0.24, 0., 1.];
PhotonsMass_bkg_8TeV_norm_cat3[ 34022 ,3.e+04,3.e+05];
PhotonsMass_bkg_8TeV_expol1_cat0[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_cat0[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_cat0[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_cat1[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_cat1[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_cat1[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_cat2[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_cat2[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_cat2[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_cat3[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_cat3[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_cat3[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_2_cat0[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_2_cat0[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_2_cat0[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_2_cat1[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_2_cat1[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_2_cat1[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_2_cat2[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_2_cat2[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_2_cat2[0.24, 0., 10.];
PhotonsMass_bkg_8TeV_expol1_2_cat3[5.7,0 ., 40.];
PhotonsMass_bkg_8TeV_expol2_2_cat3[0.000002, 0.,20];
PhotonsMass_bkg_8TeV_expol3_2_cat3[0.24, 0., 10.];
sqrtS[8000., 8000., 8000.]