ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/SusyScan.cc
Revision: 1.6
Committed: Tue Feb 22 13:59:19 2011 UTC (14 years, 2 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.5: +104 -49 lines
Log Message: 
NLO updates

File Contents

# User Rev Content
1 auterman 1.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 auterman 1.6 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 auterman 1.1 Xsection = config.read<double>("Xsection", 0);
28 auterman 1.6 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 auterman 1.1 }