/* $Id$ */
-#include "TMath.h"
-#include "TRandom.h"
+#include <TMath.h>
+#include <TRandom.h>
#include "AliGenSTRANGElib.h"
// 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
}
//=============================================================
//
- 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
// 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
//____________________________________________________________
- return PtScal(*px,2); // 2==> Kaon in the PtScal function
+ return PtScal(*px,1); // 1==> 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
//____________________________________________________________
const Double_t ka = 1000.;
- const Double_t kdy = 4.;
+ const Double_t kdy = 4.*4;
Double_t y=TMath::Abs(*py);
//============================================================================
//============================================================================
// 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
//____________________________________________________________
- return PtScal(*px,6); // 6==> Phi in the PtScal function
+ return PtScal(*px,5); // 5==> 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
//____________________________________________________________
const Double_t ka = 1000.;
- const Double_t kdy = 4.;
+ const Double_t kdy = 4.*4;
Double_t y=TMath::Abs(*py);
//===================================================================
//============================================================================
// Lambda
- Double_t AliGenSTRANGElib::PtLambda( Double_t *px, Double_t *)
+ Double_t AliGenSTRANGElib::PtLambda( const Double_t *px, const Double_t *)
{
// Lambda
// 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( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YLambda( const Double_t *py, const Double_t *)
{
// y-distribution
//____________________________________________________________
const Double_t ka = 1000.;
- const Double_t kdy = 4.;
+ const Double_t kdy = 4.*4;
Double_t y=TMath::Abs(*py);
// 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
//____________________________________________________________
- 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( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YXiMinus( const Double_t *py, const Double_t *)
{
// y-distribution
//____________________________________________________________
const Double_t ka = 1000.;
- const Double_t kdy = 4.;
+ const Double_t kdy = 4.*4;
Double_t y=TMath::Abs(*py);
// 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
//____________________________________________________________
- 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( Double_t *py, Double_t *)
+ Double_t AliGenSTRANGElib::YOmegaMinus( const Double_t *py, const Double_t *)
{
// y-distribution
//____________________________________________________________
const Double_t ka = 1000.;
- const Double_t kdy = 4.;
+ const Double_t kdy = 4.*4;
Double_t y=TMath::Abs(*py);
// 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
//____________________________________________________________
- 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( Double_t *py, Double_t *)
+Double_t AliGenSTRANGElib::YLambda1520( const Double_t *py, const Double_t *)
{
// y-distribution
//____________________________________________________________
// 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