#include <G4GIDI.hh>

Public Member Functions
	G4GIDI_target (PoPI::Database const &a_pops, MCGIDI::DomainHash const &a_domainHash, GIDI::Protare const &a_GIDI_protare, MCGIDI::Protare *a_MCGIDI_protare)
	~G4GIDI_target ()
std::string const *	getName () const
std::string const *	getFilename () const
std::string const *	getEvaluation () const
int	getZ () const
int	getA () const
int	getM () const
double	getMass () const
std::vector< int > const &	elasticIndices ()
std::vector< int > const &	captureIndices ()
std::vector< int > const &	fissionIndices ()
std::vector< int > const &	othersIndices ()
int	getNumberOfChannels () const
int	getNumberOfProductionChannels () const
channelID	getChannelsID (int channelIndex) const
std::vector< channelID > *	getChannelIDs () const
std::vector< channelID > *	getProductionChannelIDs () const
double	getTotalCrossSectionAtE (double a_energy, double a_temperature) const
double	getElasticCrossSectionAtE (double a_energy, double a_temperature) const
double	getCaptureCrossSectionAtE (double a_energy, double a_temperature) const
double	getFissionCrossSectionAtE (double a_energy, double a_temperature) const
double	getOthersCrossSectionAtE (double a_energy, double a_temperature) const
double	sumChannelCrossSectionAtE (std::vector< int > const &a_indices, double a_energy, double a_temperature) const
double	sumChannelCrossSectionAtE (int a_nIndices, int const *a_indices, double a_energy, double a_temperature) const
int	sampleChannelCrossSectionAtE (std::vector< int > const &a_indices, double a_energy, double a_temperature, double(a_rng)(void ), void *a_rngState) const
int	sampleChannelCrossSectionAtE (int a_nIndices, int const a_indices, double a_energy, double a_temperature, double(a_rng)(void ), void a_rngState) const
double	getElasticFinalState (double a_energy, double a_temperature, double(a_rng)(void ), void *a_rngState) const
std::vector< G4GIDI_Product > *	getCaptureFinalState (double a_energy, double a_temperature, double(a_rng)(void ), void *a_rngState) const
std::vector< G4GIDI_Product > *	getFissionFinalState (double a_energy, double a_temperature, double(a_rng)(void ), void *a_rngState) const
std::vector< G4GIDI_Product > *	getOthersFinalState (double a_energy, double a_temperature, double(a_rng)(void ), void *a_rngState) const
std::vector< G4GIDI_Product > *	getFinalState (std::vector< int > const &a_indices, double a_energy, double a_temperature, double(a_rng)(void ), void *a_rngState) const
std::vector< G4GIDI_Product > *	getFinalState (int a_nIndices, int const a_indices, double a_energy, double a_temperature, double(a_rng)(void ), void a_rngState) const

Detailed Description

A class to store map files for a particular projectile.

Definition at line 48 of file G4GIDI.hh.

Constructor & Destructor Documentation

◆ G4GIDI_target()

G4GIDI_target::G4GIDI_target	(	PoPI::Database const &	a_pops,
		MCGIDI::DomainHash const &	a_domainHash,
		GIDI::Protare const &	a_GIDI_protare,
		MCGIDI::Protare *	a_MCGIDI_protare )

Parameters

a_MCProtare [in] The MCGIDI protare.

