tpia_target_heated.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_heated *	tpia_target_heated_create (statusMessageReporting *smr)
int	tpia_target_heated_initialize (statusMessageReporting *, tpia_target_heated *target)
tpia_target_heated *	tpia_target_heated_createRead (statusMessageReporting *smr, const char *fileName, int checkElememtsForAccess)
tpia_target_heated *	tpia_target_heated_free (statusMessageReporting *smr, tpia_target_heated *target)
int	tpia_target_heated_release (statusMessageReporting *smr, tpia_target_heated *target)
int	tpia_target_heated_read (statusMessageReporting *smr, tpia_target_heated *target, const char *fileName, int checkElememtsForAccess)
int	tpia_target_heated_numberOfChannels (statusMessageReporting *, tpia_target_heated *target)
int	tpia_target_heated_numberOfProductionChannels (statusMessageReporting *, tpia_target_heated *target)
tpia_channel *	tpia_target_heated_getChannelAtIndex (tpia_target_heated *target, int index)
tpia_channel *	tpia_target_heated_getChannelAtIndex_smr (statusMessageReporting *smr, tpia_target_heated *target, int index)
tpia_channel *	tpia_target_heated_getProductionChannelAtIndex (tpia_target_heated *target, int index)
xData_Int	tpia_target_heated_getEnergyGrid (statusMessageReporting *, tpia_target_heated *target, double **energyGrid)
xData_Int	tpia_target_heated_getEIndex (tpia_target_heated *target, double e_in)
double	tpia_target_heated_getTotalCrossSectionAtE (statusMessageReporting *smr, tpia_target_heated *target, xData_Int, double e_in, int crossSectionType)
double	tpia_target_heated_getIndexChannelCrossSectionAtE (statusMessageReporting *smr, tpia_target_heated *target, int index, xData_Int iEg, double e_in, int crossSectionType)
int	tpia_target_heated_sampleIndexChannelProductsAtE (statusMessageReporting *smr, tpia_target_heated *target, int index, tpia_decaySamplingInfo *decaySamplingInfo, int nProductData, tpia_productOutgoingData *productDatas)

Function Documentation

tpia_target_heated* tpia_target_heated_create

(

statusMessageReporting *

smr

)

Definition at line 52 of file tpia_target_heated.cc.

References tpia_target_heated_initialize(), xData_free(), and xData_malloc2.

Referenced by tpia_target_heated_createRead().

                                                                             {

    tpia_target_heated *target;

    //if( ( target = xData_malloc2( smr, sizeof( tpia_target_heated ), 0, "target" ) ) == NULL ) return( NULL );
    if( ( target = (tpia_target_heated*) xData_malloc2( smr, sizeof( tpia_target_heated ), 0, "target" ) ) == NULL ) return( NULL );
    //if( tpia_target_heated_initialize( smr, target ) ) target = xData_free( smr, target );
    if( tpia_target_heated_initialize( smr, target ) ) target = (tpia_target_heated*) xData_free( smr, target );
    return( target );
}

tpia_target_heated* tpia_target_heated_createRead	(	statusMessageReporting *	smr,
		const char *	fileName,
		int	checkElememtsForAccess
	)

Definition at line 74 of file tpia_target_heated.cc.

References tpia_target_heated_create(), tpia_target_heated_read(), and xData_free().

Referenced by tpia_target_readHeatedTarget().

                                                                                                                                   {

    tpia_target_heated *target;

    if( ( target = tpia_target_heated_create( smr ) ) == NULL ) return( NULL );
    //if( tpia_target_heated_read( smr, target, fileName, checkElememtsForAccess ) != 0 ) target = xData_free( smr, target );
    if( tpia_target_heated_read( smr, target, fileName, checkElememtsForAccess ) != 0 ) target = (tpia_target_heated*) xData_free( smr, target );
    return( target );
}

tpia_target_heated* tpia_target_heated_free	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 86 of file tpia_target_heated.cc.

References tpia_target_heated_release(), and xData_free().

Referenced by tpia_target_release().

                                                                                                       {

    tpia_target_heated_release( smr, target );
    xData_free( smr, target );
    return( NULL );
}

tpia_channel* tpia_target_heated_getChannelAtIndex	(	tpia_target_heated *	target,
		int	index
	)

Definition at line 263 of file tpia_target_heated.cc.

References tpia_target_heated_s::channels.

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

                                                                                            {

    tpia_channel *channel = NULL;

    if( ( index >= 0 ) && ( index < target->nChannels ) ) channel = target->channels[index];
    return( channel );
}

tpia_channel* tpia_target_heated_getChannelAtIndex_smr	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		int	index
	)

Definition at line 273 of file tpia_target_heated.cc.

