models_mtot_range_m450.rs

mtot[300,1200];

mtot_sig_m0_cat0[450, 400, 500];
mtot_sig_sigma_cat0[10, 5, 40.0];
mtot_sig_alpha_cat0[-1.0, -3, -0.5];
mtot_sig_n_cat0[4.0, 0.5, 10];
mtot_sig_gsigma_cat0[40, 20., 60.0];
mtot_sig_frac_cat0[0.2, 0., 0.6];

mtot_sig_m0_cat1[455, 400, 500];
mtot_sig_sigma_cat1[10, 0.5, 20.0];
mtot_sig_alpha_cat1[-2.0, -3, 3];
mtot_sig_n_cat1[5.0, 0.4, 10];
mtot_sig_gsigma_cat1[5.0, 1.0, 10.0];
mtot_sig_frac_cat1[0.7, 0.4, 1.0];

mtotGaussSig_cat0 = Gaussian(mtot, mtot_sig_m0_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 = Gaussian(mtot, mtot_sig_m0_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.]