tpia_channel.cc File Reference

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

Go to the source code of this file.

Functions
tpia_channel *	tpia_channel_create (statusMessageReporting *smr)
int	tpia_channel_initialize (statusMessageReporting *, tpia_channel *channel)
tpia_channel *	tpia_channel_createGetFromElement (statusMessageReporting *smr, tpia_target_heated *target, xData_element *channelElement, int pointwiseRequired)
tpia_channel *	tpia_channel_free (statusMessageReporting *smr, tpia_channel *channel)
int	tpia_channel_release (statusMessageReporting *smr, tpia_channel *channel)
int	tpia_channel_getFromElement (statusMessageReporting *smr, tpia_target_heated *target, xData_element *channelElement, tpia_channel *channel, int pointwiseRequired)
tpia_product *	tpia_channel_getFirstProduct (tpia_channel *channel)
tpia_product *	tpia_channel_getProductByIndex (statusMessageReporting *, tpia_channel *channel, int index)
int	tpia_channel_numberOfProducts (statusMessageReporting *, tpia_channel *channel)
int	tpia_channel_isProduction (statusMessageReporting *, tpia_channel *channel)
double	tpia_channel_getCrossSectionAtE (statusMessageReporting *smr, tpia_channel *channel, xData_Int, double e_in, int crossSectionType)

Function Documentation

tpia_channel* tpia_channel_create

(

statusMessageReporting *

smr

)

Definition at line 50 of file tpia_channel.cc.

References tpia_channel_free(), tpia_channel_initialize(), and xData_malloc2.

Referenced by tpia_channel_createGetFromElement().

                                                                 {

    tpia_channel *channel;

    //if( ( channel = xData_malloc2( smr, sizeof( tpia_channel ), 0, "channel" ) ) == NULL ) return( NULL );
    if( ( channel = (tpia_channel*) xData_malloc2( smr, sizeof( tpia_channel ), 0, "channel" ) ) == NULL ) return( NULL );
    if( tpia_channel_initialize( smr, channel ) ) channel = tpia_channel_free( smr, channel );
    return( channel );
}

tpia_channel* tpia_channel_createGetFromElement	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		xData_element *	channelElement,
		int	pointwiseRequired
	)

Definition at line 71 of file tpia_channel.cc.

References tpia_channel_create(), tpia_channel_free(), and tpia_channel_getFromElement().

Referenced by tpia_target_heated_read().

                            {

    tpia_channel *channel;

    if( ( channel = tpia_channel_create( smr ) ) == NULL ) return( NULL );
    if( tpia_channel_getFromElement( smr, target, channelElement, channel, pointwiseRequired ) != 0 ) channel = tpia_channel_free( smr, channel );
    return( channel );
}

tpia_channel* tpia_channel_free	(	statusMessageReporting *	smr,
		tpia_channel *	channel
	)

Definition at line 83 of file tpia_channel.cc.

References tpia_channel_release(), and xData_free().

Referenced by tpia_channel_create(), tpia_channel_createGetFromElement(), and tpia_target_heated_release().

                                                                                      {

    tpia_channel_release( smr, channel );
    xData_free( smr, channel );
    return( NULL );
}

double tpia_channel_getCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_channel *	channel,
		xData_Int	,
		double	e_in,
		int	crossSectionType
	)

Definition at line 222 of file tpia_channel.cc.

References tpia_channel_s::crossSectionPointwise, tpia_target_heated_s::energyGrid, tpia_channel_s::target, tpia_crossSectionType_grouped, tpia_crossSectionType_pointwise, tpia_misc_getPointwiseCrossSectionAtE(), and tpia_target_heated_getEIndex().

Referenced by tpia_target_heated_getIndexChannelCrossSectionAtE().

                               {

    double xsec = 0.;

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

tpia_product* tpia_channel_getFirstProduct

(

tpia_channel *

channel

)

Definition at line 185 of file tpia_channel.cc.

References tpia_channel_s::decayChannel, and tpia_decayChannel_getFirstProduct().

Referenced by G4GIDI_Misc_channelProductsCompare(), and tpia_channel_getProductByIndex().

                                                                    {

    return( tpia_decayChannel_getFirstProduct( &(channel->decayChannel) ) );
}

int tpia_channel_getFromElement	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		xData_element *	channelElement,
		tpia_channel *	channel,
		int	pointwiseRequired
	)

