X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FAliCascadeVertexer.h;h=8df34874cb50d1febba198ab42e0d776ad15cd82;hb=d345a2d8f455e1290ae9d68249b0a5fe0439b791;hp=a3b7a32f488627fca00b9ed48f0b765e694a52a2;hpb=fe12e09c14250ca4a380175a67ffbd6bb844662e;p=u%2Fmrichter%2FAliRoot.git

diff --git a/STEER/AliCascadeVertexer.h b/STEER/AliCascadeVertexer.h
index a3b7a32f488..8df34874cb5 100644
--- a/STEER/AliCascadeVertexer.h
+++ b/STEER/AliCascadeVertexer.h
@@ -11,8 +11,7 @@
 
 #include "TObject.h"
 
-class AliESD;
-class TTree;
+class AliESDEvent;
 class AliESDv0;
 class AliExternalTrackParam;
 
@@ -20,77 +19,87 @@ class AliExternalTrackParam;
 class AliCascadeVertexer : public TObject {
 public:
   AliCascadeVertexer();
-  AliCascadeVertexer(const Double_t cuts[8]);
   void SetCuts(const Double_t cuts[8]);
-  void SetVertex(Double_t *vtx) { fX=vtx[0]; fY=vtx[1]; fZ=vtx[2]; }
+  static void SetDefaultCuts(const Double_t cuts[8]);
+
+  Int_t V0sTracks2CascadeVertices(AliESDEvent *event);
+  Double_t Det(Double_t a00, Double_t a01, Double_t a10, Double_t a11) const;
+  Double_t Det(Double_t a00,Double_t a01,Double_t a02,
+	       Double_t a10,Double_t a11,Double_t a12,
+	       Double_t a20,Double_t a21,Double_t a22) const;
 
-  Int_t V0sTracks2CascadeVertices(AliESD *event);
   Double_t PropagateToDCA(AliESDv0 *vtx,AliExternalTrackParam *trk,Double_t b);
 
   void GetCuts(Double_t cuts[8]) const;
-  void GetVertex(Double_t *vtx) const { vtx[0]=fX; vtx[1]=fY; vtx[2]=fZ; }
+  static void GetDefaultCuts(Double_t cuts[8]);
 
 private:
+  static
+  Double_t fgChi2max;   // maximal allowed chi2 
+  static
+  Double_t fgDV0min;    // min. allowed V0 impact parameter
+  static
+  Double_t fgMassWin;   // window around the Lambda mass
+  static
+  Double_t fgDBachMin;  // min. allowed bachelor impact parameter
+  static
+  Double_t fgDCAmax;    // maximal allowed DCA between the V0 and the track 
+  static
+  Double_t fgCPAmin;    // minimal allowed cosine of the cascade pointing angle
+  static
+  Double_t fgRmin, fgRmax;// max & min radii of the fiducial volume
+  
   Double_t fChi2max;    // maximal allowed chi2 
   Double_t fDV0min;     // min. allowed V0 impact parameter
   Double_t fMassWin;    // window around the Lambda mass
   Double_t fDBachMin;   // min. allowed bachelor impact parameter
   Double_t fDCAmax;     // maximal allowed DCA between the V0 and the track 
-  Double_t fCPAmax;     // maximal allowed cosine of the cascade pointing angle
+  Double_t fCPAmin;     // minimal allowed cosine of the cascade pointing angle
   Double_t fRmin, fRmax;// max & min radii of the fiducial volume
   
-  Double_t fX;            // X-coordinate of the primary vertex
-  Double_t fY;            // Y-coordinate of the primary vertex
-  Double_t fZ;            // Z-coordinate of the primary vertex
-
-  ClassDef(AliCascadeVertexer,1)  // cascade verterxer 
+  ClassDef(AliCascadeVertexer,3)  // cascade verterxer 
 };
 
 inline AliCascadeVertexer::AliCascadeVertexer() :
   TObject(),
-  fChi2max(33.), 
-  fDV0min(0.015),
-  fMassWin(0.05),
-  fDBachMin(0.015),
-  fDCAmax(0.01),
-  fCPAmax(0.025), 
-  fRmin(0.5),
-  fRmax(2.5), 
-  fX(0),
-  fY(0),
-  fZ(0)
-{
-}
-
-inline AliCascadeVertexer::AliCascadeVertexer(const Double_t cuts[8]) :
-  TObject(),
-  fChi2max(cuts[0]), 
-  fDV0min(cuts[1]),
-  fMassWin(cuts[2]),
-  fDBachMin(cuts[3]),
-  fDCAmax(cuts[4]),
-  fCPAmax(cuts[5]),
-  fRmin(cuts[6]),
-  fRmax(cuts[7]), 
-  fX(0),
-  fY(0),
-  fZ(0)
+  fChi2max(fgChi2max), 
+  fDV0min(fgDV0min),
+  fMassWin(fgMassWin),
+  fDBachMin(fgDBachMin),
+  fDCAmax(fgDCAmax),
+  fCPAmin(fgCPAmin), 
+  fRmin(fgRmin),
+  fRmax(fgRmax)
 {
 }
 
 inline void AliCascadeVertexer::SetCuts(const Double_t cuts[8]) {
   fChi2max=cuts[0]; 
   fDV0min=cuts[1];   fMassWin=cuts[2]; fDBachMin=cuts[3];
-  fDCAmax=cuts[4];   fCPAmax=cuts[5];
+  fDCAmax=cuts[4];   fCPAmin=cuts[5];
   fRmin=cuts[6];     fRmax=cuts[7]; 
 }
 
+inline void AliCascadeVertexer::SetDefaultCuts(const Double_t cuts[8]) {
+  fgChi2max=cuts[0]; 
+  fgDV0min=cuts[1];   fgMassWin=cuts[2]; fgDBachMin=cuts[3];
+  fgDCAmax=cuts[4];   fgCPAmin=cuts[5];
+  fgRmin=cuts[6];     fgRmax=cuts[7]; 
+}
+
 inline void AliCascadeVertexer::GetCuts(Double_t cuts[8]) const {
   cuts[0]=fChi2max; 
   cuts[1]=fDV0min;   cuts[2]=fMassWin;  cuts[3]=fDBachMin;
-  cuts[4]=fDCAmax;   cuts[5]=fCPAmax;
+  cuts[4]=fDCAmax;   cuts[5]=fCPAmin;
   cuts[6]=fRmin;     cuts[7]=fRmax; 
 }
 
+inline void AliCascadeVertexer::GetDefaultCuts(Double_t cuts[8]) {
+  cuts[0]=fgChi2max; 
+  cuts[1]=fgDV0min;   cuts[2]=fgMassWin;  cuts[3]=fgDBachMin;
+  cuts[4]=fgDCAmax;   cuts[5]=fgCPAmin;
+  cuts[6]=fgRmin;     cuts[7]=fgRmax; 
+}
+
 #endif