ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/FGolf/Tools/ScaleFactors.cc

Comparing UserCode/FGolf/Tools/ScaleFactors.cc (file contents):
Revision 1.1 by fgolf, Mon May 28 15:18:22 2012 UTC vs.
Revision 1.4 by kelley, Tue May 29 18:42:48 2012 UTC

#	Line 1 | Line 1
1		#include "ScaleFactors.h"
2	<	#include "CORE/CMS2.h"
3	<	#include "bTagEff_BTV.h"
2	>	#include "CMS2.h"
3	>	#include "DileptonHypType.h"
4
5	<	float dileptonTagAndProbeScaleFactor(int idx)
5	>	float dileptonTagAndProbeScaleFactor(int hyp_idx)
6		{
7	<	float sf1 = tagAndProbeScaleFactor(cms2.hyp_lt_id().at(idx), cms2.hyp_lt_p4().at(idx).pt(), cms2.hyp_lt_p4().at(idx).eta());
8	<	float sf2 = tagAndProbeScaleFactor(cms2.hyp_ll_id().at(idx), cms2.hyp_ll_p4().at(idx).pt(), cms2.hyp_ll_p4().at(idx).eta());
7	>	float sf1 = tagAndProbeScaleFactor(cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx).pt(), cms2.hyp_lt_p4().at(hyp_idx).eta());
8	>	float sf2 = tagAndProbeScaleFactor(cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx).pt(), cms2.hyp_ll_p4().at(hyp_idx).eta());
9
10		return (sf1 * sf2);
11		}
#	Line 16 | Line 16 | float tagAndProbeScaleFactor(int id, flo
16		float feta = fabs(eta);
17
18		if (aid == 11) {
19	<	if (feta < .1479) {
20	<	if (pt < 40. && pt > 20.)
21	<	return 0.9989;
22	<	else if (pt > 40.)
23	<	return 0.9997;
19	>	if (feta < 0.1479) {
20	>	if (pt < 40.0 && pt > 20.0)
21	>	return 0.9947;
22	>	else if (pt > 40.0)
23	>	return 0.9967;
24		else
25	<	return 0.;
25	>	return 0.0;
26		}
27		else if (feta < 2.4) {
28	<	if (pt < 40. && pt > 20.)
29	<	return 1.0335;
30	<	else if (pt > 40.)
31	<	return 1.0114;
28	>	if (pt < 40.0 && pt > 20.0)
29	>	return 0.9638;
30	>	else if (pt > 40.0)
31	>	return 0.9768;
32		else
33	<	return 0.;
33	>	return 0.0;
34		}
35	<	else return 0.;
35	>	else return 0.0;
36		}
37		else if (aid == 13) {
38	<	if (feta < .1479) {
39	<	if (pt < 40. && pt > 20.)
40	<	return 1.0116;
38	>	if (feta < 0.1479) {
39	>	if (pt < 40.0 && pt > 20.0)
40	>	return 0.9863;
41		else if (pt > 40.)
42	<	return 1.0163;
42	>	return 0.9859;
43		else
44		return 0.;
45		}
46		else if (feta < 2.4) {
47		if (pt < 40. && pt > 20.)
48	<	return 0.9965;
49	<	else if (pt > 40.)
50	<	return 1.0114;
48	>	return 0.9941;
49	>	else if (pt > 40.0)
50	>	return 0.9921;
51		else
52	<	return 0.;
52	>	return 0.0;
53		}
54	<	else return 0.;
54	>	else return 0.0;
55		}
56
57	<	return 0.;
57	>	return 0.0;
58		}
59
60	+	float dilepTriggerScaleFactor(int hyp_idx)
61	+	{
62	+	double eff_1l = 0.0;
63	+	double eff_1t = 0.0;
64	+	double eff_2l = 0.0;
65	+	double eff_2t = 0.0;
66	+	double eff_z  = 0.0;
67	+
68	+	DileptonHypType hyp_type = hyp_typeToHypType(cms2.hyp_type().at(hyp_idx));
69	+
70	+	if (hyp_type == DILEPTON_EE || hyp_type == DILEPTON_MUMU) {
71	+	eff_1l  = triggerScaleFactor(LeadDbl, cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx));
72	+	eff_1t  = triggerScaleFactor(TrailDbl, cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx));
73	+	eff_2l  = triggerScaleFactor(LeadDbl, cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx));
74	+	eff_2t  = triggerScaleFactor(TrailDbl, cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx));
75	+	if (cms2.hyp_lt_p4().at(hyp_idx).pt() > cms2.hyp_ll_p4().at(hyp_idx).pt())
76	+	eff_z = triggerScaleFactor(DzDbl, cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx));
77	+	else
78	+	eff_z = triggerScaleFactor(DzDbl, cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx));
79	+
80	+	return (eff_1l * eff_2t + eff_1t * eff_2l - eff_1l * eff_2l) * eff_z;
81	+	}
82	+	else if (hyp_type == DILEPTON_EMU) {
83	+	eff_1l  = triggerScaleFactor(LeadDbl, cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx));
84	+	eff_1t  = triggerScaleFactor(TrailDbl, cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx));
85	+	eff_2l  = triggerScaleFactor(LeadDbl, cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx));
86	+	eff_2t  = triggerScaleFactor(TrailDbl, cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx));
87	+	if (abs(cms2.hyp_ll_id().at(hyp_idx)) == 13)
88	+	eff_z = triggerScaleFactor(DzDbl, cms2.hyp_ll_id().at(hyp_idx), cms2.hyp_ll_p4().at(hyp_idx));
89	+	else
90	+	eff_z = triggerScaleFactor(DzDbl, cms2.hyp_lt_id().at(hyp_idx), cms2.hyp_lt_p4().at(hyp_idx));
91	+
92	+	return (eff_1l * eff_2t + eff_1t * eff_2l - eff_1l * eff_2l) * sqrt(eff_z);
93	+	}
94	+	return 0.0;
95	+	}
96	+
97	+	float triggerScaleFactor(enum TrigEffType trig_eff_type, int id, const LorentzVector &p4)
98	+	{
99	+	unsigned int aid  = abs(id);
100	+	float eta = fabs(p4.eta());
101	+	float pt = p4.pt();
102	+
103	+	if (aid == 11)
104	+	{
105	+	switch (trig_eff_type)
106	+	{
107	+	case DzDbl:
108	+	if (fabs(eta) < 0.000) return 0.;
109	+	if (fabs(eta) > 2.500) return 0.;
110	+
111	+	if (pt < 20.000) return 0.;
112	+	if (pt > 7000.000) return 0.;
113	+
114	+	if (fabs(eta) < 0.800000) {
115	+	if (pt < 22.500000) return 1.000;
116	+	else if (pt < 25.000000) return 1.000;
117	+	else if (pt < 27.500000) return 0.999;
118	+	else if (pt < 30.000000) return 1.000;
119	+	else if (pt < 35.000000) return 1.000;
120	+	else if (pt < 40.000000) return 1.000;
121	+	else if (pt < 50.000000) return 1.000;
122	+	else if (pt < 7000.000000) return 1.000;
123	+	}
124	+	else if (fabs(eta) < 1.479000) {
125	+	if (pt < 22.500000) return 1.000;
126	+	else if (pt < 25.000000) return 1.000;
127	+	else if (pt < 27.500000) return 1.000;
128	+	else if (pt < 30.000000) return 0.999;
129	+	else if (pt < 35.000000) return 0.999;
130	+	else if (pt < 40.000000) return 0.999;
131	+	else if (pt < 50.000000) return 0.999;
132	+	else if (pt < 7000.000000) return 1.000;
133	+	}
134	+	else if (fabs(eta) < 2.000000) {
135	+	if (pt < 22.500000) return 1.000;
136	+	else if (pt < 25.000000) return 1.000;
137	+	else if (pt < 27.500000) return 0.994;
138	+	else if (pt < 30.000000) return 0.998;
139	+	else if (pt < 35.000000) return 0.999;
140	+	else if (pt < 40.000000) return 0.999;
141	+	else if (pt < 50.000000) return 0.998;
142	+	else if (pt < 7000.000000) return 0.997;
143	+	}
144	+	else if (fabs(eta) < 2.500000) {
145	+	if (pt < 22.500000) return 1.000;
146	+	else if (pt < 25.000000) return 1.000;
147	+	else if (pt < 27.500000) return 0.991;
148	+	else if (pt < 30.000000) return 0.987;
149	+	else if (pt < 35.000000) return 0.993;
150	+	else if (pt < 40.000000) return 0.987;
151	+	else if (pt < 50.000000) return 0.987;
152	+	else if (pt < 7000.000000) return 0.985;
153	+	}
154	+	break;
155	+	case LeadDbl:
156	+	if (fabs(eta) < 0.000) return 0.;
157	+	if (fabs(eta) > 2.500) return 0.;
158	+
159	+	if (pt < 20.000) return 0.;
160	+	if (pt > 7000.000) return 0.;
161	+
162	+	if (fabs(eta) < 0.800000) {
163	+	if (pt < 22.500000) return 0.929;
164	+	else if (pt < 25.000000) return 0.964;
165	+	else if (pt < 27.500000) return 0.961;
166	+	else if (pt < 30.000000) return 0.968;
167	+	else if (pt < 35.000000) return 0.972;
168	+	else if (pt < 40.000000) return 0.974;
169	+	else if (pt < 50.000000) return 0.975;
170	+	else if (pt < 7000.000000) return 0.977;
171	+	}
172	+	else if (fabs(eta) < 1.479000) {
173	+	if (pt < 22.500000) return 0.955;
174	+	else if (pt < 25.000000) return 0.972;
175	+	else if (pt < 27.500000) return 0.985;
176	+	else if (pt < 30.000000) return 0.985;
177	+	else if (pt < 35.000000) return 0.986;
178	+	else if (pt < 40.000000) return 0.992;
179	+	else if (pt < 50.000000) return 0.995;
180	+	else if (pt < 7000.000000) return 0.996;
181	+	}
182	+	else if (fabs(eta) < 2.000000) {
183	+	if (pt < 22.500000) return 0.980;
184	+	else if (pt < 25.000000) return 0.991;
185	+	else if (pt < 27.500000) return 0.986;
186	+	else if (pt < 30.000000) return 0.983;
187	+	else if (pt < 35.000000) return 0.990;
188	+	else if (pt < 40.000000) return 0.993;
189	+	else if (pt < 50.000000) return 0.995;
190	+	else if (pt < 7000.000000) return 0.994;
191	+	}
192	+	else if (fabs(eta) < 2.500000) {
193	+	if (pt < 22.500000) return 0.888;
194	+	else if (pt < 25.000000) return 0.988;
195	+	else if (pt < 27.500000) return 0.996;
196	+	else if (pt < 30.000000) return 0.992;
197	+	else if (pt < 35.000000) return 0.986;
198	+	else if (pt < 40.000000) return 0.989;
199	+	else if (pt < 50.000000) return 0.992;
200	+	else if (pt < 7000.000000) return 0.994;
201	+	}
202	+	break;
203	+	case Sgl:
204	+	if (fabs(eta) < 0.000) return 0.;
205	+	if (fabs(eta) > 2.500) return 0.;
206	+
207	+	if (pt < 20.000) return 0.;
208	+	if (pt > 7000.000) return 0.;
209	+
210	+	if (fabs(eta) < 0.800000) {
211	+	if (pt < 22.500000) return 0.000;
212	+	else if (pt < 25.000000) return 0.000;
213	+	else if (pt < 27.500000) return 0.059;
214	+	else if (pt < 30.000000) return 0.739;
215	+	else if (pt < 35.000000) return 0.906;
216	+	else if (pt < 40.000000) return 0.927;
217	+	else if (pt < 50.000000) return 0.944;
218	+	else if (pt < 7000.000000) return 0.952;
219	+	}
220	+	else if (fabs(eta) < 1.479000) {
221	+	if (pt < 22.500000) return 0.001;
222	+	else if (pt < 25.000000) return 0.028;
223	+	else if (pt < 27.500000) return 0.110;
224	+	else if (pt < 30.000000) return 0.660;
225	+	else if (pt < 35.000000) return 0.913;
226	+	else if (pt < 40.000000) return 0.940;
227	+	else if (pt < 50.000000) return 0.956;
228	+	else if (pt < 7000.000000) return 0.969;
229	+	}
230	+	else if (fabs(eta) < 2.000000) {
231	+	if (pt < 22.500000) return 0.000;
232	+	else if (pt < 25.000000) return 0.007;
233	+	else if (pt < 27.500000) return 0.113;
234	+	else if (pt < 30.000000) return 0.456;
235	+	else if (pt < 35.000000) return 0.691;
236	+	else if (pt < 40.000000) return 0.746;
237	+	else if (pt < 50.000000) return 0.767;
238	+	else if (pt < 7000.000000) return 0.769;
239	+	}
240	+	else if (fabs(eta) < 2.500000) {
241	+	if (pt < 22.500000) return 0.000;
242	+	else if (pt < 25.000000) return 0.003;
243	+	else if (pt < 27.500000) return 0.087;
244	+	else if (pt < 30.000000) return 0.313;
245	+	else if (pt < 35.000000) return 0.533;
246	+	else if (pt < 40.000000) return 0.646;
247	+	else if (pt < 50.000000) return 0.695;
248	+	else if (pt < 7000.000000) return 0.691;
249	+	}
250	+	break;
251	+	case TrailDbl:
252	+	if (fabs(eta) < 0.000) return 0.;
253	+	if (fabs(eta) > 2.500) return 0.;
254	+
255	+	if (pt < 20.000) return 0.;
256	+	if (pt > 7000.000) return 0.;
257	+
258	+	if (fabs(eta) < 0.800000) {
259	+	if (pt < 22.500000) return 0.948;
260	+	else if (pt < 25.000000) return 0.966;
261	+	else if (pt < 27.500000) return 0.963;
262	+	else if (pt < 30.000000) return 0.968;
263	+	else if (pt < 35.000000) return 0.969;
264	+	else if (pt < 40.000000) return 0.972;
265	+	else if (pt < 50.000000) return 0.972;
266	+	else if (pt < 7000.000000) return 0.976;
267	+	}
268	+	else if (fabs(eta) < 1.479000) {
269	+	if (pt < 22.500000) return 0.969;
270	+	else if (pt < 25.000000) return 0.972;
271	+	else if (pt < 27.500000) return 0.976;
272	+	else if (pt < 30.000000) return 0.982;
273	+	else if (pt < 35.000000) return 0.985;
274	+	else if (pt < 40.000000) return 0.986;
275	+	else if (pt < 50.000000) return 0.989;
276	+	else if (pt < 7000.000000) return 0.990;
277	+	}
278	+	else if (fabs(eta) < 2.000000) {
279	+	if (pt < 22.500000) return 0.973;
280	+	else if (pt < 25.000000) return 1.000;
281	+	else if (pt < 27.500000) return 0.983;
282	+	else if (pt < 30.000000) return 0.979;
283	+	else if (pt < 35.000000) return 0.988;
284	+	else if (pt < 40.000000) return 0.987;
285	+	else if (pt < 50.000000) return 0.990;
286	+	else if (pt < 7000.000000) return 0.991;
287	+	}
288	+	else if (fabs(eta) < 2.500000) {
289	+	if (pt < 22.500000) return 0.983;
290	+	else if (pt < 25.000000) return 0.981;
291	+	else if (pt < 27.500000) return 0.996;
292	+	else if (pt < 30.000000) return 0.989;
293	+	else if (pt < 35.000000) return 0.984;
294	+	else if (pt < 40.000000) return 0.986;
295	+	else if (pt < 50.000000) return 0.990;
296	+	else if (pt < 7000.000000) return 0.992;
297	+	}
298	+	break;
299	+	case N_TrigEffType:
300	+	default:
301	+	return 0;
302	+	};
303	+	}
304	+	else if (aid == 13)
305	+	{
306	+	switch (trig_eff_type)
307	+	{
308	+	case DzDbl:
309	+	if (fabs(eta) < 0.000) return 0.;
310	+	if (fabs(eta) > 2.400) return 0.;
311	+
312	+	if (pt < 20.000) return 0.;
313	+	if (pt > 7000.000) return 0.;
314	+
315	+	if (fabs(eta) < 0.800000) {
316	+	if (pt < 22.500000) return 0.925;
317	+	else if (pt < 25.000000) return 0.931;
318	+	else if (pt < 27.500000) return 0.923;
319	+	else if (pt < 30.000000) return 0.919;
320	+	else if (pt < 35.000000) return 0.926;
321	+	else if (pt < 40.000000) return 0.932;
322	+	else if (pt < 50.000000) return 0.936;
323	+	else if (pt < 7000.000000) return 0.941;
324	+	}
325	+	else if (fabs(eta) < 1.200000) {
326	+	if (pt < 22.500000) return 0.921;
327	+	else if (pt < 25.000000) return 0.926;
328	+	else if (pt < 27.500000) return 0.922;
329	+	else if (pt < 30.000000) return 0.931;
330	+	else if (pt < 35.000000) return 0.926;
331	+	else if (pt < 40.000000) return 0.911;
332	+	else if (pt < 50.000000) return 0.910;
333	+	else if (pt < 7000.000000) return 0.897;
334	+	}
335	+	else if (fabs(eta) < 2.400000) {
336	+	if (pt < 22.500000) return 0.879;
337	+	else if (pt < 25.000000) return 0.871;
338	+	else if (pt < 27.500000) return 0.874;
339	+	else if (pt < 30.000000) return 0.877;
340	+	else if (pt < 35.000000) return 0.880;
341	+	else if (pt < 40.000000) return 0.871;
342	+	else if (pt < 50.000000) return 0.854;
343	+	else if (pt < 7000.000000) return 0.846;
344	+	}
345	+	break;
346	+	case LeadDbl:
347	+	if (fabs(eta) < 0.000) return 0.;
348	+	if (fabs(eta) > 2.400) return 0.;
349	+
350	+	if (pt < 20.000) return 0.;
351	+	if (pt > 7000.000) return 0.;
352	+
353	+	if (fabs(eta) < 0.800000) {
354	+	if (pt < 22.500000) return 0.971;
355	+	else if (pt < 25.000000) return 0.972;
356	+	else if (pt < 27.500000) return 0.980;
357	+	else if (pt < 30.000000) return 0.970;
358	+	else if (pt < 35.000000) return 0.975;
359	+	else if (pt < 40.000000) return 0.977;
360	+	else if (pt < 50.000000) return 0.976;
361	+	else if (pt < 7000.000000) return 0.975;
362	+	}
363	+	else if (fabs(eta) < 1.200000) {
364	+	if (pt < 22.500000) return 0.905;
365	+	else if (pt < 25.000000) return 0.929;
366	+	else if (pt < 27.500000) return 0.950;
367	+	else if (pt < 30.000000) return 0.934;
368	+	else if (pt < 35.000000) return 0.935;
369	+	else if (pt < 40.000000) return 0.940;
370	+	else if (pt < 50.000000) return 0.942;
371	+	else if (pt < 7000.000000) return 0.935;
372	+	}
373	+	else if (fabs(eta) < 2.400000) {
374	+	if (pt < 22.500000) return 0.886;
375	+	else if (pt < 25.000000) return 0.927;
376	+	else if (pt < 27.500000) return 0.921;
377	+	else if (pt < 30.000000) return 0.916;
378	+	else if (pt < 35.000000) return 0.918;
379	+	else if (pt < 40.000000) return 0.921;
380	+	else if (pt < 50.000000) return 0.924;
381	+	else if (pt < 7000.000000) return 0.923;
382	+	}
383	+	break;
384	+	case Sgl:
385	+	if (fabs(eta) < 0.000) return 0.;
386	+	if (fabs(eta) > 2.400) return 0.;
387	+
388	+	if (pt < 20.000) return 0.;
389	+	if (pt > 7000.000) return 0.;
390	+
391	+	if (fabs(eta) < 0.800000) {
392	+	if (pt < 22.500000) return 0.001;
393	+	else if (pt < 25.000000) return 0.396;
394	+	else if (pt < 27.500000) return 0.886;
395	+	else if (pt < 30.000000) return 0.893;
396	+	else if (pt < 35.000000) return 0.914;
397	+	else if (pt < 40.000000) return 0.923;
398	+	else if (pt < 50.000000) return 0.935;
399	+	else if (pt < 7000.000000) return 0.945;
400	+	}
401	+	else if (fabs(eta) < 1.200000) {
402	+	if (pt < 22.500000) return 0.001;
403	+	else if (pt < 25.000000) return 0.369;
404	+	else if (pt < 27.500000) return 0.817;
405	+	else if (pt < 30.000000) return 0.823;
406	+	else if (pt < 35.000000) return 0.833;
407	+	else if (pt < 40.000000) return 0.854;
408	+	else if (pt < 50.000000) return 0.862;
409	+	else if (pt < 7000.000000) return 0.857;
410	+	}
411	+	else if (fabs(eta) < 2.400000) {
412	+	if (pt < 22.500000) return 0.006;
413	+	else if (pt < 25.000000) return 0.305;
414	+	else if (pt < 27.500000) return 0.597;
415	+	else if (pt < 30.000000) return 0.617;
416	+	else if (pt < 35.000000) return 0.624;
417	+	else if (pt < 40.000000) return 0.658;
418	+	else if (pt < 50.000000) return 0.703;
419	+	else if (pt < 7000.000000) return 0.711;
420	+	}
421	+	break;
422	+	case TrailDbl:
423	+	if (fabs(eta) < 0.000) return 0.;
424	+	if (fabs(eta) > 2.400) return 0.;
425	+
426	+	if (pt < 20.000) return 0.;
427	+	if (pt > 7000.000) return 0.;
428	+
429	+	if (fabs(eta) < 0.800000) {
430	+	if (pt < 22.500000) return 0.975;
431	+	else if (pt < 25.000000) return 0.973;
432	+	else if (pt < 27.500000) return 0.983;
433	+	else if (pt < 30.000000) return 0.971;
434	+	else if (pt < 35.000000) return 0.977;
435	+	else if (pt < 40.000000) return 0.978;
436	+	else if (pt < 50.000000) return 0.978;
437	+	else if (pt < 7000.000000) return 0.977;
438	+	}
439	+	else if (fabs(eta) < 1.200000) {
440	+	if (pt < 22.500000) return 0.925;
441	+	else if (pt < 25.000000) return 0.938;
442	+	else if (pt < 27.500000) return 0.962;
443	+	else if (pt < 30.000000) return 0.942;
444	+	else if (pt < 35.000000) return 0.945;
445	+	else if (pt < 40.000000) return 0.950;
446	+	else if (pt < 50.000000) return 0.950;
447	+	else if (pt < 7000.000000) return 0.943;
448	+	}
449	+	else if (fabs(eta) < 2.400000) {
450	+	if (pt < 22.500000) return 0.927;
451	+	else if (pt < 25.000000) return 0.951;
452	+	else if (pt < 27.500000) return 0.947;
453	+	else if (pt < 30.000000) return 0.942;
454	+	else if (pt < 35.000000) return 0.941;
455	+	else if (pt < 40.000000) return 0.938;
456	+	else if (pt < 50.000000) return 0.939;
457	+	else if (pt < 7000.000000) return 0.936;
458	+	}
459	+	break;
460	+	case N_TrigEffType:
461	+	default:
462	+	return 0;
463	+	};
464	+	}
465	+	else
466	+	{
467	+	return 0.0;
468	+	}
469	+	return 0.0;
470	+	}

Diff Legend

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