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

# User Rev Content
1 auterman 1.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 auterman 1.2 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 auterman 1.1
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; }