#include <TObjArray.h>
#include <TPDGCode.h>
#include <TParticle.h>
+#include <TDatabasePDG.h>
#include <TROOT.h>
fPhiFormula(0),
fCurrentForm(0),
fPtYHist(0),
- fPartTypes(0) {
+ fPartTypes(0)
+{
//
// Default constructor
//
//////////////////////////////////////////////////////////////////////////////////
-AliGenGeVSim::AliGenGeVSim(Float_t psi, Bool_t isMultTotal) : AliGenerator(-1) {
+AliGenGeVSim::AliGenGeVSim(Float_t psi, Bool_t isMultTotal)
+ : AliGenerator(-1),
+ fModel(0),
+ fPsi(psi),
+ fIsMultTotal(isMultTotal),
+ fPtFormula(0),
+ fYFormula(0),
+ fPhiFormula(0),
+ fCurrentForm(0),
+ fPtYHist(0),
+ fPartTypes(0)
+ {
//
// Standard Constructor.
//
// checking consistancy
if (psi < 0 || psi > 360 )
- Error ("AliGenGeVSim", "Reaction plane angle ( %d )out of range [0..360]", psi);
+ Error ("AliGenGeVSim", "Reaction plane angle ( %13.3f )out of range [0..360]", psi);
fPsi = psi * TMath::Pi() / 180. ;
fIsMultTotal = isMultTotal;
- // initialization
+ // Initialization
fPartTypes = new TObjArray();
InitFormula();
// mass -> [0] , temperature -> [1] , expansion velocity -> [2]
Double_t aFormE = TMath::Sqrt(par[0]*par[0] + x[0]*x[0]) * TMath::CosH(x[1]);
- Double_t aFormG = 1 / TMath::Sqrt( 1 - par[2]*par[2] );
+ Double_t aFormG = 1 / TMath::Sqrt((1.-par[2])*(1.+par[2]));
Double_t aFormYp = par[2]*TMath::Sqrt( (par[0]*par[0] + x[0]*x[0])
- * TMath::CosH(x[1])*TMath::CosH(x[1])
- - par[0]*par[0] )
- /( par[1]*TMath::Sqrt(1-par[2]*par[2]));
+ * (TMath::CosH(x[1])-par[0])*(TMath::CosH(x[1])+par[0]))
+ /( par[1]*TMath::Sqrt((1.-par[2])*(1.+par[2])));
return x[0] * aFormE * TMath::Exp( - aFormG * aFormE / par[1])
*( TMath::SinH(aFormYp)/aFormYp