New response functions for tracking updates
[u/mrichter/AliRoot.git] / TRD / AliTRDprobdist.cxx
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
26 ClassImp(AliTRDprobdist)
27
28 //_________________________________________________________________________
29 AliTRDprobdist::AliTRDprobdist(Int_t multiplicity)
30 {
31   //
32   //  The main constructor
33   //
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();
39 }
40
41 //_________________________________________________________________________
42 Double_t  AliTRDprobdist::GetMeanPI(Int_t ip) const
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 //_________________________________________________________________________
55 Double_t  AliTRDprobdist::GetMeanEL(Int_t ip) const
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 //_________________________________________________________________________
69 Double_t  AliTRDprobdist::GetNormalizationPI(Int_t ip) const
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 //_________________________________________________________________________
81 Double_t  AliTRDprobdist::GetNormalizationEL(Int_t ip) const
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 //_________________________________________________________________________
93 Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx) const
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 
99   Double_t probability = 1.0;
100   Int_t iEnBin= ((Int_t) (dedx/fEnBinSize));
101   if(iEnBin > fNEbins-1) iEnBin = fNEbins-1;
102
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];
106   }
107   if(k==2){    // pion
108     if(mom<=fTrackMomentum[0]) probability = fProbPI[0][iEnBin];
109     if(mom>=fTrackMomentum[fNMom-1]) probability = fProbPI[fNMom-1][iEnBin];
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
117     probability= fProbEL[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
118     return probability;
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
126     probability= fProbPI[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
127     return probability;
128   }
129   return probability;
130 }
131
132
133 //_________________________________________________________________________
134 Double_t AliTRDprobdist::GetProbabilityT(Int_t k, Double_t mom, Int_t timbin) const
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 
140   Double_t probabilityT = 1.0;
141   if(timbin<=0) return 0.;
142   Int_t iTBin=timbin;
143
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];
147   }
148   if(k==2){    // pion
149     if(mom<=fTrackMomentum[0]) probabilityT = fProbPIT[0][iTBin];
150     if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbPIT[fNMom-1][iTBin];
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
158     probabilityT= fProbELT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
159     return probabilityT;
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
167     probabilityT= fProbPIT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
168     return probabilityT;
169   }
170   return probabilityT;
171 }
172
173
174 void AliTRDprobdist::FillData() 
175 {
176   //
177   // Energy loss Distributions for e and pi
178   fNMom=7;
179   Double_t trackMomentum[kNo_Mom]= {1, 1.5, 2, 3, 4, 5, 6};
180   fNEbins=250;
181   fEnBinSize=10;
182
183   Double_t energyLoss[kNo_Mom][kNo_EnBins]={  
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, 
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
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, 
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
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, 
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
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, 
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
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, 
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
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, 
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
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, 
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
372             }
373           };
374
375   Double_t probPI[kNo_Mom][kNo_EnBins]={  
376             {
377                0, 0, 0, 0, 0.000133547, 
378                0, 0, 0, 0, 0, 
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, 
416                0, 0, 0, 0, 0, 
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
427             }, 
428             {
429                0, 0, 0, 0, 0, 
430                0, 0, 0, 0, 0, 
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
479             }, 
480             {
481                0, 0, 0, 0, 0, 
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
531             }, 
532             {
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, 
581                0, 0, 0, 0, 0, 
582                0, 0, 0.000117536, 0, 0
583             }, 
584             {
585                0, 0, 0, 0, 0, 
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
635             }, 
636             {
637                0, 0, 0, 0, 0, 
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
687             }, 
688             {
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, 
733                0, 0, 0, 0, 0, 
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
739             }
740           };
741
742   Double_t probEL[kNo_Mom][kNo_EnBins]={  
743             {
744                0, 0, 0, 0, 0, 
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, 
790                0, 0, 0, 0, 0, 
791                0, 0.000141603, 0, 0, 0, 
792                0, 0.000141603, 0, 0, 0, 
793                0.000141603, 0, 0, 0.000141603, 0
794             }, 
795             {
796                0, 0, 0, 0, 0, 
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, 
837                0, 0, 0, 0, 0, 
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, 
842                0, 0, 0, 0, 0, 
843                0.000276855, 0.000138427, 0, 0, 0, 
844                0, 0, 0, 0, 0, 
845                0, 0, 0.000276855, 0.000276855, 0
846             }, 
847             {
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
898             }, 
899             {
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, 
946                0, 0, 0, 0, 0, 
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
950             }, 
951             {
952                0, 0, 0, 0, 0, 
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
1002             }, 
1003             {
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
1054             }, 
1055             {
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
1106             }
1107           };
1108
1109   // Time Bin of Max. Cluster Distributions for e and pi
1110   fNTbins=20;
1111   fTBinSize=1;
1112
1113   Int_t timBin[kNo_Mom][kNo_TBins]={  
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
1144   Double_t probPIT[kNo_Mom][kNo_TBins]={  
1145             {
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
1150             }, 
1151             {
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
1156             }, 
1157             {
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
1162             }, 
1163             {
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
1168             }, 
1169             {
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
1174             }, 
1175             {
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
1180             }, 
1181             {
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
1186             }
1187           };
1188
1189   Double_t probELT[kNo_Mom][kNo_TBins]={  
1190             {
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
1195             }, 
1196             {
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
1201             }, 
1202             {
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
1207             }, 
1208             {
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
1213             }, 
1214             {
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
1219             }, 
1220             {
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
1225             }, 
1226             {
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
1231             }
1232           };
1233
1234   for(Int_t ip=0; ip<fNMom; ip++) {
1235     fTrackMomentum[ip]= trackMomentum[ip];
1236     for(Int_t ie=0; ie<fNEbins; ie++) {
1237       fEnergyLoss[ip][ie]=energyLoss[ip][ie];
1238       fProbPI[ip][ie]=probPI[ip][ie]; 
1239       fProbEL[ip][ie]=probEL[ip][ie];
1240     }
1241   }
1242
1243   for(Int_t ip=0; ip<fNMom; ip++) {
1244     fTrackMomentum[ip]=trackMomentum[ip] ;
1245     for(Int_t it=0; it<fNTbins; it++) {
1246       fTimBin[ip][it]=timBin[ip][it];
1247       fProbPIT[ip][it]=probPIT[ip][it];
1248       fProbELT[ip][it]=probELT[ip][it];
1249     }
1250   }
1251
1252 }
1253
1254
1255 /*
1256 void 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