tpia_angularEnergy.cc File Reference

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

Go to the source code of this file.

Functions
int	tpia_angularEnergy_initialize (statusMessageReporting *smr, tpia_angularEnergy *angularEnergy)
int	tpia_angularEnergy_release (statusMessageReporting *smr, tpia_angularEnergy *angularEnergy)
int	tpia_angularEnergy_getFromElement (statusMessageReporting *smr, xData_element *angularEnergyElement, tpia_angularEnergy *angularEnergy)
int	tpia_angularEnergy_SampleEp (statusMessageReporting *smr, tpia_angularEnergy *angularEnergy, tpia_decaySamplingInfo *decaySamplingInfo)

Function Documentation

int tpia_angularEnergy_getFromElement	(	statusMessageReporting *	smr,
		xData_element *	angularEnergyElement,
		tpia_angularEnergy *	angularEnergy
	)

Definition at line 73 of file tpia_angularEnergy.cc.

References tpia_angularEnergy_s::binned, tpia_EqualProbableBinSpectra_s::dValue, xData_elementListItem_s::element, tpia_EqualProbableBinSpectra_s::energies, tpia_angularEnergyBin_s::energies, tpia_angularEnergy_s::frame, int(), xData_elementList_s::items, tpia_EqualProbableBinSpectra_s::iValue, n, xData_elementList_s::n, tpia_angularEnergyBin_s::nBins, tpia_EqualProbableBinSpectra_s::numberOfEs, tpia_angularEnergyBin_s::numberOfEs, status, tpia_frame_setFromElement(), tpia_misc_getEqualProbableBins(), tpia_misc_setMessageError_Element(), xData_addToAccessed(), xData_convertAttributeTo_xData_Int(), xData_convertAttributeToDouble(), xData_freeElementList(), xData_getElementsByTagNameAndSort(), xData_getOneElementByTagName(), and xData_malloc2.

                                                                                                                                             {

    int i, j, status = 1;
    xData_Int nBins, size, n, index;
    xData_element *epbElement, *element;
    xData_elementList *list;
    tpia_EqualProbableBinSpectra *mus;

    xData_addToAccessed( smr, angularEnergyElement, 1 );
    if( ( tpia_frame_setFromElement( smr, angularEnergyElement, 4, &(angularEnergy->frame) ) ) != 0 ) return( 1 );
    //if( ( epbElement = xData_getOneElementByTagName( smr, angularEnergyElement, "equalProbableBins", 0 ) ) == NULL ) return( 1 );
    if( ( epbElement = xData_getOneElementByTagName( smr, angularEnergyElement, (char*)"equalProbableBins", 0 ) ) == NULL ) return( 1 );
    xData_addToAccessed( smr, epbElement, 1 );
    if( xData_convertAttributeTo_xData_Int( smr, epbElement, "nBins", &nBins ) != 0 ) {
        tpia_misc_setMessageError_Element( smr, NULL, epbElement, __FILE__, __LINE__, 1, "missing or invalid nBins attribute" );
        return( 1 );
    }
    list = xData_getElementsByTagNameAndSort( smr, epbElement, "energy", NULL, NULL );
    if( list->n == 0 ) {
        tpia_misc_setMessageError_Element( smr, NULL, epbElement, __FILE__, __LINE__, 1, "bins does not contain any energy elements" ); }
    else {
        size = list->n * sizeof( tpia_EqualProbableBinSpectra );
        angularEnergy->binned.nBins = (int) nBins;
        //if( ( angularEnergy->binned.energies = xData_malloc2( smr, size, 1, "energies" ) ) != NULL ) {
        if( ( angularEnergy->binned.energies = (tpia_EqualProbableBinSpectra*) xData_malloc2( smr, size, 1, "energies" ) ) != NULL ) {
            status = 0;
            for( i = 0, mus = angularEnergy->binned.energies; i < list->n; i++, mus++ ) {
                mus->iValue = 0;
                element = list->items[i].element;
                if( xData_convertAttributeTo_xData_Int( smr, element, "index", &index ) != 0 ) {
                    tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "missing or invalid index attribute" );
                    status = 1;
                    break;
                }
                if( index != i ) {
                    tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "index = %lld is not incremental", index );
                    status = 1;
                    break;
                }
                if( ( j = xData_convertAttributeToDouble( smr, element, "value", &(mus->dValue) ) ) != 0 ) {
                    if( j == 1 ) {
                        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "element does not have value attribute" ); }
                    else {
                        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "failed to convert value attribute to double" );
                    }
                    status = 1;
                    break;
                }
                if( ( mus->energies = tpia_misc_getEqualProbableBins( smr, element, "mu", nBins, &n ) ) == NULL ) {
                    status = 1;
                    break; }
                else {
                    mus->numberOfEs = n;
                }
                angularEnergy->binned.numberOfEs++;
            }
        }
    }
    xData_freeElementList( smr, list );
    return( status );
}

