ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/SusyScan.cc
(Generate patch)

Comparing UserCode/auterman/SusyScan/PlotScript/SusyScan.cc (file contents):
Revision 1.5 by auterman, Fri Feb 4 15:42:08 2011 UTC vs.
Revision 1.6 by auterman, Tue Feb 22 13:59:19 2011 UTC

# Line 11 | Line 11 | SusyScan::SusyScan(std::string filename)
11    Azero         = config.read<double>("Azero", 0);
12    Run           = config.read<double>("Run", 0);
13    background    = config.read<double>("background", 0);
14  signal_contamination  = config.read<double>("signal.signalregion.IsoMuon", 10);
14    background_uncertainty = config.read<double>("background.uncertainty", 0);
15    data          = config.read<double>("data", 0);
16 <  signal        = config.read<double>("signal", 0);
17 <  signal_uncertainty = config.read<double>("signal.uncertainty", 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>("ExpXsecLimit", 0);
29 <  ObsXsecLimit  = config.read<double>("ObsXsecLimit", 0);
30 <  ExpNsigLimit  = config.read<double>("ExpNsigLimit", 0);
31 <  ObsNsigLimit  = config.read<double>("ObsNsigLimit", 0);
32 <  lnQ_b_exp     = config.read<double>("-2lnQ_b@exp", 0);
33 <  lnQ_b_obs     = config.read<double>("-2lnQ_b@obs", 0);
34 <  lnQ_b_xsec    = config.read<double>("-2lnQ_b@xsec", 0);
35 <  lnQ_sb_exp    = config.read<double>("-2lnQ_sb@exp", 0);
36 <  lnQ_sb_obs    = config.read<double>("-2lnQ_sb@obs", 0);
37 <  lnQ_sb_xsec   = config.read<double>("-2lnQ_sb@xsec", 0);
38 <  CLb_b_exp     = config.read<double>("CLb_b@exp", 0);
39 <  CLb_b_obs     = config.read<double>("CLb_b@obs", 0);
40 <  CLb_b_xsec    = config.read<double>("CLb_b@xsec", 0);
41 <  CLs_exp       = config.read<double>("CLs@exp", 0);
42 <  CLs_obs       = config.read<double>("CLs@obs", 0);
43 <  CLs_xsec      = config.read<double>("CLs@xsec", 0);
44 <  CLs_b_exp     = config.read<double>("CLs_b@exp", 0);
45 <  CLs_b_obs     = config.read<double>("CLs_b@obs", 0);
46 <  CLs_b_xsec    = config.read<double>("CLs_b@xsec", 0);
47 <  CLs_b_n1_exp  = config.read<double>("CLs_b_n1@exp", 0);
48 <  CLs_b_n1_obs  = config.read<double>("CLs_b_n1@obs", 0);
49 <  CLs_b_n1_xsec = config.read<double>("CLs_b_n1@xsec", 0);
50 <  CLs_b_n2_exp  = config.read<double>("CLs_b_n2@exp", 0);
51 <  CLs_b_n2_obs  = config.read<double>("CLs_b_n2@obs", 0);
52 <  CLs_b_n2_xsec = config.read<double>("CLs_b_n2@xsec", 0);
53 <  CLs_b_p1_exp  = config.read<double>("CLs_b_p1@exp", 0);
54 <  CLs_b_p1_obs  = config.read<double>("CLs_b_p1@obs", 0);
55 <  CLs_b_p1_xsec = config.read<double>("CLs_b_p1@xsec", 0);
56 <  CLs_b_p2_exp  = config.read<double>("CLs_b_p2@exp", 0);
57 <  CLs_b_p2_obs  = config.read<double>("CLs_b_p2@obs", 0);
58 <  CLs_b_p2_xsec = config.read<double>("CLs_b_p2@xsec", 0);
59 <  CLsb_b_exp    = config.read<double>("CLsb_b@exp", 0);
60 <  CLsb_b_obs    = config.read<double>("CLsb_b@obs", 0);
61 <  CLsb_b_xsec   = config.read<double>("CLsb_b@xsec", 0);
62 <  PLExpXsecLimit= config.read<double>("RooSimpleProfile.xsec.ExpUpperLimit", 0);
63 <  PLObsXsecLimit= config.read<double>("RooSimpleProfile.xsec.ObsUpperLimit", 0);
64 <  PLExpNsigLimit= config.read<double>("RooSimpleProfile.signal.ExpUpperLimit", 0);
65 <  PLObsNsigLimit= config.read<double>("RooSimpleProfile.signal.ObsUpperLimit", 0);
66 <  FCExpXsecLimit= config.read<double>("RooFC.xsec.ExpUpperLimit", 0);
67 <  FCObsXsecLimit= config.read<double>("RooFC.xsec.ObsUpperLimit", 0);
68 <  FCExpNsigLimit= config.read<double>("RooFC.signal.ExpUpperLimit", 0);
69 <  FCObsNsigLimit= config.read<double>("RooFC.signal.ObsUpperLimit", 0);
70 <  MCMCExpXsecLimit= config.read<double>("RooMCMC.xsec.ExpUpperLimit", 0);
71 <  MCMCObsXsecLimit= config.read<double>("RooMCMC.xsec.ObsUpperLimit", 0);
72 <  MCMCExpNsigLimit= config.read<double>("RooMCMC.signal.ExpUpperLimit", 0);
73 <  MCMCObsNsigLimit= config.read<double>("RooMCMC.signal.ObsUpperLimit", 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   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines