TauolaHEPEVTParticle.h

#ifndef _TauolaHEPEVTParticle_h_included_
#define _TauolaHEPEVTParticle_h_included_

#include <iostream>
#include <vector>

#include "f_Decay.h"
#include "TauolaParticle.h"
#include "TauolaHEPEVTEvent.h"

namespace Tauolapp
{

class TauolaHEPEVTEvent;

class TauolaHEPEVTParticle: public TauolaParticle {

 public:
  ~TauolaHEPEVTParticle();

  TauolaHEPEVTParticle(int pdgid, int status, double px, double py, double pz, double e, double m, int ms, int me, int ds, int de);

  void undecay();

  void setMothers(std::vector<TauolaParticle*> mothers);

  void setDaughters(std::vector<TauolaParticle*> daughters);

  std::vector<TauolaParticle*> getMothers();

  std::vector<TauolaParticle*> getDaughters();

  void checkMomentumConservation();

  TauolaHEPEVTParticle * createNewParticle(int pdg_id, int status, double mass,
                                          double px, double py,
                                          double pz, double e);

  bool isDaughterOf(TauolaHEPEVTParticle *p);

  bool isMotherOf  (TauolaHEPEVTParticle *p);

  void print();

  void setPdgID(int pdg_id);

  void setStatus(int statu);

  void setMass(double mass);

  int getPdgID();

  int getStatus();

  double getMass();

  double getPx();

  double getPy();

  double getPz();

  double getE();

  void setPx( double px );

  void setPy( double py );

  void setPz( double pz );

  void setE( double e );

  int getBarcode();

  void setBarcode(int barcode);

  void setEvent(TauolaHEPEVTEvent *event);
  
  int getFirstMotherIndex();
  
  int getSecondMotherIndex();

  int getDaughterRangeStart();

  int getDaughterRangeEnd();  

private:

  TauolaHEPEVTEvent *m_event;

  int m_barcode;

  int m_first_mother, m_second_mother;

  int m_daughter_start, m_daughter_end;

  int m_pdgid;

  int m_status;

  double m_px, m_py, m_pz, m_e;

  double m_generated_mass;

  vector<TauolaHEPEVTParticle*> cache;
};

} // namespace Tauolapp
#endif