Adding a new method calculating V0's impact parameter in 2D in XY plance
authorbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 23 Jul 2011 23:27:13 +0000 (23:27 +0000)
committerbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 23 Jul 2011 23:27:13 +0000 (23:27 +0000)
STEER/ESD/AliESDv0.cxx
STEER/ESD/AliESDv0.h

index 8cf2c30..2449759 100644 (file)
@@ -417,9 +417,22 @@ void AliESDv0::GetXYZ(Double_t &x, Double_t &y, Double_t &z) const {
   z=fPos[2]; 
 }
 
+Float_t AliESDv0::GetD(Double_t x0, Double_t y0) const {
+  //--------------------------------------------------------------------
+  // This function returns V0's impact parameter calculated in 2D in XY plane
+  //--------------------------------------------------------------------
+  Double_t x=fPos[0],y=fPos[1];
+  Double_t px=fNmom[0]+fPmom[0];
+  Double_t py=fNmom[1]+fPmom[1];
+
+  Double_t dz=(x0-x)*py - (y0-y)*px;
+  Double_t d=TMath::Sqrt(dz*dz/(px*px+py*py));
+  return d;
+}
+
 Float_t AliESDv0::GetD(Double_t x0, Double_t y0, Double_t z0) const {
   //--------------------------------------------------------------------
-  // This function returns V0's impact parameter
+  // This function returns V0's impact parameter calculated in 3D
   //--------------------------------------------------------------------
   Double_t x=fPos[0],y=fPos[1],z=fPos[2];
   Double_t px=fNmom[0]+fPmom[0];
index 09114d7..d5d5aeb 100644 (file)
@@ -77,6 +77,7 @@ public:
   void     GetNPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const;
   void     GetPPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const;
   void     GetXYZ(Double_t &x, Double_t &y, Double_t &z) const;
+  Float_t  GetD(Double_t x0,Double_t y0) const;
   Float_t  GetD(Double_t x0,Double_t y0,Double_t z0) const;
   Int_t    GetNindex() const {return fNidx;}
   Int_t    GetPindex() const {return fPidx;}