All coding rule violations except RS3 corrected
[u/mrichter/AliRoot.git] / EVGEN / AliGenScan.cxx
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$
18 Revision 1.4  1999/11/03 17:43:20  fca
19 New version from G.Martinez & A.Morsch
20
21 Revision 1.3  1999/09/29 09:24:14  fca
22 Introduction of the Copyright and cvs Log
23
24 */
25
26 #include "AliGenScan.h"
27 #include "AliRun.h"
28
29  ClassImp(AliGenScan)
30     
31  AliGenScan::AliGenScan()
32          :AliGenerator(-1)
33 {
34 // Constructor
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
49 AliGenScan::AliGenScan(Int_t npart)
50     :AliGenerator(npart)
51 {
52 // Constructor
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 //____________________________________________________________
65 AliGenScan::~AliGenScan()
66 {
67 // Destructor
68 }
69
70 void 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 {
74 // Define the grid
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 //____________________________________________________________
87 void AliGenScan::Generate()
88 {
89   //
90   // Generate one trigger
91   //
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++){
124               gMC->Rndm(random,6);
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);
134               gAlice->SetTrack(fTrackIt,-1,fIpart,p,origin,polar,0,"Primary",nt);
135           }
136       }
137   }
138 }
139
140
141
142
143
144
145
146
147