ClassImp(AliTRDprobdist)
//_________________________________________________________________________
-AliTRDprobdist::AliTRDprobdist(Int_t Multiplicity)
+AliTRDprobdist::AliTRDprobdist(Char_t *responseFile)
{
//
// The main constructor
//
- if (Multiplicity == 1) FillData();
- // if (Multiplicity == 2000) FillData2000();
- // if (Multiplicity == 4000) FillData4000();
- // if (Multiplicity == 6000) FillData6000();
- // if (Multiplicity == 8000) FillData8000();
+ Char_t *partName[5] = {"electron", "muon", "pion", "kaon", "proton"};
+ for(Int_t i=0; i<5; i++) {
+ fpartName[i] = partName[i];
+ }
+ // ADC Gain normalization
+ fADCNorm=1.0;
+ fNMom=kNMom;
+ // Set track momenta for which response functions are available
+ Double_t trackMomentum[kNMom]= {0.6, 0.8, 1, 1.5, 2, 3, 4, 5, 6, 8, 10};
+ for(Int_t imom=0; imom<fNMom; imom++) {
+ fTrackMomentum[imom] = trackMomentum[imom];
+ }
+ // Read the histogram file
+ ReadData(responseFile);
+ //
}
-//_________________________________________________________________________
-Double_t AliTRDprobdist::GetMeanPI(Int_t ip)
+Bool_t AliTRDprobdist::ReadData(Char_t *responseFile)
{
- // 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];
+ //
+ // Read the TRD dEdx histograms.
+ //
+ // Read histogram Root file
+ TFile *histFile = new TFile(responseFile);
+ if (!histFile || !histFile->IsOpen()) {
+ Error("AliTRDprobdist", "opening TRD histgram file %s failed", responseFile);
+ return kFALSE;
}
- return integral/norm;
+
+ // Read histograms
+ Char_t text[200];
+ for (Int_t imom = 0; imom < kNMom; imom++) {
+ sprintf(text,"h1dEdxEL%01d",imom+1);
+ fh1dEdxEL[imom] = (TH1F*)histFile->Get(text);
+ fh1dEdxEL[imom]->Scale(1.0/fh1dEdxEL[imom]->Integral());
+
+ sprintf(text,"h1dEdxPI%01d",imom+1);
+ fh1dEdxPI[imom] = (TH1F*)histFile->Get(text);
+ fh1dEdxPI[imom]->Scale(1.0/fh1dEdxPI[imom]->Integral());
+
+ sprintf(text,"h1dEdxMU%01d",imom+1);
+ fh1dEdxMU[imom] = (TH1F*)histFile->Get(text);
+ fh1dEdxMU[imom]->Scale(1.0/fh1dEdxMU[imom]->Integral());
+
+ sprintf(text,"h1dEdxKA%01d",imom+1);
+ fh1dEdxKA[imom] = (TH1F*)histFile->Get(text);
+ fh1dEdxKA[imom]->Scale(1.0/fh1dEdxKA[imom]->Integral());
+
+ sprintf(text,"h1dEdxPR%01d",imom+1);
+ fh1dEdxPR[imom] = (TH1F*)histFile->Get(text);
+ fh1dEdxPR[imom]->Scale(1.0/fh1dEdxPR[imom]->Integral());
+
+ sprintf(text,"h1MaxTimBinEL%01d",imom+1);
+ fh1MaxTimBinEL[imom] = (TH1F*)histFile->Get(text);
+ fh1MaxTimBinEL[imom]->Scale(1.0/fh1MaxTimBinEL[imom]->Integral());
+
+ sprintf(text,"h1MaxTimBinPI%01d",imom+1);
+ fh1MaxTimBinPI[imom] = (TH1F*)histFile->Get(text);
+ fh1MaxTimBinPI[imom]->Scale(1.0/fh1MaxTimBinPI[imom]->Integral());
+ }
+ // Number of bins and bin size
+ fNbins = fh1dEdxPI[1]->GetNbinsX();
+ fBinSize = fh1dEdxPI[1]->GetBinWidth(1);
+ return kTRUE;
}
-//_________________________________________________________________________
-Double_t AliTRDprobdist::GetMeanEL(Int_t ip)
+AliTRDprobdist::AliTRDprobdist(const AliTRDprobdist& pd):TNamed()
{
//
- // 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];
+ // Copy constructor.
+ //
+ for(Int_t i=0; i<5; i++) {
+ fpartName[i] = pd.fpartName[i];
+ }
+ // ADC Gain normalization
+ fADCNorm=pd.fADCNorm;
+ fNMom=pd.fNMom;
+ // Set track momenta for which response functions are available
+ for(Int_t imom=0; imom<fNMom; imom++) {
+ fTrackMomentum[imom] = pd.fTrackMomentum[imom];
}
- return integral/norm;
+
+ fNbins=pd.fNbins;
+ fBinSize=pd.fBinSize;
+
+ for (Int_t imom = 0; imom < kNMom; imom++) {
+ fh1dEdxEL[imom] = pd.fh1dEdxEL[imom];
+ fh1dEdxPI[imom] = pd.fh1dEdxPI[imom];
+ fh1dEdxMU[imom] = pd.fh1dEdxMU[imom];
+ fh1dEdxKA[imom] = pd.fh1dEdxKA[imom];
+ fh1dEdxPR[imom] = pd.fh1dEdxPR[imom];
+ fh1MaxTimBinEL[imom] = pd.fh1MaxTimBinEL[imom];
+ fh1MaxTimBinPI[imom] = pd.fh1MaxTimBinPI[imom];
+ }
+}
+
+
+AliTRDprobdist::~AliTRDprobdist()
+{
+ // Destructor
+
}
+
//_________________________________________________________________________
-Double_t AliTRDprobdist::GetNormalizationPI(Int_t ip)
+Double_t AliTRDprobdist::GetMean(Int_t k, 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;
+ // Gets mean of de/dx dist. of e
+ printf("Mean for particle = %s and momentum = %.2f is:\n", fpartName[k], fTrackMomentum[ip]);
+ if(k==0) return fh1dEdxEL[ip]->GetMean();
+ if(k==1) return fh1dEdxMU[ip]->GetMean();
+ if(k==2) return fh1dEdxPI[ip]->GetMean();
+ if(k==3) return fh1dEdxKA[ip]->GetMean();
+ if(k==4) return fh1dEdxPR[ip]->GetMean();
+ return fh1dEdxPR[ip]->GetMean();
}
//_________________________________________________________________________
-Double_t AliTRDprobdist::GetNormalizationEL(Int_t ip)
+Double_t AliTRDprobdist::GetNormalization(Int_t k, 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;
+
+ printf("Normalization for particle = %s and momentum = %.2f is:\n",fpartName[k], fTrackMomentum[ip]);
+ if(k==0) return fh1dEdxEL[ip]->Integral();
+ if(k==1) return fh1dEdxMU[ip]->Integral();
+ if(k==2) return fh1dEdxPI[ip]->Integral();
+ if(k==3) return fh1dEdxKA[ip]->Integral();
+ if(k==4) return fh1dEdxPR[ip]->Integral();
+ return fh1dEdxPR[ip]->Integral();
+}
+
+TH1F* AliTRDprobdist::GetHistogram(Int_t k, Int_t ip) const
+{
+ //
+ //
+ printf("Histogram for particle = %s and momentum = %.2f is:\n", fpartName[k], fTrackMomentum[ip]);
+ if(k==0) return fh1dEdxEL[ip];
+ if(k==1) return fh1dEdxMU[ip];
+ if(k==2) return fh1dEdxPI[ip];
+ if(k==3) return fh1dEdxKA[ip];
+ if(k==4) return fh1dEdxPR[ip];
+ return fh1dEdxPR[ip];
}
//_________________________________________________________________________
-Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx)
+Double_t AliTRDprobdist::GetProbability(Int_t k, Double_t mom, Double_t dedx1) 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;
+ Double_t probability = 1.0;
+ Double_t dedx = dedx1/fADCNorm;
+ Int_t iEnBin= ((Int_t) (dedx/fBinSize+1));
+ if(iEnBin > fNbins) iEnBin = fNbins;
+
+ ////Electron//////////////////////////
+ if(k==0){ // electron
+ Double_t slop;
+ // Lower limit
+ if(mom<=fTrackMomentum[0]) {
+ slop=(fh1dEdxEL[1]->GetBinContent(iEnBin)-fh1dEdxEL[0]->GetBinContent(iEnBin))/(fTrackMomentum[1] - fTrackMomentum[0]);
+ probability= fh1dEdxEL[0]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[0]);
+ return probability;
+ }
+ // Upper Limit
+ if(mom>=fTrackMomentum[fNMom-1]) {
+ slop=(fh1dEdxEL[fNMom-1]->GetBinContent(iEnBin)-fh1dEdxEL[fNMom-2]->GetBinContent(iEnBin))/(fTrackMomentum[fNMom-1] - fTrackMomentum[fNMom-2]);
+ probability= fh1dEdxEL[fNMom-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[fNMom-1]);
+ return probability;
+ }
+ // In the range
+ for(Int_t ip=1; ip<fNMom; ip++)
+ if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
+ slop=(fh1dEdxEL[ip]->GetBinContent(iEnBin)-fh1dEdxEL[ip-1]->GetBinContent(iEnBin))/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
+ // Linear Interpolation
+ probability= fh1dEdxEL[ip-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[ip-1]);
+ return probability;
+ }
+ }
- if(k==0){ // electron
- if(mom<=fTrackMomentum[0]) Probability = fProbEL[0][iEnBin];
- if(mom>=fTrackMomentum[fNMom-1]) Probability = fProbEL[fNMom-1][iEnBin];
+ ////Pion//////////////////////////
+ if(k==2){ // Pion
+ Double_t slop;
+ // Lower limit
+ if(mom<=fTrackMomentum[0]) {
+ slop=(fh1dEdxPI[1]->GetBinContent(iEnBin)-fh1dEdxPI[0]->GetBinContent(iEnBin))/(fTrackMomentum[1] - fTrackMomentum[0]);
+ probability= fh1dEdxPI[0]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[0]);
+ return probability;
+ }
+ // Upper Limit
+ if(mom>=fTrackMomentum[fNMom-1]) {
+ slop=(fh1dEdxPI[fNMom-1]->GetBinContent(iEnBin)-fh1dEdxPI[fNMom-2]->GetBinContent(iEnBin))/(fTrackMomentum[fNMom-1] - fTrackMomentum[fNMom-2]);
+ probability= fh1dEdxPI[fNMom-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[fNMom-1]);
+ return probability;
+ }
+ // In the range
+ for(Int_t ip=1; ip<fNMom; ip++)
+ if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
+ slop=(fh1dEdxPI[ip]->GetBinContent(iEnBin)-fh1dEdxPI[ip-1]->GetBinContent(iEnBin))/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
+ // Linear Interpolation
+ probability= fh1dEdxPI[ip-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[ip-1]);
+ return probability;
+ }
}
- if(k==2){ // pion
- if(mom<=fTrackMomentum[0]) Probability = fProbPI[0][iEnBin];
- if(mom>=fTrackMomentum[fNMom-1]) Probability = fProbPI[fNMom-1][iEnBin];
+
+ ////Muon//////////////////////////
+ if(k==1){ // Muon
+ Double_t slop;
+ // Lower limit
+ if(mom<=fTrackMomentum[0]) {
+ slop=(fh1dEdxMU[1]->GetBinContent(iEnBin)-fh1dEdxMU[0]->GetBinContent(iEnBin))/(fTrackMomentum[1] - fTrackMomentum[0]);
+ probability= fh1dEdxMU[0]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[0]);
+ return probability;
+ }
+ // Upper Limit
+ if(mom>=fTrackMomentum[fNMom-1]) {
+ slop=(fh1dEdxMU[fNMom-1]->GetBinContent(iEnBin)-fh1dEdxMU[fNMom-2]->GetBinContent(iEnBin))/(fTrackMomentum[fNMom-1] - fTrackMomentum[fNMom-2]);
+ probability= fh1dEdxMU[fNMom-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[fNMom-1]);
+ return probability;
+ }
+ // In the range
+ for(Int_t ip=1; ip<fNMom; ip++)
+ if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
+ slop=(fh1dEdxMU[ip]->GetBinContent(iEnBin)-fh1dEdxMU[ip-1]->GetBinContent(iEnBin))/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
+ // Linear Interpolation
+ probability= fh1dEdxMU[ip-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[ip-1]);
+ return probability;
+ }
}
- 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;
+ ////Kaon//////////////////////////
+ if(k==3){ // Kaon
+ Double_t slop;
+ // Lower limit
+ if(mom<=fTrackMomentum[0]) {
+ slop=(fh1dEdxKA[1]->GetBinContent(iEnBin)-fh1dEdxKA[0]->GetBinContent(iEnBin))/(fTrackMomentum[1] - fTrackMomentum[0]);
+ probability= fh1dEdxKA[0]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[0]);
+ return probability;
+ }
+ // Upper Limit
+ if(mom>=fTrackMomentum[fNMom-1]) {
+ slop=(fh1dEdxKA[fNMom-1]->GetBinContent(iEnBin)-fh1dEdxKA[fNMom-2]->GetBinContent(iEnBin))/(fTrackMomentum[fNMom-1] - fTrackMomentum[fNMom-2]);
+ probability= fh1dEdxKA[fNMom-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[fNMom-1]);
+ return probability;
+ }
+ // In the range
+ for(Int_t ip=1; ip<fNMom; ip++)
+ if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
+ slop=(fh1dEdxKA[ip]->GetBinContent(iEnBin)-fh1dEdxKA[ip-1]->GetBinContent(iEnBin))/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
+ // Linear Interpolation
+ probability= fh1dEdxKA[ip-1]->GetBinContent(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;
+ ////Proton//////////////////////////
+ if(k==4){ // Proton
+ Double_t slop;
+ // Lower limit
+ if(mom<=fTrackMomentum[0]) {
+ slop=(fh1dEdxPR[1]->GetBinContent(iEnBin)-fh1dEdxPR[0]->GetBinContent(iEnBin))/(fTrackMomentum[1] - fTrackMomentum[0]);
+ probability= fh1dEdxPR[0]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[0]);
+ return probability;
+ }
+ // Upper Limit
+ if(mom>=fTrackMomentum[fNMom-1]) {
+ slop=(fh1dEdxPR[fNMom-1]->GetBinContent(iEnBin)-fh1dEdxPR[fNMom-2]->GetBinContent(iEnBin))/(fTrackMomentum[fNMom-1] - fTrackMomentum[fNMom-2]);
+ probability= fh1dEdxPR[fNMom-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[fNMom-1]);
+ return probability;
+ }
+ // In the range
+ for(Int_t ip=1; ip<fNMom; ip++)
+ if((fTrackMomentum[ip-1]<= mom) && (mom<fTrackMomentum[ip])) {
+ slop=(fh1dEdxPR[ip]->GetBinContent(iEnBin)-fh1dEdxPR[ip-1]->GetBinContent(iEnBin))/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
+ // Linear Interpolation
+ probability= fh1dEdxPR[ip-1]->GetBinContent(iEnBin) + slop*(mom-fTrackMomentum[ip-1]);
+ return probability;
+ }
}
- return Probability;
+
+ return probability;
}
//_________________________________________________________________________
-Double_t AliTRDprobdist::GetProbabilityT(Int_t k, Double_t mom, Int_t timbin)
+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;
+ Double_t probabilityT = 1.0;
if(timbin<=0) return 0.;
- Int_t iTBin=timbin;
+ Int_t iTBin=timbin+1;
if(k==0){ // electron
- if(mom<=fTrackMomentum[0]) ProbabilityT = fProbELT[0][iTBin];
- if(mom>=fTrackMomentum[fNMom-1]) ProbabilityT = fProbELT[fNMom-1][iTBin];
+ if(mom<=fTrackMomentum[0]) probabilityT = fh1MaxTimBinEL[0]->GetBinContent(iTBin);
+ if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fh1MaxTimBinEL[fNMom-1]->GetBinContent(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==1||k==2||k==3||k==4){ // pion
+ if(mom<=fTrackMomentum[0]) probabilityT = fh1MaxTimBinPI[0]->GetBinContent(iTBin);
+ if(mom>=fTrackMomentum[fNMom-1]) probabilityT = fh1MaxTimBinPI[fNMom-1]->GetBinContent(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]);
+ Double_t slop=(fh1MaxTimBinEL[ip]->GetBinContent(iTBin)-fh1MaxTimBinEL[ip-1]->GetBinContent(iTBin))/(fTrackMomentum[ip] - fTrackMomentum[ip-1]);
// Linear Interpolation
- ProbabilityT= fProbELT[ip-1][iTBin] + slop*(mom-fTrackMomentum[ip-1]);
- return ProbabilityT;
+ probabilityT= fh1MaxTimBinEL[ip-1]->GetBinContent(iTBin) + slop*(mom-fTrackMomentum[ip-1]);
+ return probabilityT;
}
- if(k==2) // pion
+ if(k==1||k==2||k==3||k==4) // pion and other particles
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]);
+ Double_t slop=(fh1MaxTimBinPI[ip]->GetBinContent(iTBin)-fh1MaxTimBinPI[ip-1]->GetBinContent(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=200;
- 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
- },
- {
- 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
- },
- {
- 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
- },
- {
- 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
- },
- {
- 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
- },
- {
- 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
- },
- {
- 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
- }
- };
-
- Double_t ProbPI[kNo_Mom][kNo_EnBins]={
- {
- 0, 0, 0, 0, 0,
- 0, 0.000136166, 0.000408497, 0, 0,
- 0.000408497, 0.000816993, 0.00231481, 0.00381264, 0.00558279,
- 0.0111656, 0.0193355, 0.0270969, 0.0299564, 0.0363562,
- 0.0364924, 0.0472495, 0.0473856, 0.0430283, 0.045207,
- 0.0439815, 0.0401688, 0.03622, 0.0336329, 0.0326797,
- 0.0324074, 0.0281863, 0.0257353, 0.0228758, 0.0231481,
- 0.0204248, 0.0219227, 0.0206972, 0.0174292, 0.0142974,
- 0.0132081, 0.0108932, 0.00994009, 0.0113017, 0.0106209,
- 0.00680828, 0.00803377, 0.00708061, 0.00885076, 0.00544662,
- 0.00694444, 0.00667211, 0.00408497, 0.00653595, 0.00503813,
- 0.00558279, 0.00367647, 0.00367647, 0.0039488, 0.00258715,
- 0.00367647, 0.00217865, 0.00340414, 0.00299564, 0.00245098,
- 0.00285948, 0.00204248, 0.00285948, 0.00190632, 0.00272331,
- 0.00204248, 0.000953159, 0.00136166, 0.00231481, 0.00108932,
- 0.00204248, 0.00245098, 0.00136166, 0.00177015, 0.00149782,
- 0.00108932, 0.000953159, 0.000544662, 0.00122549, 0.00108932,
- 0.00149782, 0.00108932, 0.000816993, 0.000953159, 0.000136166,
- 0.00108932, 0.000816993, 0.00122549, 0.000816993, 0.000272331,
- 0.000680828, 0.000816993, 0.000544662, 0, 0.000544662,
- 0.000544662, 0.000544662, 0.000680828, 0.000680828, 0.000272331,
- 0.000136166, 0.000136166, 0.000544662, 0.000272331, 0.000408497,
- 0.000816993, 0.000136166, 0.000953159, 0.000136166, 0.000272331,
- 0.000272331, 0.000953159, 0.000272331, 0.000544662, 0,
- 0.000408497, 0.000544662, 0.000272331, 0.000136166, 0.000272331,
- 0.000136166, 0.000408497, 0.000136166, 0.000408497, 0.000272331,
- 0.000272331, 0.000408497, 0, 0.000136166, 0.000136166,
- 0.000272331, 0.000136166, 0.000136166, 0.000136166, 0.000136166,
- 0.000136166, 0, 0.000408497, 0.000408497, 0,
- 0.000136166, 0.000136166, 0.000680828, 0.000408497, 0.000272331,
- 0.000136166, 0.000136166, 0.000136166, 0.000136166, 0.000272331,
- 0.000136166, 0.000272331, 0.000408497, 0.000272331, 0.000136166,
- 0.000136166, 0.000408497, 0.000272331, 0.000272331, 0.000136166,
- 0, 0.000136166, 0, 0, 0,
- 0, 0.000408497, 0.000272331, 0.000272331, 0.000136166,
- 0.000136166, 0.000272331, 0, 0.000136166, 0.000136166,
- 0.000136166, 0.000136166, 0, 0, 0,
- 0.000272331, 0, 0.000272331, 0, 0.000136166,
- 0, 0, 0, 0, 0,
- 0.000136166, 0, 0.000272331, 0.000136166, 0.000136166
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0, 0,
- 0.000124008, 0.000372024, 0.00124008, 0.00198413, 0.00285218,
- 0.00806052, 0.0122768, 0.0205853, 0.0249256, 0.03063,
- 0.0375744, 0.0390625, 0.0404266, 0.0460069, 0.0455109,
- 0.0438988, 0.0414187, 0.0403026, 0.0400546, 0.0342262,
- 0.0345982, 0.0303819, 0.0300099, 0.0288938, 0.0257937,
- 0.0231895, 0.0183532, 0.0194692, 0.0182292, 0.0141369,
- 0.0138889, 0.0110367, 0.0122768, 0.0120288, 0.0102927,
- 0.00954861, 0.00892857, 0.00756448, 0.0078125, 0.00682044,
- 0.00607639, 0.00818452, 0.00545635, 0.00458829, 0.00558036,
- 0.00496032, 0.00545635, 0.00421627, 0.00359623, 0.00272817,
- 0.00372024, 0.00210813, 0.00272817, 0.00409226, 0.00334821,
- 0.00223214, 0.00347222, 0.00173611, 0.00136409, 0.00235615,
- 0.00260417, 0.00248016, 0.00235615, 0.00210813, 0.0014881,
- 0.00186012, 0.0014881, 0.00124008, 0.00124008, 0.0014881,
- 0.00111607, 0.00124008, 0.00111607, 0.000992063, 0.00062004,
- 0.00136409, 0.000744048, 0.00136409, 0.000868056, 0.00124008,
- 0.000744048, 0.000992063, 0.000744048, 0.000992063, 0.00111607,
- 0.000496032, 0.00136409, 0.000496032, 0.000372024, 0.000868056,
- 0.000744048, 0.000124008, 0.000496032, 0.00062004, 0.000124008,
- 0.00062004, 0.000372024, 0.000992063, 0.000372024, 0.000372024,
- 0.000372024, 0.000372024, 0.000744048, 0.000496032, 0.000372024,
- 0.000496032, 0.00062004, 0.000744048, 0.000496032, 0.000868056,
- 0.000496032, 0.000372024, 0.000124008, 0.00062004, 0.000372024,
- 0.000248016, 0, 0.000248016, 0.000124008, 0.000496032,
- 0.000496032, 0.000248016, 0.000372024, 0, 0,
- 0.000248016, 0.000124008, 0.000248016, 0.000248016, 0.000124008,
- 0.000124008, 0.000124008, 0.000372024, 0, 0.000124008,
- 0.000248016, 0, 0.000124008, 0.000124008, 0.000124008,
- 0.000124008, 0.000744048, 0.000248016, 0, 0.00062004,
- 0.000124008, 0.000372024, 0, 0.000124008, 0.000124008,
- 0, 0, 0.000248016, 0, 0.000248016,
- 0.000248016, 0.000124008, 0.000248016, 0.000124008, 0,
- 0.000248016, 0, 0.000248016, 0, 0,
- 0, 0.000124008, 0.000248016, 0.000124008, 0.000124008,
- 0.000124008, 0.000496032, 0.000248016, 0, 0.000124008,
- 0.000248016, 0.000124008, 0, 0.000124008, 0.000124008,
- 0, 0.000124008, 0, 0.000124008, 0.000496032,
- 0.000248016, 0, 0, 0.000124008, 0
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0, 0,
- 0, 0.000230362, 0.000460723, 0.00138217, 0.00334024,
- 0.00529832, 0.00817784, 0.014628, 0.0209629, 0.0222299,
- 0.0313292, 0.0338632, 0.039507, 0.0393918, 0.0412347,
- 0.0416955, 0.0423865, 0.0392767, 0.0403133, 0.0373186,
- 0.0336328, 0.0342087, 0.0277586, 0.0282193, 0.0244183,
- 0.0252246, 0.0214236, 0.0232665, 0.0179682, 0.0142824,
- 0.0131306, 0.0143976, 0.0137065, 0.0117484, 0.0126699,
- 0.0101359, 0.0101359, 0.00944483, 0.00956001, 0.0103663,
- 0.0084082, 0.00621977, 0.0059894, 0.00610458, 0.00679567,
- 0.00691085, 0.00472241, 0.00518314, 0.00460723, 0.00472241,
- 0.00518314, 0.00403133, 0.00380097, 0.00403133, 0.00391615,
- 0.00310988, 0.00264916, 0.00253398, 0.00276434, 0.00310988,
- 0.00218844, 0.00253398, 0.00322506, 0.00287952, 0.00207326,
- 0.00218844, 0.00218844, 0.000921447, 0.00172771, 0.00161253,
- 0.00138217, 0.00149735, 0.00172771, 0.00138217, 0.00126699,
- 0.00149735, 0.00149735, 0.00161253, 0.00103663, 0.00184289,
- 0.00103663, 0.00103663, 0.000575904, 0.000575904, 0.00115181,
- 0.000460723, 0.000806266, 0.00161253, 0.000921447, 0.00103663,
- 0.000575904, 0.00138217, 0.000806266, 0.000230362, 0.000345543,
- 0.000806266, 0.000460723, 0.000691085, 0.000921447, 0.000115181,
- 0.000230362, 0.000575904, 0.000345543, 0.000345543, 0.000345543,
- 0.000575904, 0.000345543, 0.000460723, 0.000230362, 0.000345543,
- 0.000691085, 0.000460723, 0.000115181, 0.000345543, 0.000345543,
- 0.000230362, 0.000115181, 0.000345543, 0.000230362, 0.000230362,
- 0.000691085, 0.000230362, 0.000345543, 0.000345543, 0.000230362,
- 0.000115181, 0.000230362, 0.000230362, 0, 0.000230362,
- 0, 0.000230362, 0.000345543, 0.000230362, 0.000575904,
- 0, 0.000345543, 0.000345543, 0.000230362, 0.000230362,
- 0.000115181, 0.000115181, 0.000345543, 0.000115181, 0.000115181,
- 0, 0.000115181, 0.000345543, 0.000575904, 0.000115181,
- 0, 0.000345543, 0.000115181, 0.000230362, 0,
- 0.000230362, 0.000115181, 0.000115181, 0, 0.000230362,
- 0.000230362, 0.000230362, 0.000230362, 0.000230362, 0,
- 0.000230362, 0, 0, 0, 0.000230362,
- 0, 0.000230362, 0, 0.000115181, 0,
- 0.000345543, 0.000115181, 0.000115181, 0, 0,
- 0.000115181, 0.000230362, 0, 0.000230362, 0,
- 0.000115181, 0.000115181, 0, 0, 0
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0.000108436, 0, 0.000108436, 0,
- 0, 0, 0.000108436, 0.000433745, 0.00173498,
- 0.00260247, 0.00466276, 0.0071568, 0.0107352, 0.0197354,
- 0.0228801, 0.0301453, 0.0318803, 0.0352418, 0.0386033,
- 0.040013, 0.0416396, 0.0367599, 0.0402299, 0.0412058,
- 0.0335068, 0.0358924, 0.0324225, 0.0309044, 0.0290609,
- 0.0249404, 0.0297116, 0.0242897, 0.0220126, 0.021362,
- 0.0182173, 0.0174583, 0.0147473, 0.0143136, 0.0122533,
- 0.00997614, 0.00943396, 0.00921709, 0.00997614, 0.00910865,
- 0.00824116, 0.00726524, 0.00661462, 0.00704836, 0.00628931,
- 0.00737367, 0.00672305, 0.00520494, 0.00628931, 0.00455433,
- 0.00487964, 0.00509651, 0.00314465, 0.00433745, 0.00422902,
- 0.00314465, 0.00336153, 0.00325309, 0.00412058, 0.00303622,
- 0.00227716, 0.00249404, 0.00184342, 0.00292778, 0.00249404,
- 0.00206029, 0.00227716, 0.00206029, 0.00206029, 0.00184342,
- 0.00216873, 0.00151811, 0.00173498, 0.00140967, 0.00130124,
- 0.00140967, 0.00108436, 0.0011928, 0.00108436, 0.000759054,
- 0.000650618, 0.00108436, 0.00184342, 0.00108436, 0.0011928,
- 0.000867491, 0.00108436, 0.000867491, 0.000433745, 0.000216873,
- 0.00108436, 0.000650618, 0.0011928, 0.000759054, 0.000650618,
- 0.000542182, 0.000433745, 0.00140967, 0.000542182, 0.000433745,
- 0.000650618, 0.000650618, 0.000650618, 0.000542182, 0.000650618,
- 0.000433745, 0.000216873, 0.000542182, 0.000542182, 0.000759054,
- 0.000650618, 0.000216873, 0.000433745, 0.000542182, 0.000216873,
- 0.000433745, 0.000759054, 0.000542182, 0.000325309, 0.000433745,
- 0.000216873, 0.000433745, 0.000325309, 0.000325309, 0.000108436,
- 0.000216873, 0.000216873, 0.000325309, 0.000108436, 0,
- 0.000433745, 0.000325309, 0.000216873, 0.000325309, 0,
- 0.000216873, 0.000325309, 0.000216873, 0, 0.000216873,
- 0.000108436, 0.000542182, 0.000433745, 0.000433745, 0.000216873,
- 0, 0.000216873, 0.000216873, 0.000108436, 0,
- 0.000216873, 0.000108436, 0.000433745, 0.000216873, 0.000216873,
- 0.000216873, 0, 0, 0.000108436, 0,
- 0.000108436, 0, 0, 0.000216873, 0.000108436,
- 0, 0.000108436, 0.000108436, 0.000108436, 0,
- 0, 0.000108436, 0.000216873, 0, 0.000108436,
- 0.000108436, 0, 0, 0.000325309, 0.000216873,
- 0.000108436, 0.000108436, 0, 0, 0.000216873,
- 0, 0.000108436, 0.000108436, 0.000108436, 0.000216873
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0, 0.000113611, 0, 0,
- 0, 0, 0.000113611, 0.000340832, 0.000454442,
- 0.00113611, 0.00363554, 0.0074983, 0.0096569, 0.0122699,
- 0.0222677, 0.0246535, 0.0296524, 0.034424, 0.0382868,
- 0.0378323, 0.0357873, 0.039082, 0.0399909, 0.0356737,
- 0.0388548, 0.0347648, 0.032379, 0.0328334, 0.0302204,
- 0.025676, 0.024199, 0.0220404, 0.0224949, 0.018973,
- 0.0212452, 0.0179505, 0.0140877, 0.0118155, 0.0145421,
- 0.013406, 0.0119291, 0.0117019, 0.0123835, 0.0102249,
- 0.00988412, 0.00715746, 0.0101113, 0.00852079, 0.00636219,
- 0.00670302, 0.00613497, 0.00624858, 0.00522608, 0.00556692,
- 0.00511247, 0.00352193, 0.00374915, 0.00352193, 0.00408998,
- 0.0053397, 0.00465803, 0.00443081, 0.00443081, 0.00227221,
- 0.00284026, 0.0043172, 0.00147694, 0.00227221, 0.00204499,
- 0.0031811, 0.00295387, 0.00249943, 0.00249943, 0.00170416,
- 0.00159055, 0.00193138, 0.00147694, 0.00238582, 0.00113611,
- 0.00147694, 0.000454442, 0.00181777, 0.00136333, 0.00204499,
- 0.00147694, 0.00136333, 0.00102249, 0.000795274, 0.00181777,
- 0.000908884, 0.000568053, 0.00147694, 0.000908884, 0.000454442,
- 0.000681663, 0.000681663, 0.000795274, 0.000568053, 0.000454442,
- 0.000908884, 0.00102249, 0.000454442, 0.00113611, 0.00102249,
- 0.000454442, 0.00147694, 0.000681663, 0.000681663, 0.000340832,
- 0.000568053, 0.000568053, 0.000908884, 0.000454442, 0.000795274,
- 0.000340832, 0.000454442, 0.000113611, 0.000113611, 0.000340832,
- 0.000568053, 0.000568053, 0.000340832, 0.000454442, 0.000454442,
- 0.000113611, 0.000795274, 0.000227221, 0.000340832, 0.000568053,
- 0.000340832, 0.000568053, 0.000340832, 0.000113611, 0.000227221,
- 0.000568053, 0.000340832, 0.000454442, 0, 0.000568053,
- 0.000795274, 0.000113611, 0, 0.000340832, 0.000340832,
- 0, 0.000454442, 0.000113611, 0.000113611, 0.000113611,
- 0.000340832, 0.000568053, 0.000113611, 0.000113611, 0,
- 0, 0.000227221, 0.000340832, 0.000113611, 0.000227221,
- 0.000227221, 0.000113611, 0, 0, 0,
- 0.000340832, 0.000113611, 0.000227221, 0.000227221, 0.000113611,
- 0.000340832, 0, 0.000340832, 0, 0,
- 0, 0.000113611, 0.000113611, 0.000113611, 0.000113611,
- 0.000454442, 0.000227221, 0.000227221, 0.000113611, 0.000113611,
- 0.000454442, 0.000113611, 0.000113611, 0.000113611, 0,
- 0, 0.000227221, 0, 0.000113611, 0
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0.000225836, 0.000112918, 0, 0.000112918,
- 0.000112918, 0.000225836, 0.000112918, 0.000338753, 0.000790425,
- 0.000677507, 0.00214544, 0.00361337, 0.00880759, 0.0112918,
- 0.0173893, 0.0203252, 0.0238257, 0.0310524, 0.0324074,
- 0.0368112, 0.0362466, 0.0394083, 0.0389566, 0.0423442,
- 0.0386179, 0.0363595, 0.0336495, 0.035682, 0.0290199,
- 0.0300361, 0.0265357, 0.0245032, 0.0257453, 0.0235998,
- 0.0208898, 0.0204381, 0.0160343, 0.0147922, 0.0147922,
- 0.0135501, 0.0120822, 0.0111789, 0.00982385, 0.0102755,
- 0.00677507, 0.00937218, 0.00835592, 0.00779133, 0.00609756,
- 0.00745257, 0.00722674, 0.00575881, 0.00564589, 0.00654923,
- 0.00474255, 0.0050813, 0.00417796, 0.00440379, 0.00429088,
- 0.00372629, 0.00496838, 0.00293586, 0.00293586, 0.00383921,
- 0.00361337, 0.00237127, 0.00293586, 0.0019196, 0.00169377,
- 0.00203252, 0.00225836, 0.00169377, 0.00372629, 0.00237127,
- 0.00101626, 0.00237127, 0.00169377, 0.00180668, 0.0012421,
- 0.00146793, 0.0012421, 0.00180668, 0.000790425, 0.00214544,
- 0.00135501, 0.00158085, 0.00135501, 0.000564589, 0.0012421,
- 0.00135501, 0.00101626, 0.00112918, 0.000677507, 0.000677507,
- 0.000564589, 0.00112918, 0.000564589, 0.000564589, 0.000677507,
- 0.000451671, 0.000790425, 0.000790425, 0.000677507, 0.000225836,
- 0.000903342, 0.000451671, 0.000790425, 0.00101626, 0.000564589,
- 0.000564589, 0.000564589, 0.000790425, 0.00101626, 0.000790425,
- 0.000564589, 0.000677507, 0.000564589, 0.000451671, 0.000677507,
- 0.000451671, 0.000451671, 0.000677507, 0.000564589, 0.000451671,
- 0.000451671, 0.000338753, 0.000225836, 0.000338753, 0.000451671,
- 0.000451671, 0.000338753, 0.000451671, 0.000112918, 0,
- 0.000338753, 0, 0.000225836, 0.000225836, 0.000112918,
- 0.000112918, 0.000225836, 0.000225836, 0.000225836, 0.000225836,
- 0.000112918, 0.000338753, 0.000112918, 0.000338753, 0.000564589,
- 0.000225836, 0.000338753, 0.000225836, 0.000338753, 0.000225836,
- 0.000225836, 0.000451671, 0, 0.000112918, 0,
- 0.000225836, 0.000112918, 0.000338753, 0.000112918, 0.000112918,
- 0.000225836, 0.000338753, 0.000225836, 0.000112918, 0.000112918,
- 0, 0.000338753, 0, 0.000225836, 0.000564589,
- 0, 0.000338753, 0, 0, 0,
- 0.000112918, 0.000112918, 0, 0, 0.000225836,
- 0.000112918, 0, 0.000112918, 0.000225836, 0,
- 0, 0, 0.000112918, 0.000112918, 0.000112918
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0.000106701, 0, 0.000106701, 0,
- 0, 0, 0.000106701, 0.000213402, 0.000320102,
- 0.000213402, 0.00170721, 0.00362783, 0.00416133, 0.00960307,
- 0.0128041, 0.0189927, 0.0253948, 0.028169, 0.0325437,
- 0.0368118, 0.039586, 0.0339309, 0.0373453, 0.0352113,
- 0.0364917, 0.0332907, 0.0308365, 0.0338242, 0.0290226,
- 0.0290226, 0.0256082, 0.028169, 0.022834, 0.0237943,
- 0.0224072, 0.0226206, 0.0181391, 0.0161118, 0.0157917,
- 0.0137644, 0.0132309, 0.0132309, 0.0109902, 0.0105634,
- 0.00960307, 0.00992318, 0.00960307, 0.00832266, 0.00821596,
- 0.00778916, 0.00661545, 0.00565514, 0.00554844, 0.00501494,
- 0.00597525, 0.00618865, 0.00426803, 0.00469484, 0.00437473,
- 0.00426803, 0.00597525, 0.00448143, 0.00309432, 0.00373453,
- 0.00362783, 0.00373453, 0.00330773, 0.00362783, 0.00224072,
- 0.00309432, 0.00266752, 0.00309432, 0.00245412, 0.00213402,
- 0.00202732, 0.00181391, 0.00138711, 0.00181391, 0.00245412,
- 0.00288092, 0.00170721, 0.00138711, 0.00202732, 0.00138711,
- 0.00149381, 0.00138711, 0.00192061, 0.00106701, 0.00117371,
- 0.00213402, 0.00181391, 0.00128041, 0.00117371, 0.00117371,
- 0.00117371, 0.000960307, 0.000960307, 0.000746906, 0.000853606,
- 0.000426803, 0.00106701, 0.000960307, 0.000426803, 0.000960307,
- 0.00117371, 0.00138711, 0.000853606, 0.000853606, 0.000426803,
- 0.000426803, 0.000960307, 0.000426803, 0.000640205, 0.00106701,
- 0.000746906, 0.000320102, 0.000320102, 0.000533504, 0,
- 0.000426803, 0.000426803, 0.000640205, 0.000746906, 0.000213402,
- 0.000106701, 0.000213402, 0.000426803, 0.000533504, 0.000640205,
- 0.000426803, 0.000106701, 0.000213402, 0.000320102, 0.000213402,
- 0.000213402, 0.000320102, 0.000426803, 0.000213402, 0.000106701,
- 0.000320102, 0.000213402, 0, 0.000106701, 0.000320102,
- 0.000426803, 0.000746906, 0.000533504, 0.000426803, 0.000426803,
- 0.000106701, 0, 0.000320102, 0.000213402, 0.000106701,
- 0.000320102, 0.000213402, 0, 0.000213402, 0.000320102,
- 0.000213402, 0, 0.000320102, 0.000213402, 0.000213402,
- 0.000213402, 0.000106701, 0, 0.000106701, 0.000213402,
- 0.000106701, 0.000106701, 0, 0.000213402, 0.000106701,
- 0.000106701, 0.000320102, 0.000106701, 0, 0,
- 0, 0.000106701, 0.000213402, 0.000106701, 0.000213402,
- 0, 0.000106701, 0, 0.000213402, 0,
- 0.000106701, 0.000106701, 0, 0, 0.000320102
- }
- };
-
- Double_t ProbEL[kNo_Mom][kNo_EnBins]={
- {
- 0, 0, 0, 0, 0,
- 0.000136166, 0, 0, 0, 0,
- 0, 0, 0, 0.000136166, 0,
- 0, 0.000544662, 0.000136166, 0.000816993, 0.000816993,
- 0.00190632, 0.0043573, 0.00449346, 0.00585512, 0.00844227,
- 0.00925926, 0.0119826, 0.0178377, 0.0177015, 0.0182462,
- 0.0212418, 0.0223312, 0.0249183, 0.0275054, 0.0251906,
- 0.0280501, 0.0265523, 0.0280501, 0.0337691, 0.025463,
- 0.0272331, 0.0290033, 0.0250545, 0.0268246, 0.025463,
- 0.0226035, 0.0241013, 0.0208333, 0.020561, 0.0224673,
- 0.0201525, 0.0191993, 0.0191993, 0.0163399, 0.0157952,
- 0.016476, 0.0144336, 0.0115741, 0.0115741, 0.0117102,
- 0.0115741, 0.0114379, 0.00994009, 0.0087146, 0.00748911,
- 0.00844227, 0.00708061, 0.00803377, 0.00735294, 0.00599129,
- 0.00517429, 0.00544662, 0.00626362, 0.00571895, 0.00503813,
- 0.00708061, 0.00544662, 0.00490196, 0.00367647, 0.0047658,
- 0.00340414, 0.00326797, 0.00285948, 0.00245098, 0.00367647,
- 0.00258715, 0.00340414, 0.00204248, 0.00313181, 0.00204248,
- 0.00190632, 0.00136166, 0.00204248, 0.00177015, 0.00177015,
- 0.00231481, 0.00217865, 0.00177015, 0.00163399, 0.00136166,
- 0.00149782, 0.00177015, 0.00108932, 0.00149782, 0.00149782,
- 0.00122549, 0.000953159, 0.00108932, 0.000408497, 0.000544662,
- 0.00177015, 0.00177015, 0.00122549, 0.000544662, 0.000272331,
- 0.000953159, 0.000953159, 0.00122549, 0.000953159, 0.000544662,
- 0.000680828, 0.00108932, 0.000680828, 0.000136166, 0.000544662,
- 0.000680828, 0.000680828, 0.000544662, 0.000408497, 0.000136166,
- 0, 0.000136166, 0.000544662, 0.000136166, 0.000544662,
- 0, 0.000272331, 0.000408497, 0.000272331, 0,
- 0.000408497, 0.000408497, 0.000544662, 0.000136166, 0.000272331,
- 0.000136166, 0, 0.000136166, 0.000544662, 0.000272331,
- 0.000680828, 0.000136166, 0.000136166, 0, 0.000272331,
- 0.000136166, 0.000136166, 0.000408497, 0.000136166, 0.000272331,
- 0.000272331, 0, 0.000408497, 0, 0.000136166,
- 0.000136166, 0.000680828, 0.000272331, 0.000136166, 0.000272331,
- 0.000544662, 0, 0, 0, 0,
- 0.000136166, 0.000136166, 0.000272331, 0.000136166, 0,
- 0.000272331, 0, 0, 0, 0.000408497,
- 0.000136166, 0, 0, 0.000136166, 0.000408497,
- 0, 0, 0.000272331, 0, 0,
- 0, 0, 0, 0.000136166, 0.000136166
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0, 0,
- 0, 0.000368189, 0.00012273, 0, 0.000981836,
- 0.0012273, 0.00257732, 0.00441826, 0.00454099, 0.00552283,
- 0.00871379, 0.0125184, 0.0126411, 0.0159548, 0.0162003,
- 0.0174276, 0.0197595, 0.0206186, 0.0247914, 0.0228277,
- 0.0234413, 0.0238095, 0.0250368, 0.0277369, 0.028596,
- 0.0249141, 0.0251595, 0.0268778, 0.0250368, 0.0239323,
- 0.0214777, 0.0214777, 0.0218459, 0.0260187, 0.0201276,
- 0.0212322, 0.0192685, 0.0174276, 0.0197595, 0.0182867,
- 0.0155866, 0.0169367, 0.0150957, 0.0132548, 0.0142366,
- 0.0127639, 0.0133775, 0.0115366, 0.00945017, 0.0110457,
- 0.00969563, 0.0123957, 0.00871379, 0.00711831, 0.00797742,
- 0.00785469, 0.00638193, 0.00760923, 0.00662739, 0.00613648,
- 0.00687285, 0.00466372, 0.0054001, 0.00503191, 0.00405007,
- 0.00527737, 0.00368189, 0.0033137, 0.0033137, 0.00380461,
- 0.00270005, 0.00245459, 0.0041728, 0.00282278, 0.00306824,
- 0.00220913, 0.0020864, 0.00294551, 0.00233186, 0.0020864,
- 0.0020864, 0.00184094, 0.0020864, 0.00184094, 0.00196367,
- 0.00171821, 0.0012273, 0.00110457, 0.000981836, 0.00171821,
- 0.00171821, 0.0012273, 0.00147275, 0.0012273, 0.00135002,
- 0.000981836, 0.000981836, 0.000981836, 0.00012273, 0.000859107,
- 0.00135002, 0.000859107, 0.000859107, 0.000736377, 0.000859107,
- 0.000245459, 0.000490918, 0.000613648, 0.000736377, 0.000613648,
- 0.000490918, 0.000613648, 0.000613648, 0.0012273, 0.000736377,
- 0.000368189, 0.00012273, 0.000490918, 0.000368189, 0.000245459,
- 0.000490918, 0.000490918, 0.00012273, 0.00012273, 0.00012273,
- 0.000490918, 0.000368189, 0.000368189, 0.000245459, 0.000245459,
- 0.00012273, 0.000736377, 0.000736377, 0.000368189, 0.000613648,
- 0.000613648, 0.000245459, 0.000613648, 0.000613648, 0.000245459,
- 0.000736377, 0.00012273, 0.000245459, 0, 0.000245459,
- 0.00012273, 0.000245459, 0.000245459, 0.00012273, 0.00012273,
- 0.000368189, 0, 0.000368189, 0.00012273, 0,
- 0.00012273, 0.00012273, 0.000368189, 0.000490918, 0.00012273,
- 0.000368189, 0, 0.000245459, 0.000245459, 0.000368189,
- 0.000245459, 0, 0, 0.000490918, 0.00012273,
- 0.000245459, 0.000245459, 0, 0.00012273, 0,
- 0, 0, 0.000368189, 0, 0,
- 0.00012273, 0.00012273, 0, 0, 0.00012273
- },
- {
- 0, 0, 0, 0, 0.000121477,
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0.000121477, 0,
- 0.000121477, 0, 0, 0.000364431, 0.000971817,
- 0.00133625, 0.00194363, 0.00133625, 0.0042517, 0.00655977,
- 0.00716715, 0.00947522, 0.013484, 0.0131195, 0.0136054,
- 0.0168853, 0.015792, 0.0183431, 0.0187075, 0.0202867,
- 0.0221088, 0.0253887, 0.0228377, 0.023931, 0.0247813,
- 0.0223518, 0.0250243, 0.024174, 0.0246599, 0.0224733,
- 0.0200437, 0.0208941, 0.0206511, 0.0210155, 0.0216229,
- 0.0221088, 0.0199223, 0.0196793, 0.0205296, 0.0161565,
- 0.0183431, 0.0179786, 0.0168853, 0.0154276, 0.0179786,
- 0.0154276, 0.013484, 0.0128766, 0.0115403, 0.0109329,
- 0.00971817, 0.00935374, 0.00947522, 0.00935374, 0.00911079,
- 0.00898931, 0.00741011, 0.00668124, 0.00716715, 0.00558795,
- 0.005345, 0.00862488, 0.00570943, 0.00485909, 0.00741011,
- 0.00400875, 0.00510204, 0.00498056, 0.00546647, 0.00376579,
- 0.00461613, 0.00388727, 0.00449466, 0.00400875, 0.00206511,
- 0.00315841, 0.0026725, 0.00279397, 0.00218659, 0.00315841,
- 0.00230807, 0.00194363, 0.00206511, 0.00315841, 0.00206511,
- 0.00315841, 0.00206511, 0.00194363, 0.00133625, 0.00242954,
- 0.000971817, 0.0015792, 0.00194363, 0.00109329, 0.00109329,
- 0.00133625, 0.00109329, 0.0015792, 0.000728863, 0.00170068,
- 0.000971817, 0.00121477, 0.00085034, 0.00109329, 0.00121477,
- 0.00109329, 0.000607386, 0.000607386, 0.000607386, 0.000728863,
- 0.000728863, 0.000971817, 0.00121477, 0.00121477, 0.000364431,
- 0.000485909, 0.000242954, 0.000242954, 0.000607386, 0.000242954,
- 0.000121477, 0.000607386, 0.000728863, 0.000485909, 0.000242954,
- 0.000607386, 0.000607386, 0.000242954, 0.000242954, 0.000364431,
- 0.000121477, 0.000364431, 0.000121477, 0.000364431, 0.000485909,
- 0.000728863, 0.000242954, 0.000364431, 0, 0.000121477,
- 0.000242954, 0, 0, 0.000121477, 0.000121477,
- 0.000364431, 0.000121477, 0.000364431, 0.000121477, 0.000364431,
- 0.000607386, 0.000242954, 0, 0.000242954, 0.000121477,
- 0.000242954, 0.000121477, 0.000242954, 0.000364431, 0.000121477,
- 0, 0.000121477, 0.000121477, 0.000242954, 0.000242954,
- 0.000121477, 0, 0.000242954, 0.000364431, 0.000121477,
- 0.000364431, 0, 0, 0.000121477, 0.000121477,
- 0.000242954, 0, 0.000121477, 0.000121477, 0,
- 0.000242954, 0.000242954, 0.000242954, 0.000121477, 0.000121477
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0, 0.00011558,
- 0.00011558, 0, 0, 0, 0,
- 0, 0.00011558, 0, 0.00023116, 0.000693481,
- 0.00150254, 0.00104022, 0.00184928, 0.00369857, 0.00485437,
- 0.00612575, 0.00889968, 0.0102866, 0.0117892, 0.0122515,
- 0.0146787, 0.0162968, 0.0161812, 0.0183773, 0.018724,
- 0.0190707, 0.0195331, 0.0186084, 0.0203421, 0.0226537,
- 0.0210356, 0.0219602, 0.0211512, 0.0230005, 0.0220758,
- 0.0205733, 0.0210356, 0.0208044, 0.0199954, 0.0193019,
- 0.017337, 0.0175682, 0.0198798, 0.0198798, 0.0194175,
- 0.0177994, 0.0160656, 0.0179149, 0.0158345, 0.0159501,
- 0.012945, 0.0124827, 0.0134073, 0.0142164, 0.0121359,
- 0.011558, 0.0114424, 0.00959316, 0.0101711, 0.0104022,
- 0.010749, 0.00809061, 0.0082062, 0.0082062, 0.0104022,
- 0.0087841, 0.00785945, 0.00693481, 0.00647249, 0.00705039,
- 0.00520111, 0.00531669, 0.00589459, 0.00554785, 0.00716597,
- 0.00612575, 0.00520111, 0.00508553, 0.00358299, 0.00404531,
- 0.00358299, 0.00346741, 0.00369857, 0.00369857, 0.00254276,
- 0.00416089, 0.00208044, 0.00323625, 0.00312067, 0.00208044,
- 0.00277393, 0.00277393, 0.00208044, 0.00219602, 0.00196486,
- 0.00300509, 0.00265834, 0.00208044, 0.00161812, 0.00127138,
- 0.00127138, 0.00150254, 0.0017337, 0.000924642, 0.00184928,
- 0.000809061, 0.00150254, 0.000924642, 0.0011558, 0.000577901,
- 0.00104022, 0.000924642, 0.00161812, 0.0011558, 0.00184928,
- 0.0011558, 0.000924642, 0.00150254, 0.000809061, 0.00127138,
- 0.00104022, 0.000346741, 0.000809061, 0.000462321, 0.000809061,
- 0.000924642, 0.000346741, 0.000346741, 0.000577901, 0.000809061,
- 0.0011558, 0.000346741, 0.000346741, 0.000462321, 0.000346741,
- 0.000346741, 0.000577901, 0.00023116, 0, 0.00023116,
- 0.000577901, 0.000346741, 0.00011558, 0.000346741, 0.000346741,
- 0.000346741, 0.00023116, 0.000693481, 0.00023116, 0.00023116,
- 0.000693481, 0.000346741, 0.000346741, 0.00023116, 0.00011558,
- 0.000462321, 0.00023116, 0.00011558, 0.000346741, 0.00011558,
- 0.00011558, 0.000462321, 0, 0.00023116, 0.00023116,
- 0.00023116, 0.00011558, 0.00023116, 0.00011558, 0.00011558,
- 0.00011558, 0.000462321, 0.00011558, 0.00011558, 0.00011558,
- 0.00023116, 0.00023116, 0, 0.00011558, 0,
- 0, 0.000346741, 0.00011558, 0.00011558, 0.00011558,
- 0.000346741, 0, 0.00011558, 0.000346741, 0
- },
- {
- 0, 0, 0, 0, 0,
- 0.000112309, 0, 0.000112309, 0, 0,
- 0, 0, 0, 0.000224618, 0.000112309,
- 0, 0, 0, 0.000224618, 0.000336927,
- 0.000786164, 0.00213387, 0.00213387, 0.00359389, 0.00460467,
- 0.00539084, 0.00730009, 0.00977089, 0.0104447, 0.0102201,
- 0.0134771, 0.013814, 0.0137017, 0.019991, 0.0149371,
- 0.0203279, 0.0203279, 0.0189802, 0.0203279, 0.0204403,
- 0.0188679, 0.0197664, 0.0196541, 0.021451, 0.019991,
- 0.0224618, 0.0184187, 0.0219003, 0.0203279, 0.0194295,
- 0.0196541, 0.0213387, 0.0162848, 0.0192049, 0.021788,
- 0.015611, 0.0151617, 0.0175202, 0.0162848, 0.0161725,
- 0.0134771, 0.0166217, 0.0133648, 0.0103324, 0.0132525,
- 0.0111186, 0.0128032, 0.0112309, 0.0103324, 0.0102201,
- 0.0098832, 0.0113432, 0.00876011, 0.00876011, 0.00977089,
- 0.00943396, 0.00774933, 0.00673854, 0.00797394, 0.00651393,
- 0.00673854, 0.00628931, 0.00426774, 0.0049416, 0.00561545,
- 0.00595238, 0.00561545, 0.00572776, 0.00426774, 0.00426774,
- 0.00381851, 0.0037062, 0.00336927, 0.00359389, 0.00336927,
- 0.00426774, 0.00359389, 0.00303235, 0.00269542, 0.00348158,
- 0.00404313, 0.00213387, 0.00190925, 0.00134771, 0.00213387,
- 0.00202156, 0.00213387, 0.00190925, 0.00213387, 0.00303235,
- 0.00224618, 0.00280773, 0.0024708, 0.0024708, 0.0012354,
- 0.00190925, 0.00146002, 0.00146002, 0.00101078, 0.00134771,
- 0.00112309, 0.00146002, 0.00134771, 0.00134771, 0.00179695,
- 0.0012354, 0.00101078, 0.000786164, 0.000561545, 0.0012354,
- 0.00101078, 0.000336927, 0.000673854, 0.000673854, 0.000673854,
- 0.000449236, 0.000224618, 0.000786164, 0.000336927, 0.000673854,
- 0.000561545, 0.000786164, 0.000898473, 0.000561545, 0.000561545,
- 0.000561545, 0.000449236, 0.000224618, 0.000112309, 0.00101078,
- 0.000336927, 0.000449236, 0.000449236, 0.000336927, 0.000561545,
- 0.000336927, 0.000336927, 0.000224618, 0.000112309, 0.000336927,
- 0.000336927, 0.000224618, 0.000112309, 0.000224618, 0.000673854,
- 0.000561545, 0.000336927, 0.000449236, 0.000449236, 0,
- 0.000449236, 0.000112309, 0.000224618, 0.000561545, 0.000224618,
- 0.000112309, 0.000224618, 0.000449236, 0.000112309, 0,
- 0, 0, 0.000224618, 0.000336927, 0.000561545,
- 0.000112309, 0.000112309, 0.000112309, 0.000224618, 0,
- 0, 0.000112309, 0.000112309, 0.000224618, 0,
- 0.000112309, 0.000224618, 0.000112309, 0.000112309, 0.000224618
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0, 0,
- 0, 0, 0, 0.000114077, 0,
- 0, 0.000228154, 0.000228154, 0.000228154, 0.000456308,
- 0.000912617, 0.00102669, 0.00296601, 0.00273785, 0.00456308,
- 0.00684463, 0.0065024, 0.00787132, 0.0126626, 0.0120922,
- 0.0119781, 0.0119781, 0.0138033, 0.0149441, 0.0143737,
- 0.017682, 0.016199, 0.0185946, 0.0216747, 0.0199635,
- 0.0197353, 0.0187086, 0.0208761, 0.022131, 0.0199635,
- 0.0215606, 0.0201916, 0.0187086, 0.0215606, 0.0172256,
- 0.0196213, 0.0180242, 0.0164271, 0.0183664, 0.0165412,
- 0.017682, 0.0195072, 0.0151723, 0.0141456, 0.0150582,
- 0.0150582, 0.0165412, 0.0156286, 0.0136893, 0.0125485,
- 0.0144878, 0.011864, 0.011864, 0.0110655, 0.011864,
- 0.00958248, 0.00992471, 0.0110655, 0.00935432, 0.00992471,
- 0.00821355, 0.00730094, 0.0065024, 0.00809948, 0.00581793,
- 0.00764317, 0.00570386, 0.00501939, 0.00581793, 0.00558978,
- 0.0065024, 0.00558978, 0.00479124, 0.0039927, 0.00456308,
- 0.00387862, 0.00353639, 0.00456308, 0.00376454, 0.00456308,
- 0.00296601, 0.00353639, 0.00342231, 0.00285193, 0.00296601,
- 0.00308008, 0.00228154, 0.00239562, 0.00353639, 0.00342231,
- 0.00273785, 0.0025097, 0.00216747, 0.00182523, 0.00193931,
- 0.00273785, 0.00171116, 0.00136893, 0.00228154, 0.001483,
- 0.00136893, 0.00114077, 0.001483, 0.00159708, 0.00079854,
- 0.00159708, 0.00205339, 0.001483, 0.00102669, 0.00125485,
- 0.001483, 0.00159708, 0.00114077, 0.00114077, 0.001483,
- 0.00136893, 0.00079854, 0.00159708, 0.00114077, 0.00102669,
- 0.000912617, 0.00102669, 0.00079854, 0.00079854, 0.000570386,
- 0.00102669, 0.000570386, 0.00102669, 0.000912617, 0.000456308,
- 0.000456308, 0.000912617, 0.000228154, 0.000570386, 0.000912617,
- 0.000456308, 0.000570386, 0.000570386, 0.000684463, 0.000114077,
- 0.000570386, 0.000114077, 0.000570386, 0.000228154, 0.000456308,
- 0.000228154, 0.000684463, 0.000342231, 0.000228154, 0.000570386,
- 0.000342231, 0.000114077, 0.000342231, 0.000342231, 0.000228154,
- 0.000228154, 0, 0.000114077, 0.000342231, 0.000228154,
- 0.000456308, 0.000342231, 0.000342231, 0.000114077, 0.000228154,
- 0.000456308, 0.000114077, 0, 0.000228154, 0.000228154,
- 0.000228154, 0.000114077, 0.000114077, 0.000114077, 0.000114077,
- 0, 0.000342231, 0.000228154, 0.000342231, 0.000456308,
- 0.000228154, 0, 0.000228154, 0.000342231, 0.000228154
- },
- {
- 0, 0, 0, 0, 0,
- 0, 0.000108155, 0, 0, 0.000324465,
- 0.000108155, 0, 0, 0, 0,
- 0.000108155, 0.000108155, 0, 0.000648929, 0.00021631,
- 0.000973394, 0.000865239, 0.00162232, 0.00205494, 0.00378542,
- 0.00529959, 0.00735453, 0.00735453, 0.00854424, 0.0107073,
- 0.0122215, 0.0130867, 0.0136275, 0.0144928, 0.0155743,
- 0.0146009, 0.0162232, 0.0178456, 0.0177374, 0.0186026,
- 0.0178456, 0.0186026, 0.0178456, 0.0184945, 0.0193597,
- 0.0187108, 0.0200087, 0.0173048, 0.0182782, 0.0189271,
- 0.0183863, 0.0177374, 0.0192516, 0.0201168, 0.0163314,
- 0.0175211, 0.0178456, 0.0144928, 0.0147091, 0.0171966,
- 0.0171966, 0.0137357, 0.012546, 0.0147091, 0.0146009,
- 0.0112481, 0.0124378, 0.01114, 0.0124378, 0.0114644,
- 0.010491, 0.0107073, 0.010491, 0.00897685, 0.0103829,
- 0.00789531, 0.00692191, 0.0100584, 0.00854424, 0.00724638,
- 0.00940947, 0.0067056, 0.00616483, 0.00681376, 0.0067056,
- 0.00519143, 0.0067056, 0.00746269, 0.00497512, 0.00562405,
- 0.00486697, 0.00540774, 0.00465066, 0.0045425, 0.00367727,
- 0.00378542, 0.00529959, 0.00378542, 0.00400173, 0.00302834,
- 0.00346096, 0.00259572, 0.00324465, 0.00302834, 0.0033528,
- 0.00237941, 0.00292018, 0.00183863, 0.00248756, 0.00259572,
- 0.00183863, 0.00227125, 0.0021631, 0.00259572, 0.00205494,
- 0.0011897, 0.00162232, 0.0011897, 0.00151417, 0.00183863,
- 0.00140601, 0.00129786, 0.00194679, 0.000757084, 0.000973394,
- 0.00205494, 0.00162232, 0.000757084, 0.0011897, 0.000973394,
- 0.00129786, 0.000973394, 0.000973394, 0.000865239, 0.0011897,
- 0.00108155, 0.000865239, 0.000648929, 0.0011897, 0.000540774,
- 0.00043262, 0.00108155, 0.000757084, 0.00129786, 0.00129786,
- 0.00021631, 0.000540774, 0.000540774, 0.00108155, 0.000973394,
- 0.00021631, 0.000540774, 0.00043262, 0.000648929, 0.000324465,
- 0.000757084, 0.000865239, 0.000108155, 0.00021631, 0.000324465,
- 0.00043262, 0.00043262, 0.000324465, 0.000540774, 0.00043262,
- 0.000757084, 0.000324465, 0.00043262, 0.000324465, 0.00043262,
- 0.00021631, 0.00043262, 0.000540774, 0.000324465, 0.000324465,
- 0.00021631, 0.00021631, 0.00021631, 0.00043262, 0.000324465,
- 0.000108155, 0.000324465, 0.000648929, 0.00021631, 0.000324465,
- 0.000540774, 0.000324465, 0.000324465, 0, 0.000108155,
- 0, 0.00021631, 0, 0.00021631, 0.00043262,
- 0.00021631, 0.000108155, 0.000108155, 0.000324465, 0.000108155
- }
- };
-
- // 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.109477, 0.0740741, 0.0750272, 0.0619553 ,
- 0.0642702, 0.0556917, 0.0586874, 0.057598, 0.059232 ,
- 0.0562364, 0.0551471, 0.0559641, 0.0514706, 0.0580065 ,
- 0.0656318, 0.0354031, 0.00558279, 0.000408497, 0.000136166
- },
- {
- 0, 0.113467, 0.0729167, 0.0715526, 0.0672123 ,
- 0.063244, 0.0613839, 0.0613839, 0.0561756, 0.0598958 ,
- 0.0555556, 0.0560516, 0.0533234, 0.0544395, 0.0549355 ,
- 0.0603919, 0.0336062, 0.00409226, 0.000124008, 0.000248016
- },
- {
- 0, 0.106427, 0.0816632, 0.0712969, 0.0685326 ,
- 0.0642709, 0.0590878, 0.0556323, 0.0610458, 0.0604699 ,
- 0.0542502, 0.0540198, 0.0550564, 0.0530984, 0.0551716 ,
- 0.0630039, 0.0313292, 0.00529832, 0.000345543, 0
- },
- {
- 0, 0.110388, 0.0731945, 0.0703752, 0.0597484 ,
- 0.0633268, 0.0605075, 0.06311, 0.0628931, 0.0570375 ,
- 0.0583388, 0.0530254, 0.0562785, 0.0541097, 0.0564953 ,
- 0.0612665, 0.0327478, 0.00639774, 0.000759054, 0
- },
- {
- 0, 0.116224, 0.0754374, 0.0736196, 0.066803 ,
- 0.0643036, 0.0601, 0.0591911, 0.0562372, 0.0594183 ,
- 0.05192, 0.0548739, 0.0514656, 0.0548739, 0.0566917 ,
- 0.0587367, 0.0341968, 0.00499886, 0.000908884, 0
- },
- {
- 0, 0.111563, 0.0715899, 0.0756549, 0.0667344 ,
- 0.0612014, 0.0615402, 0.0602981, 0.0605239, 0.0591689 ,
- 0.0569106, 0.0586043, 0.0539747, 0.0546522, 0.0538618 ,
- 0.0560072, 0.030262, 0.00654923, 0.000677507, 0.000225836
- },
- {
- 0, 0.111289, 0.0736236, 0.0735169, 0.0694622 ,
- 0.0606061, 0.0660478, 0.0591122, 0.0562313, 0.0601793 ,
- 0.0536705, 0.0545241, 0.0543107, 0.0565514, 0.0558045 ,
- 0.056338, 0.0315834, 0.00629535, 0.000533504, 0.000320102
- }
- };
-
- Double_t ProbELT[kNo_Mom][kNo_TBins]={
- {
- 0, 0.0834695, 0.0619553, 0.0551471, 0.0518791,
- 0.0517429, 0.0480664, 0.0510621, 0.0454793, 0.0556917,
- 0.0525599, 0.0574619, 0.0537854, 0.0690359, 0.0804739,
- 0.0946351, 0.0733932, 0.0132081, 0.000680828, 0.000272331
- },
- {
- 0, 0.079406, 0.04946, 0.0544919, 0.0527737,
- 0.0472509, 0.0490918, 0.0468827, 0.0470054, 0.0536328,
- 0.0509327, 0.0574374, 0.0587874, 0.0709377, 0.0801424,
- 0.106161, 0.0796514, 0.0155866, 0.000368189, 0
- },
- {
- 0, 0.0762877, 0.0569728, 0.0474976, 0.0420311,
- 0.0472546, 0.0462828, 0.0483479, 0.047862, 0.0521137,
- 0.050656, 0.0549077, 0.0614674, 0.0702138, 0.0860058,
- 0.108844, 0.0895287, 0.013241, 0.000364431, 0.000121477
- },
- {
- 0, 0.0736246, 0.0527046, 0.0496995, 0.0487748,
- 0.0488904, 0.0499307, 0.0485437, 0.0472723, 0.0517799,
- 0.0520111, 0.053976, 0.063338, 0.0644938, 0.0815996,
- 0.106796, 0.0917707, 0.0145631, 0.00023116, 0
- },
- {
- 0, 0.0727763, 0.0495283, 0.0528976, 0.0486298,
- 0.0440252, 0.0457098, 0.0493037, 0.0480683, 0.0513252,
- 0.0518868, 0.0626685, 0.06031, 0.0752471, 0.0800764,
- 0.105571, 0.0882749, 0.0132525, 0.000449236, 0
- },
- {
- 0, 0.0702715, 0.0480265, 0.0455168, 0.0446042,
- 0.0450605, 0.0464294, 0.0513347, 0.0487109, 0.0495095,
- 0.0563541, 0.059206, 0.0641113, 0.0743783, 0.079854,
- 0.103012, 0.0951403, 0.0181383, 0.000114077, 0.000228154
- },
- {
- 0, 0.0721393, 0.0472637, 0.0466148, 0.0472637,
- 0.0459658, 0.048886, 0.0470474, 0.048886, 0.0498594,
- 0.0508328, 0.0591607, 0.0634869, 0.0714904, 0.0825222,
- 0.112157, 0.0893359, 0.0164395, 0.000540774, 0.000108155
- }
- };
-
- 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];
- }
+ probabilityT= fh1MaxTimBinPI[ip-1]->GetBinContent(iTBin) + slop*(mom-fTrackMomentum[ip-1]);
+ return probabilityT;
}
-
+ return probabilityT;
}
-/*
-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;
-}
-
-*/
-