#include "AliGenMUONlib.h"
#include "AliMC.h"
#include "AliRun.h"
+#include "AliPDG.h"
+
#include <TDirectory.h>
+#include <TDatabasePDG.h>
#include <TFile.h>
#include <TTree.h>
#include <stdlib.h>
//____________________________________________________________
void AliGenHalo::Generate()
{
- AliMC* pMC = AliMC::GetMC();
FILE *fp = fopen(fFileName,"r");
if (fp) {
printf("\n File %s opened for reading ! \n ", fFileName);
}
//
// MARS particle codes
- const Int_t imars[12]={0,14, 13, 8, 9, 11, 12, 5, 6, 1, 3, 2};
+ // const Int_t imars[12]={0,14, 13, 8, 9, 11, 12, 5, 6, 1, 3, 2};
+ const Int_t imars[12]={0,kProton,kNeutron,kPiPlus,kPiMinus,kKPlus,kKMinus,
+ kMuonPlus,kMuonMinus,kGamma,kElectron,kPositron};
Float_t polar[3]= {0,0,0};
Float_t origin[3];
Float_t p[3], p0;
Float_t ekin, wgt, tx, ty, tz, txy;
- Float_t amass, charge, tlife;
- char name[100];
+ Float_t amass;
//
- Int_t ipart, itrtyp, ncols, nt;
+ Int_t ipart, ncols, nt;
Int_t nread=0;
origin[2]=2650;
nread++;
if (fNpart !=-1 && nread > fNpart) break;
ipart = imars[ipart];
- pMC->Gfpart(ipart, name, itrtyp,
- amass, charge, tlife);
+ amass = TDatabasePDG::Instance()->GetParticle(ipart)->Mass();
p0=sqrt(ekin*ekin + 2.*amass);
txy=TMath::Sqrt(tx*tx+ty*ty);