#include <MCGIDI_sampling.hpp>

Public Member Functions
LUPI_HOST_DEVICE	ModelDBRC_data ()
LUPI_HOST	ModelDBRC_data (double a_neutronMass, double a_targetMass, Vector< double > const &a_energies, Vector< double > const &a_crossSections, DomainHash const &a_domainHash)
LUPI_HOST_DEVICE	~ModelDBRC_data ()
LUPI_HOST_DEVICE double	evaluate (double a_energy)
LUPI_HOST_DEVICE double	targetThermalSpeed (double a_temperature)
LUPI_HOST_DEVICE double	crossSectionMax (double a_energy, double a_targetThermalSpeed)
LUPI_HOST_DEVICE void	serialize (LUPI::DataBuffer &a_buffer, LUPI::DataBuffer::Mode a_mode)

Public Attributes
double	m_neutronMass
double	m_targetMass
Vector< double >	m_energies
Vector< double >	m_crossSections
Vector< std::size_t >	m_hashIndices
MCGIDI::DomainHash	m_domainHash

Detailed Description

This class is used to store the cross section for the elastic scattering upscatter model B with Doppler Broadening Rejection Correction (DBRC) with enum MCGIDI::Sampling::Upscatter::DBRC.

Definition at line 72 of file MCGIDI_sampling.hpp.

Constructor & Destructor Documentation

◆ ModelDBRC_data() [1/2]

LUPI_HOST_DEVICE MCGIDI::Sampling::Upscatter::ModelDBRC_data::ModelDBRC_data ( )

Definition at line 91 of file MCGIDI_sampling.cc.

                                                 :
        m_neutronMass( 0.0 ),
        m_targetMass( ),
        m_energies( ),
        m_crossSections( ),
        m_hashIndices( ) {
 
}

◆ ModelDBRC_data() [2/2]

LUPI_HOST MCGIDI::Sampling::Upscatter::ModelDBRC_data::ModelDBRC_data	(	double	a_neutronMass,
		double	a_targetMass,
		Vector< double > const &	a_energies,
		Vector< double > const &	a_crossSections,
		DomainHash const &	a_domainHash )

Definition at line 104 of file MCGIDI_sampling.cc.

                                                 :
        m_neutronMass( a_neutronMass ),
        m_targetMass( a_targetMass ),
        m_energies( a_energies ),
        m_crossSections( a_crossSections ),
        m_hashIndices( a_domainHash.map( a_energies ) ),
        m_domainHash( 4000, 1e-8, 10 ) {
 
}

◆ ~ModelDBRC_data()

LUPI_HOST_DEVICE MCGIDI::Sampling::Upscatter::ModelDBRC_data::~ModelDBRC_data ( )

Definition at line 118 of file MCGIDI_sampling.cc.

                                                  {
 
}

Member Function Documentation

◆ crossSectionMax()

LUPI_HOST_DEVICE double MCGIDI::Sampling::Upscatter::ModelDBRC_data::crossSectionMax	(	double	a_energy,
		double	a_targetThermalSpeed )

This method sets a_crossSectionMin and a_crossSectionMax to the minimum and maximum cross section values for the cross section in the window (a_speed - 4 * a_targetThermalSpeed) < speed < (a_speed + 4 * a_targetThermalSpeed). If (a_speed - 4 * a_targetThermalSpeed) is less than the lowest speed in the data, then it is replaced with the lowest speed.

Parameters

a_energy	[in] The energy of the projectile (i.e., incident neutron).
a_targetThermalSpeed	[in] The thermal speed of the target.

Returns: The maximum cross section in the search window.

Definition at line 164 of file MCGIDI_sampling.cc.

                                                                                                      {
 
    double crossSectionMax2 = 0.0;
    double energyFraction;
    double a_speed = MCGIDI_particleBeta( m_neutronMass, a_energy );
 
    double speedMin = a_speed - 4 * a_targetThermalSpeed;
    if( speedMin < 0.0 ) speedMin = 0.0;
    double energyMin = 0.5 * m_neutronMass * speedMin * speedMin;
    std::size_t hashIndex = m_domainHash.index( energyMin );
    std::size_t indexMin = evaluationForHashIndex( hashIndex, m_hashIndices, energyMin, m_energies, &energyFraction );
 
    double speedMax = a_speed + 4 * a_targetThermalSpeed;
    double energyMax = 0.5 * m_neutronMass * speedMax * speedMax;
    hashIndex = m_domainHash.index( energyMax );
    std::size_t indexMax = evaluationForHashIndex( hashIndex, m_hashIndices, energyMax, m_energies, &energyFraction );
    if( indexMax < m_energies.size( ) ) ++indexMax;
 
    for( std::size_t index = indexMin; index < indexMax; ++index ) {
        if( crossSectionMax2 < m_crossSections[index] ) crossSectionMax2 = m_crossSections[index];
    }
 
    return( crossSectionMax2 );
}

