updated parametrisation for pi0, eta, omega, phi for different energies, allow select...

[u/mrichter/AliRoot.git] / EVGEN / AliGenPHOSlib.cxx

diff --git a/EVGEN/AliGenPHOSlib.cxx b/EVGEN/AliGenPHOSlib.cxx
index 1585d3151bf2b0c7307293ed6c98c17f9b4e5a4e..9a97f517b2b8fb0107ed697ffa42249d08ed9cc2 100644 (file)
--- a/EVGEN/AliGenPHOSlib.cxx
+++ b/EVGEN/AliGenPHOSlib.cxx
@@ -129,6 +129,14 @@ Int_t AliGenPHOSlib::IpPion(TRandom */*ran*/)
 //    Particle distribution IdPi0Flat  111 (pi0)
 //
 
+Double_t AliGenPHOSlib::PtPi0(const Double_t * px, const Double_t *)
+{
+//     Pion transverse momentum 
+    const Double_t kp0 =1.35; 
+    const Double_t kxn= 6.18;
+    return TMath::Power(kp0 /(kp0 + px[0]), kxn);
+}
+
 Double_t AliGenPHOSlib::PtPi0Flat(const Double_t */*px*/, const Double_t *)
 {
 //     Pion transverse momentum flat distribution 
@@ -556,6 +564,9 @@ GenFunc AliGenPHOSlib::GetPt(Int_t param, const char* /*tname*/) const
       case kPion:     
         func=PtPion;
         break;
+      case kPi0:     
+        func=PtPi0;
+        break;
       case kPi0Flat:     
         func=PtPi0Flat;
         break;
@@ -596,6 +607,7 @@ GenFunc AliGenPHOSlib::GetY(Int_t param, const char* /*tname*/) const
     case kPion:
       func=YPion;
       break;
+    case kPi0:     
     case kPi0Flat:
       func=YPi0Flat;
       break;
@@ -642,6 +654,7 @@ GenFuncIp AliGenPHOSlib::GetIp(Int_t param,  const char* /*tname*/) const
     case kChargedPion:       
       func=IpChargedPion;
       break;
+    case kPi0:     
     case kPi0Flat:       
       func=IpPi0Flat;
       break;