| 1 |
auterman |
1.1 |
#include "SusyScan.h"
|
| 2 |
|
|
#include "GeneratorMasses.h"
|
| 3 |
|
|
|
| 4 |
|
|
#include <cmath>
|
| 5 |
|
|
|
| 6 |
|
|
|
| 7 |
auterman |
1.3 |
double Luminosity = 36.3; //[pb^-1]
|
| 8 |
auterman |
1.1 |
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 |
auterman |
1.3 |
double ExpExclusion(const SusyScan* p){ return (ExpXsecLimit(p)<Xsection(p)&&ExpXsecLimit(p)>0.01?0.01:1); }
|
| 23 |
|
|
double ObsExclusion(const SusyScan* p){ return (ObsXsecLimit(p)<Xsection(p)&&ObsXsecLimit(p)>0.01?0.01:1); }
|
| 24 |
auterman |
1.1 |
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 |
auterman |
1.3 |
double XsecOverExpected(const SusyScan* p){ return (ExpXsecLimit(p)==0 ? 9999. : Xsection(p)/ExpXsecLimit(p)); }
|
| 27 |
auterman |
1.1 |
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.3 |
double PLExpNSignLimit(const SusyScan* p){ return p->PLExpNsigLimit; }
|
| 31 |
|
|
double PLObsNSignLimit(const SusyScan* p){ return p->PLExpNsigLimit; }
|
| 32 |
|
|
double PLExpXsecLimit(const SusyScan* p){ return p->PLExpXsecLimit; }
|
| 33 |
|
|
double PLObsXsecLimit(const SusyScan* p){ return p->PLObsXsecLimit; }
|
| 34 |
|
|
double PLExpExclusion(const SusyScan* p){ return (PLExpXsecLimit(p)<Xsection(p)&&PLExpXsecLimit(p)>0.01?0.01:1); }
|
| 35 |
|
|
double PLObsExclusion(const SusyScan* p){ return (PLObsXsecLimit(p)<Xsection(p)&&PLObsXsecLimit(p)>0.01?0.01:1); }
|
| 36 |
auterman |
1.1 |
|
| 37 |
|
|
|
| 38 |
|
|
|
| 39 |
|
|
double Mzero(const GeneratorMasses* p){ return p->Mzero; }
|
| 40 |
|
|
double Mhalf(const GeneratorMasses* p){ return p->Mhalf; }
|
| 41 |
|
|
double MGluino(const GeneratorMasses* p){ return p->MGL; }
|
| 42 |
|
|
double MSquarkL(const GeneratorMasses* p){ return p->MUL; }
|
| 43 |
|
|
double MSquarkR(const GeneratorMasses* p){ return p->MUR; }
|
| 44 |
|
|
double MChi1(const GeneratorMasses* p){ return p->MZ1; }
|
| 45 |
|
|
double MChi2(const GeneratorMasses* p){ return p->MZ2; }
|
| 46 |
|
|
double MChi3(const GeneratorMasses* p){ return p->MZ3; }
|
| 47 |
|
|
double MChi4(const GeneratorMasses* p){ return p->MZ4; }
|
| 48 |
|
|
double MCha1(const GeneratorMasses* p){ return p->MW1; }
|
| 49 |
|
|
double MCha2(const GeneratorMasses* p){ return p->MW2; }
|
