MCGIDI_target_heated.cc File Reference

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

Go to the source code of this file.

Functions
MCGIDI_target_heated *	MCGIDI_target_heated_free (statusMessageReporting *smr, MCGIDI_target_heated *target)
double	MCGIDI_target_heated_getIndexReactionCrossSectionAtE (statusMessageReporting *smr, MCGIDI_target_heated *target, int index, MCGIDI_quantitiesLookupModes &modes, bool sampling)
double	MCGIDI_target_heated_getIndexReactionFinalQ (statusMessageReporting *smr, MCGIDI_target_heated *target, int index, MCGIDI_quantitiesLookupModes &modes)
MCGIDI_POP *	MCGIDI_target_heated_getPOPForProjectile (statusMessageReporting *, MCGIDI_target_heated *target)
MCGIDI_POP *	MCGIDI_target_heated_getPOPForTarget (statusMessageReporting *, MCGIDI_target_heated *target)
double	MCGIDI_target_heated_getProjectileMass_MeV (statusMessageReporting *, MCGIDI_target_heated *target)
MCGIDI_reaction *	MCGIDI_target_heated_getReactionAtIndex (MCGIDI_target_heated *target, int index)
MCGIDI_reaction *	MCGIDI_target_heated_getReactionAtIndex_smr (statusMessageReporting *smr, MCGIDI_target_heated *target, int index)
int	MCGIDI_target_heated_getReactionsDomain (statusMessageReporting *, MCGIDI_target_heated *target, int index, double *EMin, double *EMax)
double	MCGIDI_target_heated_getReactionsThreshold (statusMessageReporting *, MCGIDI_target_heated *target, int index)
double	MCGIDI_target_heated_getTargetMass_MeV (statusMessageReporting *, MCGIDI_target_heated *target)
double	MCGIDI_target_heated_getTotalCrossSectionAtE (statusMessageReporting *smr, MCGIDI_target_heated *target, MCGIDI_quantitiesLookupModes &modes, bool sampling)
std::map< int, enum MCGIDI_transportability > const *	MCGIDI_target_heated_getUniqueProducts (statusMessageReporting *, MCGIDI_target_heated *target)
int	MCGIDI_target_heated_initialize (statusMessageReporting *smr, MCGIDI_target_heated *target)
MCGIDI_target_heated *	MCGIDI_target_heated_new (statusMessageReporting *smr)
MCGIDI_target_heated *	MCGIDI_target_heated_newRead (statusMessageReporting *smr, const char *fileName)
int	MCGIDI_target_heated_numberOfProductionReactions (statusMessageReporting *, MCGIDI_target_heated *)
int	MCGIDI_target_heated_numberOfReactions (statusMessageReporting *, MCGIDI_target_heated *target)
static int	MCGIDI_target_heated_parseParticle (statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, xDataTOM_element *particleAliases)
static int	MCGIDI_target_heated_parseParticleGammas (statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, char const *name)
static int	MCGIDI_target_heated_parseParticleLevel (statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, MCGIDI_POP *parent, double mass_MeV, xDataTOM_element *particleAliases)
static int	MCGIDI_target_heated_parsePOPs (statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, xDataTOM_element *particleAliases)
static int	MCGIDI_target_heated_parseReaction (statusMessageReporting *smr, xDataTOM_element *child, MCGIDI_target_heated *target, MCGIDI_POPs *pops, MCGIDI_reaction *reaction)
int	MCGIDI_target_heated_read (statusMessageReporting *smr, MCGIDI_target_heated *target, const char *fileName)
int	MCGIDI_target_heated_recast (statusMessageReporting *smr, MCGIDI_target_heated *target, GIDI_settings &settings)
int	MCGIDI_target_heated_release (statusMessageReporting *, MCGIDI_target_heated *target)
int	MCGIDI_target_heated_sampleIndexReactionProductsAtE (statusMessageReporting *smr, MCGIDI_target_heated *target, int index, MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productDatas)

Function Documentation

MCGIDI_target_heated_free()

MCGIDI_target_heated * MCGIDI_target_heated_free	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 64 of file MCGIDI_target_heated.cc.

                                                                                                             {
 
    MCGIDI_target_heated_release( smr, target );
    smr_freeMemory( (void **) &target );
    return( NULL );
}

References MCGIDI_target_heated_release(), and smr_freeMemory().

Referenced by MCGIDI_target_release().

MCGIDI_target_heated_getIndexReactionCrossSectionAtE()

double MCGIDI_target_heated_getIndexReactionCrossSectionAtE	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		int	index,
		MCGIDI_quantitiesLookupModes &	modes,
		bool	sampling
	)

Definition at line 496 of file MCGIDI_target_heated.cc.

                                                             {
 
    double xsec = 0.;
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex_smr( smr, target, index );
 
    if( reaction != NULL ) xsec = MCGIDI_reaction_getCrossSectionAtE( smr, reaction, modes, sampling );
    return( xsec );
}

References MCGIDI_reaction_getCrossSectionAtE(), and MCGIDI_target_heated_getReactionAtIndex_smr().

Referenced by MCGIDI_target_getIndexReactionCrossSectionAtE().

MCGIDI_target_heated_getIndexReactionFinalQ()

double MCGIDI_target_heated_getIndexReactionFinalQ	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		int	index,
		MCGIDI_quantitiesLookupModes &	modes
	)

