/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

//-----------------------------------------------------------------
// Class for dE/dx and Time Bin of Max. Cluster for Electrons and 
// pions in TRD. 
// It is instantiated in class AliTRDpidESD for particle identification
// in TRD
// Prashant Shukla <shukla@pi0.physi.uni-heidelberg.de>
//-----------------------------------------------------------------

#include "AliTRDprobdist.h"

ClassImp(AliTRDprobdist)

//_________________________________________________________________________
AliTRDprobdist::AliTRDprobdist(Int_t multiplicity)
{
  //
  //  The main constructor
  //
  if (multiplicity == 1) FillData();
  //  if (multiplicity == 2000) FillData2000();
  //  if (multiplicity == 4000) FillData4000();
  //  if (multiplicity == 6000) FillData6000();
  //  if (multiplicity == 8000) FillData8000();
}

//_________________________________________________________________________
Double_t  AliTRDprobdist::GetMeanPI(Int_t ip) const
{
  // Gets mean of de/dx dist. of pi
  Double_t integral=0.;
  Double_t norm=0.;
  for(Int_t ie=0; ie<fNEbins; ie++) { 
    integral+=fEnergyLoss[ip][ie]*fProbPI[ip][ie]; 
    norm+=fProbPI[ip][ie];
  }
  return integral/norm;
}

//_________________________________________________________________________
Double_t  AliTRDprobdist::GetMeanEL(Int_t ip) const
{
  //
  // Gets mean of de/dx dist. of e
  Double_t integral=0.;
  Double_t norm=0.;
  for(Int_t ie=0; ie<fNEbins; ie++) { 
    integral+=fEnergyLoss[ip][ie]*fProbEL[ip][ie]; 
    norm+=fProbEL[ip][ie];
  }
  return integral/norm;
}

//_________________________________________________________________________
Double_t  AliTRDprobdist::GetNormalizationPI(Int_t ip) const
{
  //
  // Gets Normalization of de/dx dist. of pi
  Double_t integral=0.;
  for(Int_t ie=0; ie<fNEbins; ie++) { 
    integral+=fProbPI[ip][ie]; 
  }
  return integral;
}

//_________________________________________________________________________
Double_t  AliTRDprobdist::GetNormalizationEL(Int_t ip) const
{
  //
  // Gets Normalization of de/dx dist. of e
  Double_t integral=0.;
  for(Int_t ie=0; ie<fNEbins; ie++) { 
    integral+=fProbEL[ip][ie]; 
  }
  return integral;
}

