MCGIDI::Distributions::AngularTwoBody Class Reference

#include <MCGIDI_distributions.hpp>

Inheritance diagram for MCGIDI::Distributions::AngularTwoBody:

Public Member Functions
LUPI_HOST_DEVICE	AngularTwoBody ()
LUPI_HOST	AngularTwoBody (GIDI::Distributions::AngularTwoBody const &a_angularTwoBody, SetupInfo &a_setupInfo)
LUPI_HOST_DEVICE	~AngularTwoBody ()
LUPI_HOST_DEVICE double	residualMass () const
LUPI_HOST_DEVICE double	Q () const
LUPI_HOST_DEVICE Probabilities::ProbabilityBase2d_d1 *	angular () const
template<typename RNG>
LUPI_HOST_DEVICE void	sample (double a_X, Sampling::Input &a_input, RNG &&a_rng) const
template<typename RNG>
LUPI_HOST_DEVICE double	angleBiasing (Reaction const *a_reaction, double a_temperature, double a_energy_in, double a_mu_lab, RNG &&a_rng, double &a_energy_out) const
LUPI_HOST_DEVICE void	serialize (LUPI::DataBuffer &a_buffer, LUPI::DataBuffer::Mode a_mode)
LUPI_HOST_DEVICE bool	Upscatter () const
LUPI_HOST void	setModelDBRC_data2 (Sampling::Upscatter::ModelDBRC_data *a_modelDBRC_data)
template<typename RNG>
LUPI_HOST_DEVICE double	angleBiasing (LUPI_maybeUnused Reaction const *a_reaction, LUPI_maybeUnused double a_temperature, double a_energy_in, double a_mu_lab, RNG &&a_rng, double &a_energy_out) const
Public Member Functions inherited from MCGIDI::Distributions::Distribution
LUPI_HOST_DEVICE	Distribution ()
LUPI_HOST	Distribution (Type a_type, GIDI::Distributions::Distribution const &a_distribution, SetupInfo &a_setupInfo)
LUPI_HOST	Distribution (Type a_type, GIDI::Frame a_productFrame, SetupInfo &a_setupInfo)
LUPI_HOST_DEVICE MCGIDI_VIRTUAL_FUNCTION	~Distribution () MCGIDI_TRUE_VIRTUAL
LUPI_HOST_DEVICE Type	type () const
LUPI_HOST_DEVICE GIDI::Frame	productFrame () const
LUPI_HOST_DEVICE double	projectileMass () const
LUPI_HOST_DEVICE double	targetMass () const
LUPI_HOST_DEVICE double	productMass () const
LUPI_HOST void	setModelDBRC_data (Sampling::Upscatter::ModelDBRC_data *a_modelDBRC_data)
template<typename RNG>
LUPI_HOST_DEVICE MCGIDI_VIRTUAL_FUNCTION void	sample (double a_X, Sampling::Input &a_input, RNG &&a_rng) const MCGIDI_TRUE_VIRTUAL
template<typename RNG>
LUPI_HOST_DEVICE MCGIDI_VIRTUAL_FUNCTION double	angleBiasing (Reaction const *a_reaction, double a_temperature, double a_energy_in, double a_mu_lab, RNG &&a_rng, double &a_energy_out) const MCGIDI_TRUE_VIRTUAL
LUPI_HOST_DEVICE void	serialize (LUPI::DataBuffer &a_buffer, LUPI::DataBuffer::Mode a_mode)
template<typename RNG>
LUPI_HOST_DEVICE void	sample (double a_X, MCGIDI::Sampling::Input &a_input, RNG &&a_rng) const
template<typename RNG>
LUPI_HOST_DEVICE double	angleBiasing (Reaction const *a_reaction, double a_temperature, double a_energy_in, double a_mu_lab, RNG &&a_rng, double &a_energy_out) const

Detailed Description

This class represents the distribution for an outgoing product for a two-body interaction.

Definition at line 66 of file MCGIDI_distributions.hpp.

Constructor & Destructor Documentation

AngularTwoBody() [1/2]

LUPI_HOST_DEVICE MCGIDI::Distributions::AngularTwoBody::AngularTwoBody

(

)

Base contructor.

