Coding rule violations fixed.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 25 Nov 2008 14:54:18 +0000 (14:54 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 25 Nov 2008 14:54:18 +0000 (14:54 +0000)
15 files changed:
EVGEN/AliGenGSIlib.cxx
EVGEN/AliGenGSIlib.h
EVGEN/AliGenHMPIDlib.cxx
EVGEN/AliGenHMPIDlib.h
EVGEN/AliGenLib.h
EVGEN/AliGenMUONlib.cxx
EVGEN/AliGenMUONlib.h
EVGEN/AliGenPHOSlib.cxx
EVGEN/AliGenPHOSlib.h
EVGEN/AliGenPMDlib.cxx
EVGEN/AliGenPMDlib.h
EVGEN/AliGenParam.cxx
EVGEN/AliGenParam.h
EVGEN/AliGenSTRANGElib.cxx
EVGEN/AliGenSTRANGElib.h

index af3bded..b21fbc1 100644 (file)
@@ -54,7 +54,7 @@ Int_t AliGenGSIlib::IpUpsilon(TRandom *)
   return 553;     
 
 }
-Double_t AliGenGSIlib::PtUpsilonFlat( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtUpsilonFlat( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -72,7 +72,7 @@ Double_t AliGenGSIlib::PtUpsilonFlat( Double_t *px, Double_t */*dummy*/ )
   return weight;
    
 }
-Double_t AliGenGSIlib::YUpsilonFlat(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YUpsilonFlat(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -93,7 +93,7 @@ Double_t AliGenGSIlib::YUpsilonFlat(Double_t *py, Double_t */*dummy*/)
   return yu;
 
 }
-Double_t AliGenGSIlib::PtUpsilonRitman( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtUpsilonRitman( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -110,7 +110,7 @@ Double_t AliGenGSIlib::PtUpsilonRitman( Double_t *px, Double_t */*dummy*/ )
   return x/TMath::Power(pass1,kxn);
    
 }
-Double_t AliGenGSIlib::YUpsilonRitman(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YUpsilonRitman(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -131,7 +131,7 @@ Double_t AliGenGSIlib::YUpsilonRitman(Double_t *py, Double_t */*dummy*/)
   return yu;
    
 }
-Double_t AliGenGSIlib::PtUpsilonKarel( Double_t */*px*/, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtUpsilonKarel( const Double_t */*px*/, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -143,7 +143,7 @@ Double_t AliGenGSIlib::PtUpsilonKarel( Double_t */*px*/, Double_t */*dummy*/ )
   return 0.1;   
 
 }
-Double_t AliGenGSIlib::YUpsilonKarel(Double_t */*py*/, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YUpsilonKarel(const Double_t */*py*/, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -156,7 +156,7 @@ Double_t AliGenGSIlib::YUpsilonKarel(Double_t */*py*/, Double_t */*dummy*/)
   return 0.2;  
 
 }
-Double_t AliGenGSIlib::PtUpsilonMUON( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtUpsilonMUON( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -173,7 +173,7 @@ Double_t AliGenGSIlib::PtUpsilonMUON( Double_t *px, Double_t */*dummy*/ )
   return x/TMath::Power(pass1,kxn);
 
 }
-Double_t AliGenGSIlib::YUpsilonMUON(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YUpsilonMUON(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -209,7 +209,7 @@ Int_t AliGenGSIlib::IpJpsi(TRandom *)
   return 443;     
 
 }
-Double_t AliGenGSIlib::PtJpsiFlat( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtJpsiFlat( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -227,7 +227,7 @@ Double_t AliGenGSIlib::PtJpsiFlat( Double_t *px, Double_t */*dummy*/ )
   return weight;
    
 }
-Double_t AliGenGSIlib::YJpsiFlat(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YJpsiFlat(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -248,7 +248,7 @@ Double_t AliGenGSIlib::YJpsiFlat(Double_t *py, Double_t */*dummy*/)
   return yu;
 
 }
-Double_t AliGenGSIlib::PtJpsiMUON( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtJpsiMUON( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -264,7 +264,7 @@ Double_t AliGenGSIlib::PtJpsiMUON( Double_t *px, Double_t */*dummy*/ )
   return x/TMath::Power(pass1,kxn);
    
 }
-Double_t AliGenGSIlib::PtJpsiRitman( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtJpsiRitman( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -281,7 +281,7 @@ Double_t AliGenGSIlib::PtJpsiRitman( Double_t *px, Double_t */*dummy*/ )
   return x/TMath::Power(pass1,kxn);
    
 }
-Double_t AliGenGSIlib::YJpsiMUON(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YJpsiMUON(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -363,7 +363,7 @@ Int_t AliGenGSIlib::IpCharm(TRandom *ran)
     return ip;
 
 }
-Double_t AliGenGSIlib::PtCharmFlat( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::PtCharmFlat( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -378,7 +378,7 @@ Double_t AliGenGSIlib::PtCharmFlat( Double_t *px, Double_t */*dummy*/)
   return x ;
 
 }
-Double_t AliGenGSIlib::PtCharmGSI( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::PtCharmGSI( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -398,7 +398,7 @@ Double_t AliGenGSIlib::PtCharmGSI( Double_t *px, Double_t */*dummy*/)
   return TMath::Power(x,kp1) * (pass1 + kp3 * pass2);
 
 }
-Double_t AliGenGSIlib::PtCharmMUON( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::PtCharmMUON( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -415,7 +415,7 @@ Double_t AliGenGSIlib::PtCharmMUON( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 
 }
-Double_t AliGenGSIlib::YCharm( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YCharm( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -456,7 +456,7 @@ Int_t AliGenGSIlib::IpBeauty(TRandom *ran)
     
     return ip;
 }
-Double_t AliGenGSIlib::PtBeautyFlat( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::PtBeautyFlat( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -471,7 +471,7 @@ Double_t AliGenGSIlib::PtBeautyFlat( Double_t *px, Double_t */*dummy*/)
   return x ;
 
 }
-Double_t AliGenGSIlib::PtBeautyGSI( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::PtBeautyGSI( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -493,7 +493,7 @@ Double_t AliGenGSIlib::PtBeautyGSI( Double_t *px, Double_t */*dummy*/)
   return TMath::Power(x,kp1) * (pass1 + kp3 * pass2);
 
 }
-Double_t AliGenGSIlib::PtBeautyMUON( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::PtBeautyMUON( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -510,7 +510,7 @@ Double_t AliGenGSIlib::PtBeautyMUON( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 
 }
-Double_t AliGenGSIlib::YBeauty( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YBeauty( const Double_t *px, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -538,7 +538,7 @@ Int_t AliGenGSIlib::IpEta(TRandom *)
   return 221;     
 
 }
-Double_t AliGenGSIlib::PtEtaPHOS( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtEtaPHOS( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -549,7 +549,7 @@ Double_t AliGenGSIlib::PtEtaPHOS( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,3);  //  3==> Eta in the PtScal function
    
 }
-Double_t AliGenGSIlib::YEtaPHOS(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YEtaPHOS(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -582,7 +582,7 @@ Int_t AliGenGSIlib::IpEtaprime(TRandom *)
   return 331;     
 
 }
-Double_t AliGenGSIlib::PtEtaprimePHOS( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtEtaprimePHOS( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -593,7 +593,7 @@ Double_t AliGenGSIlib::PtEtaprimePHOS( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,5);  //  5==> Etaprime in the PtScal function
    
 }
-Double_t AliGenGSIlib::YEtaprimePHOS(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YEtaprimePHOS(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -626,7 +626,7 @@ Int_t AliGenGSIlib::IpOmega(TRandom *)
   return 223;     
 
 }
-Double_t AliGenGSIlib::PtOmega( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtOmega( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -637,7 +637,7 @@ Double_t AliGenGSIlib::PtOmega( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,4);  //  4==> Omega in the PtScal function
    
 }
-Double_t AliGenGSIlib::YOmega(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YOmega(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -672,7 +672,7 @@ Int_t AliGenGSIlib::IpRho(TRandom *)
   return 113;     
 
 }
-Double_t AliGenGSIlib::PtRho( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtRho( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -683,7 +683,7 @@ Double_t AliGenGSIlib::PtRho( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,11);  //  11==> Rho in the PtScal function
    
 }
-Double_t AliGenGSIlib::YRho(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YRho(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -732,7 +732,7 @@ Int_t AliGenGSIlib::IpPionPHOS(TRandom *ran)
       }
     }
 }
-Double_t AliGenGSIlib::PtPion( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtPion( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -766,7 +766,7 @@ Double_t AliGenGSIlib::PtPion( Double_t *px, Double_t */*dummy*/ )
   return y*x;
    
 }
-Double_t AliGenGSIlib::YPion(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YPion(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -815,7 +815,7 @@ Int_t AliGenGSIlib::IpKaonPHOS(TRandom *ran)
       }
     }
 }
-Double_t AliGenGSIlib::PtKaonPHOS( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtKaonPHOS( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -826,7 +826,7 @@ Double_t AliGenGSIlib::PtKaonPHOS( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,2);  //  2==> Kaon in the PtScal function
    
 }
-Double_t AliGenGSIlib::YKaonPHOS(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YKaonPHOS(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -859,7 +859,7 @@ Int_t AliGenGSIlib::IpPhi(TRandom *)
   return 333;     
 
 }
-Double_t AliGenGSIlib::PtPhiPHOS( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtPhiPHOS( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -870,7 +870,7 @@ Double_t AliGenGSIlib::PtPhiPHOS( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,6);  //  6==> Phi in the PtScal function
    
 }
-Double_t AliGenGSIlib::YPhiPHOS(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YPhiPHOS(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -920,7 +920,7 @@ Int_t AliGenGSIlib::IpBaryons(TRandom *ran)
       }
     }
 }
-Double_t AliGenGSIlib::PtBaryons( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenGSIlib::PtBaryons( const Double_t *px, const Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -931,7 +931,7 @@ Double_t AliGenGSIlib::PtBaryons( Double_t *px, Double_t */*dummy*/ )
   return PtScal(*px,7);  //  7==> Baryon in the PtScal function
    
 }
-Double_t AliGenGSIlib::YBaryons(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenGSIlib::YBaryons(const Double_t *py, const Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -995,7 +995,7 @@ Double_t AliGenGSIlib::YBaryons(Double_t *py, Double_t */*dummy*/)
 //    
 //==========================================================================
 
-typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
+typedef Double_t (*GenFunc) (const Double_t*,  const Double_t*);
 
 typedef Int_t (*GenFuncIp) (TRandom *);
 
index 1aa81e1..f761403 100644 (file)
@@ -33,29 +33,29 @@ class AliGenGSIlib :public AliGenLib {
 // Upsilon
     static Int_t    IpUpsilon(TRandom *ran);
 // Upsilon RITMAN   
-    static Double_t PtUpsilonRitman( Double_t *px, Double_t *dummy );
-    static Double_t YUpsilonRitman(Double_t *py, Double_t *dummy);
+    static Double_t PtUpsilonRitman( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonRitman(const Double_t *py, const Double_t *dummy);
 // Upsilon FLAT   
-    static Double_t PtUpsilonFlat( Double_t *px, Double_t *dummy );
-    static Double_t YUpsilonFlat(Double_t *py, Double_t *dummy);
+    static Double_t PtUpsilonFlat( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonFlat(const Double_t *py, const Double_t *dummy);
 // Upsilon Karel
-    static Double_t PtUpsilonKarel( Double_t *px, Double_t *dummy );
-    static Double_t YUpsilonKarel(Double_t *py, Double_t *dummy);
+    static Double_t PtUpsilonKarel( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonKarel(const Double_t *py, const Double_t *dummy);
 // Upsilon MUONlib
-    static Double_t PtUpsilonMUON( Double_t *px, Double_t *dummy );
-    static Double_t YUpsilonMUON(Double_t *py, Double_t *dummy);
+    static Double_t PtUpsilonMUON( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonMUON(const Double_t *py, const Double_t *dummy);
 
 
 // JPsi 
     static Int_t    IpJpsi(TRandom *ran);
 // JPsi FLAT   
-    static Double_t PtJpsiFlat( Double_t *px, Double_t *dummy );
-    static Double_t YJpsiFlat(Double_t *py, Double_t *dummy);
+    static Double_t PtJpsiFlat( const Double_t *px, const Double_t *dummy );
+    static Double_t YJpsiFlat(const Double_t *py, const Double_t *dummy);
 // JPsi from MUONlib
-    static Double_t PtJpsiMUON( Double_t *px, Double_t *dummy );
-    static Double_t YJpsiMUON(Double_t *py, Double_t *dummy);
+    static Double_t PtJpsiMUON( const Double_t *px, const Double_t *dummy );
+    static Double_t YJpsiMUON(const Double_t *py, const Double_t *dummy);
 // JPsi from Ritman
-    static Double_t PtJpsiRitman( Double_t *px, Double_t *dummy );
+    static Double_t PtJpsiRitman( const Double_t *px, const Double_t *dummy );
 
     // JPsi from Sergei
     //    static Double_t PtJpsi( Double_t *px, Double_t *dummy );
@@ -65,61 +65,61 @@ class AliGenGSIlib :public AliGenLib {
 
 // Charm 
     static Int_t IpCharm(TRandom *ran);
-    static Double_t PtCharmFlat( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmMUON( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmGSI( Double_t *px, Double_t *dummy );
-    static Double_t YCharm(Double_t *py, Double_t *dummy);
+    static Double_t PtCharmFlat( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmMUON( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmGSI( const Double_t *px, const Double_t *dummy );
+    static Double_t YCharm(const Double_t *py, const Double_t *dummy);
 
 
 // Beauty
     static Int_t IpBeauty(TRandom *ran);
-    static Double_t PtBeautyFlat( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyMUON( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyGSI( Double_t *px, Double_t *dummy );
-    static Double_t YBeauty(Double_t *py, Double_t *dummy);
+    static Double_t PtBeautyFlat( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyMUON( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyGSI( const Double_t *px, const Double_t *dummy );
+    static Double_t YBeauty(const Double_t *py, const Double_t *dummy);
 
 
 // Eta
     static Int_t IpEta(TRandom *ran);
-    static Double_t PtEtaPHOS( Double_t *px, Double_t *dummy );
-    static Double_t YEtaPHOS(Double_t *py, Double_t *dummy);
+    static Double_t PtEtaPHOS( const Double_t *px, const Double_t *dummy );
+    static Double_t YEtaPHOS(const Double_t *py, const Double_t *dummy);
 
 
 // Etaprime
     static Int_t IpEtaprime(TRandom *ran);
-    static Double_t PtEtaprimePHOS( Double_t *px, Double_t *dummy );
-    static Double_t YEtaprimePHOS(Double_t *py, Double_t *dummy);
+    static Double_t PtEtaprimePHOS( const Double_t *px, const Double_t *dummy );
+    static Double_t YEtaprimePHOS(const Double_t *py, const Double_t *dummy);
 
 
 // Omega
     static Int_t IpOmega(TRandom *ran);
-    static Double_t PtOmega( Double_t *px, Double_t *dummy );
-    static Double_t YOmega(Double_t *py, Double_t *dummy);
+    static Double_t PtOmega( const Double_t *px, const Double_t *dummy );
+    static Double_t YOmega(const Double_t *py, const Double_t *dummy);
 
 
 // Rho
    static Int_t IpRho(TRandom *ran);
-   static Double_t PtRho( Double_t *px, Double_t *dummy );
-   static Double_t YRho(Double_t *py, Double_t *dummy);
+   static Double_t PtRho( const Double_t *px, const Double_t *dummy );
+   static Double_t YRho(const Double_t *py, const Double_t *dummy);
 
 
 
 // Kaon
     static Int_t IpKaonPHOS(TRandom *ran);
-    static Double_t PtKaonPHOS( Double_t *px, Double_t *dummy );
-    static Double_t YKaonPHOS(Double_t *py, Double_t *dummy);
+    static Double_t PtKaonPHOS( const Double_t *px, const Double_t *dummy );
+    static Double_t YKaonPHOS(const Double_t *py, const Double_t *dummy);
 
 
 // Pion
     static Int_t IpPionPHOS(TRandom *ran);
-    static Double_t PtPion( Double_t *px, Double_t *dummy );
-    static Double_t YPion(Double_t *py, Double_t *dummy);
+    static Double_t PtPion( const Double_t *px, const Double_t *dummy );
+    static Double_t YPion(const Double_t *py, const Double_t *dummy);
 
 
 // Phi
     static Int_t IpPhi(TRandom *ran);
-    static Double_t PtPhiPHOS( Double_t *px, Double_t *dummy );
-    static Double_t YPhiPHOS(Double_t *py, Double_t *dummy);
+    static Double_t PtPhiPHOS( const Double_t *px, const Double_t *dummy );
+    static Double_t YPhiPHOS(const Double_t *py, const Double_t *dummy);
 
 
 // Lambda
@@ -130,8 +130,8 @@ class AliGenGSIlib :public AliGenLib {
 
 // Baryons
     static Int_t IpBaryons(TRandom *ran);
-    static Double_t PtBaryons( Double_t *px, Double_t *dummy );
-    static Double_t YBaryons(Double_t *py, Double_t *dummy);
+    static Double_t PtBaryons( const Double_t *px, const Double_t *dummy );
+    static Double_t YBaryons(const Double_t *py, const Double_t *dummy);
 
 
 
index 67a1865..11526cf 100644 (file)
@@ -44,19 +44,19 @@ Int_t AliGenHMPIDlib::IpPiPlus(TRandom *)
   return 211;     
 }
 
-Double_t AliGenHMPIDlib::PtPiPlusFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtPiPlusFlat( const Double_t *, const Double_t *)
 {
   //PiPlus FLAT pt-distribution
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtPiPlusExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtPiPlusExp( const Double_t *x, const Double_t *)
 {
   //PiPlus     EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YPiPlusFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YPiPlusFlat(const Double_t *,const Double_t *)
 {
   //PiPlus            y-distribution 
  return 1;
@@ -72,19 +72,19 @@ Int_t AliGenHMPIDlib::IpPiMinus(TRandom *)
   return -211;     
 }
 
-Double_t AliGenHMPIDlib::PtPiMinusFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtPiMinusFlat( const Double_t *, const Double_t *)
 {
 // PiMinus FLAT pt-distribution
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtPiMinusExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtPiMinusExp( const Double_t *x, const Double_t *)
 {
 //PiMinus     EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YPiMinusFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YPiMinusFlat(const Double_t *,const Double_t *)
 {
 //PiMinus          y-distribution 
   return 1;
@@ -100,19 +100,19 @@ Int_t AliGenHMPIDlib::IpKPlus(TRandom *)
   return 321;
 }
 
-Double_t AliGenHMPIDlib::PtKPlusFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtKPlusFlat( const Double_t *, const Double_t *)
 {
 // K+ FLAT pt-distribution
   return 1;
 }
 
-Double_t AliGenHMPIDlib::PtKPlusExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtKPlusExp( const Double_t *x, const Double_t *)
 {
 // K+   EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);
 }
 
-Double_t AliGenHMPIDlib::YKPlusFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YKPlusFlat(const Double_t *,const Double_t *)
 {
 // K+             y-distribution
   return 1;
@@ -128,19 +128,19 @@ Int_t AliGenHMPIDlib::IpKMinus(TRandom *)
   return -321;
 }
 
-Double_t AliGenHMPIDlib::PtKMinusFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtKMinusFlat( const Double_t *, const Double_t *)
 {
 // K- FLAT pt-distribution
   return 1;
 }
 
-Double_t AliGenHMPIDlib::PtKMinusExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtKMinusExp( const Double_t *x, const Double_t *)
 {
 // K-   EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);
 }
 
-Double_t AliGenHMPIDlib::YKMinusFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YKMinusFlat(const Double_t *,const Double_t *)
 {
 // K-             y-distribution
   return 1;
@@ -156,19 +156,19 @@ Int_t AliGenHMPIDlib::IpK0s(TRandom *)
   return 310;
 }
 
-Double_t AliGenHMPIDlib::PtK0sFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtK0sFlat( const Double_t *, const Double_t *)
 {
 // K0s FLAT pt-distribution
   return 1;
 }
 
-Double_t AliGenHMPIDlib::PtK0sExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtK0sExp( const Double_t *x, const Double_t *)
 {
 // K0s   EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);
 }
 
-Double_t AliGenHMPIDlib::YK0sFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YK0sFlat(const Double_t *,const Double_t *)
 {
 // K0s             y-distribution
   return 1;
@@ -184,19 +184,19 @@ Int_t AliGenHMPIDlib::IpPhi(TRandom *)
   return 333;     
 }
 
-Double_t AliGenHMPIDlib::PtPhiFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtPhiFlat( const Double_t *, const Double_t *)
 {
 // Phi FLAT pt-distribution
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtPhiExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtPhiExp( const Double_t *x, const Double_t *)
 {
 //phi     EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YPhiFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YPhiFlat(const Double_t *,const Double_t *)
 {
 //phi             y-distribution 
   return 1;
@@ -212,20 +212,20 @@ Int_t AliGenHMPIDlib::IpProton(TRandom *)
   return 2122;     
 }
 
-Double_t AliGenHMPIDlib::PtProtonFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtProtonFlat( const Double_t *, const Double_t *)
 {
 // ProtonFLAT pt-distribution
 
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtProtonExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtProtonExp( const Double_t *x, const Double_t *)
 {
 //Proton    EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YProtonFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YProtonFlat(const Double_t *,const Double_t *)
 {
   //Proton            y-distribution 
   return 1;
@@ -241,20 +241,20 @@ Int_t AliGenHMPIDlib::IpProtonBar(TRandom *)
   return -2122;     
 }
 
-Double_t AliGenHMPIDlib::PtProtonBarFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtProtonBarFlat( const Double_t *, const Double_t *)
 {
 // ProtonBar FLAT pt-distribution
 
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtProtonBarExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtProtonBarExp( const Double_t *x, const Double_t *)
 {
 //ProtonBar     EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YProtonBarFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YProtonBarFlat(const Double_t *,const Double_t *)
 {
   //ProtonBar             y-distribution 
   return 1;
@@ -270,20 +270,20 @@ Int_t AliGenHMPIDlib::IpLambda(TRandom *)
   return 3122;     
 }
 
-Double_t AliGenHMPIDlib::PtLambdaFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtLambdaFlat( const Double_t *, const Double_t *)
 {
 // Lambda FLAT pt-distribution
 
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtLambdaExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtLambdaExp( const Double_t *x, const Double_t *)
 {
 //Lambda     EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YLambdaFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YLambdaFlat(const Double_t *,const Double_t *)
 {
   //Lambda             y-distribution 
   return 1;
@@ -299,20 +299,20 @@ Int_t AliGenHMPIDlib::IpLambdaBar(TRandom *)
   return -3122;     
 }
 
-Double_t AliGenHMPIDlib::PtLambdaBarFlat( Double_t *, Double_t *)
+Double_t AliGenHMPIDlib::PtLambdaBarFlat( const Double_t *, const Double_t *)
 {
 // LambdaBar FLAT pt-distribution
 
   return 1; 
 }
 
-Double_t AliGenHMPIDlib::PtLambdaBarExp( Double_t *x, Double_t *)
+Double_t AliGenHMPIDlib::PtLambdaBarExp( const Double_t *x, const Double_t *)
 {
 //LambdaBar     EXP  pt-distribution
   return x[0]*TMath::Exp(-x[0]/0.17);   
 }
 
-Double_t AliGenHMPIDlib::YLambdaBarFlat(Double_t *,Double_t *)
+Double_t AliGenHMPIDlib::YLambdaBarFlat(const Double_t *,const Double_t *)
 {
   //LambdaBar             y-distribution 
   return 1;
@@ -322,7 +322,7 @@ Double_t AliGenHMPIDlib::YLambdaBarFlat(Double_t *,Double_t *)
 
 
 
-typedef Double_t (*GenFunc)   (Double_t*,  Double_t*);
+typedef Double_t (*GenFunc)   (const Double_t*,  const Double_t*);
 typedef Int_t    (*GenFuncIp) (TRandom *);
 
 GenFunc AliGenHMPIDlib::GetPt(Int_t iPID, const char * sForm) const
index efcc4a7..3894379 100644 (file)
@@ -30,63 +30,63 @@ class AliGenHMPIDlib :public AliGenLib {
 
 //Pi+
   static Int_t IpPiPlus(TRandom *ran);
-  static Double_t PtPiPlusFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtPiPlusExp (Double_t *px, Double_t *dummy);
-  static Double_t YPiPlusFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtPiPlusFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtPiPlusExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YPiPlusFlat (const Double_t *py, const Double_t *dummy);
 
 //Pi-
   static Int_t IpPiMinus(TRandom *ran);
-  static Double_t PtPiMinusFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtPiMinusExp (Double_t *px, Double_t *dummy);
-  static Double_t YPiMinusFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtPiMinusFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtPiMinusExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YPiMinusFlat (const Double_t *py, const Double_t *dummy);
 
 //K+
   static Int_t IpKPlus(TRandom *ran);
-  static Double_t PtKPlusFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtKPlusExp (Double_t *px, Double_t *dummy);
-  static Double_t YKPlusFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtKPlusFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtKPlusExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YKPlusFlat (const Double_t *py, const Double_t *dummy);
 
 //K-
   static Int_t IpKMinus(TRandom *ran);
-  static Double_t PtKMinusFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtKMinusExp (Double_t *px, Double_t *dummy);
-  static Double_t YKMinusFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtKMinusFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtKMinusExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YKMinusFlat (const Double_t *py, const Double_t *dummy);
 
 // K0_s
   static Int_t IpK0s(TRandom *ran);
-  static Double_t PtK0sFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtK0sExp (Double_t *px, Double_t *dummy);
-  static Double_t YK0sFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtK0sFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtK0sExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YK0sFlat (const Double_t *py, const Double_t *dummy);
 
 // Phi(1020)
   static Int_t IpPhi(TRandom *ran);
-  static Double_t PtPhiFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtPhiExp (Double_t *px, Double_t *dummy);
-  static Double_t YPhiFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtPhiFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtPhiExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YPhiFlat (const Double_t *py, const Double_t *dummy);
 
 //Proton
   static Int_t IpProton(TRandom *ran);
-  static Double_t PtProtonFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtProtonExp (Double_t *px, Double_t *dummy);
-  static Double_t YProtonFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtProtonFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtProtonExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YProtonFlat (const Double_t *py, const Double_t *dummy);
 
 //ProtonBar
   static Int_t IpProtonBar(TRandom *ran);
-  static Double_t PtProtonBarFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtProtonBarExp (Double_t *px, Double_t *dummy);
-  static Double_t YProtonBarFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtProtonBarFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtProtonBarExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YProtonBarFlat (const Double_t *py, const Double_t *dummy);
 
 // Lambda
   static Int_t IpLambda(TRandom *ran);
-  static Double_t PtLambdaFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtLambdaExp (Double_t *px, Double_t *dummy);
-  static Double_t YLambdaFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtLambdaFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtLambdaExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YLambdaFlat (const Double_t *py, const Double_t *dummy);
 
 // LambdaBar
   static Int_t IpLambdaBar(TRandom *ran);
-  static Double_t PtLambdaBarFlat(Double_t *px, Double_t *dummy);
-  static Double_t PtLambdaBarExp (Double_t *px, Double_t *dummy);
-  static Double_t YLambdaBarFlat (Double_t *py, Double_t *dummy);
+  static Double_t PtLambdaBarFlat(const Double_t *px, const Double_t *dummy);
+  static Double_t PtLambdaBarExp (const Double_t *px, const Double_t *dummy);
+  static Double_t YLambdaBarFlat (const Double_t *py, const Double_t *dummy);
 
   ClassDef(AliGenHMPIDlib,0)
 };
index ff51cd9..306ab15 100644 (file)
@@ -15,7 +15,7 @@ class AliGenLib :
  public:
 //
     virtual ~AliGenLib(){}
-    typedef Double_t (*GenFunc)  (Double_t *, Double_t *);
+    typedef Double_t (*GenFunc)  (const Double_t *, const Double_t *);
     typedef Int_t    (*GenFuncIp)(TRandom *);    
     virtual GenFunc   GetPt(Int_t param, const char *tname) const   = 0;
     virtual GenFunc   GetY (Int_t param, const char *tname) const   = 0;
index c1703d7..7d03131 100644 (file)
@@ -32,7 +32,7 @@
 ClassImp(AliGenMUONlib)
 //
 //  Pions
-Double_t AliGenMUONlib::PtPion(Double_t *px, Double_t* /*dummy*/)
+Double_t AliGenMUONlib::PtPion(const Double_t *px, const Double_t* /*dummy*/)
 {
 //
 //     PT-PARAMETERIZATION CDF, PRL 61(88) 1819
@@ -61,7 +61,7 @@ Double_t AliGenMUONlib::PtPion(Double_t *px, Double_t* /*dummy*/)
 //
 // y-distribution
 //
-Double_t AliGenMUONlib::YPion( Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YPion( const Double_t *py, const Double_t */*dummy*/)
 {
 // Pion y
   Double_t y=TMath::Abs(*py);
@@ -111,7 +111,7 @@ Double_t AliGenMUONlib::PtScal(Double_t pt, Int_t np)
 //
 //                pt-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::PtKaon( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtKaon( const Double_t *px, const Double_t */*dummy*/)
 {
 // Kaon pT
   return PtScal(*px,2);
@@ -119,7 +119,7 @@ Double_t AliGenMUONlib::PtKaon( Double_t *px, Double_t */*dummy*/)
 
 // y-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::YKaon( Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YKaon( const Double_t *py, const Double_t */*dummy*/)
 {
 // Kaon y
   Double_t y=TMath::Abs(*py);
@@ -151,7 +151,7 @@ Int_t AliGenMUONlib::IpKaon(TRandom *ran)
 //
 //                pt-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::PtJpsi( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsi( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT
   const Double_t kpt0 = 4.;
@@ -162,7 +162,7 @@ Double_t AliGenMUONlib::PtJpsi( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtJpsiCDFscaled( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaled( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT
 //
@@ -178,7 +178,7 @@ Double_t AliGenMUONlib::PtJpsiCDFscaled( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtJpsiCDFscaledold( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledold( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT
 //
@@ -193,7 +193,7 @@ Double_t AliGenMUONlib::PtJpsiCDFscaledold( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtJpsiCDFscaledPP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT
 //
@@ -208,7 +208,7 @@ Double_t AliGenMUONlib::PtJpsiCDFscaledPP( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtJpsiCDFscaledPPold( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPPold( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT
 //
@@ -223,7 +223,7 @@ Double_t AliGenMUONlib::PtJpsiCDFscaledPPold( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtJpsiCDFscaledPP10( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT
 //
@@ -238,12 +238,12 @@ Double_t AliGenMUONlib::PtJpsiCDFscaledPP10( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtJpsiFlat( Double_t */*px*/, Double_t */*dummy*/ )
+Double_t AliGenMUONlib::PtJpsiFlat( const Double_t */*px*/, const Double_t */*dummy*/ )
 {
   return 1.;
 }
 
-Double_t AliGenMUONlib::PtJpsiPbPb( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiPbPb( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT spectrum
 //
@@ -270,7 +270,7 @@ Double_t AliGenMUONlib::PtJpsiPbPb( Double_t *px, Double_t */*dummy*/)
     return y;
 }
 
-Double_t AliGenMUONlib::PtJpsiBPbPb( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiBPbPb( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT spectrum
 // B -> J/Psi X
@@ -284,7 +284,7 @@ Double_t AliGenMUONlib::PtJpsiBPbPb( Double_t *px, Double_t */*dummy*/)
 }
 
 
-Double_t AliGenMUONlib::PtJpsiPP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtJpsiPP( const Double_t *px, const Double_t */*dummy*/)
 {
 // J/Psi pT spectrum
 //
@@ -311,7 +311,7 @@ Double_t AliGenMUONlib::PtJpsiPP( Double_t *px, Double_t */*dummy*/)
 //
 //               y-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::YJpsi(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YJpsi(const Double_t *py, const Double_t */*dummy*/)
 {
 // J/psi y
   const Double_t ky0 = 4.;
@@ -326,13 +326,13 @@ Double_t AliGenMUONlib::YJpsi(Double_t *py, Double_t */*dummy*/)
   return yj;
 }
 
-Double_t AliGenMUONlib::YJpsiFlat( Double_t */*py*/, Double_t */*dummy*/ )
+Double_t AliGenMUONlib::YJpsiFlat( const Double_t */*py*/, const Double_t */*dummy*/ )
 {
   return 1.;
 }
 
 
-Double_t AliGenMUONlib::YJpsiPbPb( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YJpsiPbPb( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -361,31 +361,31 @@ Double_t AliGenMUONlib::YJpsiPbPb( Double_t *px, Double_t */*dummy*/)
     return y;
 }
 
-Double_t AliGenMUONlib::YJpsiCDFscaled( Double_t *px, Double_t* dummy)
+Double_t AliGenMUONlib::YJpsiCDFscaled( const Double_t *px, const Double_t* dummy)
 {
     // J/Psi y 
     return AliGenMUONlib::YJpsiPbPb(px, dummy);
 }
 
-Double_t AliGenMUONlib::YJpsiCDFscaledold( Double_t *px, Double_t* dummy)
+Double_t AliGenMUONlib::YJpsiCDFscaledold( const Double_t *px, const Double_t* dummy)
 {
     // J/Psi y 
     return AliGenMUONlib::YJpsiPbPb(px, dummy);
 }
 
-Double_t AliGenMUONlib::YJpsiCDFscaledPP( Double_t *px, Double_t* dummy)
+Double_t AliGenMUONlib::YJpsiCDFscaledPP( const Double_t *px, const Double_t* dummy)
 {
     // J/Psi y 
     return AliGenMUONlib::YJpsiPP(px, dummy);
 }
 
-Double_t AliGenMUONlib::YJpsiCDFscaledPPold( Double_t *px, Double_t* dummy)
+Double_t AliGenMUONlib::YJpsiCDFscaledPPold( const Double_t *px, const Double_t* dummy)
 {
     // J/Psi y 
     return AliGenMUONlib::YJpsiPP(px, dummy);
 }
 
-Double_t AliGenMUONlib::YJpsiCDFscaledPP10( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YJpsiCDFscaledPP10( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -416,7 +416,7 @@ Double_t AliGenMUONlib::YJpsiCDFscaledPP10( Double_t *px, Double_t */*dummy*/)
     return y;
 }
 
-Double_t AliGenMUONlib::YJpsiPP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YJpsiPP( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -446,7 +446,7 @@ Double_t AliGenMUONlib::YJpsiPP( Double_t *px, Double_t */*dummy*/)
     return y;
 }
 
-Double_t AliGenMUONlib::YJpsiBPbPb( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YJpsiBPbPb( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -505,7 +505,7 @@ Int_t AliGenMUONlib::IpJpsiFamily(TRandom *)
 //
 //                  pt-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::PtUpsilon( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenMUONlib::PtUpsilon( const Double_t *px, const Double_t */*dummy*/ )
 {
 // Upsilon pT
   const Double_t kpt0 = 5.3;
@@ -516,7 +516,7 @@ Double_t AliGenMUONlib::PtUpsilon( Double_t *px, Double_t */*dummy*/ )
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtUpsilonCDFscaled( Double_t *px, Double_t */*dummy*/ )
+Double_t AliGenMUONlib::PtUpsilonCDFscaled( const Double_t *px, const Double_t */*dummy*/ )
 {
 // Upsilon pT
   const Double_t kpt0 = 7.753;
@@ -527,7 +527,7 @@ 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*/ )
+Double_t AliGenMUONlib::PtUpsilonCDFscaledPP( const Double_t *px, const Double_t */*dummy*/ )
 {
 // Upsilon pT
 //
@@ -543,12 +543,12 @@ Double_t AliGenMUONlib::PtUpsilonCDFscaledPP( Double_t *px, Double_t */*dummy*/
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtUpsilonFlat( Double_t */*px*/, Double_t */*dummy*/ )
+Double_t AliGenMUONlib::PtUpsilonFlat( const Double_t */*px*/, const Double_t */*dummy*/ )
 {
   return 1.;
 }
 
-Double_t AliGenMUONlib::PtUpsilonPbPb( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtUpsilonPbPb( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -577,7 +577,7 @@ Double_t AliGenMUONlib::PtUpsilonPbPb( Double_t *px, Double_t */*dummy*/)
     return y;
 }
 
-Double_t AliGenMUONlib::PtUpsilonPP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtUpsilonPP( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -609,7 +609,7 @@ Double_t AliGenMUONlib::PtUpsilonPP( Double_t *px, Double_t */*dummy*/)
 //                    y-distribution
 //
 //____________________________________________________________
-Double_t AliGenMUONlib::YUpsilon(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YUpsilon(const Double_t *py, const Double_t */*dummy*/)
 {
 // Upsilon y
   const Double_t ky0 = 3.;
@@ -625,7 +625,7 @@ Double_t AliGenMUONlib::YUpsilon(Double_t *py, Double_t */*dummy*/)
 }
 
 
-Double_t AliGenMUONlib::YUpsilonPbPb( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YUpsilonPbPb( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -649,28 +649,28 @@ Double_t AliGenMUONlib::YUpsilonPbPb( Double_t *px, Double_t */*dummy*/)
     return y;
 }
 
-Double_t AliGenMUONlib::YUpsilonCDFscaled( Double_t *px, Double_t *dummy)
+Double_t AliGenMUONlib::YUpsilonCDFscaled( const Double_t *px, const Double_t *dummy)
 {
     // Upsilon y
     return AliGenMUONlib::YUpsilonPbPb(px, dummy);
     
 }
 
-Double_t AliGenMUONlib::YUpsilonCDFscaledPP( Double_t *px, Double_t *dummy)
+Double_t AliGenMUONlib::YUpsilonCDFscaledPP( const Double_t *px, const Double_t *dummy)
 {
     // Upsilon y
     return AliGenMUONlib::YUpsilonPP(px, dummy);
     
 }
 
-Double_t AliGenMUONlib::YUpsilonFlat( Double_t */*px*/, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YUpsilonFlat( const Double_t */*px*/, const Double_t */*dummy*/)
 {
     // Upsilon y
     return 1.;
     
 }
 
-Double_t AliGenMUONlib::YUpsilonPP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YUpsilonPP( const Double_t *px, const Double_t */*dummy*/)
 {
 
 //
@@ -733,13 +733,13 @@ Int_t AliGenMUONlib::IpUpsilonFamily(TRandom *)
 //
 //    pt-distribution (by scaling of pion distribution)
 //____________________________________________________________
-Double_t AliGenMUONlib::PtPhi( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtPhi( const Double_t *px, const Double_t */*dummy*/)
 {
 // Phi pT
   return PtScal(*px,7);
 }
 //    y-distribution
-Double_t AliGenMUONlib::YPhi( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YPhi( const Double_t *px, const Double_t */*dummy*/)
 {
 // Phi y
     Double_t *dum=0;
@@ -759,13 +759,13 @@ Int_t AliGenMUONlib::IpPhi(TRandom *)
 //
 //    pt-distribution (by scaling of pion distribution)
 //____________________________________________________________
-Double_t AliGenMUONlib::PtOmega( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtOmega( const Double_t *px, const Double_t */*dummy*/)
 {
 // Omega pT
   return PtScal(*px,5);
 }
 //    y-distribution
-Double_t AliGenMUONlib::YOmega( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YOmega( const Double_t *px, const Double_t */*dummy*/)
 {
 // Omega y
     Double_t *dum=0;
@@ -786,13 +786,13 @@ Int_t AliGenMUONlib::IpOmega(TRandom *)
 //
 //    pt-distribution (by scaling of pion distribution)
 //____________________________________________________________
-Double_t AliGenMUONlib::PtEta( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtEta( const Double_t *px, const Double_t */*dummy*/)
 {
 // Eta pT
   return PtScal(*px,3);
 }
 //    y-distribution
-Double_t AliGenMUONlib::YEta( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YEta( const Double_t *px, const Double_t */*dummy*/)
 {
 // Eta y
     Double_t *dum=0;
@@ -812,7 +812,7 @@ Int_t AliGenMUONlib::IpEta(TRandom *)
 //
 //                    pt-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::PtCharm( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharm( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
   const Double_t kpt0 = 2.25;
@@ -823,7 +823,7 @@ Double_t AliGenMUONlib::PtCharm( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtCharmCentral( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmCentral( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
   const Double_t kpt0 = 2.12;
@@ -833,7 +833,7 @@ Double_t AliGenMUONlib::PtCharmCentral( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // FiMjSkPP define theoretical uncertainties around F0M0S0PP as follows:
 // PtCharmFiMjSkPP = PtCharmF0M0S0PP * (dN(i,j,k)/dpt / dN(0,0,0)/dpt)_MNR
@@ -864,7 +864,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF1M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF1M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 6.06 mb
@@ -876,7 +876,7 @@ Double_t AliGenMUONlib::PtCharmF1M0S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF2M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF2M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 6.06 mb
@@ -888,7 +888,7 @@ Double_t AliGenMUONlib::PtCharmF2M0S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M1S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M1S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 7.69 mb
@@ -900,7 +900,7 @@ Double_t AliGenMUONlib::PtCharmF0M1S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M2S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M2S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 4.81 mb
@@ -912,7 +912,7 @@ Double_t AliGenMUONlib::PtCharmF0M2S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S1PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S1PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 14.09 mb
@@ -924,7 +924,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S1PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S2PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S2PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 1.52 mb
@@ -936,7 +936,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S2PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S3PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S3PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 3.67 mb
@@ -948,7 +948,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S3PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S4PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S4PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 3.38 mb
@@ -960,7 +960,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S4PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S5PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S5PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 10.37 mb
@@ -972,7 +972,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S5PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn)*(1.+c[0]*x+c[1]*x*x)/(1.+c[2]*x*x);
 }
-Double_t AliGenMUONlib::PtCharmF0M0S6PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtCharmF0M0S6PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm pT
 // Corresponding NLO total cross section is 7.22 mb
@@ -986,7 +986,7 @@ Double_t AliGenMUONlib::PtCharmF0M0S6PP( Double_t *px, Double_t */*dummy*/)
 }
 
 //                  y-distribution
-Double_t AliGenMUONlib::YCharm( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharm( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y :: Carrer & Dainese : ALICE-INT-2003-019 v.3 (hep-ph/0311225) 
 // Pythia tuned to reproduce the distribution given by the HVQMNR program based on NLO calculations (pQCD)
@@ -1006,7 +1006,7 @@ Double_t AliGenMUONlib::YCharm( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // FiMjSkPP define theoretical uncertainties around F0M0S0PP as follows:
 // YCharmFiMjSkPP = YCharmF0M0S0PP * (dN(i,j,k)/dy / dN(0,0,0)/dy)_MNR
@@ -1042,7 +1042,7 @@ Double_t AliGenMUONlib::YCharmF0M0S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF1M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF1M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 6.06 mb
@@ -1060,7 +1060,7 @@ Double_t AliGenMUONlib::YCharmF1M0S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF2M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF2M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 6.06 mb
@@ -1078,7 +1078,7 @@ Double_t AliGenMUONlib::YCharmF2M0S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M1S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M1S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 7.69 mb
@@ -1096,7 +1096,7 @@ Double_t AliGenMUONlib::YCharmF0M1S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M2S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M2S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 4.81 mb
@@ -1114,7 +1114,7 @@ Double_t AliGenMUONlib::YCharmF0M2S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S1PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S1PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 14.09 mb
@@ -1132,7 +1132,7 @@ Double_t AliGenMUONlib::YCharmF0M0S1PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S2PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S2PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 1.52 mb
@@ -1150,7 +1150,7 @@ Double_t AliGenMUONlib::YCharmF0M0S2PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S3PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S3PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 3.67 mb
@@ -1168,7 +1168,7 @@ Double_t AliGenMUONlib::YCharmF0M0S3PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S4PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S4PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 3.38 mb
@@ -1186,7 +1186,7 @@ Double_t AliGenMUONlib::YCharmF0M0S4PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S5PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S5PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 10.37 mb
@@ -1204,7 +1204,7 @@ Double_t AliGenMUONlib::YCharmF0M0S5PP( Double_t *px, Double_t */*dummy*/)
     
     return ycharm;
 }
-Double_t AliGenMUONlib::YCharmF0M0S6PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YCharmF0M0S6PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Charm y
 // Corresponding NLO total cross section is 7.22 mb
@@ -1261,7 +1261,7 @@ Int_t AliGenMUONlib::IpCharm(TRandom *ran)
 //
 //                    pt-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::PtBeauty( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeauty( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
   const Double_t kpt0 = 6.53;
@@ -1272,7 +1272,7 @@ Double_t AliGenMUONlib::PtBeauty( Double_t *px, Double_t */*dummy*/)
   return x/TMath::Power(pass1,kxn);
 }
 
-Double_t AliGenMUONlib::PtBeautyCentral( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyCentral( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
   const Double_t kpt0 = 6.14;
@@ -1282,7 +1282,7 @@ Double_t AliGenMUONlib::PtBeautyCentral( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // FiMjSkPP define theoretical uncertainties around F0M0S0PP as follows:
 // PtBeautyFiMjSkPP = PtBeautyF0M0S0PP * (dN(i,j,k)/dpt / dN(0,0,0)/dpt)_MNR
@@ -1311,7 +1311,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF1M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF1M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.445 mb
@@ -1322,7 +1322,7 @@ Double_t AliGenMUONlib::PtBeautyF1M0S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF2M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF2M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.445 mb
@@ -1333,7 +1333,7 @@ Double_t AliGenMUONlib::PtBeautyF2M0S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M1S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M1S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.518 mb
@@ -1344,7 +1344,7 @@ Double_t AliGenMUONlib::PtBeautyF0M1S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M2S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M2S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.384 mb
@@ -1355,7 +1355,7 @@ Double_t AliGenMUONlib::PtBeautyF0M2S0PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S1PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S1PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.648 mb
@@ -1366,7 +1366,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S1PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S2PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S2PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.294 mb
@@ -1377,7 +1377,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S2PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S3PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S3PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.475 mb
@@ -1388,7 +1388,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S3PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S4PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S4PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.324 mb
@@ -1399,7 +1399,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S4PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S5PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S5PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.536 mb
@@ -1410,7 +1410,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S5PP( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtBeautyF0M0S6PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtBeautyF0M0S6PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty pT
 // Corresponding NLO total cross section is 0.420 mb
@@ -1423,7 +1423,7 @@ Double_t AliGenMUONlib::PtBeautyF0M0S6PP( Double_t *px, Double_t */*dummy*/)
 }
 
 //                     y-distribution
-Double_t AliGenMUONlib::YBeauty( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeauty( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y :: Carrer & Dainese : ALICE-INT-2003-019 v.3 (hep-ph/0311225) 
 // Pythia tuned to reproduce the distribution given by the HVQMNR program based on NLO calculations (pQCD)
@@ -1443,7 +1443,7 @@ Double_t AliGenMUONlib::YBeauty( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // FiMjSkPP define theoretical uncertainties around F0M0S0PP as follows:
 // YBeautyFiMjSkPP = YBeautyF0M0S0PP * (dN(i,j,k)/dy / dN(0,0,0)/dy)_MNR
@@ -1480,7 +1480,7 @@ Double_t AliGenMUONlib::YBeautyF0M0S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF1M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF1M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.445 mb
@@ -1498,7 +1498,7 @@ Double_t AliGenMUONlib::YBeautyF1M0S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF2M0S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF2M0S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.445 mb
@@ -1516,7 +1516,7 @@ Double_t AliGenMUONlib::YBeautyF2M0S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M1S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M1S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.518 mb
@@ -1534,7 +1534,7 @@ Double_t AliGenMUONlib::YBeautyF0M1S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M2S0PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M2S0PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.384 mb
@@ -1552,7 +1552,7 @@ Double_t AliGenMUONlib::YBeautyF0M2S0PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S1PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S1PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.648 mb
@@ -1570,7 +1570,7 @@ Double_t AliGenMUONlib::YBeautyF0M0S1PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S2PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S2PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.294 mb
@@ -1588,7 +1588,7 @@ Double_t AliGenMUONlib::YBeautyF0M0S2PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S3PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S3PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.475 mb
@@ -1606,7 +1606,7 @@ Double_t AliGenMUONlib::YBeautyF0M0S3PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S4PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S4PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.324 mb
@@ -1624,7 +1624,7 @@ Double_t AliGenMUONlib::YBeautyF0M0S4PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S5PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S5PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.536 mb
@@ -1642,7 +1642,7 @@ Double_t AliGenMUONlib::YBeautyF0M0S5PP( Double_t *px, Double_t */*dummy*/)
     
     return ybeauty;
 }
-Double_t AliGenMUONlib::YBeautyF0M0S6PP( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YBeautyF0M0S6PP( const Double_t *px, const Double_t */*dummy*/)
 {
 // Beauty y
 // Corresponding NLO total cross section is 0.420 mb
@@ -1693,7 +1693,7 @@ Int_t AliGenMUONlib::IpBeauty(TRandom *ran)
 }
 
 
-typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
+typedef Double_t (*GenFunc) (const Double_t*,  const Double_t*);
 GenFunc AliGenMUONlib::GetPt(Int_t param,  const char* tname) const
 {
 // Return pointer to pT parameterisation
@@ -1984,7 +1984,7 @@ GenFunc AliGenMUONlib::GetY(Int_t param, const char* tname) const
 //
 //                pt-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::PtChic0( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtChic0( const Double_t *px, const Double_t */*dummy*/)
 {
 // Chi_c1 pT
   const Double_t kpt0 = 4.;
@@ -1994,7 +1994,7 @@ Double_t AliGenMUONlib::PtChic0( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtChic1( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtChic1( const Double_t *px, const Double_t */*dummy*/)
 {
 // Chi_c1 pT
   const Double_t kpt0 = 4.;
@@ -2004,7 +2004,7 @@ Double_t AliGenMUONlib::PtChic1( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtChic2( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtChic2( const Double_t *px, const Double_t */*dummy*/)
 {
 // Chi_c2 pT
   const Double_t kpt0 = 4.;
@@ -2014,7 +2014,7 @@ Double_t AliGenMUONlib::PtChic2( Double_t *px, Double_t */*dummy*/)
   Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
-Double_t AliGenMUONlib::PtChic( Double_t *px, Double_t */*dummy*/)
+Double_t AliGenMUONlib::PtChic( const Double_t *px, const Double_t */*dummy*/)
 {
 // Chi_c family pT
   const Double_t kpt0 = 4.;
@@ -2028,7 +2028,7 @@ Double_t AliGenMUONlib::PtChic( Double_t *px, Double_t */*dummy*/)
 //
 //               y-distribution
 //____________________________________________________________
-Double_t AliGenMUONlib::YChic0(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YChic0(const Double_t *py, const Double_t */*dummy*/)
 {
 // Chi-1c y
   const Double_t ky0 = 4.;
@@ -2043,7 +2043,7 @@ Double_t AliGenMUONlib::YChic0(Double_t *py, Double_t */*dummy*/)
   return yj;
 }
 
-Double_t AliGenMUONlib::YChic1(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YChic1(const Double_t *py, const Double_t */*dummy*/)
 {
 // Chi-1c y
   const Double_t ky0 = 4.;
@@ -2058,7 +2058,7 @@ Double_t AliGenMUONlib::YChic1(Double_t *py, Double_t */*dummy*/)
   return yj;
 }
 
-Double_t AliGenMUONlib::YChic2(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YChic2(const Double_t *py, const Double_t */*dummy*/)
 {
 // Chi-2c y
   const Double_t ky0 = 4.;
@@ -2073,7 +2073,7 @@ Double_t AliGenMUONlib::YChic2(Double_t *py, Double_t */*dummy*/)
   return yj;
 }
 
-Double_t AliGenMUONlib::YChic(Double_t *py, Double_t */*dummy*/)
+Double_t AliGenMUONlib::YChic(const Double_t *py, const Double_t */*dummy*/)
 {
 // Chi_c family y
   const Double_t ky0 = 4.;
index 5dd9f56..c609ad9 100644 (file)
@@ -26,141 +26,141 @@ class AliGenMUONlib :
  private:
     
 // pions
-    static Double_t PtPion(Double_t *px, Double_t *dummy);
+    static Double_t PtPion(const Double_t *px, const Double_t *dummy);
     static Double_t PtScal(Double_t pt, Int_t np);
-    static Double_t YPion( Double_t *py, Double_t *dummy);
+    static Double_t YPion( const Double_t *py, const Double_t *dummy);
     static Int_t    IpPion(TRandom *ran);
 // kaons
-    static Double_t PtKaon(Double_t *px, Double_t *dummy);
-    static Double_t YKaon( Double_t *py, Double_t *dummy);
+    static Double_t PtKaon(const Double_t *px, const Double_t *dummy);
+    static Double_t YKaon( const Double_t *py, const Double_t *dummy);
     static Int_t    IpKaon(TRandom *ran);
 // Phi
-    static Double_t PtPhi( Double_t *px, Double_t *dummy);
-    static Double_t YPhi( Double_t *px, Double_t *dummy);
+    static Double_t PtPhi( const Double_t *px, const Double_t *dummy);
+    static Double_t YPhi( const  Double_t *px, const Double_t *dummy);
     static Int_t    IpPhi(TRandom *ran);
 // Omega
-    static Double_t PtOmega( Double_t *px, Double_t *dummy);
-    static Double_t YOmega( Double_t *px, Double_t *dummy);
+    static Double_t PtOmega( const Double_t *px, const Double_t *dummy);
+    static Double_t YOmega( const Double_t *px, const Double_t *dummy);
     static Int_t    IpOmega(TRandom *ran);
 // Eta
-    static Double_t PtEta( Double_t *px, Double_t *dummy);
-    static Double_t YEta( Double_t *px, Double_t *dummy);
+    static Double_t PtEta( const Double_t *px, const Double_t *dummy);
+    static Double_t YEta( const Double_t *px, const Double_t *dummy);
     static Int_t    IpEta(TRandom *ran);
 // 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 PtJpsiCDFscaledPP10( Double_t *px, Double_t *dummy);
-    static Double_t PtJpsiCDFscaledold( Double_t *px, Double_t *dummy);
-    static Double_t PtJpsiCDFscaledPPold( 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 YJpsiCDFscaledPP10( Double_t *px, Double_t *dummy);
-    static Double_t YJpsiCDFscaledold(Double_t *py, Double_t *dummy);
-    static Double_t YJpsiCDFscaledPPold( 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);
+    static Double_t PtJpsi( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiCDFscaled( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiCDFscaledPP( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiCDFscaledPP10( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiCDFscaledold( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiCDFscaledPPold( const Double_t *px, const Double_t *dummy);
+    static Double_t YJpsi(const Double_t *py, const Double_t *dummy);
+    static Double_t PtJpsiPbPb( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiBPbPb( const Double_t *px, const Double_t *dummy);
+    static Double_t YJpsiPbPb(const Double_t *py, const Double_t *dummy);
+    static Double_t YJpsiCDFscaled(const Double_t *py, const Double_t *dummy);
+    static Double_t YJpsiCDFscaledPP( const Double_t *px, const Double_t *dummy);
+    static Double_t YJpsiCDFscaledPP10( const Double_t *px, const Double_t *dummy);
+    static Double_t YJpsiCDFscaledold(const Double_t *py, const Double_t *dummy);
+    static Double_t YJpsiCDFscaledPPold( const Double_t *px, const Double_t *dummy);
+    static Double_t PtJpsiPP( const Double_t *px, const Double_t *dummy);
+    static Double_t YJpsiPP(const Double_t *py, const Double_t *dummy);
+    static Double_t YJpsiBPbPb(const Double_t *py, const Double_t *dummy);
     static Int_t    IpJpsi(TRandom *ran);
     static Int_t    IpJpsiFamily(TRandom *ran);
     static Int_t    IpPsiP(TRandom *ran);
-    static Double_t PtJpsiFlat( Double_t *px, Double_t *dummy );
-    static Double_t YJpsiFlat(Double_t *py, Double_t *dummy);
+    static Double_t PtJpsiFlat( const Double_t *px, const Double_t *dummy );
+    static Double_t YJpsiFlat(const Double_t *py, const Double_t *dummy);
 
 // 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 );
-    static Double_t YUpsilonPP(Double_t *py, Double_t *dummy);
+    static Double_t PtUpsilon( const Double_t *px, const Double_t *dummy );
+    static Double_t PtUpsilonCDFscaled( const Double_t *px, const Double_t *dummy );
+    static Double_t PtUpsilonCDFscaledPP( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilon(const Double_t *py, const Double_t *dummy);
+    static Double_t YUpsilonCDFscaled(const Double_t *py, const Double_t *dummy);
+    static Double_t YUpsilonCDFscaledPP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtUpsilonPbPb( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonPbPb(const Double_t *py, const Double_t *dummy);
+    static Double_t PtUpsilonPP( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonPP(const Double_t *py, const Double_t *dummy);
     static Int_t    IpUpsilon(TRandom *ran);
     static Int_t    IpUpsilonFamily(TRandom *ran);
     static Int_t    IpUpsilonP(TRandom *ran);
     static Int_t    IpUpsilonPP(TRandom *ran);
-    static Double_t PtUpsilonFlat( Double_t *px, Double_t *dummy );
-    static Double_t YUpsilonFlat(Double_t *py, Double_t *dummy);
+    static Double_t PtUpsilonFlat( const Double_t *px, const Double_t *dummy );
+    static Double_t YUpsilonFlat(const Double_t *py, const Double_t *dummy);
 //
 // Charm    
-    static Double_t PtCharm( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmCentral( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF1M0S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF2M0S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M1S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M2S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S1PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S2PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S3PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S4PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S5PP( Double_t *px, Double_t *dummy );
-    static Double_t PtCharmF0M0S6PP( Double_t *px, Double_t *dummy );
-    static Double_t YCharm(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF1M0S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF2M0S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M1S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M2S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S1PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S2PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S3PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S4PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S5PP(Double_t *py, Double_t *dummy);
-    static Double_t YCharmF0M0S6PP(Double_t *py, Double_t *dummy);
+    static Double_t PtCharm( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmCentral( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF1M0S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF2M0S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M1S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M2S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S1PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S2PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S3PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S4PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S5PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtCharmF0M0S6PP( const Double_t *px, const Double_t *dummy );
+    static Double_t YCharm(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF1M0S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF2M0S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M1S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M2S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S1PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S2PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S3PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S4PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S5PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YCharmF0M0S6PP(const Double_t *py, const Double_t *dummy);
     static Int_t    IpCharm(TRandom *ran);
 //
 // Beauty
-    static Double_t PtBeauty( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF1M0S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF2M0S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M1S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M2S0PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S1PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S2PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S3PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S4PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S5PP( Double_t *px, Double_t *dummy );
-    static Double_t PtBeautyF0M0S6PP( Double_t *px, Double_t *dummy );
-    static Double_t YBeauty(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF1M0S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF2M0S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M1S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M2S0PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S1PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S2PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S3PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S4PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S5PP(Double_t *py, Double_t *dummy);
-    static Double_t YBeautyF0M0S6PP(Double_t *py, Double_t *dummy);
-    static Double_t PtBeautyCentral( Double_t *px, Double_t *dummy );
+    static Double_t PtBeauty( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF1M0S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF2M0S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M1S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M2S0PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S1PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S2PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S3PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S4PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S5PP( const Double_t *px, const Double_t *dummy );
+    static Double_t PtBeautyF0M0S6PP( const Double_t *px, const Double_t *dummy );
+    static Double_t YBeauty(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF1M0S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF2M0S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M1S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M2S0PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S1PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S2PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S3PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S4PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S5PP(const Double_t *py, const Double_t *dummy);
+    static Double_t YBeautyF0M0S6PP(const Double_t *py, const Double_t *dummy);
+    static Double_t PtBeautyCentral( const Double_t *px, const Double_t *dummy );
     static Int_t    IpBeauty(TRandom *ran);
 //
 
    // Chi 1c 2c
-   static Double_t PtChic0( Double_t *px, Double_t *dummy);
-   static Double_t YChic0(Double_t *py, Double_t *dummy);
+   static Double_t PtChic0( const Double_t *px, const Double_t *dummy);
+   static Double_t YChic0(const Double_t *py, const Double_t *dummy);
    static Int_t    IpChic0(TRandom *ran);
 
-   static Double_t PtChic1( Double_t *px, Double_t *dummy);
-   static Double_t YChic1(Double_t *py, Double_t *dummy);
+   static Double_t PtChic1( const Double_t *px, const Double_t *dummy);
+   static Double_t YChic1(const Double_t *py, const Double_t *dummy);
    static Int_t    IpChic1(TRandom *ran);
 
-   static Double_t PtChic2( Double_t *px, Double_t *dummy);
-   static Double_t YChic2(Double_t *py, Double_t *dummy);
+   static Double_t PtChic2( const Double_t *px, const Double_t *dummy);
+   static Double_t YChic2(const Double_t *py, const Double_t *dummy);
    static Int_t    IpChic2(TRandom *ran);
 
-   static Double_t PtChic( Double_t *px, Double_t *dummy);
-   static Double_t YChic(Double_t *py, Double_t *dummy);
+   static Double_t PtChic( const Double_t *px, const Double_t *dummy);
+   static Double_t YChic(const Double_t *py, const Double_t *dummy);
    static Int_t    IpChic(TRandom *ran);
 
 //
index 00d50c3..4d4c985 100644 (file)
@@ -55,7 +55,7 @@ ClassImp(AliGenPHOSlib)
 //    Rapidity distribution YPion
 //    Particle distribution IdPion  111, 211 and -211 (pi0, pi+ and pi-)
 //
- Double_t AliGenPHOSlib::PtPion(Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtPion(const Double_t *px, const Double_t *)
 {
 //     Pion transverse momentum distribtuion taken 
 //     from AliGenMUONlib class, version 3.01 of aliroot
@@ -82,7 +82,7 @@ ClassImp(AliGenPHOSlib)
     y=kb*TMath::Exp(-sqrt(x*x+kxmpi2)/kt);
   return y*x;
 }
- Double_t AliGenPHOSlib::YPion( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YPion( const Double_t *py, const Double_t *)
 {
 //
 // pion y-distribution
@@ -145,7 +145,7 @@ ClassImp(AliGenPHOSlib)
 //    Particle distribution IdPi0Flat  111 (pi0)
 //
 
-Double_t AliGenPHOSlib::PtPi0Flat(Double_t */*px*/, Double_t *)
+Double_t AliGenPHOSlib::PtPi0Flat(const Double_t */*px*/, const Double_t *)
 {
 //     Pion transverse momentum flat distribution 
 
@@ -153,7 +153,7 @@ return 1;
 
 }
 
-Double_t AliGenPHOSlib::YPi0Flat( Double_t */*py*/, Double_t *)
+Double_t AliGenPHOSlib::YPi0Flat( const Double_t */*py*/, const Double_t *)
 {
 
 // pion y-distribution
@@ -201,7 +201,7 @@ Double_t AliGenPHOSlib::YPi0Flat( Double_t */*py*/, Double_t *)
 // End Scaling
 //============================================================================
 //    K  A  O  N  S
- Double_t AliGenPHOSlib::PtKaon( Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtKaon( const Double_t *px, const Double_t *)
 {
 //                kaon
 //                pt-distribution
@@ -210,7 +210,7 @@ Double_t AliGenPHOSlib::YPi0Flat( Double_t */*py*/, Double_t *)
   return PtScal(*px,2);  //  2==> Kaon in the PtScal function
 }
 
- Double_t AliGenPHOSlib::YKaon( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YKaon( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -276,7 +276,7 @@ Int_t AliGenPHOSlib::IpKaon0L(TRandom *)
 //============================================================================
 //============================================================================
 //   E  T  A  S
- Double_t AliGenPHOSlib::PtEta( Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtEta( const Double_t *px, const Double_t *)
 {
 //                etas
 //                pt-distribution
@@ -285,7 +285,7 @@ Int_t AliGenPHOSlib::IpKaon0L(TRandom *)
   return PtScal(*px,3);  //  3==> Eta in the PtScal function
 }
 
- Double_t AliGenPHOSlib::YEta( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YEta( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -316,7 +316,7 @@ Int_t AliGenPHOSlib::IpKaon0L(TRandom *)
 //    Particle distribution IdEtaFlat  111 (pi0)
 //
 
-Double_t AliGenPHOSlib::PtEtaFlat(Double_t */*px*/, Double_t *)
+Double_t AliGenPHOSlib::PtEtaFlat(const Double_t */*px*/, const Double_t *)
 {
 //     Eta transverse momentum flat distribution 
 
@@ -324,7 +324,7 @@ Double_t AliGenPHOSlib::PtEtaFlat(Double_t */*px*/, Double_t *)
 
 }
 
-Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
+Double_t AliGenPHOSlib::YEtaFlat( const Double_t */*py*/, const Double_t *)
 {
 //
 // pion y-distribution
@@ -343,7 +343,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
 //============================================================================
 //============================================================================
 //    O  M  E  G  A  S
- Double_t AliGenPHOSlib::PtOmega( Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtOmega( const Double_t *px, const Double_t *)
 {
 // omegas
 //                pt-distribution
@@ -352,7 +352,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
   return PtScal(*px,4);  //  4==> Omega in the PtScal function
 }
 
- Double_t AliGenPHOSlib::YOmega( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YOmega( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -378,7 +378,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
 //============================================================================
 //============================================================================
 //    E  T  A  P  R  I  M  E
- Double_t AliGenPHOSlib::PtEtaprime( Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtEtaprime( const Double_t *px, const Double_t *)
 {
 // etaprime
 //                pt-distribution
@@ -387,7 +387,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
   return PtScal(*px,5);  //  5==> Etaprime in the PtScal function
 }
 
- Double_t AliGenPHOSlib::YEtaprime( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YEtaprime( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -413,7 +413,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
 //===================================================================
 //============================================================================
 //    P  H  I   S
- Double_t AliGenPHOSlib::PtPhi( Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtPhi( const Double_t *px, const Double_t *)
 {
 // phi
 //                pt-distribution
@@ -422,7 +422,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
   return PtScal(*px,6);  //  6==> Phi in the PtScal function
 }
 
- Double_t AliGenPHOSlib::YPhi( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YPhi( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -448,7 +448,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
 //===================================================================
 //============================================================================
 //    B  A  R  Y  O  N  S  == protons, protonsbar, neutrons, and neutronsbars
- Double_t AliGenPHOSlib::PtBaryon( Double_t *px, Double_t *)
+ Double_t AliGenPHOSlib::PtBaryon( const Double_t *px, const Double_t *)
 {
 // baryons
 //                pt-distribution
@@ -457,7 +457,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
   return PtScal(*px,7);  //  7==> Baryon in the PtScal function
 }
 
- Double_t AliGenPHOSlib::YBaryon( Double_t *py, Double_t *)
+ Double_t AliGenPHOSlib::YBaryon( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -531,7 +531,7 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t */*py*/, Double_t *)
 
 
 
-typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
+typedef Double_t (*GenFunc) (const Double_t*,  const Double_t*);
 GenFunc AliGenPHOSlib::GetPt(Int_t param, const char* /*tname*/) const
 {
 // Return pinter to pT parameterisation
index 057ce74..3c60801 100644 (file)
@@ -30,51 +30,51 @@ public AliGenLib
     GenFuncIp GetIp(Int_t param, const char* tname=0) const;   
  private:
 // pions
-    static Double_t PtPion(Double_t *px, Double_t *dummy);
+    static Double_t PtPion(const Double_t *px, const Double_t *dummy);
     static Double_t PtScal(Double_t pt, Int_t np);
-    static Double_t YPion( Double_t *py, Double_t *dummy);
+    static Double_t YPion( const Double_t *py, const Double_t *dummy);
     static Int_t    IpPion(TRandom* ran);
     static Int_t    IpChargedPion(TRandom* ran);
 
 //  pi0 Flat Distribution
-    static Double_t PtPi0Flat(Double_t *px, Double_t *dummy);
-    static Double_t YPi0Flat( Double_t *py, Double_t *dummy);
+    static Double_t PtPi0Flat(const Double_t *px, const Double_t *dummy);
+    static Double_t YPi0Flat( const Double_t *py, const Double_t *dummy);
     static Int_t    IpPi0Flat(TRandom* ran); 
     
 // kaons
-    static Double_t PtKaon(Double_t *px, Double_t *dummy);
-    static Double_t YKaon( Double_t *py, Double_t *dummy);
+    static Double_t PtKaon(const Double_t *px, const Double_t *dummy);
+    static Double_t YKaon( const Double_t *py, const Double_t *dummy);
     static Int_t    IpKaon(TRandom* ran);
     static Int_t    IpChargedKaon(TRandom* ran);
     static Int_t    IpKaon0L(TRandom* ran);
 // etas
-    static Double_t PtEta(Double_t *px, Double_t *dummy);
-    static Double_t YEta( Double_t *py, Double_t *dummy);
+    static Double_t PtEta(const Double_t *px, const Double_t *dummy);
+    static Double_t YEta( const Double_t *py, const Double_t *dummy);
     static Int_t    IpEta(TRandom *ran);
     
 // etas Flat Distribution
-    static Double_t PtEtaFlat(Double_t *px, Double_t *dummy);
-    static Double_t YEtaFlat( Double_t *py, Double_t *dummy);
+    static Double_t PtEtaFlat(const Double_t *px, const Double_t *dummy);
+    static Double_t YEtaFlat( const Double_t *py, const Double_t *dummy);
     static Int_t    IpEtaFlat(TRandom *ran);
 
 // omegas
-    static Double_t PtOmega(Double_t *px, Double_t *dummy);
-    static Double_t YOmega( Double_t *py, Double_t *dummy);
+    static Double_t PtOmega(const Double_t *px, const Double_t *dummy);
+    static Double_t YOmega( const Double_t *py, const Double_t *dummy);
     static Int_t    IpOmega(TRandom *ran);
     
 // etaprime
-    static Double_t PtEtaprime(Double_t *px, Double_t *dummy);
-    static Double_t YEtaprime( Double_t *py, Double_t *dummy);
+    static Double_t PtEtaprime(const Double_t *px, const Double_t *dummy);
+    static Double_t YEtaprime( const Double_t *py, const Double_t *dummy);
     static Int_t    IpEtaprime(TRandom* ran);
     
 // phis
-    static Double_t PtPhi(Double_t *px, Double_t *dummy);
-    static Double_t YPhi( Double_t *py, Double_t *dummy);
+    static Double_t PtPhi(const Double_t *px, const Double_t *dummy);
+    static Double_t YPhi( const Double_t *py, const Double_t *dummy);
     static Int_t    IpPhi(TRandom* ran);
     
 // baryons
-    static Double_t PtBaryon(Double_t *px, Double_t *dummy);
-    static Double_t YBaryon( Double_t *py, Double_t *dummy);
+    static Double_t PtBaryon(const Double_t *px, const Double_t *dummy);
+    static Double_t YBaryon( const Double_t *py, const Double_t *dummy);
     static Int_t    IpBaryon(TRandom *ran);
     static Int_t    IpProton(TRandom *ran);
     static Int_t    IpAProton(TRandom *ran);
index a2df71c..0d91730 100644 (file)
@@ -32,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
@@ -62,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
@@ -114,7 +114,7 @@ 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);
@@ -122,7 +122,7 @@ Double_t AliGenPMDlib::PtEta( Double_t *px, Double_t */*dummy*/)
 
 // 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
@@ -148,7 +148,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
index 6352550..c0b587c 100644 (file)
@@ -24,13 +24,13 @@ class AliGenPMDlib :public AliGenLib
     GenFuncIp GetIp(Int_t param, const char* tname=0) const;    
  private:
 // Neutral pions
-    static Double_t PtPi0(Double_t *px, Double_t *dummy);
+    static Double_t PtPi0(const Double_t *px, const Double_t *dummy);
     static Double_t PtScal(Double_t pt, Int_t np);
-    static Double_t YPi0( Double_t *py, Double_t *dummy);
+    static Double_t YPi0( const Double_t *py, const Double_t *dummy);
     static Int_t    IpPi0(TRandom *ran);
 // Etas
-    static Double_t PtEta(Double_t *px, Double_t *dummy);
-    static Double_t YEta( Double_t *py, Double_t *dummy);
+    static Double_t PtEta(const Double_t *px, const Double_t *dummy);
+    static Double_t YEta( const Double_t *py, const Double_t *dummy);
     static Int_t    IpEta(TRandom *ran);
 //
     ClassDef(AliGenPMDlib,0) // Library providing y and pT parameterisations
index 59d117e..a72a584 100644 (file)
@@ -131,8 +131,8 @@ AliGenParam::AliGenParam(Int_t npart, Int_t param, const char* tname, const char
 //____________________________________________________________
 
 AliGenParam::AliGenParam(Int_t npart, Int_t param,
-                         Double_t (*PtPara) (Double_t*, Double_t*),
-                         Double_t (*YPara ) (Double_t* ,Double_t*),
+                         Double_t (*PtPara) (const Double_t*, const Double_t*),
+                         Double_t (*YPara ) (const Double_t* ,const Double_t*),
                         Int_t    (*IpPara) (TRandom *))                 
     :AliGenMC(npart),
      
index 4b978ba..4921476 100644 (file)
@@ -29,8 +29,8 @@ class AliGenParam : public AliGenMC
     AliGenParam(Int_t npart, AliGenLib * Library, Int_t param, char* tname = 0);
     AliGenParam(Int_t npart, Int_t param, const char* tname = 0, const char* name = 0);
     AliGenParam(Int_t npart, Int_t param,
-               Double_t (*PtPara)(Double_t*, Double_t*),
-               Double_t (*YPara )(Double_t*, Double_t*),
+               Double_t (*PtPara)(const Double_t*, const Double_t*),
+               Double_t (*YPara )(const Double_t*, const Double_t*),
                Int_t    (*IpPara)(TRandom*)           );
      
     virtual ~AliGenParam();
@@ -48,8 +48,8 @@ class AliGenParam : public AliGenMC
     Float_t GetRelativeArea(Float_t ptMin, Float_t ptMax, Float_t yMin, Float_t yMax, Float_t phiMin, Float_t phiMax);
 
  protected:
-    Double_t (*fPtParaFunc)(Double_t*, Double_t*); //! Pointer to Pt parametrisation function
-    Double_t (*fYParaFunc )(Double_t*, Double_t*); //! Pointer to Y parametrisation function
+    Double_t (*fPtParaFunc)(const Double_t*, const Double_t*); //! Pointer to Pt parametrisation function
+    Double_t (*fYParaFunc )(const Double_t*, const Double_t*); //! Pointer to Y parametrisation function
     Int_t    (*fIpParaFunc )(TRandom*);    //! Pointer to particle type parametrisation function
     TF1* fPtPara;              // Transverse momentum parameterisation
     TF1* fYPara;               // Rapidity parameterisation
index 49d5d30..6784237 100644 (file)
@@ -84,7 +84,7 @@ ClassImp(AliGenSTRANGElib)
 }
 //============================================================= 
 //
- Double_t AliGenSTRANGElib::PtPion(Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtPion(const Double_t *px, const Double_t *)
 {
 //     Pion transverse momentum distribtuion taken 
 //     from AliGenMUONlib class, version 3.01 of aliroot
@@ -114,7 +114,7 @@ ClassImp(AliGenSTRANGElib)
 // End Scaling
 //============================================================================
 //    K  A  O  N  
- Double_t AliGenSTRANGElib::PtKaon( Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtKaon( const Double_t *px, const Double_t *)
 {
 //                kaon
 //                pt-distribution
@@ -123,7 +123,7 @@ ClassImp(AliGenSTRANGElib)
   return PtScal(*px,2);  //  2==> Kaon in the PtScal function
 }
 
- Double_t AliGenSTRANGElib::YKaon( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YKaon( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -166,7 +166,7 @@ ClassImp(AliGenSTRANGElib)
 //============================================================================
 //============================================================================
 //    P  H  I   
- Double_t AliGenSTRANGElib::PtPhi( Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtPhi( const Double_t *px, const Double_t *)
 {
 // phi
 //                pt-distribution
@@ -175,7 +175,7 @@ ClassImp(AliGenSTRANGElib)
   return PtScal(*px,6);  //  6==> Phi in the PtScal function
 }
 
- Double_t AliGenSTRANGElib::YPhi( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YPhi( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -201,7 +201,7 @@ ClassImp(AliGenSTRANGElib)
 //===================================================================
 //============================================================================
 //    Lambda
- Double_t AliGenSTRANGElib::PtLambda( Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtLambda( const Double_t *px, const Double_t *)
 {
 // Lambda
 //                pt-distribution
@@ -210,7 +210,7 @@ ClassImp(AliGenSTRANGElib)
   return PtScal(*px,8);  //  8==> Lambda-antiLambda in the PtScal function
 }
 
- Double_t AliGenSTRANGElib::YLambda( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YLambda( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -240,7 +240,7 @@ ClassImp(AliGenSTRANGElib)
 // End Lambda
 //============================================================================
 //    XIminus
- Double_t AliGenSTRANGElib::PtXiMinus( Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtXiMinus( const Double_t *px, const Double_t *)
 {
 // Xi
 //                pt-distribution
@@ -249,7 +249,7 @@ ClassImp(AliGenSTRANGElib)
   return PtScal(*px,9);  //  9==> Xi-antiXi in the PtScal function
 }
 
- Double_t AliGenSTRANGElib::YXiMinus( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YXiMinus( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -279,7 +279,7 @@ ClassImp(AliGenSTRANGElib)
 // End Ximinus
 //============================================================================
 //    Omegaminus
- Double_t AliGenSTRANGElib::PtOmegaMinus( Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtOmegaMinus( const Double_t *px, const Double_t *)
 {
 // Omega
 //                pt-distribution
@@ -288,7 +288,7 @@ ClassImp(AliGenSTRANGElib)
   return PtScal(*px,10);  //  10==> Omega-antiOmega in the PtScal function
 }
 
- Double_t AliGenSTRANGElib::YOmegaMinus( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YOmegaMinus( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -319,7 +319,7 @@ ClassImp(AliGenSTRANGElib)
 // End Omegaminus
 //============================================================================
 //     Lambda(1520)
-Double_t AliGenSTRANGElib::PtLambda1520( Double_t *px, Double_t *)
+Double_t AliGenSTRANGElib::PtLambda1520( const Double_t *px, const Double_t *)
 {
 // Lambda(1520)
 //                  pt-distribution
@@ -328,7 +328,7 @@ Double_t AliGenSTRANGElib::PtLambda1520( Double_t *px, Double_t *)
   return PtScal(*px,11);   //   11=> Lambda(1520) in the PtScal function
 }
 
-Double_t AliGenSTRANGElib::YLambda1520( Double_t *py, Double_t *)
+Double_t AliGenSTRANGElib::YLambda1520( const Double_t *py, const Double_t *)
 {
 // y-distribution
 //____________________________________________________________
@@ -359,7 +359,7 @@ Int_t AliGenSTRANGElib::IpLambda1520(TRandom * ran)
 // End Lambda(1520)
 //============================================================================
 
-typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
+typedef Double_t (*GenFunc) (const Double_t*,  const Double_t*);
  GenFunc AliGenSTRANGElib::GetPt(Int_t param, const char* /*tname*/) const
 {
 // Return pinter to pT parameterisation
index 55eb671..af9f167 100644 (file)
@@ -27,31 +27,31 @@ public AliGenLib
     GenFuncIp GetIp(Int_t param, const char* tname=0) const;    
  private:
 // pions
-    static Double_t PtPion(Double_t *px, Double_t *dummy);
+    static Double_t PtPion(const Double_t *px, const Double_t *dummy);
     static Double_t PtScal(Double_t pt, Int_t np);
 // kaons
-    static Double_t PtKaon(Double_t *px, Double_t *dummy);
-    static Double_t YKaon( Double_t *py, Double_t *dummy);
+    static Double_t PtKaon(const Double_t *px, const Double_t *dummy);
+    static Double_t YKaon( const Double_t *py, const Double_t *dummy);
     static Int_t    IpKaon(TRandom* ran);
 // phis
-    static Double_t PtPhi(Double_t *px, Double_t *dummy);
-    static Double_t YPhi( Double_t *py, Double_t *dummy);
+    static Double_t PtPhi(const Double_t *px, const Double_t *dummy);
+    static Double_t YPhi( const Double_t *py, const Double_t *dummy);
     static Int_t    IpPhi(TRandom* ran);
 // lambda
-    static Double_t PtLambda(Double_t *px, Double_t *dummy);
-    static Double_t YLambda( Double_t *py, Double_t *dummy);
+    static Double_t PtLambda(const Double_t *px, const Double_t *dummy);
+    static Double_t YLambda( const Double_t *py, const Double_t *dummy);
     static Int_t    IpLambda(TRandom *ran);
 // Ximinus
-    static Double_t PtXiMinus(Double_t *px, Double_t *dummy);
-    static Double_t YXiMinus( Double_t *py, Double_t *dummy);
+    static Double_t PtXiMinus(const Double_t *px, const Double_t *dummy);
+    static Double_t YXiMinus( const Double_t *py, const Double_t *dummy);
     static Int_t    IpXiMinus(TRandom *ran);
 // Omegaminus
-    static Double_t PtOmegaMinus(Double_t *px, Double_t *dummy);
-    static Double_t YOmegaMinus( Double_t *py, Double_t *dummy);
+    static Double_t PtOmegaMinus(const Double_t *px, const Double_t *dummy);
+    static Double_t YOmegaMinus( const Double_t *py, const Double_t *dummy);
     static Int_t    IpOmegaMinus(TRandom *ran);
 // Lambda(1520)
-    static Double_t PtLambda1520(Double_t *px, Double_t *dummy);
-    static Double_t YLambda1520(Double_t *py, Double_t *dummy);
+    static Double_t PtLambda1520(const Double_t *px, const Double_t *dummy);
+    static Double_t YLambda1520(const Double_t *py, const Double_t *dummy);
     static Int_t    IpLambda1520(TRandom *ran);
     
     ClassDef(AliGenSTRANGElib,0) // Library providing y and pT parameterisations