Adding cluster energy distribution histograms

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

diff --git a/EVGEN/AliGenSTRANGElib.cxx b/EVGEN/AliGenSTRANGElib.cxx
index 572e027ae40648e1ee3c79211d5846aa94591fa2..bff4c4b6e8c3273d6b57657185c613a0fbcdcaa3 100644 (file)
--- a/EVGEN/AliGenSTRANGElib.cxx
+++ b/EVGEN/AliGenSTRANGElib.cxx
@@ -37,7 +37,6 @@
 /* $Id$ */
 
 #include <TMath.h>
-#include <TPDGCode.h>
 #include <TRandom.h>
 
 #include "AliGenSTRANGElib.h"
@@ -60,21 +59,20 @@ ClassImp(AliGenSTRANGElib)
 //     10=>Omega-antiOmega
 
   //    MASS SCALING RESPECT TO PIONS
-  //    MASS                1=>PI,  2=>K, 3=>ETA,4=>OMEGA,5=>ETA',6=>PHI 
+  //    MASS                0=>PI,  1=>K, 2=>ETA,3=>OMEGA,4=>ETA',5=>PHI 
   const Double_t khm[11] = {0.1396, 0.494,0.547, 0.782,   0.957,  1.02, 
-  //    MASS               7=>BARYON-BARYONBAR  
+  //    MASS               6=>BARYON-BARYONBAR  
                                  0.938, 
-  //    MASS               8=>Lambda-antiLambda  
+  //    MASS               7=>Lambda-antiLambda  
                                   1.1157,
-  //    MASS               9=>Xi-antiXi  
+  //    MASS               8=>Xi-antiXi  
                                   1.3213, 
-  //    MASS              10=>Omega-antiOmega  
+  //    MASS               9=>Omega-antiOmega  
                                   1.6725,
-  //    MASS              11=>Lambda(1520)
+  //    MASS              10=>Lambda(1520)
                                   1.5195};
   //     VALUE MESON/PI AT 5 GEV
   const Double_t kfmax[11]={1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.};
-  np--;
   Double_t f5=TMath::Power(((sqrt(100.018215)+2.)/(sqrt(100.+khm[np]*khm[np])+2.0)),12.3);
   Double_t kfmax2=f5/kfmax[np];
   // PIONS
@@ -121,7 +119,7 @@ ClassImp(AliGenSTRANGElib)
 //                pt-distribution
 //____________________________________________________________
 
-  return PtScal(*px,2);  //  2==> Kaon in the PtScal function
+  return PtScal(*px,1);  //  1==> Kaon in the PtScal function
 }
 
  Double_t AliGenSTRANGElib::YKaon( const Double_t *py, const Double_t *)
@@ -130,7 +128,7 @@ ClassImp(AliGenSTRANGElib)
 //____________________________________________________________
 
   const Double_t ka    = 1000.;
-  const Double_t kdy   = 4.;
+  const Double_t kdy   = 4.*4;
 
 
   Double_t y=TMath::Abs(*py);
@@ -173,7 +171,7 @@ ClassImp(AliGenSTRANGElib)
 //                pt-distribution
 //____________________________________________________________
 
-  return PtScal(*px,6);  //  6==> Phi in the PtScal function
+  return PtScal(*px,5);  //  5==> Phi in the PtScal function
 }
 
  Double_t AliGenSTRANGElib::YPhi( const Double_t *py, const Double_t *)
@@ -182,7 +180,7 @@ ClassImp(AliGenSTRANGElib)
 //____________________________________________________________
 
   const Double_t ka    = 1000.;
-  const Double_t kdy   = 4.;
+  const Double_t kdy   = 4.*4;
 
 
   Double_t y=TMath::Abs(*py);
@@ -208,7 +206,7 @@ ClassImp(AliGenSTRANGElib)
 //                pt-distribution
 //____________________________________________________________
 
-  return PtScal(*px,8);  //  8==> Lambda-antiLambda in the PtScal function
+  return PtScal(*px,7);  //  7==> Lambda-antiLambda in the PtScal function
 }
 
  Double_t AliGenSTRANGElib::YLambda( const Double_t *py, const Double_t *)
@@ -217,7 +215,7 @@ ClassImp(AliGenSTRANGElib)
 //____________________________________________________________
 
   const Double_t ka    = 1000.;
-  const Double_t kdy   = 4.;
+  const Double_t kdy   = 4.*4;
 
 
   Double_t y=TMath::Abs(*py);
@@ -247,7 +245,7 @@ ClassImp(AliGenSTRANGElib)
 //                pt-distribution
 //____________________________________________________________
 
-  return PtScal(*px,9);  //  9==> Xi-antiXi in the PtScal function
+  return PtScal(*px,8);  //  8==> Xi-antiXi in the PtScal function
 }
 
  Double_t AliGenSTRANGElib::YXiMinus( const Double_t *py, const Double_t *)
@@ -256,7 +254,7 @@ ClassImp(AliGenSTRANGElib)
 //____________________________________________________________
 
   const Double_t ka    = 1000.;
-  const Double_t kdy   = 4.;
+  const Double_t kdy   = 4.*4;
 
 
   Double_t y=TMath::Abs(*py);
@@ -286,7 +284,7 @@ ClassImp(AliGenSTRANGElib)
 //                pt-distribution
 //____________________________________________________________
 
-  return PtScal(*px,10);  //  10==> Omega-antiOmega in the PtScal function
+  return PtScal(*px,9);  //  9==> Omega-antiOmega in the PtScal function
 }
 
  Double_t AliGenSTRANGElib::YOmegaMinus( const Double_t *py, const Double_t *)
@@ -295,7 +293,7 @@ ClassImp(AliGenSTRANGElib)
 //____________________________________________________________
 
   const Double_t ka    = 1000.;
-  const Double_t kdy   = 4.;
+  const Double_t kdy   = 4.*4;
 
 
   Double_t y=TMath::Abs(*py);
@@ -326,7 +324,7 @@ Double_t AliGenSTRANGElib::PtLambda1520( const Double_t *px, const Double_t *)
 //                  pt-distribution
 //____________________________________________________________
 
-  return PtScal(*px,11);   //   11=> Lambda(1520) in the PtScal function
+  return PtScal(*px,10);   //   10=> Lambda(1520) in the PtScal function
 }
 
 Double_t AliGenSTRANGElib::YLambda1520( const Double_t *py, const Double_t *)