Paramterisation of CDF J/Psi data used. (G. Martinez)

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

diff --git a/EVGEN/AliGenGSIlib.cxx b/EVGEN/AliGenGSIlib.cxx
index 0057327635e487f3e1767fa0a6c34ac06d24aad6..799e1fb3cc3626bdaf12c256bf71843195682a8c 100644 (file)
--- a/EVGEN/AliGenGSIlib.cxx
+++ b/EVGEN/AliGenGSIlib.cxx
@@ -13,37 +13,7 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
   
-/*
-$Log$
-Revision 1.5.10.1  2002/06/10 14:57:41  hristov
-Merged with v3-08-02
-
-Revision 1.6  2002/04/16 15:53:19  morsch
-Coding rule violations corrected.
-
-Revision 1.5  2001/05/15 15:43:23  morsch
-Update of library for low mass resonances (Yiota Foka)
-
-Revision 1.4  2001/03/09 13:01:41  morsch
-- enum constants for paramterisation type (particle family) moved to AliGen*lib.h
-- use AliGenGSIlib::kUpsilon, AliGenPHOSlib::kEtaPrime to access the constants
-
-Revision 1.3  2000/12/21 16:24:06  morsch
-Coding convention clean-up
-
-Revision 1.2  2000/11/30 07:12:50  alibrary
-Introducing new Rndm and QA classes
-
-Revision 1.1  2000/06/15 08:48:43  morsch
-AliGenGSIlib with parametersations for GSI physics simulation added (YF, MI)
-
-Revision 1.7  2000/05/02 08:12:13  morsch
-Coding rule violations corrected.
-
-Revision 1.6  1999/09/29 09:24:14  fca
-Introduction of the Copyright and cvs Log
-
-*/
+/* $Id$ */
 
 /////////////////////////////////////////////////////////////////////////////
 //                                                                         //
@@ -84,7 +54,7 @@ Int_t AliGenGSIlib::IpUpsilon(TRandom *)
   return 553;     
 
 }
