1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 //-----------------------------------------------------------------
17 // Class for dE/dx and Time Bin of Max. Cluster for Electrons and
19 // It is instantiated in class AliTRDpidESD for particle identification
21 // Prashant Shukla <shukla@pi0.physi.uni-heidelberg.de>
22 //-----------------------------------------------------------------
24 #include "AliTRDprobdist.h"
26 ClassImp(AliTRDprobdist)
28 //_________________________________________________________________________
29 AliTRDprobdist::AliTRDprobdist(Int_t multiplicity)
32 // The main constructor
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();
41 //_________________________________________________________________________
42 Double_t AliTRDprobdist::GetMeanPI(Int_t ip) const
44 // Gets mean of de/dx dist. of pi
47 for(Int_t ie=0; ie<fNEbins; ie++) {
48 integral+=fEnergyLoss[ip][ie]*fProbPI[ip][ie];
49 norm+=fProbPI[ip][ie];
54 //_________________________________________________________________________
55 Double_t AliTRDprobdist::GetMeanEL(Int_t ip) const
58 // Gets mean of de/dx dist. of e
61 for(Int_t ie=0; ie<fNEbins; ie++) {
62 integral+=fEnergyLoss[ip][ie]*fProbEL[ip][ie];
63 norm+=fProbEL[ip][ie];
68 //_________________________________________________________________________
69 Double_t AliTRDprobdist::GetNormalizationPI(Int_t ip) const
72 // Gets Normalization of de/dx dist. of pi
74 for(Int_t ie=0; ie<fNEbins; ie++) {
75 integral+=fProbPI[ip][ie];
80 //_________________________________________________________________________
81 Double_t AliTRDprobdist::GetNormalizationEL(Int_t ip) const
84 // Gets Normalization of de/dx dist. of e
86 for(Int_t ie=0; ie<fNEbins; ie++) {
87 integral+=fProbEL[ip][ie];
92 //_________________________________________________________________________
93 Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx) const
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
99 Double_t probability = 1.0;
100 Int_t iEnBin= ((Int_t) (dedx/fEnBinSize));
101 if(iEnBin > fNEbins-1) iEnBin = fNEbins-1;
103 if(k==0){ // electron
104 if(mom<=fTrackMomentum[0]) probability = fProbEL[0][iEnBin];
105 if(mom>=fTrackMomentum[fNMom-1]) probability = fProbEL[fNMom-1][iEnBin];
108 if(mom<=fTrackMomentum[0]) probability = fProbPI[0][iEnBin];
109 if(mom>=fTrackMomentum[fNMom-1]) probability = fProbPI[fNMom-1][iEnBin];
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
117 probability= fProbEL[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
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
126 probability= fProbPI[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
133 //_________________________________________________________________________
134 Double_t AliTRDprobdist::GetProbabilityT(Int_t k, Double_t mom, Int_t timbin) const
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
140 Double_t probabilityT = 1.0;
141 if(timbin<=0) return 0.;
144 if(k==0){ // electron
145 if(mom<=fTrackMomentum[0]) probabilityT = fProbELT[0][iTBin];
146 if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbELT[fNMom-1][iTBin];
149 if(mom<=fTrackMomentum[0]) probabilityT = fProbPIT[0][iTBin];
150 if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbPIT[fNMom-1][iTBin];
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
158 probabilityT= fProbELT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
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
167 probabilityT= fProbPIT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
174 void AliTRDprobdist::FillData()
177 // Energy loss Distributions for e and pi
179 Double_t trackMomentum[kNo_Mom]= {1., 1.5, 2., 3., 4., 5., 6.};
183 Double_t energyLoss[kNo_Mom][kNo_EnBins]={
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,
204 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
207 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
208 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
209 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
210 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
211 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
212 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
213 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
214 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
215 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
216 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
217 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
218 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
219 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
220 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
221 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
222 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
223 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
224 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
225 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
226 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
229 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
230 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
231 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
232 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
233 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
234 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
235 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
236 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
237 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
238 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
239 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
240 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
241 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
242 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
243 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
244 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
245 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
246 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
247 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
248 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
251 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
252 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
253 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
254 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
255 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
256 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
257 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
258 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
259 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
260 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
261 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
262 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
263 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
264 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
265 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
266 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
267 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
268 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
269 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
270 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
273 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
274 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
275 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
276 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
277 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
278 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
279 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
280 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
281 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
282 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
283 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
284 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
285 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
286 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
287 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
288 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
289 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
290 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
291 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
292 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
295 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
296 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
297 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
298 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
299 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
300 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
301 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
302 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
303 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
304 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
305 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
306 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
307 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
308 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
309 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
310 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
311 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
312 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
313 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
314 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
317 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
318 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
319 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
320 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
321 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
322 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
323 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
324 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
325 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
326 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
327 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
328 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
329 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
330 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
331 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
332 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
333 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
334 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
335 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
336 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990
340 Double_t probPI[kNo_Mom][kNo_EnBins]={
344 0, 0, 0, 0.000277546, 0.000693866,
345 0.000832639, 0.00430197, 0.00402442, 0.00679989, 0.010408,
346 0.0138773, 0.0222037, 0.0255343, 0.0298362, 0.0359423,
347 0.0414932, 0.0452401, 0.0474604, 0.0392728, 0.0481543,
348 0.0431585, 0.047738, 0.0421871, 0.0370525, 0.0344158,
349 0.0338607, 0.0303913, 0.0341382, 0.025118, 0.0244241,
350 0.0234527, 0.020816, 0.017763, 0.0184568, 0.0162365,
351 0.0154038, 0.0112406, 0.011657, 0.0109631, 0.00957535,
352 0.00818762, 0.00679989, 0.00638357, 0.00582848, 0.00679989,
353 0.00527338, 0.00499584, 0.0036081, 0.00346933, 0.00263669,
354 0.0036081, 0.00305301, 0.00166528, 0.00222037, 0.00194283,
355 0.00180405, 0.0020816, 0.00138773, 0.00166528, 0.00166528,
356 0.000555093, 0.00111019, 0.00111019, 0.0020816, 0.000971413,
357 0.000832639, 0.00152651, 0.000971413, 0.000832639, 0.00111019,
358 0.000971413, 0.000832639, 0.00041632, 0.000138773, 0.000555093,
359 0.000555093, 0.00041632, 0.00041632, 0.000555093, 0.000555093,
360 0.000693866, 0.00041632, 0.000138773, 0.000277546, 0.00041632,
361 0.000555093, 0.000832639, 0, 0.00041632, 0.000693866,
362 0.000138773, 0.000555093, 0.000138773, 0.000277546, 0,
363 0.000277546, 0.000277546, 0.000138773, 0.000138773, 0.000138773,
364 0.000277546, 0.000138773, 0.000277546, 0, 0.000277546,
365 0, 0.000138773, 0.000555093, 0.000277546, 0.000277546,
366 0.00041632, 0.00041632, 0, 0.000277546, 0.000138773,
367 0.000138773, 0, 0.000138773, 0, 0,
368 0.00041632, 0.000138773, 0.000277546, 0.000277546, 0,
369 0.000277546, 0.000138773, 0.00041632, 0.000138773, 0,
370 0.000277546, 0.000138773, 0, 0.000277546, 0.00041632,
371 0, 0.000277546, 0.000138773, 0.00041632, 0.000138773,
372 0.000138773, 0, 0.000138773, 0, 0,
373 0, 0, 0.000138773, 0.000277546, 0.000277546,
374 0.000277546, 0, 0.000277546, 0.000138773, 0,
375 0.000277546, 0.000138773, 0, 0, 0.000138773,
376 0.000277546, 0, 0.000138773, 0, 0,
377 0.000138773, 0.000138773, 0, 0.000138773, 0.000138773,
378 0.000277546, 0, 0.000277546, 0.000138773, 0,
379 0, 0.000277546, 0, 0, 0.000138773,
380 0, 0, 0, 0.000138773, 0,
381 0, 0, 0.000138773, 0.000277546, 0.000138773
386 0, 0, 0.00012273, 0.000245459, 0,
387 0.000859107, 0.0012273, 0.00196367, 0.00454099, 0.00736377,
388 0.0100638, 0.0160776, 0.0177958, 0.0255277, 0.0300687,
389 0.0363279, 0.037678, 0.0385371, 0.0414826, 0.0393962,
390 0.0432008, 0.0411144, 0.0385371, 0.0380461, 0.0366961,
391 0.0352234, 0.037187, 0.0381689, 0.0281051, 0.0298233,
392 0.0265096, 0.0225822, 0.0228277, 0.0204958, 0.0196367,
393 0.0160776, 0.0142366, 0.014973, 0.0146048, 0.0132548,
394 0.0112911, 0.0109229, 0.00834561, 0.00895925, 0.00785469,
395 0.00589102, 0.00552283, 0.0041728, 0.00405007, 0.00380461,
396 0.00220913, 0.00220913, 0.00306824, 0.00245459, 0.00220913,
397 0.0020864, 0.00196367, 0.00233186, 0.00171821, 0.00135002,
398 0.0020864, 0.00171821, 0.000736377, 0.000981836, 0.000981836,
399 0.00147275, 0.000981836, 0.000736377, 0.00147275, 0.00135002,
400 0.0012273, 0.00147275, 0.000613648, 0.000859107, 0.000368189,
401 0.00012273, 0.000245459, 0.000981836, 0.000736377, 0.000613648,
402 0.000613648, 0.000490918, 0.000736377, 0.000368189, 0.000368189,
403 0.000490918, 0.000368189, 0.000368189, 0.000245459, 0.000613648,
404 0.000368189, 0.000245459, 0.000368189, 0.000245459, 0.000490918,
405 0.000245459, 0.000245459, 0.000368189, 0.000613648, 0.00012273,
406 0.000613648, 0.000368189, 0.000245459, 0.00012273, 0,
407 0, 0.000368189, 0.000368189, 0.000490918, 0,
408 0.000245459, 0.00012273, 0.000245459, 0.00012273, 0.00012273,
409 0.000245459, 0.000490918, 0.000490918, 0, 0,
410 0, 0.000368189, 0.00012273, 0.00012273, 0,
411 0.00012273, 0, 0.000245459, 0, 0.000245459,
412 0, 0.000245459, 0.00012273, 0.000490918, 0.000368189,
413 0.000490918, 0.000245459, 0, 0.00012273, 0,
414 0.000490918, 0.000245459, 0.000245459, 0.00012273, 0.00012273,
415 0.000245459, 0.000245459, 0.000368189, 0.00012273, 0,
416 0, 0, 0, 0.000245459, 0.00012273,
417 0.00012273, 0.00012273, 0.000490918, 0, 0.00012273,
418 0, 0.00012273, 0.000368189, 0, 0,
419 0.00012273, 0, 0, 0.000368189, 0,
420 0.000245459, 0, 0.00012273, 0, 0,
421 0.00012273, 0, 0, 0.000368189, 0,
422 0, 0.00012273, 0.00012273, 0.000245459, 0,
423 0, 0.000245459, 0, 0, 0
427 0, 0, 0, 0.000117123, 0,
428 0, 0, 0, 0, 0.000117123,
429 0.000234247, 0.000585617, 0.00117123, 0.00304521, 0.00527056,
430 0.00866714, 0.0124151, 0.0172171, 0.0194425, 0.0224877,
431 0.0311548, 0.036074, 0.0385336, 0.0377138, 0.0429843,
432 0.0405247, 0.0441555, 0.0402905, 0.0392364, 0.040759,
433 0.037011, 0.0340829, 0.0361911, 0.0327946, 0.0346685,
434 0.0264699, 0.0264699, 0.0243617, 0.0207309, 0.0174514,
435 0.0196767, 0.0190911, 0.0154603, 0.0148747, 0.0127665,
436 0.0127665, 0.0115952, 0.00925275, 0.00784727, 0.00714453,
437 0.00866714, 0.00632467, 0.00538768, 0.00538768, 0.00386507,
438 0.0039822, 0.00339658, 0.00468494, 0.00374795, 0.0019911,
439 0.00175685, 0.00175685, 0.000936988, 0.00128836, 0.000936988,
440 0.000819864, 0.0015226, 0.00105411, 0.000936988, 0.00105411,
441 0.000936988, 0.000936988, 0.00128836, 0.000702741, 0.000585617,
442 0.00117123, 0.000117123, 0.000702741, 0.000234247, 0.000819864,
443 0.000702741, 0.000819864, 0.00035137, 0.000585617, 0.000468494,
444 0.00035137, 0.00105411, 0.000585617, 0.000117123, 0.000585617,
445 0.000468494, 0.000117123, 0.000702741, 0.000585617, 0.00035137,
446 0, 0.000585617, 0.00035137, 0.000117123, 0.000234247,
447 0.00035137, 0.000585617, 0.00035137, 0.00035137, 0.000234247,
448 0.00035137, 0.000117123, 0.000117123, 0.000234247, 0.000468494,
449 0.00035137, 0.000117123, 0.000117123, 0.000585617, 0.000468494,
450 0.000234247, 0.000234247, 0.00035137, 0, 0.000117123,
451 0, 0.000117123, 0.000234247, 0.00035137, 0.00035137,
452 0, 0, 0.000234247, 0.000117123, 0.000234247,
453 0.000117123, 0, 0, 0, 0,
454 0.000117123, 0.000234247, 0.000234247, 0.000117123, 0.000468494,
455 0.000468494, 0.000117123, 0, 0, 0.000117123,
456 0.000117123, 0.000234247, 0, 0.000234247, 0,
457 0.000117123, 0, 0.000117123, 0.000234247, 0.000117123,
458 0.000234247, 0, 0, 0.000117123, 0.00035137,
459 0.000117123, 0.000117123, 0.000234247, 0.000117123, 0.00035137,
460 0.000234247, 0.000234247, 0.000117123, 0.000117123, 0,
461 0.000234247, 0, 0.000234247, 0, 0,
462 0.000117123, 0.00035137, 0.000117123, 0, 0.00035137,
463 0, 0, 0.000117123, 0.000234247, 0.000117123,
464 0.000117123, 0, 0, 0.000234247, 0,
465 0.000117123, 0, 0.000117123, 0, 0.000117123
470 0, 0.000329598, 0, 0, 0.000109866,
471 0.000109866, 0.000439464, 0.000769062, 0.000878928, 0.00252692,
472 0.00428477, 0.00681169, 0.0105471, 0.0134036, 0.0202153,
473 0.0220831, 0.0313118, 0.0288947, 0.0330697, 0.0377939,
474 0.0405405, 0.0396616, 0.0368051, 0.0416392, 0.0382334,
475 0.0379038, 0.0392221, 0.035267, 0.0332894, 0.0310921,
476 0.0314217, 0.0295539, 0.0274665, 0.0277961, 0.0238409,
477 0.0228521, 0.0203252, 0.0177983, 0.0175786, 0.0152714,
478 0.0126346, 0.0142826, 0.015601, 0.010657, 0.00911887,
479 0.00813008, 0.00758075, 0.00637223, 0.0054933, 0.0054933,
480 0.0045045, 0.00417491, 0.00538343, 0.00395517, 0.00351571,
481 0.00329598, 0.00285652, 0.00208745, 0.00252692, 0.00164799,
482 0.00153812, 0.00197759, 0.00197759, 0.00219732, 0.00164799,
483 0.000439464, 0.00120853, 0.00197759, 0.00120853, 0.00142826,
484 0.000988794, 0.000878928, 0.00054933, 0.000439464, 0.000329598,
485 0.000329598, 0.000659196, 0.000988794, 0.000878928, 0.000769062,
486 0.000329598, 0.000439464, 0.000109866, 0.000769062, 0.000329598,
487 0.000219732, 0.00054933, 0.000329598, 0.000439464, 0.000439464,
488 0.000988794, 0.000439464, 0.000219732, 0.000219732, 0.000439464,
489 0.000769062, 0, 0.000439464, 0.00054933, 0,
490 0.000769062, 0.000109866, 0.000109866, 0.000659196, 0.000219732,
491 0.000109866, 0.00054933, 0.000109866, 0.000329598, 0.000329598,
492 0.000219732, 0.000109866, 0.000439464, 0, 0.000109866,
493 0.000219732, 0.000329598, 0, 0.000109866, 0.000219732,
494 0, 0.000769062, 0.000329598, 0.000219732, 0.000439464,
495 0.000109866, 0.000219732, 0.000329598, 0, 0.000109866,
496 0.000109866, 0.000219732, 0.000219732, 0.000109866, 0.000109866,
497 0.000109866, 0, 0, 0.000109866, 0.000109866,
498 0, 0, 0.000109866, 0, 0,
499 0.000219732, 0, 0.000329598, 0.000219732, 0.000109866,
500 0.000109866, 0, 0, 0.000109866, 0,
501 0.000109866, 0.000109866, 0.000109866, 0, 0,
502 0, 0.000109866, 0.000329598, 0.000329598, 0,
503 0, 0, 0, 0, 0.000219732,
504 0.000219732, 0.000109866, 0.000109866, 0.000109866, 0,
505 0, 0.000219732, 0.000219732, 0, 0,
506 0.000109866, 0, 0, 0.000109866, 0.000109866,
507 0, 0, 0.000109866, 0.000109866, 0
511 0, 0, 0, 0, 0.000112007,
512 0, 0, 0.000112007, 0, 0.000112007,
513 0, 0.000112007, 0.000448029, 0.000448029, 0.00168011,
514 0.00235215, 0.00392025, 0.00728047, 0.0103047, 0.0141129,
515 0.0181452, 0.0197133, 0.0259857, 0.031362, 0.0323701,
516 0.0370744, 0.0374104, 0.0385305, 0.0393145, 0.0356183,
517 0.0375224, 0.0405466, 0.0378584, 0.0365143, 0.0348342,
518 0.0355063, 0.0293459, 0.0300179, 0.0256496, 0.0281138,
519 0.0239695, 0.0196013, 0.0228495, 0.0183692, 0.0194892,
520 0.0172491, 0.014897, 0.0126568, 0.0106407, 0.0112007,
521 0.0094086, 0.0108647, 0.00985663, 0.00638441, 0.00672043,
522 0.00481631, 0.00515233, 0.00515233, 0.00425627, 0.00504032,
523 0.00347222, 0.00347222, 0.00291219, 0.00280018, 0.00257616,
524 0.00179211, 0.0015681, 0.00123208, 0.00168011, 0.0015681,
525 0.000560036, 0.00112007, 0.00168011, 0.00212814, 0.00179211,
526 0.00123208, 0.00112007, 0.000560036, 0.00112007, 0.00078405,
527 0.00078405, 0.000448029, 0.00123208, 0.000896057, 0.00078405,
528 0.000448029, 0.000672043, 0.000112007, 0.00112007, 0.000224014,
529 0.000448029, 0.00100806, 0.000560036, 0.000448029, 0.000224014,
530 0.000896057, 0.000112007, 0.000560036, 0.000448029, 0.000336022,
531 0.000448029, 0.000448029, 0.000672043, 0, 0.000560036,
532 0.000448029, 0.000336022, 0.000336022, 0.000336022, 0.000112007,
533 0.000560036, 0.000112007, 0.000112007, 0.000336022, 0.000224014,
534 0.000112007, 0.000560036, 0.000560036, 0.000336022, 0.000448029,
535 0.000560036, 0.000112007, 0, 0.000224014, 0.000336022,
536 0.000224014, 0.000224014, 0.000336022, 0.000224014, 0.000336022,
537 0.000224014, 0.000112007, 0.000224014, 0.000112007, 0.000224014,
538 0, 0.000224014, 0.000112007, 0.000112007, 0.000112007,
539 0.000112007, 0.000336022, 0, 0, 0.000224014,
540 0.000112007, 0.000112007, 0.000224014, 0.000112007, 0.000112007,
541 0.000112007, 0, 0.000112007, 0.000112007, 0.000224014,
542 0, 0, 0, 0, 0.000224014,
543 0, 0.000112007, 0.000112007, 0, 0.000112007,
544 0.000224014, 0.000112007, 0, 0, 0,
545 0.000224014, 0.000336022, 0, 0, 0.000112007,
546 0.000112007, 0, 0, 0.000112007, 0.000112007,
547 0, 0, 0, 0, 0.000112007,
548 0, 0.000224014, 0.000112007, 0, 0.000224014,
549 0.000112007, 0.000224014, 0, 0, 0.000112007
553 0, 0, 0, 0, 0.000107805,
554 0, 0, 0.000107805, 0, 0,
555 0, 0.000323415, 0.000323415, 0.000323415, 0.000862441,
556 0.0021561, 0.00301854, 0.00614489, 0.00689953, 0.0115351,
557 0.0141225, 0.0208064, 0.0242561, 0.029323, 0.0299698,
558 0.0344976, 0.0353601, 0.0406425, 0.0389176, 0.0392411,
559 0.0416128, 0.0398879, 0.0367615, 0.0369771, 0.0330962,
560 0.0305088, 0.0328806, 0.0320181, 0.0268435, 0.0271669,
561 0.0260888, 0.0230703, 0.0212376, 0.0194049, 0.017141,
562 0.01768, 0.0146615, 0.0167098, 0.0129366, 0.0118586,
563 0.0100259, 0.0102415, 0.00743855, 0.0085166, 0.00776197,
564 0.00873221, 0.00657611, 0.00582147, 0.00366537, 0.00485123,
565 0.00388098, 0.00344976, 0.00377318, 0.00291074, 0.00258732,
566 0.00344976, 0.00161708, 0.00172488, 0.00172488, 0.00194049,
567 0.00140147, 0.000970246, 0.00107805, 0.00161708, 0.00107805,
568 0.00150927, 0.000107805, 0.000754636, 0.000323415, 0.00107805,
569 0.000754636, 0.00043122, 0.000539025, 0.00021561, 0.000539025,
570 0.000539025, 0.00043122, 0.00043122, 0.000323415, 0.000539025,
571 0.000646831, 0.000754636, 0.000539025, 0.000539025, 0.000323415,
572 0.00043122, 0.000323415, 0.000539025, 0.000539025, 0.000754636,
573 0.000539025, 0.00021561, 0, 0.000754636, 0.00021561,
574 0.000107805, 0.00043122, 0.00021561, 0.000646831, 0.000323415,
575 0.000107805, 0.00021561, 0.000323415, 0.000107805, 0.00021561,
576 0.00021561, 0.000323415, 0.00021561, 0, 0.000862441,
577 0.00021561, 0.000323415, 0.000323415, 0.00021561, 0.000107805,
578 0.000323415, 0.000323415, 0.000107805, 0.000323415, 0,
579 0.000323415, 0, 0.00043122, 0.00021561, 0.00043122,
580 0.00021561, 0, 0, 0.00021561, 0.000107805,
581 0.00021561, 0.000323415, 0.000107805, 0, 0.000107805,
582 0, 0.000323415, 0.000107805, 0, 0,
583 0, 0.000107805, 0, 0.000107805, 0,
584 0.00021561, 0, 0.000107805, 0.00021561, 0.00021561,
585 0.00021561, 0.000323415, 0.000107805, 0, 0.000107805,
586 0, 0, 0.00021561, 0.000107805, 0,
587 0.000107805, 0.000107805, 0, 0.000107805, 0.000107805,
588 0, 0.00021561, 0, 0.000107805, 0.000107805,
589 0.000107805, 0, 0.000323415, 0.000107805, 0,
590 0, 0.000107805, 0, 0.00021561, 0,
591 0, 0.000107805, 0.00021561, 0.000107805, 0.00021561
595 0, 0, 0, 0, 0.000107043,
596 0, 0.000107043, 0.000107043, 0.000107043, 0.000107043,
597 0.000107043, 0.000214087, 0.00032113, 0.000214087, 0.000856348,
598 0.00107043, 0.00299722, 0.00406765, 0.00524513, 0.00941982,
599 0.0134875, 0.0160565, 0.0190537, 0.0205523, 0.0308285,
600 0.0354314, 0.0322201, 0.0331835, 0.0358596, 0.0400343,
601 0.0436737, 0.0416399, 0.037037, 0.041854, 0.0359666,
602 0.0331835, 0.0332905, 0.0324342, 0.0260116, 0.0264397,
603 0.027082, 0.0257975, 0.0219439, 0.0203383, 0.0214087,
604 0.0194819, 0.0169129, 0.0165917, 0.0149861, 0.0125241,
605 0.00931278, 0.0103832, 0.00995504, 0.00792122, 0.00663669,
606 0.00770713, 0.00770713, 0.00620852, 0.00449583, 0.00460287,
607 0.00428174, 0.00310426, 0.00363948, 0.00235496, 0.00224791,
608 0.00224791, 0.00256904, 0.002462, 0.00224791, 0.002462,
609 0.00181974, 0.00139156, 0.00149861, 0.00149861, 0.000963391,
610 0.00107043, 0.00107043, 0.000963391, 0.00128452, 0.000535217,
611 0.000642261, 0.00128452, 0.000963391, 0.000749304, 0.000642261,
612 0.000428174, 0.00032113, 0.000856348, 0.000963391, 0.000856348,
613 0.000642261, 0.000535217, 0.000428174, 0.00032113, 0.000963391,
614 0.000428174, 0.000428174, 0.00032113, 0.000107043, 0.000428174,
615 0.000535217, 0.00032113, 0.000428174, 0.000428174, 0.000535217,
616 0.000642261, 0.000214087, 0.000214087, 0, 0.00032113,
617 0.000428174, 0.000428174, 0.000214087, 0.000214087, 0.000214087,
618 0.000428174, 0.000214087, 0.000214087, 0.000107043, 0.000214087,
619 0.00032113, 0.000107043, 0.000214087, 0.00032113, 0,
620 0.000214087, 0.00032113, 0.000214087, 0.00032113, 0.000107043,
621 0.000214087, 0.00032113, 0.000428174, 0.00032113, 0.000107043,
622 0.000214087, 0.000214087, 0.000107043, 0.000214087, 0,
623 0.000107043, 0, 0, 0, 0.00032113,
624 0.00032113, 0, 0.000214087, 0.000107043, 0.000214087,
625 0.000214087, 0.000107043, 0.000214087, 0.00032113, 0.000214087,
626 0.000107043, 0.00032113, 0.000107043, 0.000214087, 0.000214087,
627 0, 0.000107043, 0.000107043, 0.000214087, 0.000107043,
628 0.000214087, 0.000107043, 0.000107043, 0, 0,
629 0, 0, 0, 0, 0.000107043,
630 0, 0.000107043, 0, 0.000107043, 0.000107043,
631 0, 0.000107043, 0.000107043, 0.000107043, 0,
632 0, 0.00032113, 0.000107043, 0.000107043, 0.000107043,
633 0, 0, 0.000107043, 0, 0
637 Double_t probEL[kNo_Mom][kNo_EnBins]={
643 0.000292654, 0, 0.000292654, 0.000292654, 0.000731636,
644 0.00117062, 0.00175593, 0.0016096, 0.00234124, 0.00599941,
645 0.00570676, 0.00717003, 0.00921861, 0.0121452, 0.0136084,
646 0.015218, 0.0155107, 0.0166813, 0.0179982, 0.0215101,
647 0.0247293, 0.0213638, 0.0206321, 0.023266, 0.0272169,
648 0.0278022, 0.0244366, 0.0257536, 0.0244366, 0.0272169,
649 0.023705, 0.0260462, 0.0270705, 0.0216564, 0.0257536,
650 0.023266, 0.0235587, 0.0218028, 0.0225344, 0.0200468,
651 0.0203395, 0.0188762, 0.0162423, 0.0184372, 0.0149254,
652 0.0155107, 0.015657, 0.016096, 0.0125841, 0.0105356,
653 0.0102429, 0.0115598, 0.0119988, 0.00892596, 0.00877963,
654 0.00760901, 0.00746269, 0.00599941, 0.00629207, 0.00731636,
655 0.00497512, 0.00702371, 0.00438982, 0.00599941, 0.00395083,
656 0.0048288, 0.00409716, 0.00380451, 0.00409716, 0.00395083,
657 0.00278022, 0.00292654, 0.00351185, 0.00307287, 0.00278022,
658 0.0032192, 0.00307287, 0.00102429, 0.00292654, 0.00131694,
659 0.00219491, 0.0016096, 0.00175593, 0.00117062, 0.00102429,
660 0.00117062, 0.000585309, 0.00190225, 0.0016096, 0.00102429,
661 0.00190225, 0.00117062, 0.00175593, 0.000438982, 0.000731636,
662 0.00146327, 0.00102429, 0.000585309, 0.000877963, 0.000877963,
663 0.000877963, 0.000585309, 0.000438982, 0.000585309, 0.000438982,
664 0.000292654, 0.000585309, 0.000438982, 0.000438982, 0.000292654,
665 0.000438982, 0.000877963, 0.000438982, 0.000877963, 0,
666 0.000146327, 0.000146327, 0.000438982, 0.000146327, 0.000585309,
667 0.000438982, 0.000292654, 0.000292654, 0.000292654, 0.000146327,
668 0.000292654, 0.000292654, 0.000292654, 0, 0.000438982,
669 0.000146327, 0.000146327, 0.000438982, 0, 0.000292654,
670 0.000292654, 0.000146327, 0, 0.000146327, 0.000438982,
671 0.000146327, 0.000292654, 0.000438982, 0, 0,
672 0.000585309, 0.000438982, 0, 0.000146327, 0,
673 0.000146327, 0.000585309, 0.000146327, 0.000438982, 0.000146327,
674 0, 0.000146327, 0.000146327, 0, 0,
675 0, 0, 0, 0.000292654, 0,
676 0.000146327, 0.000146327, 0, 0, 0.000292654,
678 0.000146327, 0, 0, 0.000146327, 0.000292654
685 0, 0.000131961, 0.000527844, 0.000395883, 0.000263922,
686 0.000923727, 0.00145157, 0.000923727, 0.0023753, 0.00329902,
687 0.00461863, 0.00633413, 0.00936923, 0.00897334, 0.0131961,
688 0.0108208, 0.0183426, 0.0129322, 0.0179467, 0.0164951,
689 0.0199261, 0.0197941, 0.0226973, 0.0182106, 0.0213777,
690 0.0207179, 0.0225653, 0.0249406, 0.0258643, 0.0215096,
691 0.020322, 0.0224334, 0.0225653, 0.0242808, 0.0207179,
692 0.0238849, 0.023753, 0.0204539, 0.0196622, 0.0204539,
693 0.0182106, 0.0187385, 0.0188704, 0.0158353, 0.0150435,
694 0.0160992, 0.0153075, 0.0147796, 0.0163632, 0.0137239,
695 0.0129322, 0.0113486, 0.0109528, 0.00884138, 0.0112167,
696 0.0104249, 0.0091053, 0.0114806, 0.00884138, 0.00923727,
697 0.00620216, 0.00712589, 0.00765373, 0.00752177, 0.00514648,
698 0.00461863, 0.00527844, 0.00580628, 0.0060702, 0.00343098,
699 0.00461863, 0.00409079, 0.00514648, 0.00461863, 0.00343098,
700 0.00409079, 0.00356295, 0.00422275, 0.00356295, 0.00277118,
701 0.00224334, 0.0030351, 0.00263922, 0.00131961, 0.00224334,
702 0.00224334, 0.00250726, 0.00145157, 0.00171549, 0.00171549,
703 0.000791766, 0.00224334, 0.00224334, 0.00145157, 0.00131961,
704 0.00171549, 0.000659805, 0.00118765, 0.00118765, 0.00105569,
705 0.00118765, 0.000791766, 0.00118765, 0.000395883, 0.00158353,
706 0.000395883, 0.000923727, 0.000791766, 0.000659805, 0.000791766,
707 0.000263922, 0.000659805, 0.000791766, 0.00105569, 0.000395883,
708 0.000527844, 0.000263922, 0.000527844, 0.000659805, 0.000263922,
709 0.000527844, 0.000263922, 0.000395883, 0.000527844, 0.000131961,
710 0.000395883, 0, 0.000395883, 0.000131961, 0.000527844,
711 0.000527844, 0.000263922, 0.000659805, 0.000131961, 0.000263922,
712 0.000527844, 0.000131961, 0.000131961, 0.000395883, 0,
713 0.000527844, 0.000263922, 0.000263922, 0.000131961, 0.000131961,
714 0.000395883, 0.000263922, 0.000131961, 0.000131961, 0.000131961,
715 0.000131961, 0.000263922, 0.000395883, 0.000395883, 0.000131961,
716 0.000131961, 0.000263922, 0.000263922, 0, 0,
717 0, 0.000263922, 0.000263922, 0.000263922, 0,
718 0.000131961, 0.000263922, 0.000131961, 0, 0,
719 0.000131961, 0.000131961, 0, 0, 0.000131961,
720 0, 0, 0.000131961, 0.000131961, 0
726 0.000120685, 0, 0, 0, 0,
727 0, 0.000120685, 0, 0.000120685, 0.000603427,
728 0.000844798, 0.00120685, 0.00156891, 0.00229302, 0.00470673,
729 0.00398262, 0.00555153, 0.00796524, 0.00820661, 0.00893073,
730 0.0144823, 0.0112238, 0.0131547, 0.0133961, 0.016896,
731 0.0179821, 0.0189476, 0.0161719, 0.0201545, 0.0225682,
732 0.0225682, 0.0212406, 0.023413, 0.0218441, 0.021482,
733 0.0228096, 0.0197924, 0.0199131, 0.021482, 0.0219648,
734 0.0207579, 0.0208786, 0.021482, 0.0216027, 0.0176201,
735 0.0196717, 0.0164132, 0.0206372, 0.0187063, 0.0162925,
736 0.0138788, 0.0143616, 0.0166546, 0.0159305, 0.0147236,
737 0.013034, 0.0119479, 0.0124306, 0.0111031, 0.00941347,
738 0.010379, 0.00977552, 0.00989621, 0.00941347, 0.0091721,
739 0.00724113, 0.00856867, 0.00905141, 0.00699976, 0.0083273,
740 0.00856867, 0.00687907, 0.00639633, 0.00627565, 0.00518948,
741 0.00446536, 0.00591359, 0.00362056, 0.00362056, 0.00470673,
742 0.00518948, 0.00398262, 0.00289645, 0.00422399, 0.00362056,
743 0.00313782, 0.00265508, 0.00325851, 0.00265508, 0.00277577,
744 0.00289645, 0.00181028, 0.00205165, 0.0016896, 0.00193097,
745 0.00193097, 0.00193097, 0.00144823, 0.0016896, 0.00156891,
746 0.00217234, 0.000965484, 0.00108617, 0.000844798, 0.00144823,
747 0.00144823, 0.000724113, 0.000965484, 0.00181028, 0.000724113,
748 0.000844798, 0.000965484, 0.000603427, 0.00120685, 0.000965484,
749 0.00144823, 0.000724113, 0.00108617, 0.000724113, 0.000362056,
750 0.000241371, 0.000241371, 0.000844798, 0.000724113, 0.000844798,
751 0.000362056, 0.000482742, 0.000362056, 0.000724113, 0.000844798,
752 0.000362056, 0.000120685, 0.000603427, 0.000603427, 0.000241371,
753 0.000241371, 0.000603427, 0.000724113, 0.000482742, 0.000120685,
754 0.000120685, 0.000482742, 0.000241371, 0.000241371, 0.000362056,
755 0.000844798, 0, 0.000120685, 0.000724113, 0.000120685,
756 0.000362056, 0.000603427, 0.000241371, 0.000482742, 0.000120685,
757 0, 0, 0.000241371, 0.000120685, 0.000120685,
758 0.000120685, 0.000362056, 0, 0.000241371, 0.000241371,
759 0.000120685, 0.000241371, 0.000241371, 0.000241371, 0,
760 0.000120685, 0.000120685, 0, 0, 0.000120685,
761 0.000241371, 0.000120685, 0.000120685, 0.000241371, 0,
762 0.000241371, 0, 0.000120685, 0, 0.000120685
767 0, 0, 0.000113921, 0, 0,
768 0, 0, 0.000113921, 0, 0,
769 0, 0, 0, 0, 0.000455685,
770 0.000569606, 0.00102529, 0.0015949, 0.00250627, 0.00364548,
771 0.00364548, 0.00535429, 0.00501253, 0.0075188, 0.00774664,
772 0.0119617, 0.0118478, 0.0125313, 0.0126452, 0.0149237,
773 0.0140123, 0.0148098, 0.0170882, 0.0197084, 0.018797,
774 0.0194805, 0.0179995, 0.0209615, 0.018797, 0.0197084,
775 0.0200501, 0.0191388, 0.020164, 0.020278, 0.0206197,
776 0.0192527, 0.0191388, 0.0192527, 0.0225564, 0.0172021,
777 0.018797, 0.0176578, 0.0150376, 0.0189109, 0.0168603,
778 0.0150376, 0.0137845, 0.0181135, 0.0134427, 0.0192527,
779 0.0140123, 0.0132149, 0.0109364, 0.01162, 0.0119617,
780 0.010139, 0.0110504, 0.00865801, 0.00899977, 0.011506,
781 0.00877193, 0.00843017, 0.00843017, 0.00899977, 0.00922761,
782 0.00865801, 0.00740488, 0.00774664, 0.0104807, 0.00569606,
783 0.00535429, 0.00421508, 0.00535429, 0.00603782, 0.00535429,
784 0.00546822, 0.004329, 0.00364548, 0.00444293, 0.00364548,
785 0.00512645, 0.00239234, 0.00478469, 0.00353156, 0.00307587,
786 0.00341763, 0.00284803, 0.00330371, 0.00273411, 0.00262019,
787 0.00307587, 0.00227842, 0.00205058, 0.0021645, 0.00170882,
788 0.00170882, 0.00296195, 0.00125313, 0.00113921, 0.00193666,
789 0.0021645, 0.00113921, 0.00148098, 0.00136705, 0.00148098,
790 0.00170882, 0.00125313, 0.000911369, 0.00102529, 0.00148098,
791 0.000797448, 0.00125313, 0.00182274, 0.00125313, 0.000683527,
792 0.000911369, 0.000569606, 0.000683527, 0.000683527, 0.00125313,
793 0.00125313, 0.000341763, 0.000911369, 0.000113921, 0.000797448,
794 0.000341763, 0.000683527, 0.000455685, 0.000569606, 0.00125313,
795 0.000341763, 0.000683527, 0.000341763, 0.000113921, 0.000455685,
796 0.000569606, 0.000341763, 0.000341763, 0.000455685, 0.000455685,
797 0.000455685, 0.000113921, 0.000569606, 0.000683527, 0.000113921,
798 0.000455685, 0.000455685, 0.000227842, 0.000227842, 0.000341763,
799 0.000227842, 0.000341763, 0.000341763, 0, 0.000227842,
800 0.000227842, 0.000113921, 0.000341763, 0.000227842, 0.000113921,
801 0.000683527, 0.000113921, 0.000227842, 0.000341763, 0.000455685,
802 0.000227842, 0.000113921, 0.000455685, 0.000227842, 0.000455685,
803 0.000227842, 0.000227842, 0.000113921, 0.000227842, 0.000113921,
804 0.000113921, 0.000113921, 0, 0.000113921, 0.000113921
810 0, 0.000113999, 0, 0, 0,
811 0, 0.000113999, 0.000113999, 0.000569995, 0.000227998,
812 0.000341997, 0.000797994, 0.00102599, 0.00239398, 0.00273598,
813 0.00410397, 0.00421797, 0.00615595, 0.00786594, 0.00854993,
814 0.00900593, 0.0106019, 0.0101459, 0.0113999, 0.0134519,
815 0.0140219, 0.0158459, 0.0152759, 0.0163019, 0.0167579,
816 0.0188098, 0.0159599, 0.0182399, 0.0182399, 0.0199498,
817 0.0230278, 0.0196078, 0.0185819, 0.0194938, 0.0183539,
818 0.0173279, 0.0166439, 0.0191518, 0.0193798, 0.0188098,
819 0.0194938, 0.0186959, 0.0168719, 0.0153899, 0.0151619,
820 0.0139079, 0.0173279, 0.0157319, 0.0167579, 0.0155039,
821 0.0128819, 0.0142499, 0.0129959, 0.0125399, 0.0112859,
822 0.0139079, 0.0112859, 0.0100319, 0.0117419, 0.0116279,
823 0.0102599, 0.0112859, 0.00809394, 0.00797994, 0.00820793,
824 0.00946192, 0.00968992, 0.00889193, 0.00889193, 0.00706794,
825 0.00729594, 0.00615595, 0.00683995, 0.00455996, 0.00444596,
826 0.00524396, 0.00433197, 0.00478796, 0.00376197, 0.00421797,
827 0.00501596, 0.00604195, 0.00376197, 0.00341997, 0.00330597,
828 0.00353397, 0.00353397, 0.00319197, 0.00296398, 0.00205198,
829 0.00227998, 0.00193798, 0.00284998, 0.00205198, 0.00227998,
830 0.00250798, 0.00159599, 0.00239398, 0.00159599, 0.00170999,
831 0.00182399, 0.00205198, 0.00125399, 0.00205198, 0.00125399,
832 0.000911993, 0.00159599, 0.00125399, 0.00182399, 0.00113999,
833 0.00148199, 0.00102599, 0.00227998, 0.00102599, 0.000911993,
834 0.000797994, 0.00113999, 0.000911993, 0.000911993, 0.000683995,
835 0.00102599, 0.00125399, 0.00102599, 0.000797994, 0.000569995,
836 0.000683995, 0.000455996, 0.000683995, 0.000569995, 0.000341997,
837 0.000455996, 0.000683995, 0.000455996, 0.000683995, 0.000341997,
838 0.000683995, 0.000569995, 0.000569995, 0.000455996, 0.000911993,
839 0.000341997, 0.000113999, 0.000341997, 0.000455996, 0.000227998,
840 0.000227998, 0.000227998, 0.000683995, 0.000227998, 0.000227998,
841 0.000341997, 0.000227998, 0.000455996, 0.000227998, 0.000227998,
842 0.000227998, 0.000227998, 0.000227998, 0.000113999, 0.000341997,
843 0.000455996, 0.000113999, 0.000455996, 0, 0.000341997,
844 0.000113999, 0.000341997, 0.000227998, 0.000227998, 0,
845 0.000113999, 0.000113999, 0.000455996, 0.000341997, 0.000227998,
846 0, 0, 0.000113999, 0.000113999, 0
852 0, 0, 0.000107181, 0, 0,
853 0, 0.000107181, 0, 0, 0.000107181,
854 0.000321543, 0.000643087, 0.000321543, 0.0022508, 0.00203644,
855 0.00375134, 0.00460879, 0.00407288, 0.00621651, 0.00782422,
856 0.00953912, 0.0101822, 0.00986066, 0.0113612, 0.0130761,
857 0.0114684, 0.017149, 0.0148982, 0.0157556, 0.018328,
858 0.0145766, 0.01597, 0.0178992, 0.0150054, 0.0187567,
859 0.0180064, 0.0180064, 0.0186495, 0.0184352, 0.0175777,
860 0.019507, 0.0185423, 0.0191854, 0.0219721, 0.0160772,
861 0.0176849, 0.0165059, 0.0184352, 0.0178992, 0.0150054,
862 0.0148982, 0.0154341, 0.0172562, 0.0148982, 0.0145766,
863 0.0153269, 0.0142551, 0.0128617, 0.0132905, 0.0129689,
864 0.0129689, 0.0118971, 0.011254, 0.00889603, 0.0107181,
865 0.00857449, 0.0117899, 0.00932476, 0.00750268, 0.0091104,
866 0.00878885, 0.0073955, 0.00868167, 0.00814577, 0.00718114,
867 0.00664523, 0.00664523, 0.00685959, 0.00685959, 0.00750268,
868 0.0073955, 0.00610932, 0.00600214, 0.00643087, 0.00546624,
869 0.00353698, 0.00493033, 0.00450161, 0.00450161, 0.00439443,
870 0.00321543, 0.0034298, 0.00364416, 0.00321543, 0.00257235,
871 0.00321543, 0.00321543, 0.00278671, 0.00300107, 0.00257235,
872 0.00278671, 0.00246517, 0.00246517, 0.00128617, 0.00257235,
873 0.00192926, 0.00203644, 0.00160772, 0.00160772, 0.00182208,
874 0.00160772, 0.00203644, 0.00214362, 0.000643087, 0.00150054,
875 0.000750268, 0.00139335, 0.000857449, 0.000857449, 0.00096463,
876 0.00128617, 0.000535906, 0.000857449, 0.00117899, 0.000535906,
877 0.000750268, 0.000750268, 0.000643087, 0.000857449, 0.000643087,
878 0.000750268, 0.00096463, 0.000428725, 0.000535906, 0.00128617,
879 0.000428725, 0.000643087, 0.000857449, 0.000750268, 0.000428725,
880 0.000857449, 0.000750268, 0.000321543, 0.000643087, 0.000643087,
881 0.000214362, 0.000107181, 0.000750268, 0.00107181, 0.000428725,
882 0.000643087, 0.000428725, 0.000321543, 0.000643087, 0.000428725,
883 0.000321543, 0.000321543, 0.000643087, 0.000321543, 0.000107181,
884 0.000535906, 0.000428725, 0.000321543, 0.000535906, 0.000214362,
885 0, 0.000428725, 0.000214362, 0.000107181, 0.000214362,
886 0, 0.000107181, 0.000535906, 0, 0.000321543,
887 0.000107181, 0.000321543, 0.000428725, 0.000107181, 0.000321543,
888 0.000214362, 0.000214362, 0.000107181, 0.000428725, 0.000107181
895 0.000109004, 0, 0, 0, 0.000327011,
896 0.000436015, 0.00141705, 0.000981033, 0.00174406, 0.00348812,
897 0.00316111, 0.00425114, 0.00632221, 0.00664923, 0.00708524,
898 0.00632221, 0.0101373, 0.00959233, 0.0126444, 0.00981033,
899 0.0125354, 0.0141705, 0.0143885, 0.0148245, 0.0174406,
900 0.0160235, 0.0159145, 0.0158055, 0.0195117, 0.0163506,
901 0.0182036, 0.0167866, 0.0208197, 0.0177676, 0.0195117,
902 0.0180946, 0.0167866, 0.0186396, 0.0198387, 0.0164596,
903 0.0182036, 0.0165686, 0.0163506, 0.0163506, 0.0140615,
904 0.0161325, 0.0138435, 0.0148245, 0.0151515, 0.0144975,
905 0.0131894, 0.0149335, 0.0123174, 0.0139525, 0.0122084,
906 0.0116634, 0.0126444, 0.0105734, 0.0101373, 0.00937432,
907 0.0104644, 0.0100283, 0.0126444, 0.0089383, 0.00915631,
908 0.0088293, 0.00806627, 0.0089383, 0.00904731, 0.00741225,
909 0.00719424, 0.00632221, 0.00719424, 0.0057772, 0.0057772,
910 0.00457816, 0.00654022, 0.00555919, 0.00468716, 0.00468716,
911 0.00436015, 0.00425114, 0.00446915, 0.00403314, 0.00403314,
912 0.00479616, 0.00392413, 0.00370613, 0.00512317, 0.0028341,
913 0.00261609, 0.0029431, 0.00370613, 0.00218007, 0.0029431,
914 0.00272509, 0.00207107, 0.00228908, 0.00228908, 0.00174406,
915 0.00250709, 0.00337911, 0.000981033, 0.00185306, 0.00185306,
916 0.00196207, 0.00152605, 0.00130804, 0.00218007, 0.00185306,
917 0.00218007, 0.00130804, 0.00130804, 0.000981033, 0.00119904,
918 0.000981033, 0.00141705, 0.00119904, 0.00152605, 0.00109004,
919 0.00141705, 0.000763026, 0.000654022, 0.00119904, 0.00109004,
920 0.00141705, 0.000981033, 0.000327011, 0.000981033, 0.000436015,
921 0.00119904, 0.00130804, 0.00130804, 0.000981033, 0.00087203,
922 0.00087203, 0.000763026, 0.000763026, 0.000763026, 0.000545019,
923 0.000545019, 0.00109004, 0.000654022, 0.000981033, 0.000327011,
924 0.000436015, 0.000218007, 0.000545019, 0.000654022, 0.000654022,
925 0.000436015, 0.000436015, 0.000218007, 0.00130804, 0.00087203,
926 0.000327011, 0.000218007, 0.000654022, 0.000109004, 0.000436015,
927 0.000109004, 0.000436015, 0.000218007, 0.000436015, 0.000109004,
928 0.000327011, 0.000545019, 0.000436015, 0.000654022, 0.000327011,
929 0.000327011, 0, 0.000218007, 0.000436015, 0.000109004,
930 0.000109004, 0.000327011, 0.000218007, 0.000218007, 0.000109004
934 // Time Bin of Max. Cluster Distributions for e and pi
938 Int_t timBin[kNo_Mom][kNo_TBins]={
940 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
941 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
944 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
945 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
948 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
949 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
952 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
953 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
956 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
957 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
960 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
961 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
964 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
965 10, 11, 12, 13, 14, 15, 16, 17, 18, 19
969 Double_t probPIT[kNo_Mom][kNo_TBins]={
971 0, 0.177075, 0.0771579, 0.0666112, 0.0653622 ,
972 0.0595337, 0.0577297, 0.056897, 0.055093, 0.0514849 ,
973 0.0552318, 0.0463503, 0.0496808, 0.0523175, 0.0478768 ,
974 0.0512073, 0.0274771, 0.00291424, 0, 0
977 0, 0.179185, 0.0719195, 0.0729013, 0.0614875 ,
978 0.0609966, 0.0554737, 0.0551055, 0.0565783, 0.0517919 ,
979 0.0543692, 0.0544919, 0.0478645, 0.0542464, 0.0471281 ,
980 0.0459008, 0.0293324, 0.0012273, 0, 0
983 0, 0.173694, 0.0767159, 0.0631295, 0.0646521 ,
984 0.0596158, 0.0576247, 0.058796, 0.0532912, 0.0598501 ,
985 0.0489576, 0.0561021, 0.0482549, 0.0539939, 0.0460295 ,
986 0.0510658, 0.0270555, 0.00117123, 0, 0
989 0, 0.175456, 0.0840475, 0.0667985, 0.0626236 ,
990 0.0564711, 0.0563612, 0.0577895, 0.0541639, 0.0537245 ,
991 0.04922, 0.0535047, 0.0520765, 0.053285, 0.0476818 ,
992 0.0487805, 0.0269172, 0.00109866, 0, 0
995 0, 0.173835, 0.0757168, 0.0704525, 0.0609319 ,
996 0.0601478, 0.0536514, 0.0562276, 0.0521953, 0.0545475 ,
997 0.0524194, 0.0528674, 0.0506272, 0.0543235, 0.0500672 ,
998 0.0502912, 0.030802, 0.000896057, 0, 0
1001 0, 0.175938, 0.0821475, 0.0642518, 0.0600474 ,
1002 0.0620957, 0.0534713, 0.058646, 0.0560586, 0.0529323 ,
1003 0.0523933, 0.0557352, 0.0500216, 0.0502372, 0.0525011 ,
1004 0.0467874, 0.0255498, 0.00118586, 0, 0
1007 0, 0.183151, 0.0766431, 0.0645472, 0.064012 ,
1008 0.0618711, 0.0558767, 0.0557696, 0.0543781, 0.0556626 ,
1009 0.0497752, 0.0520231, 0.049347, 0.0507386, 0.0470991 ,
1010 0.0525583, 0.0260116, 0.000535217, 0, 0
1014 Double_t probELT[kNo_Mom][kNo_TBins]={
1016 0, 0.121598, 0.0534094, 0.0572139, 0.0507755,
1017 0.048288, 0.0453614, 0.0471174, 0.040679, 0.0424349,
1018 0.0463857, 0.0569213, 0.0595552, 0.063506, 0.0690664,
1019 0.0942347, 0.0896986, 0.0137548, 0, 0
1022 0, 0.116522, 0.0500132, 0.0538401, 0.0451306,
1023 0.0461863, 0.0407759, 0.0459224, 0.0482977, 0.0501452,
1024 0.0496173, 0.0560834, 0.053972, 0.0680918, 0.0695434,
1025 0.0991027, 0.0954078, 0.0113486, 0, 0
1028 0, 0.108376, 0.0517741, 0.0483949, 0.0456191,
1029 0.0482742, 0.0428434, 0.0475501, 0.0421192, 0.0469467,
1030 0.0521361, 0.0538257, 0.0562394, 0.0650495, 0.0742216,
1031 0.106445, 0.0954622, 0.0147236, 0, 0
1034 0, 0.104124, 0.0471634, 0.0464798, 0.0416951,
1035 0.0399863, 0.0443153, 0.0467077, 0.0461381, 0.0500114,
1036 0.0498975, 0.0602643, 0.0603782, 0.0682388, 0.0811119,
1037 0.106402, 0.0920483, 0.0150376, 0, 0
1040 0, 0.101573, 0.0478796, 0.0441176, 0.0411537,
1041 0.0454856, 0.0449156, 0.0463976, 0.0467396, 0.0523256,
1042 0.0512996, 0.0605335, 0.0580255, 0.0649795, 0.0802554,
1043 0.0999772, 0.0992932, 0.0150479, 0, 0
1046 0, 0.103966, 0.0463023, 0.046731, 0.0398714,
1047 0.044373, 0.0437299, 0.0435155, 0.0459807, 0.0483387,
1048 0.049732, 0.0531618, 0.0613076, 0.0767417, 0.0811361,
1049 0.101501, 0.0982851, 0.0153269, 0, 0
1052 0, 0.102245, 0.0491607, 0.0417484, 0.0443645,
1053 0.0408764, 0.0406584, 0.0448005, 0.0501417, 0.0514497,
1054 0.0537388, 0.058208, 0.0604971, 0.0671463, 0.0783737,
1055 0.101809, 0.0987574, 0.0160235, 0, 0
1059 for(Int_t ip=0; ip<fNMom; ip++) {
1060 fTrackMomentum[ip]= trackMomentum[ip];
1061 for(Int_t ie=0; ie<fNEbins; ie++) {
1062 fEnergyLoss[ip][ie]=energyLoss[ip][ie];
1063 fProbPI[ip][ie]=probPI[ip][ie];
1064 fProbEL[ip][ie]=probEL[ip][ie];
1068 for(Int_t ip=0; ip<fNMom; ip++) {
1069 fTrackMomentum[ip]=trackMomentum[ip] ;
1070 for(Int_t it=0; it<fNTbins; it++) {
1071 fTimBin[ip][it]=timBin[ip][it];
1072 fProbPIT[ip][it]=probPIT[ip][it];
1073 fProbELT[ip][it]=probELT[ip][it];
1081 void AliTRDprobdist::FillData()
1084 // Energy loss Distributions for e and pi
1086 Double_t fADCNorm=1.;
1087 ifstream filein("pp.dat");
1089 for(Int_t ip=0; ip<fNMom; ip++) {
1090 filein >> fTrackMomentum[ip] >> fNEbins;
1091 for(Int_t ie=0; ie<fNEbins; ie++) {
1092 filein >> fEnergyLoss[ip][ie]>> fProbPI[ip][ie]>> fProbEL[ip][ie];
1093 fEnergyLoss[ip][ie]=fADCNorm*fEnergyLoss[ip][ie];
1096 fEnBinSize=fEnergyLoss[0][1]-fEnergyLoss[0][0];
1098 // Time Bin of Max. Cluster Distributions for e and pi
1099 ifstream fileint("tt.dat");
1101 for(Int_t ip=0; ip<fNMom; ip++) {
1102 fileint >> fTrackMomentum[ip] >> fNTbins;
1103 for(Int_t it=0; it<fNTbins; it++) {
1104 fileint >> fTimBin[ip][it]>> fProbPIT[ip][it]>> fProbELT[ip][it];