MCGIDI::Sampling::ProductHandler Class Reference

#include <MCGIDI_sampling.hpp>

Inheritance diagram for MCGIDI::Sampling::ProductHandler:

Public Member Functions
LUPI_HOST_DEVICE	ProductHandler ()
LUPI_HOST_DEVICE	~ProductHandler ()
template<typename RNG, typename PUSHBACK>
LUPI_HOST_DEVICE void	add (double a_projectileEnergy, int a_productIntid, int a_productIndex, int a_userProductIndex, double a_productMass, Input &a_input, RNG &&a_rng, PUSHBACK &&push_back, bool isPhoton)

Detailed Description

Definition at line 242 of file MCGIDI_sampling.hpp.

Constructor & Destructor Documentation

ProductHandler()

LUPI_HOST_DEVICE MCGIDI::Sampling::ProductHandler::ProductHandler ( )

inline

Definition at line 245 of file MCGIDI_sampling.hpp.

{}

~ProductHandler()

LUPI_HOST_DEVICE MCGIDI::Sampling::ProductHandler::~ProductHandler ( )

inline

Definition at line 246 of file MCGIDI_sampling.hpp.

{}

Member Function Documentation

add()

template<typename RNG, typename PUSHBACK>

LUPI_HOST_DEVICE void MCGIDI::Sampling::ProductHandler::add	(	double	a_projectileEnergy,
		int	a_productIntid,
		int	a_productIndex,
		int	a_userProductIndex,
		double	a_productMass,
		Input &	a_input,
		RNG &&	a_rng,
		PUSHBACK &&	push_back,
		bool	isPhoton )

Definition at line 3364 of file MCGIDI_headerSource.hpp.

                                                                                                               {
 
    Product product;
 
    if( a_isPhoton && ( a_input.m_sampledType != SampledType::unspecified ) ) a_input.m_sampledType = SampledType::photon;
 
    product.m_sampledType = a_input.m_sampledType;
    product.m_isVelocity = a_input.wantVelocity( );
    product.m_productIntid = a_productIntid;
    product.m_productIndex = a_productIndex;
    product.m_userProductIndex = a_userProductIndex;
    product.m_numberOfDBRC_rejections = a_input.m_numberOfDBRC_rejections;
    product.m_productMass = a_productMass;
 
    product.m_delayedNeutronIndex = a_input.m_delayedNeutronIndex;
    product.m_delayedNeutronDecayRate = a_input.m_delayedNeutronDecayRate;
    product.m_birthTimeSec = 0.;
    if( product.m_delayedNeutronDecayRate > 0. ) {
        product.m_birthTimeSec = -log( a_rng( ) ) / product.m_delayedNeutronDecayRate;
    }
 
    if( a_input.m_sampledType == SampledType::unspecified ) {
        product.m_kineticEnergy = 0.0;
        product.m_px_vx = 0.0;
        product.m_py_vy = 0.0;
        product.m_pz_vz = 0.0; }
    else if( a_input.m_sampledType == SampledType::uncorrelatedBody ) {
        if( a_input.m_frame == GIDI::Frame::centerOfMass ) {
            a_input.m_frame = GIDI::Frame::lab;
 
            double massRatio = a_input.m_projectileMass + a_input.m_targetMass;
            massRatio = a_input.m_projectileMass * a_productMass / ( massRatio * massRatio );
            double modifiedProjectileEnergy = massRatio * a_projectileEnergy;
 
            double sqrtModifiedProjectileEnergy = sqrt( modifiedProjectileEnergy );
            double sqrtEnergyOut_com = a_input.m_mu * sqrt( a_input.m_energyOut1 );
 
            a_input.m_energyOut1 += modifiedProjectileEnergy + 2. * sqrtModifiedProjectileEnergy * sqrtEnergyOut_com;
            if( a_input.m_energyOut1 != 0 ) a_input.m_mu = ( sqrtModifiedProjectileEnergy + sqrtEnergyOut_com ) / sqrt( a_input.m_energyOut1 );
        }
 
        product.m_kineticEnergy = a_input.m_energyOut1;
 
        double p_v = sqrt( a_input.m_energyOut1 * ( a_input.m_energyOut1 + 2. * a_productMass ) );
        if( product.m_isVelocity ) p_v *= MCGIDI_speedOfLight_cm_sec / ( a_input.m_energyOut1 + a_productMass );
 
        product.m_pz_vz = p_v * a_input.m_mu;
        p_v *= sqrt( 1. - a_input.m_mu * a_input.m_mu );
        product.m_px_vx = p_v * sin( a_input.m_phi );
        product.m_py_vy = p_v * cos( a_input.m_phi ); }
    else if( a_input.m_sampledType == SampledType::firstTwoBody ) {
        product.m_kineticEnergy = a_input.m_energyOut1;
        product.m_px_vx = a_input.m_px_vx1;
        product.m_py_vy = a_input.m_py_vy1;
        product.m_pz_vz = a_input.m_pz_vz1;
        a_input.m_sampledType = SampledType::secondTwoBody; }
    else if( a_input.m_sampledType == SampledType::secondTwoBody ) {
        product.m_kineticEnergy = a_input.m_energyOut2;
        product.m_px_vx = a_input.m_px_vx2;
        product.m_py_vy = a_input.m_py_vy2;
        product.m_pz_vz = a_input.m_pz_vz2; }
    else if( a_input.m_sampledType == SampledType::photon ) {
        product.m_kineticEnergy = a_input.m_energyOut1;
 
        double pz_vz_factor = a_input.m_energyOut1;
        if( product.m_isVelocity ) pz_vz_factor = MCGIDI_speedOfLight_cm_sec;
        product.m_pz_vz = a_input.m_mu * pz_vz_factor;
 
        double v_perp = sqrt( 1.0 - a_input.m_mu * a_input.m_mu ) * pz_vz_factor;
        product.m_px_vx = cos( a_input.m_phi ) * v_perp;
        product.m_py_vy = sin( a_input.m_phi ) * v_perp; }
    else {
        product.m_kineticEnergy = a_input.m_energyOut2;
        product.m_px_vx = a_input.m_px_vx2;
        product.m_py_vy = a_input.m_py_vy2;
        product.m_pz_vz = a_input.m_pz_vz2;
    }
 
    if( a_input.m_dataInTargetFrame && ( a_input.m_sampledType != SampledType::photon ) ) upScatterModelABoostParticle( a_input, a_rng, product );
 
    a_push_back( product );
}

Referenced by MCGIDI::ProtareSingle::sampleBranchingGammas(), MCGIDI::Product::sampleFinalState(), MCGIDI::Reaction::sampleNullProducts(), MCGIDI::GRIN_capture::sampleProducts(), MCGIDI::GRIN_inelastic::sampleProducts(), MCGIDI::OutputChannel::sampleProducts(), MCGIDI::Product::sampleProducts(), and MCGIDI::Reaction::sampleProducts().

---

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

• MCGIDI_sampling.hpp
• MCGIDI_headerSource.hpp