+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);
+}
+