ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/OSNote2010/plotTrgEff/triggerSuperModel.cc
(Generate patch)

Comparing UserCode/claudioc/OSNote2010/plotTrgEff/triggerSuperModel.cc (file contents):
Revision 1.1 by claudioc, Tue Nov 9 04:39:45 2010 UTC vs.
Revision 1.2 by claudioc, Tue Nov 9 13:31:13 2010 UTC

# Line 21 | Line 21 | float triggerSuperModelEffic(int hyp, fl
21    // Plateau efficiency of muon trigger
22    // (For now we ignore the fact that in the mu9 period
23    //  it was not quite as good)
24 <  float effmu = 0.93;  
24 >  float effmu   = 0.93;  // eta<2.1
25 >  float effmu24 = 0.50;  // 2.4<eta<2.4
26  
27    // Fraction of luminosity where mu9 was not prescaled
28    // run <= 147116
# Line 30 | Line 31 | float triggerSuperModelEffic(int hyp, fl
31    // Fraction of luminosity where mu9 was prescaled and mu11
32    // was unprescaled
33    // 147196 <= run <= 148058
34 <  float f11=0.04;
34 >  float f11=0.273;
35  
36    // Fraction of luminosity where mu9 and mu11 were prescaled
37    // and mu15 was not
# Line 121 | Line 122 | float triggerSuperModelEffic(int hyp, fl
122  
123      // both at high eta
124      if (eta1>2.1 && eta2>2.1) {
125 <      float eff = effmu*effmu;
125 >      
126 >      float eff = effmu*effmu; // this neglects the trigger at high eta
127 >      if (pt2>=15) {
128 >        eff = eff + 2*(1.-effmu)*effmu24;
129 >      } else if (pt2>=11) {
130 >        eff = eff + (1.+f9+f11)*(1-effmu)*effmu24;
131 >      } else {
132 >        eff = eff + (1.+f9)*(1-effmu)*effmu24;
133 >      }
134        return eff;
135      }
136  
137      // One with pt>15 eta<2.1.  The other with eta>2.1
138      if ( (pt1>=15 && eta1<=2.1 && eta2>2.1) ||
139           (pt2>=15 && eta2<=2.1 && eta1>2.1)    ) {
140 <      float eff = effmu;
140 >      float eff = effmu;  // this neglects the trigger at high eta
141 >      if ( (pt1>=15 && eta1>2.1) || (pt2>=15 && eta2>2.1) ) {
142 >        eff = eff + (1-effmu)*effmu24;
143 >      } else if ( (pt1>=11 && eta1>2.1) || (pt2>=11 && eta2>2.1) ) {
144 >        eff = eff + (f9+f11)*(1-effmu)*effmu24;
145 >      } else {
146 >        eff = eff + f9*(1-effmu)*effmu24;
147 >      }
148        return eff;
149      }
150      
# Line 136 | Line 152 | float triggerSuperModelEffic(int hyp, fl
152      // First with 11<pt<15 eta<2.1.  Second one with eta>2.1
153      if ( (pt1>=11 && pt1<15 && eta1<=2.1 && eta2>2.1) ||
154           (pt2>=11 && pt2<15 && eta2<=2.1 && eta1>2.1)   ) {
155 <      float eff = (f9+f11)*effmu + (1-f9-f11)*effmu*effmu;
155 >      float eff = (f9+f11)*effmu + (1-f9-f11)*effmu*effmu;
156 >      eff = eff + (1-effmu)*effmu24; // allow trigger at high eta
157        return eff;
158      }
159  
# Line 144 | Line 161 | float triggerSuperModelEffic(int hyp, fl
161      if ( (pt1<11 && eta1<=2.1 && eta2>2.1) ||
162           (pt2<11 && eta2<=2.1 && eta1>2.1)   ) {
163        float eff = f9*effmu + (1-f9)*effmu*effmu;
164 +      eff = eff + (1-effmu)*effmu24;  // allow trigger at high eta
165        return eff;
166      }
167  
# Line 230 | Line 248 | float triggerSuperModelEffic(int hyp, fl
248        } else {
249          eff = e10 + e15 + e17 + e17b*eff17b + emess*eff32andUp;
250        }
251 +      // now allow for the possibility of a muon trigger
252 +      if (ptmu>=15) {
253 +        eff = eff + (1-eff)*effmu24;
254 +      } else if (ptmu>=11) {
255 +        eff = eff + (f9+f11)*(1-eff)*effmu24;
256 +      } else {
257 +        eff = eff + f9*(1-eff)*effmu24;
258 +      }
259        return eff;
260      }
261  

Diff Legend

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