Removing warnings (Sun)
[u/mrichter/AliRoot.git] / STEER / AliESDVertex.cxx
CommitLineData
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 33ClassImp(AliESDVertex)
504de69b 34
c5f0f3c1 35//--------------------------------------------------------------------------
d681bb2d 36AliESDVertex::AliESDVertex() {
c5f0f3c1 37//
38// Default Constructor, set everything to 0
39//
0fc5cc25 40 SetToZero();
c5f0f3c1 41}
42//--------------------------------------------------------------------------
d681bb2d 43AliESDVertex::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 58AliESDVertex::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 83AliESDVertex::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 105AliESDVertex::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 130void 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 152AliESDVertex::~AliESDVertex() {
c5f0f3c1 153//
154// Default Destructor
155//
0fc5cc25 156
c5f0f3c1 157}
158//--------------------------------------------------------------------------
d681bb2d 159void 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 170void 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 181void 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 195Double_t AliESDVertex::GetXv() const {
c5f0f3c1 196//
197// Return global x
198//
11ba84a4 199 return fPosition[0];
c5f0f3c1 200}
201//--------------------------------------------------------------------------
d681bb2d 202Double_t AliESDVertex::GetYv() const {
c5f0f3c1 203//
204// Return global y
205//
11ba84a4 206 return fPosition[1];
c5f0f3c1 207}
208//--------------------------------------------------------------------------
d681bb2d 209Double_t AliESDVertex::GetZv() const {
c5f0f3c1 210//
211// Return global z
212//
213 return fPosition[2];
214}
215//--------------------------------------------------------------------------
d681bb2d 216Double_t AliESDVertex::GetXRes() const {
c5f0f3c1 217//
218// Return error on global x
219//
11ba84a4 220 return TMath::Sqrt(fCovXX);
c5f0f3c1 221}
222//--------------------------------------------------------------------------
d681bb2d 223Double_t AliESDVertex::GetYRes() const {
c5f0f3c1 224//
225// Return error on global y
226//
11ba84a4 227 return TMath::Sqrt(fCovYY);
c5f0f3c1 228}
229//--------------------------------------------------------------------------
d681bb2d 230Double_t AliESDVertex::GetZRes() const {
c5f0f3c1 231//
232// Return error on global z
233//
11ba84a4 234 return TMath::Sqrt(fCovZZ);
c5f0f3c1 235}
236//--------------------------------------------------------------------------
d681bb2d 237void 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 246void 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