test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
MuonEffectiveArea.h
Go to the documentation of this file.
1 //--------------------------------------------------------------------------------------------------
2 // $Id $
3 //
4 // MuonEffectiveArea
5 //
6 // Helper Class for storing effective areas
7 //
8 // Authors: S.Xie, E. DiMarco
9 //--------------------------------------------------------------------------------------------------
10 
11 
15 #define STANDALONE // <---- this line
16 
17 #ifndef MuonEffectiveArea_H
18 #define MuonEffectiveArea_H
19 
20 #include <math.h>
21 
22 #ifndef STANDALONE
23 #endif
24 
25 using namespace std;
26 
27 class MuonEffectiveArea{
28  public:
29  MuonEffectiveArea();
30  ~MuonEffectiveArea();
31 
32  enum MuonEffectiveAreaType {
33  kMuTrkIso03,
34  kMuEcalIso03,
35  kMuHcalIso03,
36  kMuTrkIso05,
37  kMuEcalIso05,
38  kMuHcalIso05,
39  kMuChargedIso03,
40  kMuGammaIso03,
41  kMuNeutralHadronIso03,
42  kMuGammaAndNeutralHadronIso03,
43  kMuChargedIso04,
44  kMuGammaIso04,
45  kMuNeutralHadronIso04,
46  kMuGammaAndNeutralHadronIso04,
47  kMuGammaIsoDR0p0To0p1,
48  kMuGammaIsoDR0p1To0p2,
49  kMuGammaIsoDR0p2To0p3,
50  kMuGammaIsoDR0p3To0p4,
51  kMuGammaIsoDR0p4To0p5,
52  kMuNeutralHadronIsoDR0p0To0p1,
53  kMuNeutralHadronIsoDR0p1To0p2,
54  kMuNeutralHadronIsoDR0p2To0p3,
55  kMuNeutralHadronIsoDR0p3To0p4,
56  kMuNeutralHadronIsoDR0p4To0p5,
57  kMuGammaIso05,
58  kMuNeutralIso05
59  };
60 
61  enum MuonEffectiveAreaTarget {
62  kMuEANoCorr,
63  kMuEAData2011,
64  kMuEASummer11MC,
65  kMuEAFall11MC,
66  kMuEAData2012
67  };
68 
69  static const Double_t GetMuonEffectiveArea(MuonEffectiveAreaType type, Double_t SCEta,
70  MuonEffectiveAreaTarget EffectiveAreaTarget = kMuEAData2011) {
71 
72  Double_t EffectiveArea = 0;
73 
74 
75  if (EffectiveAreaTarget == kMuEANoCorr) {
76  return 0.0;
77  }
78 
79  //2012 Data Effective Areas
80  else if (EffectiveAreaTarget == kMuEAData2012) {
81  if (type == kMuGammaIso04){
82  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.50419;
83  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.30582;
84  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.19765;
85  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.28723;
86  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.52529;
87  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.48818;
88  }
89  if (type == kMuNeutralHadronIso04){
90  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.16580;
91  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25904;
92  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.24695;
93  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.22021;
94  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.34045;
95  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.21592;
96  }
97  }
98 
99  //2011 Data Effective Areas
100  else if (EffectiveAreaTarget == kMuEAData2011) {
101 
102  if (type == kMuGammaIsoDR0p0To0p1) {
103  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
104  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
105  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.002;
106  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
107  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
108  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.005;
109  }
110  if (type == kMuGammaIsoDR0p1To0p2) {
111  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.011;
112  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
113  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.005;
114  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
115  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.008;
116  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.011;
117  }
118  if (type == kMuGammaIsoDR0p2To0p3) {
119  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.023;
120  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.016;
121  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.010;
122  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.014;
123  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.017;
124  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.021;
125  }
126  if (type == kMuGammaIsoDR0p3To0p4) {
127  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
128  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
129  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
130  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.023;
131  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.028;
132  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.032;
133  }
134  if (type == kMuGammaIsoDR0p4To0p5) {
135  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
136  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.037;
137  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.028;
138  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.033;
139  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.042;
140  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.052;
141  }
142  if (type == kMuNeutralHadronIsoDR0p0To0p1) {
143  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
144  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.001;
145  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.001;
146  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.001;
147  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
148  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.007;
149  }
150  if (type == kMuNeutralHadronIsoDR0p1To0p2) {
151  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
152  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
153  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
154  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
155  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.010;
156  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.014;
157  }
158  if (type == kMuNeutralHadronIsoDR0p2To0p3) {
159  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
160  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.015;
161  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
162  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.017;
163  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.019;
164  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.024;
165  }
166  if (type == kMuNeutralHadronIsoDR0p3To0p4) {
167  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
168  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
169  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.024;
170  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.032;
171  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
172  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.038;
173  }
174  if (type == kMuNeutralHadronIsoDR0p4To0p5) {
175  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
176  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
177  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.033;
178  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
179  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
180  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.114;
181  }
182  if (type == kMuGammaIso05){
183  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.05317;
184  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.03502;
185  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.03689;
186  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.05221;
187  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.06668;
188  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.0744;
189  }
190  if (type == kMuNeutralIso05) {
191  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.06408;
192  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.07557;
193  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.08864;
194  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.11492;
195  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.13784;
196  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.18745;
197  }
198  }
199 
200  //Summer11 MC Effective Areas
201  else if (EffectiveAreaTarget == kMuEASummer11MC) {
202  if (type == kMuGammaIsoDR0p0To0p1) {
203  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.000;
204  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.000;
205  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
206  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
207  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
208  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.006;
209  }
210  if (type == kMuGammaIsoDR0p1To0p2) {
211  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
212  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
213  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
214  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
215  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.019;
216  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.015;
217  }
218  if (type == kMuGammaIsoDR0p2To0p3) {
219  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.023;
220  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.018;
221  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
222  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
223  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
224  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.036;
225  }
226  if (type == kMuGammaIsoDR0p3To0p4) {
227  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.038;
228  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
229  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.019;
230  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.033;
231  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
232  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.062;
233  }
234  if (type == kMuGammaIsoDR0p4To0p5) {
235  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.055;
236  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.038;
237  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.032;
238  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.052;
239  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
240  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.093;
241  }
242  if (type == kMuNeutralHadronIsoDR0p0To0p1) {
243  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
244  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.005;
245  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
246  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
247  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
248  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.003;
249  }
250  if (type == kMuNeutralHadronIsoDR0p1To0p2) {
251  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
252  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.006;
253  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
254  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
255  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.009;
256  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.013;
257  }
258  if (type == kMuNeutralHadronIsoDR0p2To0p3) {
259  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
260  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
261  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.015;
262  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
263  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.020;
264  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.024;
265  }
266  if (type == kMuNeutralHadronIsoDR0p3To0p4) {
267  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
268  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.019;
269  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
270  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.025;
271  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.030;
272  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.044;
273  }
274  if (type == kMuNeutralHadronIsoDR0p4To0p5) {
275  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.016;
276  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
277  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.030;
278  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
279  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.048;
280  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.118;
281  }
282  }
283 
284  //Fall11 MC Effective Areas
285  else if (EffectiveAreaTarget == kMuEAFall11MC) {
286  if (type == kMuGammaIsoDR0p0To0p1) {
287  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
288  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
289  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.003;
290  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
291  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
292  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.011;
293  }
294  if (type == kMuGammaIsoDR0p1To0p2) {
295  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
296  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
297  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
298  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.012;
299  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.019;
300  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.024;
301  }
302  if (type == kMuGammaIsoDR0p2To0p3) {
303  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.026;
304  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
305  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
306  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.022;
307  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.027;
308  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.034;
309  }
310  if (type == kMuGammaIsoDR0p3To0p4) {
311  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.042;
312  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
313  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.022;
314  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
315  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.059;
316  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.068;
317  }
318  if (type == kMuGammaIsoDR0p4To0p5) {
319  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.060;
320  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.043;
321  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.036;
322  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.055;
323  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.092;
324  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.115;
325  }
326  if (type == kMuNeutralHadronIsoDR0p0To0p1) {
327  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
328  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.004;
329  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
330  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.004;
331  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.010;
332  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.014;
333  }
334  if (type == kMuNeutralHadronIsoDR0p1To0p2) {
335  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
336  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
337  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
338  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
339  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.015;
340  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.017;
341  }
342  if (type == kMuNeutralHadronIsoDR0p2To0p3) {
343  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
344  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.015;
345  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
346  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.018;
347  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.022;
348  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.026;
349  }
350  if (type == kMuNeutralHadronIsoDR0p3To0p4) {
351  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
352  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
353  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.026;
354  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.032;
355  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.037;
356  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.042;
357  }
358  if (type == kMuNeutralHadronIsoDR0p4To0p5) {
359  if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
360  if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
361  if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.035;
362  if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.046;
363  if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.063;
364  if (fabs(SCEta) >= 2.3 ) EffectiveArea = 0.135;
365  }
366  }
367 
368 
369  return EffectiveArea;
370  }
371 };
372 
373 #endif
type
type
Definition: HCALResponse.h:21
MuonEffectiveArea::GetMuonEffectiveArea
static const Double_t GetMuonEffectiveArea(MuonEffectiveAreaType type, Double_t SCEta, MuonEffectiveAreaTarget EffectiveAreaTarget=kMuEAData2011)
Definition: MuonEffectiveArea.h:69
Generated for CMSSW Reference Manual by doxygen 1.8.5