D and B distributions updated according to PPR.

diff --git a/EVGEN/AliGenMUONlib.cxx b/EVGEN/AliGenMUONlib.cxx
index 9df3a7188dd9300b1d5fe0fd9055ed822a32a59a..fab7712c3cfa3928f9f68b04a4e1dfa307ae1c56 100644 (file)
--- a/EVGEN/AliGenMUONlib.cxx
+++ b/EVGEN/AliGenMUONlib.cxx
@@ -621,12 +621,24 @@ Int_t AliGenMUONlib::IpEta(TRandom *)
 Double_t AliGenMUONlib::PtCharm( Double_t *px, Double_t */*dummy*/)
 {
 // Charm pT
-  const Double_t kpt0 = 4.08;
-  const Double_t kxn  = 9.40;
+  const Double_t kpt0 = 2.25;
+  const Double_t kxn  = 3.17;
 
   Double_t x=*px;
   //
-  Double_t pass1 = 1.+(x/kpt0);
+  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+  return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtCharmCentral( Double_t *px, Double_t */*dummy*/)
+{
+// Charm pT
+  const Double_t kpt0 = 2.12;
+  const Double_t kxn  = 2.78;
+
+  Double_t x=*px;
+  //
+  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
 //                  y-distribution
@@ -668,8 +680,19 @@ Int_t AliGenMUONlib::IpCharm(TRandom *ran)
 Double_t AliGenMUONlib::PtBeauty( Double_t *px, Double_t */*dummy*/)
 {
 // Beauty pT
-  const Double_t kpt0 = 4.;
-  const Double_t kxn  = 3.6;
+  const Double_t kpt0 = 6.53;
+  const Double_t kxn  = 3.59;
+  Double_t x=*px;
+  //
+  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+  return x/TMath::Power(pass1,kxn);
+}
+
+Double_t AliGenMUONlib::PtBeautyCentral( Double_t *px, Double_t */*dummy*/)
+{
+// Beauty pT
+  const Double_t kpt0 = 6.14;
+  const Double_t kxn  = 2.93;
   Double_t x=*px;
   //
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
@@ -747,10 +770,18 @@ GenFunc AliGenMUONlib::GetPt(Int_t param,  const char* tname) const
        }
        break;  
     case kCharm:
-       func=PtCharm;
+       if (sname == "central") {
+           func=PtCharmCentral;
+       } else {
+           func=PtCharm;
+       }
        break;
     case kBeauty:
-       func=PtBeauty;
+       if (sname == "central") {
+           func=PtBeautyCentral;
+       } else {
+           func=PtBeauty;
+       }
        break;
     case kPion:
        func=PtPion;

diff --git a/EVGEN/AliGenMUONlib.h b/EVGEN/AliGenMUONlib.h
index 50a6efe3fc74fedd789e17c09b8d1b2a59b2a0ee..999572a8ce5c56b3222d0c26a7732d61dc51c183 100644 (file)
--- a/EVGEN/AliGenMUONlib.h
+++ b/EVGEN/AliGenMUONlib.h
@@ -69,11 +69,13 @@ class AliGenMUONlib :
 //
 // Charm    
     static Double_t PtCharm( Double_t *px, Double_t *dummy );
+    static Double_t PtCharmCentral( Double_t *px, Double_t *dummy );
     static Double_t YCharm(Double_t *py, Double_t *dummy);
     static Int_t    IpCharm(TRandom *ran);
 //
 // Beauty
     static Double_t PtBeauty( Double_t *px, Double_t *dummy );
+    static Double_t PtBeautyCentral( Double_t *px, Double_t *dummy );
     static Double_t YBeauty(Double_t *py, Double_t *dummy);
     static Int_t    IpBeauty(TRandom *ran);
 //