/* $Id$ */
+//======================================================================
+// AliGenPHOSlib class contains parameterizations of the
+// pion, kaon, eta, omega, etaprime, phi and baryon (proton,
+// antiproton, neutron and anti-neutron) particles for the
+// study of the neutral background in PHOS detector.
+// Additional particle species simulation options has been added:
+// Charged Pion, Charged Kaons, KLong Proton, Anti-Proton, Neutron,
+// Anti-Neutron --> Changes made by Gustavo Conesa in November 2004
+//======================================================================
+
class TRandom;
class AliGenPHOSlib :
public AliGenLib
{
public:
- enum constants{kPion, kKaon, kEta, kOmega, kEtaPrime, kPhi, kBaryon};
+ enum constants{kPion, kChargedPion, kPi0Flat, kKaon, kChargedKaon, kKaon0L,
+ kEta, kEtaFlat,kOmega, kEtaPrime, kPhi,
+ kBaryon, kProton, kAProton, kNeutron, kANeutron};
+ GenFunc GetPt(Int_t param, const char* tname=0) const;
+ GenFunc GetY (Int_t param, const char* tname=0) const;
+ 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 PtScal(Double_t pt, Int_t np);
static Double_t YPion( Double_t *py, 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 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 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 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 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 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 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 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 Int_t IpBaryon(TRandom *ran);
-
- GenFunc GetPt(Int_t param, const char* tname=0);
- GenFunc GetY (Int_t param, const char* tname=0);
- GenFuncIp GetIp(Int_t param, const char* tname=0);
+ static Int_t IpProton(TRandom *ran);
+ static Int_t IpAProton(TRandom *ran);
+ static Int_t IpNeutron(TRandom *ran);
+ static Int_t IpANeutron(TRandom *ran);
+
ClassDef(AliGenPHOSlib,0) // Library providing y and pT parameterisations
};
#endif