int tpia_angularEnergy_initialize	(	statusMessageReporting *	smr,
		tpia_angularEnergy *	angularEnergy
	)

Definition at line 49 of file tpia_angularEnergy.cc.

References tpia_angularEnergy_s::binned, tpia_angularEnergyBin_s::energies, tpia_angularEnergy_s::frame, tpia_angularEnergyBin_s::numberOfEs, and tpia_frame_setFromString().

Referenced by tpia_angularEnergy_release().

                                                                                                    {

    memset( angularEnergy, 0, sizeof( tpia_angularEnergy ) );
    if( tpia_frame_setFromString( smr, "", "", 0, &(angularEnergy->frame) ) ) return( 1 );
    angularEnergy->binned.numberOfEs = 0;
    angularEnergy->binned.energies = NULL;
    return( 0 );
}

int tpia_angularEnergy_release	(	statusMessageReporting *	smr,
		tpia_angularEnergy *	angularEnergy
	)

Definition at line 60 of file tpia_angularEnergy.cc.

References tpia_angularEnergy_s::binned, tpia_EqualProbableBinSpectra_s::energies, tpia_angularEnergyBin_s::energies, tpia_angularEnergyBin_s::numberOfEs, tpia_angularEnergy_initialize(), and xData_free().

Referenced by tpia_product_release().

                                                                                                 {

    int i;

    for( i = 0; i < angularEnergy->binned.numberOfEs; i++ ) xData_free( smr, angularEnergy->binned.energies[i].energies );
    //angularEnergy->binned.energies = xData_free( smr, angularEnergy->binned.energies );
    angularEnergy->binned.energies = (tpia_EqualProbableBinSpectra*) xData_free( smr, angularEnergy->binned.energies );
    tpia_angularEnergy_initialize( smr, angularEnergy );
    return( 0 );
}

int tpia_angularEnergy_SampleEp	(	statusMessageReporting *	smr,
		tpia_angularEnergy *	angularEnergy,
		tpia_decaySamplingInfo *	decaySamplingInfo
	)

Definition at line 137 of file tpia_angularEnergy.cc.

References tpia_angularEnergy_s::binned, tpia_EqualProbableBinSpectra_s::dValue, tpia_decaySamplingInfo_s::e_in, tpia_angularEnergyBin_s::energies, tpia_decaySamplingInfo_s::Ep, tpia_decaySamplingInfo_s::mu, tpia_angularEnergyBin_s::nBins, tpia_angularEnergyBin_s::numberOfEs, smr_isOk(), and tpia_misc_sampleEqualProbableBin().

Referenced by tpia_decayChannel_sampleProductsAtE().

                                                                                                                                             {

    int iE1, iE2;
    tpia_EqualProbableBinSpectra *energies = angularEnergy->binned.energies;
    double Ep, Ep1, Ep2, f, e_in = decaySamplingInfo->e_in;
/*
Currently, only equal probable binning is supported.
Need to return frame info for Ep also.
*/

    if( !smr_isOk( smr ) ) return( 1 );
    if( angularEnergy->binned.numberOfEs == 0 ) return( 1 );
    for( iE2 = 0; iE2 < angularEnergy->binned.numberOfEs; iE2++ ) if( energies[iE2].dValue >= e_in ) break;
    if( iE2 == 0 ) {
        iE1 = iE2; }
    else if( iE2 == angularEnergy->binned.numberOfEs ) {
        iE1 = iE2 = angularEnergy->binned.numberOfEs - 1; }
    else {
        iE1 = iE2 - 1;
    }
    tpia_misc_sampleEqualProbableBin( smr, decaySamplingInfo, decaySamplingInfo->mu, angularEnergy->binned.nBins, &(energies[iE1]), &Ep1 );
    if( iE1 == iE2 ) {
        Ep = Ep1; }
    else {
        tpia_misc_sampleEqualProbableBin( smr, decaySamplingInfo, decaySamplingInfo->mu, angularEnergy->binned.nBins, &(energies[iE2]), &Ep2 );
        f = ( energies[iE2].dValue - e_in ) / ( energies[iE2].dValue - energies[iE1].dValue );
        Ep = f * Ep1 + ( 1 - f ) * Ep2;
    }
    decaySamplingInfo->Ep = Ep;

    return( 0 );
}