Definition at line 542 of file MCGIDI_target_heated.cc.

                                              {
 
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex_smr( smr, target, index );
 
    if( reaction == NULL ) return( 0. );
    return( MCGIDI_reaction_getFinalQ( smr, reaction, modes ) );
}

References MCGIDI_reaction_getFinalQ(), and MCGIDI_target_heated_getReactionAtIndex_smr().

Referenced by MCGIDI_target_getIndexReactionFinalQ().

MCGIDI_target_heated_getPOPForProjectile()

MCGIDI_POP * MCGIDI_target_heated_getPOPForProjectile	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 433 of file MCGIDI_target_heated.cc.

                                                                                                                 {
 
    return( target->projectilePOP );
}

References MCGIDI_target_heated_s::projectilePOP.

Referenced by MCGIDI_KalbachMann_parseFromTOM().

MCGIDI_target_heated_getPOPForTarget()

MCGIDI_POP * MCGIDI_target_heated_getPOPForTarget	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 440 of file MCGIDI_target_heated.cc.

                                                                                                             {
 
    return( target->targetPOP );
}

References MCGIDI_target_heated_s::targetPOP.

Referenced by MCGIDI_KalbachMann_parseFromTOM().

MCGIDI_target_heated_getProjectileMass_MeV()

double MCGIDI_target_heated_getProjectileMass_MeV	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 447 of file MCGIDI_target_heated.cc.

                                                                                                              {
 
    return( MCGIDI_POP_getMass_MeV( target->projectilePOP ) );
}

References MCGIDI_POP_getMass_MeV(), and MCGIDI_target_heated_s::projectilePOP.

Referenced by MCGIDI_reaction_getProjectileMass_MeV().

MCGIDI_target_heated_getReactionAtIndex()

MCGIDI_reaction * MCGIDI_target_heated_getReactionAtIndex	(	MCGIDI_target_heated *	target,
		int	index
	)

Definition at line 401 of file MCGIDI_target_heated.cc.

                                                                                                    {
 
    if( ( index >= 0 ) && ( index < target->numberOfReactions ) ) return( &(target->reactions[index]) );
    return( NULL );
}

References MCGIDI_target_heated_s::reactions.

Referenced by G4GIDI_target::getChannelIDs(), G4GIDI_target::init(), MCGIDI_target_getReactionAtIndex(), MCGIDI_target_heated_getReactionAtIndex_smr(), MCGIDI_target_heated_getReactionsDomain(), MCGIDI_target_heated_getReactionsThreshold(), and MCGIDI_target_sampleReaction().

MCGIDI_target_heated_getReactionAtIndex_smr()

MCGIDI_reaction * MCGIDI_target_heated_getReactionAtIndex_smr	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		int	index
	)

Definition at line 409 of file MCGIDI_target_heated.cc.

                                                                                                                                     {
 
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex( target, index );
 
    if( reaction == NULL ) {
        smr_setReportError2( smr, smr_unknownID, 1, "bad reaction index = %d for %s + %s", index, target->projectilePOP->name, target->targetPOP->name ); 
    }
    return( reaction );
}

References MCGIDI_target_heated_getReactionAtIndex(), MCGIDI_POP_s::name, MCGIDI_target_heated_s::projectilePOP, smr_setReportError2, smr_unknownID, and MCGIDI_target_heated_s::targetPOP.

Referenced by G4GIDI_target::getChannelsID(), G4GIDI_target::getElasticFinalState(), MCGIDI_target_getReactionAtIndex_smr(), MCGIDI_target_getReactionTypeAtIndex(), MCGIDI_target_heated_getIndexReactionCrossSectionAtE(), MCGIDI_target_heated_getIndexReactionFinalQ(), and MCGIDI_target_heated_sampleIndexReactionProductsAtE().

MCGIDI_target_heated_getReactionsDomain()

int MCGIDI_target_heated_getReactionsDomain	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		int	index,
		double *	EMin,
		double *	EMax
	)

Definition at line 530 of file MCGIDI_target_heated.cc.

                                                                                                                                               {
 
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex( target, index );
 
    if( reaction == NULL ) return( -1 );
    *EMin = reaction->EMin;
    *EMax = reaction->EMax;
    return( 0 );
}

References MCGIDI_reaction_s::EMax, MCGIDI_reaction_s::EMin, and MCGIDI_target_heated_getReactionAtIndex().

Referenced by G4GIDI_target::getReactionsDomain(), and MCGIDI_target_getDomain().

MCGIDI_target_heated_getReactionsThreshold()

double MCGIDI_target_heated_getReactionsThreshold	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		int	index
	)

Definition at line 520 of file MCGIDI_target_heated.cc.

                                                                                                                         {
 
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex( target, index );
 
    if( reaction == NULL ) return( -1 );
    return( reaction->EMin );
}

References MCGIDI_reaction_s::EMin, and MCGIDI_target_heated_getReactionAtIndex().

Referenced by G4GIDI_target::getReactionsThreshold().

MCGIDI_target_heated_getTargetMass_MeV()

double MCGIDI_target_heated_getTargetMass_MeV	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 454 of file MCGIDI_target_heated.cc.

                                                                                                          {
 
    return( MCGIDI_POP_getMass_MeV( target->targetPOP ) );
}

References MCGIDI_POP_getMass_MeV(), and MCGIDI_target_heated_s::targetPOP.