References tpia_particle_s::name, tpia_target_heated_s::projectileID, smr_setMessageError(), tpia_target_heated_s::targetID, and tpia_target_heated_getChannelAtIndex().

Referenced by G4GIDI_target::getElasticFinalState(), tpia_target_heated_getIndexChannelCrossSectionAtE(), and tpia_target_heated_sampleIndexChannelProductsAtE().

                                                                                                                             {

    tpia_channel *channel = tpia_target_heated_getChannelAtIndex( target, index );

    if( channel == NULL ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad channel index = %d for %s + %s", index, 
            target->projectileID->name, target->targetID->name ); 
    }
    return( channel );
}

xData_Int tpia_target_heated_getEIndex	(	tpia_target_heated *	target,
		double	e_in
	)

Definition at line 305 of file tpia_target_heated.cc.

References tpia_target_heated_s::energyGrid, tpia_target_heated_s::energyGridLength, and tpia_misc_binarySearch().

Referenced by tpia_channel_getCrossSectionAtE(), and tpia_target_heated_getTotalCrossSectionAtE().

                                                                                  {

    return( tpia_misc_binarySearch( target->energyGridLength, target->energyGrid, e_in ) );
}

xData_Int tpia_target_heated_getEnergyGrid	(	statusMessageReporting *	,
		tpia_target_heated *	target,
		double **	energyGrid
	)

Definition at line 297 of file tpia_target_heated.cc.

References tpia_target_heated_s::energyGrid, and tpia_target_heated_s::energyGridLength.

Referenced by tpia_target_getEnergyGridAtTIndex().

                                                                                                                        {

    if( energyGrid != NULL ) *energyGrid = target->energyGrid;
    return( target->energyGridLength );
}

double tpia_target_heated_getIndexChannelCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		int	index,
		xData_Int	iEg,
		double	e_in,
		int	crossSectionType
	)

Definition at line 329 of file tpia_target_heated.cc.

References tpia_channel_getCrossSectionAtE(), and tpia_target_heated_getChannelAtIndex_smr().

Referenced by tpia_target_getIndexChannelCrossSectionAtE().

                               {

    double xsec = 0.;
    tpia_channel *channel = tpia_target_heated_getChannelAtIndex_smr( smr, target, index );

    if( channel != NULL ) xsec = tpia_channel_getCrossSectionAtE( smr, channel, iEg, e_in, crossSectionType );
    return( xsec );
}

tpia_channel* tpia_target_heated_getProductionChannelAtIndex	(	tpia_target_heated *	target,
		int	index
	)

Definition at line 286 of file tpia_target_heated.cc.

References tpia_target_heated_s::productionChannels.

                                                                                                      {

    tpia_channel *channel = NULL;

    if( ( index >= 0 ) && ( index < target->nProductionChannels ) ) channel = target->productionChannels[index];
    return( channel );
}

double tpia_target_heated_getTotalCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		xData_Int	,
		double	e_in,
		int	crossSectionType
	)

Definition at line 313 of file tpia_target_heated.cc.

References tpia_target_heated_s::energyGrid, tpia_target_heated_s::totalCrossSectionPointwise, tpia_crossSectionType_grouped, tpia_crossSectionType_pointwise, tpia_misc_getPointwiseCrossSectionAtE(), and tpia_target_heated_getEIndex().

Referenced by tpia_target_getTotalCrossSectionAtTAndE().

                               {

    double xsec = 0.;

    if( crossSectionType == tpia_crossSectionType_grouped ) {
        xsec = 0; }
    else if( crossSectionType == tpia_crossSectionType_pointwise ) {
        xsec = tpia_misc_getPointwiseCrossSectionAtE( smr, &(target->totalCrossSectionPointwise), target->energyGrid, 
            tpia_target_heated_getEIndex( target, e_in ), e_in );
    }
    return( xsec );
}

int tpia_target_heated_initialize	(	statusMessageReporting *	,
		tpia_target_heated *	target
	)

Definition at line 66 of file tpia_target_heated.cc.

Referenced by tpia_target_heated_create(), and tpia_target_heated_read().

                                                                                          {

    memset( target, 0, sizeof( tpia_target_heated ) );
    return( 0 );
}

int tpia_target_heated_numberOfChannels	(	statusMessageReporting *	,
		tpia_target_heated *	target
	)

Definition at line 248 of file tpia_target_heated.cc.

References tpia_target_heated_s::nChannels.

Referenced by tpia_target_numberOfChannels().

                                                                                                {

    return( target->nChannels );
}

int tpia_target_heated_numberOfProductionChannels	(	statusMessageReporting *	,
		tpia_target_heated *	target
	)

Definition at line 256 of file tpia_target_heated.cc.

References tpia_target_heated_s::nProductionChannels.

Referenced by tpia_target_numberOfProductionChannels().

                                                                                                          {

    return( target->nProductionChannels );
}