-Double_t AliGenGSIlib::PtUpsilonFlat( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtUpsilonFlat( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -102,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -123,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( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -140,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -161,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( Double_t */*px*/, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -173,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(Double_t */*py*/, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -186,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( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -203,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -239,7 +209,7 @@ Int_t AliGenGSIlib::IpJpsi(TRandom *)
   return 443;     
 
 }
-Double_t AliGenGSIlib::PtJpsiFlat( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtJpsiFlat( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -257,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -278,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( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -294,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( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -311,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -336,7 +306,7 @@ Double_t AliGenGSIlib::YJpsiMUON(Double_t *py, Double_t *dummy)
 //                  J/Psi pt-distribution by Sergei
 //
 //--------------------------------------------------------------------------
-//Double_t AliGenGSIlib::PtJpsi( Double_t *px, Double_t *dummy )
+//Double_t AliGenGSIlib::PtJpsi( Double_t *px, Double_t */*dummy*/ )
 //{
 
 //  return x = gRandom->Rndm()*10.;
@@ -393,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -408,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -428,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -445,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -486,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -501,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -523,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -540,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( Double_t *px, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -568,7 +538,7 @@ Int_t AliGenGSIlib::IpEta(TRandom *)
   return 221;     
 
 }
-Double_t AliGenGSIlib::PtEtaPHOS( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtEtaPHOS( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -579,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -612,7 +582,7 @@ Int_t AliGenGSIlib::IpEtaprime(TRandom *)
   return 331;     
 
 }
-Double_t AliGenGSIlib::PtEtaprimePHOS( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtEtaprimePHOS( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -623,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -656,7 +626,7 @@ Int_t AliGenGSIlib::IpOmega(TRandom *)
   return 223;     
 
 }
-Double_t AliGenGSIlib::PtOmega( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtOmega( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -667,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -702,7 +672,7 @@ Int_t AliGenGSIlib::IpRho(TRandom *)
   return 113;     
 
 }
-Double_t AliGenGSIlib::PtRho( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtRho( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -713,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -762,7 +732,7 @@ Int_t AliGenGSIlib::IpPionPHOS(TRandom *ran)
       }
     }
 }
-Double_t AliGenGSIlib::PtPion( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtPion( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -796,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -845,7 +815,7 @@ Int_t AliGenGSIlib::IpKaonPHOS(TRandom *ran)
       }
     }
 }
-Double_t AliGenGSIlib::PtKaonPHOS( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtKaonPHOS( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -856,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -889,7 +859,7 @@ Int_t AliGenGSIlib::IpPhi(TRandom *)
   return 333;     
 
 }
-Double_t AliGenGSIlib::PtPhiPHOS( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtPhiPHOS( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -900,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -950,7 +920,7 @@ Int_t AliGenGSIlib::IpBaryons(TRandom *ran)
       }
     }
 }
-Double_t AliGenGSIlib::PtBaryons( Double_t *px, Double_t *dummy )
+Double_t AliGenGSIlib::PtBaryons( Double_t *px, Double_t */*dummy*/ )
 {
 //--------------------------------------------------------------------------
 //
@@ -961,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(Double_t *py, Double_t */*dummy*/)
 {
 //--------------------------------------------------------------------------
 //
@@ -1039,120 +1009,120 @@ GenFunc AliGenGSIlib::GetPt(Int_t param, const char * tname) const
     {
     case kUpsilon:
       if (sname=="FLAT"){
-       func= PtUpsilonFlat;
-       break;
+        func= PtUpsilonFlat;
+        break;
       }
       if (sname=="MUON"){
-       func= PtUpsilonMUON;
-       break;
+        func= PtUpsilonMUON;
+        break;
       }
       if (sname=="RITMAN"){
-       func=PtUpsilonRitman;
-       break;
+        func=PtUpsilonRitman;
+        break;
       }
       if (sname=="KAREL"){
-       func=PtUpsilonKarel;
-       break;
+        func=PtUpsilonKarel;
+        break;
       }
+      func=0;
+      printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
       break;
-       func=0;
-        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
 
     case kJPsi:
       if (sname=="FLAT"){
-       func= PtJpsiFlat;
-       break;
+        func= PtJpsiFlat;
+        break;
       }
       if (sname=="MUON"){
-       func= PtJpsiMUON;
-       break;
+        func= PtJpsiMUON;
+        break;
       }
       //      if (sname=="SERGEI"){
       //       func= PtJpsi;
       //       break;
       //     }
+      func=0;
+      printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
       break;
-        func=0;
-        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
 
     case kCharm: 
       if (sname=="FLAT"){
-       func= PtCharmFlat;
-       break;
+        func= PtCharmFlat;
+        break;
       }
 
       if (sname=="MUON"){
-       func= PtCharmMUON;
-       break;
+        func= PtCharmMUON;
+        break;
       }
 
       if (sname=="GSI"){
-       func= PtCharmGSI;
-       break;
+        func= PtCharmGSI;
+        break;
       }
+      func=0;
+      printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
       break;
-       func=0;
-        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
 
     case kBeauty: 
       if (sname=="FLAT"){
-       func= PtBeautyFlat;
-       break;
+        func= PtBeautyFlat;
+        break;
       }
       if (sname=="MUON"){
-       func= PtBeautyMUON;
-       break;
+        func= PtBeautyMUON;
+        break;
       }
       if (sname=="GSI"){
-       func= PtBeautyGSI;
-       break;
+        func= PtBeautyGSI;
+        break;
       }
+      func=0;
+      printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
       break;
-       func=0;
-        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
 
 
     case kEta:
-         func=PtEtaPHOS;
-         break;
+      func=PtEtaPHOS;
+      break;
 
     case kEtaprime:
-         func=PtEtaprimePHOS;
-         break;
+      func=PtEtaprimePHOS;
+      break;
 
     case kOmega:
-         func=PtOmega;
-         break;
+      func=PtOmega;
+      break;
 
     case kRho:
-        func=PtRho;
-        break;
+      func=PtRho;
+      break;
 
     case kKaon:
-         func=PtKaonPHOS;
-         break;
+      func=PtKaonPHOS;
+      break;
 
     case kPion:
-         func=PtPion;
-         break;
+      func=PtPion;
+      break;
 
     case kPhi:
-         func=PtPhiPHOS;
-         break;
+      func=PtPhiPHOS;
+      break;
 
         //    case kLambda:
         //         func=PtLambda;
         //         break;
 
     case kBaryons:
-         func=PtBaryons;
-         break;
+      func=PtBaryons;
+      break;
 
     default:
-        func=0;
-        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
+      func=0;
+      printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
     }
-    return func;
+   return func;
 }
 
 GenFunc AliGenGSIlib::GetY(Int_t param, const char * tname) const