Referenced by MCGIDI_reaction_getTargetMass_MeV().

MCGIDI_target_heated_getTotalCrossSectionAtE()

double MCGIDI_target_heated_getTotalCrossSectionAtE	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		MCGIDI_quantitiesLookupModes &	modes,
		bool	sampling
	)

Definition at line 461 of file MCGIDI_target_heated.cc.

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

References MCGIDI_target_heated_s::crossSection, MCGIDI_target_heated_s::crossSectionGrouped, MCGIDI_target_heated_s::crossSectionGroupedForSampling, MCGIDI_target_heated_s::EMax, MCGIDI_target_heated_s::EMin, MCGIDI_quantitiesLookupModes::getCrossSectionMode(), MCGIDI_quantitiesLookupModes::getGroupIndex(), MCGIDI_quantitiesLookupModes::getProjectileEnergy(), MCGIDI_quantityLookupMode_grouped, MCGIDI_quantityLookupMode_pointwise, ptwX_getPointAtIndex(), ptwX_length(), ptwXY_getValueAtX(), smr_setReportError2, and smr_unknownID.

Referenced by MCGIDI_target_getTotalCrossSectionAtTAndE().

MCGIDI_target_heated_getUniqueProducts()

std::map< int, enum MCGIDI_transportability > const * MCGIDI_target_heated_getUniqueProducts	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 553 of file MCGIDI_target_heated.cc.

                                                                                                                                                    {
 
    return( target->transportabilities );
}

References MCGIDI_target_heated_s::transportabilities.

Referenced by MCGIDI_target_getUniqueProducts().

MCGIDI_target_heated_initialize()

int MCGIDI_target_heated_initialize	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 43 of file MCGIDI_target_heated.cc.

                                                                                                 {
 
    memset( target, 0, sizeof( MCGIDI_target_heated ) );
    MCGIDI_POPs_initial( smr, &(target->pops), 100 );
    target->transportabilities = new transportabilitiesMap( );
    return( 0 );
}

References MCGIDI_POPs_initial(), MCGIDI_target_heated_s::pops, and MCGIDI_target_heated_s::transportabilities.

Referenced by MCGIDI_target_heated_new().

MCGIDI_target_heated_new()

MCGIDI_target_heated * MCGIDI_target_heated_new

(

statusMessageReporting *

smr

)

Definition at line 32 of file MCGIDI_target_heated.cc.

                                                                              {
 
    MCGIDI_target_heated *target;
 
    if( ( target = (MCGIDI_target_heated *) smr_malloc2( smr, sizeof( MCGIDI_target_heated ), 0, "target" ) ) == NULL ) return( NULL );
    if( MCGIDI_target_heated_initialize( smr, target ) ) target = (MCGIDI_target_heated *) smr_freeMemory( (void **) &target );
    return( target );
}

References MCGIDI_target_heated_initialize(), smr_freeMemory(), and smr_malloc2.

Referenced by MCGIDI_target_heated_newRead().

MCGIDI_target_heated_newRead()

MCGIDI_target_heated * MCGIDI_target_heated_newRead	(	statusMessageReporting *	smr,
		const char *	fileName
	)

Definition at line 53 of file MCGIDI_target_heated.cc.

                                                                                                        {
 
    MCGIDI_target_heated *target;
 
    if( ( target = MCGIDI_target_heated_new( smr ) ) == NULL ) return( NULL );
    if( MCGIDI_target_heated_read( smr, target, fileName ) != 0 ) target = (MCGIDI_target_heated *) smr_freeMemory( (void **) &target );
    return( target );
}

References MCGIDI_target_heated_new(), MCGIDI_target_heated_read(), and smr_freeMemory().

Referenced by MCGIDI_target_readHeatedTarget().

MCGIDI_target_heated_numberOfProductionReactions()

int MCGIDI_target_heated_numberOfProductionReactions	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 394 of file MCGIDI_target_heated.cc.

                                                                                                             {
 
    return( 0 );
}

Referenced by MCGIDI_target_numberOfProductionReactions().

MCGIDI_target_heated_numberOfReactions()

int MCGIDI_target_heated_numberOfReactions	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 387 of file MCGIDI_target_heated.cc.

                                                                                                       {
 
    return( target->numberOfReactions );
}

References MCGIDI_target_heated_s::numberOfReactions.

Referenced by MCGIDI_target_numberOfReactions().

MCGIDI_target_heated_parseParticle()

static int MCGIDI_target_heated_parseParticle ( statusMessageReporting *  smr, MCGIDI_target_heated *  target, xDataTOM_element *  element, xDataTOM_element *  particleAliases  )

static

