Data members and methods to work with the mean beam position (Massimo)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Mar 2007 10:58:32 +0000 (10:58 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Mar 2007 10:58:32 +0000 (10:58 +0000)
STEER/AliVertexer.cxx
STEER/AliVertexer.h

index c5db58e962d04dc88e68f645905218b991a71c77..4b41e37d984f82aa32cf1f4a0e819aae6613385c 100644 (file)
@@ -16,6 +16,7 @@
 /* $Id$ */
 
 #include "AliLog.h"
+#include "AliESDVertex.h"
 #include "AliVertexer.h"
 
 ClassImp(AliVertexer)
@@ -32,12 +33,13 @@ AliVertexer::AliVertexer() :
   fCurrentVertex(0),
   fFirstEvent(0),
   fLastEvent(0),
-  fDebug(0),
   fMult()
 {
   //
   // Default Constructor
   //
+  SetVtxStart(0.,0.,0.);
+  for(Int_t i=0;i<6;i++)fNominalCov[i]=0.;
 }
 
 
@@ -47,7 +49,6 @@ AliVertexer::AliVertexer(const AliVertexer &vtxr) :
   fCurrentVertex(vtxr.fCurrentVertex),
   fFirstEvent(vtxr.fFirstEvent),
   fLastEvent(vtxr.fLastEvent),
-  fDebug(vtxr.fDebug),
   fMult(vtxr.fMult)
 {
   // Copy constructor
@@ -75,9 +76,14 @@ AliVertexer::~AliVertexer() {
 
 }
 
-//______________________________________________________________________
-void AliVertexer::SetDebug(Int_t debug)
-{
-  AliWarning("Don't use this method any more, use AliDebug instead");
-  fDebug = debug;
+
+//---------------------------------------------------------------------------
+void  AliVertexer::SetVtxStart(AliESDVertex *vtx) 
+{ 
+//
+// Set initial vertex knowledge
+//
+  vtx->GetXYZ(fNominalPos);
+  vtx->GetCovMatrix(fNominalCov);
+  return; 
 }
index 546acfff5cdb6c8a0fe6f1f211df185bc9f4763e..2d2a05c5f21ed296317a8036ed7133c3b2e5aee3 100644 (file)
@@ -30,9 +30,14 @@ class AliVertexer : public TObject {
     virtual void FindVertices()= 0;
     virtual AliMultiplicity* GetMultiplicity() const {return fMult;}
     virtual void PrintStatus() const = 0;
-    virtual void SetDebug(Int_t debug = 0);
     virtual void SetFirstEvent(Int_t ev){fFirstEvent = ev;}
     virtual void SetLastEvent(Int_t ev){fLastEvent = ev;}
+    virtual void SetVtxStart(Double_t x,Double_t y,Double_t z) 
+           { fNominalPos[0]=x; fNominalPos[1]=y; fNominalPos[2]=z; }
+    virtual void SetVtxStartSigma(Double_t sx,Double_t sy,Double_t sz) 
+           { fNominalCov[0]=sx*sx; fNominalCov[2]=sy*sy; fNominalCov[5]=sz*sz;
+             fNominalCov[1]=0.; fNominalCov[3]=0.; fNominalCov[4]=0.; }
+    virtual void SetVtxStart(AliESDVertex *vtx);
     virtual void WriteCurrentVertex() = 0;
 
  
@@ -47,10 +52,11 @@ class AliVertexer : public TObject {
                                    //  event
     Int_t fFirstEvent;          // First event to be processed by FindVertices
     Int_t fLastEvent;           // Last event to be processed by FindVertices 
-    Int_t fDebug;               //! debug flag - verbose printing if >0
+    Double_t  fNominalPos[3];   // initial knowledge on vertex position
+    Double_t  fNominalCov[6];   // initial knowledge on vertex position
     AliMultiplicity *fMult;     //! Multiplicity object
 
-  ClassDef(AliVertexer,2);
+  ClassDef(AliVertexer,3);
 };
 
 #endif