Definition at line 120 of file MCGIDI_distributions.cc.

                                                 :
        m_residualMass( 0.0 ),
        m_Q( 0.0 ),
        m_twoBodyThreshold( 0.0 ),
        m_Upscatter( false ),
        m_angular( nullptr ),
        m_modelDBRC_data( nullptr )  {
 
}

AngularTwoBody() [2/2]

LUPI_HOST MCGIDI::Distributions::AngularTwoBody::AngularTwoBody	(	GIDI::Distributions::AngularTwoBody const &	a_angularTwoBody,
		SetupInfo &	a_setupInfo )

Parameters

a_angularTwoBody	[in] The GIDI::Distributions::AngularTwoBody instance whose data is to be used to construct this.
a_setupInfo	[in] Used internally when constructing a Protare to pass information to other constructors.

Definition at line 135 of file MCGIDI_distributions.cc.

                                                                                                                          :
        Distribution( Type::angularTwoBody, a_angularTwoBody, a_setupInfo ),
        m_residualMass( a_setupInfo.m_product2Mass ),                           // Includes nuclear excitation energy.
        m_Q( a_setupInfo.m_Q ),
        m_twoBodyThreshold( a_setupInfo.m_reaction->twoBodyThreshold( ) ),
        m_Upscatter( false ),
        m_angular( Probabilities::parseProbability2d_d1( a_angularTwoBody.angular( ), &a_setupInfo ) ),
        m_modelDBRC_data( nullptr )  {
 
    if( a_setupInfo.m_protare.projectileIntid( ) == PoPI::Intids::neutron ) {
        m_Upscatter = a_setupInfo.m_reaction->ENDF_MT( ) == 2;
    }
}

~AngularTwoBody()

LUPI_HOST_DEVICE MCGIDI::Distributions::AngularTwoBody::~AngularTwoBody

(

)

Definition at line 152 of file MCGIDI_distributions.cc.

                                                  {
 
    delete m_angular;
    delete m_modelDBRC_data;
}

Member Function Documentation

angleBiasing() [1/2]

template<typename RNG>

LUPI_HOST_DEVICE double MCGIDI::Distributions::AngularTwoBody::angleBiasing	(	LUPI_maybeUnused Reaction const *	a_reaction,
		LUPI_maybeUnused double	a_temperature,
		double	a_energy_in,
		double	a_mu_lab,
		RNG &&	a_rng,
		double &	a_energy_out ) const

Returns the probability for a projectile with energy a_energy_in to cause a particle to be emitted at angle a_mu_lab as seen in the lab frame. a_energy_out is the sampled outgoing energy.

Parameters

a_reaction	[in] The reaction containing the particle which this distribution describes.
a_temperature	[in] The temperature of the material.
a_energy_in	[in] The energy of the incident particle.
a_mu_lab	[in] The desired mu in the lab frame for the emitted particle.
a_rng	[in] The random number generator function that returns a double in the range [0, 1.0).
a_energy_out	[out] The energy of the emitted outgoing particle.

Returns

The probability of emitting outgoing particle into lab angle a_mu_lab.

