1 #ifndef ALICASCADEVERTEXER_H
2 #define ALICASCADEVERTEXER_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 //------------------------------------------------------------------
7 // Cascade Vertexer Class
8 // Reads V0s and tracks, writes out cascade vertices
9 // Origin: Christian Kuhn, IReS, Strasbourg, christian.kuhn@ires.in2p3.fr
10 //------------------------------------------------------------------
17 class AliExternalTrackParam;
19 //_____________________________________________________________________________
20 class AliCascadeVertexer : public TObject {
23 AliCascadeVertexer(const Double_t cuts[8]);
24 void SetCuts(const Double_t cuts[8]);
25 void SetVertex(Double_t *vtx) { fX=vtx[0]; fY=vtx[1]; fZ=vtx[2]; }
27 Int_t V0sTracks2CascadeVertices(AliESD *event);
28 Double_t PropagateToDCA(AliESDv0 *vtx,AliExternalTrackParam *trk,Double_t b);
30 void GetCuts(Double_t cuts[8]) const;
31 void GetVertex(Double_t *vtx) const { vtx[0]=fX; vtx[1]=fY; vtx[2]=fZ; }
34 Double_t fChi2max; // maximal allowed chi2
35 Double_t fDV0min; // min. allowed V0 impact parameter
36 Double_t fMassWin; // window around the Lambda mass
37 Double_t fDBachMin; // min. allowed bachelor impact parameter
38 Double_t fDCAmax; // maximal allowed DCA between the V0 and the track
39 Double_t fCPAmax; // maximal allowed cosine of the cascade pointing angle
40 Double_t fRmin, fRmax;// max & min radii of the fiducial volume
42 Double_t fX; // X-coordinate of the primary vertex
43 Double_t fY; // Y-coordinate of the primary vertex
44 Double_t fZ; // Z-coordinate of the primary vertex
46 ClassDef(AliCascadeVertexer,1) // cascade verterxer
49 inline AliCascadeVertexer::AliCascadeVertexer() :
65 inline AliCascadeVertexer::AliCascadeVertexer(const Double_t cuts[8]) :
81 inline void AliCascadeVertexer::SetCuts(const Double_t cuts[8]) {
83 fDV0min=cuts[1]; fMassWin=cuts[2]; fDBachMin=cuts[3];
84 fDCAmax=cuts[4]; fCPAmax=cuts[5];
85 fRmin=cuts[6]; fRmax=cuts[7];
88 inline void AliCascadeVertexer::GetCuts(Double_t cuts[8]) const {
90 cuts[1]=fDV0min; cuts[2]=fMassWin; cuts[3]=fDBachMin;
91 cuts[4]=fDCAmax; cuts[5]=fCPAmax;
92 cuts[6]=fRmin; cuts[7]=fRmax;