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.5 by auterman, Wed Jan 26 15:23:05 2011 UTC vs.
Revision 1.8 by auterman, Wed Feb 9 08:51:32 2011 UTC

# Line 6 | Line 6
6   #include "TSpline.h"
7  
8   #include <cmath>
9 + #include <iostream>
10  
11  
12   double Luminosity = 36.3; //[pb^-1]
# Line 28 | Line 29 | double SignalUncertainty(const SusyScan*
29   double SignalRelUncertainty(const SusyScan* p){ return p->signal_uncertainty/p->signal; }
30   double ExpExclusion(const SusyScan* p){ return (ExpXsecLimit(p)<Xsection(p)&&ExpXsecLimit(p)>0.01?1:0.01); }
31   double ObsExclusion(const SusyScan* p){ return (ObsXsecLimit(p)<Xsection(p)&&ObsXsecLimit(p)>0.01?1:0.01); }
32 + double ExpExclCL(const SusyScan* p){ return (p->CLs_b_xsec<=0.05 ? 1:0.01); }
33 + double ExpExclCLm2sigma(const SusyScan* p){ return (p->CLs_b_n2_xsec<=0.05 ? 1:0.01); }
34 + double ExpExclCLm1sigma(const SusyScan* p){ return (p->CLs_b_n1_xsec<=0.05 ? 1:0.01); }
35 + double ExpExclCLp1sigma(const SusyScan* p){ return (p->CLs_b_p1_xsec<=0.05 ? 1:0.01); }
36 + double ExpExclCLp2sigma(const SusyScan* p){ return (p->CLs_b_p2_xsec<=0.05 ? 1:0.01); }
37 + double ObsExclCL(const SusyScan* p){ return (p->CLs_xsec<=0.05 ? 1:0.01); }
38   double SoverSqrtB(const SusyScan* p){ return p->signal/(sqrt(p->background)+p->background_uncertainty+p->signal_uncertainty); }
39   double XsecOverObserved(const SusyScan* p){ return (ObsXsecLimit(p)==0 ? 9999. : Xsection(p)/ObsXsecLimit(p)); }
40   double XsecOverExpected(const SusyScan* p){ return (ExpXsecLimit(p)==0 ? 9999. : Xsection(p)/ExpXsecLimit(p)); }
# Line 52 | Line 59 | double MCMCExpXsecLimit(const SusyScan*
59   double MCMCObsXsecLimit(const SusyScan* p){ return p->MCMCObsXsecLimit; }
60   double MCMCExpExclusion(const SusyScan* p){ return (MCMCExpXsecLimit(p)<Xsection(p)&&MCMCExpXsecLimit(p)>0.01?1:0.01); }
61   double MCMCObsExclusion(const SusyScan* p){ return (MCMCObsXsecLimit(p)<Xsection(p)&&MCMCObsXsecLimit(p)>0.01?1:0.01); }
62 + double SignalContamination(const SusyScan* p){return p->signal_contamination; }
63  
64   double Mzero(const GeneratorMasses* p){ return p->Mzero; }
65   double Mhalf(const GeneratorMasses* p){ return p->Mhalf; }

Diff Legend

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