SusyScan/PlotScript/SusyScan.cc

#include "SusyScan.h"
#include "ConfigFile.h"

SusyScan::SusyScan(std::string filename)
{
  ConfigFile config(filename);
  Mzero         = config.read<double>("Mzero", 0);
  Mhalf         = config.read<double>("Mhalf", 0);
  Mu            = config.read<double>("Mu", 0);
  TanBeta       = config.read<double>("TanBeta", 0);
  Azero         = config.read<double>("Azero", 0);
  Run           = config.read<double>("Run", 0);
  background    = config.read<double>("background", 0);
  signal_contamination  = config.read<double>("signal.signalregion.IsoMuon", 10);
  background_uncertainty = config.read<double>("background.uncertainty", 0);
  data          = config.read<double>("data", 0);
  signal        = config.read<double>("signal", 0);
  signal_uncertainty = config.read<double>("signal.uncertainty", 0);
  Xsection      = config.read<double>("Xsection", 0);
  ExpXsecLimit  = config.read<double>("ExpXsecLimit", 0);
  ObsXsecLimit  = config.read<double>("ObsXsecLimit", 0);
  ExpNsigLimit  = config.read<double>("ExpNsigLimit", 0);
  ObsNsigLimit  = config.read<double>("ObsNsigLimit", 0);
  lnQ_b_exp     = config.read<double>("-2lnQ_b@exp", 0);
  lnQ_b_obs     = config.read<double>("-2lnQ_b@obs", 0);
  lnQ_b_xsec    = config.read<double>("-2lnQ_b@xsec", 0);
  lnQ_sb_exp    = config.read<double>("-2lnQ_sb@exp", 0);
  lnQ_sb_obs    = config.read<double>("-2lnQ_sb@obs", 0);
  lnQ_sb_xsec   = config.read<double>("-2lnQ_sb@xsec", 0);
  CLb_b_exp     = config.read<double>("CLb_b@exp", 0);
  CLb_b_obs     = config.read<double>("CLb_b@obs", 0);
  CLb_b_xsec    = config.read<double>("CLb_b@xsec", 0);
  CLs_exp       = config.read<double>("CLs@exp", 0);
  CLs_obs       = config.read<double>("CLs@obs", 0);
  CLs_xsec      = config.read<double>("CLs@xsec", 0);
  CLs_b_exp     = config.read<double>("CLs_b@exp", 0);
  CLs_b_obs     = config.read<double>("CLs_b@obs", 0);
  CLs_b_xsec    = config.read<double>("CLs_b@xsec", 0);
  CLs_b_n1_exp  = config.read<double>("CLs_b_n1@exp", 0);
  CLs_b_n1_obs  = config.read<double>("CLs_b_n1@obs", 0);
  CLs_b_n1_xsec = config.read<double>("CLs_b_n1@xsec", 0);
  CLs_b_n2_exp  = config.read<double>("CLs_b_n2@exp", 0);
  CLs_b_n2_obs  = config.read<double>("CLs_b_n2@obs", 0);
  CLs_b_n2_xsec = config.read<double>("CLs_b_n2@xsec", 0);
  CLs_b_p1_exp  = config.read<double>("CLs_b_p1@exp", 0);
  CLs_b_p1_obs  = config.read<double>("CLs_b_p1@obs", 0);
  CLs_b_p1_xsec = config.read<double>("CLs_b_p1@xsec", 0);
  CLs_b_p2_exp  = config.read<double>("CLs_b_p2@exp", 0);
  CLs_b_p2_obs  = config.read<double>("CLs_b_p2@obs", 0);
  CLs_b_p2_xsec = config.read<double>("CLs_b_p2@xsec", 0);
  CLsb_b_exp    = config.read<double>("CLsb_b@exp", 0);
  CLsb_b_obs    = config.read<double>("CLsb_b@obs", 0);
  CLsb_b_xsec   = config.read<double>("CLsb_b@xsec", 0);
  PLExpXsecLimit= config.read<double>("RooSimpleProfile.xsec.ExpUpperLimit", 0);
  PLObsXsecLimit= config.read<double>("RooSimpleProfile.xsec.ObsUpperLimit", 0);
  PLExpNsigLimit= config.read<double>("RooSimpleProfile.signal.ExpUpperLimit", 0);
  PLObsNsigLimit= config.read<double>("RooSimpleProfile.signal.ObsUpperLimit", 0);
  FCExpXsecLimit= config.read<double>("RooFC.xsec.ExpUpperLimit", 0);
  FCObsXsecLimit= config.read<double>("RooFC.xsec.ObsUpperLimit", 0);
  FCExpNsigLimit= config.read<double>("RooFC.signal.ExpUpperLimit", 0);
  FCObsNsigLimit= config.read<double>("RooFC.signal.ObsUpperLimit", 0);
  MCMCExpXsecLimit= config.read<double>("RooMCMC.xsec.ExpUpperLimit", 0);
  MCMCObsXsecLimit= config.read<double>("RooMCMC.xsec.ObsUpperLimit", 0);
  MCMCExpNsigLimit= config.read<double>("RooMCMC.signal.ExpUpperLimit", 0);
  MCMCObsNsigLimit= config.read<double>("RooMCMC.signal.ObsUpperLimit", 0);
}
Revision:	1.5
Committed:	Fri Feb 4 15:42:08 2011 UTC (14 years, 3 months ago) by auterman
Content type:	text/plain
Branch:	MAIN
Changes since 1.4:	+1 -0 lines
Log Message:	* empty log message *
#	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	auterman	1.5	signal_contamination = config.read<double>("signal.signalregion.IsoMuon", 10);
15	auterman	1.1	background_uncertainty = config.read<double>("background.uncertainty", 0);
16			data = config.read<double>("data", 0);
17			signal = config.read<double>("signal", 0);
18			signal_uncertainty = config.read<double>("signal.uncertainty", 0);
19			Xsection = config.read<double>("Xsection", 0);
20			ExpXsecLimit = config.read<double>("ExpXsecLimit", 0);
21			ObsXsecLimit = config.read<double>("ObsXsecLimit", 0);
22	auterman	1.3	ExpNsigLimit = config.read<double>("ExpNsigLimit", 0);
23			ObsNsigLimit = config.read<double>("ObsNsigLimit", 0);
24	auterman	1.1	lnQ_b_exp = config.read<double>("-2lnQ_b@exp", 0);
25			lnQ_b_obs = config.read<double>("-2lnQ_b@obs", 0);
26	auterman	1.4	lnQ_b_xsec = config.read<double>("-2lnQ_b@xsec", 0);
27	auterman	1.1	lnQ_sb_exp = config.read<double>("-2lnQ_sb@exp", 0);
28			lnQ_sb_obs = config.read<double>("-2lnQ_sb@obs", 0);
29	auterman	1.4	lnQ_sb_xsec = config.read<double>("-2lnQ_sb@xsec", 0);
30			CLb_b_exp = config.read<double>("CLb_b@exp", 0);
31			CLb_b_obs = config.read<double>("CLb_b@obs", 0);
32			CLb_b_xsec = config.read<double>("CLb_b@xsec", 0);
33			CLs_exp = config.read<double>("CLs@exp", 0);
34			CLs_obs = config.read<double>("CLs@obs", 0);
35			CLs_xsec = config.read<double>("CLs@xsec", 0);
36			CLs_b_exp = config.read<double>("CLs_b@exp", 0);
37			CLs_b_obs = config.read<double>("CLs_b@obs", 0);
38			CLs_b_xsec = config.read<double>("CLs_b@xsec", 0);
39			CLs_b_n1_exp = config.read<double>("CLs_b_n1@exp", 0);
40			CLs_b_n1_obs = config.read<double>("CLs_b_n1@obs", 0);
41			CLs_b_n1_xsec = config.read<double>("CLs_b_n1@xsec", 0);
42			CLs_b_n2_exp = config.read<double>("CLs_b_n2@exp", 0);
43			CLs_b_n2_obs = config.read<double>("CLs_b_n2@obs", 0);
44			CLs_b_n2_xsec = config.read<double>("CLs_b_n2@xsec", 0);
45			CLs_b_p1_exp = config.read<double>("CLs_b_p1@exp", 0);
46			CLs_b_p1_obs = config.read<double>("CLs_b_p1@obs", 0);
47			CLs_b_p1_xsec = config.read<double>("CLs_b_p1@xsec", 0);
48			CLs_b_p2_exp = config.read<double>("CLs_b_p2@exp", 0);
49			CLs_b_p2_obs = config.read<double>("CLs_b_p2@obs", 0);
50			CLs_b_p2_xsec = config.read<double>("CLs_b_p2@xsec", 0);
51			CLsb_b_exp = config.read<double>("CLsb_b@exp", 0);
52			CLsb_b_obs = config.read<double>("CLsb_b@obs", 0);
53			CLsb_b_xsec = config.read<double>("CLsb_b@xsec", 0);
54	auterman	1.2	PLExpXsecLimit= config.read<double>("RooSimpleProfile.xsec.ExpUpperLimit", 0);
55			PLObsXsecLimit= config.read<double>("RooSimpleProfile.xsec.ObsUpperLimit", 0);
56			PLExpNsigLimit= config.read<double>("RooSimpleProfile.signal.ExpUpperLimit", 0);
57			PLObsNsigLimit= config.read<double>("RooSimpleProfile.signal.ObsUpperLimit", 0);
58	auterman	1.3	FCExpXsecLimit= config.read<double>("RooFC.xsec.ExpUpperLimit", 0);
59			FCObsXsecLimit= config.read<double>("RooFC.xsec.ObsUpperLimit", 0);
60			FCExpNsigLimit= config.read<double>("RooFC.signal.ExpUpperLimit", 0);
61			FCObsNsigLimit= config.read<double>("RooFC.signal.ObsUpperLimit", 0);
62			MCMCExpXsecLimit= config.read<double>("RooMCMC.xsec.ExpUpperLimit", 0);
63			MCMCObsXsecLimit= config.read<double>("RooMCMC.xsec.ObsUpperLimit", 0);
64			MCMCExpNsigLimit= config.read<double>("RooMCMC.signal.ExpUpperLimit", 0);
65			MCMCObsNsigLimit= config.read<double>("RooMCMC.signal.ObsUpperLimit", 0);
66	auterman	1.1	}