MCGIDI_misc.cc File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <cmath>
#include <ctype.h>
#include <unistd.h>
#include <ptwXY.h>
#include <xDataTOM_importXML_private.h>
#include "MCGIDI.h"
#include "MCGIDI_misc.h"
#include "MCGIDI_fromTOM.h"

Go to the source code of this file.

Data Structures
struct	ZSymbol

Functions
xDataTOM_Int	MCGIDI_misc_binarySearch (xDataTOM_Int n, double *ds, double d)
int	MCGIDI_misc_copyXMLAttributesToTOM (statusMessageReporting *smr, xDataTOM_attributionList *TOM, xDataXML_attributionList *XML)
static ptwXYPoints *	MCGIDI_misc_Data2ptwXYPointsInUnitsOf (statusMessageReporting *smr, ptwXY_interpolation interpolation, int length, double *data, char const *fromUnits[2], char const *toUnits[2])
ptwXYPoints *	MCGIDI_misc_dataFromElement2ptwXYPointsInUnitsOf (statusMessageReporting *smr, xDataTOM_element *linear, char const *toUnits[2])
ptwXYPoints *	MCGIDI_misc_dataFromXYs2ptwXYPointsInUnitsOf (statusMessageReporting *smr, xDataTOM_XYs *XYs, ptwXY_interpolation interpolation, char const *toUnits[2])
char *	MCGIDI_misc_getAbsPath (statusMessageReporting *smr, const char *fileName)
enum xDataTOM_frame	MCGIDI_misc_getProductFrame (statusMessageReporting *smr, xDataTOM_element *frameElement)
double	MCGIDI_misc_getUnitConversionFactor (statusMessageReporting *smr, char const *fromUnit, char const *toUnit)
int	MCGIDI_misc_NumberOfZSymbols (void)
char const *	MCGIDI_misc_pointerToAttributeIfAllOk (statusMessageReporting *smr, xDataXML_element *element, const char *path, int required, xDataTOM_attributionList *attributes, const char *name, const char *file, int line)
char const *	MCGIDI_misc_pointerToTOMAttributeIfAllOk (statusMessageReporting *smr, const char *path, int required, xDataTOM_attributionList *attributes, const char *name, const char *file, int line)
int	MCGIDI_misc_PQUStringToDouble (statusMessageReporting *smr, char const *str, char const *unit, double conversion, double *value)
int	MCGIDI_misc_PQUStringToDoubleInUnitOf (statusMessageReporting *smr, char const *str, char const *toUnit, double *value)
int	MCGIDI_misc_setMessageError_Element (statusMessageReporting *smr, void *userInterface, xDataXML_element *element, const char *file, int line, int code, const char *fmt,...)
int	MCGIDI_misc_symbolToZ (const char *Z)
void	MCGIDI_misc_updateTransportabilitiesMap (transportabilitiesMap *transportabilities, int PoPID, enum MCGIDI_transportability transportability)
void	MCGIDI_misc_updateTransportabilitiesMap2 (transportabilitiesMap *transportabilities, int PoPID, int transportable)
const char *	MCGIDI_misc_ZToSymbol (int iZ)
int	MCGIDI_miscNameToZAm (statusMessageReporting *smr, const char *name, int *Z, int *A, int *m, int *level)
static int	MCGIDI_miscNameToZAm_getLevel (statusMessageReporting *smr, const char *name, const char *p)

Variables
static struct ZSymbol	ZSymbols []

Function Documentation

MCGIDI_misc_binarySearch()

xDataTOM_Int MCGIDI_misc_binarySearch	(	xDataTOM_Int	n,
		double *	ds,
		double	d
	)

Definition at line 228 of file MCGIDI_misc.cc.

                                                                              {
/*
*   Returns -2 is d < first point of ds, -1 if > last point of ds and the lower index of ds otherwise.
*/
    xDataTOM_Int imin = 0, imid, imax = n - 1;
 
    if( d < ds[0] ) return( -2 );
    if( d > ds[n-1] ) return( -1 );
    while( 1 ) { // Loop checking, 11.06.2015, T. Koi
        imid = ( imin + imax ) >> 1;
        if( imid == imin ) break;
        if( d < ds[imid] ) {
            imax = imid; }
        else {
            imin = imid;
        }
    }
    return( imin );
}

References anonymous_namespace{G4QuasiElRatios.cc}::ds, imax, and CLHEP::detail::n.

Referenced by MCGIDI_sampling_doubleDistribution(), MCGIDI_sampling_sampleX_from_pdfOfX(), and MCGIDI_sampling_sampleX_from_pdfsOfXGivenW().

