tpia_angular.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_angular_initialize (statusMessageReporting *smr, tpia_angular *angular)
int	tpia_angular_release (statusMessageReporting *smr, tpia_angular *angular)
int	tpia_angular_getFromElement (statusMessageReporting *smr, xData_element *angularElement, tpia_angular *angular)
int	tpia_angular_SampleMu (statusMessageReporting *smr, tpia_angular *angular, tpia_decaySamplingInfo *decaySamplingInfo)

Function Documentation

int tpia_angular_getFromElement	(	statusMessageReporting *	smr,
		xData_element *	angularElement,
		tpia_angular *	angular
	)

Definition at line 70 of file tpia_angular.cc.

References tpia_angular_s::binned, tpia_EqualProbableBinSpectra_s::dValue, tpia_EqualProbableBinSpectra_s::energies, tpia_angular_s::frame, int(), tpia_EqualProbableBinSpectra_s::iValue, n, tpia_EqualProbableBinSpectra_s::nBins, tpia_EqualProbableBinSpectra_s::numberOfEs, tpia_frame_setFromElement(), and tpia_misc_getEqualProbableBin().

                                                                                                                     {

    xData_Int n, nBins;

    if( ( tpia_frame_setFromElement( smr, angularElement, 3, &(angular->frame) ) ) != 0 ) return( 1 );
    if( ( angular->binned.energies = tpia_misc_getEqualProbableBin( smr, angularElement, &n, &nBins ) ) != NULL ) {
        angular->binned.iValue = 0;
        angular->binned.dValue = 0;
        angular->binned.nBins = (int) nBins;
        angular->binned.numberOfEs = (int) n;
        return( 0 );
    }
    return( 1 );
}

int tpia_angular_initialize	(	statusMessageReporting *	smr,
		tpia_angular *	angular
	)

Definition at line 49 of file tpia_angular.cc.

References tpia_angular_s::binned, tpia_EqualProbableBinSpectra_s::energies, tpia_angular_s::frame, tpia_EqualProbableBinSpectra_s::numberOfEs, and tpia_frame_setFromString().

Referenced by tpia_product_initialize().

                                                                                  {

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

int tpia_angular_release	(	statusMessageReporting *	smr,
		tpia_angular *	angular
	)

Definition at line 60 of file tpia_angular.cc.

References tpia_angular_s::binned, tpia_EqualProbableBinSpectra_s::energies, tpia_EqualProbableBinSpectra_s::numberOfEs, and xData_free().

Referenced by tpia_product_release().

                                                                               {

    //angular->binned.energies = xData_free( smr, angular->binned.energies );
    angular->binned.energies = (tpia_EqualProbableBinSpectrum*) xData_free( smr, angular->binned.energies );
    angular->binned.numberOfEs = 0;
    return( 0 );
}

int tpia_angular_SampleMu	(	statusMessageReporting *	smr,
		tpia_angular *	angular,
		tpia_decaySamplingInfo *	decaySamplingInfo
	)

Definition at line 87 of file tpia_angular.cc.

References tpia_angular_s::binned, tpia_decaySamplingInfo_s::e_in, tpia_decaySamplingInfo_s::mu, tpia_EqualProbableBinSpectra_s::nBins, tpia_EqualProbableBinSpectra_s::numberOfEs, status, and tpia_misc_sampleEqualProbableBin().

Referenced by G4GIDI_target::getElasticFinalState(), and tpia_decayChannel_sampleProductsAtE().

                                                                                                                           {

    tpia_EqualProbableBinSpectra *binned = &(angular->binned);
    int nBins = binned->nBins, status = 1;
    double mu;
/*
Currently, only equal probable binning is supported.
Need to return frame info for mu also.
*/
    if( binned->numberOfEs > 0 ) {
        tpia_misc_sampleEqualProbableBin( smr, decaySamplingInfo, decaySamplingInfo->e_in, nBins, binned, &mu );
        if( mu < -1. ) mu = -1.;
        if( mu >  1. ) mu =  1.;
        decaySamplingInfo->mu = mu;
        status = 0;
    }

    return( status );
}