MCGIDI_reaction.cc File Reference

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

Include dependency graph for MCGIDI_reaction.cc:

Go to the source code of this file.

Macros
#define	nParticleChanges   6

Functions
static int	MCGIDI_reaction_initialize2 (statusMessageReporting *smr, MCGIDI_reaction *reaction)
static int	MCGIDI_reaction_particleChanges (MCGIDI_POP *projectile, MCGIDI_POP *target, MCGIDI_productsInfo *productsInfo, int n1, int *particlesChanges)
static int	MCGIDI_reaction_ParseReactionTypeAndDetermineProducts (statusMessageReporting *smr, MCGIDI_POPs *pops, MCGIDI_reaction *reaction)
static int	MCGIDI_reaction_ParseDetermineReactionProducts (statusMessageReporting *smr, MCGIDI_POPs *pops, MCGIDI_outputChannel *outputChannel, MCGIDI_productsInfo *productsInfo, MCGIDI_reaction *reaction, double *finalQ, int level)
static int	MCGIDI_reaction_addReturnProduct (statusMessageReporting *smr, MCGIDI_productsInfo *productsInfo, int ID, MCGIDI_product *product, MCGIDI_reaction *reaction, int transportable)
static int	MCGIDI_reaction_setENDL_CSNumbers (statusMessageReporting *smr, MCGIDI_reaction *reaction)
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_addReturnProduct()

static int MCGIDI_reaction_addReturnProduct ( statusMessageReporting *  smr, MCGIDI_productsInfo *  productsInfo, int  ID, MCGIDI_product *  product, MCGIDI_reaction *  reaction, int  transportable  )

static

Definition at line 279 of file MCGIDI_reaction.cc.

                                                       {

    int i1;
    enum MCGIDI_productMultiplicityType productMultiplicityType;

    MCGIDI_misc_updateTransportabilitiesMap2( reaction->transportabilities, ID, transportable );
    for( i1 = 0; i1 < productsInfo->numberOfProducts; i1++ ) {
        if( productsInfo->productInfo[i1].globalPoPsIndex == ID ) break;
    }
    if( i1 == productsInfo->numberOfProducts ) {
        if( productsInfo->numberOfProducts == productsInfo->numberOfAllocatedProducts ) {
            productsInfo->numberOfAllocatedProducts += 4;
            if( ( productsInfo->productInfo = (MCGIDI_productInfo *) smr_realloc2( smr, productsInfo->productInfo, 
                productsInfo->numberOfAllocatedProducts * sizeof( MCGIDI_productInfo ), "productsInfo->productInfo" ) ) == NULL ) return( 1 );
        }
        productsInfo->numberOfProducts++;
        productsInfo->productInfo[i1].globalPoPsIndex = ID;
        productsInfo->productInfo[i1].productMultiplicityType = MCGIDI_productMultiplicityType_unknown_e;
        productsInfo->productInfo[i1].multiplicity = 0;
        productsInfo->productInfo[i1].transportable = transportable;
    }
    if( product == NULL ) {
        productMultiplicityType = MCGIDI_productMultiplicityType_gammaBranching_e; }
    else {
        if( ( product->multiplicityVsEnergy != NULL ) || ( product->piecewiseMultiplicities != NULL ) ) {
            productMultiplicityType = MCGIDI_productMultiplicityType_energyDependent_e; }
        else {
            productsInfo->productInfo[i1].multiplicity += product->multiplicity;
            productMultiplicityType = MCGIDI_productMultiplicityType_integer_e;
        }
    }
    if( ( productsInfo->productInfo[i1].productMultiplicityType == MCGIDI_productMultiplicityType_unknown_e ) ||
        ( productsInfo->productInfo[i1].productMultiplicityType == productMultiplicityType ) ) {
        productsInfo->productInfo[i1].productMultiplicityType = productMultiplicityType; }
    else {
        productsInfo->productInfo[i1].productMultiplicityType = MCGIDI_productMultiplicityType_mixed_e;
    }
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_fixDomains()

int MCGIDI_reaction_fixDomains	(	statusMessageReporting *	,
		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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_getDomain()

int MCGIDI_reaction_getDomain	(	statusMessageReporting *	,
		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 );
}

Here is the caller graph for this function:

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 );
}

Here is the caller graph for this function:

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 *	,
		MCGIDI_reaction *	reaction,
		MCGIDI_quantitiesLookupModes &
	)

Definition at line 408 of file MCGIDI_reaction.cc.

                                                                                                                           {

    return( reaction->finalQ );
}