MCGIDI_misc_copyXMLAttributesToTOM()

int MCGIDI_misc_copyXMLAttributesToTOM	(	statusMessageReporting *	smr,
		xDataTOM_attributionList *	TOM,
		xDataXML_attributionList *	XML
	)

Definition at line 296 of file MCGIDI_misc.cc.

                                                                                                                                    {
 
    int i;
    xDataXML_attribute *attribute;
 
    xDataTOMAL_initial( smr, TOM );
    for( i = 0; ; i++ ) {
        if( ( attribute = xDataXML_attributeByIndex( XML, i ) ) == NULL ) break;
        if( xDataTOMAL_addAttribute( smr, TOM, attribute->name, attribute->value ) != 0 ) goto err;
    }
    return( 0 );
 
err:
    xDataTOMAL_release( TOM );
    return( 1 );
}

References xDataXML_attribute_s::name, xDataXML_attribute_s::value, xDataTOMAL_addAttribute(), xDataTOMAL_initial(), xDataTOMAL_release(), and xDataXML_attributeByIndex().

MCGIDI_misc_Data2ptwXYPointsInUnitsOf()

static ptwXYPoints * MCGIDI_misc_Data2ptwXYPointsInUnitsOf ( statusMessageReporting *  smr, ptwXY_interpolation  interpolation, int  length, double *  data, char const *  fromUnits[2], char const *  toUnits[2]  )

static

Definition at line 447 of file MCGIDI_misc.cc.

                                                                                     {
 
    double xFactor, yFactor;
    ptwXYPoints *ptwXY = NULL;
    nfu_status status;
 
    xFactor = MCGIDI_misc_getUnitConversionFactor( smr, fromUnits[0], toUnits[0] );
    if( !smr_isOk( smr ) ) goto err;
    yFactor = MCGIDI_misc_getUnitConversionFactor( smr, fromUnits[1], toUnits[1] );
    if( !smr_isOk( smr ) ) goto err;
 
 
    ptwXY = ptwXY_create( interpolation, NULL, 2., 1e-3, length, 10, length, data, &status, 0 );
    if( status != nfu_Okay ) {
        smr_setReportError2( smr, smr_unknownID, 1, "ptwXY_create err = %d: %s\n", status, nfu_statusMessage( status ) );
        goto err;
    }
 
    if( ( xFactor != 1. ) || ( yFactor != 1. ) ) {
        if( ( status = ptwXY_scaleOffsetXAndY( ptwXY, xFactor, 0., yFactor, 0. ) ) != nfu_Okay ) {
            smr_setReportError2( smr, smr_unknownID, 1, "ptwXY_scaleOffsetXAndY err = %d: %s\n", status, nfu_statusMessage( status ) );
            goto err;
        }
    }
 
    return( ptwXY );
 
err:
    if( ptwXY != NULL ) ptwXY_free( ptwXY );
    return( NULL );
}

References MCGIDI_misc_getUnitConversionFactor(), nfu_Okay, nfu_statusMessage(), ptwXY_create(), ptwXY_free(), ptwXY_scaleOffsetXAndY(), smr_isOk(), smr_setReportError2, and smr_unknownID.

Referenced by MCGIDI_misc_dataFromElement2ptwXYPointsInUnitsOf(), and MCGIDI_misc_dataFromXYs2ptwXYPointsInUnitsOf().

MCGIDI_misc_dataFromElement2ptwXYPointsInUnitsOf()

ptwXYPoints * MCGIDI_misc_dataFromElement2ptwXYPointsInUnitsOf	(	statusMessageReporting *	smr,
		xDataTOM_element *	linear,
		char const *	toUnits[2]
	)

Definition at line 424 of file MCGIDI_misc.cc.

                                                                                                                                               {
 
    int length;
    double *data;
    xDataTOM_axes *axes = &(linear->xDataInfo.axes);
    char const *fromUnits[2];
    ptwXY_interpolation interpolation;
 
    if( axes->numberOfAxes != 2 ) {
        smr_setReportError2( smr, smr_unknownID, 1, "axes must have 2 axis, it has %d", axes->numberOfAxes );
        return( NULL );
    }
 
    if( MCGIDI_fromTOM_interpolation( smr, linear, 0, &interpolation ) != 0 ) return( NULL );
    fromUnits[0] = axes->axis[0].unit;
    fromUnits[1] = axes->axis[1].unit;
 
    length = xDataTOM_XYs_getDataFromXDataInfo( (xDataTOM_xDataInfo *) &(linear->xDataInfo), &data );
    return( MCGIDI_misc_Data2ptwXYPointsInUnitsOf( smr, interpolation, length, data, fromUnits, toUnits ) );
}