Definition at line 114 of file tpia_channel.cc.

References xData_element_s::attributes, tpia_channel_s::attributes, tpia_channel_s::availableEnergyGrouped, tpia_target_heated_s::contents, tpia_channel_s::crossSectionFrame, tpia_channel_s::crossSectionGrouped, tpia_channel_s::crossSectionPointwise, tpia_1dData_s::data, tpia_1dData_s::end, tpia_channel_s::ENDF_MT, tpia_channel_s::ENDL_C, tpia_target_heated_s::energyGrid, tpia_channel_s::fission, tpia_channel_s::genre, int(), tpia_channel_s::outputChannel, tpia_channel_s::Q, tpia_channel_s::QIsFloat, tpia_channel_s::QString, smr_isOk(), tpia_1dData_s::start, tpia_channel_s::target, tpia_channel_isProduction(), tpia_frame_setFromElement(), tpia_misc_get2d_xShared_yHistogram_data_Grouped(), tpia_misc_get2dxindex_y_data(), tpia_misc_pointerToAttributeIfAllOk2, xData_addToAccessed(), xData_convertAttributeTo_xData_Int(), xData_copyAttributionList(), and xData_getOneElementByTagName().

Referenced by tpia_channel_createGetFromElement().

                            {

    xData_Int ll;
    char *p;
    xData_element *element, *pElement, *gElement, *eElement;

    xData_addToAccessed( smr, channelElement, 1 );
    channel->target = target;
    xData_copyAttributionList( smr, &(channel->attributes), &(channelElement->attributes) );
    channel->outputChannel = tpia_misc_pointerToAttributeIfAllOk2(smr, channelElement, 1, &(channel->attributes), "outputChannel" );
    channel->genre = tpia_misc_pointerToAttributeIfAllOk2( smr, channelElement, 1, &(channel->attributes), "genre" );
    channel->QString = tpia_misc_pointerToAttributeIfAllOk2( smr, channelElement, 1, &(channel->attributes), "Q" );
    channel->fission = tpia_misc_pointerToAttributeIfAllOk2( smr, channelElement, 0, &(channel->attributes), "fission" );
    if( smr_isOk( smr ) ) {
        ll = 0;
        if( xData_convertAttributeTo_xData_Int( smr, channelElement, "ENDL_C", &ll ) >= 0 ) channel->ENDL_C = (int) ll;
    }
    if( smr_isOk( smr ) ) {
        ll = 0;
        if( xData_convertAttributeTo_xData_Int( smr, channelElement, "ENDF_MT2", &ll ) >= 0 ) channel->ENDF_MT = (int) ll;
    }
    if( smr_isOk( smr ) ) {
        channel->QIsFloat = 1;
        channel->Q = strtod( channel->QString, &p );      /* Q string may be something like "notApplicable". */
        if( *p != 0 ) {                                     /* In that case set QIsFloat to false. */
            channel->QIsFloat = 0;
            channel->Q = 0.;
        }
        //if( ( element = xData_getOneElementByTagName( smr, channelElement, "crossSection", 1 ) ) != NULL ) {
        if( ( element = xData_getOneElementByTagName( smr, channelElement, (char*) "crossSection", 1 ) ) != NULL ) {
            if( ( tpia_frame_setFromElement( smr, element, 2, &channel->crossSectionFrame ) ) == 0 ) {
                xData_addToAccessed( smr, element, 1 );
                //if( ( pElement = xData_getOneElementByTagName( smr, element, "indexed", 1 ) ) != NULL ) {
                if( ( pElement = xData_getOneElementByTagName( smr, element, (char*) "indexed", 1 ) ) != NULL ) {
                    channel->crossSectionPointwise.data = tpia_misc_get2dxindex_y_data( smr, pElement,
                        &(channel->crossSectionPointwise.start), &(channel->crossSectionPointwise.end), target->energyGrid );
                }
                if( ( gElement = xData_getOneElementByTagName( smr, element, (char*) "grouped", 1 ) ) != NULL ) {
                    tpia_misc_get2d_xShared_yHistogram_data_Grouped( smr, gElement, &(channel->crossSectionGrouped) );
                }
                if( ( channel->crossSectionGrouped.data != NULL ) && ( ( channel->crossSectionPointwise.data != NULL ) || !pointwiseRequired ) ) {
                    if( target->contents == NULL ) {                /* Only supported "crossSection" currently. */
                        if( !tpia_channel_isProduction( smr, channel ) ) {
                            //if( ( eElement = xData_getOneElementByTagName( smr, channelElement, "availableEnergy", 1 ) ) != NULL ) {
                            if( ( eElement = xData_getOneElementByTagName( smr, channelElement, (char*) "availableEnergy", 1 ) ) != NULL ) {
                                xData_addToAccessed( smr, eElement, 1 );
                                //if( ( gElement = xData_getOneElementByTagName( smr, eElement, "grouped", 1 ) ) != NULL ) {
                                if( ( gElement = xData_getOneElementByTagName( smr, eElement, (char*) "grouped", 1 ) ) != NULL ) {
                                    tpia_misc_get2d_xShared_yHistogram_data_Grouped( smr, gElement, &(channel->availableEnergyGrouped) );
                                }
                            }
                        }
                    }
                    if( smr_isOk( smr ) ) _tpia_channel_getProductData( smr, channelElement, channel );
                }
            }
        }
    }
    return( !smr_isOk( smr ) );
}

