ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/scan.h
Revision: 1.2
Committed: Tue Dec 4 12:57:34 2007 UTC (17 years, 5 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.1: +102 -0 lines
Log Message: 
*** empty log message ***

File Contents

# Content
1 #ifndef SCAN_H
2 #define SCAN_H
3
4 #include <vector>
5
6 #include <TH2.h>
7 #include <TGraph.h>
8 #include <TGraph2D.h>
9
10 class MyTGraph : public TGraph{
11 public:
12 MyTGraph(int b):TGraph(b){};
13 void Add(double x,double y){ this->SetPoint( this->GetN(),x,y); };
14 void ls(){
15 double x,y;
16 for (int i=0; i<this->GetN(); ++i){
17 this->GetPoint(i,x,y);
18 std::cout << i << "th: "<<x<<", " <<y<<std::endl;
19 }};
20 };
21
22 class SUSY_POINT {
23 public:
24 double MZERO;
25 double MHALF;
26 double TANB;
27 double SGNMU;
28 double AZERO;
29 double MTOP;
30 double muQ;
31 double Q;
32 double M1;
33 double M2;
34 double M3;
35 double MGL;
36 double MUL;
37 double MB1;
38 double MSN;
39 double MNTAU;
40 double MZ1;
41 double MW1;
42 double MHL;
43 double MUR;
44 double MB2;
45 double MEL;
46 double MTAU1;
47 double MZ2;
48 double MW2;
49 double MHH;
50 double MDL;
51 double MT1;
52 double MER;
53 double MTAU2;
54 double MZ3;
55 double MHA;
56 double MDR;
57 double MT2;
58 double MZ4;
59 double MHp;
60 };
61
62 class SUSY_XSECS {
63 public:
64 double MZERO;
65 double MHALF;
66 double TANB;
67 double SGNMU;
68 double AZERO;
69 double MTOP;
70 double XS0;
71 double XS201;
72 double XS202;
73 double XS204;
74 double XS205;
75 double XS207;
76 double XS208;
77 double XS209;
78 double XS210;
79 double XS211;
80 double XS212;
81 double XS213;
82 double XS214;
83 double XS216;
84 double XS217;
85 double XS218;
86 double XS219;
87 double XS220;
88 double XS221;
89 double XS222;
90 double XS223;
91 double XS224;
92 double XS225;
93 double XS226;
94 double XS227;
95 double XS228;
96 double XS229;
97 double XS230;
98 double XS231;
99 double XS232;
100 double XS233;
101 double XS234;
102 double XS235;
103 double XS236;
104 double XS237;
105 double XS238;
106 double XS239;
107 double XS240;
108 double XS241;
109 double XS242;
110 double XS243;
111 double XS244;
112 double XS246;
113 double XS247;
114 double XS248;
115 double XS249;
116 double XS250;
117 double XS251;
118 double XS252;
119 double XS253;
120 double XS254;
121 double XS256;
122 double XS258;
123 double XS259;
124 double XS261;
125 double XS262;
126 double XS263;
127 double XS264;
128 double XS265;
129 double XS271;
130 double XS272;
131 double XS273;
132 double XS274;
133 double XS275;
134 double XS276;
135 double XS277;
136 double XS278;
137 double XS279;
138 double XS280;
139 double XS281;
140 double XS282;
141 double XS283;
142 double XS284;
143 double XS285;
144 double XS286;
145 double XS287;
146 double XS288;
147 double XS289;
148 double XS290;
149 double XS291;
150 double XS292;
151 double XS293;
152 double XS294;
153 double XS295;
154 double XS296;
155 double XS297;
156 double XS298;
157 double XS299;
158 double XS300;
159 double XS301;
160 };
161
162
163 class TScan {
164 public:
165 TScan();
166 ~TScan();
167
168 MyTGraph * IsoMassLine(bool(*func)(SUSY_POINT,double), double);
169 MyTGraph * GetContour(TH2F*, int flag=0);
170 TH2F * Area(bool(*func)(SUSY_POINT));
171 TH2F * Area51();
172 //static bool ul(const SUSY_POINT point, const double mass);
173
174 int DoStuff();
175
176 private:
177 void ReadGeneratedMasses(std::string file);
178 void ReadGeneratedXsects(std::string file);
179
180 int plot_id, bins_x, bins_y;
181 double min_x, max_x, min_y, max_y;
182 std::vector<SUSY_POINT> points;
183 std::vector<SUSY_XSECS> xsects;
184
185 };
186
187
188 #endif