In This Package:

DecayRates.h

Go to the documentation of this file.

#ifndef DECAYRATES_H
#define DECAYRATES_H

#include "GenDecay/NucVisitor.h"

#include "GaudiKernel/RndmGenerators.h"

#include <more/phys/ens.h>
#include <map>
#include <vector>

class IRndmGenSvc;

namespace GenDecay {

class DecayRates : public GenDecay::NucVisitor
{
public:
    typedef std::map<more::phys::nucleus,double> AbundanceMap_t;
    typedef std::pair<NucDecay*,double> DecayValuePair_t;
    typedef std::vector<DecayValuePair_t> DecayValues_t;
    typedef std::map<NucState*,double> StateValueMap_t;

    DecayRates(const AbundanceMap_t& abundance,
               bool secularEquilibrium, double corrTime, double epsTime,
               IRndmGenSvc *rgs);
    virtual ~DecayRates();

    // Return a randomly chosen decay from the chain
    NucDecay* decay();

    // Return a randomly chosen decay from a particular mother
    NucDecay* decay_state(NucState* nuc);

    // Build up map from decay to mean rate
    bool visit(NucState* mother, NucState* /*daughter*/, NucDecay* decay);

    // Access the collection of uncorrelated mothers from this chain
    const std::vector<NucState*>& mothers() const { return m_mothers; }
    double totalRate() const { return m_totalDecayRate; }

    DecayValues_t decayTimes(NucDecay* decay);

    double uni() { return m_uni(); }

private:
    void get_decay_times(NucDecay* decay, DecayValues_t& decayTime, double now);

    AbundanceMap_t m_abundance;
    bool m_secEq;
    double m_correlationTime;   // seconds
    double m_epsilonTime;        // seconds

    // secular equilibrium rate, take from chain's mother
    double m_secEqRate;

    // Cache mean decay rate of each mother
    StateValueMap_t m_meanDecayTimes;

    // uncorrelated mothers
    std::vector<NucState*> m_mothers;
    std::vector<double> m_motherBranching;
    double m_totalDecayRate;

    Rndm::Numbers m_uni;

};

}
#endif  // DECAYRATES_H

---