tpia_target.cc File Reference

#include <string.h>
#include <tpia_target.h>
#include <tpia_misc.h>

Go to the source code of this file.

Functions
tpia_target *	tpia_target_create (statusMessageReporting *smr)
int	tpia_target_initialize (statusMessageReporting *smr, tpia_target *target)
tpia_target *	tpia_target_createRead (statusMessageReporting *smr, const char *fileName)
int	tpia_target_readFromMap (statusMessageReporting *smr, tpia_target *target, tpia_map *map, const char *evaluation, const char *projectileName, const char *targetName)
tpia_target *	tpia_target_createReadFromMap (statusMessageReporting *smr, tpia_map *map, const char *evaluation, const char *projectileName, const char *targetName)
tpia_target *	tpia_target_free (statusMessageReporting *smr, tpia_target *target)
int	tpia_target_release (statusMessageReporting *smr, tpia_target *target)
int	tpia_target_read (statusMessageReporting *smr, tpia_target *target, const char *fileName)
char *	tpia_target_getAttributesValue (statusMessageReporting *, tpia_target *target, char const *name)
int	tpia_target_getTemperatures (statusMessageReporting *, tpia_target *target, double *temperatures)
int	tpia_target_readHeatedTarget (statusMessageReporting *smr, tpia_target *target, int index, int checkElememtsForAccess)
tpia_target_heated *	tpia_target_getHeatedTargetAtIndex_ReadIfNeeded (statusMessageReporting *smr, tpia_target *target, int index)
tpia_target_heated *	tpia_target_getHeatedTargetAtTIndex (statusMessageReporting *smr, tpia_target *target, int index)
int	tpia_target_numberOfChannels (statusMessageReporting *smr, tpia_target *target)
int	tpia_target_numberOfProductionChannels (statusMessageReporting *smr, tpia_target *target)
xData_Int	tpia_target_getEnergyGridAtTIndex (statusMessageReporting *smr, tpia_target *target, int index, double **energyGrid)
tpia_1dData *	tpia_target_getTotalCrossSectionAtTIndex (statusMessageReporting *smr, tpia_target *target, int index, int crossSectionType)
double	tpia_target_getTotalCrossSectionAtTAndE (statusMessageReporting *smr, tpia_target *target, double T, xData_Int iEg, double e_in, int crossSectionType)
double	tpia_target_getIndexChannelCrossSectionAtE (statusMessageReporting *smr, tpia_target *target, int index, double T, xData_Int iEg, double e_in, int crossSectionType)
int	tpia_target_sampleIndexChannelProductsAtE (statusMessageReporting *smr, tpia_target *target, int index, double, tpia_decaySamplingInfo *decaySamplingInfo, int nProductData, tpia_productOutgoingData *productData)

Function Documentation

tpia_target_create()

tpia_target * tpia_target_create

(

statusMessageReporting *

smr

)

Definition at line 49 of file tpia_target.cc.

                                                               {
 
    tpia_target *target;
 
    //if( ( target = xData_malloc2( smr, sizeof( tpia_target ), 0, "target" ) ) == NULL ) return( NULL );
    if( ( target = (tpia_target*) xData_malloc2( smr, sizeof( tpia_target ), 0, "target" ) ) == NULL ) return( NULL );
    if( tpia_target_initialize( smr, target ) ) target = tpia_target_free( smr, target );
    return( target );
}

Referenced by tpia_target_createRead().

tpia_target_createRead()

tpia_target * tpia_target_createRead	(	statusMessageReporting *	smr,
		const char *	fileName
	)

Definition at line 70 of file tpia_target.cc.

                                                                                         {
 
    tpia_target *target;
 
    if( ( target = tpia_target_create( smr ) ) == NULL ) return( NULL );
    //if( tpia_target_read( smr, target, fileName ) != 0 ) target = xData_free( smr, target );
    if( tpia_target_read( smr, target, fileName ) != 0 ) target = (tpia_target*) xData_free( smr, target );
    return( target );
}

Referenced by G4GIDI_target::init(), and tpia_target_createReadFromMap().

tpia_target_createReadFromMap()

tpia_target * tpia_target_createReadFromMap	(	statusMessageReporting *	smr,
		tpia_map *	map,
		const char *	evaluation,
		const char *	projectileName,
		const char *	targetName
	)

Definition at line 93 of file tpia_target.cc.

                                 {
 
    char *targetPath;
    tpia_target *target;
 
    targetPath = tpia_map_findTarget( smr, map, evaluation, projectileName, targetName );
    if( targetPath == NULL ) return( NULL );
    target = tpia_target_createRead( smr, targetPath );
    xData_free( smr, targetPath );
    return( target );
}

