PhotosHEPEVTParticle.h

#ifndef _PhotosHEPEVTParticle_h_included_
#define _PhotosHEPEVTParticle_h_included_

#include <iostream>
#include <vector>
#include <cmath>
#include <cstdio>

#include "Photos.h"
#include "PhotosParticle.h"
#include "PhotosHEPEVTEvent.h"

namespace Photospp
{

class PhotosHEPEVTEvent;

class PhotosHEPEVTParticle: public PhotosParticle {

 public:
  ~PhotosHEPEVTParticle();

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

  void addDaughter(PhotosParticle* daughter);

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

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

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

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

  std::vector<PhotosParticle*> getAllDecayProducts();

  bool checkMomentumConservation();

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

  void createHistoryEntry();

  void createSelfDecayVertex(PhotosParticle *out);

  bool isDaughterOf(PhotosHEPEVTParticle *p);

  bool isMotherOf  (PhotosHEPEVTParticle *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(PhotosHEPEVTEvent *event);
  
  int getFirstMotherIndex();
  
  int getSecondMotherIndex();

  int getDaughterRangeStart();

  int getDaughterRangeEnd();  

private:

  void setDaughterRangeStart(int i) { m_daughter_start=i; }
  
  void setDaughterRangeEnd(int i)   { m_daughter_end  =i; }

  PhotosHEPEVTEvent *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<PhotosHEPEVTParticle*> cache;
};

} // namespace Photospp
#endif