Definition at line 235 of file MCGIDI_target_heated.cc.

                                            {
/*
    This routine, MCGIDI_target_heated_parseParticleLevel, MCGIDI_target_heated_parseParticleGammas handle the parsing of a
    particle which can have one of the following three forms:
 
1)
    <particle name="Sc47" genre="nucleus" mass="46.9524440292728 amu"/>
2)
    <particle name="Sc46" genre="nucleus" mass="45.9551770270807 amu">
        <level name="Sc46_e0" label="0" energy="0 eV"/>
        <level name="Sc46_e1" label="1" energy="142500 eV"/></particle>
3)
    <particle name="S36" genre="nucleus" mass="35.9669735654569 amu">
        <level name="S36_e0" label="0" energy="0 eV" spin="0"/>
        <level name="S36_e1" label="1" energy="3.291e6 eV">
            <gamma finalLevel="S36_e0" probability="1.0"/></level>
        <level name="S36_e2" label="2" energy="3.346e6 eV">
            <gamma finalLevel="S36_e0" probability="1.0"/></level>
        <level name="S36_e3" label="3" energy="4191990 eV">
            <gamma finalLevel="S36_e1" probability="1.0"/></level>
        <level name="S36_e4" label="4" energy="4522990 eV">
            <gamma finalLevel="S36_e1" probability="0.248"/>
            <gamma finalLevel="S36_e0" probability="0.752"/></level>
        <level name="S36_e5" label="5" energy="4574990 eV">
            <gamma finalLevel="S36_e1" probability="1.0"/></level>
        <level name="S36_c" label="c" energy="u:4999990 eV"/></particle>
*/
    int globalParticle = 1;
    char const *name = NULL, *mass = NULL;      /* Do not free name or mass, do not own them. */
    double mass_MeV;
    xDataTOM_element *child;
    MCGIDI_POP *pop;
 
    if( ( name = xDataTOM_getAttributesValueInElement( element, "name" ) ) == NULL ) {
        smr_setReportError2p( smr, smr_unknownID, 1, "particle missing name attribute" );
        goto err;
    }
    if( ( mass = xDataTOM_getAttributesValueInElement( element, "mass" ) ) == NULL ) {
        smr_setReportError2( smr, smr_unknownID, 1, "particle '%s' missing mass attribute", name );
        goto err;
    }
    if( MCGIDI_misc_PQUStringToDouble( smr, mass, "amu", MCGIDI_AMU2MeV, &mass_MeV ) ) goto err;
    if( ( pop = MCGIDI_POPs_addParticleIfNeeded( smr, &(target->pops), name, mass_MeV, 0., NULL, globalParticle ) ) == NULL ) goto err;
 
    for( child = xDataTOME_getFirstElement( element ); child != NULL; child = xDataTOME_getNextElement( child ) ) {
        if( strcmp( child->name, "level" ) ) {
            smr_setReportError2( smr, smr_unknownID, 1, "invalid element '%s' in %s", child->name, element->name );
            goto err;
        }
        if( MCGIDI_target_heated_parseParticleLevel( smr, target, child, pop, mass_MeV, particleAliases ) ) goto err;
    }
 
    return( 0 );
 
err:
    return( 1 );
}

References MCGIDI_AMU2MeV, MCGIDI_misc_PQUStringToDouble(), MCGIDI_POPs_addParticleIfNeeded(), MCGIDI_target_heated_parseParticleLevel(), G4InuclParticleNames::name(), xDataTOM_element_s::name, MCGIDI_target_heated_s::pops, smr_setReportError2, smr_setReportError2p, smr_unknownID, xDataTOM_getAttributesValueInElement(), xDataTOME_getFirstElement(), and xDataTOME_getNextElement().

Referenced by MCGIDI_target_heated_parsePOPs().

MCGIDI_target_heated_parseParticleGammas()

static int MCGIDI_target_heated_parseParticleGammas ( statusMessageReporting *  smr, MCGIDI_target_heated *  target, xDataTOM_element *  element, char const *  name  )

static

Definition at line 331 of file MCGIDI_target_heated.cc.

                                                                                                                                                              {
 
    int gammaCounts = 0;
    MCGIDI_POP *pop = MCGIDI_POPs_findParticle( &(target->pops), name );
    xDataTOM_element *child;
    MCGIDI_GammaBranching *gammas = NULL;
    char const *finalLevelString;
    double probability;
 
    for( child = xDataTOME_getFirstElement( element ); child != NULL; child = xDataTOME_getNextElement( child ) ) {
        if( strcmp( child->name, "gamma" ) ) {
            smr_setReportError2( smr, smr_unknownID, 1, "invalid element '%s' in %s", child->name, element->name );
            goto err;
        }
        gammaCounts++;
    }
    if( gammaCounts > 0 ) {
        if( ( gammas = (MCGIDI_GammaBranching *) smr_malloc2( smr, gammaCounts * sizeof( MCGIDI_GammaBranching), 0, "gammas" ) ) == NULL ) goto err;
        for( child = xDataTOME_getFirstElement( element ); child != NULL; child = xDataTOME_getNextElement( child ) ) {
            if( ( finalLevelString = xDataTOM_getAttributesValueInElement( child, "finalLevel" ) ) == NULL ) {
                smr_setReportError2p( smr, smr_unknownID, 1, "gamma missing 'finalLevel'" );
                goto err;
            }
            if( xDataTOME_convertAttributeToDouble( smr, child, "probability", &probability ) != 0 ) {
                smr_setReportError2p( smr, smr_unknownID, 1, "gamma missing 'probability' attribute" );
                goto err;
            }
        }
    }
    pop->numberOfGammaBranchs = gammaCounts;
    pop->gammas = gammas;
 
    return( 0 );
 
err:
    if( gammas != NULL ) smr_freeMemory( (void **) &gammas );
    return( 1 );
}

