GtDecayerator.h

Go to the documentation of this file.

#ifndef GTDECAYERATOR_H
#define GTDECAYERATOR_H


#include "GenTools/IHepMCEventMutator.h"
#include "GaudiAlg/GaudiTool.h"
#include "GenDecay/NucState.h"
#include "GenDecay/NucDecay.h"

#include "HepMC/GenParticle.h"

#include <string>
#include <map>


namespace GenDecay {
    class DecayRates;
}

class GtDecayerator : public GaudiTool,
                      virtual public IHepMCEventMutator
{

public:

    GtDecayerator(const std::string& type,
                  const std::string& name,
                  const IInterface* parent);
    virtual ~GtDecayerator();

    virtual StatusCode initialize();
    virtual StatusCode finalize();

    // HepMCEventMutator interface
    virtual StatusCode mutate(HepMC::GenEvent& event);

private:
    std::string m_parentNuclide;
    

    double m_parentAbundance;


    std::map<std::string,double> m_abundanceMap;


    bool m_secularEquilibrium;


    double m_correlationTime;


    GenDecay::DecayRates* m_rates;

    HepMC::GenParticle* make_particle(GenDecay::NucState& state, int status=2);

    HepMC::GenParticle* make_particle(GenDecay::NucDecay& dk);

};


#endif  // GTDECAYERATOR_H