Keep products of decay(cascades) of primary particles on the stack.

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

diff --git a/STEER/AliAODv0.h b/STEER/AliAODv0.h
index 1b45d4dab8b63aaf18d1c365c88e5ddf62e13d0b..d3c11a94637559119b40c404cbf0a5ba566ce377 100644 (file)
--- a/STEER/AliAODv0.h
+++ b/STEER/AliAODv0.h
@@ -1,5 +1,5 @@
-#ifndef AliAODv0_H
-#define AliAODv0_H
+#ifndef ALIAODV0_H
+#define ALIAODV0_H
 /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
@@ -24,7 +24,7 @@ public:
   AliAODv0& operator=(const AliAODv0& rAliAODv0);
 
   void     Fill(AliAODVertex *rAODVertex, Double_t rDcaV0Daughters, Double_t rDcaV0ToPrimVertex,
-               const Double_t *rMomPos, const Double_t *rMomNeg, Double_t *rDcaDaughterToPrimVertex);
+               const Double_t *rMomPos, const Double_t *rMomNeg, const Double_t *rDcaDaughterToPrimVertex);
   void     ResetV0();
   void     Print(Option_t* option = "") const;
 
@@ -82,8 +82,13 @@ public:
   
   Short_t  GetPosID()       const;
   Short_t  GetNegID()       const;
-  Int_t    GetLabel() const {return -1;} // Dummy
+  Int_t    GetLabel()       const {return -1;} // Dummy
+  Int_t    PdgCode()        const {return  0;} // Dummy
+  
+  virtual Bool_t   GetPxPyPz(Double_t */*p*/) const { return kFALSE; }
+  virtual void     SetID(Short_t /*id*/) {;}
 
+  
 protected:
   Double32_t fDcaV0ToPrimVertex;    // dca of V0 to primary vertex 
   Bool_t     fOnFlyStatus;          // if kTRUE, then this V0 is recontructed
@@ -111,9 +116,9 @@ inline Double_t AliAODv0::RadiusV0() const {
 }
 
 inline Double_t AliAODv0::OpenAngleV0() const {
-  Double_t lPtot1xPtot2 = PxProng(0)*PxProng(1)+PyProng(0)*PyProng(1)+PzProng(0)*PzProng(1);
-  Double_t lPtot1Ptot2_2 = Ptot2Pos()*Ptot2Neg();
-  return ::acos(lPtot1xPtot2/::sqrt(lPtot1Ptot2_2) );
+  Double_t lScalPtot1Ptot2 = PxProng(0)*PxProng(1)+PyProng(0)*PyProng(1)+PzProng(0)*PzProng(1);
+  Double_t lPtot1xPtot2 = Ptot2Pos()*Ptot2Neg();
+  return ::acos(lScalPtot1Ptot2/::sqrt(lPtot1xPtot2) );
 }
 
 inline Double_t AliAODv0::MomPosX() const {return fPx[0];}