References MCGIDI_POP_s::gammas, MCGIDI_POPs_findParticle(), G4InuclParticleNames::name(), xDataTOM_element_s::name, MCGIDI_POP_s::numberOfGammaBranchs, MCGIDI_target_heated_s::pops, smr_freeMemory(), smr_malloc2, smr_setReportError2, smr_setReportError2p, smr_unknownID, xDataTOM_getAttributesValueInElement(), xDataTOME_convertAttributeToDouble(), xDataTOME_getFirstElement(), and xDataTOME_getNextElement().

Referenced by MCGIDI_target_heated_parseParticleLevel().

MCGIDI_target_heated_parseParticleLevel()

static int MCGIDI_target_heated_parseParticleLevel ( statusMessageReporting *  smr, MCGIDI_target_heated *  target, xDataTOM_element *  element, MCGIDI_POP *  parent, double  mass_MeV, xDataTOM_element *  particleAliases  )

static

Definition at line 296 of file MCGIDI_target_heated.cc.

                                                         {
 
    int globalParticle = 0;
    char const *name, *level, *aliasValue;     /* Do not free any of these as they are owned by called routine. */
    double level_MeV = 0.;
    xDataTOM_element *alias;
 
    if( ( name = xDataTOM_getAttributesValueInElement( element, "name" ) ) == NULL ) {
        smr_setReportError2p( smr, smr_unknownID, 1, "particle missing name attribute" );
        goto err;
    }
    if( ( level = xDataTOM_getAttributesValueInElement( element, "energy" ) ) == NULL ) {
        smr_setReportError2( smr, smr_unknownID, 1, "particle '%s' level missing energy attribute", name );
        goto err;
    }
        /* Special case for 'c' labels. Correct mass is only needed for two-body. */
    if( level[0] != 'u' ) if( MCGIDI_misc_PQUStringToDoubleInUnitOf( smr, level, "MeV", &level_MeV ) ) goto err;
    for( alias = xDataTOME_getFirstElement( particleAliases ); alias != NULL; alias = xDataTOME_getNextElement( alias ) ) {
        if( ( aliasValue = xDataTOM_getAttributesValueInElement( alias, "value" ) ) == NULL ) {
            smr_setReportError2p( smr, smr_unknownID, 1, "particle missing name attribute" );
            goto err;
        }
        if( strcmp( aliasValue, name ) == 0 ) globalParticle = 1;
    }
    if( MCGIDI_POPs_addParticleIfNeeded( smr, &(target->pops), name, mass_MeV + level_MeV, level_MeV, parent, globalParticle ) == NULL ) goto err;
 
    return( MCGIDI_target_heated_parseParticleGammas( smr, target, element, name ) );
 
err:
    return( 1 );
}

References MCGIDI_misc_PQUStringToDoubleInUnitOf(), MCGIDI_POPs_addParticleIfNeeded(), MCGIDI_target_heated_parseParticleGammas(), G4InuclParticleNames::name(), MCGIDI_target_heated_s::pops, smr_setReportError2, smr_setReportError2p, smr_unknownID, xDataTOM_getAttributesValueInElement(), xDataTOME_getFirstElement(), and xDataTOME_getNextElement().

Referenced by MCGIDI_target_heated_parseParticle().

MCGIDI_target_heated_parsePOPs()

static int MCGIDI_target_heated_parsePOPs ( statusMessageReporting *  smr, MCGIDI_target_heated *  target, xDataTOM_element *  element, xDataTOM_element *  particleAliases  )

static

Definition at line 215 of file MCGIDI_target_heated.cc.

                                            {
 
    xDataTOM_element *child;
 
    for( child = xDataTOME_getFirstElement( element ); child != NULL; child = xDataTOME_getNextElement( child ) ) {
        if( strcmp( child->name, "particle" ) ) {
            smr_setReportError2( smr, smr_unknownID, 1, "invalid element '%s' in %s", child->name, element->name );
            goto err;
        }
        if( MCGIDI_target_heated_parseParticle( smr, target, child, particleAliases ) ) goto err;
    }
    return( 0 );
 
err:
    return( 1 );
}

References MCGIDI_target_heated_parseParticle(), xDataTOM_element_s::name, smr_setReportError2, smr_unknownID, xDataTOME_getFirstElement(), and xDataTOME_getNextElement().

Referenced by MCGIDI_target_heated_read().

MCGIDI_target_heated_parseReaction()

static int MCGIDI_target_heated_parseReaction ( statusMessageReporting *  smr, xDataTOM_element *  child, MCGIDI_target_heated *  target, MCGIDI_POPs *  pops, MCGIDI_reaction *  reaction  )

static

Definition at line 372 of file MCGIDI_target_heated.cc.

                                                       {
 
    if( MCGIDI_reaction_parseFromTOM( smr, element, target, pops, reaction ) ) goto err;
    target->numberOfReactions++;
    
    return( 0 );
 
err:
    smr_setReportError2( smr, smr_unknownID, 1, "%s\n", xDataTOM_getAttributesValueInElement( element, "outputChannel" ) );
    return( 1 );
}

References MCGIDI_reaction_parseFromTOM(), MCGIDI_target_heated_s::numberOfReactions, smr_setReportError2, smr_unknownID, and xDataTOM_getAttributesValueInElement().

Referenced by MCGIDI_target_heated_read().

MCGIDI_target_heated_read()

int MCGIDI_target_heated_read	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		const char *	fileName
	)

Definition at line 92 of file MCGIDI_target_heated.cc.

                                                                                                                 {
/*
*   If a target has already been read into this target, user must have called MCGIDI_target_heated_release before calling this routine.
*   Otherwise, there will be memory leaks.
*/
    int n, ir;
    xDataTOM_TOM *doc = NULL;
    xDataTOM_element *element, *child, *particles, *particleAliases;
    char const *name, *version, *temperatureStr;
    char *e1;
    MCGIDI_reaction *reaction;
    double crossSectionInit[4] = { 0., 0., 0., 0., };
    nfu_status status;
    ptwXYPoints *crossSection;
int subtag1_Notice = 0;
 
    if( ( target->path = smr_allocateCopyString2( smr, fileName, "path" ) ) == NULL ) goto err;
    if( ( target->absPath = xDataTOMMisc_getAbsPath( smr, fileName ) ) == NULL ) goto err;
    if( ( doc = xDataTOM_importFile( smr, fileName ) ) == NULL ) goto err;
    element = xDataTOM_getDocumentsElement( doc );
    if( ( version = xDataTOM_getAttributesValueInElement( element, "version" ) ) == NULL ) {
            smr_setReportError2( smr, smr_unknownID, 1, "version attribute missing from element '%s'", element->name );
            goto err; }
    else {
        if( strcmp( version, "GND 1.3" ) != 0 ) {
            smr_setReportError2( smr, smr_unknownID, 1, "Unsupported version '%s' for element %s", version, element->name );
            goto err;
        }
    }
    if( strcmp( element->name, "reactionSuite" ) != 0 ) {
        smr_setReportError2( smr, smr_unknownID, 1, "input file's top element must be reactionSuite and not %s", element->name );
        goto err; }
    else {
        xDataTOMAL_copyAttributionList( smr, &(target->attributes), &(element->attributes) );
        particleAliases = xDataTOME_getOneElementByName( smr, element, "aliases", 0 );
        if( ( particles = xDataTOME_getOneElementByName( smr, element, "particles", 1 ) ) == NULL ) goto err;
        if( MCGIDI_target_heated_parsePOPs( smr, target, particles, particleAliases ) != 0 ) goto err;
 
        if( ( temperatureStr = MCGIDI_misc_pointerToTOMAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "temperature" ) ) == NULL ) goto err;
        target->temperature_MeV = strtod( temperatureStr, &e1 );
        while( isspace( *e1 ) ) ++e1; // Loop checking, 11.06.2015, T. Koi
        target->temperature_MeV *= MCGIDI_misc_getUnitConversionFactor( smr, e1, "MeV/k" );
        if( !smr_isOk( smr ) ) goto err;
 
        if( ( name = MCGIDI_misc_pointerToTOMAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "projectile" ) ) != NULL )
            target->projectilePOP = MCGIDI_POPs_findParticle( &(target->pops), name );
        if( !smr_isOk( smr ) ) goto err;
 
        if( ( name = MCGIDI_misc_pointerToTOMAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "target" ) ) != NULL )
        if( !smr_isOk( smr ) ) goto err;
            target->targetPOP = MCGIDI_POPs_findParticle( &(target->pops), name );
 
        n = xDataTOM_numberOfElementsByName( smr, element, "reaction" );
        if( n == 0 ) {
            smr_setReportError2( smr, smr_unknownID, 1, "target does not have any reactions: file = '%s'", fileName );
            goto err;
        }
        if( ( target->reactions = (MCGIDI_reaction *) smr_malloc2( smr, n * sizeof( MCGIDI_reaction ), 1, "target->reactions" ) ) == NULL ) goto err;
 
        for( ir = 0, child = xDataTOME_getFirstElement( element ); child != NULL; child = xDataTOME_getNextElement( child ) ) {
            if( strcmp( child->name, "particles" ) == 0 ) continue;
            if( strcmp( child->name, "styles" ) == 0 ) continue;
            if( strcmp( child->name, "documentations" ) == 0 ) continue;
            if( strcmp( child->name, "resonances" ) == 0 ) continue;
            if( strcmp( child->name, "summedReaction" ) == 0 ) continue;
            if( strcmp( child->name, "fissionComponent" ) == 0 ) continue;
            if( strcmp( child->name, "reaction" ) == 0 ) {
                double EMin, EMax;
 
                reaction =  &(target->reactions[ir]);
                if( MCGIDI_target_heated_parseReaction( smr, child, target, &(target->pops), reaction ) ) goto err;
                if( MCGIDI_reaction_getDomain( smr, reaction, &EMin, &EMax ) ) goto err;
                if( ir == 0 ) { target->EMin = EMin; target->EMax = EMax; }
                if( EMin < target->EMin ) target->EMin = EMin;
                if( EMax > target->EMax ) target->EMax = EMax;
                for( transportabilitiesMap::const_iterator iter = reaction->transportabilities->begin( ); 
                            iter != reaction->transportabilities->end( ); ++iter ) {
                    MCGIDI_misc_updateTransportabilitiesMap( target->transportabilities, iter->first, iter->second );
                }
                ir++; }
            else if( strcmp( child->name, "production" ) == 0 ) {
                continue; }
            else if( strcmp( child->name, "aliases" ) == 0 ) {
                continue; }
            else if( strcmp( child->name, "partialGammaProduction" ) == 0 ) {
                if( subtag1_Notice == 0 ) printf( "Unsupported reactionSuite sub-tag = '%s'\n", child->name );
                subtag1_Notice++; }
            else {
                printf( "Unsupported reactionSuite sub-tag = '%s'\n", child->name );
            }
        }
        crossSectionInit[0] = target->EMin;
        crossSectionInit[2] = target->EMax;
        if( ( target->crossSection = ptwXY_create( ptwXY_interpolationLinLin, NULL, 2., 1e-3, 2, 10, 2, crossSectionInit, &status, 0 ) ) == NULL ) {
            smr_setReportError2( smr, smr_unknownID, 1, "ptwXY_create err = %d: %s\n", status, nfu_statusMessage( status ) );
            goto err;
        }
        for( ir = 0; ir < target->numberOfReactions; ir++ ) {
            reaction =  &(target->reactions[ir]);
            if( MCGIDI_reaction_fixDomains( smr, reaction, target->EMin, target->EMax, &status ) ) {
                smr_setReportError2( smr, smr_unknownID, 1, "ptwXY_add_ptwXY err = %d: %s\n", status, nfu_statusMessage( status ) );
                goto err;
            }
            if( ( crossSection = ptwXY_add_ptwXY( target->crossSection, reaction->crossSection, &status ) ) == NULL ) {
                smr_setReportError2( smr, smr_unknownID, 1, "ptwXY_add_ptwXY err = %d: %s\n", status, nfu_statusMessage( status ) );
                goto err;
            }
            target->crossSection = ptwXY_free( target->crossSection );
            target->crossSection = crossSection;
        }
    }
    xDataTOM_freeTOM( smr, &doc );
    return( 0 );
 
