ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/GlobalFunctions.h
(Generate patch)

Comparing UserCode/auterman/SusyScan/PlotScript/GlobalFunctions.h (file contents):
Revision 1.2 by auterman, Thu Dec 23 12:57:53 2010 UTC vs.
Revision 1.3 by auterman, Fri Jan 7 08:54:47 2011 UTC

# Line 4 | Line 4
4   #include <cmath>
5  
6  
7 < double Luminosity = 31.2; //[pb^-1]
7 > double Luminosity = 36.3; //[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; }
# Line 19 | Line 19 | double MCha2(const SusyScan* p){ return
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); }
22 > 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   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)); }
26 > double XsecOverExpected(const SusyScan* p){ return (ExpXsecLimit(p)==0 ? 9999. : Xsection(p)/ExpXsecLimit(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 <
30 > 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  
37  
38  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines