1 #include "AliGenScan.h"
6 AliGenScan::AliGenScan()
24 AliGenScan::AliGenScan(Int_t npart)
38 //____________________________________________________________
39 AliGenScan::~AliGenScan()
42 void AliGenScan::SetRange(Int_t nx, Float_t xmin, Float_t xmax,
43 Int_t ny, Float_t ymin, Float_t ymax,
44 Int_t nz, Float_t zmin, Float_t zmax)
57 //____________________________________________________________
58 void AliGenScan::Generate()
61 // Generate one trigger
63 AliMC* pMC = AliMC::GetMC();
65 Float_t polar[3]= {0,0,0};
70 Float_t pmom, theta, phi;
93 for (Int_t ix=0; ix<fNx; ix++) {
94 for (Int_t iy=0; iy<fNy; iy++) {
95 for (Int_t iz=0; iz<fNz; iz++){
97 origin[0]=fXmin+ix*dx+2*(random[0]-0.5)*fOsigma[0];
98 origin[1]=fYmin+iy*dy+2*(random[1]-0.5)*fOsigma[1];
99 origin[2]=fZmin+iz*dz+2*(random[2]-0.5)*fOsigma[2];
100 pmom=fPMin+random[3]*(fPMax-fPMin);
101 theta=fThetaMin+random[4]*(fThetaMax-fThetaMin);
102 phi=fPhiMin+random[5]*(fPhiMax-fPhiMin);
103 p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
104 p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
105 p[2] = pmom*TMath::Cos(theta);
106 gAlice->SetTrack(1,-1,fIpart,p,origin,polar,0,"Primary",nt);