Here is the caller graph for this function:

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 ) );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_getReactionType()

enum MCGIDI_reactionType MCGIDI_reaction_getReactionType	(	statusMessageReporting *	,
		MCGIDI_reaction *	reaction
	)

Definition at line 322 of file MCGIDI_reaction.cc.

                                                                                                                  {

    return( reaction->reactionType );
}

Here is the caller graph for this function:

MCGIDI_reaction_getTargetHeated()

MCGIDI_target_heated* MCGIDI_reaction_getTargetHeated	(	statusMessageReporting *	,
		MCGIDI_reaction *	reaction
	)

Definition at line 329 of file MCGIDI_reaction.cc.

                                                                                                               {

    return( reaction->target );
}

Here is the caller graph for this function:

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 ) );
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_initialize2()

static int MCGIDI_reaction_initialize2 ( statusMessageReporting *  smr, MCGIDI_reaction *  reaction  )

static

Definition at line 51 of file MCGIDI_reaction.cc.

                                                                                                 {

    memset( reaction, 0, sizeof( MCGIDI_reaction ) );
    xDataTOMAL_initial( smr, &(reaction->attributes) );
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 );
}

Here is the call graph for this function:

MCGIDI_reaction_ParseDetermineReactionProducts()

static int MCGIDI_reaction_ParseDetermineReactionProducts ( statusMessageReporting *  smr, MCGIDI_POPs *  pops, MCGIDI_outputChannel *  outputChannel, MCGIDI_productsInfo *  productsInfo, MCGIDI_reaction *  reaction, double *  finalQ, int  level  )

static