References xDataTOM_xDataInfo_s::axes, xDataTOM_axes_s::axis, MCGIDI_fromTOM_interpolation(), MCGIDI_misc_Data2ptwXYPointsInUnitsOf(), xDataTOM_axes_s::numberOfAxes, smr_setReportError2, smr_unknownID, xDataTOM_axis_s::unit, xDataTOM_element_s::xDataInfo, and xDataTOM_XYs_getDataFromXDataInfo().

Referenced by MCGIDI_energy_parseEvaporationFromTOM(), MCGIDI_energy_parseGeneralEvaporationFromTOM(), MCGIDI_energy_parseSimpleMaxwellianFissionFromTOM(), MCGIDI_energy_parseWattFromTOM(), MCGIDI_energy_parseWeightFromTOM(), MCGIDI_product_parseFromTOM(), MCGIDI_product_parseWeightedReferenceMultiplicityFromTOM(), and MCGIDI_reaction_parseFromTOM().

MCGIDI_misc_dataFromXYs2ptwXYPointsInUnitsOf()

ptwXYPoints * MCGIDI_misc_dataFromXYs2ptwXYPointsInUnitsOf	(	statusMessageReporting *	smr,
		xDataTOM_XYs *	XYs,
		ptwXY_interpolation	interpolation,
		char const *	toUnits[2]
	)

Definition at line 405 of file MCGIDI_misc.cc.

                                                                    {
 
    int length;
    double *data;
    char const *fromUnits[2];
 
    fromUnits[0] = xDataTOM_subAxes_getUnit( smr, &(XYs->subAxes), 0 );
    if( !smr_isOk( smr ) ) return( NULL );
    fromUnits[1] = xDataTOM_subAxes_getUnit( smr, &(XYs->subAxes), 1 );
    if( !smr_isOk( smr ) ) return( NULL );
 
    length = xDataTOM_XYs_getData( XYs, &data );
 
    return( MCGIDI_misc_Data2ptwXYPointsInUnitsOf( smr, interpolation, length, data, fromUnits, toUnits ) );
}

References MCGIDI_misc_Data2ptwXYPointsInUnitsOf(), smr_isOk(), xDataTOM_XYs_s::subAxes, xDataTOM_subAxes_getUnit(), and xDataTOM_XYs_getData().

Referenced by MCGIDI_angular_parseFromTOM(), MCGIDI_angularEnergy_parsePointwiseFromTOM(), MCGIDI_energy_parseMadlandNixFromTOM(), MCGIDI_energyAngular_linear_parseFromTOM(), MCGIDI_fromTOM_pdfsOfXGivenW(), and MCGIDI_LLNL_angularEnergy_parsePointwiseFromTOM().

MCGIDI_misc_getAbsPath()

char * MCGIDI_misc_getAbsPath	(	statusMessageReporting *	smr,
		const char *	fileName
	)

Definition at line 250 of file MCGIDI_misc.cc.

                                                                                  {
/*
*   User must free returned string.
*/
    int n = (int) strlen( fileName ) + 1, nCwd = 0;
    char *absPath, cwd[4 * 1024] = "", *p, *needle;
 
    if( fileName[0] != '/' ) {
        //if( getcwd( cwd, sizeof( cwd ) + 1 ) == NULL ) {
        //TK modified above line for compiler(gcc.4.8) warning message
        if( getcwd( cwd, sizeof( cwd ) ) == NULL ) {
            smr_setReportError2p( smr, smr_unknownID, -1, "hardwired cwd too small" );
            return( NULL );
        }
        nCwd = (int) strlen( cwd );
        n += nCwd + 1;                                  /* cwd + '/'. */
    }
    if( ( absPath = (char *) smr_malloc2( smr, n, 0, "absPath" ) ) == NULL ) return( NULL );
    if( fileName[0] != '/' ) {
        strcpy( absPath, cwd );
        strcat( absPath, "/" );
        strcat( absPath, fileName ); }
    else {
        strcpy( absPath, fileName );
    }
 
    while( 1 ) {                                        /* Remove all ./ from path. */ // Loop checking, 11.06.2015, T. Koi
        if( ( needle = strstr( absPath, "/./" ) ) == NULL ) break;
        p = needle;
        for( needle += 2; *needle; p++, needle++ ) *p = *needle;
        *p = 0;
    }
    while( 1 ) {                                        /* Remove all ../ from path. */ // Loop checking, 11.06.2015, T. Koi
        if( ( needle = strstr( absPath, "/../" ) ) == NULL ) break;
        p = needle - 1;
        while( ( p > absPath ) && ( *p != '/' ) ) p--; // Loop checking, 11.06.2015, T. Koi
        if( *p != '/' ) break;                           /* This should not happen if path is legit, I think, and I do not know what to do so will leave it. */
        if( p == absPath ) break;                       /* Ditto. */
        for( needle += 3; *needle; p++, needle++ ) *p = *needle;
        *p = 0;
    }
    return( absPath );
}