Definition at line 410 of file MCGIDI_headerSource.hpp.

                                                                                                {
 
    a_energy_out = 0.0;
 
    double initialMass = projectileMass( ) + targetMass( ), finalMass = productMass( ) + m_residualMass;
    double _x = targetMass( ) * ( a_energy_in - m_twoBodyThreshold ) / ( finalMass * finalMass );
    double Kp;                      // Total kinetic energy of products in the center-of-mass.
 
    if( _x < 2e-5 ) {
        Kp = finalMass * _x * ( 1 - 0.5 * _x * ( 1 - _x ) ); }
    else {          // This is the relativistic formula derived from E^2 - (pc)^2 which is frame independent (i.e., an invariant).
        Kp = sqrt( finalMass * finalMass + 2.0 * targetMass( ) * ( a_energy_in - m_twoBodyThreshold ) ) - finalMass;
    }
    if( Kp < 0 ) Kp = 0.;           // FIXME There needs to be a better test here.
 
    double energy_product_com = 0.5 * Kp * ( Kp + 2.0 * m_residualMass ) / ( Kp + productMass( ) + m_residualMass );
 
    if( productMass( ) == 0.0 ) {
        double boostBeta = sqrt( a_energy_in * ( a_energy_in + 2. * projectileMass( ) ) ) / ( a_energy_in + initialMass );  // Good, even for projectileMass = 0.
        double one_mu_beta = 1.0 - a_mu_lab * boostBeta;
        double mu_com = ( a_mu_lab - boostBeta ) / one_mu_beta;
        double Jacobian = ( 1.0 - boostBeta * boostBeta ) / ( one_mu_beta * one_mu_beta );
 
        a_energy_out = sqrt( 1.0 - boostBeta * boostBeta ) * energy_product_com * ( 1.0 + mu_com * boostBeta );
 
        return( Jacobian * m_angular->evaluate( a_energy_in, mu_com ) );
    }
 
    double productBeta = MCGIDI_particleBeta( productMass( ), energy_product_com );
    double boostBeta = sqrt( a_energy_in * ( a_energy_in + 2. * projectileMass( ) ) ) / ( a_energy_in + initialMass );      // beta = v/c.
    double muPlus = 0.0, JacobianPlus = 0.0, muMinus = 0.0, JacobianMinus = 0.0;
 
    int numberOfMus = muCOM_From_muLab( a_mu_lab, boostBeta, productBeta, muPlus, JacobianPlus, muMinus, JacobianMinus );
 
    if( numberOfMus == 0 ) return( 0.0 );
 
    double probability = JacobianPlus * m_angular->evaluate( a_energy_in, muPlus );
 
    if( numberOfMus == 2 ) {
        double probabilityMinus = JacobianMinus * m_angular->evaluate( a_energy_in, muMinus );
        probability += probabilityMinus;
        if( probabilityMinus > a_rng( ) * probability ) {
            muPlus = muMinus;
        }
    }
 
    double productBeta2 = productBeta * productBeta;
    double productBetaLab2 = productBeta2 + boostBeta * boostBeta * ( 1.0 - productBeta2 * ( 1.0 - muPlus * muPlus ) ) + 2.0 * muPlus * productBeta * boostBeta;
    productBetaLab2 /= 1.0 - muPlus * productBeta * boostBeta;
    a_energy_out = MCGIDI::particleKineticEnergyFromBeta2( productMass( ), productBetaLab2 );
 
    return( probability );
}

angleBiasing() [2/2]

template<typename RNG>

LUPI_HOST_DEVICE double MCGIDI::Distributions::AngularTwoBody::angleBiasing	(	Reaction const *	a_reaction,
		double	a_temperature,
		double	a_energy_in,
		double	a_mu_lab,
		RNG &&	a_rng,
		double &	a_energy_out ) const

angular()

LUPI_HOST_DEVICE Probabilities::ProbabilityBase2d_d1 * MCGIDI::Distributions::AngularTwoBody::angular ( ) const

inline

Returns the value of the m_angular.

Definition at line 86 of file MCGIDI_distributions.hpp.

Referenced by AngularTwoBody().

Q()

LUPI_HOST_DEVICE double MCGIDI::Distributions::AngularTwoBody::Q ( ) const

inline

Returns the value of the m_Q.

Definition at line 85 of file MCGIDI_distributions.hpp.

residualMass()

LUPI_HOST_DEVICE double MCGIDI::Distributions::AngularTwoBody::residualMass ( ) const

inline

Returns the value of the m_residualMass.

Definition at line 84 of file MCGIDI_distributions.hpp.

sample()

template<typename RNG>

LUPI_HOST_DEVICE void MCGIDI::Distributions::AngularTwoBody::sample	(	double	a_X,
		Sampling::Input &	a_input,
		RNG &&	a_rng ) const

This method samples the outgoing product data for the two outgoing particles in a two-body outgoing channel. First, is samples mu, the cosine of the product's outgoing angle, since this is for two-body interactions, mu is in the center-of-mass frame. It then calls kinetics_COMKineticEnergy2LabEnergyAndMomentum.

Parameters

a_X	[in] The energy of the projectile in the lab frame.
a_input	[in] Sample options requested by user.
a_rng	[in] The random number generator function that returns a double in the range [0, 1.0).

