Typo fixed

[u/mrichter/AliRoot.git] / STEER / AliCascadeVertexer.h

diff --git a/STEER/AliCascadeVertexer.h b/STEER/AliCascadeVertexer.h
index c44d9a76e3a3082b6263fdd85d8b5db3499f2eef..8df34874cb50d1febba198ab42e0d776ad15cd82 100644 (file)
--- 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,60 +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() {
- fChi2max=33.; 
- fDV0min=0.015; fMassWin=0.05; fDBachMin=0.015;
- fDCAmax=0.01;  fCPAmax=0.025; 
- fRmin=0.5;     fRmax=2.5; 
+inline AliCascadeVertexer::AliCascadeVertexer() :
+  TObject(),
+  fChi2max(fgChi2max), 
+  fDV0min(fgDV0min),
+  fMassWin(fgMassWin),
+  fDBachMin(fgDBachMin),
+  fDCAmax(fgDCAmax),
+  fCPAmin(fgCPAmin), 
+  fRmin(fgRmin),
+  fRmax(fgRmax)
+{
 }
 
-inline AliCascadeVertexer::AliCascadeVertexer(const Double_t cuts[8]) {
+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]; 
-  fX=fY=fZ=0.;
 }
 
-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];
-  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