//_________________________________________________________________________
Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx) const
{
  //
  // Gets the Probability of having dedx at a given momentum (mom)
  // and particle type k (0 for e) and (2 for pi)
  // from the precalculated de/dx distributions 
  Double_t probability = 1.0;
  Int_t iEnBin= ((Int_t) (dedx/fEnBinSize));
  if(iEnBin > fNEbins-1) iEnBin = fNEbins-1;

  if(k==0){    // electron
    if(mom<=fTrackMomentum[0]) probability = fProbEL[0][iEnBin];
    if(mom>=fTrackMomentum[fNMom-1]) probability = fProbEL[fNMom-1][iEnBin];
  }
  if(k==2){    // pion
    if(mom<=fTrackMomentum[0]) probability = fProbPI[0][iEnBin];
    if(mom>=fTrackMomentum[fNMom-1]) probability = fProbPI[fNMom-1][iEnBin];
  }

  if(k==0)    // electron
  for(Int_t ip=1; ip<fNMom; ip++)  
  if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
    Double_t slop=(fProbEL[ip][iEnBin]-fProbEL[ip-1][iEnBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
    // Linear Interpolation
    probability= fProbEL[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
    return probability;
  }

  if(k==2)   //pion
  for(Int_t ip=1; ip<fNMom; ip++)
  if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
    Double_t slop=(fProbPI[ip][iEnBin]-fProbPI[ip-1][iEnBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
    // Linear Interpolation
    probability= fProbPI[ip-1][iEnBin] + slop*(mom-fTrackMomentum[ip-1]);
    return probability;
  }
  return probability;
}


//_________________________________________________________________________
Double_t AliTRDprobdist::GetProbabilityT(Int_t k, Double_t mom, Int_t timbin) const
{
  //
  // Gets the Probability of having timbin at a given momentum (mom)
  // and particle type k (0 for e) and (2 for pi)
  // from the precalculated timbin distributions 
  Double_t probabilityT = 1.0;
  if(timbin<=0) return 0.;
  Int_t iTBin=timbin;

  if(k==0){    // electron
    if(mom<=fTrackMomentum[0]) probabilityT = fProbELT[0][iTBin];
    if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbELT[fNMom-1][iTBin];
  }
  if(k==2){    // pion
    if(mom<=fTrackMomentum[0]) probabilityT = fProbPIT[0][iTBin];
    if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fProbPIT[fNMom-1][iTBin];
  }

  if(k==0)    // electron
  for(Int_t ip=1; ip<fNMom; ip++)
  if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
    Double_t slop=(fProbELT[ip][iTBin]-fProbELT[ip-1][iTBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
    // Linear Interpolation
    probabilityT= fProbELT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
    return probabilityT;
  }

  if(k==2)   // pion
  for(Int_t ip=1; ip<fNMom; ip++)
  if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
    Double_t slop=(fProbPIT[ip][iTBin]-fProbPIT[ip-1][iTBin])/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
    // Linear Interpolation
    probabilityT= fProbPIT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
    return probabilityT;
  }
  return probabilityT;
}


void AliTRDprobdist::FillData() 
{
  //
  // Energy loss Distributions for e and pi
  fNMom=7;
  Double_t trackMomentum[kNo_Mom]= {1, 1.5, 2, 3, 4, 5, 6};
  fNEbins=250;
  fEnBinSize=10;

  Double_t energyLoss[kNo_Mom][kNo_EnBins]={  
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }, 
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }, 
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }, 
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }, 
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }, 
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }, 
            {
               0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 
               100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 
               200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 
               300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 
               400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 
               500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 
               600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 
               700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 
               800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 
               900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 
               1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 
               1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 
               1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 
               1300, 1310, 1320, 1330, 1340, 1350, 1360, 1370, 1380, 1390, 
               1400, 1410, 1420, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 
               1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 
               1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 
               1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 
               1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 
               1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 
               2000, 2010, 2020, 2030, 2040, 2050, 2060, 2070, 2080, 2090, 
               2100, 2110, 2120, 2130, 2140, 2150, 2160, 2170, 2180, 2190, 
               2200, 2210, 2220, 2230, 2240, 2250, 2260, 2270, 2280, 2290, 
               2300, 2310, 2320, 2330, 2340, 2350, 2360, 2370, 2380, 2390, 
               2400, 2410, 2420, 2430, 2440, 2450, 2460, 2470, 2480, 2490
            }
          };

  Double_t probPI[kNo_Mom][kNo_EnBins]={  
            {
               0, 0, 0, 0, 0.000133547, 
               0, 0, 0, 0, 0, 
               0.000133547, 0, 0, 0, 0, 
               0.000667735, 0.000267094, 0.000934829, 0.00186966, 0.00280449, 
               0.00413996, 0.00534188, 0.0117521, 0.0132212, 0.0205662, 
               0.0206998, 0.0220353, 0.0251068, 0.0331197, 0.0335203, 
               0.0363248, 0.0371261, 0.0413996, 0.0403312, 0.0376603, 
               0.0419338, 0.0351229, 0.0345887, 0.0301816, 0.0347222, 
               0.0319177, 0.0308494, 0.0295139, 0.0247062, 0.0257746, 
               0.0220353, 0.0200321, 0.0202991, 0.0189637, 0.0153579, 
               0.0140224, 0.0150908, 0.0121528, 0.0125534, 0.0110844, 
               0.0101496, 0.00894765, 0.00734509, 0.0068109, 0.00467415, 
               0.00667735, 0.00534188, 0.00320513, 0.00320513, 0.00333868, 
               0.00400641, 0.00307158, 0.00280449, 0.00240385, 0.00333868, 
               0.00253739, 0.00240385, 0.00240385, 0.00186966, 0.00106838, 
               0.00240385, 0.000667735, 0.00200321, 0.00133547, 0.00173611, 
               0.00133547, 0.000534188, 0.00106838, 0.000534188, 0.000667735, 
               0.000801282, 0.00106838, 0.000667735, 0.00106838, 0.000934829, 
               0.000400641, 0.00133547, 0.00106838, 0.000801282, 0.000801282, 
               0.000667735, 0.000667735, 0.000534188, 0.000133547, 0.000534188, 
               0.000133547, 0.000667735, 0.000534188, 0.000400641, 0.000400641, 
               0.000400641, 0.00106838, 0.000133547, 0.000267094, 0.000534188, 
               0.000400641, 0.000133547, 0.000534188, 0.000267094, 0.000133547, 
               0.000133547, 0, 0.000133547, 0.000400641, 0.000400641, 
               0.000400641, 0.000400641, 0.000133547, 0.000133547, 0.000133547, 
               0.000534188, 0, 0.000133547, 0, 0.000267094, 
               0.000133547, 0.000133547, 0, 0.000267094, 0.000267094, 
               0.000133547, 0.000133547, 0.000801282, 0, 0.000400641, 
               0.000267094, 0.000133547, 0, 0.000267094, 0.000267094, 
               0.000267094, 0.000267094, 0, 0.000133547, 0.000133547, 
               0.000133547, 0, 0, 0, 0, 
               0.000133547, 0, 0, 0, 0.000267094, 
               0, 0.000133547, 0, 0, 0.000133547, 
               0.000400641, 0, 0.000133547, 0.000133547, 0.000133547, 
               0.000133547, 0, 0.000133547, 0.000133547, 0, 
               0.000133547, 0.000133547, 0, 0, 0.000133547, 
               0.000133547, 0.000267094, 0.000133547, 0, 0.000400641, 
               0.000133547, 0.000400641, 0, 0, 0, 
               0, 0.000267094, 0, 0, 0.000133547, 
               0, 0, 0, 0, 0, 
               0.000133547, 0.000400641, 0.000133547, 0, 0.000133547, 
               0.000133547, 0, 0.000133547, 0.000133547, 0, 
               0.000133547, 0, 0, 0.000133547, 0.000133547, 
               0, 0.000133547, 0, 0, 0, 
               0.000133547, 0, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0, 0.000133547, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0, 0.000133547, 0, 0, 0, 
               0, 0, 0, 0, 0
            }, 
            {
               0, 0, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0.000159642, 0, 0.000159642, 0.00063857, 0.000159642, 
               0, 0, 0.000159642, 0.00063857, 0.00175607, 
               0.00303321, 0.00462963, 0.00718391, 0.00846105, 0.00989783, 
               0.0164432, 0.019636, 0.0244253, 0.024106, 0.0316092, 
               0.0359195, 0.0371967, 0.0351213, 0.0387931, 0.0403895, 
               0.0386335, 0.0371967, 0.0320881, 0.0371967, 0.0320881, 
               0.0341635, 0.0333653, 0.0304917, 0.0301724, 0.0252235, 
               0.0228289, 0.0204342, 0.0231481, 0.0189974, 0.0205939, 
               0.0177203, 0.0150064, 0.0143678, 0.0119732, 0.0127714, 
               0.0116539, 0.0102171, 0.0113346, 0.00670498, 0.0063857, 
               0.00718391, 0.00718391, 0.00686462, 0.0052682, 0.0063857, 
               0.00431034, 0.00319285, 0.00399106, 0.00287356, 0.00367178, 
               0.00287356, 0.00207535, 0.00127714, 0.000957854, 0.00207535, 
               0.000957854, 0.00143678, 0.00175607, 0.00223499, 0.00175607, 
               0.0011175, 0.00127714, 0.00127714, 0.00143678, 0.000798212, 
               0.0011175, 0.000159642, 0.000478927, 0.000478927, 0.000957854, 
               0.00127714, 0.000798212, 0.000478927, 0.000319285, 0.000319285, 
               0.00063857, 0.000159642, 0.000478927, 0.000478927, 0.000478927, 
               0.00063857, 0.000319285, 0.0011175, 0.000159642, 0.000159642, 
               0.000478927, 0.000319285, 0.000159642, 0.000159642, 0.000159642, 
               0.000159642, 0.000798212, 0.000319285, 0.000319285, 0.000159642, 
               0.000319285, 0.000159642, 0.000159642, 0, 0.000478927, 
               0, 0, 0.000159642, 0.000319285, 0.000319285, 
               0.000478927, 0.000159642, 0.000319285, 0, 0.000159642, 
               0.000159642, 0.000319285, 0.000319285, 0.000159642, 0.000159642, 
               0.000159642, 0.000319285, 0.000159642, 0, 0.000319285, 
               0.000478927, 0.000319285, 0.000159642, 0, 0.000159642, 
               0.000319285, 0.000159642, 0.000159642, 0.000159642, 0.000159642, 
               0.000159642, 0.000319285, 0, 0.000159642, 0.000319285, 
               0.000159642, 0, 0.000319285, 0, 0.000478927, 
               0.000319285, 0, 0.000159642, 0, 0.000319285, 
               0, 0, 0.000159642, 0.000319285, 0, 
               0.000478927, 0.000159642, 0, 0.000159642, 0.000159642, 
               0, 0, 0.000159642, 0.000319285, 0, 
               0, 0.000319285, 0, 0.000319285, 0, 
               0, 0, 0, 0.000159642, 0, 
               0.000159642, 0.000159642, 0, 0, 0, 
               0, 0, 0.000159642, 0.000159642, 0, 
               0, 0, 0.000159642, 0, 0, 
               0.000159642, 0, 0, 0.000159642, 0, 
               0.000159642, 0, 0, 0, 0, 
               0.000159642, 0, 0, 0, 0, 
               0.000159642, 0, 0, 0, 0, 
               0, 0, 0.000159642, 0, 0, 
               0, 0.000159642, 0, 0, 0, 
               0, 0, 0, 0, 0.000159642, 
               0, 0, 0.000159642, 0, 0, 
               0.000159642, 0, 0, 0, 0
            }, 
            {
               0, 0, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0, 0.000120773, 0, 0.000120773, 0, 
               0.000120773, 0.000362319, 0.000483092, 0.000724638, 0.00120773, 
               0.0013285, 0.0031401, 0.00326087, 0.00845411, 0.00809179, 
               0.0141304, 0.018599, 0.0218599, 0.021256, 0.026087, 
               0.0288647, 0.0318841, 0.0342995, 0.0349034, 0.0369565, 
               0.0359903, 0.0363527, 0.035628, 0.0347826, 0.0375604, 
               0.0369565, 0.0324879, 0.0321256, 0.0306763, 0.0298309, 
               0.0277778, 0.0240338, 0.0225845, 0.0183575, 0.0176329, 
               0.0210145, 0.0187198, 0.0150966, 0.0146135, 0.0111111, 
               0.0119565, 0.0129227, 0.00917874, 0.00954106, 0.00821256, 
               0.00688406, 0.00785024, 0.00688406, 0.00603865, 0.00615942, 
               0.00434783, 0.00422705, 0.00410628, 0.00338164, 0.00301932, 
               0.0044686, 0.00301932, 0.00181159, 0.00289855, 0.00229469, 
               0.00181159, 0.00241546, 0.00217391, 0.0013285, 0.00169082, 
               0.000845411, 0.00108696, 0.00108696, 0.000483092, 0.00205314, 
               0.00108696, 0.0013285, 0.000241546, 0.000845411, 0.000966184, 
               0.000724638, 0.0013285, 0.000241546, 0.000603865, 0.000724638, 
               0.00169082, 0.000483092, 0.000483092, 0.000724638, 0.000603865, 
               0.000724638, 0.000483092, 0.000362319, 0.000483092, 0.000483092, 
               0.000845411, 0.000362319, 0.000362319, 0.000724638, 0.000966184, 
               0.000966184, 0.000362319, 0.000120773, 0.000483092, 0.000362319, 
               0.000241546, 0.000724638, 0.000241546, 0.000362319, 0.000483092, 
               0, 0.000483092, 0.000362319, 0, 0, 
               0.000603865, 0.000483092, 0.000120773, 0.000362319, 0, 
               0.000241546, 0.000241546, 0.000120773, 0.000120773, 0.000120773, 
               0.000120773, 0.000241546, 0.000362319, 0.000241546, 0.000120773, 
               0.000483092, 0.000120773, 0.000362319, 0.000362319, 0.000120773, 
               0.000120773, 0.000362319, 0, 0, 0.000120773, 
               0.000603865, 0.000241546, 0.000362319, 0.000120773, 0, 
               0.000120773, 0.000120773, 0.000241546, 0.000120773, 0, 
               0.000120773, 0, 0.000120773, 0.000120773, 0, 
               0.000120773, 0, 0.000241546, 0.000120773, 0.000120773, 
               0.000120773, 0.000362319, 0.000120773, 0.000241546, 0, 
               0, 0, 0, 0, 0.000120773, 
               0, 0, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0, 0.000120773, 0.000120773, 0, 0, 
               0.000120773, 0, 0, 0.000120773, 0, 
               0, 0, 0.000120773, 0, 0, 
               0.000120773, 0.000120773, 0, 0, 0.000241546, 
               0, 0.000120773, 0, 0.000241546, 0.000120773, 
               0, 0, 0, 0, 0, 
               0.000241546, 0, 0, 0, 0, 
               0.000241546, 0, 0, 0, 0, 
               0.000120773, 0, 0, 0, 0, 
               0, 0.000120773, 0, 0.000120773, 0, 
               0, 0, 0, 0, 0, 
               0, 0, 0, 0, 0
            }, 
            {
               0, 0, 0, 0, 0.000235073, 
               0.000117536, 0, 0, 0, 0.000117536, 
               0, 0.000235073, 0, 0.000117536, 0.000235073, 
               0.000117536, 0.000117536, 0.000352609, 0.000117536, 0.000235073, 
               0.00105783, 0.00141044, 0.00152797, 0.00376117, 0.00564175, 
               0.00693465, 0.0109309, 0.0105783, 0.0188058, 0.0185708, 
               0.0256229, 0.0245651, 0.0303244, 0.030677, 0.0332628, 
               0.0404325, 0.0342031, 0.0353785, 0.0349083, 0.0378467, 
               0.0352609, 0.0338505, 0.0343206, 0.0326751, 0.0330277, 
               0.0293841, 0.0287964, 0.0278561, 0.0235073, 0.0224495, 
               0.0209215, 0.0209215, 0.01622, 0.0183357, 0.0152797, 
               0.01622, 0.0124589, 0.0131641, 0.0112835, 0.00987306, 
               0.0102257, 0.00822755, 0.00716972, 0.00811001, 0.0064645, 
               0.00681711, 0.0051716, 0.00423131, 0.0038787, 0.00423131, 
               0.00434885, 0.00329102, 0.00235073, 0.00317348, 0.00235073, 
               0.00211566, 0.00211566, 0.00211566, 0.0012929, 0.00164551, 
               0.0012929, 0.0012929, 0.000940291, 0.0012929, 0.00176305, 
               0.000705219, 0.00141044, 0.000705219, 0.00105783, 0.00117536, 
               0.00105783, 0.000822755, 0.000587682, 0.000822755, 0.00117536, 
               0.000235073, 0.000470146, 0.000470146, 0.000470146, 0.000587682, 
               0.000470146, 0.000587682, 0.000352609, 0.000940291, 0.000352609, 
               0.000587682, 0.000352609, 0.000470146, 0.000470146, 0.000352609, 
               0.000352609, 0.000470146, 0.000235073, 0.000235073, 0.000587682, 
               0.000117536, 0.000235073, 0.000352609, 0.000235073, 0.000235073, 
               0.000235073, 0.000235073, 0.000352609, 0.000352609, 0.000470146, 
               0.000470146, 0.000705219, 0.000470146, 0.000470146, 0.000117536, 
               0.000470146, 0.000235073, 0.000352609, 0, 0.000235073, 
               0.000235073, 0.000470146, 0.000117536, 0, 0.000352609, 
               0, 0.000117536, 0, 0.000470146, 0, 
               0.000235073, 0.000352609, 0.000235073, 0.000235073, 0, 
               0.000235073, 0.000117536, 0, 0, 0.000117536, 
               0.000352609, 0.000470146, 0, 0.000587682, 0, 
               0, 0.000352609, 0.000352609, 0.000352609, 0.000235073, 
               0.000235073, 0.000235073, 0.000352609, 0.000235073, 0.000235073, 
               0.000117536, 0, 0.000235073, 0, 0.000352609, 
               0, 0.000117536, 0.000117536, 0.000352609, 0, 
               0, 0, 0.000235073, 0.000117536, 0, 
               0, 0.000117536, 0, 0, 0, 
               0, 0, 0.000235073, 0, 0, 
               0, 0, 0, 0.000117536, 0.000117536, 
               0, 0, 0.000117536, 0.000117536, 0.000117536, 
               0, 0.000117536, 0.000117536, 0.000117536, 0, 
               0.000235073, 0, 0.000117536, 0, 0, 
               0, 0.000117536, 0.000117536, 0, 0, 
               0, 0.000117536, 0.000117536, 0.000117536, 0.000117536, 
               0, 0.000235073, 0.000117536, 0.000117536, 0, 
               0.000117536, 0, 0, 0, 0.000117536, 
               0, 0, 0.000117536, 0.000117536, 0, 
               0, 0, 0, 0, 0, 
               0, 0, 0.000117536, 0, 0
            }, 
            {
               0, 0, 0, 0, 0, 
               0, 0, 0, 0, 0.000107112, 
               0, 0.000321337, 0, 0, 0.000214225, 
               0.000214225, 0, 0, 0.000107112, 0.000214225, 
               0.000856898, 0.000749786, 0.00160668, 0.00214225, 0.0035347, 
               0.00739075, 0.00771208, 0.0108183, 0.0138175, 0.0153171, 
               0.0194944, 0.0220651, 0.0250643, 0.0292416, 0.0302057, 
               0.0323479, 0.0397386, 0.0345973, 0.0366324, 0.0347044, 
               0.0393102, 0.0368466, 0.0351328, 0.0315981, 0.0320266, 
               0.0292416, 0.0293488, 0.0298843, 0.0244216, 0.0254927, 
               0.0247429, 0.0224936, 0.0214225, 0.018102, 0.0147815, 
               0.0128535, 0.0152099, 0.0145673, 0.0121037, 0.0123179, 
               0.00996144, 0.0103899, 0.00814053, 0.00846187, 0.00664096, 
               0.00803342, 0.00685518, 0.00685518, 0.0052485, 0.0043916, 
               0.00589117, 0.00342759, 0.00299914, 0.00278492, 0.00299914, 
               0.00321337, 0.00299914, 0.00289203, 0.00278492, 0.00160668, 
               0.00117823, 0.0017138, 0.00289203, 0.0017138, 0.000749786, 
               0.00139246, 0.00107112, 0.00139246, 0.000642674, 0.000642674, 
               0.000428449, 0.000428449, 0.000535561, 0.000321337, 0.000428449, 
               0.000642674, 0.000749786, 0.000642674, 0.000856898, 0.000214225, 
               0.000428449, 0.000214225, 0.000749786, 0.000214225, 0.000642674, 
               0.000428449, 0.000535561, 0.000428449, 0.000749786, 0.000214225, 
               0.000535561, 0.000535561, 0.000535561, 0.000428449, 0.000856898, 
               0.000535561, 0.000321337, 0.000214225, 0.000535561, 0.000321337, 
               0.000107112, 0.000107112, 0.000107112, 0.000642674, 0.000107112, 
               0.000321337, 0.000214225, 0.000214225, 0, 0.000321337, 
               0.000107112, 0.000428449, 0.000214225, 0.000107112, 0.000214225, 
               0.000321337, 0.000214225, 0.000107112, 0.000107112, 0.000321337, 
               0, 0.000214225, 0.000321337, 0.000107112, 0, 
               0.000214225, 0.000107112, 0.000214225, 0.000321337, 0.000428449, 
               0.000321337, 0.000214225, 0.000107112, 0, 0, 
               0.000321337, 0.000107112, 0.000321337, 0.000107112, 0, 
               0.000321337, 0, 0, 0.000214225, 0.000107112, 
               0.000321337, 0.000107112, 0.000214225, 0.000214225, 0, 
               0, 0.000107112, 0, 0, 0, 
               0, 0.000321337, 0.000214225, 0, 0, 
               0, 0, 0, 0.000214225, 0, 
               0, 0, 0.000107112, 0, 0.000107112, 
               0, 0.000107112, 0.000214225, 0.000107112, 0, 
               0.000107112, 0, 0, 0, 0.000107112, 
               0, 0.000107112, 0, 0.000214225, 0.000107112, 
               0.000107112, 0, 0.000107112, 0.000107112, 0, 
               0, 0, 0.000107112, 0, 0, 
               0, 0, 0, 0.000107112, 0, 
               0, 0, 0.000107112, 0.000321337, 0, 
               0, 0.000214225, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0.000107112, 0, 0, 0.000214225, 0, 
               0, 0.000107112, 0.000107112, 0, 0
            }, 
            {
               0, 0, 0, 0, 0, 
               0, 0, 0.00011558, 0, 0, 
               0, 0, 0, 0.00011558, 0.00011558, 
               0.00011558, 0.00011558, 0.000346741, 0.00023116, 0.00023116, 
               0.000346741, 0.000924642, 0.00127138, 0.00127138, 0.00300509, 
               0.00369857, 0.00612575, 0.00705039, 0.00809061, 0.0122515, 
               0.0162968, 0.0177994, 0.0238095, 0.0271613, 0.0309755, 
               0.0293574, 0.0308599, 0.0381415, 0.0365233, 0.0340962, 
               0.0351364, 0.0328248, 0.0336338, 0.0357143, 0.0330559, 
               0.0329404, 0.0292418, 0.0275081, 0.0286639, 0.0269302, 
               0.0247342, 0.0206889, 0.02092, 0.0206889, 0.0184928, 
               0.0181461, 0.018724, 0.0154877, 0.0150254, 0.0128294, 
               0.0149098, 0.00901526, 0.0105178, 0.00751271, 0.00947758, 
               0.00728155, 0.00520111, 0.00866852, 0.00554785, 0.00658807, 
               0.00566343, 0.00358299, 0.00473879, 0.00312067, 0.00473879, 
               0.00358299, 0.00277393, 0.0023116, 0.00219602, 0.00127138, 
               0.00242718, 0.00138696, 0.00196486, 0.0011558, 0.00127138, 
               0.00127138, 0.00219602, 0.00138696, 0.000462321, 0.0011558, 
               0.000924642, 0.00104022, 0.00127138, 0.00196486, 0.000577901, 
               0.000693481, 0.000693481, 0.000577901, 0.000809061, 0.000346741, 
               0.000693481, 0.000693481, 0.000577901, 0.0011558, 0.000577901, 
               0.000346741, 0.000577901, 0.000462321, 0.000346741, 0.000809061, 
               0.000346741, 0.00023116, 0.00011558, 0.000346741, 0.000577901, 
               0.00011558, 0.000346741, 0.00023116, 0.000346741, 0.000346741, 
               0.00011558, 0.00011558, 0.000577901, 0.00011558, 0.000462321, 
               0.000809061, 0.00023116, 0.00011558, 0.00011558, 0.00011558, 
               0.00011558, 0.00011558, 0.000346741, 0.00023116, 0, 
               0.000462321, 0.00011558, 0, 0.000346741, 0.00023116, 
               0.00023116, 0.000577901, 0.00011558, 0.00023116, 0, 
               0.000346741, 0.00023116, 0.00023116, 0, 0.00023116, 
               0.00011558, 0.000346741, 0.00011558, 0.00011558, 0, 
               0.00023116, 0.00023116, 0.00023116, 0.00011558, 0.00011558, 
               0.00011558, 0.00011558, 0.00023116, 0, 0.00023116, 
               0, 0.00011558, 0.00011558, 0, 0, 
               0.00011558, 0.00011558, 0.00011558, 0.00023116, 0.000346741, 
               0.000346741, 0.00011558, 0, 0.00023116, 0.00011558, 
               0.00011558, 0, 0, 0.00023116, 0, 
               0.00011558, 0, 0, 0, 0, 
               0, 0.00011558, 0, 0, 0.00011558, 
               0.00023116, 0.00011558, 0, 0, 0.00011558, 
               0.00011558, 0, 0, 0, 0, 
               0.00011558, 0.00011558, 0.00011558, 0, 0, 
               0.00011558, 0.00011558, 0.00023116, 0, 0, 
               0, 0.00023116, 0, 0, 0.00011558, 
               0.00011558, 0.00011558, 0.00011558, 0, 0, 
               0, 0, 0.00011558, 0, 0, 
               0, 0, 0.00011558, 0, 0, 
               0, 0, 0, 0, 0.00011558, 
               0, 0.00011558, 0, 0, 0.00011558, 
               0, 0, 0, 0, 0
            }, 
            {
               0, 0, 0, 0.00010582, 0.00010582, 
               0, 0, 0.00010582, 0, 0, 
               0, 0.00010582, 0, 0, 0, 
               0.00010582, 0, 0, 0.00021164, 0.00021164, 
               0.00010582, 0.000529101, 0.00010582, 0.00169312, 0.00116402, 
               0.00169312, 0.00465608, 0.00560847, 0.00878307, 0.0114286, 
               0.0122751, 0.0183069, 0.0201058, 0.0207407, 0.0266667, 
               0.0273016, 0.0316402, 0.0283598, 0.0331217, 0.0380952, 
               0.0342857, 0.0362963, 0.0361905, 0.035873, 0.0356614, 
               0.0320635, 0.0330159, 0.0326984, 0.0288889, 0.0251852, 
               0.0245503, 0.0274074, 0.0249735, 0.0197884, 0.0201058, 
               0.0180952, 0.0201058, 0.0169312, 0.0156614, 0.0143915, 
               0.0126984, 0.0116402, 0.00867725, 0.0108995, 0.0108995, 
               0.00814815, 0.00867725, 0.00740741, 0.00592593, 0.00497354, 
               0.00582011, 0.00497354, 0.00560847, 0.00412698, 0.00380952, 
               0.0031746, 0.00275132, 0.00232804, 0.00296296, 0.00232804, 
               0.00179894, 0.00190476, 0.00179894, 0.00275132, 0.0021164, 
               0.0010582, 0.00116402, 0.00148148, 0.00137566, 0.00137566, 
               0.00116402, 0.000846561, 0.00031746, 0.000529101, 0.000846561, 
               0.0010582, 0.00042328, 0.000634921, 0.000740741, 0.00031746, 
               0.000529101, 0.00021164, 0.000634921, 0.00021164, 0.00021164, 
               0.00042328, 0.00010582, 0.000529101, 0.00010582, 0.000634921, 
               0.00021164, 0.00042328, 0.00031746, 0.000529101, 0.000952381, 
               0.00021164, 0, 0.00031746, 0.00010582, 0.00031746, 
               0.00042328, 0.00031746, 0.00031746, 0.000529101, 0.00010582, 
               0.000529101, 0.000529101, 0.000529101, 0.00010582, 0.00021164, 
               0.00010582, 0.00021164, 0.00010582, 0.00010582, 0.00021164, 
               0.00010582, 0.00031746, 0.00021164, 0.00021164, 0, 
               0.000529101, 0, 0, 0.00010582, 0.00010582, 
               0.000634921, 0, 0.00010582, 0, 0.00031746, 
               0, 0.00010582, 0, 0, 0, 
               0.00031746, 0.00010582, 0.00021164, 0.000529101, 0.00021164, 
               0.00021164, 0.00010582, 0.00010582, 0, 0.00010582, 
               0.00021164, 0, 0, 0, 0, 
               0.00010582, 0.00021164, 0.00031746, 0, 0.00021164, 
               0.00010582, 0.00010582, 0, 0, 0, 
               0.00010582, 0.00010582, 0.00010582, 0, 0, 
               0, 0.00021164, 0.00010582, 0.00021164, 0.00021164, 
               0.00021164, 0, 0, 0.00010582, 0.00010582, 
               0, 0, 0.00021164, 0, 0, 
               0, 0.00010582, 0.00021164, 0.00010582, 0.00010582, 
               0, 0, 0, 0.00010582, 0.00010582, 
               0.00010582, 0, 0.00021164, 0, 0, 
               0.00010582, 0.00010582, 0, 0.00010582, 0, 
               0, 0, 0, 0, 0, 
               0, 0, 0, 0.00010582, 0, 
               0.00010582, 0.00031746, 0, 0, 0, 
               0, 0.00010582, 0, 0, 0.00010582, 
               0, 0.00010582, 0, 0, 0, 
               0, 0, 0, 0, 0.00010582
            }
          };

  Double_t probEL[kNo_Mom][kNo_EnBins]={  
            {
               0, 0, 0, 0, 0, 
               0.000141603, 0, 0, 0, 0, 
               0, 0.000141603, 0.000283206, 0.000141603, 0, 
               0, 0, 0.000141603, 0.000141603, 0, 
               0, 0.000141603, 0, 0, 0.000141603, 
               0.000141603, 0.000708015, 0.000424809, 0.000849618, 0.000849618, 
               0.00155763, 0.00240725, 0.000991221, 0.00354007, 0.00580572, 
               0.00608893, 0.00736335, 0.0107618, 0.011045, 0.0104786, 
               0.0127443, 0.0184084, 0.0167091, 0.0155763, 0.0209572, 
               0.01855, 0.0239309, 0.0225149, 0.0244973, 0.0252053, 
               0.0247805, 0.0242141, 0.0291702, 0.0237893, 0.0278958, 
               0.0226565, 0.026763, 0.0230813, 0.0223733, 0.0257717, 
               0.0230813, 0.0215236, 0.0210988, 0.0198244, 0.0196828, 
               0.0191164, 0.0175588, 0.0179836, 0.0198244, 0.017134, 
               0.0155763, 0.0128859, 0.0147267, 0.0130275, 0.0127443, 
               0.0128859, 0.0120363, 0.00764656, 0.00877938, 0.00991221, 
               0.00920419, 0.00778816, 0.00708015, 0.00623053, 0.00679694, 
               0.00538091, 0.00566412, 0.00693854, 0.0049561, 0.00538091, 
               0.00566412, 0.00424809, 0.00438969, 0.00523931, 0.00254885, 
               0.0046729, 0.00354007, 0.00254885, 0.00311526, 0.00269046, 
               0.00226565, 0.00354007, 0.00212404, 0.00212404, 0.00254885, 
               0.00226565, 0.00240725, 0.00184084, 0.00226565, 0.00198244, 
               0.000991221, 0.00113282, 0.000566412, 0.00127443, 0.00155763, 
               0.00127443, 0.000849618, 0.00113282, 0.000849618, 0.000566412, 
               0.00141603, 0.00113282, 0.00155763, 0.000566412, 0.00113282, 
               0.000991221, 0.000566412, 0.000424809, 0.000424809, 0.000283206, 
               0.000849618, 0.000283206, 0.000991221, 0.000141603, 0.000566412, 
               0.000708015, 0.000708015, 0.000283206, 0.000283206, 0.000424809, 
               0.000283206, 0.000708015, 0.000424809, 0.000708015, 0.000141603, 
               0.000566412, 0, 0.000708015, 0.000283206, 0.000283206, 
               0.000141603, 0.000566412, 0.000424809, 0, 0.000424809, 
               0.000141603, 0, 0.000283206, 0.000424809, 0.000283206, 
               0.000141603, 0.000283206, 0.000141603, 0, 0.000708015, 
               0.000283206, 0, 0.000283206, 0.000141603, 0.000141603, 
               0.000141603, 0.000283206, 0, 0.000283206, 0.000141603, 
               0.000141603, 0.000141603, 0.000141603, 0.000424809, 0.000141603, 
               0.000141603, 0.000283206, 0.000141603, 0, 0.000283206, 
               0.000141603, 0.000141603, 0.000424809, 0, 0.000141603, 
               0.000283206, 0, 0.000283206, 0, 0, 
               0, 0, 0.000283206, 0.000283206, 0.000283206, 
               0, 0, 0.000141603, 0, 0.000141603, 
               0, 0, 0, 0, 0.000141603, 
               0.000141603, 0, 0.000141603, 0.000141603, 0.000141603, 
               0, 0, 0.000141603, 0, 0, 
               0, 0, 0.000283206, 0, 0, 
               0, 0, 0, 0, 0.000141603, 
               0, 0, 0, 0, 0, 
               0, 0.000141603, 0, 0, 0, 
               0, 0.000141603, 0, 0, 0, 
               0.000141603, 0, 0, 0.000141603, 0
            }, 
            {
               0, 0, 0, 0, 0, 
               0, 0.000138427, 0.000138427, 0, 0, 
               0, 0.000138427, 0.00055371, 0, 0, 
               0, 0, 0.000138427, 0, 0.000276855, 
               0.000276855, 0, 0.000138427, 0.00055371, 0.000138427, 
               0, 0.000276855, 0.000276855, 0.00055371, 0.00124585, 
               0.00110742, 0.00221484, 0.00373754, 0.0030454, 0.00484496, 
               0.00609081, 0.00567553, 0.00733666, 0.00816722, 0.01232, 
               0.0110742, 0.0160576, 0.0138427, 0.0160576, 0.0159192, 
               0.0186877, 0.021041, 0.021041, 0.0184109, 0.0221484, 
               0.0202104, 0.019103, 0.020072, 0.0236711, 0.0211794, 
               0.0235327, 0.0213178, 0.021041, 0.0224252, 0.0235327, 
               0.0227021, 0.0164729, 0.0218715, 0.0213178, 0.0184109, 
               0.0203488, 0.0202104, 0.0164729, 0.0135659, 0.0145349, 
               0.0153654, 0.014258, 0.0146733, 0.0163344, 0.0141196, 
               0.0130122, 0.0141196, 0.0131506, 0.0102436, 0.00899779, 
               0.00968992, 0.0101052, 0.0127353, 0.0095515, 0.0102436, 
               0.0095515, 0.00927464, 0.00816722, 0.00636766, 0.0055371, 
               0.00609081, 0.00636766, 0.00484496, 0.00498339, 0.00442968, 
               0.00263012, 0.0055371, 0.00442968, 0.00332226, 0.00359911, 
               0.00318383, 0.00373754, 0.00373754, 0.00193798, 0.00263012, 
               0.00263012, 0.00332226, 0.00276855, 0.00221484, 0.00207641, 
               0.00346069, 0.00166113, 0.00207641, 0.00221484, 0.00166113, 
               0.00207641, 0.00166113, 0.00193798, 0.00138427, 0.00166113, 
               0.000968992, 0.000692137, 0.00124585, 0.00124585, 0.0015227, 
               0.00207641, 0.0015227, 0.0015227, 0.00124585, 0.000830565, 
               0.000415282, 0.000692137, 0.00055371, 0.000692137, 0.00055371, 
               0.000968992, 0.000968992, 0.00110742, 0.000968992, 0.000692137, 
               0.000415282, 0.000692137, 0.000968992, 0.00055371, 0.000830565, 
               0.000415282, 0.000415282, 0.000968992, 0.000968992, 0.000415282, 
               0.000692137, 0.000276855, 0.000138427, 0.00055371, 0.000138427, 
               0.00055371, 0.000415282, 0.000276855, 0.000138427, 0, 
               0.00055371, 0.000138427, 0.000415282, 0.000415282, 0.000415282, 
               0.000968992, 0.000276855, 0.000415282, 0.000415282, 0.000138427, 
               0.000138427, 0.000415282, 0.000276855, 0.000138427, 0, 
               0.000276855, 0, 0.000138427, 0.000276855, 0.000830565, 
               0.000415282, 0.000415282, 0, 0.000138427, 0.000138427, 
               0, 0.000276855, 0.000138427, 0.000138427, 0, 
               0.000138427, 0.000138427, 0.000276855, 0, 0.000138427, 
               0, 0.000138427, 0, 0.000138427, 0.000138427, 
               0, 0.000138427, 0.000138427, 0.000138427, 0, 
               0, 0, 0, 0, 0, 
               0.000138427, 0, 0.000138427, 0, 0, 
               0, 0.000138427, 0.000276855, 0.000138427, 0, 
               0, 0.000276855, 0, 0.000138427, 0.000138427, 
               0.000138427, 0, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0.000276855, 0.000138427, 0, 0, 0, 
               0, 0, 0, 0, 0, 
               0, 0, 0.000276855, 0.000276855, 0
            }, 
            {
               0, 0, 0.000228781, 0, 0.00011439, 
               0, 0, 0.00011439, 0, 0.00011439, 
               0, 0, 0.000228781, 0.00011439, 0, 
               0, 0, 0.00011439, 0, 0.00011439, 
               0, 0.000343171, 0, 0.000228781, 0.00011439, 
               0, 0.000228781, 0.000800732, 0.000571951, 0.000686342, 
               0.00102951, 0.00125829, 0.0024022, 0.00205903, 0.00537634, 
               0.00514756, 0.00583391, 0.00754976, 0.00674903, 0.0117822, 
               0.0133837, 0.0108671, 0.0104095, 0.0122398, 0.0153283, 
               0.0179593, 0.0217342, 0.0171585, 0.0186456, 0.0181881, 
               0.0232212, 0.021391, 0.0237932, 0.0224205, 0.0226493, 
               0.0256234, 0.0202471, 0.0194464, 0.0235644, 0.0186456, 
               0.021391, 0.0195607, 0.0210478, 0.0202471, 0.0178449, 
               0.0155571, 0.0220773, 0.0171585, 0.0173873, 0.0175017, 
               0.0173873, 0.0152139, 0.0141844, 0.0147563, 0.0136124, 
               0.0134981, 0.0128117, 0.0125829, 0.0101807, 0.0112102, 
               0.00972318, 0.0101807, 0.00869366, 0.00926561, 0.0082361, 
               0.00926561, 0.0070922, 0.00663464, 0.00743537, 0.00629147, 
               0.00789293, 0.00846488, 0.00732098, 0.00571951, 0.00526195, 
               0.0059483, 0.0059483, 0.00549073, 0.00526195, 0.00491878, 
               0.0024022, 0.00343171, 0.00423244, 0.00400366, 0.00469, 
               0.00388927, 0.00308854, 0.00331732, 0.00320293, 0.00251659, 
               0.00205903, 0.00251659, 0.00205903, 0.00171585, 0.00171585, 
               0.00148707, 0.00263098, 0.00137268, 0.00274537, 0.0011439, 
               0.00183024, 0.00137268, 0.00102951, 0.000800732, 0.00102951, 
               0.000571951, 0.00160146, 0.00194464, 0.00137268, 0.0011439, 
               0.00148707, 0.000915122, 0.000800732, 0.00137268, 0.000800732, 
               0.000686342, 0.000457561, 0.000800732, 0.000571951, 0.0011439, 
               0.000343171, 0.000800732, 0.000800732, 0.000571951, 0.000686342, 
               0.000343171, 0.000457561, 0.000457561, 0.000571951, 0.000571951, 
               0.000571951, 0.000457561, 0.000571951, 0.00011439, 0.000457561, 
               0.000571951, 0.000571951, 0.000228781, 0.000343171, 0.000228781, 
               0.000686342, 0.000457561, 0.00011439, 0.000228781, 0.000571951, 
               0.000571951, 0.00011439, 0.000228781, 0.000571951, 0.000343171, 
               0.00011439, 0.000228781, 0.000228781, 0.000228781, 0.000228781, 
               0.000228781, 0.000457561, 0.000343171, 0.00011439, 0.00011439, 
               0, 0.00011439, 0, 0.000228781, 0.000343171, 
               0, 0.000228781, 0.00011439, 0, 0, 
               0.00011439, 0.00011439, 0.00011439, 0, 0.000228781, 
               0.00011439, 0, 0.000343171, 0.00011439, 0.000228781, 
               0.00011439, 0.00011439, 0, 0.00011439, 0, 
               0.000228781, 0.000343171, 0, 0.00011439, 0.00011439, 
               0, 0.00011439, 0.000343171, 0.00011439, 0.00011439, 
               0, 0.00011439, 0.00011439, 0, 0, 
               0, 0, 0.000343171, 0, 0, 
               0, 0, 0, 0.00011439, 0, 
               0, 0, 0.00011439, 0.00011439, 0, 
               0, 0.00011439, 0, 0.00011439, 0, 
               0.000228781, 0, 0, 0.00011439, 0.00011439, 
               0, 0, 0.000228781, 0, 0
            }, 
            {
               0, 0, 0, 0.00010144, 0, 
               0.00010144, 0.00010144, 0, 0, 0.00010144, 
               0.000202881, 0, 0.00010144, 0.00010144, 0.000202881, 
               0.00010144, 0, 0.000304321, 0.000405762, 0.00010144, 
               0, 0.000202881, 0.000202881, 0.000507202, 0.000304321, 
               0, 0.000202881, 0.000405762, 0.000811524, 0.000507202, 
               0.00142017, 0.00152161, 0.00273889, 0.00223169, 0.00334753, 
               0.00436194, 0.00547778, 0.00618787, 0.00821668, 0.0095354, 
               0.0095354, 0.0118685, 0.0121729, 0.0121729, 0.014506, 
               0.0142017, 0.0165348, 0.0185636, 0.0177521, 0.0202881, 
               0.0185636, 0.0192737, 0.0208967, 0.0200852, 0.0222155, 
               0.0185636, 0.0181578, 0.0194766, 0.0207953, 0.0201867, 
               0.0180564, 0.017955, 0.017955, 0.0184622, 0.0191722, 
               0.0174478, 0.0174478, 0.017955, 0.0157233, 0.0155204, 
               0.017955, 0.0141002, 0.0164334, 0.0155204, 0.017042, 
               0.0146074, 0.0131873, 0.0141002, 0.0104484, 0.0129844, 
               0.0121729, 0.0114628, 0.0113613, 0.00963684, 0.00872388, 
               0.0108541, 0.00841956, 0.00892676, 0.00811524, 0.00831812, 
               0.00507202, 0.00811524, 0.00730371, 0.00689795, 0.00689795, 
               0.00730371, 0.00497058, 0.00568067, 0.00466626, 0.00486914, 
               0.00466626, 0.00405762, 0.00507202, 0.00355042, 0.00415906, 
               0.00284033, 0.00385474, 0.00334753, 0.00334753, 0.00243457, 
               0.00334753, 0.00344898, 0.00385474, 0.00182593, 0.00243457, 
               0.00223169, 0.00263745, 0.00273889, 0.00213025, 0.00213025, 
               0.00243457, 0.00152161, 0.00152161, 0.00162305, 0.00172449, 
               0.00111584, 0.00142017, 0.00152161, 0.00111584, 0.00162305, 
               0.00121729, 0.00111584, 0.000507202, 0.00121729, 0.000710083, 
               0.0010144, 0.000405762, 0.00121729, 0.000507202, 0.00131873, 
               0.0010144, 0.000608643, 0.000710083, 0.000507202, 0.000710083, 
               0.000608643, 0.000507202, 0.000507202, 0.000304321, 0.000608643, 
               0.000912964, 0.000608643, 0.000710083, 0.000405762, 0.000507202, 
               0.000811524, 0.000912964, 0.000507202, 0.000304321, 0.000304321, 
               0.000405762, 0.000507202, 0.000811524, 0.000405762, 0.000507202, 
               0.000608643, 0.00010144, 0.000405762, 0.000202881, 0.000304321, 
               0.000608643, 0.000608643, 0.00010144, 0.000507202, 0.00010144, 
               0.000304321, 0.000507202, 0.000304321, 0.000202881, 0.000304321, 
               0.000405762, 0.000304321, 0.000304321, 0.000405762, 0.000202881, 
               0.00010144, 0.000507202, 0, 0.000405762, 0, 
               0.000507202, 0.000405762, 0.000304321, 0.000405762, 0, 
               0.000507202, 0.000304321, 0.000405762, 0.000304321, 0.000507202, 
               0.000202881, 0.00010144, 0.000202881, 0, 0.000405762, 
               0.000202881, 0.000202881, 0.000405762, 0, 0.000202881, 
               0.000202881, 0, 0.000202881, 0.000202881, 0.000202881, 
               0.00010144, 0, 0.00010144, 0.000202881, 0.000202881, 
               0.000202881, 0.00010144, 0.000202881, 0, 0, 
               0.00010144, 0.00010144, 0.00010144, 0.00010144, 0, 
               0, 0, 0, 0, 0, 
               0, 0.00010144, 0, 0, 0, 
               0.00010144, 0.000202881, 0.00010144, 0.00010144, 0, 
               0, 0.00010144, 0, 0.00010144, 0.00010144
            }, 
            {
               0, 0, 0, 0, 0, 
               0, 0, 0, 0.000316056, 0.000105352, 
               0.000210704, 0.000210704, 0.000316056, 0, 0, 
               0, 0.000105352, 0.000210704, 0.000105352, 0.000105352, 
               0.000210704, 0.000105352, 0, 0.000316056, 0.000105352, 
               0.000632111, 0.000210704, 0.000526759, 0.000737463, 0.00115887, 
               0.00105352, 0.00115887, 0.00147493, 0.00273915, 0.00337126, 
               0.0054783, 0.00495154, 0.00579435, 0.00758534, 0.00790139, 
               0.00842815, 0.00937632, 0.0106405, 0.013169, 0.0144332, 
               0.0142225, 0.015276, 0.0168563, 0.0150653, 0.0172777, 
               0.0184366, 0.0175938, 0.0189633, 0.0200169, 0.0198062, 
               0.0185419, 0.0187526, 0.0189633, 0.018858, 0.0198062, 
               0.0202276, 0.0166456, 0.0164349, 0.0198062, 0.0207543, 
               0.0154867, 0.0179098, 0.0189633, 0.0169617, 0.0147493, 
               0.0169617, 0.0155921, 0.0161188, 0.0163295, 0.0148546, 
               0.0154867, 0.0126422, 0.0112727, 0.0139064, 0.0120101, 
               0.0130636, 0.0117994, 0.0085335, 0.0110619, 0.0102191, 
               0.00842815, 0.00906026, 0.0104298, 0.0085335, 0.00874421, 
               0.00716393, 0.00769069, 0.00821745, 0.00716393, 0.00737463, 
               0.00684787, 0.00495154, 0.00453013, 0.00579435, 0.00558365, 
               0.00558365, 0.00368732, 0.00453013, 0.00400337, 0.00442478, 
               0.00474083, 0.00463548, 0.00431943, 0.00379267, 0.00400337, 
               0.00421408, 0.00273915, 0.0028445, 0.0028445, 0.0030552, 
               0.0026338, 0.0028445, 0.00273915, 0.0026338, 0.00294985, 
               0.00168563, 0.00189633, 0.00158028, 0.00105352, 0.00179098, 
               0.00316056, 0.00168563, 0.00126422, 0.00221239, 0.00242309, 
               0.00136957, 0.00147493, 0.00136957, 0.000948167, 0.00168563, 
               0.00126422, 0.00158028, 0.00126422, 0.00105352, 0.000948167, 
               0.00147493, 0.000421408, 0.000948167, 0.000948167, 0.00115887, 
               0.000842815, 0.000210704, 0.00126422, 0.000105352, 0.000632111, 
               0.000421408, 0.000948167, 0.000526759, 0.000526759, 0.000421408, 
               0.000526759, 0.000421408, 0.000737463, 0.000526759, 0.000632111, 
               0.000210704, 0.000737463, 0.000526759, 0.000210704, 0.000632111, 
               0.000526759, 0.000632111, 0.000526759, 0.000316056, 0.000316056, 
               0.000526759, 0.000526759, 0.000105352, 0.000948167, 0.000526759, 
               0.000210704, 0.000632111, 0.000632111, 0.000421408, 0.000421408, 
               0.000421408, 0, 0.000316056, 0.000210704, 0, 
               0.000421408, 0.000421408, 0.000316056, 0.000105352, 0.000105352, 
               0.000316056, 0.000105352, 0.000105352, 0.000210704, 0.000210704, 
               0.000105352, 0.000210704, 0.000105352, 0, 0.000842815, 
               0.000105352, 0, 0, 0.000421408, 0.000105352, 
               0.000105352, 0.000105352, 0, 0.000210704, 0, 
               0.000105352, 0.000210704, 0.000105352, 0.000105352, 0.000210704, 
               0.000105352, 0, 0.000105352, 0.000421408, 0.000210704, 
               0.000210704, 0.000210704, 0, 0, 0.000210704, 
               0, 0.000105352, 0.000210704, 0, 0.000210704, 
               0, 0, 0.000105352, 0.000105352, 0, 
               0, 0, 0.000316056, 0.000105352, 0.000105352, 
               0, 0.000105352, 0, 0.000105352, 0, 
               0.000105352, 0, 0.000105352, 0, 0.000210704
            }, 
            {
               0, 0, 0, 0, 0.000105552, 
               0, 0.000105552, 0.000316656, 0.000105552, 0, 
               0, 0.000105552, 0.000105552, 0.000105552, 0.000211104, 
               0.000422208, 0.000422208, 0.000211104, 0.000105552, 0.000211104, 
               0.000211104, 0.000105552, 0.000316656, 0.000211104, 0.00052776, 
               0.00052776, 0.000316656, 0.000422208, 0.000422208, 0.00126662, 
               0.000738864, 0.00147773, 0.00221659, 0.00147773, 0.00295546, 
               0.00358877, 0.00369432, 0.00591091, 0.00749419, 0.00812751, 
               0.00992189, 0.0102385, 0.011083, 0.0112941, 0.0164661, 
               0.0150939, 0.017205, 0.0142495, 0.0137218, 0.0160439, 
               0.0184716, 0.0170994, 0.016255, 0.0184716, 0.0196327, 
               0.0180494, 0.0165717, 0.0189994, 0.0186827, 0.0168883, 
               0.0200549, 0.0197382, 0.0183661, 0.0192105, 0.016255, 
               0.0185772, 0.0183661, 0.018155, 0.0158328, 0.0166772, 
               0.0151995, 0.014144, 0.0155161, 0.0145662, 0.0151995, 
               0.0124551, 0.0132996, 0.0155161, 0.0128773, 0.0134051, 
               0.0117163, 0.0113996, 0.0102385, 0.00876082, 0.0107663, 
               0.0117163, 0.00802195, 0.00876082, 0.00918303, 0.0108719, 
               0.0105552, 0.00907748, 0.0077053, 0.00643867, 0.00717754, 
               0.00696643, 0.00675533, 0.00591091, 0.0050665, 0.00496095, 
               0.00474984, 0.00538315, 0.00443319, 0.00453874, 0.00443319, 
               0.00591091, 0.00390543, 0.00432763, 0.00422208, 0.00443319, 
               0.00327211, 0.00401098, 0.00401098, 0.00327211, 0.00221659, 
               0.00295546, 0.00369432, 0.00253325, 0.00274435, 0.00200549, 
               0.00274435, 0.00232214, 0.00253325, 0.00179438, 0.00147773, 
               0.00168883, 0.00158328, 0.00158328, 0.00158328, 0.00168883, 
               0.00200549, 0.000949968, 0.00137218, 0.00105552, 0.00147773, 
               0.000949968, 0.000844416, 0.000949968, 0.00105552, 0.000633312, 
               0.000422208, 0.000633312, 0.000949968, 0.00116107, 0.000738864, 
               0.000844416, 0.00052776, 0.000844416, 0.000422208, 0.00105552, 
               0.000844416, 0.000633312, 0.000738864, 0.00105552, 0.000949968, 
               0.00137218, 0.000633312, 0.00052776, 0.000633312, 0.000422208, 
               0.000738864, 0.000316656, 0.000211104, 0.000633312, 0.000949968, 
               0.000633312, 0.00052776, 0.000633312, 0.00052776, 0.00105552, 
               0.000211104, 0.000422208, 0.000422208, 0.000422208, 0.000105552, 
               0.000316656, 0.000105552, 0.000316656, 0.000211104, 0.000422208, 
               0.000211104, 0, 0.000316656, 0, 0.000211104, 
               0.000211104, 0.000422208, 0, 0.000211104, 0.000105552, 
               0.000316656, 0.000316656, 0.000211104, 0.000211104, 0.000211104, 
               0.000211104, 0.000211104, 0.00052776, 0.000211104, 0.000316656, 
               0.000105552, 0.000316656, 0, 0, 0, 
               0.000211104, 0.000316656, 0, 0, 0.000105552, 
               0.000105552, 0.000105552, 0.000105552, 0.000422208, 0.000105552, 
               0.000422208, 0.000105552, 0.000105552, 0.000105552, 0, 
               0, 0, 0.000422208, 0, 0, 
               0.000211104, 0, 0.000105552, 0.000316656, 0, 
               0, 0, 0, 0.000211104, 0.000316656, 
               0.000211104, 0, 0, 0, 0.000105552, 
               0.000105552, 0.000422208, 0, 0.000211104, 0.000105552, 
               0.000211104, 0, 0, 0, 0
            }, 
            {
               0, 0, 0, 0.000105686, 0.000105686, 
               0, 0, 0, 0.000105686, 0.000105686, 
               0.000105686, 0, 0.000317058, 0.000105686, 0.000105686, 
               0.000105686, 0.000422744, 0, 0.000105686, 0.000105686, 
               0.000105686, 0.000211372, 0.000211372, 0, 0.000105686, 
               0, 0.000422744, 0.000105686, 0.000211372, 0.00052843, 
               0.000951173, 0.00126823, 0.00158529, 0.0022194, 0.00264215, 
               0.00338195, 0.00401606, 0.00570704, 0.00581272, 0.00686958, 
               0.00951173, 0.0102515, 0.00887762, 0.0114141, 0.014796, 
               0.0126823, 0.0121539, 0.0155358, 0.0158529, 0.0177552, 
               0.0171211, 0.0150074, 0.0193405, 0.0175439, 0.020186, 
               0.0198689, 0.0190235, 0.0174382, 0.0179666, 0.0171211, 
               0.0158529, 0.0190235, 0.0189178, 0.018495, 0.0190235, 
               0.0174382, 0.0180723, 0.0176495, 0.0154301, 0.0183893, 
               0.0129994, 0.0131051, 0.0131051, 0.0132107, 0.0152188, 
               0.0142676, 0.0152188, 0.0140562, 0.0128937, 0.0119425, 
               0.0133164, 0.0138449, 0.0124709, 0.00919467, 0.0100402, 
               0.011097, 0.00940605, 0.0105686, 0.00877193, 0.00803213, 
               0.00845487, 0.0082435, 0.00718664, 0.0075037, 0.00729233, 
               0.00581272, 0.00634115, 0.00623547, 0.00739801, 0.00560135, 
               0.0060241, 0.00549567, 0.00517861, 0.0060241, 0.00412175, 
               0.00475587, 0.00443881, 0.00285352, 0.00433312, 0.00348763, 
               0.00401606, 0.00401606, 0.00285352, 0.00295921, 0.00317058, 
               0.00454449, 0.00359332, 0.00359332, 0.00253646, 0.00327626, 
               0.00359332, 0.00264215, 0.00285352, 0.00348763, 0.0014796, 
               0.00158529, 0.00211372, 0.00200803, 0.00179666, 0.00232509, 
               0.00253646, 0.00190235, 0.000739801, 0.0022194, 0.00158529, 
               0.00179666, 0.00200803, 0.00169097, 0.000739801, 0.00126823, 
               0.000845487, 0.000845487, 0.0014796, 0.00052843, 0.00105686, 
               0.000739801, 0.00116254, 0.000739801, 0.00126823, 0.000951173, 
               0.000634115, 0.000317058, 0.00052843, 0.000951173, 0.000951173, 
               0.000739801, 0.000739801, 0.000951173, 0.000845487, 0.000951173, 
               0.00052843, 0.000422744, 0.000739801, 0.000422744, 0.000739801, 
               0.000105686, 0.000317058, 0.00052843, 0.00052843, 0.00052843, 
               0.000845487, 0.000105686, 0.000317058, 0.000634115, 0.00052843, 
               0.000317058, 0.000634115, 0.000317058, 0.000211372, 0.000422744, 
               0.000211372, 0.000422744, 0.000105686, 0.000422744, 0.000105686, 
               0.000105686, 0.000105686, 0.000105686, 0.00052843, 0.000211372, 
               0.000317058, 0.00052843, 0.000105686, 0.000422744, 0.000211372, 
               0.000105686, 0, 0.000211372, 0.000211372, 0.000211372, 
               0.000634115, 0.000105686, 0, 0.000211372, 0.000317058, 
               0.000105686, 0.000211372, 0.000211372, 0.000211372, 0.000317058, 
               0, 0, 0, 0.00052843, 0.000105686, 
               0.000105686, 0.000105686, 0.000211372, 0.000105686, 0.000211372, 
               0, 0.000211372, 0.000105686, 0.000105686, 0.000211372, 
               0, 0.000105686, 0, 0, 0.000317058, 
               0, 0.000105686, 0.000211372, 0.000105686, 0, 
               0, 0, 0.000422744, 0, 0, 
               0, 0, 0.000105686, 0, 0, 
               0, 0.000105686, 0, 0, 0.000105686
            }
          };

  // Time Bin of Max. Cluster Distributions for e and pi
  fNTbins=20;
  fTBinSize=1;

  Int_t timBin[kNo_Mom][kNo_TBins]={  
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }, 
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }, 
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }, 
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }, 
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }, 
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }, 
            {
               0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
               10, 11, 12, 13, 14, 15, 16, 17, 18, 19
            }
          };

  Double_t probPIT[kNo_Mom][kNo_TBins]={  
            {
               0, 0.252938, 0.0801282, 0.0559562, 0.0532853 ,
               0.0491453, 0.0436699, 0.0431357, 0.0455395, 0.0422009 ,
               0.045406, 0.0432692, 0.0396635, 0.0419338, 0.0404647 ,
               0.0409989, 0.0389957, 0.0432692, 0, 0
            }, 
            {
               0, 0.265805, 0.076788, 0.057152, 0.0491699 ,
               0.041507, 0.0451788, 0.0506066, 0.0458174, 0.0439017 ,
               0.0431034, 0.0434227, 0.0431034, 0.0426245, 0.0391124 ,
               0.0408685, 0.0322478, 0.0395913, 0, 0
            }, 
            {
               0, 0.278986, 0.0757246, 0.0601449, 0.0496377 ,
               0.044686, 0.0422705, 0.0433575, 0.0388889, 0.0396135 ,
               0.042029, 0.0402174, 0.0427536, 0.0433575, 0.0376812 ,
               0.0386473, 0.0379227, 0.0440821, 0, 0
            }, 
            {
               0, 0.269746, 0.0749882, 0.0590033, 0.0465444 ,
               0.0470146, 0.0429008, 0.0429008, 0.0429008, 0.0440762 ,
               0.0394922, 0.0419605, 0.0446638, 0.0413728, 0.0400799 ,
               0.0393747, 0.0393747, 0.043606, 0, 0
            }, 
            {
               0, 0.269387, 0.0710154, 0.0571979, 0.0470223 ,
               0.0498072, 0.0501285, 0.0433805, 0.0460583, 0.0426307 ,
               0.043916, 0.0443445, 0.0435947, 0.0392031, 0.0344901 ,
               0.0396315, 0.0357755, 0.0424165, 0, 0
            }, 
            {
               0, 0.268377, 0.0723532, 0.0542071, 0.0520111 ,
               0.0447295, 0.0455386, 0.0420712, 0.0473879, 0.0453074 ,
               0.0420712, 0.0438049, 0.0420712, 0.0449607, 0.0394129 ,
               0.0372168, 0.0340962, 0.0443828, 0, 0
            }, 
            {
               0, 0.280106, 0.0703704, 0.0536508, 0.0493122 ,
               0.042328, 0.0434921, 0.0385185, 0.0421164, 0.0437037 ,
               0.0393651, 0.0404233, 0.0442328, 0.047619, 0.0390476 ,
               0.0386243, 0.0409524, 0.0461376, 0, 0
            }
          };

  Double_t probELT[kNo_Mom][kNo_TBins]={  
            {
               0, 0.190173, 0.0727839, 0.0512603, 0.0444633, 
               0.0446049, 0.0404984, 0.0346927, 0.0376664, 0.0379496, 
               0.0407816, 0.0397904, 0.0421977, 0.0468706, 0.049561, 
               0.0560748, 0.068819, 0.101813, 0, 0
            }, 
            {
               0, 0.195321, 0.0639535, 0.0449889, 0.0412514, 
               0.040144, 0.0390365, 0.0398671, 0.0366833, 0.0368217, 
               0.0336379, 0.0408361, 0.0460963, 0.053433, 0.0510797, 
               0.0617386, 0.0661683, 0.108942, 0, 0
            }, 
            {
               0, 0.175017, 0.0619995, 0.0464425, 0.039579, 
               0.0377488, 0.0386639, 0.0359186, 0.0412949, 0.0367193, 
               0.0385495, 0.0431251, 0.0419812, 0.0465569, 0.0553649, 
               0.0630291, 0.0720659, 0.125944, 0, 0
            }, 
            {
               0, 0.160783, 0.0604585, 0.0420978, 0.0420978, 
               0.035707, 0.0361128, 0.0348955, 0.0390546, 0.0412863, 
               0.0404747, 0.0438223, 0.0463583, 0.0521404, 0.0569081, 
               0.0619801, 0.0722256, 0.133597, 0, 0
            }, 
            {
               0, 0.159292, 0.054783, 0.0401391, 0.0371892, 
               0.0404551, 0.0385588, 0.0414033, 0.0363464, 0.0360303, 
               0.0441424, 0.0442478, 0.0510957, 0.0517278, 0.0552044, 
               0.0694269, 0.0711125, 0.128845, 0, 0
            }, 
            {
               0, 0.15495, 0.0529871, 0.040532, 0.0393709, 
               0.037471, 0.0407431, 0.0376821, 0.0385265, 0.0403209, 
               0.0401098, 0.0397931, 0.0525649, 0.0545704, 0.0575259, 
               0.0670255, 0.0754697, 0.130357, 0, 0
            }, 
            {
               0, 0.157366, 0.0504122, 0.0432255, 0.0394208, 
               0.0384697, 0.0374128, 0.0374128, 0.0394208, 0.0370958, 
               0.0430142, 0.0455506, 0.0467132, 0.0526316, 0.0583386, 
               0.0616149, 0.0783133, 0.133587, 0, 0
            }
          };

  for(Int_t ip=0; ip<fNMom; ip++) {
    fTrackMomentum[ip]= trackMomentum[ip];
    for(Int_t ie=0; ie<fNEbins; ie++) {
      fEnergyLoss[ip][ie]=energyLoss[ip][ie];
      fProbPI[ip][ie]=probPI[ip][ie]; 
      fProbEL[ip][ie]=probEL[ip][ie];
    }
  }

  for(Int_t ip=0; ip<fNMom; ip++) {
    fTrackMomentum[ip]=trackMomentum[ip] ;
    for(Int_t it=0; it<fNTbins; it++) {
      fTimBin[ip][it]=timBin[ip][it];
      fProbPIT[ip][it]=probPIT[ip][it];
      fProbELT[ip][it]=probELT[ip][it];
    }
  }

}


/*
void AliTRDprobdist::FillData()
{
  //
  // Energy loss Distributions for e and pi
  Double_t temp;
  Double_t fADCNorm=1.;
  ifstream filein("pp.dat");
  filein >> fNMom;
  for(Int_t ip=0; ip<fNMom; ip++) {
    filein >> fTrackMomentum[ip] >> fNEbins;
    for(Int_t ie=0; ie<fNEbins; ie++) {
      filein >> fEnergyLoss[ip][ie]>> fProbPI[ip][ie]>> fProbEL[ip][ie];
      fEnergyLoss[ip][ie]=fADCNorm*fEnergyLoss[ip][ie];
    }
  }
  fEnBinSize=fEnergyLoss[0][1]-fEnergyLoss[0][0];

  // Time Bin of Max. Cluster Distributions for e and pi
  ifstream fileint("tt.dat");
  fileint >> fNMom;
  for(Int_t ip=0; ip<fNMom; ip++) {
    fileint >> fTrackMomentum[ip] >> fNTbins;
    for(Int_t it=0; it<fNTbins; it++) {
      fileint >> fTimBin[ip][it]>> fProbPIT[ip][it]>> fProbELT[ip][it];
    }
  }
  fTBinSize=1;
}

*/