err:
    smr_setReportError2( smr, smr_unknownID, 1, "Sub-error while reading file '%s'", fileName );
    if( doc != NULL ) xDataTOM_freeTOM( smr, &doc );
    MCGIDI_target_heated_release( smr, target );
    return( 1 );
}

References MCGIDI_target_heated_s::absPath, MCGIDI_target_heated_s::attributes, xDataTOM_element_s::attributes, MCGIDI_reaction_s::crossSection, MCGIDI_target_heated_s::crossSection, e1, MCGIDI_target_heated_s::EMax, MCGIDI_target_heated_s::EMin, MCGIDI_misc_getUnitConversionFactor(), MCGIDI_misc_pointerToTOMAttributeIfAllOk3, MCGIDI_misc_updateTransportabilitiesMap(), MCGIDI_POPs_findParticle(), MCGIDI_reaction_fixDomains(), MCGIDI_reaction_getDomain(), MCGIDI_target_heated_parsePOPs(), MCGIDI_target_heated_parseReaction(), MCGIDI_target_heated_release(), CLHEP::detail::n, G4InuclParticleNames::name(), xDataTOM_element_s::name, nfu_statusMessage(), MCGIDI_target_heated_s::numberOfReactions, MCGIDI_target_heated_s::path, MCGIDI_target_heated_s::pops, MCGIDI_target_heated_s::projectilePOP, ptwXY_add_ptwXY(), ptwXY_create(), ptwXY_free(), ptwXY_interpolationLinLin, MCGIDI_target_heated_s::reactions, smr_allocateCopyString2, smr_isOk(), smr_malloc2, smr_setReportError2, smr_unknownID, MCGIDI_target_heated_s::targetPOP, MCGIDI_target_heated_s::temperature_MeV, MCGIDI_reaction_s::transportabilities, MCGIDI_target_heated_s::transportabilities, xDataTOM_freeTOM(), xDataTOM_getAttributesValueInElement(), xDataTOM_getDocumentsElement(), xDataTOM_importFile(), xDataTOM_numberOfElementsByName(), xDataTOMAL_copyAttributionList(), xDataTOME_getFirstElement(), xDataTOME_getNextElement(), xDataTOME_getOneElementByName(), and xDataTOMMisc_getAbsPath().

Referenced by MCGIDI_target_heated_newRead().

MCGIDI_target_heated_recast()

int MCGIDI_target_heated_recast	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		GIDI_settings &	settings
	)

Definition at line 560 of file MCGIDI_target_heated.cc.

                                                                                                                      {
 
    int ir, projectilePoPID = target->projectilePOP->globalPoPsIndex;
    ptwXPoints *totalGroupedCrossSection = NULL;
    GIDI_settings_particle const *projectileSettings = settings.getParticle( projectilePoPID );
    nfu_status status_nf;
 
    if( projectileSettings == NULL ) {
        smr_setReportError2( smr, smr_unknownID, 1, "Settings missing for projectile %s", target->projectilePOP->name ); 
        return( 1 );
    }
    target->crossSectionGrouped = ptwX_free( target->crossSectionGrouped );
    target->crossSectionGroupedForSampling = ptwX_free( target->crossSectionGroupedForSampling );
    if( projectileSettings->isEnergyMode_grouped( ) ) {
        int64_t numberOfGroups = projectileSettings->getNumberOfGroups( );
 
        if( ( totalGroupedCrossSection = ptwX_createLine( numberOfGroups, numberOfGroups, 0, 0, &status_nf ) ) == NULL ) {
            smr_setReportError2( smr, smr_unknownID, 1, "totalGroupedCrossSection allocation failed: status_nf = %d, '%s'", 
                    status_nf, nfu_statusMessage( status_nf ) ); 
            goto err;
        }
    }
 
    for( ir = 0; ir < target->numberOfReactions; ++ir ) {
        if( MCGIDI_reaction_recast( smr, &(target->reactions[ir]), settings, projectileSettings, target->temperature_MeV, totalGroupedCrossSection ) != 0 ) goto err;
    }
    if( projectileSettings->isEnergyMode_grouped( ) ) {
        if( ( target->crossSectionGroupedForSampling = ptwX_clone( totalGroupedCrossSection, &status_nf ) ) == NULL ) {
            smr_setReportError2( smr, smr_unknownID, 1, "totalGroupedCrossSection allocation failed: status_nf = %d, '%s'", 
                    status_nf, nfu_statusMessage( status_nf ) ); 
            goto err;
        }
        for( ir = 0; ir < target->numberOfReactions; ++ir ) {
            int index = target->reactions[ir].thresholdGroupIndex;
 
            if( index > -1 ) {
                double xSec = target->reactions[ir].thresholdGroupedDeltaCrossSection +
                        ptwX_getPointAtIndex_Unsafely( target->crossSectionGroupedForSampling, index );
 
                ptwX_setPointAtIndex( target->crossSectionGroupedForSampling, index, xSec );
            }
        }
    }
    target->crossSectionGrouped = totalGroupedCrossSection;
    totalGroupedCrossSection = NULL;
 
    return( 0 );
 
err:
    ptwX_free( totalGroupedCrossSection );
    target->crossSectionGroupedForSampling = ptwX_free( target->crossSectionGroupedForSampling );
    return( 1 );
}

