4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 //-------------------------------------------------------------------------
10 // ESD V0 Vertex Class
11 // This class is part of the Event Summary Data set of classes
12 // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
13 // Modified by: Marian Ivanov, CERN, Marian.Ivanov@cern.ch
14 // and Boris Hippolyte,IPHC, hippolyt@in2p3.fr
15 //-------------------------------------------------------------------------
19 #include "AliESDV0Params.h"
20 #include "AliExternalTrackParam.h"
22 class AliESDv0 : public TObject {
25 AliESDv0(const AliExternalTrackParam &t1, Int_t i1,
26 const AliExternalTrackParam &t2, Int_t i2);
28 AliESDv0(const AliESDv0&);
31 Double_t ChangeMassHypothesis(Int_t code=kK0Short);
33 Int_t GetPdgCode() const {return fPdgCode;}
34 Float_t GetEffMass(UInt_t p1, UInt_t p2);
35 Float_t GetEffMass() const {return fEffMass;}
36 Float_t GetChi2V0() const {return fChi2V0;}
37 void GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const;
38 void GetNPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const;
39 void GetPPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const;
40 void GetXYZ(Double_t &x, Double_t &y, Double_t &z) const;
41 Float_t GetD(Double_t x0=0.,Double_t y0=0.,Double_t z0=0.) const;
42 Int_t GetNindex() const {return fNidx;}
43 Int_t GetPindex() const {return fPidx;}
44 void SetDcaV0Daughters(Double_t rDcaV0Daughters=0.);
45 Float_t GetDcaV0Daughters() {return fDcaV0Daughters;}
46 Float_t GetV0CosineOfPointingAngle(Double_t&, Double_t&, Double_t&) const;
47 Float_t GetV0CosineOfPointingAngle() const {return fPointAngle;}
48 void SetV0CosineOfPointingAngle(Double_t cpa) {fPointAngle=cpa;}
49 void SetOnFlyStatus(Bool_t status){fOnFlyStatus=status;}
50 Bool_t GetOnFlyStatus() const {return fOnFlyStatus;}
51 const AliExternalTrackParam *GetParamP() const {return &fParamP;}
52 const AliExternalTrackParam *GetParamN() const {return &fParamN;}
56 // **** The following member functions need to be revised ***
58 void GetPosCov(Double_t cov[6])const ; // getter for the covariance matrix of the V0 position
59 Double_t GetSigmaY(); // sigma of y coordinate at vertex posistion
60 Double_t GetSigmaZ(); // sigma of z coordinate at vertex posistion
61 Double_t GetSigmaAP0(); // calculate sigma of Point angle resolution at vertex pos.
62 Double_t GetSigmaD0(); // calculate sigma of position resolution at vertex pos.
63 Double_t GetEffectiveSigmaAP0(); // calculate sigma of point angle resolution at vertex pos. effecive parameterization
64 Double_t GetEffectiveSigmaD0(); // calculate sigma of position resolution at vertex pos.
65 Double_t GetMinimaxSigmaAP0(); // calculate mini-max sigma of point angle resolution
66 Double_t GetMinimaxSigmaD0(); // calculate mini-max sigma of dca resolution
67 Double_t GetLikelihoodAP(Int_t mode0, Int_t mode1); // get likelihood for point angle
68 Double_t GetLikelihoodD(Int_t mode0, Int_t mode1); // get likelihood for DCA
69 Double_t GetLikelihoodC(Int_t mode0, Int_t mode1); // get likelihood for Causality
72 static const AliESDV0Params & GetParameterization(){return fgkParams;}
73 void SetParamP(const AliExternalTrackParam & paramP) {fParamP = paramP;}
74 void SetParamN(const AliExternalTrackParam & paramN) {fParamN = paramN;}
75 void SetStatus(Int_t status){fStatus=status;}
76 Int_t GetStatus() const {return fStatus;}
77 Int_t GetIndex(Int_t i) const {return (i==0) ? fNidx : fPidx;}
78 void SetIndex(Int_t i, Int_t ind) {(i==0) ? (fNidx=ind) : (fPidx=ind);}
79 Double_t *GetAnglep() {return fAngle;}
80 Double_t GetRr() const {return fRr;}
81 Double_t GetDistSigma() const {return fDistSigma;}
82 void SetDistSigma(Double_t ds) {fDistSigma=ds;}
83 Float_t GetChi2Before() const {return fChi2Before;}
84 void SetChi2Before(Float_t cb) {fChi2Before=cb;}
85 Float_t GetChi2After() const {return fChi2After;}
86 void SetChi2After(Float_t ca) {fChi2After=ca;}
87 Float_t GetNAfter() const {return fNAfter;}
88 void SetNAfter(Float_t na) {fNAfter=na;}
89 Float_t GetNBefore() const {return fNBefore;}
90 void SetNBefore(Float_t nb) {fNBefore=nb;}
91 void SetCausality(Float_t pb0, Float_t pb1, Float_t pa0, Float_t pa1);
92 const Float_t * GetCausalityP() const {return fCausality;}
93 void SetClusters(Int_t *clp, Int_t *clm);
94 const Int_t * GetClusters(Int_t i) const {return fClusters[i];}
95 void SetNormDCAPrim(Float_t nd0, Float_t nd1){fNormDCAPrim[0] = nd0; fNormDCAPrim[1]=nd1;}
96 const Float_t *GetNormDCAPrimP() const {return fNormDCAPrim;}
99 Bool_t fOnFlyStatus; // if kTRUE, then this V0 is recontructed
100 // "on fly" during the tracking
102 Int_t fPdgCode; // reconstructed V0's type (PDG code)
103 Float_t fEffMass; // reconstructed V0's effective mass
104 Float_t fDcaV0Daughters; // dca between V0's daughters
105 Float_t fPointAngle; //cosine of the pointing angle
106 Float_t fChi2V0; // V0's chi2 value
108 Double32_t fPos[3]; // V0's position (global)
109 Double32_t fPosCov[6]; // covariance matrix of the vertex position
111 Int_t fNidx; // index of the negative daughter
112 Double32_t fNmom[3]; // momentum of the negative daughter (global)
113 AliExternalTrackParam fParamN; // external parameters of negative particle
114 Int_t fPidx; // index of the positive daughter
115 Double32_t fPmom[3]; // momentum of the positive daughter (global)
116 AliExternalTrackParam fParamP; // external parameters of positive particle
119 // **** The following data members need to be revised ***
121 Int_t fClusters[2][6]; //! its clusters
123 Float_t fNormDCAPrim[2]; // normalize distance to the priary vertex
125 Double32_t fAngle[3]; //three angles
126 Float_t fRr; //rec position of the vertex
127 Int_t fStatus; //status
128 Float_t fDistSigma; //sigma of distance
129 Float_t fCausality[4]; // causality information - see comments in SetCausality
130 Float_t fChi2Before; //chi2 of the tracks before V0
131 Float_t fNBefore; // number of possible points before V0
132 Float_t fChi2After; // chi2 of the tracks after V0
133 Float_t fNAfter; // number of possible points after V0
134 Float_t fPointAngleFi; //point angle fi
135 Float_t fPointAngleTh; //point angle theta
137 // parameterization coefficients
138 static AliESDV0Params fgkParams; // resolution and likelihood parameterization
141 AliESDv0& operator=(const AliESDv0&);
143 ClassDef(AliESDv0,3) // ESD V0 vertex
147 void AliESDv0::GetNPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const {
148 px=fNmom[0]; py=fNmom[1]; pz=fNmom[2];
152 void AliESDv0::GetPPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const {
153 px=fPmom[0]; py=fPmom[1]; pz=fPmom[2];
157 void AliESDv0::SetDcaV0Daughters(Double_t rDcaV0Daughters){
158 fDcaV0Daughters=rDcaV0Daughters;