tune p0 default value

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

diff --git a/EVGEN/AliGenPMDlib.cxx b/EVGEN/AliGenPMDlib.cxx
index fcf3242209f03307593ce30df5f29ebafc885fa9..38b2e6e2604502ed2582fe765455acf931afa81a 100644 (file)
--- a/EVGEN/AliGenPMDlib.cxx
+++ b/EVGEN/AliGenPMDlib.cxx
@@ -13,19 +13,26 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-$Log$
-*/
+/* $Id$ */
+
+//
+// Library of generators for PMD
+// providing y and pt parametrisation
+// for generated tracks
+// Specific for PMD needs
+// Author: PMD Offline Group
+//
+
+#include <TMath.h>
+#include <TPDGCode.h>
 
 #include "AliGenPMDlib.h"
-#include "AliMC.h"
-#include "AliPDG.h"
 
 ClassImp(AliGenPMDlib)
 //
 //  Neutral Pions
 
-Double_t AliGenPMDlib::PtPi0(Double_t *px, Double_t *dummy)
+Double_t AliGenPMDlib::PtPi0(const Double_t *px, const Double_t */*dummy*/)
 {
 //
 //     PT-PARAMETERIZATION CDF, PRL 61(88) 1819
@@ -55,7 +62,7 @@ Double_t AliGenPMDlib::PtPi0(Double_t *px, Double_t *dummy)
 //
 // y-distribution
 //
-Double_t AliGenPMDlib::YPi0( Double_t *py, Double_t *dummy)
+Double_t AliGenPMDlib::YPi0( const Double_t *py, const Double_t */*dummy*/)
 {
   //
   // y parametrisation for pi0
@@ -75,7 +82,7 @@ Double_t AliGenPMDlib::YPi0( Double_t *py, Double_t *dummy)
 
 //                 particle composition
 //
-Int_t AliGenPMDlib::IpPi0()
+Int_t AliGenPMDlib::IpPi0(TRandom *)
 {
 // Pi0
     return kPi0;
@@ -92,7 +99,6 @@ Double_t AliGenPMDlib::PtScal(Double_t pt, Int_t np)
   const Double_t khm[10] = {.13957,.493,.5488,.769,.7826,.958,1.02,0,0,0};
   //     VALUE MESON/PI AT 5 GEV
   const Double_t kfmax[10]={1.,0.3,0.55,1.0,1.0,1.0,1.0,0,0,0};
-  np--;
   Double_t f5=TMath::Power(((sqrt(100.018215)+2.)/(sqrt(100.+khm[np]*khm[np])+2.0)),12.3);
   Double_t fmax2=f5/kfmax[np];
   // PIONS
@@ -107,15 +113,15 @@ Double_t AliGenPMDlib::PtScal(Double_t pt, Int_t np)
 //                pt-distribution
 //____________________________________________________________
 
-Double_t AliGenPMDlib::PtEta( Double_t *px, Double_t *dummy)
+Double_t AliGenPMDlib::PtEta( const Double_t *px, const Double_t */*dummy*/)
 {
 // Kaon pT
-  return PtScal(*px,3);
+  return PtScal(*px, 1);
 }
 
 // y-distribution
 //____________________________________________________________
-Double_t AliGenPMDlib::YEta( Double_t *py, Double_t *dummy)
+Double_t AliGenPMDlib::YEta( const Double_t *py, const Double_t */*dummy*/)
 {
     //
     // y parametrisation for etas
@@ -135,23 +141,23 @@ Double_t AliGenPMDlib::YEta( Double_t *py, Double_t *dummy)
 
 //                 particle composition
 //
-Int_t AliGenPMDlib::IpEta()
+Int_t AliGenPMDlib::IpEta(TRandom *)
 {
     return 221;
 }
 
 
-typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
-GenFunc AliGenPMDlib::GetPt(Param_t param,  const char* tname)
+typedef Double_t (*GenFunc) (const Double_t*,  const Double_t*);
+GenFunc AliGenPMDlib::GetPt(Int_t param,  const char* /*tname*/) const
 {
 // Return pointer to pT parameterisation
     GenFunc func=NULL;
     switch (param) 
     {
-    case Pion:
+    case kPion:
        func=PtPi0;
        break;
-    case Eta:
+    case kEta:
        func=PtEta;
        break;
     default:
@@ -161,16 +167,16 @@ GenFunc AliGenPMDlib::GetPt(Param_t param,  const char* tname)
     return func;
 }
 
-GenFunc AliGenPMDlib::GetY(Param_t param, const char* tname)
+GenFunc AliGenPMDlib::GetY(Int_t param, const char* /*tname*/) const
 {
 // Return pointer to y- parameterisation
     GenFunc func=NULL;
     switch (param) 
     {
-    case Pion:
+    case kPion:
        func=YPi0;
        break;
-    case Eta:
+    case kEta:
        func=YEta;
        break;
     default:
@@ -180,17 +186,17 @@ GenFunc AliGenPMDlib::GetY(Param_t param, const char* tname)
     return func;
 
 }
-typedef Int_t (*GenFuncIp) ();
-GenFuncIp AliGenPMDlib::GetIp(Param_t param,  const char* tname)
+typedef Int_t (*GenFuncIp) (TRandom *);
+GenFuncIp AliGenPMDlib::GetIp(Int_t param,  const char* /*tname*/) const
 {
 // Return pointer to particle type parameterisation
     GenFuncIp func=NULL;
     switch (param) 
     {
-    case Pion:
+    case kPion:
        func=IpPi0;
        break;
-    case Eta:
+    case kEta:
        func=IpEta;
        break;
     default: