MCGIDI_reaction.cc File Reference

#include <string.h>
#include <cmath>
#include <PoPs.h>
#include "MCGIDI.h"
#include "MCGIDI_misc.h"
#include "MCGIDI_private.h"

Go to the source code of this file.

Macros
#define	nParticleChanges   6

Functions
MCGIDI_reaction *	MCGIDI_reaction_new (statusMessageReporting *smr)
int	MCGIDI_reaction_initialize (statusMessageReporting *smr, MCGIDI_reaction *reaction)
MCGIDI_reaction *	MCGIDI_reaction_free (statusMessageReporting *smr, MCGIDI_reaction *reaction)
int	MCGIDI_reaction_release (statusMessageReporting *smr, MCGIDI_reaction *reaction)
int	MCGIDI_reaction_parseFromTOM (statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_target_heated *target, MCGIDI_POPs *pops, MCGIDI_reaction *reaction)
enum MCGIDI_reactionType	MCGIDI_reaction_getReactionType (statusMessageReporting *, MCGIDI_reaction *reaction)
MCGIDI_target_heated *	MCGIDI_reaction_getTargetHeated (statusMessageReporting *, MCGIDI_reaction *reaction)
double	MCGIDI_reaction_getProjectileMass_MeV (statusMessageReporting *smr, MCGIDI_reaction *reaction)
double	MCGIDI_reaction_getTargetMass_MeV (statusMessageReporting *smr, MCGIDI_reaction *reaction)
int	MCGIDI_reaction_getDomain (statusMessageReporting *, MCGIDI_reaction *reaction, double *EMin, double *EMax)
int	MCGIDI_reaction_fixDomains (statusMessageReporting *, MCGIDI_reaction *reaction, double EMin, double EMax, nfu_status *status)
double	MCGIDI_reaction_getCrossSectionAtE (statusMessageReporting *smr, MCGIDI_reaction *reaction, MCGIDI_quantitiesLookupModes &modes, bool sampling)
double	MCGIDI_reaction_getFinalQ (statusMessageReporting *, MCGIDI_reaction *reaction, MCGIDI_quantitiesLookupModes &)
int	MCGIDI_reaction_getENDF_MTNumber (MCGIDI_reaction *reaction)
int	MCGIDI_reaction_getENDL_CSNumbers (MCGIDI_reaction *reaction, int *S)
MCGIDI_productsInfo *	MCGIDI_reaction_getProductsInfo (MCGIDI_reaction *reaction)
int	MCGIDI_reaction_recast (statusMessageReporting *smr, MCGIDI_reaction *reaction, GIDI_settings &, GIDI_settings_particle const *projectileSettings, double temperature_MeV, ptwXPoints *totalGroupedCrossSection)
int	MCGIDI_productsInfo_getNumberOfUniqueProducts (MCGIDI_productsInfo *productsInfo)
int	MCGIDI_productsInfo_getPoPsIndexAtIndex (MCGIDI_productsInfo *productsInfo, int index)
enum MCGIDI_productMultiplicityType	MCGIDI_productsInfo_getMultiplicityTypeAtIndex (MCGIDI_productsInfo *productsInfo, int index)
int	MCGIDI_productsInfo_getIntegerMultiplicityAtIndex (MCGIDI_productsInfo *productsInfo, int index)
int	MCGIDI_productsInfo_getTransportableAtIndex (MCGIDI_productsInfo *productsInfo, int index)

Macro Definition Documentation

nParticleChanges

#define nParticleChanges   6

Definition at line 18 of file MCGIDI_reaction.cc.

Function Documentation

MCGIDI_productsInfo_getIntegerMultiplicityAtIndex()

int MCGIDI_productsInfo_getIntegerMultiplicityAtIndex	(	MCGIDI_productsInfo *	productsInfo,
		int	index )

Definition at line 550 of file MCGIDI_reaction.cc.

                                                                                                      {
 
    if( ( index < 0 ) || ( index >= productsInfo->numberOfProducts ) ) return( -1 );
    return( productsInfo->productInfo[index].multiplicity );
}

MCGIDI_productsInfo_getMultiplicityTypeAtIndex()

enum MCGIDI_productMultiplicityType MCGIDI_productsInfo_getMultiplicityTypeAtIndex	(	MCGIDI_productsInfo *	productsInfo,
		int	index )

Definition at line 542 of file MCGIDI_reaction.cc.

                                                                                                                                   {
 
    if( ( index < 0 ) || ( index >= productsInfo->numberOfProducts ) ) return( MCGIDI_productMultiplicityType_invalid_e );
    return( productsInfo->productInfo[index].productMultiplicityType );
}

MCGIDI_productsInfo_getNumberOfUniqueProducts()

int MCGIDI_productsInfo_getNumberOfUniqueProducts

(

MCGIDI_productsInfo *

productsInfo

)