Referenced by tpia_depot_addTargetFromMap().

tpia_target_free()

tpia_target * tpia_target_free	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 108 of file tpia_target.cc.

                                                                                  {
 
    tpia_target_release( smr, target );
    xData_free( smr, target );
    return( NULL );
}

Referenced by tpia_depot_release(), tpia_target_create(), and G4GIDI_target::~G4GIDI_target().

tpia_target_getAttributesValue()

char * tpia_target_getAttributesValue	(	statusMessageReporting *	smr,
		tpia_target *	target,
		char const *	name
	)

Definition at line 248 of file tpia_target.cc.

                                                                                                        {
 
    return( xData_getAttributesValue( &(target->attributes), name ) );
}

tpia_target_getEnergyGridAtTIndex()

xData_Int tpia_target_getEnergyGridAtTIndex	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		double **	energyGrid
	)

Definition at line 331 of file tpia_target.cc.

                                                                                                                                {
 
    tpia_target_heated *heated = tpia_target_getHeatedTargetAtTIndex( smr, target, index );
 
    if( !smr_isOk( smr ) ) return( -1 );
    return( tpia_target_heated_getEnergyGrid( smr, heated, energyGrid ) );
}

Referenced by G4GIDI_target::getEnergyGridAtTIndex().

tpia_target_getHeatedTargetAtIndex_ReadIfNeeded()

tpia_target_heated * tpia_target_getHeatedTargetAtIndex_ReadIfNeeded	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index
	)

Definition at line 290 of file tpia_target.cc.

                                                                                                                                   {
 
    if( ( index < 0 ) || ( index >= target->nHeatedTargets ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "temperature index = %d out of range (0 <= index < %d", index, target->nHeatedTargets );
        return( NULL );
    }
    if( target->heatedTargets[index].heatedTarget == NULL ) tpia_target_readHeatedTarget( smr, target, index, 0 );
    return( target->heatedTargets[index].heatedTarget );
}

tpia_target_getHeatedTargetAtTIndex()

tpia_target_heated * tpia_target_getHeatedTargetAtTIndex	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index
	)

Definition at line 302 of file tpia_target.cc.

                                                                                                                       {
 
    if( ( index < 0 ) || ( index >= target->nHeatedTargets ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "temperature index = %d out of range (0 <= index < %d", index, target->nHeatedTargets );
        return( NULL );
    }
    if( target->heatedTargets[index].heatedTarget == NULL ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "temperature index = %d not read in", index );
        return( NULL );
    }
    return( target->heatedTargets[index].heatedTarget );
}

Referenced by tpia_target_getEnergyGridAtTIndex(), and tpia_target_getTotalCrossSectionAtTIndex().

tpia_target_getIndexChannelCrossSectionAtE()

double tpia_target_getIndexChannelCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		double	T,
		xData_Int	iEg,
		double	e_in,
		int	crossSectionType
	)

Definition at line 380 of file tpia_target.cc.

                           {
 
    int i;
    double xsec = 0., xsec1, xsec2;
 
    for( i = 0; i < target->nReadHeatedTargets; i++ ) if( target->readHeatedTargets[i]->temperature > T ) break;
    if( i == 0 ) {
        xsec = tpia_target_heated_getIndexChannelCrossSectionAtE( smr, target->readHeatedTargets[0]->heatedTarget, index, iEg, e_in, crossSectionType ); }
    else if( i == target->nReadHeatedTargets ) {
        xsec = tpia_target_heated_getIndexChannelCrossSectionAtE( smr, target->readHeatedTargets[i-1]->heatedTarget, index, iEg, e_in, crossSectionType ); }
    else {
        xsec1 = tpia_target_heated_getIndexChannelCrossSectionAtE(smr, target->readHeatedTargets[i-1]->heatedTarget, index, iEg, e_in, crossSectionType);
        xsec2 = tpia_target_heated_getIndexChannelCrossSectionAtE(smr, target->readHeatedTargets[i  ]->heatedTarget, index, iEg, e_in, crossSectionType);
        xsec = ( ( target->readHeatedTargets[i]->temperature - T ) * xsec1 + ( T - target->readHeatedTargets[i-1]->temperature ) * xsec2 ) / 
            ( target->readHeatedTargets[i]->temperature - target->readHeatedTargets[i-1]->temperature );
    }
 
    return( xsec );
}

Referenced by G4GIDI_target::sampleChannelCrossSectionAtE(), and G4GIDI_target::sumChannelCrossSectionAtE().

tpia_target_getTemperatures()

int tpia_target_getTemperatures	(	statusMessageReporting *	smr,
		tpia_target *	target,
		double *	temperatures
	)

