ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/GlobalFunctions.h
Revision: 1.2
Committed: Thu Dec 23 12:57:53 2010 UTC (14 years, 4 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.1: +2 -1 lines
Log Message: 
*** empty log message ***

File Contents

# Content
1 #include "SusyScan.h"
2 #include "GeneratorMasses.h"
3
4 #include <cmath>
5
6
7 double Luminosity = 31.2; //[pb^-1]
8 double Mzero(const SusyScan* p){ return p->Mzero; }
9 double Mhalf(const SusyScan* p){ return p->Mhalf; }
10 double MGluino(const SusyScan* p){ return p->MGL; }
11 double MSquarkL(const SusyScan* p){ return p->MUL; }
12 double MSquarkR(const SusyScan* p){ return p->MUR; }
13 double MChi1(const SusyScan* p){ return p->MZ1; }
14 double MChi2(const SusyScan* p){ return p->MZ2; }
15 double MChi3(const SusyScan* p){ return p->MZ3; }
16 double MChi4(const SusyScan* p){ return p->MZ4; }
17 double MCha1(const SusyScan* p){ return p->MW1; }
18 double MCha2(const SusyScan* p){ return p->MW2; }
19 double Xsection(const SusyScan* p){ return p->Xsection; }
20 double ExpXsecLimit(const SusyScan* p){ return p->ExpXsecLimit; }
21 double ObsXsecLimit(const SusyScan* p){ return p->ObsXsecLimit; }
22 double ExpExclusion(const SusyScan* p){ return (ExpXsecLimit(p)<Xsection(p)&&ExpXsecLimit(p)>0.01?1:-1); }
23 double ObsExclusion(const SusyScan* p){ return (ObsXsecLimit(p)<Xsection(p)&&ObsXsecLimit(p)>0.01?1:-1); }
24 double SoverSqrtB(const SusyScan* p){ return p->signal/(sqrt(p->background)+p->background_uncertainty+p->signal_uncertainty); }
25 double XsecOverObserved(const SusyScan* p){ return (ObsXsecLimit(p)==0 ? 9999. : Xsection(p)/ObsXsecLimit(p)); }
26 double XsecOverExpected(const SusyScan* p){ return (ObsXsecLimit(p)==0 ? 9999. : Xsection(p)/ObsXsecLimit(p)); }
27 double SignalAcceptance(const SusyScan* p){ return p->signal / (Luminosity*Xsection(p)); }
28 double ExpNSignLimit(const SusyScan* p){ return p->signal * ExpXsecLimit(p)/Xsection(p); }
29 double ObsNSignLimit(const SusyScan* p){ return p->signal * ObsXsecLimit(p)/Xsection(p); }
30
31
32
33
34
35 double Mzero(const GeneratorMasses* p){ return p->Mzero; }
36 double Mhalf(const GeneratorMasses* p){ return p->Mhalf; }
37 double MGluino(const GeneratorMasses* p){ return p->MGL; }
38 double MSquarkL(const GeneratorMasses* p){ return p->MUL; }
39 double MSquarkR(const GeneratorMasses* p){ return p->MUR; }
40 double MChi1(const GeneratorMasses* p){ return p->MZ1; }
41 double MChi2(const GeneratorMasses* p){ return p->MZ2; }
42 double MChi3(const GeneratorMasses* p){ return p->MZ3; }
43 double MChi4(const GeneratorMasses* p){ return p->MZ4; }
44 double MCha1(const GeneratorMasses* p){ return p->MW1; }
45 double MCha2(const GeneratorMasses* p){ return p->MW2; }