tpia_product.cc File Reference

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

Go to the source code of this file.

Functions
tpia_product *	tpia_product_create (statusMessageReporting *smr)
int	tpia_product_initialize (statusMessageReporting *smr, tpia_product *product)
tpia_product *	tpia_product_createGetFromElement (statusMessageReporting *smr, tpia_channel *channel, tpia_product *parentProduct, xData_element *productElement)
tpia_product *	tpia_product_free (statusMessageReporting *smr, tpia_product *product)
int	tpia_product_release (statusMessageReporting *smr, tpia_product *product)
int	tpia_product_getFromElement (statusMessageReporting *smr, tpia_channel *channel, tpia_product *parentProduct, xData_element *productElement, tpia_product *product)
int	tpia_product_getDecayChannelFromElement (statusMessageReporting *smr, xData_element *parentElement, tpia_channel *channel, tpia_product *parentProduct, tpia_product **priorProductNext)
long	tpia_product_dataRequired (statusMessageReporting *, tpia_product *product)
tpia_product *	tpia_product_getFirstProduct (tpia_product *product)
tpia_product *	tpia_product_getProductByIndex (statusMessageReporting *, tpia_product *product, int index)
int	tpia_product_doesDecay (statusMessageReporting *, tpia_product *product)
int	tpia_product_numberOfProducts (statusMessageReporting *, tpia_product *product)
int	tpia_product_isDataPresent (statusMessageReporting *, tpia_product *product, int b_data)
int	tpia_product_sampleMultiplicity (statusMessageReporting *, tpia_product *product, double e_in, double r)

Variables
const char *	tpia_productGenre_unknown = "unknown"
const char *	tpia_productGenre_twoBody_angular = "twoBody_angular"
const char *	tpia_productGenre_twoBody_formFactor = "twoBody_formFactor"
const char *	tpia_productGenre_NBody_Legendre = "NBody_Legendre"
const char *	tpia_productGenre_NBody_angular_energy = "NBody_angular_energy"
const char *	tpia_productGenre_NBody_uncorrelate_Legendre = "NBody_uncorrelate_Legendre"
const char *	tpia_productGenre_NBody_pairProduction = "NBody_pairProduction"

Function Documentation

tpia_product* tpia_product_create

(

statusMessageReporting *

smr

)

Definition at line 71 of file tpia_product.cc.

References tpia_product_free(), tpia_product_initialize(), and xData_malloc2.

Referenced by tpia_product_createGetFromElement().

                                                                 {

    tpia_product *product;

    //if( ( product = xData_malloc2( smr, sizeof( tpia_product ), 0, "product" ) ) == NULL ) return( NULL );
    if( ( product = (tpia_product*) xData_malloc2( smr, sizeof( tpia_product ), 0, "product" ) ) == NULL ) return( NULL );
    if( tpia_product_initialize( smr, product ) ) product = tpia_product_free( smr, product );
    return( product );
}

tpia_product* tpia_product_createGetFromElement	(	statusMessageReporting *	smr,
		tpia_channel *	channel,
		tpia_product *	parentProduct,
		xData_element *	productElement
	)

Definition at line 93 of file tpia_product.cc.

References tpia_product_create(), tpia_product_free(), and tpia_product_getFromElement().

Referenced by tpia_product_getDecayChannelFromElement().

                                    {

    tpia_product *product;

    if( ( product = tpia_product_create( smr ) ) == NULL ) return( NULL );
    if( tpia_product_getFromElement( smr, channel, parentProduct, productElement, product ) != 0 ) product = tpia_product_free( smr, product );
    return( product );
}

long tpia_product_dataRequired	(	statusMessageReporting *	,
		tpia_product *	product
	)

Definition at line 359 of file tpia_product.cc.

References tpia_product_s::b_dataRequired.

                                                                                  {

    return( product->b_dataRequired );
}

int tpia_product_doesDecay	(	statusMessageReporting *	,
		tpia_product *	product
	)

Definition at line 387 of file tpia_product.cc.

References tpia_product_s::decayChannel, and tpia_decayChannel_s::products.

                                                                              {

    return( product->decayChannel.products != NULL );
}

tpia_product* tpia_product_free	(	statusMessageReporting *	smr,
		tpia_product *	product
	)

Definition at line 105 of file tpia_product.cc.

