New code from Piergiorgio added
[u/mrichter/AliRoot.git] / EVGEN / AliGenScan.cxx
CommitLineData
4c039060 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/*
17$Log$
f87cfe57 18Revision 1.4 1999/11/03 17:43:20 fca
19New version from G.Martinez & A.Morsch
20
886b6f73 21Revision 1.3 1999/09/29 09:24:14 fca
22Introduction of the Copyright and cvs Log
23
4c039060 24*/
25
b0418df4 26#include "AliGenScan.h"
b0418df4 27#include "AliRun.h"
f87cfe57 28
b0418df4 29 ClassImp(AliGenScan)
30
31 AliGenScan::AliGenScan()
32 :AliGenerator(-1)
33{
f87cfe57 34// Constructor
b0418df4 35 fXmin=0;
36 fXmax=0;
37 fNx=1;
38 fYmin=0;
39 fYmax=0;
40 fNy=1;
41 fZmin=0;
42 fZmax=0;
43 fNz=1;
44//
45// Read all particles
46 fNpart=-1;
47}
48
49AliGenScan::AliGenScan(Int_t npart)
50 :AliGenerator(npart)
51{
f87cfe57 52// Constructor
b0418df4 53 fXmin=0;
54 fXmax=0;
55 fNx=1;
56 fYmin=0;
57 fYmax=0;
58 fNy=1;
59 fZmin=0;
60 fZmax=0;
61 fNz=1;
62}
63
64//____________________________________________________________
65AliGenScan::~AliGenScan()
f87cfe57 66{
67// Destructor
68}
b0418df4 69
70void AliGenScan::SetRange(Int_t nx, Float_t xmin, Float_t xmax,
71 Int_t ny, Float_t ymin, Float_t ymax,
72 Int_t nz, Float_t zmin, Float_t zmax)
73{
f87cfe57 74// Define the grid
b0418df4 75 fXmin=xmin;
76 fXmax=xmax;
77 fNx=nx;
78 fYmin=ymin;
79 fYmax=ymax;
80 fNy=ny;
81 fZmin=zmin;
82 fZmax=zmax;
83 fNz=nz;
84}
85
86//____________________________________________________________
87void AliGenScan::Generate()
88{
89 //
90 // Generate one trigger
91 //
b0418df4 92
93 Float_t polar[3]= {0,0,0};
94 //
95 Float_t origin[3];
96 Float_t p[3];
97 Int_t nt;
98 Float_t pmom, theta, phi;
99 //
100 Float_t random[6];
101 Float_t dx,dy,dz;
102
103 //
104 if (fNy > 0) {
105 dx=(fXmax-fXmin)/fNx;
106 } else {
107 dx=1e10;
108 }
109
110 if (fNy > 0) {
111 dy=(fYmax-fYmin)/fNy;
112 } else {
113 dy=1e10;
114 }
115
116 if (fNz > 0) {
117 dz=(fZmax-fZmin)/fNz;
118 } else {
119 dz=1e10;
120 }
121 for (Int_t ix=0; ix<fNx; ix++) {
122 for (Int_t iy=0; iy<fNy; iy++) {
123 for (Int_t iz=0; iz<fNz; iz++){
cfce8870 124 gMC->Rndm(random,6);
b0418df4 125 origin[0]=fXmin+ix*dx+2*(random[0]-0.5)*fOsigma[0];
126 origin[1]=fYmin+iy*dy+2*(random[1]-0.5)*fOsigma[1];
127 origin[2]=fZmin+iz*dz+2*(random[2]-0.5)*fOsigma[2];
128 pmom=fPMin+random[3]*(fPMax-fPMin);
129 theta=fThetaMin+random[4]*(fThetaMax-fThetaMin);
130 phi=fPhiMin+random[5]*(fPhiMax-fPhiMin);
131 p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
132 p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
133 p[2] = pmom*TMath::Cos(theta);
886b6f73 134 gAlice->SetTrack(fTrackIt,-1,fIpart,p,origin,polar,0,"Primary",nt);
b0418df4 135 }
136 }
137 }
138}
139
140
141
142
143
144
145
146
147