int tpia_target_heated_read	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		const char *	fileName,
		int	checkElememtsForAccess
	)

Definition at line 126 of file tpia_target_heated.cc.

References tpia_target_heated_s::absPath, xData_element_s::attributes, tpia_target_heated_s::attributes, tpia_target_heated_s::channels, tpia_target_heated_s::contents, tpia_channel_s::crossSectionGrouped, tpia_channel_s::crossSectionPointwise, tpia_1dData_s::data, int(), xData_element_s::name, tpia_target_heated_s::nChannels, tpia_target_heated_s::nGroups, tpia_target_heated_s::nProductionChannels, tpia_target_heated_s::path, tpia_target_heated_s::productionChannels, tpia_target_heated_s::projectileID, smr_isOk(), tpia_1dData_s::start, tpia_target_heated_s::targetID, tpia_target_heated_s::totalCrossSectionGrouped, tpia_target_heated_s::totalCrossSectionPointwise, tpia_channel_createGetFromElement(), tpia_misc_pointerToAttributeIfAllOk3, tpia_misc_setMessageError_Element(), tpia_particle_getInternalID(), tpia_target_heated_initialize(), tpia_target_heated_release(), xData_addToAccessed(), xData_convertAttributeTo_xData_Int(), xData_copyAttributionList(), xData_getAttributesValue(), xData_getDocumentsElement(), xData_getFirstElement(), xData_getNextElement(), xData_getOneElementByTagName(), xData_malloc2, xData_numberOfElementsByTagName(), xData_parseFree(), xData_parseReadFile(), xDataMisc_allocateCopyString2, and xDataMisc_getAbsPath().

Referenced by tpia_target_heated_createRead().

                                                                                                                                         {
/*
*   If a target has already been read into this target, user must have called tpia_target_heated_release before calling this routine.
*   Otherwise, there will be memory leaks.
*/
    xData_document *doc = NULL;
    xData_element *element, *channelElement, *channels;
    int nChannels;
    tpia_channel *channel;
    char *name;
    xData_Int i, j;

    tpia_target_heated_initialize( smr, target );
    if( ( target->path = xDataMisc_allocateCopyString2( smr, fileName, "path" ) ) == NULL ) return( 1 );
    if( ( target->absPath = xDataMisc_getAbsPath( smr, fileName ) ) == NULL ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
    if( ( doc = xData_parseReadFile( smr, fileName, NULL, NULL ) ) == NULL ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
    element = xData_getDocumentsElement( doc );
    xData_addToAccessed( smr, element, 1 );
    if( xData_convertAttributeTo_xData_Int( smr, element, "nGroups", &i ) != 0 ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
    target->nGroups = (int) i;
    if( strcmp( element->name, "xTargetHeated" ) != 0 ) {
        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "input file's top element must be xTargetHeated and not %s", 
            element->name ); }
    else {
        xData_copyAttributionList( smr, &(target->attributes), &(element->attributes) );
        if( smr_isOk( smr ) ) target->contents = xData_getAttributesValue( &(target->attributes), "contents" );
        if( ( name = tpia_misc_pointerToAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "projectile" ) ) != NULL )
            target->projectileID = tpia_particle_getInternalID( smr, name );
        if( ( name = tpia_misc_pointerToAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "target" ) ) != NULL )
            target->targetID = tpia_particle_getInternalID( smr, name );
        if( smr_isOk( smr ) ) _tpia_target_heated_getEnergyGridAndAllocateTotalCrossSections( smr, target, element );
        if( smr_isOk( smr ) ) {            /* Get channels. */
            //if( ( channels = xData_getOneElementByTagName( smr, element, "channels", 1 ) ) == NULL ) 
            if( ( channels = xData_getOneElementByTagName( smr, element, (char*)"channels", 1 ) ) == NULL ) 
                return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
            xData_addToAccessed( smr, channels, 1 );
            if( ( nChannels = xData_numberOfElementsByTagName( smr, channels, "channel" ) ) > 0 ) {
                //if( ( target->channels = xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) == NULL ) 
                if( ( target->channels = (tpia_channel**) xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) == NULL ) 
                    return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
                for( channelElement = xData_getFirstElement( channels ); channelElement != NULL; channelElement = xData_getNextElement( channelElement ) ) {
                    if( !strcmp( channelElement->name, "channel" ) ) {
                        if( ( channel = tpia_channel_createGetFromElement( smr, target, channelElement, 1 ) ) == NULL ) break;
                        target->channels[target->nChannels] = channel;
                        target->nChannels++;
                        for( i = channel->crossSectionPointwise.start, j = 0; i < channel->crossSectionPointwise.end; i++, j++ ) 
                            target->totalCrossSectionPointwise.data[i] += channel->crossSectionPointwise.data[j];
                        for( i = channel->crossSectionGrouped.start, j = 0; i < channel->crossSectionGrouped.end; i++, j++ ) 
                            target->totalCrossSectionGrouped.data[i] += channel->crossSectionGrouped.data[j];
                    }
                }
            }
        }
        if( smr_isOk( smr ) ) {            /* Get production channels. */
            //if( ( channels = xData_getOneElementByTagName( smr, element, "productionChannels", 0 ) ) == NULL ) {
            if( ( channels = xData_getOneElementByTagName( smr, element, (char*) "productionChannels", 0 ) ) == NULL ) {
                if( !smr_isOk( smr ) ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) ); }
            else {
                xData_addToAccessed( smr, channels, 1 );
                if( ( nChannels = xData_numberOfElementsByTagName( smr, channels, "channel" ) ) > 0 ) {
                    //if( ( target->productionChannels = xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) != NULL ) {
                    if( ( target->productionChannels = (tpia_channel**) xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) != NULL ) {
                        for( channelElement = xData_getFirstElement(channels); channelElement != NULL; channelElement = xData_getNextElement(channelElement) ) {
                            if( !strcmp( channelElement->name, "channel" ) ) {
                                channel = tpia_channel_createGetFromElement( smr, target, channelElement, 1 );
                                if( channel == NULL ) break;
                                target->productionChannels[target->nProductionChannels] = channel;
                                target->nProductionChannels++;
                            }
                        }
                    }
                }
            }
        }
    }
    if( smr_isOk( smr ) && checkElememtsForAccess ) _tpia_target_heated_checkElememtsForAccess( smr, doc );
    xData_parseFree( smr, doc );
    if( !smr_isOk( smr ) ) tpia_target_heated_release( smr, target );
    return( !smr_isOk( smr ) );
}

