Updated clusteriser parameters
[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=200;
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             }, 
206             {
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
227             }, 
228             {
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
249             }, 
250             {
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
271             }, 
272             {
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
293             }, 
294             {
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
315             }, 
316             {
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
337             }
338           };
339
340   Double_t probPI[kNo_Mom][kNo_EnBins]={  
341             {
342                0, 0, 0, 0, 0, 
343                0, 0, 0, 0, 0, 
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
382             }, 
383             {
384                0, 0, 0, 0, 0, 
385                0, 0, 0, 0, 0, 
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
424             }, 
425             {
426                0, 0, 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
466             }, 
467             {
468                0, 0, 0, 0, 0, 
469                0, 0, 0, 0, 0, 
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
508             }, 
509             {
510                0, 0, 0, 0, 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
550             }, 
551             {
552                0, 0, 0, 0, 0, 
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
592             }, 
593             {
594                0, 0, 0, 0, 0, 
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
634             }
635           };
636
637   Double_t probEL[kNo_Mom][kNo_EnBins]={  
638             {
639                0, 0, 0, 0, 0, 
640                0, 0, 0, 0, 0, 
641                0, 0, 0, 0, 0, 
642                0, 0, 0, 0, 0, 
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, 
677                0, 0, 0, 0, 0, 
678                0.000146327, 0, 0, 0.000146327, 0.000292654
679             }, 
680             {
681                0, 0, 0, 0, 0, 
682                0, 0, 0, 0, 0, 
683                0, 0, 0, 0, 0, 
684                0, 0, 0, 0, 0, 
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
721             }, 
722             {
723                0, 0, 0, 0, 0, 
724                0, 0, 0, 0, 0, 
725                0, 0, 0, 0, 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
763             }, 
764             {
765                0, 0, 0, 0, 0, 
766                0, 0, 0, 0, 0, 
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
805             }, 
806             {
807                0, 0, 0, 0, 0, 
808                0, 0, 0, 0, 0, 
809                0, 0, 0, 0, 0, 
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
847             }, 
848             {
849                0, 0, 0, 0, 0, 
850                0, 0, 0, 0, 0, 
851                0, 0, 0, 0, 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
889             }, 
890             {
891                0, 0, 0, 0, 0, 
892                0, 0, 0, 0, 0, 
893                0, 0, 0, 0, 0, 
894                0, 0, 0, 0, 0, 
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
931             }
932           };
933
934   // Time Bin of Max. Cluster Distributions for e and pi
935   fNTbins=20;
936   fTBinSize=1;
937
938   Int_t timBin[kNo_Mom][kNo_TBins]={  
939             {
940                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
941                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
942             }, 
943             {
944                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
945                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
946             }, 
947             {
948                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
949                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
950             }, 
951             {
952                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
953                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
954             }, 
955             {
956                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
957                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
958             }, 
959             {
960                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
961                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
962             }, 
963             {
964                0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
965                10, 11, 12, 13, 14, 15, 16, 17, 18, 19
966             }
967           };
968
969   Double_t probPIT[kNo_Mom][kNo_TBins]={  
970             {
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
975             }, 
976             {
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
981             }, 
982             {
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
987             }, 
988             {
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
993             }, 
994             {
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
999             }, 
1000             {
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
1005             }, 
1006             {
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
1011             }
1012           };
1013
1014   Double_t probELT[kNo_Mom][kNo_TBins]={  
1015             {
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
1020             }, 
1021             {
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
1026             }, 
1027             {
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
1032             }, 
1033             {
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
1038             }, 
1039             {
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
1044             }, 
1045             {
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
1050             }, 
1051             {
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
1056             }
1057           };
1058
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];
1065     }
1066   }
1067
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];
1074     }
1075   }
1076
1077 }
1078
1079
1080 /*
1081 void AliTRDprobdist::FillData()
1082 {
1083   //
1084   // Energy loss Distributions for e and pi
1085   Double_t temp;
1086   Double_t fADCNorm=1.;
1087   ifstream filein("pp.dat");
1088   filein >> fNMom;
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];
1094     }
1095   }
1096   fEnBinSize=fEnergyLoss[0][1]-fEnergyLoss[0][0];
1097
1098   // Time Bin of Max. Cluster Distributions for e and pi
1099   ifstream fileint("tt.dat");
1100   fileint >> fNMom;
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];
1105     }
1106   }
1107   fTBinSize=1;
1108 }
1109
1110 */
1111