◆ evaluate()

LUPI_HOST_DEVICE double MCGIDI::Sampling::Upscatter::ModelDBRC_data::evaluate ( double a_energy )

This method returns the cross section evaluate at the projectile speed a_speed.

Parameters

a_temperature [in] The temperature of the target.

Returns: The thermal speed of the target.

Definition at line 130 of file MCGIDI_sampling.cc.

                                                                  {
 
    double energyFraction;
    std::size_t hashIndex = m_domainHash.index( a_energy );
    std::size_t index = evaluationForHashIndex( hashIndex, m_hashIndices, a_energy, m_energies, &energyFraction );
 
    return( energyFraction * m_crossSections[index] + ( 1.0 - energyFraction ) * m_crossSections[index+1] );
}

◆ serialize()

LUPI_HOST_DEVICE void MCGIDI::Sampling::Upscatter::ModelDBRC_data::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 197 of file MCGIDI_sampling.cc.

                                                                                                     {
 
    DATA_MEMBER_DOUBLE( m_neutronMass, a_buffer, a_mode );
    DATA_MEMBER_DOUBLE( m_targetMass, a_buffer, a_mode );
    DATA_MEMBER_VECTOR_DOUBLE( m_energies, a_buffer, a_mode );
    DATA_MEMBER_VECTOR_DOUBLE( m_crossSections, a_buffer, a_mode );
    DATA_MEMBER_VECTOR_SIZE_T( m_hashIndices, a_buffer, a_mode );
    m_domainHash.serialize( a_buffer, a_mode );
}

Referenced by MCGIDI::Sampling::Upscatter::serializeModelDBRC_data().

◆ targetThermalSpeed()

LUPI_HOST_DEVICE double MCGIDI::Sampling::Upscatter::ModelDBRC_data::targetThermalSpeed ( double a_temperature )

This method returns the thermal speed of the target with temperature a_temperature.

Parameters

a_temperature [in] The temperature of the target.

Returns: The thermal speed of the target.

Definition at line 147 of file MCGIDI_sampling.cc.

                                                                                 {
 
    return( sqrt( 2.0 * a_temperature / m_targetMass ) );
}

Member Data Documentation

◆ m_crossSections

Vector<double> MCGIDI::Sampling::Upscatter::ModelDBRC_data::m_crossSections

The cross sections corresponding to the energy grid.

Definition at line 78 of file MCGIDI_sampling.hpp.

Referenced by crossSectionMax(), evaluate(), ModelDBRC_data(), ModelDBRC_data(), and serialize().

◆ m_domainHash

MCGIDI::DomainHash MCGIDI::Sampling::Upscatter::ModelDBRC_data::m_domainHash

The hash "function".

Definition at line 80 of file MCGIDI_sampling.hpp.

Referenced by crossSectionMax(), evaluate(), ModelDBRC_data(), and serialize().

◆ m_energies

Vector<double> MCGIDI::Sampling::Upscatter::ModelDBRC_data::m_energies

The energy grid for the cross section.

Definition at line 77 of file MCGIDI_sampling.hpp.

Referenced by crossSectionMax(), evaluate(), ModelDBRC_data(), ModelDBRC_data(), and serialize().

◆ m_hashIndices

Vector<std::size_t> MCGIDI::Sampling::Upscatter::ModelDBRC_data::m_hashIndices

The indicies for the energy hash function.

Definition at line 79 of file MCGIDI_sampling.hpp.

Referenced by crossSectionMax(), evaluate(), ModelDBRC_data(), ModelDBRC_data(), and serialize().

◆ m_neutronMass

double MCGIDI::Sampling::Upscatter::ModelDBRC_data::m_neutronMass

The mass of the neutron.

Definition at line 75 of file MCGIDI_sampling.hpp.

Referenced by crossSectionMax(), ModelDBRC_data(), ModelDBRC_data(), and serialize().

◆ m_targetMass

double MCGIDI::Sampling::Upscatter::ModelDBRC_data::m_targetMass

The mass of the target.

Definition at line 76 of file MCGIDI_sampling.hpp.

Referenced by ModelDBRC_data(), ModelDBRC_data(), serialize(), and targetThermalSpeed().

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

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ModelDBRC_data() [1/2]

◆ ModelDBRC_data() [2/2]

◆ ~ModelDBRC_data()

Member Function Documentation

◆ crossSectionMax()

◆ evaluate()

◆ serialize()

◆ targetThermalSpeed()

Member Data Documentation

◆ m_crossSections

◆ m_domainHash

◆ m_energies

◆ m_hashIndices

◆ m_neutronMass

◆ m_targetMass