| 1 |
#include "SusyScan.h"
|
| 2 |
#include "ConfigFile.h"
|
| 3 |
|
| 4 |
SusyScan::SusyScan(std::string filename)
|
| 5 |
{
|
| 6 |
ConfigFile config(filename);
|
| 7 |
Mzero = config.read<double>("Mzero", 0);
|
| 8 |
Mhalf = config.read<double>("Mhalf", 0);
|
| 9 |
Mu = config.read<double>("Mu", 0);
|
| 10 |
TanBeta = config.read<double>("TanBeta", 0);
|
| 11 |
Azero = config.read<double>("Azero", 0);
|
| 12 |
Run = config.read<double>("Run", 0);
|
| 13 |
background = config.read<double>("background", 0);
|
| 14 |
background_uncertainty = config.read<double>("background.uncertainty", 0);
|
| 15 |
data = config.read<double>("data", 0);
|
| 16 |
signal_JEC_UP = config.read<double>("signal.JEC_.UP", 0);
|
| 17 |
signal_JEC_DN = config.read<double>("signal.JEC_.DN", 0);
|
| 18 |
signal_MuIso_UP = config.read<double>("signal.MuIso_.UP", 0);
|
| 19 |
signal_MuIso_DN = config.read<double>("signal.MuIso_.DN", 0);
|
| 20 |
signal_kfactor_UP = config.read<double>("signal.kFactorUP", 0);
|
| 21 |
signal_kfactor_DN = config.read<double>("signal.kFactorDN", 0);
|
| 22 |
signal_kfactor = config.read<double>("signal.kFactor", 0);
|
| 23 |
|
| 24 |
signal_contamination = config.read<double>("signal.LO.signalregion.IsoMuon", 10);
|
| 25 |
signal = config.read<double>("signal.LO", 0);
|
| 26 |
signal_uncertainty = config.read<double>("signal.LO.uncertainty", 0);
|
| 27 |
Xsection = config.read<double>("Xsection", 0);
|
| 28 |
ExpXsecLimit = config.read<double>("LO.ExpXsecLimit", 0);
|
| 29 |
ObsXsecLimit = config.read<double>("LO.ObsXsecLimit", 0);
|
| 30 |
ExpNsigLimit = config.read<double>("LO.ExpNsigLimit", 0);
|
| 31 |
ObsNsigLimit = config.read<double>("LO.ObsNsigLimit", 0);
|
| 32 |
lnQ_b_exp = config.read<double>("LO.-2lnQ_b@exp", 0);
|
| 33 |
lnQ_b_obs = config.read<double>("LO.-2lnQ_b@obs", 0);
|
| 34 |
lnQ_b_xsec = config.read<double>("LO.-2lnQ_b@xsec", 0);
|
| 35 |
lnQ_sb_exp = config.read<double>("LO.-2lnQ_sb@exp", 0);
|
| 36 |
lnQ_sb_obs = config.read<double>("LO.-2lnQ_sb@obs", 0);
|
| 37 |
lnQ_sb_xsec = config.read<double>("LO.-2lnQ_sb@xsec", 0);
|
| 38 |
CLb_b_exp = config.read<double>("LO.CLb_b@exp", 0);
|
| 39 |
CLb_b_obs = config.read<double>("LO.CLb_b@obs", 0);
|
| 40 |
CLb_b_xsec = config.read<double>("LO.CLb_b@xsec", 0);
|
| 41 |
CLs_exp = config.read<double>("LO.CLs@exp", 0);
|
| 42 |
CLs_obs = config.read<double>("LO.CLs@obs", 0);
|
| 43 |
CLs_xsec = config.read<double>("LO.CLs@xsec", 0);
|
| 44 |
CLs_b_exp = config.read<double>("LO.CLs_b@exp", 0);
|
| 45 |
CLs_b_obs = config.read<double>("LO.CLs_b@obs", 0);
|
| 46 |
CLs_b_xsec = config.read<double>("LO.CLs_b@xsec", 0);
|
| 47 |
CLs_b_n1_exp = config.read<double>("LO.CLs_b_n1@exp", 0);
|
| 48 |
CLs_b_n1_obs = config.read<double>("LO.CLs_b_n1@obs", 0);
|
| 49 |
CLs_b_n1_xsec = config.read<double>("LO.CLs_b_n1@xsec", 0);
|
| 50 |
CLs_b_n2_exp = config.read<double>("LO.CLs_b_n2@exp", 0);
|
| 51 |
CLs_b_n2_obs = config.read<double>("LO.CLs_b_n2@obs", 0);
|
| 52 |
CLs_b_n2_xsec = config.read<double>("LO.CLs_b_n2@xsec", 0);
|
| 53 |
CLs_b_p1_exp = config.read<double>("LO.CLs_b_p1@exp", 0);
|
| 54 |
CLs_b_p1_obs = config.read<double>("LO.CLs_b_p1@obs", 0);
|
| 55 |
CLs_b_p1_xsec = config.read<double>("LO.CLs_b_p1@xsec", 0);
|
| 56 |
CLs_b_p2_exp = config.read<double>("LO.CLs_b_p2@exp", 0);
|
| 57 |
CLs_b_p2_obs = config.read<double>("LO.CLs_b_p2@obs", 0);
|
| 58 |
CLs_b_p2_xsec = config.read<double>("LO.CLs_b_p2@xsec", 0);
|
| 59 |
CLsb_b_exp = config.read<double>("LO.CLsb_b@exp", 0);
|
| 60 |
CLsb_b_obs = config.read<double>("LO.CLsb_b@obs", 0);
|
| 61 |
CLsb_b_xsec = config.read<double>("LO.CLsb_b@xsec", 0);
|
| 62 |
PLExpXsecLimit= config.read<double>("RooSimpleProfile.xsec.LO.ExpUpperLimit", 0);
|
| 63 |
PLObsXsecLimit= config.read<double>("RooSimpleProfile.xsec.LO.ObsUpperLimit", 0);
|
| 64 |
PLExpNsigLimit= config.read<double>("RooSimpleProfile.signal.LO.ExpUpperLimit", 0);
|
| 65 |
PLObsNsigLimit= config.read<double>("RooSimpleProfile.signal.LO.ObsUpperLimit", 0);
|
| 66 |
FCExpXsecLimit= config.read<double>("RooFC.xsec.LO.ExpUpperLimit", 0);
|
| 67 |
FCObsXsecLimit= config.read<double>("RooFC.xsec.LO.ObsUpperLimit", 0);
|
| 68 |
FCExpNsigLimit= config.read<double>("RooFC.signal.LO.ExpUpperLimit", 0);
|
| 69 |
FCObsNsigLimit= config.read<double>("RooFC.signal.LO.ObsUpperLimit", 0);
|
| 70 |
MCMCExpXsecLimit= config.read<double>("RooMCMC.xsec.LO.ExpUpperLimit", 0);
|
| 71 |
MCMCObsXsecLimit= config.read<double>("RooMCMC.xsec.LO.ObsUpperLimit", 0);
|
| 72 |
MCMCExpNsigLimit= config.read<double>("RooMCMC.signal.LO.ExpUpperLimit", 0);
|
| 73 |
MCMCObsNsigLimit= config.read<double>("RooMCMC.signal.LO.ObsUpperLimit", 0);
|
| 74 |
|
| 75 |
NLO_ExpXsecLimit = config.read<double>("NLO.ExpXsecLimit", 0);
|
| 76 |
NLO_ObsXsecLimit = config.read<double>("NLO.ObsXsecLimit", 0);
|
| 77 |
NLO_ExpNsigLimit = config.read<double>("NLO.ExpNsigLimit", 0);
|
| 78 |
NLO_ObsNsigLimit = config.read<double>("NLO.ObsNsigLimit", 0);
|
| 79 |
NLO_lnQ_b_exp = config.read<double>("NLO.-2lnQ_b@exp", 0);
|
| 80 |
NLO_lnQ_b_obs = config.read<double>("NLO.-2lnQ_b@obs", 0);
|
| 81 |
NLO_lnQ_b_xsec = config.read<double>("NLO.-2lnQ_b@xsec", 0);
|
| 82 |
NLO_lnQ_sb_exp = config.read<double>("NLO.-2lnQ_sb@exp", 0);
|
| 83 |
NLO_lnQ_sb_obs = config.read<double>("NLO.-2lnQ_sb@obs", 0);
|
| 84 |
NLO_lnQ_sb_xsec = config.read<double>("NLO.-2lnQ_sb@xsec", 0);
|
| 85 |
NLO_CLb_b_exp = config.read<double>("NLO.CLb_b@exp", 0);
|
| 86 |
NLO_CLb_b_obs = config.read<double>("NLO.CLb_b@obs", 0);
|
| 87 |
NLO_CLb_b_xsec = config.read<double>("NLO.CLb_b@xsec", 0);
|
| 88 |
NLO_CLs_exp = config.read<double>("NLO.CLs@exp", 0);
|
| 89 |
NLO_CLs_obs = config.read<double>("NLO.CLs@obs", 0);
|
| 90 |
NLO_CLs_xsec = config.read<double>("NLO.CLs@xsec", 0);
|
| 91 |
NLO_CLs_b_exp = config.read<double>("NLO.CLs_b@exp", 0);
|
| 92 |
NLO_CLs_b_obs = config.read<double>("NLO.CLs_b@obs", 0);
|
| 93 |
NLO_CLs_b_xsec = config.read<double>("NLO.CLs_b@xsec", 0);
|
| 94 |
NLO_CLs_b_n1_exp = config.read<double>("NLO.CLs_b_n1@exp", 0);
|
| 95 |
NLO_CLs_b_n1_obs = config.read<double>("NLO.CLs_b_n1@obs", 0);
|
| 96 |
NLO_CLs_b_n1_xsec = config.read<double>("NLO.CLs_b_n1@xsec", 0);
|
| 97 |
NLO_CLs_b_n2_exp = config.read<double>("NLO.CLs_b_n2@exp", 0);
|
| 98 |
NLO_CLs_b_n2_obs = config.read<double>("NLO.CLs_b_n2@obs", 0);
|
| 99 |
NLO_CLs_b_n2_xsec = config.read<double>("NLO.CLs_b_n2@xsec", 0);
|
| 100 |
NLO_CLs_b_p1_exp = config.read<double>("NLO.CLs_b_p1@exp", 0);
|
| 101 |
NLO_CLs_b_p1_obs = config.read<double>("NLO.CLs_b_p1@obs", 0);
|
| 102 |
NLO_CLs_b_p1_xsec = config.read<double>("NLO.CLs_b_p1@xsec", 0);
|
| 103 |
NLO_CLs_b_p2_exp = config.read<double>("NLO.CLs_b_p2@exp", 0);
|
| 104 |
NLO_CLs_b_p2_obs = config.read<double>("NLO.CLs_b_p2@obs", 0);
|
| 105 |
NLO_CLs_b_p2_xsec = config.read<double>("NLO.CLs_b_p2@xsec", 0);
|
| 106 |
NLO_CLsb_b_exp = config.read<double>("NLO.CLsb_b@exp", 0);
|
| 107 |
NLO_CLsb_b_obs = config.read<double>("NLO.CLsb_b@obs", 0);
|
| 108 |
NLO_CLsb_b_xsec = config.read<double>("NLO.CLsb_b@xsec", 0);
|
| 109 |
NLO_PLExpXsecLimit= config.read<double>("RooSimpleProfile.xsec.NLO.ExpUpperLimit", 0);
|
| 110 |
NLO_PLObsXsecLimit= config.read<double>("RooSimpleProfile.xsec.NLO.ObsUpperLimit", 0);
|
| 111 |
NLO_PLExpNsigLimit= config.read<double>("RooSimpleProfile.signal.NLO.ExpUpperLimit", 0);
|
| 112 |
NLO_PLObsNsigLimit= config.read<double>("RooSimpleProfile.signal.NLO.ObsUpperLimit", 0);
|
| 113 |
NLO_FCExpXsecLimit= config.read<double>("RooFC.xsec.NLO.ExpUpperLimit", 0);
|
| 114 |
NLO_FCObsXsecLimit= config.read<double>("RooFC.xsec.NLO.ObsUpperLimit", 0);
|
| 115 |
NLO_FCExpNsigLimit= config.read<double>("RooFC.signal.NLO.ExpUpperLimit", 0);
|
| 116 |
NLO_FCObsNsigLimit= config.read<double>("RooFC.signal.NLO.ObsUpperLimit", 0);
|
| 117 |
NLO_MCMCExpXsecLimit= config.read<double>("RooMCMC.xsec.NLO.ExpUpperLimit", 0);
|
| 118 |
NLO_MCMCObsXsecLimit= config.read<double>("RooMCMC.xsec.NLO.ObsUpperLimit", 0);
|
| 119 |
NLO_MCMCExpNsigLimit= config.read<double>("RooMCMC.signal.NLO.ExpUpperLimit", 0);
|
| 120 |
NLO_MCMCObsNsigLimit= config.read<double>("RooMCMC.signal.NLO.ObsUpperLimit", 0);
|
| 121 |
}
|
