//
// pt-distribution
//____________________________________________________________
+Double_t AliGenMUONlib::PtJpsiPP7000( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi pT
+//
+// pp 7 TeV
+// using ALICE data at 2.5<y<4, see arXiv:1103.2394
+
+ const Double_t kpt0 = 2.44;
+ const Double_t kxn = 3.9;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+0.36*(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtJpsiPP2760( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi pT
+//
+// pp 2.76 TeV
+// from the fit of RHIC + LHC data, see arXiv:1103.2394
+
+ const Double_t kpt0 = 2.31;
+ const Double_t kxn = 3.9;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+0.36*(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtJpsiPbPb2760( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi pT
+//
+// PbPb 2.76 TeV, for EKS98 with minimum bias shadowing factor 0.66
+//
+ Double_t c[5] = {6.01022e-01, 4.70988e-02, -2.27917e-03, 3.09885e-05, 1.31955e-06};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtJpsiPP2760(px,dummy);
+}
+
Double_t AliGenMUONlib::PtJpsi( const Double_t *px, const Double_t */*dummy*/)
{
// J/Psi pT
return x/TMath::Power(pass1,kxn);
}
-Double_t AliGenMUONlib::PtJpsiCDFscaledold( const Double_t *px, const Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP( const Double_t *px, const Double_t */*dummy*/)
{
// J/Psi pT
//
-// PbPb 5.5 TeV
-// scaled from CDF "old" data at 2 TeV
+// pp 14 TeV
+// scaled from CDF data at 2 TeV
- const Double_t kpt0 = 4.703;
- const Double_t kxn = 3.826;
+ const Double_t kpt0 = 5.630;
+ const Double_t kxn = 4.071;
Double_t x=*px;
//
Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
return x/TMath::Power(pass1,kxn);
}
-Double_t AliGenMUONlib::PtJpsiCDFscaledPP( const Double_t *px, const Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
{
// J/Psi pT
//
-// pp 14 TeV
+// pp 10 TeV
// scaled from CDF data at 2 TeV
- const Double_t kpt0 = 5.630;
+ const Double_t kpt0 = 5.334;
const Double_t kxn = 4.071;
Double_t x=*px;
//
return x/TMath::Power(pass1,kxn);
}
-Double_t AliGenMUONlib::PtJpsiCDFscaledPPold( const Double_t *px, const Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP9( const Double_t *px, const Double_t */*dummy*/)
{
// J/Psi pT
//
-// pp 14 TeV
-// scaled from CDF "old" data at 2 TeV
+// pp 8.8 TeV
+// scaled from CDF data at 2 TeV
+//
+ const Double_t kpt0 = 5.245;
+ const Double_t kxn = 4.071;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP7( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi pT
+//
+// pp 7 TeV
+// scaled from CDF data at 2 TeV
- const Double_t kpt0 = 5.355;
- const Double_t kxn = 3.821;
+ const Double_t kpt0 = 5.072;
+ const Double_t kxn = 4.071;
Double_t x=*px;
//
Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
return x/TMath::Power(pass1,kxn);
}
-Double_t AliGenMUONlib::PtJpsiCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP4( const Double_t *px, const Double_t */*dummy*/)
{
// J/Psi pT
//
-// pp 10 TeV
+// pp 3.94 TeV
// scaled from CDF data at 2 TeV
+//
+ const Double_t kpt0 = 4.647;
+ const Double_t kxn = 4.071;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
- const Double_t kpt0 = 5.334;
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP3( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi pT
+//
+// pp 2.76 TeV
+// scaled from CDF data at 1.9 TeV
+//
+ const Double_t kpt0 = 4.435;
+ const Double_t kxn = 4.071;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP2( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi pT
+//
+// pp 1.9 TeV
+// fit of the CDF data at 1.9 TeV
+//
+ const Double_t kpt0 = 4.233;
const Double_t kxn = 4.071;
Double_t x=*px;
//
return x/TMath::Power(pass1,kxn);
}
+Double_t AliGenMUONlib::PtJpsiCDFscaledPPb9( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi pT
+//
+// pPb 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.80
+//
+ Double_t c[5] = {6.42774e-01, 1.86168e-02, -6.77296e-04, 8.93512e-06, 1.31586e-07};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtJpsiCDFscaledPP9(px,dummy);
+}
+
+Double_t AliGenMUONlib::PtJpsiCDFscaledPbP9( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi pT
+//
+// Pbp 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.80
+//
+ Double_t c[5] = {8.58557e-01, 5.39791e-02, -4.75180e-03, 2.49463e-04, 5.52396e-05};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtJpsiCDFscaledPP9(px,dummy);
+}
+
+Double_t AliGenMUONlib::PtJpsiCDFscaledPbPb4( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi pT
+//
+// PbPb 3.94 TeV, for EKS98 with minimum bias shadowing factor 0.66
+//
+ Double_t c[5] = {6.01022e-01, 4.70988e-02, -2.27917e-03, 3.09885e-05, 1.31955e-06};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtJpsiCDFscaledPP4(px,dummy);
+}
+
Double_t AliGenMUONlib::PtJpsiFlat( const Double_t */*px*/, const Double_t */*dummy*/ )
{
return 1.;
//
// y-distribution
//____________________________________________________________
+Double_t AliGenMUONlib::YJpsiPP7000( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pp 7 TeV
+// from the fit of RHIC + LHC data, see arXiv:1103.2394
+//
+ Double_t x = px[0]/7.72;
+ x = x*x;
+ Double_t y = TMath::Exp(-x/0.383/0.383/2);
+ if(x > 1) y=0;
+ return y;
+}
+
+Double_t AliGenMUONlib::YJpsiPP2760( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pp 2.76 TeV
+// from the fit of RHIC + LHC data, see arXiv:1103.2394
+//
+ Double_t x = px[0]/6.79;
+ x = x*x;
+ Double_t y = TMath::Exp(-x/0.383/0.383/2);
+ if(x > 1) y=0;
+ return y;
+}
+
+Double_t AliGenMUONlib::YJpsiPbPb2760( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi y
+//
+// PbPb 2.76 TeV, for EKS98 with minimum bias shadowing factor 0.66
+//
+ Double_t c[4] = {5.95228e-01, 9.45069e-03, 2.44710e-04, -1.32894e-05};
+ Double_t x = px[0]*px[0];
+ Double_t y;
+ Int_t j;
+ y = c[j = 3];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YJpsiPP2760(px,dummy);
+}
+
Double_t AliGenMUONlib::YJpsi(const Double_t *py, const Double_t */*dummy*/)
{
// J/psi y
return AliGenMUONlib::YJpsiPbPb(px, dummy);
}
-Double_t AliGenMUONlib::YJpsiCDFscaledold( const Double_t *px, const Double_t* dummy)
-{
- // J/Psi y
- return AliGenMUONlib::YJpsiPbPb(px, dummy);
-}
-
Double_t AliGenMUONlib::YJpsiCDFscaledPP( const Double_t *px, const Double_t* dummy)
{
// J/Psi y
return AliGenMUONlib::YJpsiPP(px, dummy);
}
-Double_t AliGenMUONlib::YJpsiCDFscaledPPold( const Double_t *px, const Double_t* dummy)
-{
- // J/Psi y
- return AliGenMUONlib::YJpsiPP(px, dummy);
-}
-
Double_t AliGenMUONlib::YJpsiCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
{
-
-//
// J/Psi y
//
// pp 10 TeV
return y;
}
+Double_t AliGenMUONlib::YJpsiCDFscaledPP9( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pp 8.8 TeV
+// rescaling of YJpsiPP(14 TeV) using 8.8 TeV / 14 TeV ratio of y-spectra in LO QCD
+//
+ Double_t c[5] = {3.33882e+02, -1.30980e+02, 2.59082e+01, -3.08935e+00, 1.56375e-01};
+ Double_t x = TMath::Abs(px[0]);
+ Double_t y;
+
+ if (x < 3.7) {
+ y = 99.236 - 1.5498 * x * x;
+ } else if (x < 7.4) {
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ } else {
+ y =0.;
+ }
+
+ if(y<0) y=0;
+
+ return y;
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPP9dummy(Double_t px)
+{
+ return AliGenMUONlib::YJpsiCDFscaledPP9(&px, (Double_t*) 0);
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPP7( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pp 7 TeV
+// scaled from YJpsiPP(14 TeV) using 7 TeV / 14 TeV ratio of y-spectra in LO pQCD.
+//
+
+ Double_t c[5] = {6.71181e+02, -3.69240e+02, 8.89644e+01, -1.04937e+01, 4.80959e-01};
+
+ Double_t x = TMath::Abs(px[0]);
+ Double_t y;
+
+ if (x < 4.0) {
+ y = 100.78 - 1.8353 * x * x;
+ } else if (x < 7.3) {
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ } else {
+ y =0.;
+ }
+
+ if(y<0) y=0;
+
+ return y;
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPP4( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pp 3.94 TeV
+// rescaling of YJpsiPP(14 TeV) using 3.94 TeV / 14 TeV ratio of y-spectra in LO QCD
+//
+ Double_t c[5] = {4.00785e+02, -1.41159e+01, -3.28599e+01, 5.53048e+00, -2.45151e-01};
+ Double_t x = TMath::Abs(px[0]);
+ Double_t y;
+
+ if (x < 5.5) {
+ y = 107.389 - 2.7454 * x * x;
+ } else if (x < 7.0) {
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ } else {
+ y =0.;
+ }
+
+ if(y<0) y=0;
+
+ return y;
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPP3( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi y
+ return AliGenMUONlib::YJpsiPP2760(px, dummy);
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPP2( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pp 1.9 TeV
+// from the fit of RHIC + LHC data, see arXiv:1103.2394
+//
+ Double_t x = px[0]/6.42;
+ x = x*x;
+ Double_t y = TMath::Exp(-x/0.383/0.383/2);
+ if(x > 1) y=0;
+ return y;
+}
+
Double_t AliGenMUONlib::YJpsiPP( const Double_t *px, const Double_t */*dummy*/)
{
return y;
}
+Double_t AliGenMUONlib::YJpsiCDFscaledPPb9( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// pPb 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.80
+//
+ Double_t c[7] = {7.52296e-01, 2.49917e-02, 3.36500e-03, 1.91187e-03, 2.92154e-04,
+ -4.16509e-05,-7.62709e-06};
+ Double_t y;
+ Double_t x = px[0] + 0.47; // rapidity shift
+ Int_t j;
+ y = c[j = 6];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YJpsiCDFscaledPP9dummy(x);
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPbP9( const Double_t *px, const Double_t */*dummy*/)
+{
+// J/Psi y
+//
+// Pbp 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.80
+//
+ Double_t c[7] = {7.52296e-01, 2.49917e-02, 3.36500e-03, 1.91187e-03, 2.92154e-04,
+ -4.16509e-05,-7.62709e-06};
+ Double_t y;
+ Double_t x = -px[0] + 0.47; // rapidity shift
+ Int_t j;
+ y = c[j = 6];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YJpsiCDFscaledPP9dummy(x);
+}
+
+Double_t AliGenMUONlib::YJpsiCDFscaledPbPb4( const Double_t *px, const Double_t *dummy)
+{
+// J/Psi y
+//
+// PbPb 3.94 TeV, for EKS98 with minimum bias shadowing factor 0.66
+//
+ Double_t c[4] = {5.95228e-01, 9.45069e-03, 2.44710e-04, -1.32894e-05};
+ Double_t x = px[0]*px[0];
+ Double_t y;
+ Int_t j;
+ y = c[j = 3];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YJpsiCDFscaledPP4(px,dummy);
+}
+
Double_t AliGenMUONlib::YJpsiBPbPb( const Double_t *px, const Double_t */*dummy*/)
{
return x/TMath::Power(pass1,kxn);
}
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon pT
+//
+// pp 10 TeV
+//
+// scaled from CDF data at 2 TeV
+
+ const Double_t kpt0 = 8.235;
+ const Double_t kxn = 3.051;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPP9( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon pT
+//
+// pp 8.8 TeV
+// scaled from CDF data at 2 TeV
+//
+ const Double_t kpt0 = 8.048;
+ const Double_t kxn = 3.051;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPP7( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon pT
+//
+// pp 7 TeV
+//
+// scaled from CDF data at 2 TeV
+
+ const Double_t kpt0 = 7.817;
+ const Double_t kxn = 3.051;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPP4( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon pT
+//
+// pp 3.94 TeV
+// scaled from CDF data at 2 TeV
+//
+ const Double_t kpt0 = 7.189;
+ const Double_t kxn = 3.051;
+ Double_t x=*px;
+ //
+ Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+ return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPPb9( const Double_t *px, const Double_t *dummy)
+{
+// Upsilon pT
+//
+// pPb 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.90
+//
+ Double_t c[5] = {7.64952e-01, 1.12501e-04, 4.96038e-04, -3.03198e-05, 3.74035e-06};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtUpsilonCDFscaledPP9(px,dummy);
+}
+
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPbP9( const Double_t *px, const Double_t *dummy)
+{
+// Upsilon pT
+//
+// Pbp 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.90
+//
+ Double_t c[5] = {1.09881e+00, 3.08329e-03, -2.00356e-04, 8.28991e-06, 2.52576e-06};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtUpsilonCDFscaledPP9(px,dummy);
+}
+
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPbPb4( const Double_t *px, const Double_t *dummy)
+{
+// Upsilon pT
+//
+// PbPb 3.94 TeV, for EKS98 with minimum bias shadowing factor 0.85
+//
+ Double_t c[5] = {8.65872e-01, 2.05465e-03, 2.56063e-04, -1.65598e-05, 2.29209e-06};
+ Double_t x=*px;
+ Double_t y;
+ Int_t j;
+ y = c[j = 4];
+ while (j > 0) y = y * x + c[--j];
+ //
+ Double_t d = 1.+c[4]*TMath::Power(x,4);
+ return y/d * AliGenMUONlib::PtUpsilonCDFscaledPP4(px,dummy);
+}
+
Double_t AliGenMUONlib::PtUpsilonFlat( const Double_t */*px*/, const Double_t */*dummy*/ )
{
return 1.;
Double_t AliGenMUONlib::PtUpsilonPbPb( const Double_t *px, const Double_t */*dummy*/)
{
-
//
// Upsilon pT
//
Double_t AliGenMUONlib::PtUpsilonPP( const Double_t *px, const Double_t */*dummy*/)
{
-
//
// Upsilon pT
//
Double_t c[7] = {3.40036e-01, -3.98882e-07, -4.48398e-03, 8.46411e-08, -6.10854e-04,
-2.99753e-09, 1.28895e-05};
-
- Double_t x = px[0];
- if (TMath::Abs(x) > 5.55) return 0.;
+ Double_t x = TMath::Abs(px[0]);
+ if (x > 5.55) return 0.;
Int_t j;
Double_t y = c[j = 6];
while (j > 0) y = y * x +c[--j];
}
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon y
+//
+// pp 10 TeV
+// scaled from YUpsilonPP(14 TeV) using 10 TeV / 14 TeV ratio of y-spectra in LO pQCD.
+// see S.Grigoryan, PWG3 Meeting, 27th Oct 2008
+//
+ Double_t c[4] = {1., -2.17877e-02, -6.52830e-04, 1.40578e-05};
+ Double_t x = TMath::Abs(px[0]);
+ if (x > 6.1) return 0.;
+ Int_t j;
+ Double_t y = c[j = 3];
+ while (j > 0) y = y * x*x +c[--j];
+ return y;
+}
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP9( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon y
+//
+// pp 8.8 TeV
+// rescaling of YUpsilonPP(14 TeV) using 8.8 TeV / 14 TeV ratio of y-spectra in LO QCD
+//
+ Double_t c[4] = {1., -2.37621e-02, -6.29610e-04, 1.47976e-05};
+ Double_t x = TMath::Abs(px[0]);
+ if (x > 6.1) return 0.;
+ Int_t j;
+ Double_t y = c[j = 3];
+ while (j > 0) y = y * x*x +c[--j];
+ return y;
+}
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP9dummy(Double_t px)
+{
+ return AliGenMUONlib::YUpsilonCDFscaledPP9(&px, (Double_t*) 0);
+}
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP7( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon y
+//
+// pp 7 TeV
+// scaled from YUpsilonPP(14 TeV) using 7 TeV / 14 TeV ratio of y-spectra in LO pQCD.
+//
+ Double_t c[4] = {1., -2.61009e-02, -6.83937e-04, 1.78451e-05};
+ Double_t x = TMath::Abs(px[0]);
+ if (x > 6.0) return 0.;
+ Int_t j;
+ Double_t y = c[j = 3];
+ while (j > 0) y = y * x*x +c[--j];
+ return y;
+}
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP4( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon y
+//
+// pp 3.94 TeV
+// rescaling of YUpsilonPP(14 TeV) using 3.94 TeV / 14 TeV ratio of y-spectra in LO QCD
+//
+ Double_t c[4] = {1., -3.91924e-02, -4.26184e-04, 2.10914e-05};
+ Double_t x = TMath::Abs(px[0]);
+ if (x > 5.7) return 0.;
+ Int_t j;
+ Double_t y = c[j = 3];
+ while (j > 0) y = y * x*x +c[--j];
+
+ return y;
+}
+
Double_t AliGenMUONlib::YUpsilonPP( const Double_t *px, const Double_t */*dummy*/)
{
//
Double_t c[7] = {8.91936e-01, -6.46645e-07, -1.52774e-02, 4.28677e-08, -7.01517e-04,
-6.20539e-10, 1.29943e-05};
-
- Double_t x = px[0];
- if (TMath::Abs(x) > 6.2) return 0.;
+ Double_t x = TMath::Abs(px[0]);
+ if (x > 6.2) return 0.;
Int_t j;
Double_t y = c[j = 6];
while (j > 0) y = y * x +c[--j];
return y;
}
+Double_t AliGenMUONlib::YUpsilonCDFscaledPPb9( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon y
+//
+// pPb 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.90
+//
+ Double_t c[7] = {8.71829e-01, 4.77467e-02, 8.09671e-03, 6.45294e-04, -2.15730e-04,
+ -4.67538e-05,-2.11683e-06};
+ Double_t y;
+ Double_t x = px[0] + 0.47; // rapidity shift
+ Int_t j;
+ y = c[j = 6];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YUpsilonCDFscaledPP9dummy(x);
+}
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPbP9( const Double_t *px, const Double_t */*dummy*/)
+{
+// Upsilon y
+//
+// Pbp 8.8 TeV, for EKS98 with minimum bias shadowing factor 0.90
+//
+ Double_t c[7] = {8.71829e-01, 4.77467e-02, 8.09671e-03, 6.45294e-04, -2.15730e-04,
+ -4.67538e-05,-2.11683e-06};
+ Double_t y;
+ Double_t x = -px[0] + 0.47; // rapidity shift
+ Int_t j;
+ y = c[j = 6];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YUpsilonCDFscaledPP9dummy(x);
+}
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPbPb4( const Double_t *px, const Double_t *dummy)
+{
+// Upsilon y
+//
+// PbPb 3.94 TeV, for EKS98 with minimum bias shadowing factor 0.85
+//
+ Double_t c[4] = {8.27837e-01, 1.70115e-02, -1.26046e-03, 1.52091e-05};
+ Double_t x = px[0]*px[0];
+ Double_t y;
+ Int_t j;
+ y = c[j = 3];
+ while (j > 0) y = y * x + c[--j];
+ if(y<0) y=0;
+
+ return y * AliGenMUONlib::YUpsilonCDFscaledPP4(px,dummy);
+}
+
+
// particle composition
//
Int_t AliGenMUONlib::IpUpsilon(TRandom *)
break;
case kJpsiFamily:
case kPsiP:
+ case kChic1:
+ case kChic2:
case kJpsi:
if (sname == "Vogt" || sname == "Vogt PbPb") {
func=PtJpsiPbPb;
} else if (sname == "Vogt pp") {
func=PtJpsiPP;
+ } else if (sname == "pp 7") {
+ func=PtJpsiPP7000;
+ } else if (sname == "pp 2.76") {
+ func=PtJpsiPP2760;
+ } else if (sname == "PbPb 2.76") {
+ func=PtJpsiPbPb2760;
} else if (sname == "CDF scaled") {
func=PtJpsiCDFscaled;
} else if (sname == "CDF pp") {
func=PtJpsiCDFscaledPP;
} else if (sname == "CDF pp 10") {
func=PtJpsiCDFscaledPP10;
- } else if (sname == "CDF scaled old") {
- func=PtJpsiCDFscaledold;
- } else if (sname == "CDF pp old") {
- func=PtJpsiCDFscaledPPold;
- } else if (sname == "Flat") {
+ } else if (sname == "CDF pp 8.8") {
+ func=PtJpsiCDFscaledPP9;
+ } else if (sname == "CDF pp 7" || sname == "CDF pp 7 flat y") {
+ func=PtJpsiCDFscaledPP7;
+ } else if (sname == "CDF pp 3.94") {
+ func=PtJpsiCDFscaledPP4;
+ } else if (sname == "CDF pp 2.76") {
+ func=PtJpsiCDFscaledPP3;
+ } else if (sname == "CDF pp 1.9") {
+ func=PtJpsiCDFscaledPP2;
+ } else if (sname == "CDF pPb 8.8") {
+ func=PtJpsiCDFscaledPPb9;
+ } else if (sname == "CDF Pbp 8.8") {
+ func=PtJpsiCDFscaledPbP9;
+ } else if (sname == "CDF PbPb 3.94") {
+ func=PtJpsiCDFscaledPbPb4;
+ } else if (sname == "Flat" || sname == "CDF pp 7 flat pt") {
func=PtJpsiFlat;
} else {
func=PtJpsi;
func=PtUpsilonCDFscaled;
} else if (sname == "CDF pp") {
func=PtUpsilonCDFscaledPP;
+ } else if (sname == "CDF pp 10") {
+ func=PtUpsilonCDFscaledPP10;
+ } else if (sname == "CDF pp 8.8") {
+ func=PtUpsilonCDFscaledPP9;
+ } else if (sname == "CDF pp 7") {
+ func=PtUpsilonCDFscaledPP7;
+ } else if (sname == "CDF pp 3.94") {
+ func=PtUpsilonCDFscaledPP4;
+ } else if (sname == "CDF pPb 8.8") {
+ func=PtUpsilonCDFscaledPPb9;
+ } else if (sname == "CDF Pbp 8.8") {
+ func=PtUpsilonCDFscaledPbP9;
+ } else if (sname == "CDF PbPb 3.94") {
+ func=PtUpsilonCDFscaledPbPb4;
} else if (sname == "Flat") {
func=PtUpsilonFlat;
} else {
case kChic0:
func=PtChic0;
break;
- case kChic1:
- func=PtChic1;
- break;
- case kChic2:
- func=PtChic2;
- break;
case kChic:
func=PtChic;
break;
break;
case kJpsiFamily:
case kPsiP:
+ case kChic1:
+ case kChic2:
case kJpsi:
if (sname == "Vogt" || sname == "Vogt PbPb") {
func=YJpsiPbPb;
} else if (sname == "Vogt pp"){
func=YJpsiPP;
+ } else if (sname == "pp 7") {
+ func=YJpsiPP7000;
+ } else if (sname == "pp 2.76") {
+ func=YJpsiPP2760;
+ } else if (sname == "PbPb 2.76") {
+ func=YJpsiPbPb2760;
} else if (sname == "CDF scaled") {
func=YJpsiCDFscaled;
} else if (sname == "CDF pp") {
func=YJpsiCDFscaledPP;
} else if (sname == "CDF pp 10") {
func=YJpsiCDFscaledPP10;
- } else if (sname == "CDF scaled old") {
- func=YJpsiCDFscaledold;
- } else if (sname == "CDF pp old") {
- func=YJpsiCDFscaledPPold;
- } else if (sname == "Flat") {
+ } else if (sname == "CDF pp 8.8") {
+ func=YJpsiCDFscaledPP9;
+ } else if (sname == "CDF pp 7" || sname == "CDF pp 7 flat pt") {
+ func=YJpsiCDFscaledPP7;
+ } else if (sname == "CDF pp 3.94") {
+ func=YJpsiCDFscaledPP4;
+ } else if (sname == "CDF pp 2.76") {
+ func=YJpsiCDFscaledPP3;
+ } else if (sname == "CDF pp 1.9") {
+ func=YJpsiCDFscaledPP2;
+ } else if (sname == "CDF pPb 8.8") {
+ func=YJpsiCDFscaledPPb9;
+ } else if (sname == "CDF Pbp 8.8") {
+ func=YJpsiCDFscaledPbP9;
+ } else if (sname == "CDF PbPb 3.94") {
+ func=YJpsiCDFscaledPbPb4;
+ } else if (sname == "Flat" || sname == "CDF pp 7 flat y") {
func=YJpsiFlat;
} else {
func=YJpsi;
func=YUpsilonCDFscaled;
} else if (sname == "CDF pp") {
func=YUpsilonCDFscaledPP;
+ } else if (sname == "CDF pp 10") {
+ func=YUpsilonCDFscaledPP10;
+ } else if (sname == "CDF pp 8.8") {
+ func=YUpsilonCDFscaledPP9;
+ } else if (sname == "CDF pp 7") {
+ func=YUpsilonCDFscaledPP7;
+ } else if (sname == "CDF pp 3.94") {
+ func=YUpsilonCDFscaledPP4;
+ } else if (sname == "CDF pPb 8.8") {
+ func=YUpsilonCDFscaledPPb9;
+ } else if (sname == "CDF Pbp 8.8") {
+ func=YUpsilonCDFscaledPbP9;
+ } else if (sname == "CDF PbPb 3.94") {
+ func=YUpsilonCDFscaledPbPb4;
} else if (sname == "Flat") {
func=YUpsilonFlat;
} else {
case kChic0:
func=YChic0;
break;
- case kChic1:
- func=YChic1;
- break;
- case kChic2:
- func=YChic2;
- break;
case kChic:
func=YChic;
break;