Definition at line 256 of file tpia_target.cc.

                                                                                                       {
 
    int i;
 
    if( temperatures != NULL ) for( i = 0; i < target->nHeatedTargets; i++ ) temperatures[i] = target->heatedTargets[i].temperature;
    return( target->nHeatedTargets );
}

Referenced by G4GIDI_target::getTemperatures().

tpia_target_getTotalCrossSectionAtTAndE()

double tpia_target_getTotalCrossSectionAtTAndE	(	statusMessageReporting *	smr,
		tpia_target *	target,
		double	T,
		xData_Int	iEg,
		double	e_in,
		int	crossSectionType
	)

Definition at line 357 of file tpia_target.cc.

                           {
 
    int i;
    double xsec = 0., xsec1, xsec2;
 
    for( i = 0; i < target->nReadHeatedTargets; i++ ) if( target->readHeatedTargets[i]->temperature > T ) break;
    if( i == 0 ) {
        xsec = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[0]->heatedTarget, iEg, e_in, crossSectionType ); }
    else if( i == target->nReadHeatedTargets ) {
        xsec = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[i-1]->heatedTarget, iEg, e_in, crossSectionType ); }
    else {
        xsec1 = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[i-1]->heatedTarget, iEg, e_in, crossSectionType );
        xsec2 = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[i  ]->heatedTarget, iEg, e_in, crossSectionType );
        xsec = ( ( target->readHeatedTargets[i]->temperature - T ) * xsec1 + ( T - target->readHeatedTargets[i-1]->temperature ) * xsec2 ) / 
            ( target->readHeatedTargets[i]->temperature - target->readHeatedTargets[i-1]->temperature );
    }
 
    return( xsec );
}

Referenced by G4GIDI_target::getTotalCrossSectionAtE().

tpia_target_getTotalCrossSectionAtTIndex()

tpia_1dData * tpia_target_getTotalCrossSectionAtTIndex	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		int	crossSectionType
	)

Definition at line 341 of file tpia_target.cc.

                                                                                                                                           {
 
    tpia_target_heated *heated = tpia_target_getHeatedTargetAtTIndex( smr, target, index );
 
    if( !smr_isOk( smr ) ) return( NULL );
    if( crossSectionType == tpia_crossSectionType_grouped ) {
        return( &(heated->totalCrossSectionGrouped) ); }
    else if( crossSectionType == tpia_crossSectionType_pointwise ) {
        return( &(heated->totalCrossSectionPointwise) );
    }
    smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Invalue crossSectionType = %d", crossSectionType );
    return( NULL );
}

tpia_target_initialize()

int tpia_target_initialize	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 61 of file tpia_target.cc.

                                                                               {
 
    memset( target, 0, sizeof( tpia_target ) );
    tpia_samplingMethods_initialize( smr, &(target->samplingMethods) );
    return( 0 );
}

Referenced by tpia_target_create(), tpia_target_read(), and tpia_target_release().

tpia_target_numberOfChannels()

int tpia_target_numberOfChannels	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 317 of file tpia_target.cc.

                                                                                     {
 
    return( tpia_target_heated_numberOfChannels( smr, target->baseHeatedTarget ) );
}

Referenced by G4GIDI_target::getChannelIDs(), G4GIDI_target::getNumberOfChannels(), and G4GIDI_target::init().

tpia_target_numberOfProductionChannels()

int tpia_target_numberOfProductionChannels	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 324 of file tpia_target.cc.

                                                                                               {
 
    return( tpia_target_heated_numberOfProductionChannels( smr, target->baseHeatedTarget ) );
}

Referenced by G4GIDI_target::getNumberOfProductionChannels(), and G4GIDI_target::getProductionChannelIDs().

tpia_target_read()

int tpia_target_read	(	statusMessageReporting *	smr,
		tpia_target *	target,
		const char *	fileName
	)