Definition at line 202 of file MCGIDI_reaction.cc.

                                                                                              {
/*
*   This function determines all products that can be returned during sampling for this outputChannel. Note, products like 'U238_c' and
*   'U238_e3' are not returned during sampling as both are decay to the groud state (unless a meta-stable is encountered). 
*   Some examples for projectile 'n' and target 'U238' are:
*       outputChannel                       products returned during sampling.
*       'n + U238'                          n, U238
*       'n + U238 + gamma'                  n, U238, gamma
*       'n + U238_c'                        n, U238     (even if no gammas are give, return ground state of residual.
*       'n + (U238_e3 -> U238 + gamma)'     n, U238, gamma
*/
    int iProduct, nProducts = MCGIDI_outputChannel_numberOfProducts( outputChannel ), globalPoPsIndex, productIsTrackable;
    int twoBodyProductsWithData = 0;
    MCGIDI_product *product;
    MCGIDI_POP *residual;

    if( ( level == 0 ) && ( outputChannel->genre == MCGIDI_channelGenre_twoBody_e ) ) {
        for( iProduct = 0; iProduct < nProducts; iProduct++ ) {
            product = MCGIDI_outputChannel_getProductAtIndex( smr, outputChannel, iProduct );
            if( product->pop->globalPoPsIndex < 0 ) {
                twoBodyProductsWithData = -1; }
            else if( product->distribution.type == MCGIDI_distributionType_angular_e ) {
                if( twoBodyProductsWithData >= 0 ) twoBodyProductsWithData = 1;
            }
        }
    }
    if( twoBodyProductsWithData < 0 ) twoBodyProductsWithData = 0;
    *finalQ += MCGIDI_outputChannel_getQ_MeV( smr, outputChannel, 0 );
    for( iProduct = 0; iProduct < nProducts; iProduct++ ) {
        productIsTrackable = twoBodyProductsWithData;
        product = MCGIDI_outputChannel_getProductAtIndex( smr, outputChannel, iProduct );
        globalPoPsIndex = product->pop->globalPoPsIndex;
        if( ( product->distribution.type != MCGIDI_distributionType_none_e ) && ( product->distribution.type != MCGIDI_distributionType_unknown_e ) ) {
            productIsTrackable = 1;
            if( globalPoPsIndex < 0 ) {
                if( product->distribution.angular != NULL ) {
                    if( product->distribution.angular->type == MCGIDI_angularType_recoil ) productIsTrackable = 0;
                }
                if( productIsTrackable ) {
                    int len = (int) strlen( product->pop->name );

                    if( len > 2 ) {     /* Special case for continuum reactions with data for residual (e.g., n + U233 -> n + U233_c). */
                        if( ( product->pop->name[len-2] == '_' ) && ( product->pop->name[len-1] == 'c' ) ) {
                            for( residual = product->pop; residual->globalPoPsIndex < 0; residual = residual->parent ) ;
                            productIsTrackable = 1;
                            globalPoPsIndex = residual->globalPoPsIndex;
                        }
                    }
                    if( globalPoPsIndex < 0 ) {
                        smr_setReportError2( smr, smr_unknownID, 1, "product determination for '%s' cannot be determined", product->pop->name );
                        return( 1 );
                    }
                }
            }
        }
        if( productIsTrackable ) {
            if( MCGIDI_reaction_addReturnProduct( smr, productsInfo, globalPoPsIndex, product, reaction, 1 ) != 0 ) return( 1 ); }
        else {
            if( product->decayChannel.genre != MCGIDI_channelGenre_undefined_e ) {
                if( MCGIDI_reaction_ParseDetermineReactionProducts( smr, pops, &(product->decayChannel), productsInfo, reaction, finalQ, level + 1 ) != 0 ) return( 1 ); }
            else {
                *finalQ += product->pop->level_MeV;
                for( residual = product->pop; residual->globalPoPsIndex < 0; residual = residual->parent ) ;
                if( MCGIDI_reaction_addReturnProduct( smr, productsInfo, residual->globalPoPsIndex, product, reaction, 0 ) != 0 ) return( 1 );
                if( product->pop->numberOfGammaBranchs != 0 ) {
                    int gammaIndex = PoPs_particleIndex( "gamma" );
                    if( MCGIDI_reaction_addReturnProduct( smr, productsInfo, gammaIndex, NULL, reaction, 1 ) != 0 ) return( 1 );
                }
            }
        }
    }
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_ParseReactionTypeAndDetermineProducts()

static int MCGIDI_reaction_ParseReactionTypeAndDetermineProducts ( statusMessageReporting *  smr, MCGIDI_POPs *  pops, MCGIDI_reaction *  reaction  )

static

Definition at line 129 of file MCGIDI_reaction.cc.

                                                                                                                                              {

    MCGIDI_outputChannel *outputChannel = &(reaction->outputChannel);
    int MT;
    int particlesChanges[nParticleChanges], numberOfChanges;
    double finalQ = 0.;

    if( MCGIDI_reaction_ParseDetermineReactionProducts( smr, pops, outputChannel, &(reaction->productsInfo), reaction, &finalQ, 0 ) != 0 ) return( 1 );
    reaction->finalQ = finalQ;
    MT = MCGIDI_reaction_getENDF_MTNumber( reaction );
    switch( MT ) {
    case 2 :
        reaction->reactionType = MCGIDI_reactionType_elastic_e;
        break;
    case 18 : case 19 : case 20 : case 21 : case 38 :
        reaction->reactionType = MCGIDI_reactionType_fission_e;
        break;
    case 102 :
        reaction->reactionType = MCGIDI_reactionType_capture_e;
        break;
    case 5 :
        reaction->reactionType = MCGIDI_reactionType_sumOfRemainingOutputChannels_e;
        break;
    default :
        numberOfChanges = MCGIDI_reaction_particleChanges( reaction->target->projectilePOP, reaction->target->targetPOP, &(reaction->productsInfo), 
            nParticleChanges, particlesChanges );

        reaction->reactionType = MCGIDI_reactionType_unknown_e;
        if( numberOfChanges == 0 ) {
            reaction->reactionType = MCGIDI_reactionType_scattering_e; }
        else {
            reaction->reactionType = MCGIDI_reactionType_nuclearIsomerTransmutation_e;
        }

/*
    Currently, these are not handled properly:
        MCGIDI_reactionType_nuclearLevelTransition_e
        MCGIDI_reactionType_atomic_e
*/
        break;
    }

    MCGIDI_reaction_setENDL_CSNumbers( smr, reaction );
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_particleChanges()

static int MCGIDI_reaction_particleChanges ( MCGIDI_POP *  projectile, MCGIDI_POP *  target, MCGIDI_productsInfo *  productsInfo, int  n1, int *  particlesChanges  )

static

Definition at line 177 of file MCGIDI_reaction.cc.

                                                                                                                                                           {

    int projectileGlobalIndex = projectile->globalPoPsIndex, targetGlobalIndex = target->globalPoPsIndex, i1, i2 = 0;
    int gammaIndex = PoPs_particleIndex( "gamma" );

    if( projectileGlobalIndex != gammaIndex ) {
        for( i1 = 0; i1 < productsInfo->numberOfProducts; i1++ ) if( projectileGlobalIndex == productsInfo->productInfo[i1].globalPoPsIndex ) break;
        if( i1 == productsInfo->numberOfProducts ) particlesChanges[i2++] = projectileGlobalIndex;
    }

    for( i1 = 0; i1 < productsInfo->numberOfProducts; i1++ ) if( targetGlobalIndex == productsInfo->productInfo[i1].globalPoPsIndex ) break;
    if( i1 == productsInfo->numberOfProducts ) particlesChanges[i2++] = targetGlobalIndex;

    for( i1 = 0; i1 < productsInfo->numberOfProducts; i1++ ) {
        if( i2 == n1 ) break;
        if( /*(*/ projectileGlobalIndex == productsInfo->productInfo[i1].globalPoPsIndex /*)*/ ) continue;
        if( /*(*/ targetGlobalIndex == productsInfo->productInfo[i1].globalPoPsIndex /*)*/ ) continue;
        if( /*(*/ gammaIndex == productsInfo->productInfo[i1].globalPoPsIndex /*)*/ ) continue;
        particlesChanges[i2++] = productsInfo->productInfo[i1].globalPoPsIndex;
    }
    return( i2 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_recast()

int MCGIDI_reaction_recast	(	statusMessageReporting *	smr,
		MCGIDI_reaction *	reaction,
		GIDI_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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MCGIDI_reaction_setENDL_CSNumbers()

static int MCGIDI_reaction_setENDL_CSNumbers ( statusMessageReporting *  smr, MCGIDI_reaction *  reaction  )

static

Definition at line 430 of file MCGIDI_reaction.cc.

                                                                                                      {

    int MT = MCGIDI_reaction_getENDF_MTNumber( reaction );
    int MT1_50ToC[] = { 1,   10,  -3,   -4,   -5,    0,    0,    0,    0,  -10,
                       32,    0,   0,    0,    0,   12,   13,   15,   15,   15,
                       15,   26,  36,   33,  -25,    0,  -27,   20,   27,  -30,
                        0,   22,  24,   25,  -35,  -36,   14,   15,    0,    0,
                       29,   16,   0,   17,   34,    0,    0,    0,    0 };
    int MT100_200ToC[] = { -101,   46,   40,   41,   42,   44,   45,   37, -109,    0,
                             18,   48, -113, -114,   19,   39,   47,    0,    0,    0,
                              0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
                              0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
                              0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
                              0, -152, -153, -154,   43, -156, -157,   23,   31, -160,
                           -161, -162, -163, -164, -165, -166, -167, -168, -169, -170,
                           -171, -172, -173, -174, -175, -176, -177, -178, -179, -180,
                           -181, -182, -183, -184, -185, -186, -187, -188,   28, -190,
                           -191, -192,   38, -194, -195, -196, -197, -198, -199, -200 };

    reaction->ENDL_C = 0;
    reaction->ENDL_S = 0;
    if( MT <= 0 ) return( 1 );
    if( MT > 891 ) return( 1 );
    if( MT < 50 ) {
        reaction->ENDL_C = MT1_50ToC[MT - 1]; }
    else if( MT <= 91 ) {
        reaction->ENDL_C = 11;
        if( MT != 91 ) reaction->ENDL_S = 1; }
    else if( ( MT > 100 ) && ( MT <= 200 ) ) {
        reaction->ENDL_C = MT100_200ToC[MT - 101]; }
    else if( ( MT == 452 ) || ( MT == 455 ) || ( MT == 456 ) || ( MT == 458 ) ) {
        reaction->ENDL_C = 15;
        if( MT == 455 ) reaction->ENDL_S = 7; }
    else if( MT >= 600 ) {
        if( MT < 650 ) {
            reaction->ENDL_C = 40;
            if( MT != 649 ) reaction->ENDL_S = 1; }
        else if( MT < 700 ) {
            reaction->ENDL_C = 41;
            if( MT != 699 ) reaction->ENDL_S = 1; }
        else if( MT < 750 ) {
            reaction->ENDL_C = 42;
            if( MT != 749 ) reaction->ENDL_S = 1; }
        else if( MT < 800 ) {
            reaction->ENDL_C = 44;
            if( MT != 799 ) reaction->ENDL_S = 1; }
        else if( MT < 850 ) {
            reaction->ENDL_C = 45;
            if( MT != 849 ) reaction->ENDL_S = 1; }
        else if( ( MT >= 875 ) && ( MT <= 891 ) ) {
            reaction->ENDL_C = 12;
            if( MT != 891 ) reaction->ENDL_S = 1;
        }
    }
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function: