]>
Commit | Line | Data |
---|---|---|
c5f0f3c1 | 1 | /************************************************************************** |
0fc5cc25 | 2 | * Copyright(c) 1998-2003, ALICE Experiment at CERN, All rights reserved. * |
c5f0f3c1 | 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 | // Implementation of the Primary Vertex class | |
af7ba10c | 18 | // for the Event Data Summary class |
19 | // This class contains the Primary Vertex | |
20 | // of the event coming from reconstruction | |
c5f0f3c1 | 21 | // Origin: A.Dainese, Padova, andrea.dainese@pd.infn.it |
22 | //----------------------------------------------------------------- | |
23 | ||
24 | //---- standard headers ---- | |
8ba665e9 | 25 | #include "Riostream.h" |
c5f0f3c1 | 26 | //---- Root headers -------- |
504de69b | 27 | #include <TMath.h> |
504de69b | 28 | #include <TROOT.h> |
c5f0f3c1 | 29 | //---- AliRoot headers ----- |
d681bb2d | 30 | #include "AliESDVertex.h" |
504de69b | 31 | |
32 | ||
d681bb2d | 33 | ClassImp(AliESDVertex) |
504de69b | 34 | |
c5f0f3c1 | 35 | //-------------------------------------------------------------------------- |
d681bb2d | 36 | AliESDVertex::AliESDVertex() { |
c5f0f3c1 | 37 | // |
38 | // Default Constructor, set everything to 0 | |
39 | // | |
0fc5cc25 | 40 | SetToZero(); |
c5f0f3c1 | 41 | } |
42 | //-------------------------------------------------------------------------- | |
d681bb2d | 43 | AliESDVertex::AliESDVertex(Double_t positionZ,Double_t sigmaZ, |
a27d7b02 | 44 | Int_t nContributors,const Char_t *vtxName) { |
0fc5cc25 | 45 | // |
46 | // Constructor for vertex Z from pixels | |
47 | // | |
c5f0f3c1 | 48 | |
0fc5cc25 | 49 | SetToZero(); |
c5f0f3c1 | 50 | |
0fc5cc25 | 51 | fPosition[2] = positionZ; |
52 | fCovZZ = sigmaZ*sigmaZ; | |
53 | fNContributors = nContributors; | |
54 | SetName(vtxName); | |
c5f0f3c1 | 55 | |
c5f0f3c1 | 56 | } |
57 | //------------------------------------------------------------------------- | |
d681bb2d | 58 | AliESDVertex::AliESDVertex(Double_t position[3],Double_t covmatrix[6], |
a27d7b02 | 59 | Double_t chi2,Int_t nContributors, const Char_t *vtxName) { |
c5f0f3c1 | 60 | // |
61 | // Constructor for vertex in 3D from tracks | |
62 | // | |
0fc5cc25 | 63 | |
64 | SetToZero(); | |
c5f0f3c1 | 65 | fPosition[0] = position[0]; |
66 | fPosition[1] = position[1]; | |
67 | fPosition[2] = position[2]; | |
68 | fCovXX = covmatrix[0]; | |
69 | fCovXY = covmatrix[1]; | |
70 | fCovYY = covmatrix[2]; | |
71 | fCovXZ = covmatrix[3]; | |
72 | fCovYZ = covmatrix[4]; | |
73 | fCovZZ = covmatrix[5]; | |
74 | ||
c5f0f3c1 | 75 | |
c5f0f3c1 | 76 | fChi2 = chi2; |
77 | fNContributors = nContributors; | |
78 | ||
79 | SetName(vtxName); | |
80 | ||
cd77595e | 81 | } |
c5f0f3c1 | 82 | //-------------------------------------------------------------------------- |
d681bb2d | 83 | AliESDVertex::AliESDVertex(Double_t position[3],Double_t sigma[3], |
a27d7b02 | 84 | const Char_t *vtxName) { |
c5f0f3c1 | 85 | // |
86 | // Constructor for smearing of true position | |
87 | // | |
0fc5cc25 | 88 | |
89 | SetToZero(); | |
c5f0f3c1 | 90 | fPosition[0] = position[0]; |
91 | fPosition[1] = position[1]; | |
92 | fPosition[2] = position[2]; | |
93 | fCovXX = sigma[0]*sigma[0]; | |
94 | fCovXY = 0; | |
95 | fCovYY = sigma[1]*sigma[1]; | |
96 | fCovXZ = 0; | |
97 | fCovYZ = 0; | |
98 | fCovZZ = sigma[2]*sigma[2]; | |
99 | ||
c5f0f3c1 | 100 | |
101 | SetName(vtxName); | |
102 | ||
c5f0f3c1 | 103 | } |
104 | //-------------------------------------------------------------------------- | |
d681bb2d | 105 | AliESDVertex::AliESDVertex(Double_t position[3],Double_t sigma[3], |
a27d7b02 | 106 | Double_t snr[3], const Char_t *vtxName) { |
0fc5cc25 | 107 | // |
108 | // Constructor for Pb-Pb | |
109 | // | |
110 | ||
111 | SetToZero(); | |
112 | fPosition[0] = position[0]; | |
113 | fPosition[1] = position[1]; | |
114 | fPosition[2] = position[2]; | |
115 | fCovXX = sigma[0]*sigma[0]; | |
116 | fCovXY = 0; | |
117 | fCovYY = sigma[1]*sigma[1]; | |
118 | fCovXZ = 0; | |
119 | fCovYZ = 0; | |
120 | fCovZZ = sigma[2]*sigma[2]; | |
121 | ||
122 | fSNR[0] = snr[0]; | |
123 | fSNR[1] = snr[1]; | |
124 | fSNR[2] = snr[2]; | |
125 | ||
126 | SetName(vtxName); | |
c5f0f3c1 | 127 | |
0fc5cc25 | 128 | } |
129 | //-------------------------------------------------------------------------- | |
d681bb2d | 130 | void AliESDVertex::SetToZero() { |
0fc5cc25 | 131 | // |
132 | // Set some data members to 0. Used by constructors | |
133 | // | |
134 | for(Int_t i=0; i<3; i++){ | |
135 | fPosition[i] = 0.; | |
136 | fTruePos[i] = 0; | |
137 | fSNR[i] = 0.; | |
138 | } | |
2257f27e | 139 | fCovXX = 0.005*0.005; |
0fc5cc25 | 140 | fCovXY = 0; |
2257f27e | 141 | fCovYY = 0.005*0.005; |
0fc5cc25 | 142 | fCovXZ = 0; |
143 | fCovYZ = 0; | |
2257f27e | 144 | fCovZZ = 5.3*5.3; |
0fc5cc25 | 145 | |
0fc5cc25 | 146 | fChi2 = 0; |
147 | fNContributors = 0; | |
148 | ||
149 | SetDebug(); | |
c5f0f3c1 | 150 | } |
151 | //-------------------------------------------------------------------------- | |
d681bb2d | 152 | AliESDVertex::~AliESDVertex() { |
c5f0f3c1 | 153 | // |
154 | // Default Destructor | |
155 | // | |
0fc5cc25 | 156 | |
c5f0f3c1 | 157 | } |
158 | //-------------------------------------------------------------------------- | |
d681bb2d | 159 | void AliESDVertex::GetXYZ(Double_t position[3]) const { |
c5f0f3c1 | 160 | // |
161 | // Return position of the vertex in global frame | |
c5f0f3c1 | 162 | // |
163 | position[0] = fPosition[0]; | |
164 | position[1] = fPosition[1]; | |
165 | position[2] = fPosition[2]; | |
166 | ||
167 | return; | |
168 | } | |
169 | //-------------------------------------------------------------------------- | |
d681bb2d | 170 | void AliESDVertex::GetSigmaXYZ(Double_t sigma[3]) const { |
c5f0f3c1 | 171 | // |
c5f0f3c1 | 172 | // Return errors on vertex position in thrust frame |
173 | // | |
174 | sigma[0] = TMath::Sqrt(fCovXX); | |
175 | sigma[1] = TMath::Sqrt(fCovYY); | |
176 | sigma[2] = TMath::Sqrt(fCovZZ); | |
177 | ||
178 | return; | |
cd77595e | 179 | } |
c5f0f3c1 | 180 | //-------------------------------------------------------------------------- |
d681bb2d | 181 | void AliESDVertex::GetCovMatrix(Double_t covmatrix[6]) const { |
c5f0f3c1 | 182 | // |
183 | // Return covariance matrix of the vertex | |
c5f0f3c1 | 184 | // |
185 | covmatrix[0] = fCovXX; | |
186 | covmatrix[1] = fCovXY; | |
187 | covmatrix[2] = fCovYY; | |
188 | covmatrix[3] = fCovXZ; | |
189 | covmatrix[4] = fCovYZ; | |
190 | covmatrix[5] = fCovZZ; | |
191 | ||
192 | return; | |
cd77595e | 193 | } |
c5f0f3c1 | 194 | //-------------------------------------------------------------------------- |
d681bb2d | 195 | Double_t AliESDVertex::GetXv() const { |
c5f0f3c1 | 196 | // |
197 | // Return global x | |
198 | // | |
11ba84a4 | 199 | return fPosition[0]; |
c5f0f3c1 | 200 | } |
201 | //-------------------------------------------------------------------------- | |
d681bb2d | 202 | Double_t AliESDVertex::GetYv() const { |
c5f0f3c1 | 203 | // |
204 | // Return global y | |
205 | // | |
11ba84a4 | 206 | return fPosition[1]; |
c5f0f3c1 | 207 | } |
208 | //-------------------------------------------------------------------------- | |
d681bb2d | 209 | Double_t AliESDVertex::GetZv() const { |
c5f0f3c1 | 210 | // |
211 | // Return global z | |
212 | // | |
213 | return fPosition[2]; | |
214 | } | |
215 | //-------------------------------------------------------------------------- | |
d681bb2d | 216 | Double_t AliESDVertex::GetXRes() const { |
c5f0f3c1 | 217 | // |
218 | // Return error on global x | |
219 | // | |
11ba84a4 | 220 | return TMath::Sqrt(fCovXX); |
c5f0f3c1 | 221 | } |
222 | //-------------------------------------------------------------------------- | |
d681bb2d | 223 | Double_t AliESDVertex::GetYRes() const { |
c5f0f3c1 | 224 | // |
225 | // Return error on global y | |
226 | // | |
11ba84a4 | 227 | return TMath::Sqrt(fCovYY); |
c5f0f3c1 | 228 | } |
229 | //-------------------------------------------------------------------------- | |
d681bb2d | 230 | Double_t AliESDVertex::GetZRes() const { |
c5f0f3c1 | 231 | // |
232 | // Return error on global z | |
233 | // | |
11ba84a4 | 234 | return TMath::Sqrt(fCovZZ); |
c5f0f3c1 | 235 | } |
236 | //-------------------------------------------------------------------------- | |
d681bb2d | 237 | void AliESDVertex::GetSNR(Double_t snr[3]) const { |
c5f0f3c1 | 238 | // |
239 | // Return S/N ratios | |
240 | // | |
241 | for(Int_t i=0;i<3;i++) snr[i] = fSNR[i]; | |
242 | ||
243 | return; | |
244 | } | |
245 | //-------------------------------------------------------------------------- | |
d681bb2d | 246 | void AliESDVertex::PrintStatus() const { |
c5f0f3c1 | 247 | // |
248 | // Print out information on all data members | |
249 | // | |
c5f0f3c1 | 250 | printf(" Vertex position:\n"); |
251 | printf(" x = %f +- %f\n",fPosition[0],TMath::Sqrt(fCovXX)); | |
252 | printf(" y = %f +- %f\n",fPosition[1],TMath::Sqrt(fCovYY)); | |
253 | printf(" z = %f +- %f\n",fPosition[2],TMath::Sqrt(fCovZZ)); | |
254 | printf(" Covariance matrix:\n"); | |
255 | printf(" %12.10f %12.10f %12.10f\n %12.10f %12.10f %12.10f\n %12.10f %12.10f %12.10f\n",fCovXX,fCovXY,fCovXZ,fCovXY,fCovYY,fCovYZ,fCovXZ,fCovYZ,fCovZZ); | |
256 | printf(" S/N = (%f, %f, %f)\n",fSNR[0],fSNR[1],fSNR[2]); | |
257 | printf(" chi2 = %f\n",fChi2); | |
258 | printf(" # tracks (or tracklets) = %d\n",fNContributors); | |
259 | ||
0fc5cc25 | 260 | printf(" True vertex position - for comparison: %12.10f %12.10f %12.10f\n ",fTruePos[0],fTruePos[1],fTruePos[2]); |
261 | ||
c5f0f3c1 | 262 | return; |
263 | } | |
264 | ||
265 | ||
266 | ||
6b88f180 | 267 |