//  
// Author: Sophia Schade <mailto: schade@physik.hu-berlin.de>
// Update: $Id: AtlPhiDecayKK.h,v 1.3 2010/05/18 10:07:01 schade Exp $
// Copyright: 2010 (C) Sophia Schade
//
#ifndef ATLAS_AtlPhiDecayKK
#define ATLAS_AtlPhiDecayKK
#ifndef HEP_Hep2BodyDecay
#include <Hep2BodyDecay.h>
#endif
#ifndef ATLAS_AtlIDTrack
#include <AtlIDTrack.h>
#endif
#ifndef HEP_HepVertex
#include <HepVertex.h>
#endif


class AtlPhiDecayKK : public Hep2BodyDecay {
    
  private:
    
    Float_t fAngle;        // Angle between KPlus and KMinus

  public:
    AtlPhiDecayKK();
    AtlPhiDecayKK(Int_t Id, Float_t Px, Float_t Py, Float_t Pz,
		  Float_t E, AtlIDTrack *KPlus, AtlIDTrack *KMinus,
		  HepVertex *DecayVtx,HepVertex *PrimaryVtx,
		  HepParticle ReFitKPlus, HepParticle ReFitKMinus);
    virtual ~AtlPhiDecayKK();
    virtual void Clear(Option_t *option = "");
    void Print(Option_t *option = "");
    static void PrintHeader();
    static void PrintFooter();
    
    inline AtlIDTrack* GetKPlus()         { return (AtlIDTrack*)GetDaughter1();           }
    inline AtlIDTrack* GetKMinus()        { return (AtlIDTrack*)GetDaughter2();           }
    inline Float_t     GetAngleKK() const { return  fAngle;                               }
    inline HepVertex*  GetVtx()           { return (HepVertex* )GetVtx();                 }
    inline HepParticle GetReFitKPlus()    { return GetReFitDaughter1();                   }
    inline HepParticle GetReFitKMinus()   { return GetReFitDaughter2();                   }
    
    ClassDef(AtlPhiDecayKK,1) // Phi ->K-Plus + K-Minus  - decay
};
#endif