From: morsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 13 Jul 2006 11:57:17 +0000 (+0000)
Subject: Quarkonia for pp add. (D. Stocco)
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=9ff138499e0da64e412c2443dcb3a94b8d955a63;hp=de0b6220a3907b24c91863b2b891301b686a5145

Quarkonia for pp add. (D. Stocco)
---

diff --git a/EVGEN/AliGenMUONlib.cxx b/EVGEN/AliGenMUONlib.cxx
index 4757a5e7f0e..be21e4d1fba 100644
--- a/EVGEN/AliGenMUONlib.cxx
+++ b/EVGEN/AliGenMUONlib.cxx
@@ -173,6 +173,22 @@ Double_t AliGenMUONlib::PtJpsiCDFscaled( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP( Double_t *px, Double_t */*dummy*/)
+{
+// J/Psi pT
+//
+// pp 14 TeV
+//
+// scaled from CDF data at 2 TeV
+
+  const Double_t kpt0 = 5.355;
+  const Double_t kxn  = 3.821;
+  Double_t x=*px;
+  //
+  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+  return x/TMath::Power(pass1,kxn);
+}
+
 Double_t AliGenMUONlib::PtJpsiFlat( Double_t */*px*/, Double_t */*dummy*/ )
 {
   return 1.;
@@ -302,6 +318,11 @@ Double_t AliGenMUONlib::YJpsiCDFscaled( Double_t *px, Double_t* dummy)
     return AliGenMUONlib::YJpsiPbPb(px, dummy);
 }
 
+Double_t AliGenMUONlib::YJpsiCDFscaledPP( Double_t *px, Double_t* dummy)
+{
+    // J/Psi y 
+    return AliGenMUONlib::YJpsiPP(px, dummy);
+}
 
 Double_t AliGenMUONlib::YJpsiPP( Double_t *px, Double_t */*dummy*/)
 {
@@ -414,6 +435,22 @@ Double_t AliGenMUONlib::PtUpsilonCDFscaled( Double_t *px, Double_t */*dummy*/ )
   return x/TMath::Power(pass1,kxn);
 }
 
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPP( Double_t *px, Double_t */*dummy*/ )
+{
+// Upsilon pT
+//
+// pp 14 TeV
+//
+// scaled from CDF data at 2 TeV
+
+  const Double_t kpt0 = 8.610;
+  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::PtUpsilonFlat( Double_t */*px*/, Double_t */*dummy*/ )
 {
   return 1.;
@@ -526,6 +563,14 @@ Double_t AliGenMUONlib::YUpsilonCDFscaled( Double_t *px, Double_t *dummy)
     return AliGenMUONlib::YUpsilonPbPb(px, dummy);
     
 }
+
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP( Double_t *px, Double_t *dummy)
+{
+    // Upsilon y
+    return AliGenMUONlib::YUpsilonPP(px, dummy);
+    
+}
+
 Double_t AliGenMUONlib::YUpsilonFlat( Double_t */*px*/, Double_t */*dummy*/)
 {
     // Upsilon y
@@ -860,6 +905,8 @@ GenFunc AliGenMUONlib::GetPt(Int_t param,  const char* tname) const
 	    func=PtJpsiPP;
 	} else if (sname == "CDF scaled") {
 	    func=PtJpsiCDFscaled;
+	} else if (sname == "CDF pp") {
+	    func=PtJpsiCDFscaledPP;
 	} else if (sname == "Flat") {
 	    func=PtJpsiFlat;
 	} else {
@@ -879,6 +926,8 @@ GenFunc AliGenMUONlib::GetPt(Int_t param,  const char* tname) const
 	    func=PtUpsilonPP;
 	} else if (sname == "CDF scaled") {
 	    func=PtUpsilonCDFscaled;
+	} else if (sname == "CDF pp") {
+	    func=PtUpsilonCDFscaledPP;
 	} else if (sname == "Flat") {
 	    func=PtUpsilonFlat;
 	} else {
@@ -939,6 +988,8 @@ GenFunc AliGenMUONlib::GetY(Int_t param, const char* tname) const
 	    func=YJpsiPP;
 	} else if (sname == "CDF scaled") {
 	    func=YJpsiCDFscaled;
+	} else if (sname == "CDF pp") {
+	    func=YJpsiCDFscaledPP;
 	} else if (sname == "Flat") {
 	    func=YJpsiFlat;
 	} else {
@@ -958,6 +1009,8 @@ GenFunc AliGenMUONlib::GetY(Int_t param, const char* tname) const
 	    func = YUpsilonPP;
 	} else if (sname == "CDF scaled") {
 	    func=YUpsilonCDFscaled;
+	} else if (sname == "CDF pp") {
+	    func=YUpsilonCDFscaledPP;
 	} else if (sname == "Flat") {
 	    func=YUpsilonFlat;
 	} else {
diff --git a/EVGEN/AliGenMUONlib.h b/EVGEN/AliGenMUONlib.h
index d2e24359614..0e2f74dbd7a 100644
--- a/EVGEN/AliGenMUONlib.h
+++ b/EVGEN/AliGenMUONlib.h
@@ -49,11 +49,13 @@ class AliGenMUONlib :
 // J/Psi     
     static Double_t PtJpsi( Double_t *px, Double_t *dummy);
     static Double_t PtJpsiCDFscaled( Double_t *px, Double_t *dummy);
+    static Double_t PtJpsiCDFscaledPP( Double_t *px, Double_t *dummy);
     static Double_t YJpsi(Double_t *py, Double_t *dummy);
     static Double_t PtJpsiPbPb( Double_t *px, Double_t *dummy);
     static Double_t PtJpsiBPbPb( Double_t *px, Double_t *dummy);
     static Double_t YJpsiPbPb(Double_t *py, Double_t *dummy);
     static Double_t YJpsiCDFscaled(Double_t *py, Double_t *dummy);
+    static Double_t YJpsiCDFscaledPP( Double_t *px, Double_t *dummy);
     static Double_t PtJpsiPP( Double_t *px, Double_t *dummy);
     static Double_t YJpsiPP(Double_t *py, Double_t *dummy);
     static Double_t YJpsiBPbPb(Double_t *py, Double_t *dummy);
@@ -66,8 +68,10 @@ class AliGenMUONlib :
 // Upsilon    
     static Double_t PtUpsilon( Double_t *px, Double_t *dummy );
     static Double_t PtUpsilonCDFscaled( Double_t *px, Double_t *dummy );
+    static Double_t PtUpsilonCDFscaledPP( Double_t *px, Double_t *dummy );
     static Double_t YUpsilon(Double_t *py, Double_t *dummy);
     static Double_t YUpsilonCDFscaled(Double_t *py, Double_t *dummy);
+    static Double_t YUpsilonCDFscaledPP( Double_t *px, Double_t *dummy );
     static Double_t PtUpsilonPbPb( Double_t *px, Double_t *dummy );
     static Double_t YUpsilonPbPb(Double_t *py, Double_t *dummy);
     static Double_t PtUpsilonPP( Double_t *px, Double_t *dummy );