Definition at line 143 of file tpia_target.cc.

                                                                                               {
/*
*   If a target has already been read into this target, user must have called tpia_target_release before calling this routine.
*   Otherwise, there will be memory leaks.
*/
    xData_document *doc;
    xData_element *element, *child;
    int i, iHeated, nHeated = 0, status = 1;
    double temperature;
    //fix for gcc4.6 warings 110602
    //char *pReturnValue, *name;
    char *name;
    char const *contents;
 
    tpia_target_initialize( smr, target );
    if( ( target->path = xDataMisc_allocateCopyString2( smr, fileName, "path" ) ) == NULL ) return( status );
    if( ( target->absPath = xDataMisc_getAbsPath( smr, fileName ) ) == NULL ) return( _tpia_target_releaseAndReturnOne( smr, target ) );
    if( ( doc = xData_parseReadFile( smr, fileName, NULL, NULL ) ) == NULL ) return( _tpia_target_releaseAndReturnOne( smr, target ) );
    element = xData_getDocumentsElement( doc );
    if( strcmp( element->name, "xTarget" ) != 0 ) {
        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "input file's top element must be xTarget and not %s", element->name ); }
    else {
        //pReturnValue = ( xData_copyAttributionList( smr, &(target->attributes), &(element->attributes) ) ) ? NULL : target->path;
        //fix for gcc4.6 warings 110602
        xData_copyAttributionList( smr, &(target->attributes),&(element->attributes) );
        name = tpia_misc_pointerToAttributeIfAllOk2( smr, element, 1, &(target->attributes), "projectile" );
        if( smr_isOk( smr ) ) target->projectileID = tpia_particle_getInternalID( smr, name );
        if( smr_isOk( smr ) && ( name = tpia_misc_pointerToAttributeIfAllOk2( smr, element, 1, &(target->attributes), "target" ) ) != NULL ) {
            if( smr_isOk( smr ) && ( target->targetID = tpia_particle_getInternalID( smr, name ) ) != NULL ) {
                status = 0;
                for( nHeated = 0, child = xData_getFirstElement( element ); child != NULL; nHeated++, child = xData_getNextElement( child ) ) {
                    if( strcmp( child->name, "target" ) != 0 ) {
                        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "element can only have target sub-elements%s", 
                                element->name );
                        status = 1;
                        break;
                    }
                }
                if( status == 0 ) {
                    //if( ( target->heatedTargets = xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info ), 1, "heatedTargets" ) ) == NULL ) {
                    if( ( target->heatedTargets = (tpia_target_heated_info*) xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info ), 1, "heatedTargets" ) ) == NULL ) {
                        status = 1; }
                    else {
                        //if( ( target->readHeatedTargets = xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info * ), 1, "heatedTargets" ) ) == NULL ) 
                        if( ( target->readHeatedTargets = (tpia_target_heated_info**) xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info * ), 1, "heatedTargets" ) ) == NULL ) 
                                status = 1;
                    }
                    for( nHeated = 0, child = xData_getFirstElement( element ); ( status == 0 ) && ( child != NULL ); 
                            nHeated++, child = xData_getNextElement( child ) ) {
                        if( ( i = xData_convertAttributeToDouble( smr, child, "temperature", &temperature ) ) != 0 ) {
                            if( i > 0 ) smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "target does not have a temperature attribute" );
                            status = 1;
                            break;
                        }
                        for( iHeated = 0; iHeated < nHeated; iHeated++ ) if( target->heatedTargets[iHeated].temperature > temperature ) break;
                        if( iHeated < nHeated ) for( i = nHeated; i >= iHeated; i-- ) target->heatedTargets[i+1] = target->heatedTargets[i];
                        target->heatedTargets[iHeated].temperature = temperature;
                        target->heatedTargets[iHeated].path = NULL;
                        target->heatedTargets[iHeated].contents = NULL;
                        target->heatedTargets[iHeated].heatedTarget = NULL;
                        if( ( contents = xData_getAttributesValueInElement( child, "contents" ) ) != NULL ) {
                            if( ( target->heatedTargets[iHeated].contents = xDataMisc_allocateCopyString2( smr, contents, "contents" ) ) == NULL ) {
                                status = 1;
                                break;
                            }
                        }
                        if( ( contents = xData_getAttributesValueInElement( child, "file" ) ) == NULL ) {
                            status = 1;
                            break;
                        }
                        //if((target->heatedTargets[iHeated].path = xData_malloc2(smr, strlen( target->absPath ) + strlen( contents ) + 2, 0, "path")) == NULL) {
                        if((target->heatedTargets[iHeated].path = (char*) xData_malloc2(smr, strlen( target->absPath ) + strlen( contents ) + 2, 0, "path")) == NULL) {
                            status = 1;
                            break;
                        }
                        strcpy( target->heatedTargets[iHeated].path, target->absPath );
                        *strrchr( target->heatedTargets[iHeated].path, '/' ) = 0;
                        strcat( target->heatedTargets[iHeated].path, "/" );
                        strcat( target->heatedTargets[iHeated].path, contents );
                        target->nHeatedTargets++;
                    }
                }
            }
        }
    }
    xData_parseFree( smr, doc );
    if( status == 0 ) {
        for( i = 0; i < nHeated; i++ ) target->heatedTargets[i].ordinal = i;
        for( i = 0; i < nHeated; i++ ) if( target->heatedTargets[i].contents == NULL ) break;
        if( i == nHeated ) i = 0;                                           /* All heated targets are crossSection only. */
        if( tpia_target_readHeatedTarget( smr, target, i, 0 ) == 0 ) {
            target->baseHeatedTarget = target->heatedTargets[i].heatedTarget; }
        else {
            tpia_target_release( NULL, target );
            status = 1;
        } }
    else {
        tpia_target_release( smr, target );
    }
    return( status );
}

