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

# User Rev Content
1 auterman 1.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 auterman 1.2 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 auterman 1.1
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 auterman 1.2 void ReadGeneratedXsects(std::string file);
179 auterman 1.1
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 auterman 1.2 std::vector<SUSY_XSECS> xsects;
184 auterman 1.1
185     };
186    
187    
188     #endif