Definition at line 37 of file G4GIDI_target.cc.

                                                  :
        m_MCGIDI_protare( a_MCGIDI_protare ),
        m_target( a_GIDI_protare.target( ).ID( ) ),
        m_fileName( a_GIDI_protare.fileName( ) ),
        m_evaluation( a_GIDI_protare.evaluation( ) ),
        m_targetZ( 0 ),
        m_targetA( 0 ),
        m_targetM( 0 ),
        m_targetMass( 0.0 ),
        m_domainHash( a_domainHash ),
        m_elasticAngular( nullptr ) {
 
    PoPI::Base const *targetAsBase = &a_pops.get<PoPI::Base const>( m_target );
    PoPI::Base const *targetAsBase2 = targetAsBase;
    if( targetAsBase->isAlias( ) ) {
        targetAsBase2 = &a_pops.get<PoPI::Base const>( a_pops.final( m_target ) );
    }
    m_targetZ = PoPI::particleZ( *targetAsBase2, true );
    m_targetA = PoPI::particleA( *targetAsBase2, true );
    m_targetM = PoPI::particleMetaStableIndex( *targetAsBase );
    if( targetAsBase2->isParticle( ) ) {
        PoPI::Particle const *targetAsParticle = static_cast<PoPI::Particle const *>( targetAsBase2 );
        m_targetMass = targetAsParticle->massValue( "amu" );
    }
 
    for( std::size_t reactionIndex = 0; reactionIndex < a_MCGIDI_protare->numberOfReactions( ); ++reactionIndex ) {
        MCGIDI::Reaction const *reaction = a_MCGIDI_protare->reaction( reactionIndex );
 
        if( reaction->ENDF_MT( ) == 2 ) {
            m_elasticIndices.push_back( static_cast<int>( reactionIndex ) ); }
        else if( reaction->ENDF_MT( ) == 18 ) {
            m_fissionIndices.push_back( static_cast<int>( reactionIndex ) ); }
        else if( reaction->ENDF_MT( ) == 102 ) {
            m_captureIndices.push_back( static_cast<int>( reactionIndex ) ); }
        else {
            m_othersIndices.push_back( static_cast<int>( reactionIndex ) );
        }
    }
 
    if( m_elasticIndices.size( ) > 0 ) {
        MCGIDI::Reaction const *elastic = a_MCGIDI_protare->reaction( m_elasticIndices[0] );
        MCGIDI::Product const *firstProduct = elastic->product( 0 );
        MCGIDI::Distributions::AngularTwoBody const *angularTwoBody = static_cast<MCGIDI::Distributions::AngularTwoBody const *>( firstProduct->distribution( ) );
        m_elasticAngular = angularTwoBody->angular( );
    }
}

◆ ~G4GIDI_target()

G4GIDI_target::~G4GIDI_target ( )

Definition at line 88 of file G4GIDI_target.cc.

                               {
 
    delete m_MCGIDI_protare;
}

Member Function Documentation

◆ captureIndices()

std::vector< int > const & G4GIDI_target::captureIndices ( )

inline

Definition at line 87 of file G4GIDI.hh.

87{ return( m_captureIndices ); }

◆ elasticIndices()

std::vector< int > const & G4GIDI_target::elasticIndices ( )

inline

Definition at line 86 of file G4GIDI.hh.

86{ return( m_elasticIndices ); }

◆ fissionIndices()

std::vector< int > const & G4GIDI_target::fissionIndices ( )

inline

Definition at line 88 of file G4GIDI.hh.

88{ return( m_fissionIndices ); }

◆ getA()

int G4GIDI_target::getA ( ) const

inline

Definition at line 76 of file G4GIDI.hh.

76{ return( m_targetA ); }

◆ getCaptureCrossSectionAtE()

double G4GIDI_target::getCaptureCrossSectionAtE	(	double	a_energy,
		double	a_temperature ) const

Definition at line 167 of file G4GIDI_target.cc.

                                                                                             {
 
    return( sumChannelCrossSectionAtE( m_captureIndices, a_energy, a_temperature ) );
}

◆ getCaptureFinalState()

