forked from acarvalh/Limit_codes
-
Notifications
You must be signed in to change notification settings - Fork 1
/
models_mtot_range.rs
33 lines (26 loc) · 1.22 KB
/
models_mtot_range.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
mtot[300,1200];
mtot_sig_m0_cat0[1100, 1000, 1200];
mtot_sig_sigma_cat0[30, 5, 120.0];
mtot_sig_alpha_cat0[1.0, 0.5, 3];
mtot_sig_n_cat0[4.0, 0.5, 10];
mtot_sig_g_cat0[1.0, 0.5, 3];
mtot_sig_gsigma_cat0[20, 10., 120.0];
mtot_sig_frac_cat0[0.4, 0.0, 0.8];
mtot_sig_m0_cat1[1100, 1000, 1200];
mtot_sig_sigma_cat1[20, 5, 120.0];
mtot_sig_alpha_cat1[1.0, 0.5, 3];
mtot_sig_n_cat1[4.0, 0.5, 10];
mtot_sig_g_cat1[1.0, 0.5, 3];
mtot_sig_gsigma_cat1[30.0, 10.0, 120.0];
mtot_sig_frac_cat1[0.8, 0.0, 1];
mtotGaussSig_cat0 = Voigtian(mtot, mtot_sig_m0_cat0,mtot_sig_g_cat0, mtot_sig_gsigma_cat0);
mtotCBSig_cat0 = CBShape(mtot, mtot_sig_m0_cat0, mtot_sig_sigma_cat0, mtot_sig_alpha_cat0, mtot_sig_n_cat0);
mtotSig_cat0 = AddPdf(mtotGaussSig_cat0, mtotCBSig_cat0, mtot_sig_frac_cat0);
mtotGaussSig_cat1 = Voigtian(mtot, mtot_sig_m0_cat1,mtot_sig_g_cat1, mtot_sig_gsigma_cat1);
mtotCBSig_cat1 = CBShape(mtot, mtot_sig_m0_cat1, mtot_sig_sigma_cat1, mtot_sig_alpha_cat1, mtot_sig_n_cat1);
mtotSig_cat1 = AddPdf(mtotGaussSig_cat1, mtotCBSig_cat1, mtot_sig_frac_cat1);
mtot_bkg_8TeV_slope1_cat0[1.9, -10., 10.];
mtot_bkg_8TeV_slope1_cat1[1.9, -10., 10.];
mtot_bkg_8TeV_slope2_cat1[1000,800, 1200];
wei[1,0,10];
sqrtS[8000., 8000., 8000., 8000.]