AnalysisFramework/DistributedModelCalculations/last_configuration.C

#include <iostream>
#include <vector>
#ifndef SampleClassLoaded
#include "SampleClass.C"
#endif
#define SetupLoaded
#ifndef ResultLibraryClassLoaded
#include "ResultLibraryClass.C"
#endif

using namespace std;

namespace PlottingSetup { 
string datajzb="datajzb_ERROR";
string mcjzb="mcjzb_ERROR";
vector<float>jzb_cuts;
float MCPeakError=-999;
}

void read_config() {
datajzb="(jzb[1]+0.06*pt-3.58612)";
mcjzb="(jzb[1]+0.04*pt-2.32211)";


MCPeakError=0.0990011;

jzb_cuts.push_back(50); // JZB cut at 50
jzb_cuts.push_back(100); // JZB cut at 100
jzb_cuts.push_back(150); // JZB cut at 150
jzb_cuts.push_back(200); // JZB cut at 200
jzb_cuts.push_back(250); // JZB cut at 250


Nobs.push_back(408); // JZB cut at 50
Nobs.push_back(88); // JZB cut at 100
Nobs.push_back(21); // JZB cut at 150
Nobs.push_back(5); // JZB cut at 200
Nobs.push_back(3); // JZB cut at 250
Npred.push_back(407.667); // JZB cut at 50
Npred.push_back(89.3333); // JZB cut at 100
Npred.push_back(22); // JZB cut at 150
Npred.push_back(8); // JZB cut at 200
Npred.push_back(2); // JZB cut at 250
Nprederr.push_back(61.1314); // JZB cut at 50
Nprederr.push_back(13.6275); // JZB cut at 100
Nprederr.push_back(4.25245); // JZB cut at 150
Nprederr.push_back(2.10159); // JZB cut at 200
Nprederr.push_back(0.874007); // JZB cut at 250


flippedNobs.push_back(149); // JZB cut at 50
flippedNobs.push_back(11); // JZB cut at 100
flippedNobs.push_back(3); // JZB cut at 150
flippedNobs.push_back(1); // JZB cut at 200
flippedNobs.push_back(0); // JZB cut at 250
flippedNpred.push_back(149.333); // JZB cut at 50
flippedNpred.push_back(9.66667); // JZB cut at 100
flippedNpred.push_back(2); // JZB cut at 150
flippedNpred.push_back(-2); // JZB cut at 200
flippedNpred.push_back(1); // JZB cut at 250
flippedNprederr.push_back(114.072); // JZB cut at 50
flippedNprederr.push_back(27.1872); // JZB cut at 100
flippedNprederr.push_back(7.9425); // JZB cut at 150
flippedNprederr.push_back(3.14907); // JZB cut at 200
flippedNprederr.push_back(2.08); // JZB cut at 250


luminosity=4980;


cout << "Configuration successfully loaded!" << endl; 
 
 } 
 
Revision:	1.17
Committed:	Wed Mar 21 16:13:09 2012 UTC (13 years, 1 month ago) by buchmann
Content type:	text/plain
Branch:	MAIN
CVS Tags:	cbaf_4_98ifb_paper, HEAD
Changes since 1.16:	+1 -1 lines
Log Message:	Updated luminosity and luminosity uncertainty (as per HN phys announcement 1531)
#	User	Rev	Content
1	buchmann	1.1	#include <iostream>
2			#include <vector>
3	buchmann	1.11	#ifndef SampleClassLoaded
4			#include "SampleClass.C"
5			#endif
6			#define SetupLoaded
7			#ifndef ResultLibraryClassLoaded
8			#include "ResultLibraryClass.C"
9			#endif
10	buchmann	1.1
11			using namespace std;
12
13			namespace PlottingSetup {
14			string datajzb="datajzb_ERROR";
15			string mcjzb="mcjzb_ERROR";
16			vector<float>jzb_cuts;
17			float MCPeakError=-999;
18			}
19
20			void read_config() {
21	buchmann	1.16	datajzb="(jzb[1]+0.06*pt-3.58612)";
22			mcjzb="(jzb[1]+0.04*pt-2.32211)";
23	buchmann	1.1
24	buchmann	1.2
25
26	buchmann	1.16	MCPeakError=0.0990011;
27	buchmann	1.1
28	buchmann	1.10	jzb_cuts.push_back(50); // JZB cut at 50
29			jzb_cuts.push_back(100); // JZB cut at 100
30			jzb_cuts.push_back(150); // JZB cut at 150
31			jzb_cuts.push_back(200); // JZB cut at 200
32			jzb_cuts.push_back(250); // JZB cut at 250
33
34
35	buchmann	1.13	Nobs.push_back(408); // JZB cut at 50
36			Nobs.push_back(88); // JZB cut at 100
37	buchmann	1.12	Nobs.push_back(21); // JZB cut at 150
38	buchmann	1.11	Nobs.push_back(5); // JZB cut at 200
39			Nobs.push_back(3); // JZB cut at 250
40	buchmann	1.13	Npred.push_back(407.667); // JZB cut at 50
41			Npred.push_back(89.3333); // JZB cut at 100
42	buchmann	1.12	Npred.push_back(22); // JZB cut at 150
43			Npred.push_back(8); // JZB cut at 200
44			Npred.push_back(2); // JZB cut at 250
45	buchmann	1.13	Nprederr.push_back(61.1314); // JZB cut at 50
46			Nprederr.push_back(13.6275); // JZB cut at 100
47	buchmann	1.12	Nprederr.push_back(4.25245); // JZB cut at 150
48			Nprederr.push_back(2.10159); // JZB cut at 200
49			Nprederr.push_back(0.874007); // JZB cut at 250
50
51
52	buchmann	1.13	flippedNobs.push_back(149); // JZB cut at 50
53	buchmann	1.12	flippedNobs.push_back(11); // JZB cut at 100
54			flippedNobs.push_back(3); // JZB cut at 150
55			flippedNobs.push_back(1); // JZB cut at 200
56	buchmann	1.10	flippedNobs.push_back(0); // JZB cut at 250
57	buchmann	1.13	flippedNpred.push_back(149.333); // JZB cut at 50
58			flippedNpred.push_back(9.66667); // JZB cut at 100
59	buchmann	1.12	flippedNpred.push_back(2); // JZB cut at 150
60			flippedNpred.push_back(-2); // JZB cut at 200
61			flippedNpred.push_back(1); // JZB cut at 250
62	buchmann	1.13	flippedNprederr.push_back(114.072); // JZB cut at 50
63			flippedNprederr.push_back(27.1872); // JZB cut at 100
64	buchmann	1.12	flippedNprederr.push_back(7.9425); // JZB cut at 150
65			flippedNprederr.push_back(3.14907); // JZB cut at 200
66			flippedNprederr.push_back(2.08); // JZB cut at 250
67	buchmann	1.2
68
69	buchmann	1.17	luminosity=4980;
70	buchmann	1.4
71	buchmann	1.1
72			cout << "Configuration successfully loaded!" << endl;
73
74			}
75