References tpia_product_release(), and xData_free().

Referenced by tpia_channel_release(), tpia_product_create(), tpia_product_createGetFromElement(), and tpia_product_release().

                                                                                      {

    tpia_product_release( smr, product );
    xData_free( smr, product );
    return( NULL );
}

int tpia_product_getDecayChannelFromElement	(	statusMessageReporting *	smr,
		xData_element *	parentElement,
		tpia_channel *	channel,
		tpia_product *	parentProduct,
		tpia_product **	priorProductNext
	)

Definition at line 268 of file tpia_product.cc.

References tpia_product_s::decayChannel, tpia_channel_s::decayChannel, xData_elementListItem_s::element, tpia_particle_s::fullMass_MeV, xData_elementList_s::items, tpia_decayChannel_s::m1_fullMass_MeV, tpia_decayChannel_s::m2_fullMass_MeV, xData_elementList_s::n, tpia_product_s::next, tpia_decayChannel_s::numberOfProducts, tpia_product_s::productID, tpia_target_heated_s::projectileID, status, tpia_channel_s::target, tpia_target_heated_s::targetID, tpia_product_createGetFromElement(), xData_freeElementList(), and xData_getElementsByTagName().

                                      {

    xData_elementList *list;
    tpia_product *product;
    int i, status = 0;

    list = xData_getElementsByTagName( smr, parentElement, "product" );
    for( i = 0; i < list->n; i++ ) {
        if( ( product = tpia_product_createGetFromElement( smr, channel, parentProduct, list->items[i].element ) ) == NULL ) {
            status = 1;
            break;
        }
        if( parentProduct == NULL ) {
            channel->decayChannel.m1_fullMass_MeV = channel->target->projectileID->fullMass_MeV;
            channel->decayChannel.m2_fullMass_MeV = channel->target->targetID->fullMass_MeV;
            channel->decayChannel.numberOfProducts++; }
        else {
            channel->decayChannel.m1_fullMass_MeV = parentProduct->productID->fullMass_MeV;
            channel->decayChannel.m2_fullMass_MeV = 0.;
            parentProduct->decayChannel.numberOfProducts++;
        }
        *priorProductNext = product;
        priorProductNext = &(product->next);
    }
    xData_freeElementList( smr, list );
    return( status );
}

tpia_product* tpia_product_getFirstProduct

(

tpia_product *

product

)

Definition at line 366 of file tpia_product.cc.

References tpia_product_s::decayChannel, and tpia_decayChannel_getFirstProduct().