std::vector< G4GIDI_Product > * G4GIDI_target::getCaptureFinalState	(	double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 257 of file G4GIDI_target.cc.

                                                                                                                                                         {
 
    return( getFinalState( m_captureIndices, a_energy, a_temperature, a_rng, a_rngState ) );
}

Referenced by G4LENDCapture::ApplyYourself().

◆ getChannelIDs()

std::vector< channelID > * G4GIDI_target::getChannelIDs ( ) const

Definition at line 119 of file G4GIDI_target.cc.

                                                          {
 
    std::vector<channelID> *channelIDs = new std::vector<channelID>( 0 );
 
    for( std::size_t reactionIndex = 0; reactionIndex < m_MCGIDI_protare->numberOfReactions( ); ++reactionIndex ) {
        MCGIDI::Reaction const *reaction = m_MCGIDI_protare->reaction( reactionIndex );
        channelIDs->push_back( reaction->label( ).c_str( ) );
    }
 
    return( channelIDs );
}

◆ getChannelsID()

channelID G4GIDI_target::getChannelsID ( int channelIndex ) const

Definition at line 111 of file G4GIDI_target.cc.

                                                               {
 
    return( channelID( m_MCGIDI_protare->reaction( channelIndex )->label( ).c_str( ) ) );
}

◆ getElasticCrossSectionAtE()

double G4GIDI_target::getElasticCrossSectionAtE	(	double	a_energy,
		double	a_temperature ) const

Definition at line 159 of file G4GIDI_target.cc.

                                                                                             {
 
    return( sumChannelCrossSectionAtE( m_elasticIndices, a_energy, a_temperature ) );
}

◆ getElasticFinalState()

double G4GIDI_target::getElasticFinalState	(	double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 249 of file G4GIDI_target.cc.

                                                                                                                                                      {
 
    return( m_elasticAngular->sample( a_energy, a_rng( a_rngState ), [&]() -> double { return a_rng( a_rngState ); } ) );
}

Referenced by G4LENDElastic::ApplyYourself(), and G4LENDModel::ApplyYourself().

◆ getEvaluation()

std::string const * G4GIDI_target::getEvaluation ( ) const

inline

Definition at line 74 of file G4GIDI.hh.

74{ return( &m_evaluation ); }

◆ getFilename()

std::string const * G4GIDI_target::getFilename ( ) const

inline

Definition at line 73 of file G4GIDI.hh.

73{ return( &m_fileName ); }

◆ getFinalState() [1/2]

std::vector< G4GIDI_Product > * G4GIDI_target::getFinalState	(	int	a_nIndices,
		int const *	a_indices,
		double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 334 of file G4GIDI_target.cc.

                                                                            {
 
    std::vector<int> indices( a_nIndices );
    for( int index = 0; index < a_nIndices; ++index ) indices[index] = a_indices[index];
 
    return( getFinalState( indices, a_energy, a_temperature, a_rng, a_rngState ) );
}

◆ getFinalState() [2/2]

std::vector< G4GIDI_Product > * G4GIDI_target::getFinalState	(	std::vector< int > const &	a_indices,
		double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 281 of file G4GIDI_target.cc.

                                                                            {
 
    int reactionIndex;
 
    if( a_indices.size( ) == 0 ) return( NULL );
    if( a_indices.size( ) == 1 ) {
        reactionIndex = a_indices[0]; }
    else {
        reactionIndex = sampleChannelCrossSectionAtE( a_indices, a_energy, a_temperature, a_rng, a_rngState );
    }
 
    MCGIDI::Sampling::StdVectorProductHandler productHandler;
    MCGIDI::Sampling::Input input( false, MCGIDI::Sampling::Upscatter::Model::none );
    input.setTemperatureAndEnergy( a_temperature, a_energy );
 
    MCGIDI::Reaction const *reaction = m_MCGIDI_protare->reaction( reactionIndex );
    reaction->sampleProducts( m_MCGIDI_protare, input, [&]() -> double { return a_rng( a_rngState ); },
        [&] (MCGIDI::Sampling::Product &a_product) -> void { productHandler.push_back( a_product ); }, productHandler );
 
    std::vector<G4GIDI_Product> *products = new std::vector<G4GIDI_Product>( productHandler.size( ) );
 
    for( std::size_t index = 0; index < productHandler.size( ); ++index ) {
        MCGIDI::Sampling::Product &productIn = productHandler[index];
        G4GIDI_Product &productOut = (*products)[index];
 
        if( productIn.m_productIntid  == 1020000000 ) {             // Neutron.
            productOut.A = 1;
            productOut.Z = 0;
            productOut.m = 0; }
        else if( productIn.m_productIntid  == 1000000000 ) {        // Photon.
            productOut.A = 0;
            productOut.Z = 0;
            productOut.m = 0; }
        else {
            PoPI::ParseIntidInfo parseIntidInfo( productIn.m_productIntid );
            productOut.A = parseIntidInfo.AAA( );
            productOut.Z = parseIntidInfo.ZZZ( );
            productOut.m = parseIntidInfo.metaStableIndex( );
        }
 
        productOut.kineticEnergy = productIn.m_kineticEnergy;
        productOut.px = productIn.m_px_vx;
        productOut.py = productIn.m_py_vy;
        productOut.pz = productIn.m_pz_vz;
    }
 
    return( products );
}

Referenced by getCaptureFinalState(), getFinalState(), getFissionFinalState(), and getOthersFinalState().

◆ getFissionCrossSectionAtE()

double G4GIDI_target::getFissionCrossSectionAtE	(	double	a_energy,
		double	a_temperature ) const

Definition at line 175 of file G4GIDI_target.cc.

                                                                                             {
 
    return( sumChannelCrossSectionAtE( m_fissionIndices, a_energy, a_temperature ) );
}

◆ getFissionFinalState()

std::vector< G4GIDI_Product > * G4GIDI_target::getFissionFinalState	(	double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 265 of file G4GIDI_target.cc.

                                                                                                                                                         {
 
    return( getFinalState( m_fissionIndices, a_energy, a_temperature, a_rng, a_rngState ) );
}

Referenced by G4LENDFission::ApplyYourself().

◆ getM()

int G4GIDI_target::getM ( ) const

inline

Definition at line 77 of file G4GIDI.hh.

77{ return( m_targetM ); }

◆ getMass()

double G4GIDI_target::getMass ( ) const

inline

Definition at line 78 of file G4GIDI.hh.

78{ return( m_targetMass ); }

◆ getName()

std::string const * G4GIDI_target::getName ( ) const

inline

Definition at line 72 of file G4GIDI.hh.

72{ return( &m_target ); }

◆ getNumberOfChannels()

int G4GIDI_target::getNumberOfChannels ( ) const

Definition at line 96 of file G4GIDI_target.cc.

                                              {
 
    return( static_cast<int>( m_MCGIDI_protare->numberOfReactions( ) ) );
}

◆ getNumberOfProductionChannels()

int G4GIDI_target::getNumberOfProductionChannels ( ) const

Definition at line 103 of file G4GIDI_target.cc.

                                                        {
 
    return( static_cast<int>( m_MCGIDI_protare->numberOfOrphanProducts( ) ) );
}

◆ getOthersCrossSectionAtE()

double G4GIDI_target::getOthersCrossSectionAtE	(	double	a_energy,
		double	a_temperature ) const

Definition at line 183 of file G4GIDI_target.cc.

                                                                                            {
 
    return( sumChannelCrossSectionAtE( m_othersIndices, a_energy, a_temperature ) );
}

◆ getOthersFinalState()

std::vector< G4GIDI_Product > * G4GIDI_target::getOthersFinalState	(	double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 273 of file G4GIDI_target.cc.

                                                                                                                                                        {
 
    return( getFinalState( m_othersIndices, a_energy, a_temperature, a_rng, a_rngState ) );
}

Referenced by G4LENDInelastic::ApplyYourself().

◆ getProductionChannelIDs()

std::vector< channelID > * G4GIDI_target::getProductionChannelIDs ( ) const

Definition at line 134 of file G4GIDI_target.cc.

                                                                    {
 
    std::vector<channelID> *channelIDs = new std::vector<channelID>( 0 );
 
    for( std::size_t reactionIndex = 0; reactionIndex < m_MCGIDI_protare->numberOfOrphanProducts( ); ++reactionIndex ) {
        MCGIDI::Reaction const *reaction = m_MCGIDI_protare->orphanProduct( reactionIndex );
        channelIDs->push_back( reaction->label( ).c_str( ) );
    }
 
    return( channelIDs );
}

◆ getTotalCrossSectionAtE()

double G4GIDI_target::getTotalCrossSectionAtE	(	double	a_energy,
		double	a_temperature ) const

Definition at line 149 of file G4GIDI_target.cc.

                                                                                           {
 
    std::size_t hashIndex = m_domainHash.index( a_energy );
 
    return( m_MCGIDI_protare->crossSection( m_URR_protareInfos, hashIndex, a_temperature, a_energy ) );
}

◆ getZ()

int G4GIDI_target::getZ ( ) const

inline

Definition at line 75 of file G4GIDI.hh.

75{ return( m_targetZ ); }

◆ othersIndices()

std::vector< int > const & G4GIDI_target::othersIndices ( )

inline

Definition at line 89 of file G4GIDI.hh.

89{ return( m_othersIndices ); }

◆ sampleChannelCrossSectionAtE() [1/2]

int G4GIDI_target::sampleChannelCrossSectionAtE	(	int	a_nIndices,
		int const *	a_indices,
		double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 237 of file G4GIDI_target.cc.

                                                                    {
 
    std::vector<int> indices( a_nIndices );
    for( int index = 0; index < a_nIndices; ++index ) indices[index] = a_indices[index];
 
    return( sampleChannelCrossSectionAtE( indices, a_energy, a_temperature, a_rng, a_rngState ) );
}

◆ sampleChannelCrossSectionAtE() [2/2]

int G4GIDI_target::sampleChannelCrossSectionAtE	(	std::vector< int > const &	a_indices,
		double	a_energy,
		double	a_temperature,
		double(*	a_rng )(void *),
		void *	a_rngState ) const

Definition at line 217 of file G4GIDI_target.cc.

                                                                    {
 
    int index = 0;
    int nIndices = static_cast<int>( a_indices.size( ) );
    double crossSectionSample = sumChannelCrossSectionAtE( a_indices, a_energy, a_temperature );
    double crossSectionSum = 0.0;
 
    crossSectionSample *= a_rng( a_rngState );
    for( ; index < nIndices; ++index ) {
        crossSectionSum += sumChannelCrossSectionAtE( 1, &a_indices[index], a_energy, a_temperature );
        if( crossSectionSum >= crossSectionSample ) break;
    }
    if( index == nIndices ) --index;          // This should really be an error.
    
    return( a_indices[index] );
}

Referenced by getFinalState(), and sampleChannelCrossSectionAtE().

◆ sumChannelCrossSectionAtE() [1/2]

double G4GIDI_target::sumChannelCrossSectionAtE	(	int	a_nIndices,
		int const *	a_indices,
		double	a_energy,
		double	a_temperature ) const

Definition at line 206 of file G4GIDI_target.cc.

                                                                                                                                   {
 
    std::vector<int> indices( a_nIndices );
    for( int index = 0; index < a_nIndices; ++index ) indices[index] = a_indices[index];
 
    return( sumChannelCrossSectionAtE( indices, a_energy, a_temperature ) );
}

◆ sumChannelCrossSectionAtE() [2/2]

double G4GIDI_target::sumChannelCrossSectionAtE	(	std::vector< int > const &	a_indices,
		double	a_energy,
		double	a_temperature ) const

Definition at line 191 of file G4GIDI_target.cc.

                                                                                                                              {
 
    std::size_t hashIndex = m_domainHash.index( a_energy );
    double crossSection = 0.0;
 
    for( auto indexIter = a_indices.begin( ); indexIter != a_indices.end( ); ++indexIter ) {
        crossSection += m_MCGIDI_protare->reactionCrossSection( *indexIter, m_URR_protareInfos, hashIndex, a_temperature, a_energy );
    }
 
    return( crossSection );
}

Referenced by getCaptureCrossSectionAtE(), getElasticCrossSectionAtE(), getFissionCrossSectionAtE(), getOthersCrossSectionAtE(), sampleChannelCrossSectionAtE(), and sumChannelCrossSectionAtE().

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ G4GIDI_target()

◆ ~G4GIDI_target()

Member Function Documentation

◆ captureIndices()

◆ elasticIndices()

◆ fissionIndices()

◆ getA()

◆ getCaptureCrossSectionAtE()

◆ getCaptureFinalState()

◆ getChannelIDs()

◆ getChannelsID()

◆ getElasticCrossSectionAtE()

◆ getElasticFinalState()

◆ getEvaluation()

◆ getFilename()

◆ getFinalState() [1/2]

◆ getFinalState() [2/2]

◆ getFissionCrossSectionAtE()

◆ getFissionFinalState()

◆ getM()

◆ getMass()

◆ getName()

◆ getNumberOfChannels()

◆ getNumberOfProductionChannels()

◆ getOthersCrossSectionAtE()

◆ getOthersFinalState()

◆ getProductionChannelIDs()

◆ getTotalCrossSectionAtE()

◆ getZ()

◆ othersIndices()

◆ sampleChannelCrossSectionAtE() [1/2]

◆ sampleChannelCrossSectionAtE() [2/2]

◆ sumChannelCrossSectionAtE() [1/2]

◆ sumChannelCrossSectionAtE() [2/2]