New response functions for tracking updates
[u/mrichter/AliRoot.git] / TRD / AliTRDprobdist.cxx
CommitLineData
33a848f6 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16//-----------------------------------------------------------------
17// Class for dE/dx and Time Bin of Max. Cluster for Electrons and
18// pions in TRD.
19// It is instantiated in class AliTRDpidESD for particle identification
20// in TRD
21// Prashant Shukla <shukla@pi0.physi.uni-heidelberg.de>
22//-----------------------------------------------------------------
23
24#include "AliTRDprobdist.h"
25
26ClassImp(AliTRDprobdist)
27
28//_________________________________________________________________________
1642e650 29AliTRDprobdist::AliTRDprobdist(Int_t multiplicity)
33a848f6 30{
31 //
32 // The main constructor
33 //
1642e650 34 if (multiplicity == 1) FillData();
35 // if (multiplicity == 2000) FillData2000();
36 // if (multiplicity == 4000) FillData4000();
37 // if (multiplicity == 6000) FillData6000();
38 // if (multiplicity == 8000) FillData8000();
33a848f6 39}
40
41//_________________________________________________________________________
1642e650 42Double_t AliTRDprobdist::GetMeanPI(Int_t ip) const
33a848f6 43{
44 // Gets mean of de/dx dist. of pi
45 Double_t integral=0.;
46 Double_t norm=0.;
47 for(Int_t ie=0; ie<fNEbins; ie++) {
48 integral+=fEnergyLoss[ip][ie]*fProbPI[ip][ie];
49 norm+=fProbPI[ip][ie];
50 }
51 return integral/norm;
52}
53
54//_________________________________________________________________________
1642e650 55Double_t AliTRDprobdist::GetMeanEL(Int_t ip) const
33a848f6 56{
57 //
58 // Gets mean of de/dx dist. of e
59 Double_t integral=0.;
60 Double_t norm=0.;
61 for(Int_t ie=0; ie<fNEbins; ie++) {
62 integral+=fEnergyLoss[ip][ie]*fProbEL[ip][ie];
63 norm+=fProbEL[ip][ie];
64 }
65 return integral/norm;
66}
67
68//_________________________________________________________________________
1642e650 69Double_t AliTRDprobdist::GetNormalizationPI(Int_t ip) const
33a848f6 70{
71 //
72 // Gets Normalization of de/dx dist. of pi
73 Double_t integral=0.;
74 for(Int_t ie=0; ie<fNEbins; ie++) {
75 integral+=fProbPI[ip][ie];
76 }
77 return integral;
78}
79
80//_________________________________________________________________________
1642e650 81Double_t AliTRDprobdist::GetNormalizationEL(Int_t ip) const
33a848f6 82{
83 //
84 // Gets Normalization of de/dx dist. of e
85 Double_t integral=0.;
86 for(Int_t ie=0; ie<fNEbins; ie++) {
87 integral+=fProbEL[ip][ie];
88 }
89 return integral;
90}
91
92//_________________________________________________________________________
1642e650 93Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx) const
33a848f6 94{
95 //
96 // Gets the Probability of having dedx at a given momentum (mom)
97 // and particle type k (0 for e) and (2 for pi)
98 // from the precalculated de/dx distributions
1642e650 99 Double_t probability = 1.0;
33a848f6 100 Int_t iEnBin= ((Int_t) (dedx/fEnBinSize));
101 if(iEnBin > fNEbins-1) iEnBin = fNEbins-1;
102
103 if(k==0){ // electron
1642e650 104 if(mom<=fTrackMomentum[0]) probability = fProbEL[0][iEnBin];
105 if(mom>=fTrackMomentum[fNMom-1]) probability = fProbEL[fNMom-1][iEnBin];
33a848f6 106 }
107 if(k==2){ // pion
1642e650 108 if(mom<=fTrackMomentum[0]) probability = fProbPI[0][iEnBin];
109 if(mom>=fTrackMomentum[fNMom-1]) probability = fProbPI[fNMom-1][iEnBin];
33a848f6 110 }
111
112 if(k==0) // electron
113 for(Int_t ip=1; ip<fNMom; ip++)
114 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
115 Double_t slop=(fProbEL[ip][iEnBin]-fProbEL[ip-1][iEnBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
116 // Linear Interpolation
1642e650 117 probability= fProbEL[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
118 return probability;
33a848f6 119 }
120
121 if(k==2) //pion
122 for(Int_t ip=1; ip<fNMom; ip++)
123 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
124 Double_t slop=(fProbPI[ip][iEnBin]-fProbPI[ip-1][iEnBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
125 // Linear Interpolation
1642e650 126 probability= fProbPI[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
127 return probability;
33a848f6 128 }
1642e650 129 return probability;
33a848f6 130}
131
132
133//_________________________________________________________________________
1642e650 134Double_t AliTRDprobdist::GetProbabilityT(Int_t k, Double_t mom, Int_t timbin) const
33a848f6 135{
136 //
137 // Gets the Probability of having timbin at a given momentum (mom)
138 // and particle type k (0 for e) and (2 for pi)
139 // from the precalculated timbin distributions
1642e650 140 Double_t probabilityT = 1.0;
33a848f6 141 if(timbin<=0) return 0.;
142 Int_t iTBin=timbin;
143
144 if(k==0){ // electron
1642e650 145 if(mom<=fTrackMomentum[0]) probabilityT = fProbELT[0][iTBin];
146 if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbELT[fNMom-1][iTBin];
33a848f6 147 }
148 if(k==2){ // pion
1642e650 149 if(mom<=fTrackMomentum[0]) probabilityT = fProbPIT[0][iTBin];
150 if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbPIT[fNMom-1][iTBin];
33a848f6 151 }
152
153 if(k==0) // electron
154 for(Int_t ip=1; ip<fNMom; ip++)
155 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
156 Double_t slop=(fProbELT[ip][iTBin]-fProbELT[ip-1][iTBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
157 // Linear Interpolation
1642e650 158 probabilityT= fProbELT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
159 return probabilityT;
33a848f6 160 }
161
162 if(k==2) // pion
163 for(Int_t ip=1; ip<fNMom; ip++)
164 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
165 Double_t slop=(fProbPIT[ip][iTBin]-fProbPIT[ip-1][iTBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
166 // Linear Interpolation
1642e650 167 probabilityT= fProbPIT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
168 return probabilityT;
33a848f6 169 }
1642e650 170 return probabilityT;
33a848f6 171}
172
173
33a848f6 174void AliTRDprobdist::FillData()
175{
176 //
177 // Energy loss Distributions for e and pi
178 fNMom=7;
32333a5e 179 Double_t trackMomentum[kNo_Mom]= {1, 1.5, 2, 3, 4, 5, 6};
180 fNEbins=250;
1642e650 181 fEnBinSize=10;
33a848f6 182
1642e650 183 Double_t energyLoss[kNo_Mom][kNo_EnBins]={
33a848f6 184 {
185 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
186 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
187 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
188 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
189 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
190 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
191 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
192 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
193 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
194 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
195 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
196 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
197 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
198 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
199 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
200 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
201 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
202 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
203 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 204 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
205 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
206 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
207 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
208 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
209 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 210 },
211 {
212 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
213 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
214 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
215 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
216 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
217 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
218 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
219 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
220 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
221 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
222 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
223 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
224 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
225 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
226 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
227 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
228 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
229 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
230 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 231 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
232 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
233 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
234 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
235 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
236 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 237 },
238 {
239 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
240 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
241 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
242 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
243 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
244 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
245 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
246 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
247 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
248 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
249 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
250 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
251 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
252 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
253 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
254 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
255 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
256 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
257 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 258 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
259 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
260 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
261 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
262 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
263 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 264 },
265 {
266 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
267 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
268 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
269 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
270 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
271 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
272 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
273 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
274 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
275 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
276 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
277 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
278 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
279 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
280 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
281 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
282 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
283 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
284 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 285 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
286 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
287 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
288 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
289 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
290 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 291 },
292 {
293 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
294 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
295 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
296 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
297 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
298 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
299 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
300 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
301 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
302 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
303 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
304 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
305 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
306 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
307 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
308 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
309 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
310 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
311 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 312 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
313 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
314 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
315 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
316 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
317 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 318 },
319 {
320 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
321 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
322 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
323 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
324 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
325 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
326 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
327 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
328 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
329 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
330 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
331 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
332 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
333 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
334 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
335 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
336 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
337 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
338 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 339 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
340 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
341 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
342 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
343 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
344 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 345 },
346 {
347 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
348 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
349 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
350 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
351 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
352 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
353 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
354 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
355 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
356 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
357 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
358 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
359 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
360 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
361 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
362 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
363 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
364 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
365 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
32333a5e 366 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990,
367 2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090,
368 2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190,
369 2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290,
370 2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390,
371 2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
33a848f6 372 }
373 };
374
1642e650 375 Double_t probPI[kNo_Mom][kNo_EnBins]={
33a848f6 376 {
32333a5e 377 0, 0, 0, 0, 0.000133547,
33a848f6 378 0, 0, 0, 0, 0,
32333a5e 379 0.000133547, 0, 0, 0, 0,
380 0.000667735, 0.000267094, 0.000934829, 0.00186966, 0.00280449,
381 0.00413996, 0.00534188, 0.0117521, 0.0132212, 0.0205662,
382 0.0206998, 0.0220353, 0.0251068, 0.0331197, 0.0335203,
383 0.0363248, 0.0371261, 0.0413996, 0.0403312, 0.0376603,
384 0.0419338, 0.0351229, 0.0345887, 0.0301816, 0.0347222,
385 0.0319177, 0.0308494, 0.0295139, 0.0247062, 0.0257746,
386 0.0220353, 0.0200321, 0.0202991, 0.0189637, 0.0153579,
387 0.0140224, 0.0150908, 0.0121528, 0.0125534, 0.0110844,
388 0.0101496, 0.00894765, 0.00734509, 0.0068109, 0.00467415,
389 0.00667735, 0.00534188, 0.00320513, 0.00320513, 0.00333868,
390 0.00400641, 0.00307158, 0.00280449, 0.00240385, 0.00333868,
391 0.00253739, 0.00240385, 0.00240385, 0.00186966, 0.00106838,
392 0.00240385, 0.000667735, 0.00200321, 0.00133547, 0.00173611,
393 0.00133547, 0.000534188, 0.00106838, 0.000534188, 0.000667735,
394 0.000801282, 0.00106838, 0.000667735, 0.00106838, 0.000934829,
395 0.000400641, 0.00133547, 0.00106838, 0.000801282, 0.000801282,
396 0.000667735, 0.000667735, 0.000534188, 0.000133547, 0.000534188,
397 0.000133547, 0.000667735, 0.000534188, 0.000400641, 0.000400641,
398 0.000400641, 0.00106838, 0.000133547, 0.000267094, 0.000534188,
399 0.000400641, 0.000133547, 0.000534188, 0.000267094, 0.000133547,
400 0.000133547, 0, 0.000133547, 0.000400641, 0.000400641,
401 0.000400641, 0.000400641, 0.000133547, 0.000133547, 0.000133547,
402 0.000534188, 0, 0.000133547, 0, 0.000267094,
403 0.000133547, 0.000133547, 0, 0.000267094, 0.000267094,
404 0.000133547, 0.000133547, 0.000801282, 0, 0.000400641,
405 0.000267094, 0.000133547, 0, 0.000267094, 0.000267094,
406 0.000267094, 0.000267094, 0, 0.000133547, 0.000133547,
407 0.000133547, 0, 0, 0, 0,
408 0.000133547, 0, 0, 0, 0.000267094,
409 0, 0.000133547, 0, 0, 0.000133547,
410 0.000400641, 0, 0.000133547, 0.000133547, 0.000133547,
411 0.000133547, 0, 0.000133547, 0.000133547, 0,
412 0.000133547, 0.000133547, 0, 0, 0.000133547,
413 0.000133547, 0.000267094, 0.000133547, 0, 0.000400641,
414 0.000133547, 0.000400641, 0, 0, 0,
415 0, 0.000267094, 0, 0, 0.000133547,
33a848f6 416 0, 0, 0, 0, 0,
32333a5e 417 0.000133547, 0.000400641, 0.000133547, 0, 0.000133547,
418 0.000133547, 0, 0.000133547, 0.000133547, 0,
419 0.000133547, 0, 0, 0.000133547, 0.000133547,
420 0, 0.000133547, 0, 0, 0,
421 0.000133547, 0, 0, 0, 0,
422 0, 0, 0, 0, 0,
423 0, 0.000133547, 0, 0, 0,
424 0, 0, 0, 0, 0,
425 0, 0.000133547, 0, 0, 0,
426 0, 0, 0, 0, 0
33a848f6 427 },
428 {
429 0, 0, 0, 0, 0,
430 0, 0, 0, 0, 0,
32333a5e 431 0.000159642, 0, 0.000159642, 0.00063857, 0.000159642,
432 0, 0, 0.000159642, 0.00063857, 0.00175607,
433 0.00303321, 0.00462963, 0.00718391, 0.00846105, 0.00989783,
434 0.0164432, 0.019636, 0.0244253, 0.024106, 0.0316092,
435 0.0359195, 0.0371967, 0.0351213, 0.0387931, 0.0403895,
436 0.0386335, 0.0371967, 0.0320881, 0.0371967, 0.0320881,
437 0.0341635, 0.0333653, 0.0304917, 0.0301724, 0.0252235,
438 0.0228289, 0.0204342, 0.0231481, 0.0189974, 0.0205939,
439 0.0177203, 0.0150064, 0.0143678, 0.0119732, 0.0127714,
440 0.0116539, 0.0102171, 0.0113346, 0.00670498, 0.0063857,
441 0.00718391, 0.00718391, 0.00686462, 0.0052682, 0.0063857,
442 0.00431034, 0.00319285, 0.00399106, 0.00287356, 0.00367178,
443 0.00287356, 0.00207535, 0.00127714, 0.000957854, 0.00207535,
444 0.000957854, 0.00143678, 0.00175607, 0.00223499, 0.00175607,
445 0.0011175, 0.00127714, 0.00127714, 0.00143678, 0.000798212,
446 0.0011175, 0.000159642, 0.000478927, 0.000478927, 0.000957854,
447 0.00127714, 0.000798212, 0.000478927, 0.000319285, 0.000319285,
448 0.00063857, 0.000159642, 0.000478927, 0.000478927, 0.000478927,
449 0.00063857, 0.000319285, 0.0011175, 0.000159642, 0.000159642,
450 0.000478927, 0.000319285, 0.000159642, 0.000159642, 0.000159642,
451 0.000159642, 0.000798212, 0.000319285, 0.000319285, 0.000159642,
452 0.000319285, 0.000159642, 0.000159642, 0, 0.000478927,
453 0, 0, 0.000159642, 0.000319285, 0.000319285,
454 0.000478927, 0.000159642, 0.000319285, 0, 0.000159642,
455 0.000159642, 0.000319285, 0.000319285, 0.000159642, 0.000159642,
456 0.000159642, 0.000319285, 0.000159642, 0, 0.000319285,
457 0.000478927, 0.000319285, 0.000159642, 0, 0.000159642,
458 0.000319285, 0.000159642, 0.000159642, 0.000159642, 0.000159642,
459 0.000159642, 0.000319285, 0, 0.000159642, 0.000319285,
460 0.000159642, 0, 0.000319285, 0, 0.000478927,
461 0.000319285, 0, 0.000159642, 0, 0.000319285,
462 0, 0, 0.000159642, 0.000319285, 0,
463 0.000478927, 0.000159642, 0, 0.000159642, 0.000159642,
464 0, 0, 0.000159642, 0.000319285, 0,
465 0, 0.000319285, 0, 0.000319285, 0,
466 0, 0, 0, 0.000159642, 0,
467 0.000159642, 0.000159642, 0, 0, 0,
468 0, 0, 0.000159642, 0.000159642, 0,
469 0, 0, 0.000159642, 0, 0,
470 0.000159642, 0, 0, 0.000159642, 0,
471 0.000159642, 0, 0, 0, 0,
472 0.000159642, 0, 0, 0, 0,
473 0.000159642, 0, 0, 0, 0,
474 0, 0, 0.000159642, 0, 0,
475 0, 0.000159642, 0, 0, 0,
476 0, 0, 0, 0, 0.000159642,
477 0, 0, 0.000159642, 0, 0,
478 0.000159642, 0, 0, 0, 0
33a848f6 479 },
480 {
481 0, 0, 0, 0, 0,
32333a5e 482 0, 0, 0, 0, 0,
483 0, 0.000120773, 0, 0.000120773, 0,
484 0.000120773, 0.000362319, 0.000483092, 0.000724638, 0.00120773,
485 0.0013285, 0.0031401, 0.00326087, 0.00845411, 0.00809179,
486 0.0141304, 0.018599, 0.0218599, 0.021256, 0.026087,
487 0.0288647, 0.0318841, 0.0342995, 0.0349034, 0.0369565,
488 0.0359903, 0.0363527, 0.035628, 0.0347826, 0.0375604,
489 0.0369565, 0.0324879, 0.0321256, 0.0306763, 0.0298309,
490 0.0277778, 0.0240338, 0.0225845, 0.0183575, 0.0176329,
491 0.0210145, 0.0187198, 0.0150966, 0.0146135, 0.0111111,
492 0.0119565, 0.0129227, 0.00917874, 0.00954106, 0.00821256,
493 0.00688406, 0.00785024, 0.00688406, 0.00603865, 0.00615942,
494 0.00434783, 0.00422705, 0.00410628, 0.00338164, 0.00301932,
495 0.0044686, 0.00301932, 0.00181159, 0.00289855, 0.00229469,
496 0.00181159, 0.00241546, 0.00217391, 0.0013285, 0.00169082,
497 0.000845411, 0.00108696, 0.00108696, 0.000483092, 0.00205314,
498 0.00108696, 0.0013285, 0.000241546, 0.000845411, 0.000966184,
499 0.000724638, 0.0013285, 0.000241546, 0.000603865, 0.000724638,
500 0.00169082, 0.000483092, 0.000483092, 0.000724638, 0.000603865,
501 0.000724638, 0.000483092, 0.000362319, 0.000483092, 0.000483092,
502 0.000845411, 0.000362319, 0.000362319, 0.000724638, 0.000966184,
503 0.000966184, 0.000362319, 0.000120773, 0.000483092, 0.000362319,
504 0.000241546, 0.000724638, 0.000241546, 0.000362319, 0.000483092,
505 0, 0.000483092, 0.000362319, 0, 0,
506 0.000603865, 0.000483092, 0.000120773, 0.000362319, 0,
507 0.000241546, 0.000241546, 0.000120773, 0.000120773, 0.000120773,
508 0.000120773, 0.000241546, 0.000362319, 0.000241546, 0.000120773,
509 0.000483092, 0.000120773, 0.000362319, 0.000362319, 0.000120773,
510 0.000120773, 0.000362319, 0, 0, 0.000120773,
511 0.000603865, 0.000241546, 0.000362319, 0.000120773, 0,
512 0.000120773, 0.000120773, 0.000241546, 0.000120773, 0,
513 0.000120773, 0, 0.000120773, 0.000120773, 0,
514 0.000120773, 0, 0.000241546, 0.000120773, 0.000120773,
515 0.000120773, 0.000362319, 0.000120773, 0.000241546, 0,
516 0, 0, 0, 0, 0.000120773,
517 0, 0, 0, 0, 0,
518 0, 0, 0, 0, 0,
519 0, 0.000120773, 0.000120773, 0, 0,
520 0.000120773, 0, 0, 0.000120773, 0,
521 0, 0, 0.000120773, 0, 0,
522 0.000120773, 0.000120773, 0, 0, 0.000241546,
523 0, 0.000120773, 0, 0.000241546, 0.000120773,
524 0, 0, 0, 0, 0,
525 0.000241546, 0, 0, 0, 0,
526 0.000241546, 0, 0, 0, 0,
527 0.000120773, 0, 0, 0, 0,
528 0, 0.000120773, 0, 0.000120773, 0,
529 0, 0, 0, 0, 0,
530 0, 0, 0, 0, 0
33a848f6 531 },
532 {
32333a5e 533 0, 0, 0, 0, 0.000235073,
534 0.000117536, 0, 0, 0, 0.000117536,
535 0, 0.000235073, 0, 0.000117536, 0.000235073,
536 0.000117536, 0.000117536, 0.000352609, 0.000117536, 0.000235073,
537 0.00105783, 0.00141044, 0.00152797, 0.00376117, 0.00564175,
538 0.00693465, 0.0109309, 0.0105783, 0.0188058, 0.0185708,
539 0.0256229, 0.0245651, 0.0303244, 0.030677, 0.0332628,
540 0.0404325, 0.0342031, 0.0353785, 0.0349083, 0.0378467,
541 0.0352609, 0.0338505, 0.0343206, 0.0326751, 0.0330277,
542 0.0293841, 0.0287964, 0.0278561, 0.0235073, 0.0224495,
543 0.0209215, 0.0209215, 0.01622, 0.0183357, 0.0152797,
544 0.01622, 0.0124589, 0.0131641, 0.0112835, 0.00987306,
545 0.0102257, 0.00822755, 0.00716972, 0.00811001, 0.0064645,
546 0.00681711, 0.0051716, 0.00423131, 0.0038787, 0.00423131,
547 0.00434885, 0.00329102, 0.00235073, 0.00317348, 0.00235073,
548 0.00211566, 0.00211566, 0.00211566, 0.0012929, 0.00164551,
549 0.0012929, 0.0012929, 0.000940291, 0.0012929, 0.00176305,
550 0.000705219, 0.00141044, 0.000705219, 0.00105783, 0.00117536,
551 0.00105783, 0.000822755, 0.000587682, 0.000822755, 0.00117536,
552 0.000235073, 0.000470146, 0.000470146, 0.000470146, 0.000587682,
553 0.000470146, 0.000587682, 0.000352609, 0.000940291, 0.000352609,
554 0.000587682, 0.000352609, 0.000470146, 0.000470146, 0.000352609,
555 0.000352609, 0.000470146, 0.000235073, 0.000235073, 0.000587682,
556 0.000117536, 0.000235073, 0.000352609, 0.000235073, 0.000235073,
557 0.000235073, 0.000235073, 0.000352609, 0.000352609, 0.000470146,
558 0.000470146, 0.000705219, 0.000470146, 0.000470146, 0.000117536,
559 0.000470146, 0.000235073, 0.000352609, 0, 0.000235073,
560 0.000235073, 0.000470146, 0.000117536, 0, 0.000352609,
561 0, 0.000117536, 0, 0.000470146, 0,
562 0.000235073, 0.000352609, 0.000235073, 0.000235073, 0,
563 0.000235073, 0.000117536, 0, 0, 0.000117536,
564 0.000352609, 0.000470146, 0, 0.000587682, 0,
565 0, 0.000352609, 0.000352609, 0.000352609, 0.000235073,
566 0.000235073, 0.000235073, 0.000352609, 0.000235073, 0.000235073,
567 0.000117536, 0, 0.000235073, 0, 0.000352609,
568 0, 0.000117536, 0.000117536, 0.000352609, 0,
569 0, 0, 0.000235073, 0.000117536, 0,
570 0, 0.000117536, 0, 0, 0,
571 0, 0, 0.000235073, 0, 0,
572 0, 0, 0, 0.000117536, 0.000117536,
573 0, 0, 0.000117536, 0.000117536, 0.000117536,
574 0, 0.000117536, 0.000117536, 0.000117536, 0,
575 0.000235073, 0, 0.000117536, 0, 0,
576 0, 0.000117536, 0.000117536, 0, 0,
577 0, 0.000117536, 0.000117536, 0.000117536, 0.000117536,
578 0, 0.000235073, 0.000117536, 0.000117536, 0,
579 0.000117536, 0, 0, 0, 0.000117536,
580 0, 0, 0.000117536, 0.000117536, 0,
33a848f6 581 0, 0, 0, 0, 0,
32333a5e 582 0, 0, 0.000117536, 0, 0
33a848f6 583 },
584 {
585 0, 0, 0, 0, 0,
32333a5e 586 0, 0, 0, 0, 0.000107112,
587 0, 0.000321337, 0, 0, 0.000214225,
588 0.000214225, 0, 0, 0.000107112, 0.000214225,
589 0.000856898, 0.000749786, 0.00160668, 0.00214225, 0.0035347,
590 0.00739075, 0.00771208, 0.0108183, 0.0138175, 0.0153171,
591 0.0194944, 0.0220651, 0.0250643, 0.0292416, 0.0302057,
592 0.0323479, 0.0397386, 0.0345973, 0.0366324, 0.0347044,
593 0.0393102, 0.0368466, 0.0351328, 0.0315981, 0.0320266,
594 0.0292416, 0.0293488, 0.0298843, 0.0244216, 0.0254927,
595 0.0247429, 0.0224936, 0.0214225, 0.018102, 0.0147815,
596 0.0128535, 0.0152099, 0.0145673, 0.0121037, 0.0123179,
597 0.00996144, 0.0103899, 0.00814053, 0.00846187, 0.00664096,
598 0.00803342, 0.00685518, 0.00685518, 0.0052485, 0.0043916,
599 0.00589117, 0.00342759, 0.00299914, 0.00278492, 0.00299914,
600 0.00321337, 0.00299914, 0.00289203, 0.00278492, 0.00160668,
601 0.00117823, 0.0017138, 0.00289203, 0.0017138, 0.000749786,
602 0.00139246, 0.00107112, 0.00139246, 0.000642674, 0.000642674,
603 0.000428449, 0.000428449, 0.000535561, 0.000321337, 0.000428449,
604 0.000642674, 0.000749786, 0.000642674, 0.000856898, 0.000214225,
605 0.000428449, 0.000214225, 0.000749786, 0.000214225, 0.000642674,
606 0.000428449, 0.000535561, 0.000428449, 0.000749786, 0.000214225,
607 0.000535561, 0.000535561, 0.000535561, 0.000428449, 0.000856898,
608 0.000535561, 0.000321337, 0.000214225, 0.000535561, 0.000321337,
609 0.000107112, 0.000107112, 0.000107112, 0.000642674, 0.000107112,
610 0.000321337, 0.000214225, 0.000214225, 0, 0.000321337,
611 0.000107112, 0.000428449, 0.000214225, 0.000107112, 0.000214225,
612 0.000321337, 0.000214225, 0.000107112, 0.000107112, 0.000321337,
613 0, 0.000214225, 0.000321337, 0.000107112, 0,
614 0.000214225, 0.000107112, 0.000214225, 0.000321337, 0.000428449,
615 0.000321337, 0.000214225, 0.000107112, 0, 0,
616 0.000321337, 0.000107112, 0.000321337, 0.000107112, 0,
617 0.000321337, 0, 0, 0.000214225, 0.000107112,
618 0.000321337, 0.000107112, 0.000214225, 0.000214225, 0,
619 0, 0.000107112, 0, 0, 0,
620 0, 0.000321337, 0.000214225, 0, 0,
621 0, 0, 0, 0.000214225, 0,
622 0, 0, 0.000107112, 0, 0.000107112,
623 0, 0.000107112, 0.000214225, 0.000107112, 0,
624 0.000107112, 0, 0, 0, 0.000107112,
625 0, 0.000107112, 0, 0.000214225, 0.000107112,
626 0.000107112, 0, 0.000107112, 0.000107112, 0,
627 0, 0, 0.000107112, 0, 0,
628 0, 0, 0, 0.000107112, 0,
629 0, 0, 0.000107112, 0.000321337, 0,
630 0, 0.000214225, 0, 0, 0,
631 0, 0, 0, 0, 0,
632 0, 0, 0, 0, 0,
633 0.000107112, 0, 0, 0.000214225, 0,
634 0, 0.000107112, 0.000107112, 0, 0
33a848f6 635 },
636 {
637 0, 0, 0, 0, 0,
32333a5e 638 0, 0, 0.00011558, 0, 0,
639 0, 0, 0, 0.00011558, 0.00011558,
640 0.00011558, 0.00011558, 0.000346741, 0.00023116, 0.00023116,
641 0.000346741, 0.000924642, 0.00127138, 0.00127138, 0.00300509,
642 0.00369857, 0.00612575, 0.00705039, 0.00809061, 0.0122515,
643 0.0162968, 0.0177994, 0.0238095, 0.0271613, 0.0309755,
644 0.0293574, 0.0308599, 0.0381415, 0.0365233, 0.0340962,
645 0.0351364, 0.0328248, 0.0336338, 0.0357143, 0.0330559,
646 0.0329404, 0.0292418, 0.0275081, 0.0286639, 0.0269302,
647 0.0247342, 0.0206889, 0.02092, 0.0206889, 0.0184928,
648 0.0181461, 0.018724, 0.0154877, 0.0150254, 0.0128294,
649 0.0149098, 0.00901526, 0.0105178, 0.00751271, 0.00947758,
650 0.00728155, 0.00520111, 0.00866852, 0.00554785, 0.00658807,
651 0.00566343, 0.00358299, 0.00473879, 0.00312067, 0.00473879,
652 0.00358299, 0.00277393, 0.0023116, 0.00219602, 0.00127138,
653 0.00242718, 0.00138696, 0.00196486, 0.0011558, 0.00127138,
654 0.00127138, 0.00219602, 0.00138696, 0.000462321, 0.0011558,
655 0.000924642, 0.00104022, 0.00127138, 0.00196486, 0.000577901,
656 0.000693481, 0.000693481, 0.000577901, 0.000809061, 0.000346741,
657 0.000693481, 0.000693481, 0.000577901, 0.0011558, 0.000577901,
658 0.000346741, 0.000577901, 0.000462321, 0.000346741, 0.000809061,
659 0.000346741, 0.00023116, 0.00011558, 0.000346741, 0.000577901,
660 0.00011558, 0.000346741, 0.00023116, 0.000346741, 0.000346741,
661 0.00011558, 0.00011558, 0.000577901, 0.00011558, 0.000462321,
662 0.000809061, 0.00023116, 0.00011558, 0.00011558, 0.00011558,
663 0.00011558, 0.00011558, 0.000346741, 0.00023116, 0,
664 0.000462321, 0.00011558, 0, 0.000346741, 0.00023116,
665 0.00023116, 0.000577901, 0.00011558, 0.00023116, 0,
666 0.000346741, 0.00023116, 0.00023116, 0, 0.00023116,
667 0.00011558, 0.000346741, 0.00011558, 0.00011558, 0,
668 0.00023116, 0.00023116, 0.00023116, 0.00011558, 0.00011558,
669 0.00011558, 0.00011558, 0.00023116, 0, 0.00023116,
670 0, 0.00011558, 0.00011558, 0, 0,
671 0.00011558, 0.00011558, 0.00011558, 0.00023116, 0.000346741,
672 0.000346741, 0.00011558, 0, 0.00023116, 0.00011558,
673 0.00011558, 0, 0, 0.00023116, 0,
674 0.00011558, 0, 0, 0, 0,
675 0, 0.00011558, 0, 0, 0.00011558,
676 0.00023116, 0.00011558, 0, 0, 0.00011558,
677 0.00011558, 0, 0, 0, 0,
678 0.00011558, 0.00011558, 0.00011558, 0, 0,
679 0.00011558, 0.00011558, 0.00023116, 0, 0,
680 0, 0.00023116, 0, 0, 0.00011558,
681 0.00011558, 0.00011558, 0.00011558, 0, 0,
682 0, 0, 0.00011558, 0, 0,
683 0, 0, 0.00011558, 0, 0,
684 0, 0, 0, 0, 0.00011558,
685 0, 0.00011558, 0, 0, 0.00011558,
686 0, 0, 0, 0, 0
33a848f6 687 },
688 {
32333a5e 689 0, 0, 0, 0.00010582, 0.00010582,
690 0, 0, 0.00010582, 0, 0,
691 0, 0.00010582, 0, 0, 0,
692 0.00010582, 0, 0, 0.00021164, 0.00021164,
693 0.00010582, 0.000529101, 0.00010582, 0.00169312, 0.00116402,
694 0.00169312, 0.00465608, 0.00560847, 0.00878307, 0.0114286,
695 0.0122751, 0.0183069, 0.0201058, 0.0207407, 0.0266667,
696 0.0273016, 0.0316402, 0.0283598, 0.0331217, 0.0380952,
697 0.0342857, 0.0362963, 0.0361905, 0.035873, 0.0356614,
698 0.0320635, 0.0330159, 0.0326984, 0.0288889, 0.0251852,
699 0.0245503, 0.0274074, 0.0249735, 0.0197884, 0.0201058,
700 0.0180952, 0.0201058, 0.0169312, 0.0156614, 0.0143915,
701 0.0126984, 0.0116402, 0.00867725, 0.0108995, 0.0108995,
702 0.00814815, 0.00867725, 0.00740741, 0.00592593, 0.00497354,
703 0.00582011, 0.00497354, 0.00560847, 0.00412698, 0.00380952,
704 0.0031746, 0.00275132, 0.00232804, 0.00296296, 0.00232804,
705 0.00179894, 0.00190476, 0.00179894, 0.00275132, 0.0021164,
706 0.0010582, 0.00116402, 0.00148148, 0.00137566, 0.00137566,
707 0.00116402, 0.000846561, 0.00031746, 0.000529101, 0.000846561,
708 0.0010582, 0.00042328, 0.000634921, 0.000740741, 0.00031746,
709 0.000529101, 0.00021164, 0.000634921, 0.00021164, 0.00021164,
710 0.00042328, 0.00010582, 0.000529101, 0.00010582, 0.000634921,
711 0.00021164, 0.00042328, 0.00031746, 0.000529101, 0.000952381,
712 0.00021164, 0, 0.00031746, 0.00010582, 0.00031746,
713 0.00042328, 0.00031746, 0.00031746, 0.000529101, 0.00010582,
714 0.000529101, 0.000529101, 0.000529101, 0.00010582, 0.00021164,
715 0.00010582, 0.00021164, 0.00010582, 0.00010582, 0.00021164,
716 0.00010582, 0.00031746, 0.00021164, 0.00021164, 0,
717 0.000529101, 0, 0, 0.00010582, 0.00010582,
718 0.000634921, 0, 0.00010582, 0, 0.00031746,
719 0, 0.00010582, 0, 0, 0,
720 0.00031746, 0.00010582, 0.00021164, 0.000529101, 0.00021164,
721 0.00021164, 0.00010582, 0.00010582, 0, 0.00010582,
722 0.00021164, 0, 0, 0, 0,
723 0.00010582, 0.00021164, 0.00031746, 0, 0.00021164,
724 0.00010582, 0.00010582, 0, 0, 0,
725 0.00010582, 0.00010582, 0.00010582, 0, 0,
726 0, 0.00021164, 0.00010582, 0.00021164, 0.00021164,
727 0.00021164, 0, 0, 0.00010582, 0.00010582,
728 0, 0, 0.00021164, 0, 0,
729 0, 0.00010582, 0.00021164, 0.00010582, 0.00010582,
730 0, 0, 0, 0.00010582, 0.00010582,
731 0.00010582, 0, 0.00021164, 0, 0,
732 0.00010582, 0.00010582, 0, 0.00010582, 0,
33a848f6 733 0, 0, 0, 0, 0,
32333a5e 734 0, 0, 0, 0.00010582, 0,
735 0.00010582, 0.00031746, 0, 0, 0,
736 0, 0.00010582, 0, 0, 0.00010582,
737 0, 0.00010582, 0, 0, 0,
738 0, 0, 0, 0, 0.00010582
33a848f6 739 }
740 };
741
1642e650 742 Double_t probEL[kNo_Mom][kNo_EnBins]={
33a848f6 743 {
744 0, 0, 0, 0, 0,
32333a5e 745 0.000141603, 0, 0, 0, 0,
746 0, 0.000141603, 0.000283206, 0.000141603, 0,
747 0, 0, 0.000141603, 0.000141603, 0,
748 0, 0.000141603, 0, 0, 0.000141603,
749 0.000141603, 0.000708015, 0.000424809, 0.000849618, 0.000849618,
750 0.00155763, 0.00240725, 0.000991221, 0.00354007, 0.00580572,
751 0.00608893, 0.00736335, 0.0107618, 0.011045, 0.0104786,
752 0.0127443, 0.0184084, 0.0167091, 0.0155763, 0.0209572,
753 0.01855, 0.0239309, 0.0225149, 0.0244973, 0.0252053,
754 0.0247805, 0.0242141, 0.0291702, 0.0237893, 0.0278958,
755 0.0226565, 0.026763, 0.0230813, 0.0223733, 0.0257717,
756 0.0230813, 0.0215236, 0.0210988, 0.0198244, 0.0196828,
757 0.0191164, 0.0175588, 0.0179836, 0.0198244, 0.017134,
758 0.0155763, 0.0128859, 0.0147267, 0.0130275, 0.0127443,
759 0.0128859, 0.0120363, 0.00764656, 0.00877938, 0.00991221,
760 0.00920419, 0.00778816, 0.00708015, 0.00623053, 0.00679694,
761 0.00538091, 0.00566412, 0.00693854, 0.0049561, 0.00538091,
762 0.00566412, 0.00424809, 0.00438969, 0.00523931, 0.00254885,
763 0.0046729, 0.00354007, 0.00254885, 0.00311526, 0.00269046,
764 0.00226565, 0.00354007, 0.00212404, 0.00212404, 0.00254885,
765 0.00226565, 0.00240725, 0.00184084, 0.00226565, 0.00198244,
766 0.000991221, 0.00113282, 0.000566412, 0.00127443, 0.00155763,
767 0.00127443, 0.000849618, 0.00113282, 0.000849618, 0.000566412,
768 0.00141603, 0.00113282, 0.00155763, 0.000566412, 0.00113282,
769 0.000991221, 0.000566412, 0.000424809, 0.000424809, 0.000283206,
770 0.000849618, 0.000283206, 0.000991221, 0.000141603, 0.000566412,
771 0.000708015, 0.000708015, 0.000283206, 0.000283206, 0.000424809,
772 0.000283206, 0.000708015, 0.000424809, 0.000708015, 0.000141603,
773 0.000566412, 0, 0.000708015, 0.000283206, 0.000283206,
774 0.000141603, 0.000566412, 0.000424809, 0, 0.000424809,
775 0.000141603, 0, 0.000283206, 0.000424809, 0.000283206,
776 0.000141603, 0.000283206, 0.000141603, 0, 0.000708015,
777 0.000283206, 0, 0.000283206, 0.000141603, 0.000141603,
778 0.000141603, 0.000283206, 0, 0.000283206, 0.000141603,
779 0.000141603, 0.000141603, 0.000141603, 0.000424809, 0.000141603,
780 0.000141603, 0.000283206, 0.000141603, 0, 0.000283206,
781 0.000141603, 0.000141603, 0.000424809, 0, 0.000141603,
782 0.000283206, 0, 0.000283206, 0, 0,
783 0, 0, 0.000283206, 0.000283206, 0.000283206,
784 0, 0, 0.000141603, 0, 0.000141603,
785 0, 0, 0, 0, 0.000141603,
786 0.000141603, 0, 0.000141603, 0.000141603, 0.000141603,
787 0, 0, 0.000141603, 0, 0,
788 0, 0, 0.000283206, 0, 0,
789 0, 0, 0, 0, 0.000141603,
1642e650 790 0, 0, 0, 0, 0,
32333a5e 791 0, 0.000141603, 0, 0, 0,
792 0, 0.000141603, 0, 0, 0,
793 0.000141603, 0, 0, 0.000141603, 0
33a848f6 794 },
795 {
796 0, 0, 0, 0, 0,
32333a5e 797 0, 0.000138427, 0.000138427, 0, 0,
798 0, 0.000138427, 0.00055371, 0, 0,
799 0, 0, 0.000138427, 0, 0.000276855,
800 0.000276855, 0, 0.000138427, 0.00055371, 0.000138427,
801 0, 0.000276855, 0.000276855, 0.00055371, 0.00124585,
802 0.00110742, 0.00221484, 0.00373754, 0.0030454, 0.00484496,
803 0.00609081, 0.00567553, 0.00733666, 0.00816722, 0.01232,
804 0.0110742, 0.0160576, 0.0138427, 0.0160576, 0.0159192,
805 0.0186877, 0.021041, 0.021041, 0.0184109, 0.0221484,
806 0.0202104, 0.019103, 0.020072, 0.0236711, 0.0211794,
807 0.0235327, 0.0213178, 0.021041, 0.0224252, 0.0235327,
808 0.0227021, 0.0164729, 0.0218715, 0.0213178, 0.0184109,
809 0.0203488, 0.0202104, 0.0164729, 0.0135659, 0.0145349,
810 0.0153654, 0.014258, 0.0146733, 0.0163344, 0.0141196,
811 0.0130122, 0.0141196, 0.0131506, 0.0102436, 0.00899779,
812 0.00968992, 0.0101052, 0.0127353, 0.0095515, 0.0102436,
813 0.0095515, 0.00927464, 0.00816722, 0.00636766, 0.0055371,
814 0.00609081, 0.00636766, 0.00484496, 0.00498339, 0.00442968,
815 0.00263012, 0.0055371, 0.00442968, 0.00332226, 0.00359911,
816 0.00318383, 0.00373754, 0.00373754, 0.00193798, 0.00263012,
817 0.00263012, 0.00332226, 0.00276855, 0.00221484, 0.00207641,
818 0.00346069, 0.00166113, 0.00207641, 0.00221484, 0.00166113,
819 0.00207641, 0.00166113, 0.00193798, 0.00138427, 0.00166113,
820 0.000968992, 0.000692137, 0.00124585, 0.00124585, 0.0015227,
821 0.00207641, 0.0015227, 0.0015227, 0.00124585, 0.000830565,
822 0.000415282, 0.000692137, 0.00055371, 0.000692137, 0.00055371,
823 0.000968992, 0.000968992, 0.00110742, 0.000968992, 0.000692137,
824 0.000415282, 0.000692137, 0.000968992, 0.00055371, 0.000830565,
825 0.000415282, 0.000415282, 0.000968992, 0.000968992, 0.000415282,
826 0.000692137, 0.000276855, 0.000138427, 0.00055371, 0.000138427,
827 0.00055371, 0.000415282, 0.000276855, 0.000138427, 0,
828 0.00055371, 0.000138427, 0.000415282, 0.000415282, 0.000415282,
829 0.000968992, 0.000276855, 0.000415282, 0.000415282, 0.000138427,
830 0.000138427, 0.000415282, 0.000276855, 0.000138427, 0,
831 0.000276855, 0, 0.000138427, 0.000276855, 0.000830565,
832 0.000415282, 0.000415282, 0, 0.000138427, 0.000138427,
833 0, 0.000276855, 0.000138427, 0.000138427, 0,
834 0.000138427, 0.000138427, 0.000276855, 0, 0.000138427,
835 0, 0.000138427, 0, 0.000138427, 0.000138427,
836 0, 0.000138427, 0.000138427, 0.000138427, 0,
33a848f6 837 0, 0, 0, 0, 0,
32333a5e 838 0.000138427, 0, 0.000138427, 0, 0,
839 0, 0.000138427, 0.000276855, 0.000138427, 0,
840 0, 0.000276855, 0, 0.000138427, 0.000138427,
841 0.000138427, 0, 0, 0, 0,
33a848f6 842 0, 0, 0, 0, 0,
32333a5e 843 0.000276855, 0.000138427, 0, 0, 0,
1642e650 844 0, 0, 0, 0, 0,
32333a5e 845 0, 0, 0.000276855, 0.000276855, 0
33a848f6 846 },
847 {
32333a5e 848 0, 0, 0.000228781, 0, 0.00011439,
849 0, 0, 0.00011439, 0, 0.00011439,
850 0, 0, 0.000228781, 0.00011439, 0,
851 0, 0, 0.00011439, 0, 0.00011439,
852 0, 0.000343171, 0, 0.000228781, 0.00011439,
853 0, 0.000228781, 0.000800732, 0.000571951, 0.000686342,
854 0.00102951, 0.00125829, 0.0024022, 0.00205903, 0.00537634,
855 0.00514756, 0.00583391, 0.00754976, 0.00674903, 0.0117822,
856 0.0133837, 0.0108671, 0.0104095, 0.0122398, 0.0153283,
857 0.0179593, 0.0217342, 0.0171585, 0.0186456, 0.0181881,
858 0.0232212, 0.021391, 0.0237932, 0.0224205, 0.0226493,
859 0.0256234, 0.0202471, 0.0194464, 0.0235644, 0.0186456,
860 0.021391, 0.0195607, 0.0210478, 0.0202471, 0.0178449,
861 0.0155571, 0.0220773, 0.0171585, 0.0173873, 0.0175017,
862 0.0173873, 0.0152139, 0.0141844, 0.0147563, 0.0136124,
863 0.0134981, 0.0128117, 0.0125829, 0.0101807, 0.0112102,
864 0.00972318, 0.0101807, 0.00869366, 0.00926561, 0.0082361,
865 0.00926561, 0.0070922, 0.00663464, 0.00743537, 0.00629147,
866 0.00789293, 0.00846488, 0.00732098, 0.00571951, 0.00526195,
867 0.0059483, 0.0059483, 0.00549073, 0.00526195, 0.00491878,
868 0.0024022, 0.00343171, 0.00423244, 0.00400366, 0.00469,
869 0.00388927, 0.00308854, 0.00331732, 0.00320293, 0.00251659,
870 0.00205903, 0.00251659, 0.00205903, 0.00171585, 0.00171585,
871 0.00148707, 0.00263098, 0.00137268, 0.00274537, 0.0011439,
872 0.00183024, 0.00137268, 0.00102951, 0.000800732, 0.00102951,
873 0.000571951, 0.00160146, 0.00194464, 0.00137268, 0.0011439,
874 0.00148707, 0.000915122, 0.000800732, 0.00137268, 0.000800732,
875 0.000686342, 0.000457561, 0.000800732, 0.000571951, 0.0011439,
876 0.000343171, 0.000800732, 0.000800732, 0.000571951, 0.000686342,
877 0.000343171, 0.000457561, 0.000457561, 0.000571951, 0.000571951,
878 0.000571951, 0.000457561, 0.000571951, 0.00011439, 0.000457561,
879 0.000571951, 0.000571951, 0.000228781, 0.000343171, 0.000228781,
880 0.000686342, 0.000457561, 0.00011439, 0.000228781, 0.000571951,
881 0.000571951, 0.00011439, 0.000228781, 0.000571951, 0.000343171,
882 0.00011439, 0.000228781, 0.000228781, 0.000228781, 0.000228781,
883 0.000228781, 0.000457561, 0.000343171, 0.00011439, 0.00011439,
884 0, 0.00011439, 0, 0.000228781, 0.000343171,
885 0, 0.000228781, 0.00011439, 0, 0,
886 0.00011439, 0.00011439, 0.00011439, 0, 0.000228781,
887 0.00011439, 0, 0.000343171, 0.00011439, 0.000228781,
888 0.00011439, 0.00011439, 0, 0.00011439, 0,
889 0.000228781, 0.000343171, 0, 0.00011439, 0.00011439,
890 0, 0.00011439, 0.000343171, 0.00011439, 0.00011439,
891 0, 0.00011439, 0.00011439, 0, 0,
892 0, 0, 0.000343171, 0, 0,
893 0, 0, 0, 0.00011439, 0,
894 0, 0, 0.00011439, 0.00011439, 0,
895 0, 0.00011439, 0, 0.00011439, 0,
896 0.000228781, 0, 0, 0.00011439, 0.00011439,
897 0, 0, 0.000228781, 0, 0
33a848f6 898 },
899 {
32333a5e 900 0, 0, 0, 0.00010144, 0,
901 0.00010144, 0.00010144, 0, 0, 0.00010144,
902 0.000202881, 0, 0.00010144, 0.00010144, 0.000202881,
903 0.00010144, 0, 0.000304321, 0.000405762, 0.00010144,
904 0, 0.000202881, 0.000202881, 0.000507202, 0.000304321,
905 0, 0.000202881, 0.000405762, 0.000811524, 0.000507202,
906 0.00142017, 0.00152161, 0.00273889, 0.00223169, 0.00334753,
907 0.00436194, 0.00547778, 0.00618787, 0.00821668, 0.0095354,
908 0.0095354, 0.0118685, 0.0121729, 0.0121729, 0.014506,
909 0.0142017, 0.0165348, 0.0185636, 0.0177521, 0.0202881,
910 0.0185636, 0.0192737, 0.0208967, 0.0200852, 0.0222155,
911 0.0185636, 0.0181578, 0.0194766, 0.0207953, 0.0201867,
912 0.0180564, 0.017955, 0.017955, 0.0184622, 0.0191722,
913 0.0174478, 0.0174478, 0.017955, 0.0157233, 0.0155204,
914 0.017955, 0.0141002, 0.0164334, 0.0155204, 0.017042,
915 0.0146074, 0.0131873, 0.0141002, 0.0104484, 0.0129844,
916 0.0121729, 0.0114628, 0.0113613, 0.00963684, 0.00872388,
917 0.0108541, 0.00841956, 0.00892676, 0.00811524, 0.00831812,
918 0.00507202, 0.00811524, 0.00730371, 0.00689795, 0.00689795,
919 0.00730371, 0.00497058, 0.00568067, 0.00466626, 0.00486914,
920 0.00466626, 0.00405762, 0.00507202, 0.00355042, 0.00415906,
921 0.00284033, 0.00385474, 0.00334753, 0.00334753, 0.00243457,
922 0.00334753, 0.00344898, 0.00385474, 0.00182593, 0.00243457,
923 0.00223169, 0.00263745, 0.00273889, 0.00213025, 0.00213025,
924 0.00243457, 0.00152161, 0.00152161, 0.00162305, 0.00172449,
925 0.00111584, 0.00142017, 0.00152161, 0.00111584, 0.00162305,
926 0.00121729, 0.00111584, 0.000507202, 0.00121729, 0.000710083,
927 0.0010144, 0.000405762, 0.00121729, 0.000507202, 0.00131873,
928 0.0010144, 0.000608643, 0.000710083, 0.000507202, 0.000710083,
929 0.000608643, 0.000507202, 0.000507202, 0.000304321, 0.000608643,
930 0.000912964, 0.000608643, 0.000710083, 0.000405762, 0.000507202,
931 0.000811524, 0.000912964, 0.000507202, 0.000304321, 0.000304321,
932 0.000405762, 0.000507202, 0.000811524, 0.000405762, 0.000507202,
933 0.000608643, 0.00010144, 0.000405762, 0.000202881, 0.000304321,
934 0.000608643, 0.000608643, 0.00010144, 0.000507202, 0.00010144,
935 0.000304321, 0.000507202, 0.000304321, 0.000202881, 0.000304321,
936 0.000405762, 0.000304321, 0.000304321, 0.000405762, 0.000202881,
937 0.00010144, 0.000507202, 0, 0.000405762, 0,
938 0.000507202, 0.000405762, 0.000304321, 0.000405762, 0,
939 0.000507202, 0.000304321, 0.000405762, 0.000304321, 0.000507202,
940 0.000202881, 0.00010144, 0.000202881, 0, 0.000405762,
941 0.000202881, 0.000202881, 0.000405762, 0, 0.000202881,
942 0.000202881, 0, 0.000202881, 0.000202881, 0.000202881,
943 0.00010144, 0, 0.00010144, 0.000202881, 0.000202881,
944 0.000202881, 0.00010144, 0.000202881, 0, 0,
945 0.00010144, 0.00010144, 0.00010144, 0.00010144, 0,
33a848f6 946 0, 0, 0, 0, 0,
32333a5e 947 0, 0.00010144, 0, 0, 0,
948 0.00010144, 0.000202881, 0.00010144, 0.00010144, 0,
949 0, 0.00010144, 0, 0.00010144, 0.00010144
33a848f6 950 },
951 {
952 0, 0, 0, 0, 0,
32333a5e 953 0, 0, 0, 0.000316056, 0.000105352,
954 0.000210704, 0.000210704, 0.000316056, 0, 0,
955 0, 0.000105352, 0.000210704, 0.000105352, 0.000105352,
956 0.000210704, 0.000105352, 0, 0.000316056, 0.000105352,
957 0.000632111, 0.000210704, 0.000526759, 0.000737463, 0.00115887,
958 0.00105352, 0.00115887, 0.00147493, 0.00273915, 0.00337126,
959 0.0054783, 0.00495154, 0.00579435, 0.00758534, 0.00790139,
960 0.00842815, 0.00937632, 0.0106405, 0.013169, 0.0144332,
961 0.0142225, 0.015276, 0.0168563, 0.0150653, 0.0172777,
962 0.0184366, 0.0175938, 0.0189633, 0.0200169, 0.0198062,
963 0.0185419, 0.0187526, 0.0189633, 0.018858, 0.0198062,
964 0.0202276, 0.0166456, 0.0164349, 0.0198062, 0.0207543,
965 0.0154867, 0.0179098, 0.0189633, 0.0169617, 0.0147493,
966 0.0169617, 0.0155921, 0.0161188, 0.0163295, 0.0148546,
967 0.0154867, 0.0126422, 0.0112727, 0.0139064, 0.0120101,
968 0.0130636, 0.0117994, 0.0085335, 0.0110619, 0.0102191,
969 0.00842815, 0.00906026, 0.0104298, 0.0085335, 0.00874421,
970 0.00716393, 0.00769069, 0.00821745, 0.00716393, 0.00737463,
971 0.00684787, 0.00495154, 0.00453013, 0.00579435, 0.00558365,
972 0.00558365, 0.00368732, 0.00453013, 0.00400337, 0.00442478,
973 0.00474083, 0.00463548, 0.00431943, 0.00379267, 0.00400337,
974 0.00421408, 0.00273915, 0.0028445, 0.0028445, 0.0030552,
975 0.0026338, 0.0028445, 0.00273915, 0.0026338, 0.00294985,
976 0.00168563, 0.00189633, 0.00158028, 0.00105352, 0.00179098,
977 0.00316056, 0.00168563, 0.00126422, 0.00221239, 0.00242309,
978 0.00136957, 0.00147493, 0.00136957, 0.000948167, 0.00168563,
979 0.00126422, 0.00158028, 0.00126422, 0.00105352, 0.000948167,
980 0.00147493, 0.000421408, 0.000948167, 0.000948167, 0.00115887,
981 0.000842815, 0.000210704, 0.00126422, 0.000105352, 0.000632111,
982 0.000421408, 0.000948167, 0.000526759, 0.000526759, 0.000421408,
983 0.000526759, 0.000421408, 0.000737463, 0.000526759, 0.000632111,
984 0.000210704, 0.000737463, 0.000526759, 0.000210704, 0.000632111,
985 0.000526759, 0.000632111, 0.000526759, 0.000316056, 0.000316056,
986 0.000526759, 0.000526759, 0.000105352, 0.000948167, 0.000526759,
987 0.000210704, 0.000632111, 0.000632111, 0.000421408, 0.000421408,
988 0.000421408, 0, 0.000316056, 0.000210704, 0,
989 0.000421408, 0.000421408, 0.000316056, 0.000105352, 0.000105352,
990 0.000316056, 0.000105352, 0.000105352, 0.000210704, 0.000210704,
991 0.000105352, 0.000210704, 0.000105352, 0, 0.000842815,
992 0.000105352, 0, 0, 0.000421408, 0.000105352,
993 0.000105352, 0.000105352, 0, 0.000210704, 0,
994 0.000105352, 0.000210704, 0.000105352, 0.000105352, 0.000210704,
995 0.000105352, 0, 0.000105352, 0.000421408, 0.000210704,
996 0.000210704, 0.000210704, 0, 0, 0.000210704,
997 0, 0.000105352, 0.000210704, 0, 0.000210704,
998 0, 0, 0.000105352, 0.000105352, 0,
999 0, 0, 0.000316056, 0.000105352, 0.000105352,
1000 0, 0.000105352, 0, 0.000105352, 0,
1001 0.000105352, 0, 0.000105352, 0, 0.000210704
33a848f6 1002 },
1003 {
32333a5e 1004 0, 0, 0, 0, 0.000105552,
1005 0, 0.000105552, 0.000316656, 0.000105552, 0,
1006 0, 0.000105552, 0.000105552, 0.000105552, 0.000211104,
1007 0.000422208, 0.000422208, 0.000211104, 0.000105552, 0.000211104,
1008 0.000211104, 0.000105552, 0.000316656, 0.000211104, 0.00052776,
1009 0.00052776, 0.000316656, 0.000422208, 0.000422208, 0.00126662,
1010 0.000738864, 0.00147773, 0.00221659, 0.00147773, 0.00295546,
1011 0.00358877, 0.00369432, 0.00591091, 0.00749419, 0.00812751,
1012 0.00992189, 0.0102385, 0.011083, 0.0112941, 0.0164661,
1013 0.0150939, 0.017205, 0.0142495, 0.0137218, 0.0160439,
1014 0.0184716, 0.0170994, 0.016255, 0.0184716, 0.0196327,
1015 0.0180494, 0.0165717, 0.0189994, 0.0186827, 0.0168883,
1016 0.0200549, 0.0197382, 0.0183661, 0.0192105, 0.016255,
1017 0.0185772, 0.0183661, 0.018155, 0.0158328, 0.0166772,
1018 0.0151995, 0.014144, 0.0155161, 0.0145662, 0.0151995,
1019 0.0124551, 0.0132996, 0.0155161, 0.0128773, 0.0134051,
1020 0.0117163, 0.0113996, 0.0102385, 0.00876082, 0.0107663,
1021 0.0117163, 0.00802195, 0.00876082, 0.00918303, 0.0108719,
1022 0.0105552, 0.00907748, 0.0077053, 0.00643867, 0.00717754,
1023 0.00696643, 0.00675533, 0.00591091, 0.0050665, 0.00496095,
1024 0.00474984, 0.00538315, 0.00443319, 0.00453874, 0.00443319,
1025 0.00591091, 0.00390543, 0.00432763, 0.00422208, 0.00443319,
1026 0.00327211, 0.00401098, 0.00401098, 0.00327211, 0.00221659,
1027 0.00295546, 0.00369432, 0.00253325, 0.00274435, 0.00200549,
1028 0.00274435, 0.00232214, 0.00253325, 0.00179438, 0.00147773,
1029 0.00168883, 0.00158328, 0.00158328, 0.00158328, 0.00168883,
1030 0.00200549, 0.000949968, 0.00137218, 0.00105552, 0.00147773,
1031 0.000949968, 0.000844416, 0.000949968, 0.00105552, 0.000633312,
1032 0.000422208, 0.000633312, 0.000949968, 0.00116107, 0.000738864,
1033 0.000844416, 0.00052776, 0.000844416, 0.000422208, 0.00105552,
1034 0.000844416, 0.000633312, 0.000738864, 0.00105552, 0.000949968,
1035 0.00137218, 0.000633312, 0.00052776, 0.000633312, 0.000422208,
1036 0.000738864, 0.000316656, 0.000211104, 0.000633312, 0.000949968,
1037 0.000633312, 0.00052776, 0.000633312, 0.00052776, 0.00105552,
1038 0.000211104, 0.000422208, 0.000422208, 0.000422208, 0.000105552,
1039 0.000316656, 0.000105552, 0.000316656, 0.000211104, 0.000422208,
1040 0.000211104, 0, 0.000316656, 0, 0.000211104,
1041 0.000211104, 0.000422208, 0, 0.000211104, 0.000105552,
1042 0.000316656, 0.000316656, 0.000211104, 0.000211104, 0.000211104,
1043 0.000211104, 0.000211104, 0.00052776, 0.000211104, 0.000316656,
1044 0.000105552, 0.000316656, 0, 0, 0,
1045 0.000211104, 0.000316656, 0, 0, 0.000105552,
1046 0.000105552, 0.000105552, 0.000105552, 0.000422208, 0.000105552,
1047 0.000422208, 0.000105552, 0.000105552, 0.000105552, 0,
1048 0, 0, 0.000422208, 0, 0,
1049 0.000211104, 0, 0.000105552, 0.000316656, 0,
1050 0, 0, 0, 0.000211104, 0.000316656,
1051 0.000211104, 0, 0, 0, 0.000105552,
1052 0.000105552, 0.000422208, 0, 0.000211104, 0.000105552,
1053 0.000211104, 0, 0, 0, 0
33a848f6 1054 },
1055 {
32333a5e 1056 0, 0, 0, 0.000105686, 0.000105686,
1057 0, 0, 0, 0.000105686, 0.000105686,
1058 0.000105686, 0, 0.000317058, 0.000105686, 0.000105686,
1059 0.000105686, 0.000422744, 0, 0.000105686, 0.000105686,
1060 0.000105686, 0.000211372, 0.000211372, 0, 0.000105686,
1061 0, 0.000422744, 0.000105686, 0.000211372, 0.00052843,
1062 0.000951173, 0.00126823, 0.00158529, 0.0022194, 0.00264215,
1063 0.00338195, 0.00401606, 0.00570704, 0.00581272, 0.00686958,
1064 0.00951173, 0.0102515, 0.00887762, 0.0114141, 0.014796,
1065 0.0126823, 0.0121539, 0.0155358, 0.0158529, 0.0177552,
1066 0.0171211, 0.0150074, 0.0193405, 0.0175439, 0.020186,
1067 0.0198689, 0.0190235, 0.0174382, 0.0179666, 0.0171211,
1068 0.0158529, 0.0190235, 0.0189178, 0.018495, 0.0190235,
1069 0.0174382, 0.0180723, 0.0176495, 0.0154301, 0.0183893,
1070 0.0129994, 0.0131051, 0.0131051, 0.0132107, 0.0152188,
1071 0.0142676, 0.0152188, 0.0140562, 0.0128937, 0.0119425,
1072 0.0133164, 0.0138449, 0.0124709, 0.00919467, 0.0100402,
1073 0.011097, 0.00940605, 0.0105686, 0.00877193, 0.00803213,
1074 0.00845487, 0.0082435, 0.00718664, 0.0075037, 0.00729233,
1075 0.00581272, 0.00634115, 0.00623547, 0.00739801, 0.00560135,
1076 0.0060241, 0.00549567, 0.00517861, 0.0060241, 0.00412175,
1077 0.00475587, 0.00443881, 0.00285352, 0.00433312, 0.00348763,
1078 0.00401606, 0.00401606, 0.00285352, 0.00295921, 0.00317058,
1079 0.00454449, 0.00359332, 0.00359332, 0.00253646, 0.00327626,
1080 0.00359332, 0.00264215, 0.00285352, 0.00348763, 0.0014796,
1081 0.00158529, 0.00211372, 0.00200803, 0.00179666, 0.00232509,
1082 0.00253646, 0.00190235, 0.000739801, 0.0022194, 0.00158529,
1083 0.00179666, 0.00200803, 0.00169097, 0.000739801, 0.00126823,
1084 0.000845487, 0.000845487, 0.0014796, 0.00052843, 0.00105686,
1085 0.000739801, 0.00116254, 0.000739801, 0.00126823, 0.000951173,
1086 0.000634115, 0.000317058, 0.00052843, 0.000951173, 0.000951173,
1087 0.000739801, 0.000739801, 0.000951173, 0.000845487, 0.000951173,
1088 0.00052843, 0.000422744, 0.000739801, 0.000422744, 0.000739801,
1089 0.000105686, 0.000317058, 0.00052843, 0.00052843, 0.00052843,
1090 0.000845487, 0.000105686, 0.000317058, 0.000634115, 0.00052843,
1091 0.000317058, 0.000634115, 0.000317058, 0.000211372, 0.000422744,
1092 0.000211372, 0.000422744, 0.000105686, 0.000422744, 0.000105686,
1093 0.000105686, 0.000105686, 0.000105686, 0.00052843, 0.000211372,
1094 0.000317058, 0.00052843, 0.000105686, 0.000422744, 0.000211372,
1095 0.000105686, 0, 0.000211372, 0.000211372, 0.000211372,
1096 0.000634115, 0.000105686, 0, 0.000211372, 0.000317058,
1097 0.000105686, 0.000211372, 0.000211372, 0.000211372, 0.000317058,
1098 0, 0, 0, 0.00052843, 0.000105686,
1099 0.000105686, 0.000105686, 0.000211372, 0.000105686, 0.000211372,
1100 0, 0.000211372, 0.000105686, 0.000105686, 0.000211372,
1101 0, 0.000105686, 0, 0, 0.000317058,
1102 0, 0.000105686, 0.000211372, 0.000105686, 0,
1103 0, 0, 0.000422744, 0, 0,
1104 0, 0, 0.000105686, 0, 0,
1105 0, 0.000105686, 0, 0, 0.000105686
33a848f6 1106 }
1107 };
1108
1109 // Time Bin of Max. Cluster Distributions for e and pi
1110 fNTbins=20;
1111 fTBinSize=1;
1112
1642e650 1113 Int_t timBin[kNo_Mom][kNo_TBins]={
33a848f6 1114 {
1115 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1116 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1117 },
1118 {
1119 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1120 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1121 },
1122 {
1123 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1124 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1125 },
1126 {
1127 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1128 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1129 },
1130 {
1131 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1132 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1133 },
1134 {
1135 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1136 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1137 },
1138 {
1139 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
1140 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
1141 }
1142 };
1143
1642e650 1144 Double_t probPIT[kNo_Mom][kNo_TBins]={
33a848f6 1145 {
32333a5e 1146 0, 0.252938, 0.0801282, 0.0559562, 0.0532853 ,
1147 0.0491453, 0.0436699, 0.0431357, 0.0455395, 0.0422009 ,
1148 0.045406, 0.0432692, 0.0396635, 0.0419338, 0.0404647 ,
1149 0.0409989, 0.0389957, 0.0432692, 0, 0
33a848f6 1150 },
1151 {
32333a5e 1152 0, 0.265805, 0.076788, 0.057152, 0.0491699 ,
1153 0.041507, 0.0451788, 0.0506066, 0.0458174, 0.0439017 ,
1154 0.0431034, 0.0434227, 0.0431034, 0.0426245, 0.0391124 ,
1155 0.0408685, 0.0322478, 0.0395913, 0, 0
33a848f6 1156 },
1157 {
32333a5e 1158 0, 0.278986, 0.0757246, 0.0601449, 0.0496377 ,
1159 0.044686, 0.0422705, 0.0433575, 0.0388889, 0.0396135 ,
1160 0.042029, 0.0402174, 0.0427536, 0.0433575, 0.0376812 ,
1161 0.0386473, 0.0379227, 0.0440821, 0, 0
33a848f6 1162 },
1163 {
32333a5e 1164 0, 0.269746, 0.0749882, 0.0590033, 0.0465444 ,
1165 0.0470146, 0.0429008, 0.0429008, 0.0429008, 0.0440762 ,
1166 0.0394922, 0.0419605, 0.0446638, 0.0413728, 0.0400799 ,
1167 0.0393747, 0.0393747, 0.043606, 0, 0
33a848f6 1168 },
1169 {
32333a5e 1170 0, 0.269387, 0.0710154, 0.0571979, 0.0470223 ,
1171 0.0498072, 0.0501285, 0.0433805, 0.0460583, 0.0426307 ,
1172 0.043916, 0.0443445, 0.0435947, 0.0392031, 0.0344901 ,
1173 0.0396315, 0.0357755, 0.0424165, 0, 0
33a848f6 1174 },
1175 {
32333a5e 1176 0, 0.268377, 0.0723532, 0.0542071, 0.0520111 ,
1177 0.0447295, 0.0455386, 0.0420712, 0.0473879, 0.0453074 ,
1178 0.0420712, 0.0438049, 0.0420712, 0.0449607, 0.0394129 ,
1179 0.0372168, 0.0340962, 0.0443828, 0, 0
33a848f6 1180 },
1181 {
32333a5e 1182 0, 0.280106, 0.0703704, 0.0536508, 0.0493122 ,
1183 0.042328, 0.0434921, 0.0385185, 0.0421164, 0.0437037 ,
1184 0.0393651, 0.0404233, 0.0442328, 0.047619, 0.0390476 ,
1185 0.0386243, 0.0409524, 0.0461376, 0, 0
33a848f6 1186 }
1187 };
1188
1642e650 1189 Double_t probELT[kNo_Mom][kNo_TBins]={
33a848f6 1190 {
32333a5e 1191 0, 0.190173, 0.0727839, 0.0512603, 0.0444633,
1192 0.0446049, 0.0404984, 0.0346927, 0.0376664, 0.0379496,
1193 0.0407816, 0.0397904, 0.0421977, 0.0468706, 0.049561,
1194 0.0560748, 0.068819, 0.101813, 0, 0
33a848f6 1195 },
1196 {
32333a5e 1197 0, 0.195321, 0.0639535, 0.0449889, 0.0412514,
1198 0.040144, 0.0390365, 0.0398671, 0.0366833, 0.0368217,
1199 0.0336379, 0.0408361, 0.0460963, 0.053433, 0.0510797,
1200 0.0617386, 0.0661683, 0.108942, 0, 0
33a848f6 1201 },
1202 {
32333a5e 1203 0, 0.175017, 0.0619995, 0.0464425, 0.039579,
1204 0.0377488, 0.0386639, 0.0359186, 0.0412949, 0.0367193,
1205 0.0385495, 0.0431251, 0.0419812, 0.0465569, 0.0553649,
1206 0.0630291, 0.0720659, 0.125944, 0, 0
33a848f6 1207 },
1208 {
32333a5e 1209 0, 0.160783, 0.0604585, 0.0420978, 0.0420978,
1210 0.035707, 0.0361128, 0.0348955, 0.0390546, 0.0412863,
1211 0.0404747, 0.0438223, 0.0463583, 0.0521404, 0.0569081,
1212 0.0619801, 0.0722256, 0.133597, 0, 0
33a848f6 1213 },
1214 {
32333a5e 1215 0, 0.159292, 0.054783, 0.0401391, 0.0371892,
1216 0.0404551, 0.0385588, 0.0414033, 0.0363464, 0.0360303,
1217 0.0441424, 0.0442478, 0.0510957, 0.0517278, 0.0552044,
1218 0.0694269, 0.0711125, 0.128845, 0, 0
33a848f6 1219 },
1220 {
32333a5e 1221 0, 0.15495, 0.0529871, 0.040532, 0.0393709,
1222 0.037471, 0.0407431, 0.0376821, 0.0385265, 0.0403209,
1223 0.0401098, 0.0397931, 0.0525649, 0.0545704, 0.0575259,
1224 0.0670255, 0.0754697, 0.130357, 0, 0
33a848f6 1225 },
1226 {
32333a5e 1227 0, 0.157366, 0.0504122, 0.0432255, 0.0394208,
1228 0.0384697, 0.0374128, 0.0374128, 0.0394208, 0.0370958,
1229 0.0430142, 0.0455506, 0.0467132, 0.0526316, 0.0583386,
1230 0.0616149, 0.0783133, 0.133587, 0, 0
33a848f6 1231 }
1232 };
1233
1234 for(Int_t ip=0; ip<fNMom; ip++) {
1642e650 1235 fTrackMomentum[ip]= trackMomentum[ip];
33a848f6 1236 for(Int_t ie=0; ie<fNEbins; ie++) {
1642e650 1237 fEnergyLoss[ip][ie]=energyLoss[ip][ie];
1238 fProbPI[ip][ie]=probPI[ip][ie];
1239 fProbEL[ip][ie]=probEL[ip][ie];
33a848f6 1240 }
1241 }
1242
1243 for(Int_t ip=0; ip<fNMom; ip++) {
1642e650 1244 fTrackMomentum[ip]=trackMomentum[ip] ;
33a848f6 1245 for(Int_t it=0; it<fNTbins; it++) {
1642e650 1246 fTimBin[ip][it]=timBin[ip][it];
1247 fProbPIT[ip][it]=probPIT[ip][it];
1248 fProbELT[ip][it]=probELT[ip][it];
33a848f6 1249 }
1250 }
1251
1252}
1253
1254
1255/*
1256void AliTRDprobdist::FillData()
1257{
1258 //
1259 // Energy loss Distributions for e and pi
1260 Double_t temp;
1261 Double_t fADCNorm=1.;
1262 ifstream filein("pp.dat");
1263 filein >> fNMom;
1264 for(Int_t ip=0; ip<fNMom; ip++) {
1265 filein >> fTrackMomentum[ip] >> fNEbins;
1266 for(Int_t ie=0; ie<fNEbins; ie++) {
1267 filein >> fEnergyLoss[ip][ie]>> fProbPI[ip][ie]>> fProbEL[ip][ie];
1268 fEnergyLoss[ip][ie]=fADCNorm*fEnergyLoss[ip][ie];
1269 }
1270 }
1271 fEnBinSize=fEnergyLoss[0][1]-fEnergyLoss[0][0];
1272
1273 // Time Bin of Max. Cluster Distributions for e and pi
1274 ifstream fileint("tt.dat");
1275 fileint >> fNMom;
1276 for(Int_t ip=0; ip<fNMom; ip++) {
1277 fileint >> fTrackMomentum[ip] >> fNTbins;
1278 for(Int_t it=0; it<fNTbins; it++) {
1279 fileint >> fTimBin[ip][it]>> fProbPIT[ip][it]>> fProbELT[ip][it];
1280 }
1281 }
1282 fTBinSize=1;
1283}
1284
1285*/
1286