Definition at line 527 of file MCGIDI_reaction.cc.

                                                                                       {
 
    return( productsInfo->numberOfProducts );
}

MCGIDI_productsInfo_getPoPsIndexAtIndex()

int MCGIDI_productsInfo_getPoPsIndexAtIndex	(	MCGIDI_productsInfo *	productsInfo,
		int	index )

Definition at line 534 of file MCGIDI_reaction.cc.

                                                                                            {
 
    if( ( index < 0 ) || ( index >= productsInfo->numberOfProducts ) ) return( -1 );
    return( productsInfo->productInfo[index].globalPoPsIndex );
}

MCGIDI_productsInfo_getTransportableAtIndex()

int MCGIDI_productsInfo_getTransportableAtIndex	(	MCGIDI_productsInfo *	productsInfo,
		int	index )

Definition at line 558 of file MCGIDI_reaction.cc.

                                                                                                {
 
    if( ( index < 0 ) || ( index >= productsInfo->numberOfProducts ) ) return( -1 );
    return( productsInfo->productInfo[index].transportable );
}

MCGIDI_reaction_fixDomains()

int MCGIDI_reaction_fixDomains	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction,
		double	EMin,
		double	EMax,
		nfu_status *	status )

Definition at line 366 of file MCGIDI_reaction.cc.

                                                                                                                                      {
 
    double lowerEps = 1e-14, upperEps = -1e-14;
 
    if( reaction->EMin == EMin ) lowerEps = 0.;
    if( reaction->EMax == EMax ) upperEps = 0.;
    if( ( lowerEps == 0. ) && ( upperEps == 0. ) ) return( 0 );
 
    *status = ptwXY_dullEdges( reaction->crossSection, lowerEps, upperEps, 1 );
    return( *status != nfu_Okay );
}

Referenced by MCGIDI_target_heated_read().

MCGIDI_reaction_free()

MCGIDI_reaction * MCGIDI_reaction_free	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 60 of file MCGIDI_reaction.cc.

                                                                                                {
 
    MCGIDI_reaction_release( smr, reaction );
    smr_freeMemory( (void **) &reaction );
    return( NULL );
}

Referenced by MCGIDI_reaction_new().

MCGIDI_reaction_getCrossSectionAtE()

double MCGIDI_reaction_getCrossSectionAtE	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction,
		MCGIDI_quantitiesLookupModes &	modes,
		bool	sampling )

Definition at line 380 of file MCGIDI_reaction.cc.

                        {
 
    double e_in = modes.getProjectileEnergy( ), xsec;
 
    if( modes.getCrossSectionMode( ) == MCGIDI_quantityLookupMode_pointwise ) {
        if( e_in < reaction->EMin ) e_in = reaction->EMin;
        if( e_in > reaction->EMax ) e_in = reaction->EMax;
        ptwXY_getValueAtX( reaction->crossSection, e_in, &xsec ); }
    else if( modes.getCrossSectionMode( ) == MCGIDI_quantityLookupMode_grouped ) {
        int index = modes.getGroupIndex( );
        double *xSecP = ptwX_getPointAtIndex( reaction->crossSectionGrouped, index );
 
        if( xSecP != NULL ) {
            xsec = *xSecP;
            if( sampling && ( index == reaction->thresholdGroupIndex ) ) xsec += reaction->thresholdGroupedDeltaCrossSection; }
        else {
            xsec = 0.;
            smr_setReportError2( smr, smr_unknownID, 1, "Invalid cross section group index %d", index );
        } }
    else {
        xsec = 0.;
    }
    return( xsec );
}

Referenced by MCGIDI_target_heated_getIndexReactionCrossSectionAtE().

MCGIDI_reaction_getDomain()

int MCGIDI_reaction_getDomain	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction,
		double *	EMin,
		double *	EMax )

Definition at line 350 of file MCGIDI_reaction.cc.

                                                                                                                   {
/*
*   Return value
*       <  0    No cross section data.
*       == 0    Okay and EMin and EMax set.
*       >  0    error, EMin and EMax undefined.
*/
 
    if( !reaction->domainValuesPresent ) return( -1 );
    *EMin = reaction->EMin;
    *EMax = reaction->EMax;
    return( 0 );
}

Referenced by MCGIDI_outputChannel_getDomain(), and MCGIDI_target_heated_read().

MCGIDI_reaction_getENDF_MTNumber()

int MCGIDI_reaction_getENDF_MTNumber

(

MCGIDI_reaction *

reaction

)

Definition at line 415 of file MCGIDI_reaction.cc.

                                                                  {
 
    return( reaction->ENDF_MT );
}

Referenced by G4GIDI_target::init().

MCGIDI_reaction_getENDL_CSNumbers()

int MCGIDI_reaction_getENDL_CSNumbers	(	MCGIDI_reaction *	reaction,
		int *	S )