int tpia_target_heated_release	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 95 of file tpia_target_heated.cc.

References tpia_target_heated_s::absPath, tpia_target_heated_s::attributes, tpia_target_heated_s::channels, tpia_1dData_s::data, tpia_target_heated_s::energyGrid, tpia_target_heated_s::energyGridLength, tpia_target_heated_s::kerma, tpia_target_heated_s::nChannels, tpia_target_heated_s::nProductionChannels, tpia_target_heated_s::path, tpia_target_heated_s::productionChannels, tpia_target_heated_s::totalCrossSectionGrouped, tpia_target_heated_s::totalCrossSectionPointwise, tpia_channel_free(), xData_free(), and xData_releaseAttributionList().

Referenced by tpia_target_heated_free(), and tpia_target_heated_read().

                                                                                          {

    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 );
    target->energyGridLength = 0;
    //target->energyGrid = xData_free( smr, target->energyGrid );
    target->energyGrid = (double*) xData_free( smr, target->energyGrid );
    //target->kerma = xData_free( smr, target->kerma );
    target->kerma = (double*) xData_free( smr, target->kerma );
    //target->totalCrossSectionPointwise.data = xData_free( smr, target->totalCrossSectionPointwise.data );
    target->totalCrossSectionPointwise.data = (double*) xData_free( smr, target->totalCrossSectionPointwise.data );
    //target->totalCrossSectionGrouped.data = xData_free( smr, target->totalCrossSectionGrouped.data );
    target->totalCrossSectionGrouped.data = (double*) xData_free( smr, target->totalCrossSectionGrouped.data );
    xData_releaseAttributionList( smr, &(target->attributes) );
    for( i = 0; i < target->nChannels; i++ ) tpia_channel_free( smr, target->channels[i] );
    target->nChannels = 0;
    //target->channels = xData_free( smr, target->channels );
    target->channels = (tpia_channel**) xData_free( smr, target->channels );
    for( i = 0; i < target->nProductionChannels; i++ ) tpia_channel_free( smr, target->productionChannels[i] );
    target->nProductionChannels = 0;
    //target->productionChannels = xData_free( smr, target->productionChannels );
    target->productionChannels = (tpia_channel**) xData_free( smr, target->productionChannels );
    return( 0 );
}

int tpia_target_heated_sampleIndexChannelProductsAtE	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		int	index,
		tpia_decaySamplingInfo *	decaySamplingInfo,
		int	nProductData,
		tpia_productOutgoingData *	productDatas
	)

Definition at line 341 of file tpia_target_heated.cc.

References tpia_channel_s::decayChannel, tpia_decayChannel_sampleProductsAtE(), and tpia_target_heated_getChannelAtIndex_smr().

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

                                                                                                              {

    tpia_channel *channel = tpia_target_heated_getChannelAtIndex_smr( smr, target, index );

    if( channel == NULL ) return( -1 );
    return( tpia_decayChannel_sampleProductsAtE( smr, &(channel->decayChannel), decaySamplingInfo, nProductData, productDatas ) );
}