Update of SSD simulation and reconstruction code by Boris and Enrico.
[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$
38f1bd58 18Revision 1.7 2000/11/30 07:12:50 alibrary
19Introducing new Rndm and QA classes
20
65fb704d 21Revision 1.6 2000/10/02 21:28:06 fca
22Removal of useless dependecies via forward declarations
23
94de3818 24Revision 1.5 2000/06/09 20:37:20 morsch
25All coding rule violations except RS3 corrected
26
f87cfe57 27Revision 1.4 1999/11/03 17:43:20 fca
28New version from G.Martinez & A.Morsch
29
886b6f73 30Revision 1.3 1999/09/29 09:24:14 fca
31Introduction of the Copyright and cvs Log
32
4c039060 33*/
34
b0418df4 35#include "AliGenScan.h"
b0418df4 36#include "AliRun.h"
f87cfe57 37
b0418df4 38 ClassImp(AliGenScan)
39
40 AliGenScan::AliGenScan()
41 :AliGenerator(-1)
42{
f87cfe57 43// Constructor
38f1bd58 44 fXCmin=0;
45 fXCmax=0;
b0418df4 46 fNx=1;
38f1bd58 47 fYCmin=0;
48 fYCmax=0;
b0418df4 49 fNy=1;
50 fZmin=0;
51 fZmax=0;
52 fNz=1;
53//
54// Read all particles
55 fNpart=-1;
56}
57
58AliGenScan::AliGenScan(Int_t npart)
59 :AliGenerator(npart)
60{
f87cfe57 61// Constructor
38f1bd58 62 fXCmin=0;
63 fXCmax=0;
b0418df4 64 fNx=1;
38f1bd58 65 fYCmin=0;
66 fYCmax=0;
b0418df4 67 fNy=1;
68 fZmin=0;
69 fZmax=0;
70 fNz=1;
71}
72
73//____________________________________________________________
74AliGenScan::~AliGenScan()
f87cfe57 75{
76// Destructor
77}
b0418df4 78
79void AliGenScan::SetRange(Int_t nx, Float_t xmin, Float_t xmax,
80 Int_t ny, Float_t ymin, Float_t ymax,
81 Int_t nz, Float_t zmin, Float_t zmax)
82{
f87cfe57 83// Define the grid
38f1bd58 84 fXCmin=xmin;
85 fXCmax=xmax;
b0418df4 86 fNx=nx;
38f1bd58 87 fYCmin=ymin;
88 fYCmax=ymax;
b0418df4 89 fNy=ny;
90 fZmin=zmin;
91 fZmax=zmax;
92 fNz=nz;
93}
94
95//____________________________________________________________
96void AliGenScan::Generate()
97{
98 //
99 // Generate one trigger
100 //
b0418df4 101
102 Float_t polar[3]= {0,0,0};
103 //
104 Float_t origin[3];
105 Float_t p[3];
106 Int_t nt;
107 Float_t pmom, theta, phi;
108 //
109 Float_t random[6];
110 Float_t dx,dy,dz;
111
112 //
113 if (fNy > 0) {
38f1bd58 114 dx=(fXCmax-fXCmin)/fNx;
b0418df4 115 } else {
116 dx=1e10;
117 }
118
119 if (fNy > 0) {
38f1bd58 120 dy=(fYCmax-fYCmin)/fNy;
b0418df4 121 } else {
122 dy=1e10;
123 }
124
125 if (fNz > 0) {
126 dz=(fZmax-fZmin)/fNz;
127 } else {
128 dz=1e10;
129 }
130 for (Int_t ix=0; ix<fNx; ix++) {
131 for (Int_t iy=0; iy<fNy; iy++) {
132 for (Int_t iz=0; iz<fNz; iz++){
65fb704d 133 Rndm(random,6);
38f1bd58 134 origin[0]=fXCmin+ix*dx+2*(random[0]-0.5)*fOsigma[0];
135 origin[1]=fYCmin+iy*dy+2*(random[1]-0.5)*fOsigma[1];
b0418df4 136 origin[2]=fZmin+iz*dz+2*(random[2]-0.5)*fOsigma[2];
137 pmom=fPMin+random[3]*(fPMax-fPMin);
138 theta=fThetaMin+random[4]*(fThetaMax-fThetaMin);
139 phi=fPhiMin+random[5]*(fPhiMax-fPhiMin);
140 p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
141 p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
142 p[2] = pmom*TMath::Cos(theta);
65fb704d 143 gAlice->SetTrack(fTrackIt,-1,fIpart,p,origin,polar,0,kPPrimary,nt);
b0418df4 144 }
145 }
146 }
147}
148
149
150
151
152
153
154
155
156