Definition at line 366 of file MCGIDI_headerSource.hpp.

                                                                                                                              {
 
    double initialMass = projectileMass( ) + targetMass( ), finalMass = productMass( ) + m_residualMass;
    double beta = sqrt( a_X * ( a_X + 2. * projectileMass( ) ) ) / ( a_X + initialMass );      // beta = v/c.
    double _x = targetMass( ) * ( a_X - m_twoBodyThreshold ) / ( finalMass * finalMass );
    double Kp;                          // Kp is the total kinetic energy for m3 and m4 in the COM frame.
 
    a_input.setSampledType( Sampling::SampledType::firstTwoBody );
 
    if( m_Upscatter ) {
       if( ( a_input.m_upscatterModel == Sampling::Upscatter::Model::B ) || ( a_input.m_upscatterModel == Sampling::Upscatter::Model::BSnLimits )
                || ( a_input.m_upscatterModel == Sampling::Upscatter::Model::DBRC ) ) {
            if( upscatterModelB( a_X, a_input, a_rng ) ) return;
        }
    }
 
    if( _x < 2e-5 ) {
        Kp = finalMass * _x * ( 1 - 0.5 * _x * ( 1 - _x ) ); }
    else {          // This is the relativistic formula derived from E^2 - (pc)^2 is frame independent.
        Kp = sqrt( finalMass * finalMass + 2 * targetMass( ) * ( a_X - m_twoBodyThreshold ) ) - finalMass;
    }
    if( Kp < 0 ) Kp = 0.;           // FIXME There needs to be a better test here.
 
    a_input.m_mu = m_angular->sample( a_X, a_rng( ), a_rng );
    a_input.m_phi = 2. * M_PI * a_rng( );
    kinetics_COMKineticEnergy2LabEnergyAndMomentum( beta, Kp, productMass( ), m_residualMass, a_input );
}

serialize()

LUPI_HOST_DEVICE void MCGIDI::Distributions::AngularTwoBody::serialize	(	LUPI::DataBuffer &	a_buffer,
		LUPI::DataBuffer::Mode	a_mode )

This method serializes this for broadcasting as needed for MPI and GPUs. The method can count the number of required bytes, pack this or unpack this depending on a_mode.

Parameters

a_buffer	[in] The buffer to read or write data to depending on a_mode.
a_mode	[in] Specifies the action of this method.

Definition at line 166 of file MCGIDI_distributions.cc.

                                                                                                     {
 
    Distribution::serialize( a_buffer, a_mode );
    DATA_MEMBER_DOUBLE( m_residualMass, a_buffer, a_mode );
    DATA_MEMBER_DOUBLE( m_Q, a_buffer, a_mode );
    DATA_MEMBER_DOUBLE( m_twoBodyThreshold, a_buffer, a_mode );
    DATA_MEMBER_INT( m_Upscatter, a_buffer, a_mode );
 
    m_angular = serializeProbability2d_d1( a_buffer, a_mode, m_angular );
    m_modelDBRC_data = serializeModelDBRC_data( a_buffer, a_mode, m_modelDBRC_data );
}

setModelDBRC_data2()

LUPI_HOST void MCGIDI::Distributions::AngularTwoBody::setModelDBRC_data2

(

Sampling::Upscatter::ModelDBRC_data *

a_modelDBRC_data

)

This method sets this m_modelDBRC_data to a_modelDBRC_data. It also deletes the current m_modelDBRC_data member.

Parameters

a_modelDBRC_data

[in] The instance storing data needed to treat the DRRC upscatter mode.

Definition at line 184 of file MCGIDI_distributions.cc.

                                                                                                     {
 
    delete m_modelDBRC_data;
    m_modelDBRC_data = a_modelDBRC_data;
}

Upscatter()

LUPI_HOST_DEVICE bool MCGIDI::Distributions::AngularTwoBody::Upscatter ( ) const

inline

Returns the value of the m_Upscatter.

Definition at line 93 of file MCGIDI_distributions.hpp.

---

The documentation for this class was generated from the following files:

• MCGIDI_distributions.hpp
• MCGIDI_headerSource.hpp
• MCGIDI_distributions.cc