TRD/AliTRDprobdist.cxx

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