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.10 by auterman, Wed Feb 9 21:00:29 2011 UTC vs.
Revision 1.12 by auterman, Wed Mar 2 13:31:32 2011 UTC

# Line 25 | Line 25 | double MChi3(const SusyScan* p){ return
25   double MChi4(const SusyScan* p){ return p->MZ4; }
26   double MCha1(const SusyScan* p){ return p->MW1; }
27   double MCha2(const SusyScan* p){ return p->MW2; }
28 + double SignalUncertKfactor(const SusyScan* p){return fabs(p->signal_kfactor_UP-p->signal_kfactor_DN)/(2.0*p->signal_kfactor); }
29 + double SignalUncertJEC(const SusyScan* p){ return (fabs(p->signal_JEC_UP)+fabs(p->signal_JEC_DN))/(2.0*p->signal); }
30 + double SignalUncertMuIso(const SusyScan* p){ return (fabs(p->signal_MuIso_UP)+fabs(p->signal_MuIso_DN))/(2.0*p->signal); }
31 + double SignalKfactor(const SusyScan* p){return p->signal_kfactor; }
32 +
33   double Xsection(const SusyScan* p){ return p->Xsection; }
34   double ExpXsecLimit(const SusyScan* p){ return p->ExpXsecLimit; }
35   double ObsXsecLimit(const SusyScan* p){ return p->ObsXsecLimit; }
# Line 65 | Line 70 | double MCMCExpExclusion(const SusyScan*
70   double MCMCObsExclusion(const SusyScan* p){ return (MCMCObsXsecLimit(p)<Xsection(p)&&MCMCObsXsecLimit(p)>0.01?1:0.01); }
71   double SignalContamination(const SusyScan* p){return p->signal_contamination; }
72  
73 + double NLOXsection(const SusyScan* p){ return p->Xsection * p->signal_kfactor; }
74 + double NLOExpXsecLimit(const SusyScan* p){ return p->NLO_ExpXsecLimit; }
75 + double NLOObsXsecLimit(const SusyScan* p){ return p->NLO_ObsXsecLimit; }
76 + double NLOSignal(const SusyScan* p){ return p->NLO_signal; }
77 + double NLOSignalUncertainty(const SusyScan* p){ return p->NLO_signal_uncertainty; }
78 + double NLOSignalRelUncertainty(const SusyScan* p){ return p->NLO_signal_uncertainty/p->NLO_signal; }
79 + double NLOExpExclusion(const SusyScan* p){ return (NLOExpXsecLimit(p)<NLOXsection(p)&&NLOExpXsecLimit(p)>0.01?1:0.01); }
80 + double NLOObsExclusion(const SusyScan* p){ return (NLOObsXsecLimit(p)<NLOXsection(p)&&NLOObsXsecLimit(p)>0.01?1:0.01); }
81 + double NLOExpExclCL(const SusyScan* p){ return (p->NLO_CLs_b_xsec<=0.05 ? 1:0.01); }
82 + double NLOExpExclCLm2sigma(const SusyScan* p){ return (p->NLO_CLs_b_n2_xsec<=0.05 ? 1:0.01); }
83 + double NLOExpExclCLm1sigma(const SusyScan* p){ return (p->NLO_CLs_b_n1_xsec<=0.05 ? 1:0.01); }
84 + double NLOExpExclCLp1sigma(const SusyScan* p){ return (p->NLO_CLs_b_p1_xsec<=0.05 ? 1:0.01); }
85 + double NLOExpExclCLp2sigma(const SusyScan* p){ return (p->NLO_CLs_b_p2_xsec<=0.05 ? 1:0.01); }
86 + double NLOObsExclCL(const SusyScan* p){ return (p->NLO_CLs_xsec<=0.05 ? 1:0.01); }
87 + double NLOSoverSqrtB(const SusyScan* p){ return p->NLO_signal/(sqrt(p->background)+p->background_uncertainty+p->NLO_signal_uncertainty); }
88 + double NLOXsecOverObserved(const SusyScan* p){ return (NLOObsXsecLimit(p)==0 ? 9999. : NLOXsection(p)/NLOObsXsecLimit(p)); }
89 + double NLOXsecOverExpected(const SusyScan* p){ return (NLOExpXsecLimit(p)==0 ? 9999. : NLOXsection(p)/NLOExpXsecLimit(p)); }
90 + double NLOSignalAcceptance(const SusyScan* p){ return  p->NLO_signal / (Luminosity*NLOXsection(p)); }
91 + double NLOExpNSignLimit(const SusyScan* p){ return  p->NLO_ExpNsigLimit; }
92 + double NLOObsNSignLimit(const SusyScan* p){ return  p->NLO_ObsNsigLimit; }
93 + double NLOPLExpNSignLimit(const SusyScan* p){ return  p->NLO_PLExpNsigLimit; }
94 + double NLOPLObsNSignLimit(const SusyScan* p){ return  p->NLO_PLObsNsigLimit; }
95 + double NLOPLExpXsecLimit(const SusyScan* p){ return p->NLO_PLExpXsecLimit; }
96 + double NLOPLObsXsecLimit(const SusyScan* p){ return p->NLO_PLObsXsecLimit; }
97 + double NLOPLExpExclusion(const SusyScan* p){ return (NLOPLExpXsecLimit(p)<NLOXsection(p)&&NLOPLExpXsecLimit(p)>0.01?1:0.01); }
98 + double NLOPLObsExclusion(const SusyScan* p){ return (NLOPLObsXsecLimit(p)<NLOXsection(p)&&NLOPLObsXsecLimit(p)>0.01?1:0.01); }
99 + double NLOFCExpNSignLimit(const SusyScan* p){ return  p->NLO_FCExpNsigLimit; }
100 + double NLOFCObsNSignLimit(const SusyScan* p){ return  p->NLO_FCObsNsigLimit; }
101 + double NLOFCExpXsecLimit(const SusyScan* p){ return p->NLO_FCExpXsecLimit; }
102 + double NLOFCObsXsecLimit(const SusyScan* p){ return p->NLO_FCObsXsecLimit; }
103 + double NLOFCExpExclusion(const SusyScan* p){ return (NLOFCExpXsecLimit(p)<NLOXsection(p)&&NLOFCExpXsecLimit(p)>0.01?1:0.01); }
104 + double NLOFCObsExclusion(const SusyScan* p){ return (NLOFCObsXsecLimit(p)<NLOXsection(p)&&NLOFCObsXsecLimit(p)>0.01?1:0.01); }
105 + double NLOMCMCExpNSignLimit(const SusyScan* p){ return  p->NLO_MCMCExpNsigLimit; }
106 + double NLOMCMCObsNSignLimit(const SusyScan* p){ return  p->NLO_MCMCObsNsigLimit; }
107 + double NLOMCMCExpXsecLimit(const SusyScan* p){ return p->NLO_MCMCExpXsecLimit; }
108 + double NLOMCMCObsXsecLimit(const SusyScan* p){ return p->NLO_MCMCObsXsecLimit; }
109 + double NLOMCMCExpExclusion(const SusyScan* p){ return (NLOMCMCExpXsecLimit(p)<NLOXsection(p)&&NLOMCMCExpXsecLimit(p)>0.01?1:0.01); }
110 + double NLOMCMCObsExclusion(const SusyScan* p){ return (NLOMCMCObsXsecLimit(p)<NLOXsection(p)&&NLOMCMCObsXsecLimit(p)>0.01?1:0.01); }
111 + double NLOSignalContamination(const SusyScan* p){return p->NLO_signal_contamination; }
112 +
113   double Mzero(const GeneratorMasses* p){ return p->Mzero; }
114   double Mhalf(const GeneratorMasses* p){ return p->Mhalf; }
115   double MGluino(const GeneratorMasses* p){ return p->MGL; }
# Line 104 | Line 149 | TGraph* set_sneutrino_d0_2(Int_t tanBeta
149  
150   TGraph* set_lep_ch_tanBeta10(){
151  
152 <  double ch_m0[11];
153 <  double ch_m12[11];
152 >  double ch_m0[12];
153 >  double ch_m12[12];
154  
155    ch_m0[0] = 0;
156    ch_m0[1] = 100;
# Line 116 | Line 161 | TGraph* set_lep_ch_tanBeta10(){
161    ch_m0[6] = 600;
162    ch_m0[7] = 700;
163    ch_m0[8] = 800;
164 <  ch_m0[9] = 800;
165 <  ch_m0[10] = 0;
164 >  ch_m0[9] = 1000;
165 >  ch_m0[10] = 1000;
166 >  ch_m0[11] = 0;
167  
168    ch_m12[0] = 163;
169    ch_m12[1] = 162;
# Line 128 | Line 174 | TGraph* set_lep_ch_tanBeta10(){
174    ch_m12[6] = 157;
175    ch_m12[7] = 156;
176    ch_m12[8] = 155.4;
177 <  ch_m12[9] = 0;
177 >  ch_m12[9] = 154;
178    ch_m12[10] = 0;
179 +  ch_m12[11] = 0;
180    
181    
182 <  TGraph* ch_gr = new TGraph(11,ch_m0,ch_m12);
182 >  TGraph* ch_gr = new TGraph(12,ch_m0,ch_m12);
183  
184    ch_gr->SetFillColor(3);
185    ch_gr->SetLineColor(3);
# Line 414 | Line 461 | TGraph* set_tev_stau(Int_t tanBeta){
461  
462   }
463  
464 + TGraph* sq_LEP(){
465 +    double sq[] = {0,0,100,100};
466 +    double gl[] = {0,2000,2000,0};  
467 +    TGraph* res = new TGraph(4,sq,gl);
468 +    res->SetFillColor(kBlue);
469 +    return res;
470 + }
471 +
472 + TGraph* gl_TEV(){
473 +    double sq[] = {0,2000,2000,0};
474 +    double gl[] = {0,0,190,190};  
475 +    TGraph* res = new TGraph(4,sq,gl);
476 +    res->SetFillColor(kGreen+2);
477 +    return res;
478 + }
479 +
480 + TGraph* gl_CDF(){
481 +    double sq[] = {0,2000,2000,0};
482 +    double gl[] = {190,190,230,230};  
483 +    TGraph* res = new TGraph(4,sq,gl);
484 +    res->SetFillColor(kOrange+5);
485 +    return res;
486 + }
487 +
488 + TGraph* gl_DEZ(){
489 +    double sq[] = {0,2000,2000,0};
490 +    double gl[] = {230,230,255,255};  
491 +    TGraph* res = new TGraph(4,sq,gl);
492 +    res->SetFillColor(kYellow-5);
493 +    return res;
494 + }
495 +
496 + TGraph* gl_WHT(){
497 +    double sq[] = {101,2000,2000,101};
498 +    double gl[] = {256,256,400,400};  
499 +    TGraph* res = new TGraph(4,sq,gl);
500 +    res->SetFillColor(kWhite);
501 +    return res;
502 + }
503  
504  
505  

Diff Legend

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