ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/OSUT3Analysis/AnaTools/src/SFWeight.cc
(Generate patch)

Comparing UserCode/OSUT3Analysis/AnaTools/src/SFWeight.cc (file contents):
Revision 1.1 by ahart, Tue Mar 12 20:43:04 2013 UTC vs.
Revision 1.6 by lantonel, Mon Aug 19 15:51:14 2013 UTC

# Line 2 | Line 2
2  
3   MuonSFWeight::MuonSFWeight (const string &sfFile, const string &dataOverMC)
4   {
5 +
6    TFile *fin = TFile::Open (sfFile.c_str ());
7    TGraphAsymmErrors *data = (TGraphAsymmErrors *) fin->Get (dataOverMC.c_str ());
8  
# Line 81 | Line 82 | ElectronSFWeight::ElectronSFWeight (cons
82   double
83   ElectronSFWeight::at (const double &eta, const double &pt)
84   {
85 +  //cout << eta << " " << pt << endl;
86 +  //cout << "doing it\n";
87    if (cmsswRelease_ == "53X")
88      {
89 +      if (id_ == "loose")
90 +        {
91 +          if (fabs (eta) < 0.8)
92 +            {
93 +              if (pt < 15.0 && pt >= 10.0)
94 +                return 0.855;
95 +              if (pt < 20.0 && pt >= 15.0)
96 +                return 0.962;
97 +              if (pt < 30.0 && pt >= 20.0)
98 +                return 1.005;
99 +              if (pt < 40.0 && pt >= 30.0)
100 +                return 1.004;
101 +              if (pt < 50.0 && pt >= 40.0)
102 +                return 1.008;
103 +              if (pt >= 50.0)
104 +                return 1.008;
105 +            }
106 +          if (fabs (eta) < 1.442 && fabs (eta) >= 0.8)
107 +            {
108 +              if (pt < 15.0 && pt >= 10.0)
109 +                return 0.858;
110 +              if (pt < 20.0 && pt >= 15.0)
111 +                return 0.962;
112 +              if (pt < 30.0 && pt >= 20.0)
113 +                return 0.981;
114 +              if (pt < 40.0 && pt >= 30.0)
115 +                return 0.991;
116 +              if (pt < 50.0 && pt >= 40.0)
117 +                return 0.994;
118 +              if (pt >= 50.0)
119 +                return 0.999;
120 +            }
121 +          if (fabs (eta) < 1.556 && fabs (eta) >= 1.442)
122 +            {
123 +              if (pt < 15.0 && pt >= 10.0)
124 +                return 1.109;
125 +              if (pt < 20.0 && pt >= 15.0)
126 +                return 0.903;
127 +              if (pt < 30.0 && pt >= 20.0)
128 +                return 1.044;
129 +              if (pt < 40.0 && pt >= 30.0)
130 +                return 0.998;
131 +              if (pt < 50.0 && pt >= 40.0)
132 +                return 0.989;
133 +              if (pt >= 50.0)
134 +                return 0.994;
135 +            }
136 +          if (fabs (eta) < 2.0 && fabs (eta) >= 1.556)
137 +            {
138 +              if (pt < 15.0 && pt >= 10.0)
139 +                return 0.838;
140 +              if (pt < 20.0 && pt >= 15.0)
141 +                return 0.939;
142 +              if (pt < 30.0 && pt >= 20.0)
143 +                return 0.980;
144 +              if (pt < 40.0 && pt >= 30.0)
145 +                return 0.992;
146 +              if (pt < 50.0 && pt >= 40.0)
147 +                return 1.004;
148 +              if (pt >= 50.0)
149 +                return 1.006;
150 +            }
151 +          if (fabs (eta) < 2.5 && fabs (eta) >= 2.0)
152 +            {
153 +              if (pt < 15.0 && pt >= 10.0)
154 +                return 1.034;
155 +              if (pt < 20.0 && pt >= 15.0)
156 +                return 0.970;
157 +              if (pt < 30.0 && pt >= 20.0)
158 +                return 1.017;
159 +              if (pt < 40.0 && pt >= 30.0)
160 +                return 1.019;
161 +              if (pt < 50.0 && pt >= 40.0)
162 +                return 1.005;
163 +              if (pt >= 50.0)
164 +                return 1.009;
165 +            }
166 +        }
167        if (id_ == "tight")
168          {
169            if (fabs (eta) < 0.8)
# Line 97 | Line 178 | ElectronSFWeight::at (const double &eta,
178                  return 0.979;
179                if (pt < 50.0 && pt >= 40.0)
180                  return 0.984;
181 <              if (pt < 200.0 && pt >= 50.0)
181 >              if (pt >= 50.0)
182                  return 0.983;
183              }
184            if (fabs (eta) < 1.442 && fabs (eta) >= 0.8)
# Line 112 | Line 193 | ElectronSFWeight::at (const double &eta,
193                  return 0.961;
194                if (pt < 50.0 && pt >= 40.0)
195                  return 0.972;
196 <              if (pt < 200.0 && pt >= 50.0)
196 >              if (pt >= 50.0)
197                  return 0.977;
198              }
199            if (fabs (eta) < 1.556 && fabs (eta) >= 1.442)
# Line 127 | Line 208 | ElectronSFWeight::at (const double &eta,
208                  return 0.983;
209                if (pt < 50.0 && pt >= 40.0)
210                  return 0.957;
211 <              if (pt < 200.0 && pt >= 50.0)
211 >              if (pt >= 50.0)
212                  return 0.978;
213              }
214            if (fabs (eta) < 2.0 && fabs (eta) >= 1.556)
# Line 142 | Line 223 | ElectronSFWeight::at (const double &eta,
223                  return 0.962;
224                if (pt < 50.0 && pt >= 40.0)
225                  return 0.985;
226 <              if (pt < 200.0 && pt >= 50.0)
226 >              if (pt >= 50.0)
227                  return 0.986;
228              }
229            if (fabs (eta) < 2.5 && fabs (eta) >= 2.0)
# Line 157 | Line 238 | ElectronSFWeight::at (const double &eta,
238                  return 1.002;
239                if (pt < 50.0 && pt >= 40.0)
240                  return 0.999;
241 <              if (pt < 200.0 && pt >= 50.0)
241 >              if (pt >= 50.0)
242                  return 0.995;
243              }
244          }
245 +      if (id_ == "mvaNonTrig")
246 +        {
247 +          //cout << "mvaNonTrig\n";
248 +          if (fabs (eta) < 0.8)
249 +            {
250 +              if (pt < 20.0)
251 +                return 1.0;
252 +              if (pt < 30.0 && pt >= 20.0)
253 +                return 0.994;
254 +              if (pt < 40.0 && pt >= 30.0)
255 +                return 0.974;
256 +              if (pt < 50.0 && pt >= 40.0)
257 +                return 0.988;
258 +              if (pt >= 50.0)
259 +                return 0.999;
260 +            }
261 +          if (fabs (eta) < 1.478 && fabs (eta) >= 0.8)
262 +            {
263 +              if (pt < 20.0)
264 +                return 1.0;
265 +              if (pt < 30.0 && pt >= 20.0)
266 +                return 0.971;
267 +              if (pt < 40.0 && pt >= 30.0)
268 +                return 0.967;
269 +              if (pt < 50.0 && pt >= 40.0)
270 +                return 0.984;
271 +              if (pt >= 50.0)
272 +                return 0.998;
273 +            }
274 +          if (fabs (eta) <= 2.5 && fabs (eta) >= 1.478)
275 +            {
276 +              if (pt < 20.0)
277 +                return 1.0;
278 +              if (pt < 30.0 && pt >= 20.0)
279 +                return 0.921;
280 +              if (pt < 40.0 && pt >= 30.0)
281 +                return 0.932;
282 +              if (pt < 50.0 && pt >= 40.0)
283 +                return 0.968;
284 +              if (pt >= 50.0)
285 +                return 0.983;
286 +            }
287 +        }
288      }
289  
290    return 0.0;

Diff Legend

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