Definition at line 422 of file MCGIDI_reaction.cc.

                                                                           {
 
    if( S != NULL ) *S = reaction->ENDL_S;
    return( reaction->ENDL_C );
}

MCGIDI_reaction_getFinalQ()

double MCGIDI_reaction_getFinalQ	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction,
		MCGIDI_quantitiesLookupModes &	modes )

Definition at line 408 of file MCGIDI_reaction.cc.

                                                                                                                           {
 
    return( reaction->finalQ );
}

Referenced by MCGIDI_target_heated_getIndexReactionFinalQ().

MCGIDI_reaction_getProductsInfo()

MCGIDI_productsInfo * MCGIDI_reaction_getProductsInfo

(

MCGIDI_reaction *

reaction

)

Definition at line 489 of file MCGIDI_reaction.cc.

                                                                                  {
 
    return( &(reaction->productsInfo) );
}

MCGIDI_reaction_getProjectileMass_MeV()

double MCGIDI_reaction_getProjectileMass_MeV	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 336 of file MCGIDI_reaction.cc.

                                                                                                       {
 
    return( MCGIDI_target_heated_getProjectileMass_MeV( smr, reaction->target ) );
}

Referenced by MCGIDI_outputChannel_getProjectileMass_MeV(), MCGIDI_outputChannel_parseFromTOM(), and MCGIDI_outputChannel_sampleProductsAtE().

MCGIDI_reaction_getReactionType()

enum MCGIDI_reactionType MCGIDI_reaction_getReactionType	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 322 of file MCGIDI_reaction.cc.

                                                                                                                  {
 
    return( reaction->reactionType );
}

Referenced by MCGIDI_target_getReactionTypeAtIndex().

MCGIDI_reaction_getTargetHeated()

MCGIDI_target_heated * MCGIDI_reaction_getTargetHeated	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 329 of file MCGIDI_reaction.cc.

                                                                                                               {
 
    return( reaction->target );
}

Referenced by MCGIDI_outputChannel_getTargetHeated().

MCGIDI_reaction_getTargetMass_MeV()

double MCGIDI_reaction_getTargetMass_MeV	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 343 of file MCGIDI_reaction.cc.

                                                                                                   {
 
    return( MCGIDI_target_heated_getTargetMass_MeV( smr, reaction->target ) );
}

Referenced by MCGIDI_outputChannel_getTargetMass_MeV(), MCGIDI_outputChannel_parseFromTOM(), and MCGIDI_outputChannel_sampleProductsAtE().

MCGIDI_reaction_initialize()

int MCGIDI_reaction_initialize	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 42 of file MCGIDI_reaction.cc.

                                                                                         {
 
    if( MCGIDI_reaction_initialize2( smr, reaction ) != 0 ) return( 1 );
    reaction->transportabilities = new transportabilitiesMap( );
    return( 0 );
}

Referenced by MCGIDI_reaction_new(), and MCGIDI_reaction_parseFromTOM().

MCGIDI_reaction_new()

MCGIDI_reaction * MCGIDI_reaction_new

(

statusMessageReporting *

smr

)

Definition at line 31 of file MCGIDI_reaction.cc.

                                                                    {
 
    MCGIDI_reaction *reaction;
 
    if( ( reaction = (MCGIDI_reaction *) smr_malloc2( smr, sizeof( MCGIDI_reaction ), 0, "reaction" ) ) == NULL ) return( NULL );
    if( MCGIDI_reaction_initialize( smr, reaction ) ) reaction = MCGIDI_reaction_free( smr, reaction );
    return( reaction );
}

MCGIDI_reaction_parseFromTOM()

int MCGIDI_reaction_parseFromTOM	(	statusMessageReporting *	smr,
		xDataTOM_element *	element,
		MCGIDI_target_heated *	target,
		MCGIDI_POPs *	pops,
		MCGIDI_reaction *	reaction )

