STEER/AliESDcascade.h

   1 #ifndef ALIESDCASCADE_H
   2 #define ALIESDCASCADE_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 /* $Id$ */
   8
   9 //-------------------------------------------------------------------------
  10 //                        ESD Cascade Vertex Class
  11 //               Implementation of the cascade vertex class
  12 //    Origin: Christian Kuhn, IReS, Strasbourg, christian.kuhn@ires.in2p3.fr
  13 //-------------------------------------------------------------------------
  14
  15 #include <TObject.h>
  16 #include <TPDGCode.h>
  17 #include "AliESDv0.h"
  18
  19 class AliExternalTrackParam;
  20
  21
  22 class AliESDcascade : public AliESDv0 {
  23
  24   enum  { kOmegaPlusBar =  -3334,
  25           kXiPlusBar    =  -3312,
  26           kXiMinus      =   3312,
  27           kOmegaMinus   =   3334};
  28
  29 public:
  30   AliESDcascade();
  31   AliESDcascade(const AliESDcascade&);
  32   AliESDcascade(const AliESDv0 &v0,
  33                 const AliExternalTrackParam &t, Int_t i);
  34   ~AliESDcascade();
  35   AliESDcascade& operator=(const AliESDcascade&);
  36   virtual void Copy(TObject &obj) const;
  37
  38 // Start with AliVParticle functions
  39   virtual Double_t Px() const { return fNmom[0]+fPmom[0]+fBachMom[0]; }
  40   virtual Double_t Py() const { return fNmom[1]+fPmom[1]+fBachMom[1]; }
  41   virtual Double_t Pz() const { return fNmom[2]+fPmom[2]+fBachMom[2]; }
  42   virtual Double_t Pt() const { return TMath::Sqrt(Px()*Px()+Py()*Py()); }
  43   virtual Double_t P()  const {
  44      return TMath::Sqrt(Px()*Px()+Py()*Py()+Pz()*Pz());
  45   }
  46   virtual Bool_t   PxPyPz(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; }
  47   virtual Double_t Xv() const { return fPosXi[0]; }
  48   virtual Double_t Yv() const { return fPosXi[1]; }
  49   virtual Double_t Zv() const { return fPosXi[2]; }
  50   virtual Bool_t   XvYvZv(Double_t x[3]) const { x[0] = Xv(); x[1] = Yv(); x[2] = Zv(); return kTRUE; }
  51   virtual Double_t OneOverPt() const { return (Pt() != 0.) ? 1./Pt() : -999.; }
  52   virtual Double_t Phi() const {return TMath::Pi()+TMath::ATan2(-Py(),-Px()); }
  53   virtual Double_t Theta() const {return 0.5*TMath::Pi()-TMath::ATan(Pz()/(Pt()+1.e-13)); }
  54   virtual Double_t E() const; // default is Xis but can be changed via ChangeMassHypothesis (defined in the .cxx)
  55   virtual Double_t M() const { return GetEffMassXi(); }
  56   virtual Double_t Eta() const { return 0.5*TMath::Log((P()+Pz())/(P()-Pz()+1.e-13)); }
  57   virtual Double_t Y() const { return 0.5*TMath::Log((E()+Pz())/(E()-Pz()+1.e-13)); }
  58   virtual Short_t  Charge() const { return 0; }
  59   virtual Int_t    GetLabel() const { return -1; }  // temporary
  60   virtual const Double_t *PID() const { return 0; } // return PID object ? (to be discussed!)
  61
  62   // Then the older functions
  63   Double_t ChangeMassHypothesis(Double_t &v0q, Int_t code=kXiMinus);
  64
  65   Int_t    GetPdgCodeXi() const {return fPdgCodeXi;}
  66   Double_t GetEffMassXi() const {return fEffMassXi;}
  67   Double_t GetChi2Xi()  const {return fChi2Xi;}
  68   void     GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const;
  69   void     GetXYZcascade(Double_t &x, Double_t &y, Double_t &z) const;
  70   Double_t GetDcascade(Double_t x0=0.,Double_t y0=0.,Double_t z0=0.) const;
  71
  72   void     GetBPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const {
  73      px=fBachMom[0]; py=fBachMom[1]; pz=fBachMom[2];
  74   }
  75
  76   Int_t    GetBindex() const {return fBachIdx;}
  77   void     SetIndex(Int_t i) {fBachIdx=i;}        //for the consistency with V0
  78   Int_t    GetIndex() const {return GetBindex();} //for the consistency with V0
  79   void     SetDcaXiDaughters(Double_t rDcaXiDaughters=0.);
  80   Double_t GetDcaXiDaughters() const {return fDcaXiDaughters;}
  81   Double_t GetCascadeCosineOfPointingAngle(Double_t&, Double_t&, Double_t&) const;
  82
  83   void GetPosCovXi(Double_t cov[6]) const;
  84
  85 protected:
  86
  87   Double32_t fEffMassXi;      // reconstructed cascade effective mass
  88   Double32_t fChi2Xi;         // chi2 value
  89   Double32_t fDcaXiDaughters; // dca between Xi's daughters
  90   Double32_t fPosXi[3];       // cascade vertex position (global)
  91   Double32_t fPosCovXi[6];    // covariance matrix of the vertex position
  92   Double32_t fBachMom[3];     // bachelor momentum (global)
  93   Double32_t fBachMomCov[6];  // covariance matrix of the bachelor momentum.
  94   Int_t      fPdgCodeXi;      // reconstructed cascade type (PDG code)
  95   Int_t      fBachIdx;        // label of the bachelor track
  96
  97
  98 private:
  99
 100
 101   ClassDef(AliESDcascade,5) // reconstructed cascade vertex
 102 };
 103
 104 inline
 105 void AliESDcascade::SetDcaXiDaughters(Double_t rDcaXiDaughters){
 106   fDcaXiDaughters=rDcaXiDaughters;
 107 }
 108
 109 #endif