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