PHHepMCParticleSelectorDecayProductChain.h

Go to the documentation of this file.

#ifndef PHHepMCParticleSelectorDecayProductChain_H__
#define PHHepMCParticleSelectorDecayProductChain_H__

#include <fun4all/SubsysReco.h>

#include <HepMC/GenEvent.h>
#include <HepMC/GenParticle.h>

#include <vector>

class PHHepMCParticleSelectorDecayProductChain: public SubsysReco
{
 public:
  PHHepMCParticleSelectorDecayProductChain(const std::string &name="PARTICLESELECTOR");
  virtual ~PHHepMCParticleSelectorDecayProductChain(){}

  int InitRun(PHCompositeNode *topNode);
  int process_event(PHCompositeNode *topNode);

  virtual void SetParticle(const int pid);

  virtual void AddAncestor(const int pid);

  virtual void AddDaughter(const int pid);

  int get_embedding_id() const { return _embedding_id; }
  //
  void set_embedding_id(int id) { _embedding_id = id; }
 protected:

HepMC::GenParticle*  GetParent(HepMC::GenParticle* p, HepMC::GenEvent* event);

  int _theParticle;
  std::vector<int> _theDaughters;
  std::vector<int> _theAncestors;

  int _embedding_id;
};

#endif