References g4zmq::cwd(), CLHEP::detail::n, smr_malloc2, smr_setReportError2p, and smr_unknownID.

Referenced by MCGIDI_target_read().

MCGIDI_misc_getProductFrame()

enum xDataTOM_frame MCGIDI_misc_getProductFrame	(	statusMessageReporting *	smr,
		xDataTOM_element *	frameElement
	)

Definition at line 315 of file MCGIDI_misc.cc.

                                                                                                               {
 
    char const *frameString;
    enum xDataTOM_frame frame = xDataTOM_frame_invalid;
 
    if( ( frameString = xDataTOM_getAttributesValueInElement( frameElement, MCGIDI_token_productFrame ) ) != NULL ) {
        if( ( frame = xDataTOM_axis_stringToFrame( smr, frameString ) ) == xDataTOM_frame_invalid ) {
            smr_setReportError2( smr, smr_unknownID, 1, "Invalid frame = '%s'", frameString );
        }
    }
    return( frame );
}

References MCGIDI_token_productFrame, smr_setReportError2, smr_unknownID, xDataTOM_axis_stringToFrame(), xDataTOM_frame_invalid, and xDataTOM_getAttributesValueInElement().

Referenced by MCGIDI_angular_parseFromTOM(), MCGIDI_angularEnergy_parsePointwiseFromTOM(), MCGIDI_energy_parseFromTOM(), MCGIDI_energyAngular_linear_parseFromTOM(), MCGIDI_KalbachMann_parseFromTOM(), and MCGIDI_LLNL_angularEnergy_parsePointwiseFromTOM().

MCGIDI_misc_getUnitConversionFactor()

double MCGIDI_misc_getUnitConversionFactor	(	statusMessageReporting *	smr,
		char const *	fromUnit,
		char const *	toUnit
	)

Definition at line 381 of file MCGIDI_misc.cc.

                                                                                                                    {
/*
*   This is a kludge until units are better supported.
*/
    if( strcmp( fromUnit, toUnit ) == 0 ) return( 1.0 );
 
    if( strcmp( fromUnit, "eV" ) == 0 ) {
        if( strcmp( toUnit, "MeV" ) == 0 ) return( 1e-6 ); }
    else if( strcmp( fromUnit, "MeV" ) == 0 ) {
        if( strcmp( toUnit, "eV" ) == 0 ) return( 1e+6 ); }
    else if( strcmp( fromUnit, "1/eV" ) == 0 ) {
        if( strcmp( toUnit, "1/MeV" ) == 0 ) return( 1e+6 ); }
    else if( strcmp( fromUnit, "1/MeV" ) == 0 ) {
        if( strcmp( toUnit, "1/eV" ) == 0 ) return( 1e-6 ); }
    else if( strcmp( fromUnit, "K" ) == 0 ) {
        if( strcmp( toUnit, "MeV/k" ) == 0 ) return( 8.617343183775137e-11 );
    }
 
    smr_setReportError2( smr, smr_unknownID, 1, "Cannot convert unit '%s' to unit '%s'", fromUnit, toUnit );
    return( 1.0 );
}

References smr_setReportError2, and smr_unknownID.

Referenced by MCGIDI_angular_parseFromTOM(), MCGIDI_angularEnergy_parsePointwiseFromTOM(), MCGIDI_energyAngular_linear_parseFromTOM(), MCGIDI_fromTOM_pdfsOfXGivenW(), MCGIDI_KalbachMann_parseFromTOM(), MCGIDI_misc_Data2ptwXYPointsInUnitsOf(), MCGIDI_misc_PQUStringToDoubleInUnitOf(), MCGIDI_product_parsePolynomialMultiplicity(), and MCGIDI_target_heated_read().

MCGIDI_misc_NumberOfZSymbols()

int MCGIDI_misc_NumberOfZSymbols

(

void

)

Definition at line 54 of file MCGIDI_misc.cc.

                                         {
 
    return( sizeof( ZSymbols ) / sizeof( struct ZSymbol ) );
}