Definition at line 84 of file MCGIDI_reaction.cc.

                                                       {
 
    xDataTOM_element *child, *linear, *outputChannel;
    enum xDataTOM_interpolationFlag independent, dependent;
    enum xDataTOM_interpolationQualifier qualifier;
    char const *outputChannelStr, *crossSectionUnits[2] = { "MeV", "b" };
 
    MCGIDI_reaction_initialize( smr, reaction );
 
    reaction->target = target;
    reaction->reactionType = MCGIDI_reactionType_unknown_e;
    if( xDataTOME_copyAttributionList( smr, &(reaction->attributes), element ) ) goto err;
    if( xDataTOME_convertAttributeToInteger( smr, element, "ENDF_MT", &(reaction->ENDF_MT) ) ) goto err;
    if( ( outputChannelStr = xDataTOM_getAttributesValueInElement( element, "outputChannel" ) ) == NULL ) goto err;
    if( ( reaction->outputChannelStr = smr_allocateCopyString2( smr, outputChannelStr, "reaction->outputChannelStr" ) ) == NULL ) goto err;
 
    if( ( child = xDataTOME_getOneElementByName( smr, element, "crossSection", 1 ) ) == NULL ) goto err;
    if( ( linear = xDataTOME_getOneElementByName( smr, child, "linear", 0 ) ) == NULL ) {
        if( ( linear = xDataTOME_getOneElementByName( smr, child, "pointwise", 1 ) ) == NULL ) goto err;
    }
    if( xDataTOME_getInterpolation( smr, linear, 0, &independent, &dependent, &qualifier ) ) goto err;
    if( ( independent != xDataTOM_interpolationFlag_linear ) || ( dependent != xDataTOM_interpolationFlag_linear ) ) {
        smr_setReportError2( smr, smr_unknownID, 1, "cross section interpolation (%d,%d) is not linear-linear", independent, dependent );
        goto err;
    }
    if( ( reaction->crossSection = MCGIDI_misc_dataFromElement2ptwXYPointsInUnitsOf( smr, linear, crossSectionUnits ) ) == NULL ) goto err;
    reaction->domainValuesPresent = 1;
    reaction->EMin = ptwXY_getXMin( reaction->crossSection );
    reaction->EMax = ptwXY_getXMax( reaction->crossSection );
 
    if( ( outputChannel = xDataTOME_getOneElementByName( smr, element, "outputChannel", 1 ) ) == NULL ) goto err;
    if( MCGIDI_outputChannel_parseFromTOM( smr, outputChannel, pops, &(reaction->outputChannel), reaction, NULL ) ) goto err;
 
    if( MCGIDI_reaction_ParseReactionTypeAndDetermineProducts( smr, pops, reaction ) != 0 ) goto err;
 
    return( 0 );
 
err:
    MCGIDI_reaction_release( smr, reaction );
    return( 1 );
}

MCGIDI_reaction_recast()

int MCGIDI_reaction_recast	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction,
		GIDI_settings &	settings,
		GIDI_settings_particle const *	projectileSettings,
		double	temperature_MeV,
		ptwXPoints *	totalGroupedCrossSection )

Definition at line 496 of file MCGIDI_reaction.cc.

                                                                                                                         {
 
    if( totalGroupedCrossSection != NULL ) {
        nfu_status status_nf;
        GIDI_settings_group group( projectileSettings->getGroup( ) );
 
        if( reaction->crossSectionGrouped != NULL ) reaction->crossSectionGrouped = ptwX_free( reaction->crossSectionGrouped );
        if( ( reaction->crossSectionGrouped = projectileSettings->groupFunction( smr, reaction->crossSection, temperature_MeV, 0 ) ) == NULL ) return( 1 );
        if( ( status_nf = ptwX_add_ptwX( totalGroupedCrossSection, reaction->crossSectionGrouped ) ) != nfu_Okay ) return( 1 );
 
        reaction->thresholdGroupDomain = reaction->thresholdGroupedDeltaCrossSection = 0.;
        reaction->thresholdGroupIndex = group.getGroupIndexFromEnergy( reaction->EMin, false );
        if( reaction->thresholdGroupIndex > -1 ) {
            reaction->thresholdGroupDomain = group[reaction->thresholdGroupIndex+1] - reaction->EMin;
            if( reaction->thresholdGroupDomain > 0 ) {
                                                            /* factor 2 for linear reject in bin but above threshold. */
                reaction->thresholdGroupedDeltaCrossSection = *ptwX_getPointAtIndex( reaction->crossSectionGrouped, reaction->thresholdGroupIndex ) *
                        ( 2 * ( group[reaction->thresholdGroupIndex+1] - group[reaction->thresholdGroupIndex] ) / reaction->thresholdGroupDomain - 1 );
            }
        }
    }
    return( 0 );
}

Referenced by MCGIDI_target_heated_recast().

MCGIDI_reaction_release()

int MCGIDI_reaction_release	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction )

Definition at line 69 of file MCGIDI_reaction.cc.

                                                                                      {
 
    ptwXY_free( reaction->crossSection );
    ptwX_free( reaction->crossSectionGrouped );
    MCGIDI_outputChannel_release( smr, &(reaction->outputChannel) );
    xDataTOMAL_release( &(reaction->attributes) );
    smr_freeMemory( (void **) &(reaction->outputChannelStr) );
    if( reaction->productsInfo.productInfo != NULL ) smr_freeMemory( (void **) &(reaction->productsInfo.productInfo) );
    delete reaction->transportabilities;
    MCGIDI_reaction_initialize2( smr, reaction );
    return( 0 );
}

Referenced by MCGIDI_reaction_free(), MCGIDI_reaction_parseFromTOM(), and MCGIDI_target_heated_release().