tpia_product* tpia_channel_getProductByIndex	(	statusMessageReporting *	,
		tpia_channel *	channel,
		int	index
	)

Definition at line 193 of file tpia_channel.cc.

References tpia_channel_getFirstProduct(), and tpia_decayChannel_getNextProduct().

                                                                                                           {

    int i = 0;
    tpia_product *p;

    if( index < 0 ) return( NULL );
    for( p = tpia_channel_getFirstProduct( channel ); ( p != NULL ) && ( i < index ); p = tpia_decayChannel_getNextProduct( p ), i++ ) ;
    return( p );
}

int tpia_channel_initialize	(	statusMessageReporting *	,
		tpia_channel *	channel
	)

Definition at line 63 of file tpia_channel.cc.

Referenced by tpia_channel_create().

                                                                               {

    memset( channel, 0, sizeof( tpia_channel ) );
    return( 0 );
}

int tpia_channel_isProduction	(	statusMessageReporting *	,
		tpia_channel *	channel
	)

Definition at line 214 of file tpia_channel.cc.

References tpia_channel_s::genre.

Referenced by tpia_channel_getFromElement().

                                                                                 {

    return( strcmp( channel->genre, "production" ) == 0 );
}

int tpia_channel_numberOfProducts	(	statusMessageReporting *	,
		tpia_channel *	channel
	)

Definition at line 206 of file tpia_channel.cc.

References tpia_channel_s::decayChannel, and tpia_decayChannel_s::numberOfProducts.

                                                                                     {

    return( channel->decayChannel.numberOfProducts );
}

int tpia_channel_release	(	statusMessageReporting *	smr,
		tpia_channel *	channel
	)

Definition at line 92 of file tpia_channel.cc.

References tpia_channel_s::attributes, tpia_channel_s::availableEnergyGrouped, tpia_channel_s::crossSectionGrouped, tpia_channel_s::crossSectionPointwise, tpia_1dData_s::data, tpia_channel_s::decayChannel, tpia_product_s::next, tpia_decayChannel_s::numberOfProducts, tpia_decayChannel_s::products, tpia_product_free(), xData_free(), and xData_releaseAttributionList().

Referenced by tpia_channel_free().

                                                                               {

    tpia_product *product, *nextProduct;

    xData_releaseAttributionList( smr, &(channel->attributes) );
    //channel->crossSectionPointwise.data = xData_free( smr, channel->crossSectionPointwise.data );
    channel->crossSectionPointwise.data = (double*) xData_free( smr, channel->crossSectionPointwise.data );
    //channel->crossSectionGrouped.data = xData_free( smr, channel->crossSectionGrouped.data );
    channel->crossSectionGrouped.data = (double*) xData_free( smr, channel->crossSectionGrouped.data );
    //channel->availableEnergyGrouped.data = xData_free( smr, channel->availableEnergyGrouped.data );
    channel->availableEnergyGrouped.data = (double*) xData_free( smr, channel->availableEnergyGrouped.data );
    for( product = channel->decayChannel.products; product != NULL; product = nextProduct ) {
        nextProduct = product->next;
        tpia_product_free( smr, product );
    }
    channel->decayChannel.numberOfProducts = 0;
    channel->decayChannel.products = NULL;
    return( 0 );
}