Referenced by tpia_product_getProductByIndex().

                                                                    {

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

int tpia_product_getFromElement	(	statusMessageReporting *	smr,
		tpia_channel *	channel,
		tpia_product *	parentProduct,
		xData_element *	productElement,
		tpia_product *	product
	)

Definition at line 142 of file tpia_product.cc.

References xData_element_s::attributes, tpia_product_s::attributes, tpia_product_s::b_dataRequired, tpia_product_s::channel, tpia_product_s::genre, tpia_product_s::multiplicity, tpia_product_s::parentProduct, tpia_product_s::productID, smr_isOk(), tpia_misc_pointerToAttributeIfAllOk2, tpia_misc_setMessageError_Element(), tpia_particle_getInternalID(), tpia_productGenre_NBody_angular_energy, tpia_productGenre_NBody_Legendre, tpia_productGenre_NBody_pairProduction, tpia_productGenre_NBody_uncorrelate_Legendre, tpia_productGenre_twoBody_angular, tpia_productGenre_twoBody_formFactor, tpia_productGenre_unknown, xData_addToAccessed(), and xData_copyAttributionList().

Referenced by tpia_product_createGetFromElement().

                            {

    char const *productGenre;
    char *name, *multiplicity, *e;

    xData_addToAccessed( smr, productElement, 1 );
    product->channel = channel;
    product->parentProduct = parentProduct;
    if( xData_copyAttributionList( smr, &(product->attributes), &(productElement->attributes) ) != 0 ) return( 0 );
    name = tpia_misc_pointerToAttributeIfAllOk2( smr, productElement, 1, &(product->attributes), "particle" );
    if( name != NULL ) {
        product->productID = tpia_particle_getInternalID( smr, name );
        multiplicity = tpia_misc_pointerToAttributeIfAllOk2( smr, productElement, 1, &(product->attributes), "multiplicity" );
        if( multiplicity != NULL ) {
            if( strcmp( multiplicity, "energyDependent" ) && strcmp( multiplicity, "partialProduction" ) ) {    /* Must be an integer. */
                product->multiplicity = strtol( multiplicity, &e, 10 );
                while( isspace( *e ) ) e++;
                if( *e != 0 ) tpia_misc_setMessageError_Element( smr, NULL, productElement, __FILE__, __LINE__, 1, "bad multiplicity = %s", multiplicity );
            }
        }
    }
    if( ( productGenre = tpia_misc_pointerToAttributeIfAllOk2( smr, productElement, 1, &(product->attributes), "genre" ) ) != NULL ) {
        if( strcmp( productGenre, tpia_productGenre_unknown ) == 0 ) {
            product->b_dataRequired = 0;
            product->genre = tpia_productGenre_unknown; }
        else if( strcmp( productGenre, tpia_productGenre_twoBody_angular ) == 0 ) {
            product->b_dataRequired = tpia_b_twoBody_angular;
            product->genre = tpia_productGenre_twoBody_angular; }
        else if( strcmp( productGenre, tpia_productGenre_twoBody_formFactor ) == 0 ) {
            product->b_dataRequired = tpia_b_twoBody_formFactor;
            product->genre = tpia_productGenre_twoBody_formFactor; }
        else if( strcmp( productGenre, tpia_productGenre_NBody_Legendre ) == 0 ) {
            product->b_dataRequired = tpia_b_NBody_Legendre;
            product->genre = tpia_productGenre_NBody_Legendre; }
        else if( strcmp( productGenre, tpia_productGenre_NBody_angular_energy ) == 0 ) {
            product->b_dataRequired = tpia_b_NBody_angular_energy;
            product->genre = tpia_productGenre_NBody_angular_energy; }
        else if( strcmp( productGenre, tpia_productGenre_NBody_uncorrelate_Legendre ) == 0 ) {
            product->b_dataRequired = tpia_b_NBody_uncorrelate_Legendre;
            product->genre = tpia_productGenre_NBody_uncorrelate_Legendre; }
        else if( strcmp( productGenre, tpia_productGenre_NBody_pairProduction ) == 0 ) {
            product->b_dataRequired = tpia_b_NBody_pairProduction;
            product->genre = tpia_productGenre_NBody_pairProduction; }
        else {
            tpia_misc_setMessageError_Element( smr, NULL, productElement, __FILE__, __LINE__, 1, "unsupported product genre = %s", productGenre );
        }
        if( smr_isOk( smr ) ) _tpia_product_getProductOutgoingData( smr, productElement, product );
    }
    return( !smr_isOk( smr ) );
}

tpia_product* tpia_product_getProductByIndex	(	statusMessageReporting *	,
		tpia_product *	product,
		int	index
	)

Definition at line 374 of file tpia_product.cc.

References tpia_decayChannel_getNextProduct(), and tpia_product_getFirstProduct().

                                                                                                           {

    int i = 0;
    tpia_product *p;

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

int tpia_product_initialize	(	statusMessageReporting *	smr,
		tpia_product *	product
	)

Definition at line 83 of file tpia_product.cc.

References tpia_product_s::angular, tpia_product_s::Legendre, tpia_angular_initialize(), and tpia_Legendre_initialize().

Referenced by tpia_product_create().

                                                                                  {

    memset( product, 0, sizeof( tpia_product ) );
    if( tpia_angular_initialize( smr, &(product->angular) ) ) return( 1 );
    if( tpia_Legendre_initialize( smr, &(product->Legendre) ) ) return( 1 );
    return( 0 );
}

int tpia_product_isDataPresent	(	statusMessageReporting *	,
		tpia_product *	product,
		int	b_data
	)

Definition at line 403 of file tpia_product.cc.

References tpia_product_s::b_dataPresent.

                                                                                              {

    return( product->b_dataPresent && b_data );
}

int tpia_product_numberOfProducts	(	statusMessageReporting *	,
		tpia_product *	product
	)

Definition at line 395 of file tpia_product.cc.

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

                                                                                     {

    return( product->decayChannel.numberOfProducts );
}

int tpia_product_release	(	statusMessageReporting *	smr,
		tpia_product *	product
	)

Definition at line 114 of file tpia_product.cc.

References tpia_product_s::angular, tpia_product_s::angularEnergy, tpia_product_s::attributes, tpia_1dData_s::data, tpia_product_s::decayChannel, tpia_product_s::delayedNeutronMultiplicityVsEnergy, tpia_product_s::depositionEnergyGrouped, tpia_product_s::Legendre, tpia_product_s::multiplicityVsEnergy, tpia_multiplicity_s::next, tpia_product_s::next, tpia_decayChannel_s::numberOfProducts, tpia_decayChannel_s::products, tpia_angular_release(), tpia_angularEnergy_release(), tpia_Legendre_release(), tpia_multiplicity_free(), tpia_product_free(), xData_free(), and xData_releaseAttributionList().

Referenced by tpia_product_free().

                                                                              {

    tpia_multiplicity *multiplicity, *multiplicity_next;
    tpia_product *decayProduct, *nextProduct;

    xData_releaseAttributionList( smr, &(product->attributes) );
    //product->depositionEnergyGrouped.data = xData_free( smr, product->depositionEnergyGrouped.data );
    product->depositionEnergyGrouped.data = (double*) xData_free( smr, product->depositionEnergyGrouped.data );

    if( product->multiplicityVsEnergy != NULL ) tpia_multiplicity_free( smr, product->multiplicityVsEnergy );
    for( multiplicity = product->delayedNeutronMultiplicityVsEnergy; multiplicity != NULL; multiplicity = multiplicity_next ) {
        multiplicity_next = multiplicity->next;
        tpia_multiplicity_free( smr, multiplicity );
    }
    tpia_angular_release( smr, &(product->angular) );
    tpia_Legendre_release( smr, &(product->Legendre ) );
    tpia_angularEnergy_release( smr, &(product->angularEnergy) );
    for( decayProduct = product->decayChannel.products; decayProduct != NULL; decayProduct = nextProduct ) {
        nextProduct = decayProduct->next;
        tpia_product_free( smr, decayProduct );
    }
    product->decayChannel.numberOfProducts = 0;
    product->decayChannel.products = NULL;
    return( 0 );
}

int tpia_product_sampleMultiplicity	(	statusMessageReporting *	,
		tpia_product *	product,
		double	e_in,
		double	r
	)

Definition at line 411 of file tpia_product.cc.

References int(), tpia_product_s::multiplicityVsEnergy, tpia_multiplicity_s::numberOfPointwise, and tpia_multiplicity_s::pointwise.

Referenced by tpia_decayChannel_sampleProductsAtE().

                                                                                                              {

    int i, multiplicity;
    tpia_multiplicity *multiplicityVsEnergy = product->multiplicityVsEnergy;
    double *p = multiplicityVsEnergy->pointwise, dMult;

    if( e_in <= p[0] ) {
        dMult = p[1]; }
    else if( e_in >= p[2 * ( multiplicityVsEnergy->numberOfPointwise - 1 )] ) {
        dMult = p[2 * multiplicityVsEnergy->numberOfPointwise - 1]; }
    else {
        for( i = 0; i < multiplicityVsEnergy->numberOfPointwise - 1; i++, p += 2 ) if( e_in < p[2] ) break;
        dMult = ( e_in - p[0] ) / ( p[2] - p[0] );
        dMult = dMult * p[3] + ( 1. - dMult ) * p[1];
    }
    multiplicity = (int) dMult;
    if( r < ( dMult - multiplicity ) ) multiplicity++;

    return( multiplicity );
}

Variable Documentation

const char * tpia_productGenre_NBody_angular_energy = "NBody_angular_energy"

Definition at line 60 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

const char * tpia_productGenre_NBody_Legendre = "NBody_Legendre"

Definition at line 59 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

const char * tpia_productGenre_NBody_pairProduction = "NBody_pairProduction"

Definition at line 62 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

const char * tpia_productGenre_NBody_uncorrelate_Legendre = "NBody_uncorrelate_Legendre"

Definition at line 61 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

const char * tpia_productGenre_twoBody_angular = "twoBody_angular"

Definition at line 57 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

const char * tpia_productGenre_twoBody_formFactor = "twoBody_formFactor"

Definition at line 58 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

const char* tpia_productGenre_unknown = "unknown"

Definition at line 56 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().