]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/AliTRDprobdist.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
New code from Prashant
[u/mrichter/AliRoot.git] / TRD / AliTRDprobdist.cxx
CommitLineData
33a848f6 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16//-----------------------------------------------------------------
17// Class for dE/dx and Time Bin of Max. Cluster for Electrons and
18// pions in TRD.
19// It is instantiated in class AliTRDpidESD for particle identification
20// in TRD
21// Prashant Shukla <shukla@pi0.physi.uni-heidelberg.de>
22//-----------------------------------------------------------------
23
24#include "AliTRDprobdist.h"
25
26ClassImp(AliTRDprobdist)
27
28//_________________________________________________________________________
1642e650 29AliTRDprobdist::AliTRDprobdist(Int_t multiplicity)
33a848f6 30{
31 //
32 // The main constructor
33 //
1642e650 34 if (multiplicity == 1) FillData();
35 // if (multiplicity == 2000) FillData2000();
36 // if (multiplicity == 4000) FillData4000();
37 // if (multiplicity == 6000) FillData6000();
38 // if (multiplicity == 8000) FillData8000();
33a848f6 39}
40
41//_________________________________________________________________________
1642e650 42Double_t AliTRDprobdist::GetMeanPI(Int_t ip) const
33a848f6 43{
44 // Gets mean of de/dx dist. of pi
45 Double_t integral=0.;
46 Double_t norm=0.;
47 for(Int_t ie=0; ie<fNEbins; ie++) {
48 integral+=fEnergyLoss[ip][ie]*fProbPI[ip][ie];
49 norm+=fProbPI[ip][ie];
50 }
51 return integral/norm;
52}
53
54//_________________________________________________________________________
1642e650 55Double_t AliTRDprobdist::GetMeanEL(Int_t ip) const
33a848f6 56{
57 //
58 // Gets mean of de/dx dist. of e
59 Double_t integral=0.;
60 Double_t norm=0.;
61 for(Int_t ie=0; ie<fNEbins; ie++) {
62 integral+=fEnergyLoss[ip][ie]*fProbEL[ip][ie];
63 norm+=fProbEL[ip][ie];
64 }
65 return integral/norm;
66}
67
68//_________________________________________________________________________
1642e650 69Double_t AliTRDprobdist::GetNormalizationPI(Int_t ip) const
33a848f6 70{
71 //
72 // Gets Normalization of de/dx dist. of pi
73 Double_t integral=0.;
74 for(Int_t ie=0; ie<fNEbins; ie++) {
75 integral+=fProbPI[ip][ie];
76 }
77 return integral;
78}
79
80//_________________________________________________________________________
1642e650 81Double_t AliTRDprobdist::GetNormalizationEL(Int_t ip) const
33a848f6 82{
83 //
84 // Gets Normalization of de/dx dist. of e
85 Double_t integral=0.;
86 for(Int_t ie=0; ie<fNEbins; ie++) {
87 integral+=fProbEL[ip][ie];
88 }
89 return integral;
90}
91
92//_________________________________________________________________________
1642e650 93Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx) const
33a848f6 94{
95 //
96 // Gets the Probability of having dedx at a given momentum (mom)
97 // and particle type k (0 for e) and (2 for pi)
98 // from the precalculated de/dx distributions
1642e650 99 Double_t probability = 1.0;
33a848f6 100 Int_t iEnBin= ((Int_t) (dedx/fEnBinSize));
101 if(iEnBin > fNEbins-1) iEnBin = fNEbins-1;
102
103 if(k==0){ // electron
1642e650 104 if(mom<=fTrackMomentum[0]) probability = fProbEL[0][iEnBin];
105 if(mom>=fTrackMomentum[fNMom-1]) probability = fProbEL[fNMom-1][iEnBin];
33a848f6 106 }
107 if(k==2){ // pion
1642e650 108 if(mom<=fTrackMomentum[0]) probability = fProbPI[0][iEnBin];
109 if(mom>=fTrackMomentum[fNMom-1]) probability = fProbPI[fNMom-1][iEnBin];
33a848f6 110 }
111
112 if(k==0) // electron
113 for(Int_t ip=1; ip<fNMom; ip++)
114 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
115 Double_t slop=(fProbEL[ip][iEnBin]-fProbEL[ip-1][iEnBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
116 // Linear Interpolation
1642e650 117 probability= fProbEL[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
118 return probability;
33a848f6 119 }
120
121 if(k==2) //pion
122 for(Int_t ip=1; ip<fNMom; ip++)
123 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
124 Double_t slop=(fProbPI[ip][iEnBin]-fProbPI[ip-1][iEnBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
125 // Linear Interpolation
1642e650 126 probability= fProbPI[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
127 return probability;
33a848f6 128 }
1642e650 129 return probability;
33a848f6 130}
131
132
133//_________________________________________________________________________
1642e650 134Double_t AliTRDprobdist::GetProbabilityT(Int_t k, Double_t mom, Int_t timbin) const
33a848f6 135{
136 //
137 // Gets the Probability of having timbin at a given momentum (mom)
138 // and particle type k (0 for e) and (2 for pi)
139 // from the precalculated timbin distributions
1642e650 140 Double_t probabilityT = 1.0;
33a848f6 141 if(timbin<=0) return 0.;
142 Int_t iTBin=timbin;
143
144 if(k==0){ // electron
1642e650 145 if(mom<=fTrackMomentum[0]) probabilityT = fProbELT[0][iTBin];
146 if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbELT[fNMom-1][iTBin];
33a848f6 147 }
148 if(k==2){ // pion
1642e650 149 if(mom<=fTrackMomentum[0]) probabilityT = fProbPIT[0][iTBin];
150 if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbPIT[fNMom-1][iTBin];
33a848f6 151 }
152
153 if(k==0) // electron
154 for(Int_t ip=1; ip<fNMom; ip++)
155 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
156 Double_t slop=(fProbELT[ip][iTBin]-fProbELT[ip-1][iTBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
157 // Linear Interpolation
1642e650 158 probabilityT= fProbELT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
159 return probabilityT;
33a848f6 160 }
161
162 if(k==2) // pion
163 for(Int_t ip=1; ip<fNMom; ip++)
164 if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
165 Double_t slop=(fProbPIT[ip][iTBin]-fProbPIT[ip-1][iTBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
166 // Linear Interpolation
1642e650 167 probabilityT= fProbPIT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
168 return probabilityT;
33a848f6 169 }
1642e650 170 return probabilityT;
33a848f6 171}
172
173
33a848f6 174void AliTRDprobdist::FillData()
175{
176 //
177 // Energy loss Distributions for e and pi
178 fNMom=7;
1642e650 179 Double_t trackMomentum[kNo_Mom]= {1., 1.5, 2., 3., 4., 5., 6.};
33a848f6 180 fNEbins=200;
1642e650 181 fEnBinSize=10;
33a848f6 182
1642e650 183 Double_t energyLoss[kNo_Mom][kNo_EnBins]={
33a848f6 184 {
185 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
186 100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
187 200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
188 300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
189 400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
190 500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
191 600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
192 700, 710, 720, 730, 740, 750, 760, 770, 780, 790,
193 800, 810, 820, 830, 840, 850, 860, 870, 880, 890,
194 900, 910, 920, 930, 940, 950, 960, 970, 980, 990,
195 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090,
196 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190,
197 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290,
198 1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390,
199 1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490,
200 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590,
201 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690,
202 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790,
203 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890,
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
1642e650 340 Double_t probPI[kNo_Mom][kNo_EnBins]={
33a848f6 341 {
342 0, 0, 0, 0, 0,
33a848f6 343 0, 0, 0, 0, 0,
1642e650 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
33a848f6 382 },
383 {
384 0, 0, 0, 0, 0,
385 0, 0, 0, 0, 0,
1642e650 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
33a848f6 424 },
425 {
426 0, 0, 0, 0, 0,
1642e650 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
33a848f6 466 },
467 {
468 0, 0, 0, 0, 0,
1642e650 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
33a848f6 508 },
509 {
510 0, 0, 0, 0, 0,
1642e650 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
33a848f6 550 },
551 {
552 0, 0, 0, 0, 0,
1642e650 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
33a848f6 592 },
593 {
594 0, 0, 0, 0, 0,
1642e650 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
33a848f6 634 }
635 };
636
1642e650 637 Double_t probEL[kNo_Mom][kNo_EnBins]={
33a848f6 638 {
639 0, 0, 0, 0, 0,
1642e650 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
33a848f6 679 },
680 {
681 0, 0, 0, 0, 0,
682 0, 0, 0, 0, 0,
683 0, 0, 0, 0, 0,
1642e650 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
33a848f6 721 },
722 {
33a848f6 723 0, 0, 0, 0, 0,
1642e650 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
33a848f6 763 },
764 {
765 0, 0, 0, 0, 0,
1642e650 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
33a848f6 805 },
806 {
807 0, 0, 0, 0, 0,
1642e650 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
33a848f6 847 },
848 {
849 0, 0, 0, 0, 0,
850 0, 0, 0, 0, 0,
1642e650 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
33a848f6 889 },
890 {
891 0, 0, 0, 0, 0,
1642e650 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
33a848f6 931 }
932 };
933
934 // Time Bin of Max. Cluster Distributions for e and pi
935 fNTbins=20;
936 fTBinSize=1;
937
1642e650 938 Int_t timBin[kNo_Mom][kNo_TBins]={
33a848f6 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
1642e650 969 Double_t probPIT[kNo_Mom][kNo_TBins]={
33a848f6 970 {
1642e650 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
33a848f6 975 },
976 {
1642e650 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
33a848f6 981 },
982 {
1642e650 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
33a848f6 987 },
988 {
1642e650 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
33a848f6 993 },
994 {
1642e650 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
33a848f6 999 },
1000 {
1642e650 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
33a848f6 1005 },
1006 {
1642e650 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
33a848f6 1011 }
1012 };
1013
1642e650 1014 Double_t probELT[kNo_Mom][kNo_TBins]={
33a848f6 1015 {
1642e650 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
33a848f6 1020 },
1021 {
1642e650 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
33a848f6 1026 },
1027 {
1642e650 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
33a848f6 1032 },
1033 {
1642e650 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
33a848f6 1038 },
1039 {
1642e650 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
33a848f6 1044 },
1045 {
1642e650 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
33a848f6 1050 },
1051 {
1642e650 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
33a848f6 1056 }
1057 };
1058
1059 for(Int_t ip=0; ip<fNMom; ip++) {
1642e650 1060 fTrackMomentum[ip]= trackMomentum[ip];
33a848f6 1061 for(Int_t ie=0; ie<fNEbins; ie++) {
1642e650 1062 fEnergyLoss[ip][ie]=energyLoss[ip][ie];
1063 fProbPI[ip][ie]=probPI[ip][ie];
1064 fProbEL[ip][ie]=probEL[ip][ie];
33a848f6 1065 }
1066 }
1067
1068 for(Int_t ip=0; ip<fNMom; ip++) {
1642e650 1069 fTrackMomentum[ip]=trackMomentum[ip] ;
33a848f6 1070 for(Int_t it=0; it<fNTbins; it++) {
1642e650 1071 fTimBin[ip][it]=timBin[ip][it];
1072 fProbPIT[ip][it]=probPIT[ip][it];
1073 fProbELT[ip][it]=probELT[ip][it];
33a848f6 1074 }
1075 }
1076
1077}
1078
1079
1080/*
1081void 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