TRD/AliTRDprobdist.cxx

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