#ifndef ROOT_TObject
#include <TObject.h>
#endif

class pytpart : public TObject {

 private:
  
  // See pages 57-60, Pythia manual
  
  Int_t   Flavor;      // particle flavor code (KF)
  Int_t   Status;      // particle status code (KS)
  Int_t   Parent;      // parent particle KF code
  Int_t   Charge;      // charge

  // Particle's kinematics:
  Double_t E;
  Double_t M;
  Double_t Pt;
  Double_t Phi;
  Double_t Theta;
  // Parent's kinematics:
  Double_t EP;
  Double_t MP;
  Double_t PtP;
  Double_t PhiP;
  Double_t ThetaP;
  
 public:

  pytpart(){;}              // constructor
  virtual ~pytpart(){;}     // destructor

  Int_t   GetFlavor()      {return Flavor;};         
  Int_t   GetStatus()      {return Status;};     
  Int_t   GetParent()      {return Parent;};  
  Int_t   GetCharge()      {return Charge;};     
  
  Double_t GetE()           {return E;}; 
  Double_t GetM()           {return M;}; 
  Double_t GetPt()          {return Pt;}; 
  Double_t GetPhi()         {return Phi;}; 
  Double_t GetTheta()       {return Theta;};
  
  Double_t GetEP()          {return EP;}; 
  Double_t GetMP()          {return MP;}; 
  Double_t GetPtP()         {return PtP;}; 
  Double_t GetPhiP()        {return PhiP;}; 
  Double_t GetThetaP()      {return ThetaP;};

  void    SetFlavor(Int_t v)          {Flavor = v;};     
  void    SetStatus(Int_t v)          {Status = v;};     
  void    SetParent(Int_t v)          {Parent = v;};         
  void    SetCharge(Int_t v)          {Charge = v;};     
  
  void    SetE(Double_t v)             {E = v;};     
  void    SetM(Double_t v)             {M = v;};     
  void    SetPt(Double_t v)            {Pt = v;};     
  void    SetPhi(Double_t v)           {Phi = v;};     
  void    SetTheta(Double_t v)         {Theta = v;};
  
  void    SetEP(Double_t v)            {EP = v;};     
  void    SetMP(Double_t v)            {MP = v;};     
  void    SetPtP(Double_t v)           {PtP = v;};     
  void    SetPhiP(Double_t v)          {PhiP = v;};     
  void    SetThetaP(Double_t v)        {ThetaP = v;};     

  ClassDef(pytpart,2)

};