ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/ScansICHEP2012/code/Functions.cc
Revision: 1.2
Committed: Mon Jul 2 04:42:05 2012 UTC (12 years, 10 months ago) by claudioc
Content type: text/plain
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +82 -4 lines
Log Message: 
lipstick

File Contents

# User Rev Content
1 claudioc 1.1
2    
3    
4     double observedLimitOnNumberOfEvents( int mySR, float accUncertainty, bool Posteriori ){
5    
6    
7     int SR ;
8 claudioc 1.2 float err,obs,exp, exp_PlusSigma, exp_MinusSigma, myUL, myMinErr = 9999;
9 claudioc 1.1 TString MyFile = "empty";
10    
11    
12     //input files:
13     if( Posteriori )
14     MyFile = "LandsOutput_Posterior.txt";
15     else
16     MyFile = "LandsOutput_Apriori.txt";
17    
18     ifstream fLimits(""+MyFile+"", ifstream::in);
19     if( !fLimits.is_open() ){
20     cerr << "[ERROR] Limit file is null " << endl;
21     return 1;
22     }
23    
24 claudioc 1.2 while( fLimits >> SR >> err >> exp >> obs >> exp_PlusSigma >> exp_MinusSigma ){
25 claudioc 1.1
26     if( SR == mySR ){
27    
28     // To Find the UL for a given accUncert!!
29     float dAncc = TMath::Abs( err - accUncertainty );
30     if( dAncc < myMinErr ){
31     myMinErr = dAncc ;
32     myUL = obs;
33     }
34     }
35     }
36    
37     return myUL;
38     }
39    
40    
41    
42     double expectedLimitOnNumberOfEvents( int mySR, float accUncertainty, bool Posteriori ){
43    
44    
45     int SR ;
46 claudioc 1.2 float err,obs,exp, exp_PlusSigma, exp_MinusSigma, myUL, myMinErr = 9999;
47 claudioc 1.1 TString MyFile = "empty";
48    
49    
50     //input files:
51     if( Posteriori )
52     MyFile = "LandsOutput_Posterior.txt";
53     else
54     MyFile = "LandsOutput_Apriori.txt";
55    
56     ifstream fLimits(""+MyFile+"", ifstream::in);
57     if( !fLimits.is_open() ){
58     cerr << "[ERROR] Limit file is null " << endl;
59     return 1;
60     }
61    
62 claudioc 1.2 while( fLimits >> SR >> err >> exp >> obs >> exp_PlusSigma >> exp_MinusSigma ){
63 claudioc 1.1
64     if( SR == mySR ){
65    
66     // To Find the UL for a given accUncert!!
67     float dAncc = TMath::Abs( err - accUncertainty );
68     if( dAncc < myMinErr ){
69     myMinErr = dAncc ;
70     myUL = exp;
71     }
72     }
73     }
74    
75     return myUL;
76     }
77    
78 claudioc 1.2
79    
80     double expectedMinusOneSigmaLimitOnNumberOfEvents( int mySR, float accUncertainty, bool Posteriori ){
81    
82    
83     int SR ;
84     float err,obs,exp, exp_PlusSigma, exp_MinusSigma, myUL, myMinErr = 9999;
85     TString MyFile = "empty";
86    
87    
88     //input files:
89     if( Posteriori )
90     MyFile = "LandsOutput_Posterior.txt";
91     else
92     MyFile = "LandsOutput_Apriori.txt";
93    
94     ifstream fLimits(""+MyFile+"", ifstream::in);
95     if( !fLimits.is_open() ){
96     cerr << "[ERROR] Limit file is null " << endl;
97     return 1;
98     }
99    
100     while( fLimits >> SR >> err >> exp >> obs >> exp_PlusSigma >> exp_MinusSigma ){
101    
102     if( SR == mySR ){
103    
104     // To Find the UL for a given accUncert!!
105     float dAncc = TMath::Abs( err - accUncertainty );
106     if( dAncc < myMinErr ){
107     myMinErr = dAncc ;
108     myUL = exp_MinusSigma ;
109     }
110     }
111     }
112    
113     return myUL;
114     }
115    
116    
117    
118    
119    
120     double expectedPlusOneSigmaLimitOnNumberOfEvents( int mySR, float accUncertainty, bool Posteriori ){
121    
122    
123     int SR ;
124     float err,obs,exp, exp_PlusSigma, exp_MinusSigma, myUL, myMinErr = 9999;
125     TString MyFile = "empty";
126    
127    
128     //input files:
129     if( Posteriori )
130     MyFile = "LandsOutput_Posterior.txt";
131     else
132     MyFile = "LandsOutput_Apriori.txt";
133    
134     ifstream fLimits(""+MyFile+"", ifstream::in);
135     if( !fLimits.is_open() ){
136     cerr << "[ERROR] Limit file is null " << endl;
137     return 1;
138     }
139    
140     while( fLimits >> SR >> err >> exp >> obs >> exp_PlusSigma >> exp_MinusSigma ){
141    
142     if( SR == mySR ){
143    
144     // To Find the UL for a given accUncert!!
145     float dAncc = TMath::Abs( err - accUncertainty );
146     if( dAncc < myMinErr ){
147     myMinErr = dAncc ;
148     myUL = exp_PlusSigma ;
149     }
150     }
151     }
152    
153     return myUL;
154     }
155