]>
Commit | Line | Data |
---|---|---|
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 | /* $Id$ */ | |
17 | ||
18 | #include "TMevSimPartTypeParams.h" | |
19 | ||
20 | ClassImp(TMevSimPartTypeParams) | |
21 | ||
22 | ||
23 | //______________________________________________________________________________ | |
24 | TMevSimPartTypeParams::TMevSimPartTypeParams() | |
25 | { | |
26 | // Default constructor | |
27 | ||
28 | for (Int_t i = 0; i < 4; i++) | |
29 | for (Int_t j = 0; j < NFLOWTERMS; j++) { | |
30 | fVnMean[j][i] = 0.0; | |
31 | fVnStDev[j][i] = 0.0; | |
32 | } | |
33 | } | |
34 | ||
35 | //__________________________________________________________________________ | |
36 | TMevSimPartTypeParams::TMevSimPartTypeParams(Int_t agpid, Int_t amultmean, Int_t amultvc, | |
37 | Float_t atempmean, Float_t atempstdev, Float_t asigmamean, | |
38 | Float_t asigmastdev, Float_t aexpvelmean, Float_t aexpvelstdev) | |
39 | { | |
40 | // Construct the particle type parametrs class. Use the values provide | |
41 | // by the user to initialize the internal variables. For the meaning of | |
42 | // the parametrs see the TMevSim class documentation. | |
43 | ||
44 | fGPid = agpid; | |
45 | fMultMean = amultmean; | |
46 | fMultVarianceControl = amultvc; | |
47 | fTempMean = atempmean; | |
48 | fTempStDev = atempstdev; | |
49 | fSigmaMean = asigmamean; | |
50 | fSigmaStDev = asigmastdev; | |
51 | fExpVelMean = aexpvelmean; | |
52 | fExpVelStDev = aexpvelstdev; | |
53 | for (Int_t i = 0; i < 4; i++) | |
54 | for (Int_t j = 0; j < NFLOWTERMS; j++) { | |
55 | fVnMean[j][i] = 0.0; | |
56 | fVnStDev[j][i] = 0.0; | |
57 | } | |
58 | ||
59 | } | |
60 | ||
61 | //______________________________________________________________________________ | |
62 | TMevSimPartTypeParams::~TMevSimPartTypeParams() | |
63 | { | |
64 | //dtor | |
65 | } | |
66 | ||
67 | //______________________________________________________________________________ | |
68 | TMevSimPartTypeParams::TMevSimPartTypeParams (const TMevSimPartTypeParams& pars) : TObject(pars) { | |
69 | // The copy constructor | |
70 | ||
71 | this->fGPid = pars.GetGPid(); | |
72 | this->fMultMean = pars.GetMultMean(); | |
73 | this->fMultVarianceControl = pars.GetMultVarianceControl(); | |
74 | this->fTempMean = pars.GetTempMean(); | |
75 | this->fTempStDev = pars.GetTempStDev(); | |
76 | this->fSigmaMean = pars.GetSigmaMean(); | |
77 | this->fSigmaStDev = pars.GetSigmaStDev(); | |
78 | this->fExpVelMean = pars.GetExpVelMean(); | |
79 | this->fExpVelStDev = pars.GetExpVelStDev(); | |
80 | for (Int_t i = 0; i < 4; i++) | |
81 | for (Int_t j = 0; j < NFLOWTERMS; j++) { | |
82 | this->fVnMean[j][i] = pars.GetVnMeanComponent(j, i); | |
83 | this->fVnStDev[j][i] = pars.GetVnStDevComponent(j, i); | |
84 | } | |
85 | } | |
86 | ||
87 | //______________________________________________________________________________ | |
88 | TMevSimPartTypeParams& TMevSimPartTypeParams::operator=(const TMevSimPartTypeParams& pars) { | |
89 | // The assignment operator | |
90 | ||
91 | this->fGPid = pars.GetGPid(); | |
92 | this->fMultMean = pars.GetMultMean(); | |
93 | this->fMultVarianceControl = pars.GetMultVarianceControl(); | |
94 | this->fTempMean = pars.GetTempMean(); | |
95 | this->fTempStDev = pars.GetTempStDev(); | |
96 | this->fSigmaMean = pars.GetSigmaMean(); | |
97 | this->fSigmaStDev = pars.GetSigmaStDev(); | |
98 | this->fExpVelMean = pars.GetExpVelMean(); | |
99 | this->fExpVelStDev = pars.GetExpVelStDev(); | |
100 | for (Int_t i = 0; i < 4; i++) | |
101 | for (Int_t j = 0; j < NFLOWTERMS; j++) { | |
102 | this->fVnMean[j][i] = GetVnMeanComponent(j, i); | |
103 | this->fVnStDev[j][i] = GetVnStDevComponent(j, i); | |
104 | } | |
105 | return (*this); | |
106 | } | |
107 | ||
108 | //______________________________________________________________________________ | |
109 | void TMevSimPartTypeParams::SetGPid(Int_t gpid) { | |
110 | //Sets pid | |
111 | fGPid = gpid; | |
112 | } | |
113 | //______________________________________________________________________________ | |
114 | Int_t TMevSimPartTypeParams::GetGPid() const { | |
115 | //returns pid | |
116 | return fGPid; | |
117 | } | |
118 | //______________________________________________________________________________ | |
119 | void TMevSimPartTypeParams::SetMultMean(Int_t multmean) { | |
120 | //sets multiplicity mean | |
121 | fMultMean = multmean; | |
122 | } | |
123 | //______________________________________________________________________________ | |
124 | Int_t TMevSimPartTypeParams::GetMultMean() const { | |
125 | //returns multiplicity mean | |
126 | return fMultMean; | |
127 | } | |
128 | //______________________________________________________________________________ | |
129 | void TMevSimPartTypeParams::SetMultVarianceControl(Int_t multvc) { | |
130 | //sets multiplicity variance | |
131 | fMultVarianceControl = multvc; | |
132 | } | |
133 | //______________________________________________________________________________ | |
134 | Int_t TMevSimPartTypeParams::GetMultVarianceControl() const { | |
135 | //returns multiplicity variance | |
136 | return fMultVarianceControl; | |
137 | } | |
138 | //______________________________________________________________________________ | |
139 | void TMevSimPartTypeParams::SetTempParams(Float_t tempmean, Float_t tempsdev) { | |
140 | //Sets mean and variance of temperature | |
141 | fTempMean = tempmean; | |
142 | fTempStDev = tempsdev; | |
143 | } | |
144 | //______________________________________________________________________________ | |
145 | Float_t TMevSimPartTypeParams::GetTempMean() const { | |
146 | //returns mean temperature | |
147 | return fTempMean; | |
148 | } | |
149 | //______________________________________________________________________________ | |
150 | Float_t TMevSimPartTypeParams::GetTempStDev() const { | |
151 | //returns temperature variance | |
152 | return fTempStDev; | |
153 | } | |
154 | //______________________________________________________________________________ | |
155 | void TMevSimPartTypeParams::SetSigmaPrams(Float_t sigmamean, Float_t sigmastdev) { | |
156 | //Sets mean and variance of rapidity | |
157 | fSigmaMean = sigmamean; | |
158 | fSigmaStDev = sigmastdev; | |
159 | } | |
160 | //______________________________________________________________________________ | |
161 | Float_t TMevSimPartTypeParams::GetSigmaMean() const { | |
162 | //returns mean rapidity | |
163 | return fSigmaMean; | |
164 | } | |
165 | //______________________________________________________________________________ | |
166 | Float_t TMevSimPartTypeParams::GetSigmaStDev() const { | |
167 | //returns rapidity variance | |
168 | return fSigmaStDev; | |
169 | } | |
170 | //______________________________________________________________________________ | |
171 | void TMevSimPartTypeParams::SetExpVelParams(Float_t expvelmean, Float_t expvelstdev) { | |
172 | //Sets mean and variance of Expansion velocity ala Scott Pratt (in units of c) | |
173 | fExpVelMean = expvelmean; | |
174 | fExpVelStDev = expvelstdev; | |
175 | } | |
176 | //______________________________________________________________________________ | |
177 | Float_t TMevSimPartTypeParams::GetExpVelMean() const { | |
178 | //Sets mean of Expansion velocity ala Scott Pratt (in units of c) | |
179 | return fExpVelMean; | |
180 | } | |
181 | //______________________________________________________________________________ | |
182 | Float_t TMevSimPartTypeParams::GetExpVelStDev() const { | |
183 | //Sets variance of Expansion velocity ala Scott Pratt (in units of c) | |
184 | return fExpVelStDev; | |
185 | } | |
186 | //______________________________________________________________________________ | |
187 | ||
188 | void TMevSimPartTypeParams::SetVnMeanComponent | |
189 | (Int_t nComponent, Float_t mean1, Float_t mean2, Float_t mean3, Float_t mean4) { | |
190 | //Sets VnMean: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6 | |
191 | // mean include all 6 sets of parameters, set them to 0 if not used | |
192 | ||
193 | if (nComponent < 0 || nComponent > NFLOWTERMS ) | |
194 | Error("SetVnMeanComponent", "Wrong Vn component n = %d (must be [%d , %d])", | |
195 | nComponent, 0, NFLOWTERMS); | |
196 | ||
197 | fVnMean[nComponent][0] = mean1; | |
198 | fVnMean[nComponent][1] = mean2; | |
199 | fVnMean[nComponent][2] = mean3; | |
200 | fVnMean[nComponent][3] = mean4; | |
201 | } | |
202 | //______________________________________________________________________________ | |
203 | ||
204 | void TMevSimPartTypeParams::SetVnStDevComponent | |
205 | (Int_t nComponent, Float_t stdev1, Float_t stdev2,Float_t stdev3, Float_t stdev4) { | |
206 | //Sets VnStDev: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6 | |
207 | // standard deviation include all 6 sets of parameters, set them to 0 if not used | |
208 | ||
209 | if (nComponent < 0 || nComponent > NFLOWTERMS ) | |
210 | Error("SetVnStDevComponent", "Wrong Vn component n = %d (must be [%d , %d])", | |
211 | nComponent, 0, NFLOWTERMS); | |
212 | ||
213 | fVnStDev[nComponent][0] = stdev1; | |
214 | fVnStDev[nComponent][1] = stdev2; | |
215 | fVnStDev[nComponent][2] = stdev3; | |
216 | fVnStDev[nComponent][3] = stdev4; | |
217 | } | |
218 | ||
219 | //______________________________________________________________________________ | |
220 | ||
221 | void TMevSimPartTypeParams::SetVnMeanComponent (Int_t nComponent, Float_t mean[4]) { | |
222 | //Sets VnMean: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6 | |
223 | // mean include all 6 sets of parameters, set them to 0 if not used | |
224 | if (nComponent < 0 || nComponent > NFLOWTERMS ) | |
225 | Error("SetVnMeanComponent", "Wrong Vn component (%d) (must be [%d,%d])", | |
226 | nComponent, 0, NFLOWTERMS); | |
227 | ||
228 | Int_t n = 4; | |
229 | for (Int_t i=0; i<n; i++) | |
230 | fVnMean[nComponent][i] = mean[i]; | |
231 | ||
232 | } | |
233 | ||
234 | //______________________________________________________________________________ | |
235 | ||
236 | void TMevSimPartTypeParams::SetVnStDevComponent (Int_t nComponent, Float_t stdev[4]) { | |
237 | //Sets VnStDev: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6 | |
238 | // standard deviation include all 6 sets of parameters, set them to 0 if not used | |
239 | ||
240 | if (nComponent < 0 || nComponent > NFLOWTERMS ) | |
241 | Error("SetVnStDevComponent", "Wrong Vn component n = %d (must be [%d , %d])", | |
242 | nComponent, 0, NFLOWTERMS); | |
243 | ||
244 | Int_t n = 4; | |
245 | for (Int_t i=0; i<n; i++) | |
246 | fVnStDev[nComponent][i] = stdev[i]; | |
247 | } | |
248 | ||
249 | //______________________________________________________________________________ | |
250 | ||
251 | Float_t TMevSimPartTypeParams::GetVnMeanComponent(Int_t nComponent, Int_t nMean) const { | |
252 | //Returns VnMean: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6 | |
253 | // mean include all 6 sets of parameters, set them to 0 if not used | |
254 | ||
255 | if ((nComponent < 0) || (nComponent>NFLOWTERMS) || (nMean < 0) || (nMean > 3)) return 0.0; | |
256 | return fVnMean[nComponent][nMean]; | |
257 | } | |
258 | //______________________________________________________________________________ | |
259 | Float_t TMevSimPartTypeParams::GetVnStDevComponent(Int_t nComponent, Int_t nStDev) const { | |
260 | //Returns VnStDev: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6 | |
261 | // standard deviation include all 6 sets of parameters, set them to 0 if not used | |
262 | ||
263 | if ((nComponent < 0) || (nComponent>NFLOWTERMS) || (nStDev < 0) || (nStDev > 3)) return 0.0; | |
264 | return fVnMean[nComponent][nStDev]; | |
265 | } | |
266 | ||
267 | //______________________________________________________________________________ | |
268 | ||
269 | ||
270 | ||
271 | ||
272 |