ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/FGolf/Tools/fliprate_egun.cc (file contents):
Revision 1.1 by fgolf, Fri Sep 10 21:36:36 2010 UTC vs.
Revision 1.3 by fgolf, Mon May 23 23:13:14 2011 UTC

#	Line 15 | Line 15
15		#include <iostream>
16		#include <stdio.h>
17		#include <math.h>
18	+	#include <cmath>
19
20		using namespace std;
21
#	Line 22 | Line 23 | double getSingleEleNum(double el_pt, dou
23
24		el_eta = fabs(el_eta);
25
26	<	if( el_eta < 1.28 ){
27	<	if( el_pt > 70 ) return 31;
28	<	if( el_pt > 50 ) return 20;
29	<	if( el_pt > 40 ) return 5;
26	>	if( el_eta < 0.5 ){
27	>	if( el_pt > 95 ) return 9;
28	>	if( el_pt > 90 ) return 5;
29	>	if( el_pt > 85 ) return 4;
30	>	if( el_pt > 80 ) return 6;
31	>	if( el_pt > 75 ) return 7;
32	>	if( el_pt > 70 ) return 5;
33	>	if( el_pt > 65 ) return 8;
34	>	if( el_pt > 60 ) return 4;
35	>	if( el_pt > 55 ) return 4;
36	>	if( el_pt > 50 ) return 3;
37	>	if( el_pt > 45 ) return 6;
38	>	if( el_pt > 40 ) return 6;
39	>	if( el_pt > 35 ) return 6;
40	>	if( el_pt > 30 ) return 3;
41	>	if( el_pt > 25 ) return 5;
42	>	if( el_pt > 20 ) return 0;
43	>	if( el_pt > 10 ) return 1;
44	>	return 0.0;
45	>	}
46	>	if( el_eta < 1 ){
47	>	if( el_pt > 95 ) return 17;
48	>	if( el_pt > 90 ) return 25;
49	>	if( el_pt > 85 ) return 16;
50	>	if( el_pt > 80 ) return 9;
51	>	if( el_pt > 75 ) return 18;
52	>	if( el_pt > 70 ) return 10;
53	>	if( el_pt > 65 ) return 17;
54	>	if( el_pt > 60 ) return 15;
55	>	if( el_pt > 55 ) return 13;
56	>	if( el_pt > 50 ) return 21;
57	>	if( el_pt > 45 ) return 14;
58	>	if( el_pt > 40 ) return 7;
59	>	if( el_pt > 35 ) return 8;
60		if( el_pt > 30 ) return 9;
61	<	if( el_pt > 10 ) return 4;
61	>	if( el_pt > 25 ) return 12;
62	>	if( el_pt > 20 ) return 4;
63	>	if( el_pt > 10 ) return 2;
64		return 0.0;
65		}
66	<	if( el_eta < 1.56 ){
67	<	if( el_pt > 70 ) return 48;
68	<	if( el_pt > 50 ) return 30;
69	<	if( el_pt > 40 ) return 12;
70	<	if( el_pt > 30 ) return 8;
71	<	if( el_pt > 10 ) return 5;
66	>	if( el_eta < 1.479 ){
67	>	if( el_pt > 95 ) return 88;
68	>	if( el_pt > 90 ) return 76;
69	>	if( el_pt > 85 ) return 95;
70	>	if( el_pt > 80 ) return 73;
71	>	if( el_pt > 75 ) return 84;
72	>	if( el_pt > 70 ) return 72;
73	>	if( el_pt > 65 ) return 75;
74	>	if( el_pt > 60 ) return 75;
75	>	if( el_pt > 55 ) return 77;
76	>	if( el_pt > 50 ) return 65;
77	>	if( el_pt > 45 ) return 45;
78	>	if( el_pt > 40 ) return 39;
79	>	if( el_pt > 35 ) return 48;
80	>	if( el_pt > 30 ) return 33;
81	>	if( el_pt > 25 ) return 25;
82	>	if( el_pt > 20 ) return 19;
83	>	if( el_pt > 10 ) return 8;
84		return 0.0;
85		}
86	<	if( el_eta < 1.84 ){
87	<	if( el_pt > 70 ) return 92;
88	<	if( el_pt > 50 ) return 67;
89	<	if( el_pt > 40 ) return 35;
90	<	if( el_pt > 30 ) return 20;
91	<	if( el_pt > 10 ) return 13;
86	>	if( el_eta < 1.8 ){
87	>	if( el_pt > 95 ) return 144;
88	>	if( el_pt > 90 ) return 138;
89	>	if( el_pt > 85 ) return 129;
90	>	if( el_pt > 80 ) return 116;
91	>	if( el_pt > 75 ) return 122;
92	>	if( el_pt > 70 ) return 103;
93	>	if( el_pt > 65 ) return 111;
94	>	if( el_pt > 60 ) return 108;
95	>	if( el_pt > 55 ) return 90;
96	>	if( el_pt > 50 ) return 93;
97	>	if( el_pt > 45 ) return 91;
98	>	if( el_pt > 40 ) return 69;
99	>	if( el_pt > 35 ) return 46;
100	>	if( el_pt > 30 ) return 44;
101	>	if( el_pt > 25 ) return 30;
102	>	if( el_pt > 20 ) return 16;
103	>	if( el_pt > 10 ) return 9;
104		return 0.0;
105		}
106	<	if( el_eta < 2.12 ){
107	<	if( el_pt > 70 ) return 72;
108	<	if( el_pt > 50 ) return 29;
109	<	if( el_pt > 40 ) return 16;
110	<	if( el_pt > 30 ) return 13;
106	>	if( el_eta < 2 ){
107	>	if( el_pt > 95 ) return 70;
108	>	if( el_pt > 90 ) return 60;
109	>	if( el_pt > 85 ) return 70;
110	>	if( el_pt > 80 ) return 60;
111	>	if( el_pt > 75 ) return 56;
112	>	if( el_pt > 70 ) return 74;
113	>	if( el_pt > 65 ) return 58;
114	>	if( el_pt > 60 ) return 54;
115	>	if( el_pt > 55 ) return 59;
116	>	if( el_pt > 50 ) return 60;
117	>	if( el_pt > 45 ) return 51;
118	>	if( el_pt > 40 ) return 43;
119	>	if( el_pt > 35 ) return 25;
120	>	if( el_pt > 30 ) return 28;
121	>	if( el_pt > 25 ) return 17;
122	>	if( el_pt > 20 ) return 10;
123		if( el_pt > 10 ) return 6;
124		return 0.0;
125		}
126	<	if( el_eta < = 2.5 ){
127	<	if( el_pt > 70 ) return 93;
128	<	if( el_pt > 50 ) return 42;
126	>	if( el_eta < 2.1 ){
127	>	if( el_pt > 95 ) return 40;
128	>	if( el_pt > 90 ) return 33;
129	>	if( el_pt > 85 ) return 33;
130	>	if( el_pt > 80 ) return 29;
131	>	if( el_pt > 75 ) return 28;
132	>	if( el_pt > 70 ) return 23;
133	>	if( el_pt > 65 ) return 27;
134	>	if( el_pt > 60 ) return 21;
135	>	if( el_pt > 55 ) return 27;
136	>	if( el_pt > 50 ) return 15;
137	>	if( el_pt > 45 ) return 18;
138		if( el_pt > 40 ) return 17;
139	<	if( el_pt > 30 ) return 9;
140	<	if( el_pt > 10 ) return 20;
139	>	if( el_pt > 35 ) return 11;
140	>	if( el_pt > 30 ) return 12;
141	>	if( el_pt > 25 ) return 8;
142	>	if( el_pt > 20 ) return 8;
143	>	if( el_pt > 10 ) return 6;
144		return 0.0;
145		}
146	<	std::cout << "Error: eta > 2.5 value found" << endl;
146	>	if( el_eta < 2.2 ){
147	>	if( el_pt > 95 ) return 35;
148	>	if( el_pt > 90 ) return 35;
149	>	if( el_pt > 85 ) return 27;
150	>	if( el_pt > 80 ) return 29;
151	>	if( el_pt > 75 ) return 36;
152	>	if( el_pt > 70 ) return 20;
153	>	if( el_pt > 65 ) return 24;
154	>	if( el_pt > 60 ) return 21;
155	>	if( el_pt > 55 ) return 19;
156	>	if( el_pt > 50 ) return 32;
157	>	if( el_pt > 45 ) return 21;
158	>	if( el_pt > 40 ) return 19;
159	>	if( el_pt > 35 ) return 12;
160	>	if( el_pt > 30 ) return 10;
161	>	if( el_pt > 25 ) return 13;
162	>	if( el_pt > 20 ) return 5;
163	>	if( el_pt > 10 ) return 9;
164	>	return 0.0;
165	>	}
166	>	if( el_eta <= 2.4 ){
167	>	if( el_pt > 95 ) return 64;
168	>	if( el_pt > 90 ) return 84;
169	>	if( el_pt > 85 ) return 68;
170	>	if( el_pt > 80 ) return 64;
171	>	if( el_pt > 75 ) return 56;
172	>	if( el_pt > 70 ) return 57;
173	>	if( el_pt > 65 ) return 54;
174	>	if( el_pt > 60 ) return 47;
175	>	if( el_pt > 55 ) return 38;
176	>	if( el_pt > 50 ) return 46;
177	>	if( el_pt > 45 ) return 39;
178	>	if( el_pt > 40 ) return 29;
179	>	if( el_pt > 35 ) return 23;
180	>	if( el_pt > 30 ) return 25;
181	>	if( el_pt > 25 ) return 14;
182	>	if( el_pt > 20 ) return 20;
183	>	if( el_pt > 10 ) return 14;
184	>	return 0.0;
185	>	}
186	>	std::cout << "Error: eta > 2.4 value found" << endl;
187		return 0.0;
188		}
189
190
191		double getSingleEleDenom(double el_pt, double el_eta) {
192
193	<	el_eta = fabs(el_eta);
193	>	el_eta = std::fabs(el_eta);
194
195	<	if( el_eta < 1.28 ){
196	<	if( el_pt > 70 ) return 117286;
197	<	if( el_pt > 50 ) return 78397;
198	<	if( el_pt > 40 ) return 39090;
199	<	if( el_pt > 30 ) return 38858;
200	<	if( el_pt > 10 ) return 69323;
201	<	return 0.0;
202	<	}
203	<	if( el_eta < 1.56 ){
204	<	if( el_pt > 70 ) return 18452;
205	<	if( el_pt > 50 ) return 12276;
206	<	if( el_pt > 40 ) return 6087;
207	<	if( el_pt > 30 ) return 5793;
208	<	if( el_pt > 10 ) return 9374;
209	<	return 0.0;
210	<	}
211	<	if( el_eta < 1.84 ){
212	<	if( el_pt > 70 ) return 20720;
213	<	if( el_pt > 50 ) return 14174;
214	<	if( el_pt > 40 ) return 7021;
215	<	if( el_pt > 30 ) return 6852;
216	<	if( el_pt > 10 ) return 10145;
217	<	return 0.0;
218	<	}
219	<	if( el_eta < 2.12 ){
220	<	if( el_pt > 70 ) return 21436;
221	<	if( el_pt > 50 ) return 14537;
222	<	if( el_pt > 40 ) return 7401;
223	<	if( el_pt > 30 ) return 7241;
224	<	if( el_pt > 10 ) return 12178;
225	<	return 0.0;
226	<	}
227	<	if( el_eta < = 2.5 ){
228	<	if( el_pt > 70 ) return 25029;
229	<	if( el_pt > 50 ) return 16802;
230	<	if( el_pt > 40 ) return 8462;
231	<	if( el_pt > 30 ) return 8277;
232	<	if( el_pt > 10 ) return 15646;
195	>	if( el_eta < 0.5 ){
196	>	if( el_pt > 95 ) return 80174;
197	>	if( el_pt > 90 ) return 79910;
198	>	if( el_pt > 85 ) return 79820;
199	>	if( el_pt > 80 ) return 79722;
200	>	if( el_pt > 75 ) return 79890;
201	>	if( el_pt > 70 ) return 79653;
202	>	if( el_pt > 65 ) return 79892;
203	>	if( el_pt > 60 ) return 80077;
204	>	if( el_pt > 55 ) return 80011;
205	>	if( el_pt > 50 ) return 79597;
206	>	if( el_pt > 45 ) return 80053;
207	>	if( el_pt > 40 ) return 79263;
208	>	if( el_pt > 35 ) return 79239;
209	>	if( el_pt > 30 ) return 78273;
210	>	if( el_pt > 25 ) return 77894;
211	>	if( el_pt > 20 ) return 75100;
212	>	if( el_pt > 10 ) return 114108;
213	>	return 0.0;
214	>	}
215	>	if( el_eta < 1 ){
216	>	if( el_pt > 95 ) return 79315;
217	>	if( el_pt > 90 ) return 79561;
218	>	if( el_pt > 85 ) return 79406;
219	>	if( el_pt > 80 ) return 79842;
220	>	if( el_pt > 75 ) return 79654;
221	>	if( el_pt > 70 ) return 79671;
222	>	if( el_pt > 65 ) return 79640;
223	>	if( el_pt > 60 ) return 79599;
224	>	if( el_pt > 55 ) return 79540;
225	>	if( el_pt > 50 ) return 79343;
226	>	if( el_pt > 45 ) return 79944;
227	>	if( el_pt > 40 ) return 79623;
228	>	if( el_pt > 35 ) return 78557;
229	>	if( el_pt > 30 ) return 78303;
230	>	if( el_pt > 25 ) return 77378;
231	>	if( el_pt > 20 ) return 75091;
232	>	if( el_pt > 10 ) return 122322;
233	>	return 0.0;
234	>	}
235	>	if( el_eta < 1.479 ){
236	>	if( el_pt > 95 ) return 71523;
237	>	if( el_pt > 90 ) return 70668;
238	>	if( el_pt > 85 ) return 71662;
239	>	if( el_pt > 80 ) return 71434;
240	>	if( el_pt > 75 ) return 71467;
241	>	if( el_pt > 70 ) return 71204;
242	>	if( el_pt > 65 ) return 70951;
243	>	if( el_pt > 60 ) return 71044;
244	>	if( el_pt > 55 ) return 70664;
245	>	if( el_pt > 50 ) return 70559;
246	>	if( el_pt > 45 ) return 70218;
247	>	if( el_pt > 40 ) return 69528;
248	>	if( el_pt > 35 ) return 68507;
249	>	if( el_pt > 30 ) return 67672;
250	>	if( el_pt > 25 ) return 65776;
251	>	if( el_pt > 20 ) return 61672;
252	>	if( el_pt > 10 ) return 93412;
253	>	return 0.0;
254	>	}
255	>	if( el_eta < 1.8 ){
256	>	if( el_pt > 95 ) return 37802;
257	>	if( el_pt > 90 ) return 38002;
258	>	if( el_pt > 85 ) return 37636;
259	>	if( el_pt > 80 ) return 38035;
260	>	if( el_pt > 75 ) return 38432;
261	>	if( el_pt > 70 ) return 37862;
262	>	if( el_pt > 65 ) return 37989;
263	>	if( el_pt > 60 ) return 37686;
264	>	if( el_pt > 55 ) return 37738;
265	>	if( el_pt > 50 ) return 37321;
266	>	if( el_pt > 45 ) return 37092;
267	>	if( el_pt > 40 ) return 36538;
268	>	if( el_pt > 35 ) return 35582;
269	>	if( el_pt > 30 ) return 34886;
270	>	if( el_pt > 25 ) return 32819;
271	>	if( el_pt > 20 ) return 30064;
272	>	if( el_pt > 10 ) return 36945;
273	>	return 0.0;
274	>	}
275	>	if( el_eta < 2 ){
276	>	if( el_pt > 95 ) return 27683;
277	>	if( el_pt > 90 ) return 27689;
278	>	if( el_pt > 85 ) return 27736;
279	>	if( el_pt > 80 ) return 27917;
280	>	if( el_pt > 75 ) return 27924;
281	>	if( el_pt > 70 ) return 27763;
282	>	if( el_pt > 65 ) return 28018;
283	>	if( el_pt > 60 ) return 27977;
284	>	if( el_pt > 55 ) return 28120;
285	>	if( el_pt > 50 ) return 27549;
286	>	if( el_pt > 45 ) return 27622;
287	>	if( el_pt > 40 ) return 27251;
288	>	if( el_pt > 35 ) return 26803;
289	>	if( el_pt > 30 ) return 26001;
290	>	if( el_pt > 25 ) return 25559;
291	>	if( el_pt > 20 ) return 23249;
292	>	if( el_pt > 10 ) return 29668;
293	>	return 0.0;
294	>	}
295	>	if( el_eta < 2.1 ){
296	>	if( el_pt > 95 ) return 13372;
297	>	if( el_pt > 90 ) return 13612;
298	>	if( el_pt > 85 ) return 13518;
299	>	if( el_pt > 80 ) return 13317;
300	>	if( el_pt > 75 ) return 13707;
301	>	if( el_pt > 70 ) return 13562;
302	>	if( el_pt > 65 ) return 13637;
303	>	if( el_pt > 60 ) return 13714;
304	>	if( el_pt > 55 ) return 13635;
305	>	if( el_pt > 50 ) return 13628;
306	>	if( el_pt > 45 ) return 13481;
307	>	if( el_pt > 40 ) return 13384;
308	>	if( el_pt > 35 ) return 13437;
309	>	if( el_pt > 30 ) return 13341;
310	>	if( el_pt > 25 ) return 13083;
311	>	if( el_pt > 20 ) return 12750;
312	>	if( el_pt > 10 ) return 17362;
313	>	return 0.0;
314	>	}
315	>	if( el_eta < 2.2 ){
316	>	if( el_pt > 95 ) return 12545;
317	>	if( el_pt > 90 ) return 12730;
318	>	if( el_pt > 85 ) return 12914;
319	>	if( el_pt > 80 ) return 12765;
320	>	if( el_pt > 75 ) return 12899;
321	>	if( el_pt > 70 ) return 12897;
322	>	if( el_pt > 65 ) return 12758;
323	>	if( el_pt > 60 ) return 12824;
324	>	if( el_pt > 55 ) return 13093;
325	>	if( el_pt > 50 ) return 12773;
326	>	if( el_pt > 45 ) return 12724;
327	>	if( el_pt > 40 ) return 12958;
328	>	if( el_pt > 35 ) return 12976;
329	>	if( el_pt > 30 ) return 12877;
330	>	if( el_pt > 25 ) return 12626;
331	>	if( el_pt > 20 ) return 12046;
332	>	if( el_pt > 10 ) return 17226;
333	>	return 0.0;
334	>	}
335	>	if( el_eta <= 2.4 ){
336	>	if( el_pt > 95 ) return 23896;
337	>	if( el_pt > 90 ) return 24146;
338	>	if( el_pt > 85 ) return 24045;
339	>	if( el_pt > 80 ) return 24087;
340	>	if( el_pt > 75 ) return 24197;
341	>	if( el_pt > 70 ) return 24369;
342	>	if( el_pt > 65 ) return 24486;
343	>	if( el_pt > 60 ) return 24379;
344	>	if( el_pt > 55 ) return 24598;
345	>	if( el_pt > 50 ) return 24470;
346	>	if( el_pt > 45 ) return 24531;
347	>	if( el_pt > 40 ) return 24214;
348	>	if( el_pt > 35 ) return 24295;
349	>	if( el_pt > 30 ) return 24178;
350	>	if( el_pt > 25 ) return 23611;
351	>	if( el_pt > 20 ) return 23022;
352	>	if( el_pt > 10 ) return 31178;
353		return 0.0;
354		}
355		std::cout << "Error: eta > 2.5 value found" << endl;

Diff Legend

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