References ZSymbols.

Referenced by MCGIDI_misc_symbolToZ(), and MCGIDI_misc_ZToSymbol().

MCGIDI_misc_pointerToAttributeIfAllOk()

char const * MCGIDI_misc_pointerToAttributeIfAllOk	(	statusMessageReporting *	smr,
		xDataXML_element *	element,
		const char *	path,
		int	required,
		xDataTOM_attributionList *	attributes,
		const char *	name,
		const char *	file,
		int	line
	)

Definition at line 184 of file MCGIDI_misc.cc.

                                                                                             {
 
    char const *value;
 
    if( !smr_isOk( smr ) ) return( NULL );
    if( ( value = xDataTOMAL_getAttributesValue( attributes, name ) ) == NULL ) {
        if( required ) {
            if( element != NULL ) {
                MCGIDI_misc_setMessageError_Element( smr, NULL, element, file, line, 1, "element does not have attribute named %s", name ); }
            else {
                smr_setReportError( smr, NULL, file, line, __func__, smr_unknownID, 1, "element does not have attribute named %s for file = %d", name, path );
            }
        }
    }
    return( value );
}

References geant4_check_module_cycles::file, MCGIDI_misc_setMessageError_Element(), G4InuclParticleNames::name(), geant4_check_module_cycles::required, smr_isOk(), smr_setReportError(), smr_unknownID, and xDataTOMAL_getAttributesValue().

MCGIDI_misc_pointerToTOMAttributeIfAllOk()

char const * MCGIDI_misc_pointerToTOMAttributeIfAllOk	(	statusMessageReporting *	smr,
		const char *	path,
		int	required,
		xDataTOM_attributionList *	attributes,
		const char *	name,
		const char *	file,
		int	line
	)

Definition at line 168 of file MCGIDI_misc.cc.

                                                                                             {
 
    char const *value;
 
    if( !smr_isOk( smr ) ) return( NULL );
    if( ( value = xDataTOMAL_getAttributesValue( attributes, name ) ) == NULL ) {
        if( required ) {
            smr_setReportError( smr, NULL, file, line, __func__, smr_unknownID, 1, "element does not have attribute named %s for file = %d", name, path );
        }
    }
    return( value );
}

References geant4_check_module_cycles::file, G4InuclParticleNames::name(), geant4_check_module_cycles::required, smr_isOk(), smr_setReportError(), smr_unknownID, and xDataTOMAL_getAttributesValue().

MCGIDI_misc_PQUStringToDouble()

int MCGIDI_misc_PQUStringToDouble	(	statusMessageReporting *	smr,
		char const *	str,
		char const *	unit,
		double	conversion,
		double *	value
	)

Definition at line 330 of file MCGIDI_misc.cc.

                                                                                                                                      {
/*
*   Currently, white spaces are not allowed after the unit.
*
*   Examples of allowed strings are: "2.39e6 eV", " 2.39e6eV" and " 2.39e6 eV".
*/
    char const *s = str;
    char *e;
 
 
    while( isspace( *s ) ) s++; // Loop checking, 11.06.2015, T. Koi
    *value = strtod( s, &e ) * conversion;
    if( e == s ) {
        smr_setReportError2( smr, smr_unknownID, 1, "no number at start of string = <%s>", str );
        return( 1 );
    }
    while( isspace( *e ) ) e++; // Loop checking, 11.06.2015, T. Koi
    if( strcmp( e, unit ) != 0 ) {
        smr_setReportError2( smr, smr_unknownID, 1, "unit = '%s' not '%s' in '%s'", e, unit, str );
        return( 1 );
    }
    return( 0 );
}

References s, smr_setReportError2, and smr_unknownID.

Referenced by MCGIDI_energy_parseNBodyPhaseSpaceFromTOM(), and MCGIDI_target_heated_parseParticle().

MCGIDI_misc_PQUStringToDoubleInUnitOf()

int MCGIDI_misc_PQUStringToDoubleInUnitOf	(	statusMessageReporting *	smr,
		char const *	str,
		char const *	toUnit,
		double *	value
	)

Definition at line 356 of file MCGIDI_misc.cc.

                                                                                                                             {
/*
*   Currently, white spaces are not allowed after the unit.
*
*   Examples of allowed strings are: "2.39e6 eV", " 2.39e6eV" and " 2.39e6 eV".
*/
    char const *s1 = str;
    char *e1;
    double factor;
 
    while( isspace( *s1 ) ) s1++; // Loop checking, 11.06.2015, T. Koi
    *value = strtod( s1, &e1 );
    if( e1 == s1 ) {
        smr_setReportError2( smr, smr_unknownID, 1, "no number at start of string = <%s>", str );
        return( 1 );
    }
    while( isspace( *e1 ) ) e1++; // Loop checking, 11.06.2015, T. Koi
 
    factor = MCGIDI_misc_getUnitConversionFactor( smr, e1, toUnit );
    *value *= factor;
    return( !smr_isOk( smr ) );
}

