* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.1 2000/09/08 15:37:58 morsch
-pt and y-parameterisations for PMD physics simulation.
-THIS VERSION IS JUST A TEMPLATE THE PT AND Y DISTRIBUTIONS HAVE TO BE CORRECTED !!!!.
+/* $Id$ */
-*/
+//
+// Library of generators for PMD
+// providing y and pt parametrisation
+// for generated tracks
+// Specific for PMD needs
+// Author: PMD Offline Group
+//
+
+#include <TMath.h>
+#include <TPDGCode.h>
#include "AliGenPMDlib.h"
-#include "AliPDG.h"
-#include "TMath.h"
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
//
// 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
const Double_t khm[10] = {.13957,.493,.5488,.769,.7826,.958,1.02,0,0,0};
// VALUE MESON/PI AT 5 GEV
const Double_t kfmax[10]={1.,0.3,0.55,1.0,1.0,1.0,1.0,0,0,0};
- np--;
Double_t f5=TMath::Power(((sqrt(100.018215)+2.)/(sqrt(100.+khm[np]*khm[np])+2.0)),12.3);
Double_t fmax2=f5/kfmax[np];
// PIONS
// 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);
+ return PtScal(*px, 1);
}
// 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
}
-typedef Double_t (*GenFunc) (Double_t*, Double_t*);
-GenFunc AliGenPMDlib::GetPt(Param_t param, const char* tname)
+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
GenFunc func=NULL;
switch (param)
{
- case Pion:
+ case kPion:
func=PtPi0;
break;
- case Eta:
+ case kEta:
func=PtEta;
break;
default:
return func;
}
-GenFunc AliGenPMDlib::GetY(Param_t param, const char* tname)
+GenFunc AliGenPMDlib::GetY(Int_t param, const char* /*tname*/) const
{
// Return pointer to y- parameterisation
GenFunc func=NULL;
switch (param)
{
- case Pion:
+ case kPion:
func=YPi0;
break;
- case Eta:
+ case kEta:
func=YEta;
break;
default:
}
typedef Int_t (*GenFuncIp) (TRandom *);
-GenFuncIp AliGenPMDlib::GetIp(Param_t param, const char* tname)
+GenFuncIp AliGenPMDlib::GetIp(Int_t param, const char* /*tname*/) const
{
// Return pointer to particle type parameterisation
GenFuncIp func=NULL;
switch (param)
{
- case Pion:
+ case kPion:
func=IpPi0;
break;
- case Eta:
+ case kEta:
func=IpEta;
break;
default: