tpia_Legendre.cc File Reference

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

Go to the source code of this file.

Functions
int	tpia_Legendre_initialize (statusMessageReporting *smr, tpia_Legendre *Legendre)
int	tpia_Legendre_release (statusMessageReporting *smr, tpia_Legendre *Legendre)
int	tpia_Legendre_getFromElement (statusMessageReporting *smr, xData_element *LegendreElement, tpia_Legendre *Legendre)
int	tpia_Legendre_SampleEp (statusMessageReporting *smr, tpia_Legendre *Legendre, int sampleMu, tpia_decaySamplingInfo *decaySamplingInfo)

Function Documentation

int tpia_Legendre_getFromElement	(	statusMessageReporting *	smr,
		xData_element *	LegendreElement,
		tpia_Legendre *	Legendre
	)

Definition at line 78 of file tpia_Legendre.cc.

References tpia_Legendre_s::binned, xData_elementListItem_s::element, tpia_Legendre_s::frame, xData_elementList_s::items, tpia_LegendreBin_s::ls, xData_elementList_s::n, tpia_LegendreBin_s::numberOfLs, status, tpia_frame_setFromElement(), tpia_misc_setMessageError_Element(), xData_addToAccessed(), xData_freeElementList(), xData_getElementsByTagNameAndSort(), and xData_malloc2.

                                                                                                                         {

    int i, status = 0;
    xData_elementList *list;

    xData_addToAccessed( smr, LegendreElement, 1 );
    if( ( tpia_frame_setFromElement( smr, LegendreElement, 4, &(Legendre->frame) ) ) != 0 ) return( 1 );
    list = xData_getElementsByTagNameAndSort( smr, LegendreElement, "l", NULL, NULL );
    if( list->n == 0 ) {
        status = 1;
        tpia_misc_setMessageError_Element( smr, NULL, LegendreElement, __FILE__, __LINE__, 1, "Legendre element does not contain any l elements" ); }
    else {
        //if( ( Legendre->binned.ls = xData_malloc2( smr, list->n * sizeof( tpia_EqualProbableBinSpectra ), 1, "ls" ) ) != NULL ) {
        if( ( Legendre->binned.ls = (tpia_EqualProbableBinSpectra*) xData_malloc2( smr, list->n * sizeof( tpia_EqualProbableBinSpectra ), 1, "ls" ) ) != NULL ) {
            Legendre->binned.numberOfLs = 0;
            for( i = 0; i < list->n; i++ ) {
                if( ( status = _tpia_Legendre_getOrder( smr, list->items[i].element, Legendre, &(Legendre->binned.ls[i]) ) ) != 0 ) break;
                Legendre->binned.numberOfLs++;
            }
        }
    }
    xData_freeElementList( smr, list );
    return( status );
}

int tpia_Legendre_initialize	(	statusMessageReporting *	smr,
		tpia_Legendre *	Legendre
	)

Definition at line 54 of file tpia_Legendre.cc.

References tpia_Legendre_s::binned, tpia_Legendre_s::frame, tpia_LegendreBin_s::ls, tpia_LegendreBin_s::numberOfLs, and tpia_frame_setFromString().

Referenced by tpia_Legendre_release(), and tpia_product_initialize().

                                                                                     {

    memset( Legendre, 0, sizeof( tpia_Legendre ) );
    if( tpia_frame_setFromString( smr, "", "", 0, &(Legendre->frame) ) ) return( 1 );
    Legendre->binned.numberOfLs = 0;
    Legendre->binned.ls = NULL;
    return( 0 );
}

int tpia_Legendre_release	(	statusMessageReporting *	smr,
		tpia_Legendre *	Legendre
	)

Definition at line 65 of file tpia_Legendre.cc.

References tpia_Legendre_s::binned, tpia_EqualProbableBinSpectra_s::energies, tpia_LegendreBin_s::ls, tpia_LegendreBin_s::numberOfLs, tpia_Legendre_initialize(), and xData_free().

Referenced by tpia_product_release().

                                                                                  {

    int i;

    for( i = 0; i < Legendre->binned.numberOfLs; i++ ) xData_free( smr, Legendre->binned.ls[i].energies );
    //Legendre->binned.ls = xData_free( smr, Legendre->binned.ls );
    Legendre->binned.ls = (tpia_EqualProbableBinSpectra*) xData_free( smr, Legendre->binned.ls );
    tpia_Legendre_initialize( smr, Legendre );
    return( 0 );
}

int tpia_Legendre_SampleEp	(	statusMessageReporting *	smr,
		tpia_Legendre *	Legendre,
		int	sampleMu,
		tpia_decaySamplingInfo *	decaySamplingInfo
	)

Definition at line 127 of file tpia_Legendre.cc.

References tpia_Legendre_s::binned, tpia_decaySamplingInfo_s::e_in, tpia_decaySamplingInfo_s::Ep, tpia_LegendreBin_s::ls, tpia_decaySamplingInfo_s::mu, tpia_EqualProbableBinSpectra_s::nBins, tpia_LegendreBin_s::numberOfLs, tpia_decaySamplingInfo_s::rng, tpia_decaySamplingInfo_s::rngState, tpia_misc_drng(), and tpia_misc_sampleEqualProbableBin().

Referenced by tpia_decayChannel_sampleProductsAtE().

                                                                                                                                            {

    tpia_EqualProbableBinSpectra *binned_l0 = Legendre->binned.ls;
    double Ep;

/*
Currently, only l = 0, equal probable binning is supported.
Need to return frame info for Ep, mu, also.
*/
    if( Legendre->binned.numberOfLs > 0 ) {
        if( sampleMu ) decaySamplingInfo->mu = 2. * tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState ) - 1.;
#ifndef WIN32
        if( decaySamplingInfo->mu <= -1 ) decaySamplingInfo->mu = nextafter( -1., 0. );
#endif
#ifdef WIN32
        if( decaySamplingInfo->mu <= -1 ) decaySamplingInfo->mu = _nextafter( -1., 0. );
#endif

        tpia_misc_sampleEqualProbableBin( smr, decaySamplingInfo, decaySamplingInfo->e_in, binned_l0->nBins, binned_l0, &Ep );
        /* ??? Need to check that e_in > E_Threshold */
        decaySamplingInfo->Ep = Ep; }
    else {
        return( 1 );
    }

    return( 0 );
}