References e1, MCGIDI_misc_getUnitConversionFactor(), smr_isOk(), smr_setReportError2, and smr_unknownID.

Referenced by MCGIDI_distribution_parseFromTOM(), MCGIDI_energy_parseEvaporationFromTOM(), MCGIDI_energy_parseMadlandNixFromTOM(), MCGIDI_energy_parseSimpleMaxwellianFissionFromTOM(), MCGIDI_energy_parseWattFromTOM(), MCGIDI_outputChannel_parseFromTOM(), MCGIDI_product_parseFromTOM(), and MCGIDI_target_heated_parseParticleLevel().

MCGIDI_misc_setMessageError_Element()

int MCGIDI_misc_setMessageError_Element	(	statusMessageReporting *	smr,
		void *	userInterface,
		xDataXML_element *	element,
		const char *	file,
		int	line,
		int	code,
		const char *	fmt,
			...
	)

Definition at line 204 of file MCGIDI_misc.cc.

                           {
 
    int status = 0;
    va_list args;
    char *msg;
 
    va_start( args, fmt );
    msg = smr_vallocateFormatMessage( fmt, &args );
    va_end( args );
    if( msg == NULL ) {
        status = 1;
        va_start( args, fmt );
        smr_vsetReportError( smr, userInterface, file, line, __func__, smr_unknownID, code, fmt, &args );
        va_end( args ); }
    else {
        status = smr_setReportError( smr, userInterface, file, line, __func__, smr_unknownID, code, "%s for element %s", msg, element->name );
        smr_freeMemory( (void **) &msg );
    }
    return( status );
}

References args, geant4_check_module_cycles::file, xDataXML_element_s::name, smr_freeMemory(), smr_setReportError(), smr_unknownID, smr_vallocateFormatMessage(), and smr_vsetReportError().

Referenced by MCGIDI_misc_pointerToAttributeIfAllOk(), and MCGIDI_target_read().

MCGIDI_misc_symbolToZ()

int MCGIDI_misc_symbolToZ

(

const char *

Z

)

Definition at line 69 of file MCGIDI_misc.cc.

                                           {
 
    int i, n = MCGIDI_misc_NumberOfZSymbols( );
 
    for( i = 0; i < n; i++ ) {
        if( strcmp( Z, ZSymbols[i].symbol ) == 0 ) return( ZSymbols[i].Z );
    }
    return( -1 );
}

References MCGIDI_misc_NumberOfZSymbols(), CLHEP::detail::n, ZSymbol::symbol, Z, and ZSymbols.

Referenced by MCGIDI_miscNameToZAm().

MCGIDI_misc_updateTransportabilitiesMap()

void MCGIDI_misc_updateTransportabilitiesMap	(	transportabilitiesMap *	transportabilities,
		int	PoPID,
		enum MCGIDI_transportability	transportability
	)

Definition at line 482 of file MCGIDI_misc.cc.

                                                                                                                                                    {
 
    transportabilitiesMap::iterator iter = transportabilities->find( PoPID );
 
    if( iter != transportabilities->end( ) ) {
        switch ( iter->second ) {
        case MCGIDI_transportability_unknown :
            break;
        case MCGIDI_transportability_none :
            switch( transportability ) {
                case MCGIDI_transportability_unknown :
                case MCGIDI_transportability_none :
                    transportability = MCGIDI_transportability_none;
                    break;
                case MCGIDI_transportability_partial :
                    break;
                case MCGIDI_transportability_full :
                    transportability = MCGIDI_transportability_partial;
                    break;
            }
            break;
        case MCGIDI_transportability_partial :
            transportability = MCGIDI_transportability_partial;
            break;
        case MCGIDI_transportability_full :
            switch( transportability ) {
                case MCGIDI_transportability_none :
                case MCGIDI_transportability_partial :
                    transportability = MCGIDI_transportability_partial;
                    break;
                case MCGIDI_transportability_unknown :
                case MCGIDI_transportability_full :
                    break;
            }
            break;
        }
    }
    (*transportabilities)[PoPID] = transportability;
}

References MCGIDI_transportability_full, MCGIDI_transportability_none, MCGIDI_transportability_partial, and MCGIDI_transportability_unknown.

