HLT/global/AliFlatESDVertex.h

   1 #ifndef ALIFLATESDVERTEX_H
   2 #define ALIFLATESDVERTEX_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               *
   6  * Primary Authors : Sergey Gorbunov, Jochen Thaeder, Chiara Zampolli     */
   7
   8 /**
   9  * >> Flat structure representing a ESD vertex <<
  10  */
  11
  12 #include "Rtypes.h"
  13 #include "AliVVMisc.h"
  14 #include "AliESDVertex.h"
  15
  16 struct AliFlatESDVertex
  17 {
  18   // -- Constructor / Destructors
  19
  20   AliFlatESDVertex();
  21  ~AliFlatESDVertex() {}
  22
  23   // constructor and method for reinitialisation of virtual table
  24   AliFlatESDVertex( AliVVConstructorReinitialisationFlag );
  25   void Reinitialize() const {} // no virtual table - do nothing
  26
  27   //--
  28
  29   Double32_t fPosition[3];    // vertex position
  30   Double32_t fCov[6];  // vertex covariance matrix
  31   Int_t    fNContributors;  // # of tracklets/tracks used for the estimate
  32   Double32_t fChi2;  // chi2 of vertex fit
  33   /*
  34     Double32_t fSigma;          // track dispersion around found vertex
  35     Char_t fID;       // ID of this vertex within an ESD event
  36     Char_t fBCID;     // BC ID assigned to vertex
  37   */
  38
  39
  40   void Set(const AliESDVertex &v );
  41
  42   Double32_t GetX() const { return fPosition[0]; }
  43   Double32_t GetY() const { return fPosition[1]; }
  44   Double32_t GetZ() const { return fPosition[2]; }
  45
  46   void GetXYZ(Double_t pos[3]) const { for(Int_t j=0; j<3; j++) pos[j]=fPosition[j]; }
  47   void SetXYZ(Double_t pos[3]) { for(Int_t j=0; j<3; j++) fPosition[j]=pos[j]; }
  48
  49   void   SetX(Double_t xVert) {fPosition[0]=xVert; }
  50   void   SetY(Double_t yVert) {fPosition[1]=yVert; }
  51   void   SetZ(Double_t zVert) {fPosition[2]=zVert; }
  52   void   SetNContributors(Int_t nContr) {fNContributors=nContr; }
  53
  54
  55   Int_t    GetNContributors() const { return fNContributors; }
  56
  57   /*
  58   void     GetCovarianceMatrix(Double_t covmatrix[6]) const;
  59   void     SetCovarianceMatrix(const Double_t *) {}
  60
  61   Double_t GetChi2perNDF() const {return -999.;}
  62   Double_t GetChi2() const {return -999.;}
  63   void     SetChi2(Double_t ) {}
  64   Int_t    GetNDF() const {return -999;}
  65
  66   void     GetSigmaXYZ(Double_t sigma[3]) const;
  67   void     GetCovMatrix(Double_t covmatrix[6]) const;
  68   void     GetCovarianceMatrix(Double_t covmatrix[6]) const
  69                     {GetCovMatrix(covmatrix);}
  70   void     GetSNR(Double_t snr[3]) const;
  71   void     SetCovarianceMatrix(const Double_t *cov);
  72
  73   Double_t GetXRes() const {return TMath::Sqrt(fCovXX);}
  74   Double_t GetYRes() const {return TMath::Sqrt(fCovYY);}
  75   Double_t GetZRes() const {return TMath::Sqrt(fCovZZ);}
  76   Double_t GetXSNR() const { return fSNR[0]; }
  77   Double_t GetYSNR() const { return fSNR[1]; }
  78   Double_t GetZSNR() const { return fSNR[2]; }
  79   void     SetSNR(double snr, int i) {if (i<3 && i>=0) fSNR[i] = snr;}
  80
  81   Double_t GetChi2() const { return fChi2; }
  82   void     SetChi2(Double_t chi) { fChi2 = chi; }
  83   Double_t GetChi2toNDF() const
  84     { return fChi2/(2.*(Double_t)fNContributors-3.); }
  85   Double_t GetChi2perNDF() const { return GetChi2toNDF();}
  86   Int_t    GetNDF() const {return (2*fNContributors-3);}
  87
  88   void     Print(Option_t* option = "") const;
  89   void     PrintStatus() const {Print();}
  90
  91   void     Reset() { SetToZero(); SetName("Vertex"); }
  92
  93   void     SetID(Char_t id) {fID=id;}
  94   Char_t   GetID() const {return fID;}
  95   //
  96   Double_t GetWDist(const AliESDVertex* v) const;
  97   */
  98
  99   static size_t GetSize() { return sizeof(AliFlatESDVertex); }
 100
 101 };
 102
 103 inline AliFlatESDVertex::AliFlatESDVertex() :
 104 fNContributors(0),
 105 fChi2(0)
 106 {
 107   for( int i=0; i<3; i++) fPosition[i] = -9999;
 108   for( int i=0; i<6; i++) fCov[i] = -9999;
 109 }
 110
 111 inline AliFlatESDVertex::AliFlatESDVertex( AliVVConstructorReinitialisationFlag ):
 112 fNContributors(fNContributors),
 113 fChi2(fChi2)
 114 {
 115   // do nothing
 116   for( int i=0; i<3; i++) fPosition[i] = fPosition[i];
 117   for( int i=0; i<6; i++) fCov[i] = fCov[i];
 118 }
 119
 120 inline void AliFlatESDVertex::Set(const AliESDVertex &v )
 121 {
 122   fPosition[0] = v.GetX();
 123   fPosition[1] = v.GetY();
 124   fPosition[2] = v.GetZ();
 125   Double_t c[6];
 126   v.GetCovarianceMatrix( c );
 127   for( int i=0; i<6; i++) fCov[i] = c[i];
 128   fNContributors = v.GetNContributors();
 129   fChi2 = v.GetChi2();
 130 }
 131
 132 typedef struct AliFlatESDVertex AliFlatESDVertex;
 133
 134 #endif