References MCGIDI_target_heated_s::crossSectionGrouped, MCGIDI_target_heated_s::crossSectionGroupedForSampling, GIDI_settings_particle::getNumberOfGroups(), GIDI_settings::getParticle(), MCGIDI_POP_s::globalPoPsIndex, GIDI_settings_particle::isEnergyMode_grouped(), MCGIDI_reaction_recast(), MCGIDI_POP_s::name, nfu_statusMessage(), MCGIDI_target_heated_s::numberOfReactions, MCGIDI_target_heated_s::projectilePOP, ptwX_clone(), ptwX_createLine(), ptwX_free(), ptwX_getPointAtIndex_Unsafely(), ptwX_setPointAtIndex(), MCGIDI_target_heated_s::reactions, smr_setReportError2, smr_unknownID, MCGIDI_target_heated_s::temperature_MeV, MCGIDI_reaction_s::thresholdGroupedDeltaCrossSection, and MCGIDI_reaction_s::thresholdGroupIndex.

Referenced by MCGIDI_target_recast().

MCGIDI_target_heated_release()

int MCGIDI_target_heated_release	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target
	)

Definition at line 73 of file MCGIDI_target_heated.cc.

                                                                                             {
 
    int ir;
 
    ptwXY_free( target->crossSection );
    ptwX_free( target->crossSectionGrouped );
    ptwX_free( target->crossSectionGroupedForSampling );
    for( ir = 0; ir < target->numberOfReactions; ir++ ) MCGIDI_reaction_release( NULL, &(target->reactions[ir]) );
    smr_freeMemory( (void **) &(target->reactions) );
    MCGIDI_POPs_release( &(target->pops) );
    smr_freeMemory( (void **) &(target->path) );
    smr_freeMemory( (void **) &(target->absPath) );
    xDataTOMAL_release( &(target->attributes) );
    delete target->transportabilities;
    return( 0 );
}

References MCGIDI_target_heated_s::absPath, MCGIDI_target_heated_s::attributes, MCGIDI_target_heated_s::crossSection, MCGIDI_target_heated_s::crossSectionGrouped, MCGIDI_target_heated_s::crossSectionGroupedForSampling, MCGIDI_POPs_release(), MCGIDI_reaction_release(), MCGIDI_target_heated_s::numberOfReactions, MCGIDI_target_heated_s::path, MCGIDI_target_heated_s::pops, ptwX_free(), ptwXY_free(), MCGIDI_target_heated_s::reactions, smr_freeMemory(), MCGIDI_target_heated_s::transportabilities, and xDataTOMAL_release().

Referenced by MCGIDI_target_heated_free(), and MCGIDI_target_heated_read().

MCGIDI_target_heated_sampleIndexReactionProductsAtE()

int MCGIDI_target_heated_sampleIndexReactionProductsAtE	(	statusMessageReporting *	smr,
		MCGIDI_target_heated *	target,
		int	index,
		MCGIDI_quantitiesLookupModes &	modes,
		MCGIDI_decaySamplingInfo *	decaySamplingInfo,
		MCGIDI_sampledProductsDatas *	productDatas
	)

Definition at line 508 of file MCGIDI_target_heated.cc.

                                                                                                                                      {
 
    MCGIDI_reaction *reaction = MCGIDI_target_heated_getReactionAtIndex_smr( smr, target, index );
 
    productDatas->numberOfProducts = 0;
    if( reaction == NULL ) return( -1 );
    return( MCGIDI_outputChannel_sampleProductsAtE( smr, &(reaction->outputChannel), modes, decaySamplingInfo, productDatas, NULL ) );
}

References MCGIDI_outputChannel_sampleProductsAtE(), MCGIDI_target_heated_getReactionAtIndex_smr(), MCGIDI_sampledProductsDatas_s::numberOfProducts, and MCGIDI_reaction_s::outputChannel.

Referenced by G4GIDI_target::getFinalState(), and MCGIDI_target_sampleIndexReactionProductsAtE().