Referenced by MCGIDI_misc_updateTransportabilitiesMap2(), and MCGIDI_target_heated_read().

MCGIDI_misc_updateTransportabilitiesMap2()

void MCGIDI_misc_updateTransportabilitiesMap2	(	transportabilitiesMap *	transportabilities,
		int	PoPID,
		int	transportable
	)

Definition at line 524 of file MCGIDI_misc.cc.

                                                                                                                         {
 
    MCGIDI_misc_updateTransportabilitiesMap( transportabilities, PoPID, ( transportable ? MCGIDI_transportability_full : MCGIDI_transportability_none ) );
}

References MCGIDI_misc_updateTransportabilitiesMap(), MCGIDI_transportability_full, and MCGIDI_transportability_none.

Referenced by MCGIDI_reaction_addReturnProduct().

MCGIDI_misc_ZToSymbol()

const char * MCGIDI_misc_ZToSymbol

(

int

iZ

)

Definition at line 61 of file MCGIDI_misc.cc.

                                            {
 
    if( ( iZ < 0 ) || ( iZ >= MCGIDI_misc_NumberOfZSymbols( ) ) ) return( NULL );
    return( ZSymbols[iZ].symbol );
}

References MCGIDI_misc_NumberOfZSymbols(), ZSymbol::symbol, and ZSymbols.

Referenced by G4GIDI_Misc_Z_A_m_ToName().

MCGIDI_miscNameToZAm()

int MCGIDI_miscNameToZAm	(	statusMessageReporting *	smr,
		const char *	name,
		int *	Z,
		int *	A,
		int *	m,
		int *	level
	)

Definition at line 81 of file MCGIDI_misc.cc.

                                                                                                              {
 
    const char *p;
    char s[1024] = "", *q, *e;   /* Note 1) routine will fail when parts of a particle name can be longer than 1024. */
 
    if( strlen( name ) >= ( sizeof( s ) - 1 ) ) {
        smr_setReportError2( smr, smr_unknownID, 0, "particle name too long: '%s'", name );
        return( 1 );
    }
 
    *Z = *A = *m = *level = 0;
    if( ( !strncmp( "FissionProduct", name, 14 ) ) || !strncmp( "99120", name, 5 ) ) {
        *Z = 99;
        *A = 120;
        return( 0 );
    }
    if( strcmp( "gamma", name ) == 0 ) return( 0 );
    if( strcmp( "n", name ) == 0 ) { *A = 1; return( 0 ); }
 
    for( p = name, q = s; ( *p != 0 ) && !isdigit( *p ) && ( *p != '_' ); p++, q++ ) *q = *p;   /* '_' only for "natural". */
    if( *p == 0 ) {
        smr_setReportError2( smr, smr_unknownID, 0, "unsupported particle name = '%s'", name );
        return( 1 );
    }
    *q = 0;
    if( ( *Z = MCGIDI_misc_symbolToZ( s ) ) < 0 ) {
        smr_setReportError2( smr, smr_unknownID, 1, "Particle %s's symbol = '%s' not found", name, s ); }
    else {                  /* Getting here implies that *p is a digit. */
        if( *p == '_' ) {
            if( strncmp( p, "_natural", 8 ) == 0 ) {
                p += 8;
                if( *p ) *level = MCGIDI_miscNameToZAm_getLevel( smr, name, p ); }
            else {
                smr_setReportError2( smr, smr_unknownID, 0, "expecting 'natural': %s", name );
            } }
        else {
            for( q = s; isdigit( *p ); p++, q++ ) *q = *p;
            *q = 0;
            if( strcmp( s, "natural" ) == 0 ) {
                e = s;
                while( *e ) e++; /* Loop checking, 11.06.2015, T. Koi*/ }
            else {
                *A = (int) strtol( s, &e, 10 );
            }
            if( *e != 0 ) {
                smr_setReportError2( smr, smr_unknownID, 1, "Failed to convert A to integer in particle name %s", name ); }
            else {          /* Getting here implies that *p == '_' or 0. */
                if( *p ) *level = MCGIDI_miscNameToZAm_getLevel( smr, name, p );
            }
        }
    }
 
    return( !smr_isOk( smr ) );
}

References A, m, MCGIDI_misc_symbolToZ(), MCGIDI_miscNameToZAm_getLevel(), G4InuclParticleNames::name(), s, smr_isOk(), smr_setReportError2, smr_unknownID, and Z.

Referenced by G4GIDI_Misc_channelCompound(), MCGIDI_particle_getInternalID(), and MCGIDI_POP_new().

