Added possibilit to revertex before refiltering

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

diff --git a/EVGEN/AliGenPMDlib.cxx b/EVGEN/AliGenPMDlib.cxx
index 48030a27c22b314c914348ad4d82735bd6664228..38b2e6e2604502ed2582fe765455acf931afa81a 100644 (file)
--- a/EVGEN/AliGenPMDlib.cxx
+++ b/EVGEN/AliGenPMDlib.cxx
@@ -15,6 +15,14 @@
 
 /* $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>
 
@@ -24,7 +32,7 @@ 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
@@ -54,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
@@ -91,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
@@ -106,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
@@ -140,7 +147,7 @@ Int_t AliGenPMDlib::IpEta(TRandom *)
 }
 
 
-typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
+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