102 {
103
107
108 std::string residualState( "" );
109 Decay
const &
decay = m_decayPath[0];
110 Suite<Product, Decay>
const &products =
decay.products( );
111 for( std::size_t i1 = 0; i1 < products.size( ); ++i1 ) {
112 Product
const &
product = products[i1];
113
115 }
116
117 Particle const &initialState = a_pops.get<Particle>( a_nuclideGammaBranchStateInfo.state( ) );
118 Particle const &finalState = a_pops.get<Particle>( residualState );
119 double gammaEnergy =
PoPI_AMU2MeV_c2 * ( initialState.massValue(
"amu" ) - finalState.massValue(
"amu" ) );
120
121 NuclideGammaBranchInfo nuclideGammaBranchInfo( _probability, _photonEmissionProbabilities, gammaEnergy, residualState );
122 a_nuclideGammaBranchStateInfo.add( nuclideGammaBranchInfo );
123 }
124}
#define PoPI_decayModeElectroMagnetic
PQ_suite const & photonEmissionProbabilities() const
PQ_suite const & probability() const
ParticleList decay(Cluster *const c)
Carries out a cluster decay.
double getPhysicalQuantityOfSuiteAsDouble(PQ_suite const &a_suite, bool a_allowEmpty=false, double a_emptyValue=0.0)
static std::string const photon