MCGIDI_miscNameToZAm_getLevel()

static int MCGIDI_miscNameToZAm_getLevel ( statusMessageReporting *  smr, const char *  name, const char *  p  )

static

Definition at line 138 of file MCGIDI_misc.cc.

                                                                                                         {
 
    int level = 0;
    char *e;
 
    if( *p == '_' ) {
        p++;
        switch( *p ) {
        case 'e' :
            p++;
            level = (int) strtol( p, &e, 10 );
            if( *e != 0 ) smr_setReportError2( smr, smr_unknownID, 1, "Failed to convert level to integer in particle name %s", name );
            break;
        case 'c' :
            level = MCGIDI_particleLevel_continuum;
            break;
        case 's' :
            level = MCGIDI_particleLevel_sum;
            break;
        default :
            smr_setReportError2( smr, smr_unknownID, 0, "invalid 'natural': %s", name );
        } }
    else {
        smr_setReportError2( smr, smr_unknownID, 0, "invalid level specifier: %s", name );
    }
    return( level );
}

References MCGIDI_particleLevel_continuum, MCGIDI_particleLevel_sum, G4InuclParticleNames::name(), smr_setReportError2, and smr_unknownID.

Referenced by MCGIDI_miscNameToZAm().

Variable Documentation

ZSymbols

struct ZSymbol ZSymbols[]

static

Initial value:

= {    {   0, "n"  },  {   1, "H"  },  {   2, "He" },  {   3, "Li" },  {   4, "Be" },  {   5, "B"  },  {   6, "C"  },
        {   7, "N"  },  {   8, "O"  },  {   9, "F"  },  {  10, "Ne" },  {  11, "Na" },  {  12, "Mg" },  {  13, "Al" },  {  14, "Si" },  {  15, "P"  },
        {  16, "S"  },  {  17, "Cl" },  {  18, "Ar" },  {  19, "K"  },  {  20, "Ca" },  {  21, "Sc" },  {  22, "Ti" },  {  23, "V"  },  {  24, "Cr" },
        {  25, "Mn" },  {  26, "Fe" },  {  27, "Co" },  {  28, "Ni" },  {  29, "Cu" },  {  30, "Zn" },  {  31, "Ga" },  {  32, "Ge" },  {  33, "As" },
        {  34, "Se" },  {  35, "Br" },  {  36, "Kr" },  {  37, "Rb" },  {  38, "Sr" },  {  39, "Y"  },  {  40, "Zr" },  {  41, "Nb" },  {  42, "Mo" },
        {  43, "Tc" },  {  44, "Ru" },  {  45, "Rh" },  {  46, "Pd" },  {  47, "Ag" },  {  48, "Cd" },  {  49, "In" },  {  50, "Sn" },  {  51, "Sb" },
        {  52, "Te" },  {  53, "I"  },  {  54, "Xe" },  {  55, "Cs" },  {  56, "Ba" },  {  57, "La" },  {  58, "Ce" },  {  59, "Pr" },  {  60, "Nd" },
        {  61, "Pm" },  {  62, "Sm" },  {  63, "Eu" },  {  64, "Gd" },  {  65, "Tb" },  {  66, "Dy" },  {  67, "Ho" },  {  68, "Er" },  {  69, "Tm" },
        {  70, "Yb" },  {  71, "Lu" },  {  72, "Hf" },  {  73, "Ta" },  {  74, "W"  },  {  75, "Re" },  {  76, "Os" },  {  77, "Ir" },  {  78, "Pt" },
        {  79, "Au" },  {  80, "Hg" },  {  81, "Tl" },  {  82, "Pb" },  {  83, "Bi" },  {  84, "Po" },  {  85, "At" },  {  86, "Rn" },  {  87, "Fr" },
        {  88, "Ra" },  {  89, "Ac" },  {  90, "Th" },  {  91, "Pa" },  {  92, "U"  },  {  93, "Np" },  {  94, "Pu" },  {  95, "Am" },  {  96, "Cm" },
        {  97, "Bk" },  {  98, "Cf" },  {  99, "Es" },  { 100, "Fm" },  { 101, "Md" },  { 102, "No" },  { 103, "Lr" },  { 104, "Rf" },  { 105, "Db" },
        { 106, "Sg" }, { 107, "Bh" },  { 108, "Hs" },  { 109, "Mt" } }

Definition at line 34 of file MCGIDI_misc.cc.

Referenced by MCGIDI_misc_NumberOfZSymbols(), MCGIDI_misc_symbolToZ(), and MCGIDI_misc_ZToSymbol().