Referenced by tpia_target_createRead(), and tpia_target_readFromMap().

tpia_target_readFromMap()

int tpia_target_readFromMap	(	statusMessageReporting *	smr,
		tpia_target *	target,
		tpia_map *	map,
		const char *	evaluation,
		const char *	projectileName,
		const char *	targetName
	)

Definition at line 82 of file tpia_target.cc.

                                 {
 
    char *targetPath;
 
    if( ( targetPath = tpia_map_findTarget( smr, map, evaluation, projectileName, targetName ) ) == NULL ) return( 1 );
    return( tpia_target_read( smr, target, targetPath ) );
}

tpia_target_readHeatedTarget()

int tpia_target_readHeatedTarget	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		int	checkElememtsForAccess
	)

Definition at line 266 of file tpia_target.cc.

                                                                                                                            {
 
    int i;
 
    if( ( index < 0 ) || ( index >= target->nHeatedTargets ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "temperature index = %d out of range (0 <= index < %d", index, target->nHeatedTargets );
        return( -1 );
    }
    if( target->heatedTargets[index].heatedTarget != NULL ) return( 1 );
    target->heatedTargets[index].heatedTarget = tpia_target_heated_createRead( smr, target->heatedTargets[index].path, checkElememtsForAccess );
    if( target->heatedTargets[index].heatedTarget != NULL ) {
        target->heatedTargets[index].heatedTarget->ordinal = target->heatedTargets[index].ordinal;
        for( i = target->nReadHeatedTargets; i > 0; i-- ) {
            if( target->readHeatedTargets[i-1]->temperature < target->heatedTargets[index].temperature ) break;
            target->readHeatedTargets[i] = target->readHeatedTargets[i-1];
        }
        target->readHeatedTargets[i] = &(target->heatedTargets[i]);
        target->nReadHeatedTargets++;
    }
    return( ( target->heatedTargets[index].heatedTarget == NULL ? -1 : 0 ) );
}

Referenced by G4GIDI_target::readTemperature(), tpia_target_getHeatedTargetAtIndex_ReadIfNeeded(), and tpia_target_read().

tpia_target_release()

int tpia_target_release	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 117 of file tpia_target.cc.

                                                                            {
 
    int i;
 
    //target->path = xData_free( smr, target->path );
    target->path = (char*) xData_free( smr, target->path );
    //target->absPath = xData_free( smr, target->absPath );
    target->absPath = (char*) xData_free( smr, target->absPath );
    xData_releaseAttributionList( smr, &(target->attributes) );
    for( i = 0; i < target->nHeatedTargets; i++ ) {
        //target->heatedTargets[i].path = xData_free( smr, target->heatedTargets[i].path );
        target->heatedTargets[i].path = (char*) xData_free( smr, target->heatedTargets[i].path );
        //target->heatedTargets[i].contents = xData_free( smr, target->heatedTargets[i].contents );
        target->heatedTargets[i].contents = (char*) xData_free( smr, target->heatedTargets[i].contents );
            if( target->heatedTargets[i].heatedTarget != NULL ) tpia_target_heated_free( smr, target->heatedTargets[i].heatedTarget );
    }
    //target->heatedTargets = xData_free( smr, target->heatedTargets );
    target->heatedTargets = (tpia_target_heated_info*) xData_free( smr, target->heatedTargets );
    //target->readHeatedTargets = xData_free( smr, target->readHeatedTargets );
    target->readHeatedTargets = (tpia_target_heated_info**) xData_free( smr, target->readHeatedTargets );
    tpia_target_initialize( smr, target );
    return( 0 );
}

Referenced by tpia_target_free(), and tpia_target_read().

tpia_target_sampleIndexChannelProductsAtE()

int tpia_target_sampleIndexChannelProductsAtE	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		double	T,
		tpia_decaySamplingInfo *	decaySamplingInfo,
		int	nProductData,
		tpia_productOutgoingData *	productData
	)

Definition at line 404 of file tpia_target.cc.

                                                                                                             {
 
    return( tpia_target_heated_sampleIndexChannelProductsAtE( smr, target->baseHeatedTarget, index, decaySamplingInfo,
        nProductData, productData ) );
}