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.2 by lantonel, Tue Apr 30 22:22:08 2013 UTC vs.
Revision 1.6 by lantonel, Mon Aug 19 15:51:14 2013 UTC

# Line 82 | 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";
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 100 | 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 115 | 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 130 | 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 145 | 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 160 | 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";
247 >          //cout << "mvaNonTrig\n";
248            if (fabs (eta) < 0.8)
249              {
250                if (pt < 20.0)
# Line 177 | Line 255 | ElectronSFWeight::at (const double &eta,
255                  return 0.974;
256                if (pt < 50.0 && pt >= 40.0)
257                  return 0.988;
258 <              if (pt < 200.0 && pt >= 50.0)
258 >              if (pt >= 50.0)
259                  return 0.999;
260              }
261            if (fabs (eta) < 1.478 && fabs (eta) >= 0.8)
# Line 190 | Line 268 | ElectronSFWeight::at (const double &eta,
268                  return 0.967;
269                if (pt < 50.0 && pt >= 40.0)
270                  return 0.984;
271 <              if (pt < 200.0 && pt >= 50.0)
271 >              if (pt >= 50.0)
272                  return 0.998;
273              }
274            if (fabs (eta) <= 2.5 && fabs (eta) >= 1.478)
# Line 203 | Line 281 | ElectronSFWeight::at (const double &eta,
281                  return 0.932;
282                if (pt < 50.0 && pt >= 40.0)
283                  return 0.968;
284 <              if (pt < 200.0 && pt >= 50.0)
284 >              if (pt >= 50.0)
285                  return 0.983;
286              }
287          }

Diff Legend

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