tpia_target.h File Reference

#include <xData.h>
#include <tpia_defs.h>
#include <tpia_map.h>
#include <tpi_IDs.h>
#include <tpia_mass.h>

Go to the source code of this file.

Data Structures
struct	tpia_particle_s
struct	tpia_samplingMethods_s
struct	tpia_data_frame_s
struct	tpia_1dData_s
struct	tpia_decaySamplingInfo_s
struct	tpia_productOutgoingData_s
struct	tpia_multiplicity_s
struct	tpia_EqualProbableBinSpectrum_s
struct	tpia_EqualProbableBinSpectra_s
struct	tpia_LegendreBin_s
struct	tpia_angularEnergyBin_s
struct	tpia_decayChannel_s
struct	tpia_angular_s
struct	tpia_Legendre_s
struct	tpia_angularEnergy_s
struct	tpia_product_s
struct	tpia_channel_s
struct	tpia_target_heated_s
struct	tpia_target_heated_info_s
struct	tpia_target_s

Macros
#define	tpia_mode_MonteCarlo   1
#define	tpia_mode_Pn   2
#define	tpia_maxNumberOfFrames   8
#define	tpia_referenceFrame_None   0
#define	tpia_referenceFrame_COM   1
#define	tpia_referenceFrame_lab   2
#define	tpia_referenceFrame_Max   tpia_referenceFrame_lab
#define	tpia_crossSectionType_grouped   1
#define	tpia_crossSectionType_pointwise   2
#define	tpia_m_depositionEnergy   1
#define	tpia_m_multiplicity   ( 1 << 1 )
#define	tpia_m_decayChannel   ( 1 << 2 )
#define	tpia_m_commonShift   3
#define	tpia_m_angular   ( 1 << 3 )
#define	tpia_m_formFactor   ( 1 << 4 )
#define	tpia_m_Legendre   ( 1 << 5 )
#define	tpia_m_angular_energy   ( 1 << 6 )
#define	tpia_speedOfLight_cm_sec   2.99792458e10
#define	tpia_AMU2MeV   931.494028
#define	tpia_samplingMethods_isConstant(method)   ( strcmp ( method , "constant" ) == 0 )
#define	tpia_samplingMethods_isLinear(method)   ( strcmp ( method , "linear" ) == 0 )

Typedefs
typedef struct tpia_particle_s	tpia_particle
typedef struct tpia_target_s	tpia_target
typedef struct tpia_target_heated_info_s	tpia_target_heated_info
typedef struct tpia_target_heated_sorted_s	tpia_target_heated_sorted
typedef struct tpia_target_heated_s	tpia_target_heated
typedef struct tpia_decaySamplingInfo_s	tpia_decaySamplingInfo
typedef struct tpia_productOutgoingData_s	tpia_productOutgoingData
typedef struct tpia_multiplicity_s	tpia_multiplicity
typedef struct tpia_samplingMethods_s	tpia_samplingMethods
typedef struct tpia_EqualProbableBinSpectra_s	tpia_EqualProbableBinSpectra
typedef struct tpia_EqualProbableBinSpectrum_s	tpia_EqualProbableBinSpectrum
typedef struct tpia_LegendreBin_s	tpia_LegendreBin
typedef struct tpia_angularEnergyBin_s	tpia_angularEnergyBin
typedef struct tpia_angular_s	tpia_angular
typedef struct tpia_Legendre_s	tpia_Legendre
typedef struct tpia_angularEnergy_s	tpia_angularEnergy
typedef struct tpia_decayChannel_s	tpia_decayChannel
typedef struct tpia_product_s	tpia_product
typedef struct tpia_1dData_s	tpia_1dData
typedef struct tpia_2dData_s	tpia_2dData
typedef struct tpia_channel_s	tpia_channel
typedef struct tpia_data_frame_s	tpia_data_frame

Functions
tpia_target *	tpia_target_create (statusMessageReporting *smr)
int	tpia_target_initialize (statusMessageReporting *smr, tpia_target *target)
tpia_target *	tpia_target_createRead (statusMessageReporting *smr, const char *fileName)
int	tpia_target_readFromMap (statusMessageReporting *smr, tpia_target *target, tpia_map *map, const char *evaluation, const char *projectileName, const char *targetName)
tpia_target *	tpia_target_createReadFromMap (statusMessageReporting *smr, tpia_map *map, const char *evaluation, const char *projectileName, const char *targetName)
tpia_target *	tpia_target_free (statusMessageReporting *smr, tpia_target *target)
int	tpia_target_release (statusMessageReporting *smr, tpia_target *target)
int	tpia_target_read (statusMessageReporting *smr, tpia_target *target, const char *fileName)
char *	tpia_target_getAttributesValue (statusMessageReporting *smr, tpia_target *target, char const *name)
int	tpia_target_getTemperatures (statusMessageReporting *smr, tpia_target *target, double *temperatures)
int	tpia_target_readHeatedTarget (statusMessageReporting *smr, tpia_target *target, int index, int checkElememtsForAccess)
tpia_target_heated *	tpia_target_getHeatedTargetAtIndex_ReadIfNeeded (statusMessageReporting *smr, tpia_target *target, int index)
int	tpia_target_numberOfChannels (statusMessageReporting *smr, tpia_target *target)
int	tpia_target_numberOfProductionChannels (statusMessageReporting *smr, tpia_target *target)
xData_Int	tpia_target_getEnergyGridAtTIndex (statusMessageReporting *smr, tpia_target *target, int index, double **energyGrid)
tpia_1dData *	tpia_target_getTotalCrossSectionAtTIndex (statusMessageReporting *smr, tpia_target *target, int index, int crossSectionType)
double	tpia_target_getTotalCrossSectionAtTAndE (statusMessageReporting *smr, tpia_target *target, double T, xData_Int iEg, double e_in, int crossSectionType)
double	tpia_target_getIndexChannelCrossSectionAtE (statusMessageReporting *smr, tpia_target *target, int index, double T, xData_Int iEg, double e_in, int crossSectionType)
int	tpia_target_sampleIndexChannelProductsAtE (statusMessageReporting *smr, tpia_target *target, int index, double T, tpia_decaySamplingInfo *decaySamplingInfo, int nProductData, tpia_productOutgoingData *productData)
tpia_target_heated *	tpia_target_heated_create (statusMessageReporting *smr)
int	tpia_target_heated_initialize (statusMessageReporting *smr, tpia_target_heated *target)
tpia_target_heated *	tpia_target_heated_createRead (statusMessageReporting *smr, const char *fileName, int checkElememtsForAccess)
tpia_target_heated *	tpia_target_heated_free (statusMessageReporting *smr, tpia_target_heated *target)
int	tpia_target_heated_release (statusMessageReporting *smr, tpia_target_heated *target)
int	tpia_target_heated_read (statusMessageReporting *smr, tpia_target_heated *target, const char *fileName, int checkElememtsForAccess)
int	tpia_target_heated_numberOfChannels (statusMessageReporting *smr, tpia_target_heated *target)
int	tpia_target_heated_numberOfProductionChannels (statusMessageReporting *smr, tpia_target_heated *target)
tpia_channel *	tpia_target_heated_getChannelAtIndex (tpia_target_heated *target, int index)
tpia_channel *	tpia_target_heated_getChannelAtIndex_smr (statusMessageReporting *smr, tpia_target_heated *target, int index)
tpia_channel *	tpia_target_heated_getProductionChannelAtIndex (tpia_target_heated *target, int index)
xData_Int	tpia_target_heated_getEnergyGrid (statusMessageReporting *smr, tpia_target_heated *target, double **energyGrid)
xData_Int	tpia_target_heated_getEIndex (tpia_target_heated *target, double e_in)
double	tpia_target_heated_getTotalCrossSectionAtE (statusMessageReporting *smr, tpia_target_heated *target, xData_Int gE, double e_in, int crossSectionType)
double	tpia_target_heated_getIndexChannelCrossSectionAtE (statusMessageReporting *smr, tpia_target_heated *target, int index, xData_Int iEg, double e_in, int crossSectionType)
int	tpia_target_heated_sampleIndexChannelProductsAtE (statusMessageReporting *smr, tpia_target_heated *target, int index, tpia_decaySamplingInfo *decaySamplingInfo, int nProductData, tpia_productOutgoingData *productData)
tpia_channel *	tpia_channel_create (statusMessageReporting *smr)
int	tpia_channel_initialize (statusMessageReporting *smr, 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 *smr, tpia_channel *channel, int index)
int	tpia_channel_numberOfProducts (statusMessageReporting *smr, tpia_channel *channel)
int	tpia_channel_isProduction (statusMessageReporting *smr, tpia_channel *channel)
double	tpia_channel_getCrossSectionAtE (statusMessageReporting *smr, tpia_channel *channel, xData_Int iEg, double e_in, int crossSectionType)
tpia_particle *	tpia_particle_create (statusMessageReporting *smr)
int	tpia_particle_initialize (statusMessageReporting *smr, tpia_particle *particle)
tpia_particle *	tpia_particle_free (statusMessageReporting *smr, tpia_particle *particle)
int	tpia_particle_release (statusMessageReporting *smr, tpia_particle *particle)
int	tpia_particle_freeInternalList (statusMessageReporting *smr)
tpia_particle *	tpia_particle_getInternalID (statusMessageReporting *smr, const char *const name)
int	tpia_particle_printInternalSortedList (statusMessageReporting *smr)
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 *smr, tpia_product *product)
tpia_product *	tpia_product_getFirstProduct (tpia_product *product)
tpia_product *	tpia_product_getProductByIndex (statusMessageReporting *smr, tpia_product *product, int index)
int	tpia_product_doesDecay (statusMessageReporting *smr, tpia_product *product)
int	tpia_product_numberOfProducts (statusMessageReporting *smr, tpia_product *product)
int	tpia_product_isDataPresent (statusMessageReporting *smr, tpia_product *product, int b_data)
int	tpia_product_sampleMultiplicity (statusMessageReporting *smr, tpia_product *product, double e_in, double r)
tpia_product *	tpia_decayChannel_getFirstProduct (tpia_decayChannel *decayChannel)
tpia_product *	tpia_decayChannel_getNextProduct (tpia_product *product)
int	tpia_decayChannel_sampleProductsAtE (statusMessageReporting *smr, tpia_decayChannel *decayChannel, tpia_decaySamplingInfo *decaySamplingInfo, int nProductData, tpia_productOutgoingData *productData)
int	tpia_kinetics_2BodyReaction (statusMessageReporting *smr, tpia_decayChannel *decayChannel, double K, double mu, double phi, tpia_productOutgoingData *outgoingData)
int	tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum (statusMessageReporting *smr, double beta, double e_kinetic_com, double mu, double phi, double m3cc, double m4cc, tpia_productOutgoingData *outgoingData)
int	tpia_frame_clear (statusMessageReporting *smr, tpia_data_frame *frame)
int	tpia_frame_setFromElement (statusMessageReporting *smr, xData_element *element, int dimension, tpia_data_frame *frame)
int	tpia_frame_setFromString (statusMessageReporting *smr, const char *forItem, const char *value, int dimension, tpia_data_frame *frame)
int	tpia_frame_getDimensions (statusMessageReporting *smr, tpia_data_frame *frame)
char *	tpia_frame_toString (statusMessageReporting *smr, const char *fromItem, tpia_data_frame *frame)
int	tpia_frame_setColumns (statusMessageReporting *smr, tpia_data_frame *frame, int nColumns, int *values)
int	tpia_frame_setColumn (statusMessageReporting *smr, tpia_data_frame *frame, int column, int value)
int	tpia_frame_getColumn (statusMessageReporting *smr, tpia_data_frame *frame, int column)
tpia_multiplicity *	tpia_multiplicity_create (statusMessageReporting *smr)
int	tpia_multiplicity_initialize (statusMessageReporting *smr, tpia_multiplicity *multiplicity)
tpia_multiplicity *	tpia_multiplicity_free (statusMessageReporting *smr, tpia_multiplicity *multiplicity)
int	tpia_multiplicity_release (statusMessageReporting *smr, tpia_multiplicity *multiplicity)
tpia_multiplicity *	tpia_multiplicity_createGetFromElement (statusMessageReporting *smr, xData_element *multiplicityElement, int nGroups)
int	tpia_multiplicity_getFromElement (statusMessageReporting *smr, xData_element *multiplicityElement, tpia_multiplicity *multiplicity, int nGroups)
int	tpia_multiplicity_getTimeScaleFromElement (statusMessageReporting *smr, xData_element *element, const char **timeScale, int *isDelayedNeutrons, double *dTimeScale)
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)
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)
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)
int	tpia_samplingMethods_initialize (statusMessageReporting *smr, tpia_samplingMethods *samplingMethods)
int	tpia_misc_NumberOfZSymbols (void)
const char *	tpia_misc_ZToSymbol (int iZ)
int	tpia_misc_symbolToZ (const char *Z)
int	tpia_miscNameToZAm (statusMessageReporting *smr, const char *name, int *Z, int *A, int *m)
char *	tpia_misc_pointerToAttributeIfAllOk (statusMessageReporting *smr, xData_element *element, const char *path, int required, xData_attributionList *attributes, const char *name, const char *file, int line)
double *	tpia_misc_get2dx_y_data (statusMessageReporting *smr, xData_element *element, xData_Int *length)
double *	tpia_misc_get2dxindex_y_data (statusMessageReporting *smr, xData_element *element, xData_Int *start, xData_Int *end, double *xValues)
double *	tpia_misc_get2d_xShared_yHistogram_data (statusMessageReporting *smr, xData_element *element, xData_Int *start, xData_Int *end, xData_Int *length)
int	tpia_misc_get2d_xShared_yHistogram_data_Grouped (statusMessageReporting *smr, xData_element *element, tpia_1dData *group)
double	tpia_misc_getPointwiseCrossSectionAtE (statusMessageReporting *smr, tpia_1dData *crossSection, double *energyGrid, xData_Int index, double e_in)
double	tpia_misc_drng (double(*rng)(void *), void *rngState)
int	tpia_misc_sampleEqualProbableBin (statusMessageReporting *smr, tpia_decaySamplingInfo *decaySamplingInfo, double e_in, int nBins, tpia_EqualProbableBinSpectra *binned, double *value_)

Variables
DLL_LEND const char *	tpia_productGenre_unknown
DLL_LEND const char *	tpia_productGenre_twoBody_angular
DLL_LEND const char *	tpia_productGenre_twoBody_formFactor
DLL_LEND const char *	tpia_productGenre_NBody_Legendre
DLL_LEND const char *	tpia_productGenre_NBody_angular_energy
DLL_LEND const char *	tpia_productGenre_NBody_uncorrelate_Legendre
DLL_LEND const char *	tpia_productGenre_NBody_pairProduction
DLL_LEND const char *	tpia_samplingMethods_constant
DLL_LEND const char *	tpia_samplingMethods_linear

Macro Definition Documentation

#define tpia_AMU2MeV   931.494028

Definition at line 101 of file tpia_target.h.

Referenced by tpia_particle_getInternalID().

#define tpia_crossSectionType_grouped   1

Definition at line 88 of file tpia_target.h.

Referenced by tpia_channel_getCrossSectionAtE(), tpia_target_getTotalCrossSectionAtTIndex(), and tpia_target_heated_getTotalCrossSectionAtE().

#define tpia_crossSectionType_pointwise   2

Definition at line 89 of file tpia_target.h.

Referenced by G4GIDI_target::getTotalCrossSectionAtE(), G4GIDI_target::sampleChannelCrossSectionAtE(), G4GIDI_target::sumChannelCrossSectionAtE(), tpia_channel_getCrossSectionAtE(), tpia_target_getTotalCrossSectionAtTIndex(), and tpia_target_heated_getTotalCrossSectionAtE().

#define tpia_m_angular   ( 1 << 3 )

Definition at line 95 of file tpia_target.h.

#define tpia_m_angular_energy   ( 1 << 6 )

Definition at line 98 of file tpia_target.h.

#define tpia_m_commonShift   3

Definition at line 94 of file tpia_target.h.

#define tpia_m_decayChannel   ( 1 << 2 )

Definition at line 93 of file tpia_target.h.

#define tpia_m_depositionEnergy   1

Definition at line 91 of file tpia_target.h.

#define tpia_m_formFactor   ( 1 << 4 )

Definition at line 96 of file tpia_target.h.

#define tpia_m_Legendre   ( 1 << 5 )

Definition at line 97 of file tpia_target.h.

#define tpia_m_multiplicity   ( 1 << 1 )

Definition at line 92 of file tpia_target.h.

#define tpia_maxNumberOfFrames   8

Definition at line 82 of file tpia_target.h.

Referenced by tpia_frame_getColumn(), tpia_frame_getDimensions(), tpia_frame_setColumn(), and tpia_frame_setFromString().

#define tpia_mode_MonteCarlo   1

Definition at line 80 of file tpia_target.h.

#define tpia_mode_Pn   2

Definition at line 81 of file tpia_target.h.

#define tpia_referenceFrame_COM   1

Definition at line 84 of file tpia_target.h.

Referenced by tpia_frame_setFromString(), and tpia_frame_toString().

#define tpia_referenceFrame_lab   2

Definition at line 85 of file tpia_target.h.

Referenced by G4GIDI_target::getElasticFinalState(), G4GIDI_target::getFinalState(), tpia_decayChannel_sampleProductsAtE(), tpia_frame_setFromString(), tpia_frame_toString(), and tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum().

#define tpia_referenceFrame_Max   tpia_referenceFrame_lab

Definition at line 86 of file tpia_target.h.

Referenced by tpia_frame_setColumn().

#define tpia_referenceFrame_None   0

Definition at line 83 of file tpia_target.h.

Referenced by tpia_frame_setColumn(), and tpia_frame_setFromString().

#define tpia_samplingMethods_isConstant

(

method

)

( strcmp ( method , "constant" ) == 0 )

Definition at line 120 of file tpia_target.h.

#define tpia_samplingMethods_isLinear

(

method

)

( strcmp ( method , "linear" ) == 0 )

Definition at line 121 of file tpia_target.h.

Referenced by tpia_misc_sampleEqualProbableBin().

#define tpia_speedOfLight_cm_sec   2.99792458e10

Definition at line 100 of file tpia_target.h.

Referenced by tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum().

Typedef Documentation

typedef struct tpia_1dData_s tpia_1dData

Definition at line 70 of file tpia_target.h.

typedef struct tpia_2dData_s tpia_2dData

Definition at line 71 of file tpia_target.h.

typedef struct tpia_angular_s tpia_angular

Definition at line 65 of file tpia_target.h.

typedef struct tpia_angularEnergy_s tpia_angularEnergy

Definition at line 67 of file tpia_target.h.

typedef struct tpia_angularEnergyBin_s tpia_angularEnergyBin

Definition at line 64 of file tpia_target.h.

typedef struct tpia_channel_s tpia_channel

Definition at line 72 of file tpia_target.h.

typedef struct tpia_data_frame_s tpia_data_frame

Definition at line 73 of file tpia_target.h.

typedef struct tpia_decayChannel_s tpia_decayChannel

Definition at line 68 of file tpia_target.h.

typedef struct tpia_decaySamplingInfo_s tpia_decaySamplingInfo

Definition at line 57 of file tpia_target.h.

typedef struct tpia_EqualProbableBinSpectra_s tpia_EqualProbableBinSpectra

Definition at line 61 of file tpia_target.h.

typedef struct tpia_EqualProbableBinSpectrum_s tpia_EqualProbableBinSpectrum

Definition at line 62 of file tpia_target.h.

typedef struct tpia_Legendre_s tpia_Legendre

Definition at line 66 of file tpia_target.h.

typedef struct tpia_LegendreBin_s tpia_LegendreBin

Definition at line 63 of file tpia_target.h.

typedef struct tpia_multiplicity_s tpia_multiplicity

Definition at line 59 of file tpia_target.h.

typedef struct tpia_particle_s tpia_particle

Definition at line 52 of file tpia_target.h.

typedef struct tpia_product_s tpia_product

Definition at line 69 of file tpia_target.h.

typedef struct tpia_productOutgoingData_s tpia_productOutgoingData

Definition at line 58 of file tpia_target.h.

typedef struct tpia_samplingMethods_s tpia_samplingMethods

Definition at line 60 of file tpia_target.h.

typedef struct tpia_target_s tpia_target

Definition at line 53 of file tpia_target.h.

typedef struct tpia_target_heated_s tpia_target_heated

Definition at line 56 of file tpia_target.h.

typedef struct tpia_target_heated_info_s tpia_target_heated_info

Definition at line 54 of file tpia_target.h.

typedef struct tpia_target_heated_sorted_s tpia_target_heated_sorted

Definition at line 55 of file tpia_target.h.

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 );
}

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 );
}

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	iEg,
		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 *	smr,
		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 *	smr,
		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 *	smr,
		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 *	smr,
		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 );
}

tpia_product* tpia_decayChannel_getFirstProduct

(

tpia_decayChannel *

decayChannel

)

Definition at line 51 of file tpia_decayChannel.cc.

References tpia_decayChannel_s::products.

Referenced by G4GIDI_target::getElasticFinalState(), tpia_channel_getFirstProduct(), tpia_decayChannel_sampleProductsAtE(), tpia_kinetics_2BodyReaction(), and tpia_product_getFirstProduct().

                                                                                   {

    return( decayChannel->products );
}

tpia_product* tpia_decayChannel_getNextProduct

(

tpia_product *

product

)

Definition at line 58 of file tpia_decayChannel.cc.

References tpia_product_s::next.

Referenced by G4GIDI_Misc_channelProductsCompare(), tpia_channel_getProductByIndex(), tpia_decayChannel_sampleProductsAtE(), tpia_kinetics_2BodyReaction(), and tpia_product_getProductByIndex().

                                                                        {

    return( product->next );
}

int tpia_decayChannel_sampleProductsAtE	(	statusMessageReporting *	smr,
		tpia_decayChannel *	decayChannel,
		tpia_decaySamplingInfo *	decaySamplingInfo,
		int	nProductData,
		tpia_productOutgoingData *	productData
	)

Definition at line 65 of file tpia_decayChannel.cc.

References tpia_productOutgoingData_s::decayChannel, tpia_decaySamplingInfo_s::e_in, tpia_decaySamplingInfo_s::Ep, tpia_decaySamplingInfo_s::frame, tpia_productOutgoingData_s::frame, tpia_decaySamplingInfo_s::genre, tpia_productOutgoingData_s::genre, tpia_product_s::genre, tpia_decaySamplingInfo_s::isVelocity, tpia_productOutgoingData_s::isVelocity, tpia_productOutgoingData_s::kineticEnergy, M_PI, tpia_decaySamplingInfo_s::mu, n, tpia_decayChannel_s::numberOfProducts, tpia_decaySamplingInfo_s::productID, tpia_productOutgoingData_s::productID, tpia_product_s::productID, tpia_productOutgoingData_s::px_vx, tpia_productOutgoingData_s::py_vy, tpia_productOutgoingData_s::pz_vz, tpia_decaySamplingInfo_s::rng, tpia_decaySamplingInfo_s::rngState, smr_isOk(), smr_setMessageError(), tpia_angular_SampleMu(), tpia_angularEnergy_SampleEp(), tpia_decayChannel_getFirstProduct(), tpia_decayChannel_getNextProduct(), tpia_frame_setColumns(), tpia_kinetics_2BodyReaction(), tpia_Legendre_SampleEp(), tpia_misc_drng(), tpia_product_sampleMultiplicity(), and tpia_referenceFrame_lab.

Referenced by tpia_target_heated_sampleIndexChannelProductsAtE().

                                                                   {

    int i, n = 0, multiplicity, secondTwoBody = 0, labFrame = tpia_referenceFrame_lab;
    tpia_product *product, *nextProduct;
    double phi, p;

    if( nProductData < decayChannel->numberOfProducts ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "nProductData = %d < decayChannel->numberOfProducts = %d", nProductData,
            decayChannel->numberOfProducts );
        return( -1 );
    }
    for( i = 0, product = tpia_decayChannel_getFirstProduct( decayChannel ); product != NULL; i++, product = tpia_decayChannel_getNextProduct( product ) ) {
        nextProduct = tpia_decayChannel_getNextProduct( product );
        if( !secondTwoBody ) {
            if( ( multiplicity = product->multiplicity ) == 0 ) multiplicity = tpia_product_sampleMultiplicity( smr, product, decaySamplingInfo->e_in,
                tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState ) );
            while( multiplicity > 0 ) {
                if( n >= nProductData ) break;          /* This needs work. */
                multiplicity--;
                decaySamplingInfo->genre = product->genre;
                decaySamplingInfo->productID = product->productID;
                decaySamplingInfo->mu = 0;
                decaySamplingInfo->Ep = 0;
                productDatas[n].genre = product->genre;
                productDatas[n].isVelocity = decaySamplingInfo->isVelocity;
                tpia_frame_setColumns( smr, &(productDatas[n].frame), 1, &labFrame );
                productDatas[n].productID = product->productID;
                productDatas[n].decayChannel = &(product->decayChannel);
                if( strcmp( product->genre, "twoBody_angular" ) == 0 ) {
                    secondTwoBody = 1;
                    productDatas[n+1].productID = nextProduct->productID;
                    productDatas[n].genre = product->genre;
                    tpia_angular_SampleMu( smr, &(product->angular), decaySamplingInfo );   /* Need to test for success. */
                    if( smr_isOk( smr ) ) {
                        phi = 2. * M_PI * tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState );
                        productDatas[n].isVelocity = decaySamplingInfo->isVelocity;
                        productDatas[n].frame = decaySamplingInfo->frame;
                        tpia_kinetics_2BodyReaction( smr, decayChannel, decaySamplingInfo->e_in, decaySamplingInfo->mu, phi, &(productDatas[n]) );
                    } }
                else if( strcmp( product->genre, "NBody_Legendre" ) == 0 ) {
                    tpia_Legendre_SampleEp( smr, &(product->Legendre), 1, decaySamplingInfo ); }
                else if( strcmp( product->genre, "NBody_angular_energy" ) == 0 ) {
                    tpia_angular_SampleMu( smr, &(product->angular), decaySamplingInfo );   /* Need to test for success. */
                    tpia_angularEnergy_SampleEp( smr, &(product->angularEnergy), decaySamplingInfo ); }
                else if( strcmp( product->genre, "NBody_uncorrelate_Legendre" ) == 0 ) {
                    tpia_angular_SampleMu( smr, &(product->angular), decaySamplingInfo );   /* Need to test for success. */
                    tpia_Legendre_SampleEp( smr, &(product->Legendre), 0, decaySamplingInfo ); }
                else if( strcmp( product->genre, "unknown" ) == 0 ) {
                    }
                else {
                    printf( "Unknown spectral data form product name = %s, genre = %s\n", product->productID->name, product->genre );
                }
                if( !smr_isOk( smr ) ) return( -1 );
                if( secondTwoBody ) {
                    n++;
                    productDatas[n].productID = nextProduct->productID;
                    productDatas[n].genre = nextProduct->genre; }
                else {
                    productDatas[n].kineticEnergy = decaySamplingInfo->Ep;
                    p = std::sqrt( decaySamplingInfo->Ep * ( decaySamplingInfo->Ep + 2. * product->productID->fullMass_MeV ) );
                    productDatas[n].pz_vz = p * decaySamplingInfo->mu;
                    p = std::sqrt( 1. - decaySamplingInfo->mu * decaySamplingInfo->mu ) * p;
                    phi = 2. * M_PI * tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState );
                    productDatas[n].px_vx = p * std::sin( phi );
                    productDatas[n].py_vy = p * std::cos( phi );
                }
                n++;
            }
        }
    }
    return( n );
}

int tpia_frame_clear	(	statusMessageReporting *	smr,
		tpia_data_frame *	frame
	)

Definition at line 54 of file tpia_frame.cc.

References tpia_data_frame_s::frames.

Referenced by tpia_frame_setColumns().

                                                                         {

    frame->frames = 0;
    return( 0 );
}

int tpia_frame_getColumn	(	statusMessageReporting *	smr,
		tpia_data_frame *	frame,
		int	column
	)

Definition at line 214 of file tpia_frame.cc.

References tpia_data_frame_s::frames, smr_setMessageError(), and tpia_maxNumberOfFrames.

Referenced by tpia_frame_toString(), and tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum().

                                                                                            {

    unsigned int value;

    if( ( column < 0 ) || ( column >= tpia_maxNumberOfFrames ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad column = %d value for getting frame (0 <= column < %d)", column, tpia_maxNumberOfFrames );
        return( -1 ); }
    else {
        value = frame->frames >> ( 2 * column );
    }
    return( (int) ( value & 3 ) );
}

int tpia_frame_getDimensions	(	statusMessageReporting *	smr,
		tpia_data_frame *	frame
	)

Definition at line 122 of file tpia_frame.cc.

References tpia_data_frame_s::frames, and tpia_maxNumberOfFrames.

Referenced by tpia_frame_toString().

                                                                                 {

    int i, dimension = 0;
    unsigned value = frame->frames;

    for( i = 0; i < tpia_maxNumberOfFrames; i++ ) {
        if( ( value & 3 ) == 0 ) break;
        dimension++;
        value = value >> 2;
    }
    return( dimension );
}

int tpia_frame_setColumn	(	statusMessageReporting *	smr,
		tpia_data_frame *	frame,
		int	column,
		int	value
	)

Definition at line 191 of file tpia_frame.cc.

References tpia_data_frame_s::frames, smr_setMessageError(), tpia_maxNumberOfFrames, tpia_referenceFrame_Max, and tpia_referenceFrame_None.

Referenced by G4GIDI_target::getElasticFinalState(), G4GIDI_target::getFinalState(), tpia_frame_setColumns(), and tpia_frame_setFromString().

                                                                                                       {

    int i;

    if( ( column < 0 ) || ( column >= tpia_maxNumberOfFrames ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad column = %d value for setting frame (0 <= column < %d)", column, tpia_maxNumberOfFrames );
        return( 1 ); }
    else {
        if( ( value < tpia_referenceFrame_None ) || ( value > tpia_referenceFrame_Max ) ) {
            smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad frame value = %d for column = %d", value, column );
            return( 1 ); }
        else {
            i = 3 << ( 2 * column );
            frame->frames = frame->frames & (~i);
            value = value << ( 2 * column );
            frame->frames = frame->frames | value;
        }
    }
    return( 0 );
}

int tpia_frame_setColumns	(	statusMessageReporting *	smr,
		tpia_data_frame *	frame,
		int	nColumns,
		int *	values
	)

Definition at line 179 of file tpia_frame.cc.

References tpia_frame_clear(), and tpia_frame_setColumn().

Referenced by tpia_decayChannel_sampleProductsAtE().

                                                                                                            {

    int i;

    tpia_frame_clear( smr, frame );
    for( i = 0; i < nColumns; i++ ) if( tpia_frame_setColumn( smr, frame, i, values[i] ) != 0 ) return( 1 );

    return( 0 );
}

int tpia_frame_setFromElement	(	statusMessageReporting *	smr,
		xData_element *	element,
		int	dimension,
		tpia_data_frame *	frame
	)

Definition at line 62 of file tpia_frame.cc.

References xData_element_s::fullName, tpia_frame_setFromString(), tpia_misc_setMessageError_Element(), and xData_getAttributesValueInElement().

Referenced by tpia_angular_getFromElement(), tpia_angularEnergy_getFromElement(), tpia_channel_getFromElement(), tpia_Legendre_getFromElement(), and tpia_multiplicity_getFromElement().

                                                                                                                            {

    char const *value = xData_getAttributesValueInElement( element, "frame" );

    if( value == NULL ) {
        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "element is missing frame attribute" );
        return( 1 );
    }
    return( tpia_frame_setFromString( smr, element->fullName, value, dimension, frame ) );
}

int tpia_frame_setFromString	(	statusMessageReporting *	smr,
		const char *	forItem,
		const char *	value,
		int	dimension,
		tpia_data_frame *	frame
	)

Definition at line 75 of file tpia_frame.cc.

References smr_isOk(), smr_setMessageError(), status, tpia_frame_setColumn(), tpia_maxNumberOfFrames, tpia_referenceFrame_COM, tpia_referenceFrame_lab, and tpia_referenceFrame_None.

Referenced by tpia_angular_initialize(), tpia_angularEnergy_initialize(), tpia_frame_setFromElement(), tpia_Legendre_initialize(), and tpia_multiplicity_initialize().

                                                                                                                                           {

    char const *e;
    int status = 1, i;

    for( i = 0; i < tpia_maxNumberOfFrames; i++ ) tpia_frame_setColumn( smr, frame, i, tpia_referenceFrame_None );
    if( dimension > tpia_maxNumberOfFrames ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "dimension argument = %d is greater than tpia_maxNumberOfFrames = %d", dimension,
            tpia_maxNumberOfFrames );
        return( status );
    }
    for( i = 0, e = value; ( i < dimension ) && ( *e != 0 ); i++ ) {
        if( strstr( e, lab_Str ) == e ) {
            tpia_frame_setColumn( smr, frame, i, tpia_referenceFrame_lab );
            e += 3; }
        else if( strstr( e, COM_Str ) == e ) {
            tpia_frame_setColumn( smr, frame, i, tpia_referenceFrame_COM );
            e += 12; }
        else {
            smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad frame '%s' for %s", value, forItem );
            break;
        }
        if( *e != 0 ) {
            if( *e != ',' ) {
                smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad separater for frame '%s' for %s", value, forItem );
                break;
            }
            e++;
        }
    }
    if( smr_isOk( smr ) ) {
        if( i == dimension ) {
            if( *e == 0 ) {
                status = 0; }
            else {
                smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "extra values for frame '%s' for %s", value, forItem );
            } }
        else {
            smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "missing values for frame '%s' for %s", value, forItem );
        }
    }
    return( status );
}

char* tpia_frame_toString	(	statusMessageReporting *	smr,
		const char *	fromItem,
		tpia_data_frame *	frame
	)

Definition at line 138 of file tpia_frame.cc.

References n, smr_setMessageError(), tpia_frame_getColumn(), tpia_frame_getDimensions(), tpia_referenceFrame_COM, tpia_referenceFrame_lab, and xData_malloc2.

                                                                                               {

    int i, n = tpia_frame_getDimensions( smr, frame ), value, nStr = 0;
    char *str = NULL, *p;

    if( nlab_Str == 0 ) {
        nlab_Str = strlen( lab_Str );
        nCOM_Str = strlen( COM_Str );
    }
    for( i = 0; i < n; i++ ) {
        value = tpia_frame_getColumn( smr, frame, i );
        if( value == tpia_referenceFrame_COM ) {
            nStr += nCOM_Str + 1; }
        else if( value == tpia_referenceFrame_lab ) {
            nStr += nlab_Str + 1; }
        else {
            smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad frame value = %d for column = %d", value, i );
            return( NULL );
        }
    }
    if( nStr == 0 ) nStr = 1;
    //if( ( str = xData_malloc2( smr, nStr, 1, "str" ) ) == NULL ) return( NULL );
    if( ( str = (char*) xData_malloc2( smr, nStr, 1, "str" ) ) == NULL ) return( NULL );
    for( i = 0, p = str - 1; i < n; i++ ) {
        p++;
        value = tpia_frame_getColumn( smr, frame, i );
        if( value == tpia_referenceFrame_COM ) {
            strcpy( p, COM_Str );
            p += nCOM_Str; }
        else {
            strcpy( p, lab_Str );
            p += nlab_Str;
        }
        *p = ',';
    }
    *p = 0;
    return( str );
}

int tpia_kinetics_2BodyReaction	(	statusMessageReporting *	smr,
		tpia_decayChannel *	decayChannel,
		double	K,
		double	mu,
		double	phi,
		tpia_productOutgoingData *	outgoingData
	)

Definition at line 48 of file tpia_kinetics.cc.

References tpia_productOutgoingData_s::decayChannel, tpia_product_s::decayChannel, tpia_particle_s::fullMass_MeV, tpia_productOutgoingData_s::genre, tpia_decayChannel_s::m1_fullMass_MeV, python.hepunit::m2, tpia_decayChannel_s::m2_fullMass_MeV, python.hepunit::m3, tpia_productOutgoingData_s::productID, tpia_product_s::productID, tpia_decayChannel_getFirstProduct(), tpia_decayChannel_getNextProduct(), tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum(), and test::x.

Referenced by tpia_decayChannel_sampleProductsAtE().

                                                 {

    tpia_product *pp3 = tpia_decayChannel_getFirstProduct( decayChannel ), *pp4;
    double m1 = decayChannel->m1_fullMass_MeV, m2 = decayChannel->m2_fullMass_MeV, m3, m4, mi, mf, Kp, x, beta;

    pp4 = tpia_decayChannel_getNextProduct( pp3 );
    m3 = pp3->productID->fullMass_MeV;
    m4 = pp4->productID->fullMass_MeV;
    mi = m1 + m2;
    mf = m3 + m4;
    beta = std::sqrt( K * ( K + 2. * m1 ) ) / ( K + mi );
    x = K * m2 / ( mi * mi );
    if( x < 2e-5 ) {                                        /* Kp is the total kinetic energy for m3 and m4 in the COM frame. */
        Kp = mi - mf + K * m2 / mi * ( 1 - 0.5 * x * ( 1 - x ) ); }
    else {
        Kp = std::sqrt( mi * mi + 2 * K * m2 ) - mf;
    }
    if( Kp < 0 ) Kp = 0.;           /* ???? There needs to be a better test here. */
    outgoingData[0].decayChannel = &(pp3->decayChannel);
    outgoingData[1].genre = outgoingData[0].genre;
    outgoingData[1].productID = pp4->productID;
    outgoingData[1].decayChannel = &(pp4->decayChannel);
    return( tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum( smr, beta, Kp, mu, phi, m3, m4, outgoingData ) );
}

int tpia_kinetics_COMKineticEnergy2LabEnergyAndMomentum	(	statusMessageReporting *	smr,
		double	beta,
		double	e_kinetic_com,
		double	mu,
		double	phi,
		double	m3cc,
		double	m4cc,
		tpia_productOutgoingData *	outgoingData
	)

Definition at line 77 of file tpia_kinetics.cc.

References tpia_productOutgoingData_s::frame, tpia_productOutgoingData_s::isVelocity, tpia_productOutgoingData_s::kineticEnergy, tpia_productOutgoingData_s::px_vx, tpia_productOutgoingData_s::py_vy, tpia_productOutgoingData_s::pz_vz, tpia_frame_getColumn(), tpia_referenceFrame_lab, tpia_speedOfLight_cm_sec, and test::x.

Referenced by tpia_kinetics_2BodyReaction().

                                                                           {
/*
*   beta            the velocity/speedOflight of the com frame relative to the lab frame.
*   e_kinetic_com   Total kinetic energy (K1 + K2) in the COM frame.
*   mu              std::cos( theta ) in the COM frame.
*/
    double x, v_p, p, pp3, pp4, px3, py3, pz3, pz4, pz, p_perp2, E3, E4, gamma, m3cc2 = m3cc * m3cc, m4cc2 = m4cc * m4cc;

    p = std::sqrt( e_kinetic_com * ( e_kinetic_com + 2. * m3cc ) * ( e_kinetic_com + 2. * m4cc )  * ( e_kinetic_com + 2. * ( m3cc + m4cc ) ) ) /
            ( 2. * ( e_kinetic_com + m3cc + m4cc ) );
    py3 = p * std::sqrt( 1 - mu * mu );
    px3 = py3 * std::cos( phi );
    py3 *= std::sin( phi );
    pz = p * mu;
    if( tpia_frame_getColumn( NULL, &(outgoingData[0].frame), 0 ) == tpia_referenceFrame_lab ) {
        E3 = std::sqrt( p * p + m3cc2 );
        E4 = std::sqrt( p * p + m4cc2 );
        gamma = std::sqrt( 1. / ( 1. - beta * beta ) );
        pz3 = gamma * (  pz + beta * E3 );
        pz4 = gamma * ( -pz + beta * E4 ); }
    else {
        pz3 = pz;
        pz4 = -pz;
    }
    outgoingData[1].isVelocity = outgoingData[0].isVelocity;
    outgoingData[1].frame = outgoingData[0].frame;

    p_perp2 = px3 * px3 + py3 * py3;

    outgoingData[0].px_vx = px3;
    outgoingData[0].py_vy = py3;
    outgoingData[0].pz_vz = pz3;
    pp3 = p_perp2 + pz3 * pz3;
    x = pp3 / ( 2 * m3cc2 );
    if( x < 1e-5 ) {
        outgoingData[0].kineticEnergy = m3cc * x  * ( 1 - 0.5 * x * ( 1 - x ) ); }
    else {
        outgoingData[0].kineticEnergy = std::sqrt( m3cc2 + pp3 ) - m3cc;
    }
    outgoingData[1].px_vx = -px3;
    outgoingData[1].py_vy = -py3;
    outgoingData[1].pz_vz = pz4;
    pp4 = p_perp2 + pz4 * pz4;
    x = pp4 / ( 2 * m4cc2 );
    if( x < 1e-5 ) {
        outgoingData[1].kineticEnergy = m4cc * x  * ( 1 - 0.5 * x * ( 1 - x ) ); }
    else {
        outgoingData[1].kineticEnergy = std::sqrt( m4cc2 + pp4 ) - m4cc;
    }

    if( outgoingData[0].isVelocity ) {
        v_p = tpia_speedOfLight_cm_sec / std::sqrt( pp3 + m3cc2 );
        outgoingData[0].px_vx *= v_p;
        outgoingData[0].py_vy *= v_p;
        outgoingData[0].pz_vz *= v_p;

        v_p = tpia_speedOfLight_cm_sec / std::sqrt( pp4 + m4cc2 );
        outgoingData[1].px_vx *= v_p;
        outgoingData[1].py_vy *= v_p;
        outgoingData[1].pz_vz *= v_p;
    }

    return( 0 );
}

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 );
}

double tpia_misc_drng	(	double(*)(void *)	rng,
		void *	rngState
	)

Definition at line 403 of file tpia_misc.cc.

References CLHEP::HepRandomEngine::flat(), and CLHEP::HepRandom::getTheEngine().

Referenced by G4GIDI_target::sampleChannelCrossSectionAtE(), tpia_decayChannel_sampleProductsAtE(), tpia_Legendre_SampleEp(), and tpia_misc_sampleEqualProbableBin().

                                                                 {

    double r;

    if( rng != NULL ) {
        r = rng( rngState ); }
    else {
        //r = drand48( );
        r = CLHEP::HepRandom::getTheEngine()->flat();
          
    }
    return( r );
}

double* tpia_misc_get2d_xShared_yHistogram_data	(	statusMessageReporting *	smr,
		xData_element *	element,
		xData_Int *	start,
		xData_Int *	end,
		xData_Int *	length
	)

Definition at line 268 of file tpia_misc.cc.

References xDataType_s::end, xDataType_s::length, xDataType_s::start, xData_2d_xShared_yHistogram_copyData(), xData_addToAccessed(), xData_getOneElementByTagName(), and xData_element_s::xDataTypeInfo.

Referenced by tpia_misc_get2d_xShared_yHistogram_data_Grouped().

                                                                                                                                                            {

    xData_Int i;
    xData_element *xDataElement;
    double *d = NULL;

    xData_addToAccessed( smr, element, 1 );
    //if( ( xDataElement = xData_getOneElementByTagName( smr, element, "xData", 1 ) ) != NULL ) {
    if( ( xDataElement = xData_getOneElementByTagName( smr, element, (char*) "xData", 1 ) ) != NULL ) {
        xData_addToAccessed( smr, xDataElement, 1 );
        if( ( d = xData_2d_xShared_yHistogram_copyData( smr, xDataElement, &i ) ) != NULL ) {
            if( start != NULL ) *start = xDataElement->xDataTypeInfo.start;
            if( end != NULL ) *end = xDataElement->xDataTypeInfo.end;
            if( length != NULL ) *length = xDataElement->xDataTypeInfo.length;
        }
    }
    return( d );
}

int tpia_misc_get2d_xShared_yHistogram_data_Grouped	(	statusMessageReporting *	smr,
		xData_element *	element,
		tpia_1dData *	group
	)

Definition at line 289 of file tpia_misc.cc.

References tpia_1dData_s::data, tpia_1dData_s::end, tpia_1dData_s::length, tpia_1dData_s::start, and tpia_misc_get2d_xShared_yHistogram_data().

Referenced by tpia_channel_getFromElement(), and tpia_multiplicity_getFromElement().

                                                                                                                               {

    if( ( group->data = tpia_misc_get2d_xShared_yHistogram_data( smr, element, &(group->start), &(group->end), &(group->length) ) ) == NULL ) return( 1 );
    return( 0 );
}

double* tpia_misc_get2dx_y_data	(	statusMessageReporting *	smr,
		xData_element *	element,
		xData_Int *	length
	)

Definition at line 228 of file tpia_misc.cc.

References xData_2d_xy_allocateCopyData(), xData_addToAccessed(), xData_getOneElementByTagName(), xData_is_2d_xy(), and xData_element_s::xDataTypeInfo.

Referenced by tpia_multiplicity_getFromElement().

                                                                                                          {

    xData_element *xDataElement;
    double *d = NULL;
    xData_Int length_;

    xData_addToAccessed( smr, element, 1 );
    //if( ( xDataElement = xData_getOneElementByTagName( smr, element, "xData", 1 ) ) != NULL ) {
    if( ( xDataElement = xData_getOneElementByTagName( smr, element, (char*)"xData", 1 ) ) != NULL ) {
        xData_addToAccessed( smr, xDataElement, 1 );
        if( xData_is_2d_xy( smr, &(xDataElement->xDataTypeInfo), 1 ) ) {
            d = xData_2d_xy_allocateCopyData( smr, xDataElement, &length_ );
            *length = (xData_Int) length_;
        }
    }
    return( d );
}

double* tpia_misc_get2dxindex_y_data	(	statusMessageReporting *	smr,
		xData_element *	element,
		xData_Int *	start,
		xData_Int *	end,
		double *	xValues
	)

Definition at line 248 of file tpia_misc.cc.

References xDataType_s::end, xDataType_s::start, xData_2d_xindex_y_toFilledYs(), xData_addToAccessed(), xData_getOneElementByTagName(), xData_is_2d_xindex_y(), and xData_element_s::xDataTypeInfo.

Referenced by tpia_channel_getFromElement().

                                                                                                                                               {

    xData_element *xDataElement;
    double *d = NULL;

    xData_addToAccessed( smr, element, 1 );
    //if( ( xDataElement = xData_getOneElementByTagName( smr, element, "xData", 1 ) ) != NULL ) {
    if( ( xDataElement = xData_getOneElementByTagName( smr, element, (char*) "xData", 1 ) ) != NULL ) {
        xData_addToAccessed( smr, xDataElement, 1 );
        if( xData_is_2d_xindex_y( smr, &(xDataElement->xDataTypeInfo), 1 ) ) {
            if( start != NULL ) *start = xDataElement->xDataTypeInfo.start;
            if( end != NULL ) *end = xDataElement->xDataTypeInfo.end;
            d = xData_2d_xindex_y_toFilledYs( smr, xDataElement, xValues );
        }
    }
    return( d );
}

double tpia_misc_getPointwiseCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_1dData *	crossSection,
		double *	energyGrid,
		xData_Int	index,
		double	e_in
	)

Definition at line 298 of file tpia_misc.cc.

References tpia_1dData_s::data, and tpia_1dData_s::start.

Referenced by tpia_channel_getCrossSectionAtE(), and tpia_target_heated_getTotalCrossSectionAtE().

                                                                                                                                                      {

    double xsec = 0.0, e1, e2;

    if( ( index >= crossSection->start ) && ( index < crossSection->end ) ) {
        e1 = energyGrid[index];
        e2 = energyGrid[index + 1];
        index -= crossSection->start;
        if( e1 == e2 ) {                            /* Allow for future where step function may be allowed. */
            xsec = 0.5 * ( crossSection->data[index] + crossSection->data[index + 1] ); }
        else {
            xsec = ( crossSection->data[index] * ( e2 - e_in ) + crossSection->data[index + 1] * ( e_in - e1 ) ) / ( e2 - e1 );
        }
    }
    return( xsec );
}

int tpia_misc_NumberOfZSymbols

(

void

)

Definition at line 76 of file tpia_misc.cc.

Referenced by tpia_misc_symbolToZ(), and tpia_misc_ZToSymbol().

                                       {

    return( sizeof( ZSymbols ) / sizeof( struct ZSymbol ) );
}

char* tpia_misc_pointerToAttributeIfAllOk	(	statusMessageReporting *	smr,
		xData_element *	element,
		const char *	path,
		int	required,
		xData_attributionList *	attributes,
		const char *	name,
		const char *	file,
		int	line
	)

Definition at line 163 of file tpia_misc.cc.

References smr_isOk(), smr_setMessageError(), tpia_misc_setMessageError_Element(), and xData_getAttributesValue().

                                                                                          {

    char *value;

    if( !smr_isOk( smr ) ) return( NULL );
    if( ( value = xData_getAttributesValue( attributes, name ) ) == NULL ) {
        if( required ) {
            if( element != NULL ) {
                tpia_misc_setMessageError_Element( smr, NULL, element, file, line, 1, "element does not have attribute named %s", name ); }
            else {
                smr_setMessageError( smr, NULL, file, line, 1, "element does not have attribute named %s for file = %d", name, path );
            }
        }
    }
    return( value );
}

int tpia_misc_sampleEqualProbableBin	(	statusMessageReporting *	smr,
		tpia_decaySamplingInfo *	decaySamplingInfo,
		double	e_in,
		int	nBins,
		tpia_EqualProbableBinSpectra *	binned,
		double *	value_
	)

Definition at line 420 of file tpia_misc.cc.

References tpia_samplingMethods_s::angular_equalProbableBinMethod, tpia_EqualProbableBinSpectrum_s::bins, tpia_EqualProbableBinSpectra_s::energies, int(), tpia_EqualProbableBinSpectra_s::numberOfEs, tpia_decaySamplingInfo_s::rng, tpia_decaySamplingInfo_s::rngState, tpia_decaySamplingInfo_s::samplingMethods, tpia_misc_drng(), tpia_samplingMethods_isLinear, and tpia_EqualProbableBinSpectrum_s::value.

Referenced by tpia_angular_SampleMu(), tpia_angularEnergy_SampleEp(), and tpia_Legendre_SampleEp().

                                                               {

    int i, j, index1, index2, method = 0;
    double fE = 1., r, value1, value2, value3, P12, P23, value = -2.;

    if( e_in <= binned->energies[0].value ) {
                index1 = 0;
                index2 = 0; }
    else if( e_in >= binned->energies[binned->numberOfEs-1].value ) {
                index1 = binned->numberOfEs - 1;
                index2 = binned->numberOfEs - 1; }
    else {
            for( i = 0; i < binned->numberOfEs - 1; i++ ) {
                if( e_in <= binned->energies[i].value ) break;
            }
            i--;
            index1 = i;
            index2 = i;
            if( e_in != binned->energies[i].value ) {
                index2++;
                fE = ( e_in - binned->energies[i].value ) / ( binned->energies[i+1].value - binned->energies[i].value );
            }
    }
    r = tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState );
    j = (int) (r * nBins);
    if( j >= nBins ) j = nBins - 1;
    if( j < 0 ) j = 0;
    r = tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState );          // Do not change r until after Point1 below.
    if( tpia_samplingMethods_isLinear( decaySamplingInfo->samplingMethods->angular_equalProbableBinMethod ) ) {
        method = 1;
        if( ( ( j == 0 ) && ( r <= 0.5 ) ) || ( j == ( nBins - 1 ) && r > 0.5 ) ) method = 0;
    }
    if( method == 0 ) {                 /* Constant interpolaton. */
        value1 = ( 1. - fE ) * binned->energies[index1].bins[j] + fE * binned->energies[index2].bins[j];
        value2 = ( 1. - fE ) * binned->energies[index1].bins[j+1] + fE * binned->energies[index2].bins[j+1];
        value = ( 1. - r ) * value1 + r * value2; }
    else {                              /* Linear interpolation. */
        if( r <= 0.5 ) j--;             /* Point1: Above test insures that j is not 0 (nBins-1) if r <= 0.5 (> 0.5); */
        value1 = ( 1. - fE ) * binned->energies[index1].bins[j] + fE * binned->energies[index2].bins[j];
        value2 = ( 1. - fE ) * binned->energies[index1].bins[j+1] + fE * binned->energies[index2].bins[j+1];
        value3 = ( 1. - fE ) * binned->energies[index1].bins[j+2] + fE * binned->energies[index2].bins[j+2];
        P12 = 1. / ( value2 - value1 );
        P23 = 1. / ( value3 - value2 );
        r = tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState );
        if( 0.25 * ( 1.0 + 2.0 * ( value2 - value1 ) / ( value3 - value1 ) ) > r ) {
            P23 = 2. / ( value3 - value1 );
            value3 = value2; }
        else {
            P12 = 2. / ( value3 - value1 );
            value1 = value2;
        }
        r = tpia_misc_drng( decaySamplingInfo->rng, decaySamplingInfo->rngState );
        if( P23 != P12 ) r = ( -P12 + std::sqrt( P12 * P12 * ( 1. - r ) + r * P23 * P23 ) ) / ( P23 - P12 );
        value = 0.5 * ( value1 + value2 + r * ( value3 - value1 ) );
    }
    *value_ = value;
    return( 0 );
}

int tpia_misc_symbolToZ

(

const char *

Z

)

Definition at line 92 of file tpia_misc.cc.

References n, symbol, and tpia_misc_NumberOfZSymbols().

Referenced by tpia_miscNameToZAm().

                                         {

    int i, n = tpia_misc_NumberOfZSymbols( );

    for( i = 0; i < n; i++ ) {
        //if( !strcmp( Z, ZSymbols[i].symbol ) ) return( ZSymbols[i].Z );
        if( !strcmp( Z, ZSymbols[i].symbol.c_str() ) ) return( ZSymbols[i].Z );
    }
    return( -1 );
}

const char* tpia_misc_ZToSymbol

(

int

iZ

)

Definition at line 83 of file tpia_misc.cc.

References symbol, and tpia_misc_NumberOfZSymbols().

Referenced by G4GIDI_Misc_Z_A_m_ToName().

                                          {

    if( ( iZ < 0 ) || ( iZ >= tpia_misc_NumberOfZSymbols( ) ) ) return( NULL );
    //return( ZSymbols[iZ].symbol );
    return( ZSymbols[iZ].symbol.c_str() );
}

int tpia_miscNameToZAm	(	statusMessageReporting *	smr,
		const char *	name,
		int *	Z,
		int *	A,
		int *	m
	)

Definition at line 105 of file tpia_misc.cc.

References int(), n, smr_isOk(), smr_setMessageError(), and tpia_misc_symbolToZ().

Referenced by G4GIDI_Misc_channelCompound(), and tpia_particle_getInternalID().

                                                                                                {

    int i, n;
    const char *p;
    char s[1024] = "", *q, *e;   /* Note 1) routine will fail when parts of a particle name can be longer than 1024. */

    n = sizeof( s ) - 1;

    *Z = *A = *m = 0;
    if( !strncmp( "FissionProduct", name, 14 ) ) {
        *Z = 99;
        *A = 120;
        return( 0 );
    }
    if( !strcmp( "gamma", name ) ) return( 0 );
    for( p = name, q = s, i = 0; ( *p != '_' ) && ( i != n ); p++, q++, i++ ) *q = *p;
    if( i == n ) {              /* See note 1 above. */
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Failed to find first '_' in particle name %s", name ); }
    else {
        *q = 0;
        if( ( *Z = tpia_misc_symbolToZ( s ) ) < 0 ) {
            smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Particle %s's symbol = '%s' not found", name, s ); }
        else {                  /* Getting here implies that *p == '_'. */
            for( p++, q = s; ( *p != '_' ) && ( *p != 0 ) && ( i != n ); p++, q++, i++ ) *q = *p;
            if( i == n ) {      /* See note 1 above. */
                smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Failed to find second '_' in particle name %s", name ); }
            else {
                *q = 0;
                if( strcmp( s, "natural" ) == 0 ) {
                    e = s;
                    while( *e ) e++; }
                else {
                    *A = (int) strtol( s, &e, 10 );
                }
                if( *e != 0 ) {
                    smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Failed to convert A to integer in particle name %s", name ); }
                else {          /* Getting here implies that *p == '_' or 0. */
                    *m = 0;
                    if( *p == '_' ) {
                        p++;
                        if( *p != 'm' ) {
                            smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Particle name %s missing meta-stable label 'm'", name ); }
                        else {
                            p++;
                            *m = (int) strtol( p, &e, 10 );
                            if( *e != 0 ) smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Failed to convert m to integer in particle name %s", name );
                        }
                    }
                }
            }
        }
    }

    return( !smr_isOk( smr ) );
}

tpia_multiplicity* tpia_multiplicity_create

(

statusMessageReporting *

smr

)

Definition at line 50 of file tpia_multiplicity.cc.

References tpia_multiplicity_free(), tpia_multiplicity_initialize(), and xData_malloc2.

Referenced by tpia_multiplicity_createGetFromElement().

                                                                           {

    tpia_multiplicity *multiplicity;

    //if( ( multiplicity = xData_malloc2( smr, sizeof( tpia_multiplicity ), 0, "multiplicity" ) ) == NULL ) return( NULL );
    if( ( multiplicity = (tpia_multiplicity*) xData_malloc2( smr, sizeof( tpia_multiplicity ), 0, "multiplicity" ) ) == NULL ) return( NULL );
    if( tpia_multiplicity_initialize( smr, multiplicity ) ) multiplicity = tpia_multiplicity_free( smr, multiplicity );
    return( multiplicity );
}

tpia_multiplicity* tpia_multiplicity_createGetFromElement	(	statusMessageReporting *	smr,
		xData_element *	multiplicityElement,
		int	nGroups
	)

Definition at line 92 of file tpia_multiplicity.cc.

References tpia_multiplicity_create(), tpia_multiplicity_free(), and tpia_multiplicity_getFromElement().

                                                                                                                                          {

    tpia_multiplicity *multiplicity;

    if( ( multiplicity = tpia_multiplicity_create( smr ) ) == NULL ) return( NULL );
    if( tpia_multiplicity_getFromElement( smr, multiplicityElement, multiplicity, nGroups ) != 0 ) multiplicity = tpia_multiplicity_free( smr, multiplicity  );
    return( multiplicity );
}

tpia_multiplicity* tpia_multiplicity_free	(	statusMessageReporting *	smr,
		tpia_multiplicity *	multiplicity
	)

Definition at line 71 of file tpia_multiplicity.cc.

References tpia_multiplicity_release(), and xData_free().

Referenced by tpia_multiplicity_create(), tpia_multiplicity_createGetFromElement(), and tpia_product_release().

                                                                                                          {

    tpia_multiplicity_release( smr, multiplicity );
    xData_free( smr, multiplicity );
    return( NULL );
}

int tpia_multiplicity_getFromElement	(	statusMessageReporting *	smr,
		xData_element *	multiplicityElement,
		tpia_multiplicity *	multiplicity,
		int	nGroups
	)

Definition at line 104 of file tpia_multiplicity.cc.

References tpia_multiplicity_s::frame, tpia_multiplicity_s::grouped, xData_element_s::name, tpia_multiplicity_s::numberOfPointwise, tpia_multiplicity_s::pointwise, tpia_multiplicity_s::timeScale, tpia_frame_setFromElement(), tpia_misc_get2d_xShared_yHistogram_data_Grouped(), tpia_misc_get2dx_y_data(), tpia_misc_setMessageError_Element(), tpia_multiplicity_getTimeScaleFromElement(), xData_addToAccessed(), xData_getFirstElement(), and xData_getNextElement().

Referenced by tpia_multiplicity_createGetFromElement().

                                                                                                                                              {

    const char *timeScale;
    int isDelayedNeutrons;
    xData_element *data;

    xData_addToAccessed( smr, multiplicityElement, 1 );
    if( ( tpia_frame_setFromElement( smr, multiplicityElement, 2, &(multiplicity->frame) ) ) != 0 ) return( 1 );
    if( tpia_multiplicity_getTimeScaleFromElement( smr, multiplicityElement, &timeScale, &isDelayedNeutrons, &(multiplicity->timeScale) ) ) return( 1 );
    for( data = xData_getFirstElement( multiplicityElement ); data != NULL; data = xData_getNextElement( data ) ) {
        if( strcmp( data->name, "grouped" ) == 0 ) {
            if( tpia_misc_get2d_xShared_yHistogram_data_Grouped( smr, data, &(multiplicity->grouped) ) ) return( 1 ); }
        else if( strcmp( data->name, "pointwise" ) == 0 ) {
            if( ( multiplicity->pointwise = tpia_misc_get2dx_y_data( smr, data, &multiplicity->numberOfPointwise ) ) == NULL ) return( 1 ); }
        else {
            tpia_misc_setMessageError_Element( smr, NULL, multiplicityElement, __FILE__, __LINE__, 1, "unsupported multiplicity type = %s", data->name );
            return( 1 );
        }
    }
    return( 0 );
}

int tpia_multiplicity_getTimeScaleFromElement	(	statusMessageReporting *	smr,
		xData_element *	element,
		const char **	timeScale,
		int *	isDelayedNeutrons,
		double *	dTimeScale
	)

Definition at line 128 of file tpia_multiplicity.cc.

References xData_element_s::attributes, tpia_misc_setMessageError_Element(), and xData_getAttributesValue().

Referenced by tpia_multiplicity_getFromElement().

                             {

    const char *p;
    char *e;

    *isDelayedNeutrons = 0;
    *dTimeScale = 0.;
    *timeScale = xData_getAttributesValue( &(element->attributes), "timeScale" );
    if( *timeScale != NULL ) {
        if( strcmp( *timeScale, "prompt" ) ) {
            if( strncmp( *timeScale, "delayed", 7 ) ) {
                tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "invalid timeScale attribute = %s", *timeScale );
                return( 1 ); }
            else {
                for( p = &((*timeScale)[7]); *p; p++ ) if( !isspace( *p ) ) break;
                if( *p != '=' ) {
                    tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "timeScale attribute '%s' missing '='", *timeScale );
                    return( 1 );
                }
                p++;
                *dTimeScale = strtod( p, &e );
                if( *e != 0 ) {
                    tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "could not convert timeScale attribute '%s' to double",
                        *timeScale );
                    return( 1 );
                }
                *isDelayedNeutrons = 1;
            }
        }
    }
    return( 0 );
}

int tpia_multiplicity_initialize	(	statusMessageReporting *	smr,
		tpia_multiplicity *	multiplicity
	)

Definition at line 62 of file tpia_multiplicity.cc.

References tpia_multiplicity_s::frame, and tpia_frame_setFromString().

Referenced by tpia_multiplicity_create(), and tpia_multiplicity_release().

                                                                                                 {

    memset( multiplicity, 0, sizeof( tpia_multiplicity ) );
    if( tpia_frame_setFromString( smr, "", "", 0, &(multiplicity->frame) ) ) return( 1 );
    return( 0 );
}

int tpia_multiplicity_release	(	statusMessageReporting *	smr,
		tpia_multiplicity *	multiplicity
	)

Definition at line 80 of file tpia_multiplicity.cc.

References tpia_1dData_s::data, tpia_multiplicity_s::grouped, tpia_multiplicity_s::pointwise, tpia_multiplicity_initialize(), and xData_free().

Referenced by tpia_multiplicity_free().

                                                                                              {

    //multiplicity->grouped.data = xData_free( smr, multiplicity->grouped.data );
    multiplicity->grouped.data = (double*) xData_free( smr, multiplicity->grouped.data );
    //multiplicity->pointwise = xData_free( smr, multiplicity->pointwise );
    multiplicity->pointwise = (double*) xData_free( smr, multiplicity->pointwise );
    tpia_multiplicity_initialize( smr, multiplicity );
    return( 0 );
}

tpia_particle* tpia_particle_create

(

statusMessageReporting *

smr

)

Definition at line 52 of file tpia_particle.cc.

References tpia_particle_initialize(), and xData_malloc2.

Referenced by tpia_particle_getInternalID().

                                                                   {

    //tpia_particle *particle = xData_malloc2( smr, sizeof( tpia_particle ), 0, "particle" );
    tpia_particle *particle = (tpia_particle*) xData_malloc2( smr, sizeof( tpia_particle ), 0, "particle" );

    if( particle == NULL ) return( NULL );
    tpia_particle_initialize( smr, particle );
    return( particle );
}

tpia_particle* tpia_particle_free	(	statusMessageReporting *	smr,
		tpia_particle *	particle
	)

Definition at line 73 of file tpia_particle.cc.

References tpia_particle_s::next, tpia_particle_s::prior, tpia_particle_release(), and xData_free().

Referenced by tpia_particle_freeInternalList(), and tpia_particle_getInternalID().

                                                                                          {

    int i, j;
    tpia_particle **p;

    for( i = 0, p = particleSortedList; i < numberOfParticles; i++, p++ ) {
        if( *p == particle ) {
            numberOfParticles--;
            for( j = i; j < numberOfParticles; j++, p++ ) *p = p[1];
            break;
        }
    }
    if( particle == particleListEnd ) particleListEnd = particle->prior;
    if( particle == particleList ) particleList = particle->next;
    if( particle->prior != NULL ) particle->prior->next = particle->next;
    if( particle->next != NULL ) particle->next->prior = particle->prior;
    tpia_particle_release( smr, particle );
    xData_free( smr, particle );
    return( NULL );
}

int tpia_particle_freeInternalList

(

statusMessageReporting *

smr

)

Definition at line 104 of file tpia_particle.cc.

References tpia_particle_free(), and xData_free().

                                                                  {

    while( particleList != NULL ) tpia_particle_free( smr, particleList );
    //particleSortedList = xData_free( smr, particleSortedList );
    particleSortedList = (tpia_particle**) xData_free( smr, particleSortedList );
    return( 0 );
}

tpia_particle* tpia_particle_getInternalID	(	statusMessageReporting *	smr,
		const char *const	name
	)

Definition at line 114 of file tpia_particle.cc.

References tpia_particle_s::A, tpia_particle_s::fullMass_MeV, tpi_spectralID_s::level, python.hepunit::m, tpia_particle_s::m, tpia_particle_s::mass, G4INCL::Math::max(), G4INCL::Math::min(), tpi_spectralID_s::name, tpia_particle_s::name, tpia_particle_s::next, tpia_particle_s::ordinal, tpia_particle_s::prior, smr_isOk(), tpia_particle_s::spectralID, tpi_spectralID_parse(), tpia_AMU2MeV, tpia_miscNameToZAm(), tpia_particle_create(), tpia_particle_free(), tpia_particleMass_AMU(), xData_realloc2, and tpia_particle_s::Z.

Referenced by tpia_product_getFromElement(), tpia_target_heated_read(), and tpia_target_read().

                                                                                                   {

    int i, iCmp, min, mid, max, Z, A, m;
    tpia_particle *particle;
    char *EOP;

    iCmp = 0;
    min = mid = 0;
    max = numberOfParticles;
    while( min != max ) {
        mid = ( min + max ) / 2;
        iCmp = strcmp( name, particleSortedList[mid]->name );
        if( iCmp == 0 ) return( particleSortedList[mid] );
        if( iCmp < 0 ) {
            max = mid - 1;
            if( mid == 0 ) max = 0; }
        else {
            min = mid + 1;
            if( min > max ) min = max;
        }
    }
    mid = min;
    if( numberOfParticles > 0 ) {
        iCmp = strcmp( name, particleSortedList[mid]->name );
        if( iCmp == 0 ) return( particleSortedList[mid] );
        if( ( iCmp < 0 ) && ( mid != 0 ) ) {
            mid--;
            iCmp = strcmp( name, particleSortedList[mid]->name );
        }
    }

    if( ( particle = tpia_particle_create( smr ) ) == NULL ) return( NULL );
    if( ( particle->spectralID = tpi_spectralID_parse( smr, name, &(EOP) ) ) == NULL ) return( tpia_particle_free( smr, particle ) );
    particle->name = particle->spectralID->name;
    if( tpia_miscNameToZAm( smr, particle->name, &Z, &A, &m ) != 0 ) return( tpia_particle_free( smr, particle ) );
    particle->prior = NULL;
    particle->next = NULL;
    particle->Z = Z;
    particle->A = A;
    particle->m = m;
    particle->mass = tpia_particleMass_AMU( smr, particle->name );
    if( !smr_isOk( smr ) ) return( tpia_particle_free( smr, particle ) );
    particle->fullMass_MeV = tpia_AMU2MeV * particle->mass + particle->spectralID->level;

    if( sizeOfParticleSortedList < ( numberOfParticles + 1 ) ) {
        if( sizeOfParticleSortedList == 0 ) {
            sizeOfParticleSortedList = initialSizeOfList; }
        else {
            sizeOfParticleSortedList += incrementalSizeOfList;
        }
        //if( ( particleSortedList = xData_realloc2( smr, particleSortedList, sizeOfParticleSortedList * sizeof( tpia_particle * ), 
        if( ( particleSortedList = (tpia_particle** ) xData_realloc2( smr, particleSortedList, sizeOfParticleSortedList * sizeof( tpia_particle * ), 
            "particleSortedList" ) ) == NULL ) return( tpia_particle_free( smr, particle ) );
    }

    if( particleList == NULL ) {
        particle->ordinal = 0;
        particleListEnd = particleList = particle; }
    else {
        particle->ordinal = particleListEnd->ordinal + 1;
        particle->prior = particleListEnd;
        particleListEnd->next = particle;
        particleListEnd = particle;
    }

    if( ( mid != 0 ) || ( iCmp > 0 ) ) mid++;
    for( i = numberOfParticles; i > mid; i-- ) particleSortedList[i] = particleSortedList[i-1];
    particleSortedList[mid] = particle;
    numberOfParticles++;

    return( particle );
}

int tpia_particle_initialize	(	statusMessageReporting *	smr,
		tpia_particle *	particle
	)

Definition at line 65 of file tpia_particle.cc.

Referenced by tpia_particle_create().

                                                                                  {

    memset( particle, 0, sizeof( tpia_particle ) );
    return( 0 );
}

int tpia_particle_printInternalSortedList

(

statusMessageReporting *

smr

)

Definition at line 190 of file tpia_particle.cc.

                                                                      {

    int i;

    for( i = 0; i < numberOfParticles; i++ ) printf( "%s\n", particleSortedList[i]->name );
    return( 0 );
}

int tpia_particle_release	(	statusMessageReporting *	smr,
		tpia_particle *	particle
	)

Definition at line 96 of file tpia_particle.cc.

References tpia_particle_s::spectralID, and tpi_spectralID_free().

Referenced by tpia_particle_free().

                                                                                  {

    if( particle->spectralID != NULL ) tpi_spectralID_free( smr, particle->spectralID );
    return( 0 );
}

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 *	smr,
		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 *	smr,
		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 *	smr,
		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 *	smr,
		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 *	smr,
		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 *	smr,
		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 );
}

int tpia_samplingMethods_initialize	(	statusMessageReporting *	smr,
		tpia_samplingMethods *	samplingMethods
	)

Definition at line 53 of file tpia_samplingMethods.cc.

References tpia_samplingMethods_s::angular_equalProbableBinMethod, and tpia_samplingMethods_linear.

Referenced by tpia_target_initialize().

                                                                                                       {

    memset( samplingMethods, 0, sizeof( tpia_samplingMethods ) );
    samplingMethods->angular_equalProbableBinMethod = tpia_samplingMethods_linear;
    return( 0 );
}

tpia_target* tpia_target_create

(

statusMessageReporting *

smr

)

Definition at line 49 of file tpia_target.cc.

References tpia_target_free(), tpia_target_initialize(), and xData_malloc2.

Referenced by tpia_target_createRead().

                                                               {

    tpia_target *target;

    //if( ( target = xData_malloc2( smr, sizeof( tpia_target ), 0, "target" ) ) == NULL ) return( NULL );
    if( ( target = (tpia_target*) xData_malloc2( smr, sizeof( tpia_target ), 0, "target" ) ) == NULL ) return( NULL );
    if( tpia_target_initialize( smr, target ) ) target = tpia_target_free( smr, target );
    return( target );
}

tpia_target* tpia_target_createRead	(	statusMessageReporting *	smr,
		const char *	fileName
	)

Definition at line 70 of file tpia_target.cc.

References tpia_target_create(), tpia_target_read(), and xData_free().

Referenced by G4GIDI_target::init(), and tpia_target_createReadFromMap().

                                                                                         {

    tpia_target *target;

    if( ( target = tpia_target_create( smr ) ) == NULL ) return( NULL );
    //if( tpia_target_read( smr, target, fileName ) != 0 ) target = xData_free( smr, target );
    if( tpia_target_read( smr, target, fileName ) != 0 ) target = (tpia_target*) xData_free( smr, target );
    return( target );
}

tpia_target* tpia_target_createReadFromMap	(	statusMessageReporting *	smr,
		tpia_map *	map,
		const char *	evaluation,
		const char *	projectileName,
		const char *	targetName
	)

Definition at line 93 of file tpia_target.cc.

References tpia_map_findTarget(), tpia_target_createRead(), and xData_free().

Referenced by tpia_depot_addTargetFromMap().

                                 {

    char *targetPath;
    tpia_target *target;

    targetPath = tpia_map_findTarget( smr, map, evaluation, projectileName, targetName );
    if( targetPath == NULL ) return( NULL );
    target = tpia_target_createRead( smr, targetPath );
    xData_free( smr, targetPath );
    return( target );
}

tpia_target* tpia_target_free	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 108 of file tpia_target.cc.

References tpia_target_release(), and xData_free().

Referenced by tpia_depot_release(), tpia_target_create(), and G4GIDI_target::~G4GIDI_target().

                                                                                  {

    tpia_target_release( smr, target );
    xData_free( smr, target );
    return( NULL );
}

char* tpia_target_getAttributesValue	(	statusMessageReporting *	smr,
		tpia_target *	target,
		char const *	name
	)

Definition at line 248 of file tpia_target.cc.

References tpia_target_s::attributes, and xData_getAttributesValue().

                                                                                                        {

    return( xData_getAttributesValue( &(target->attributes), name ) );
}

xData_Int tpia_target_getEnergyGridAtTIndex	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		double **	energyGrid
	)

Definition at line 331 of file tpia_target.cc.

References smr_isOk(), tpia_target_getHeatedTargetAtTIndex(), and tpia_target_heated_getEnergyGrid().

Referenced by G4GIDI_target::getEnergyGridAtTIndex().

                                                                                                                                {

    tpia_target_heated *heated = tpia_target_getHeatedTargetAtTIndex( smr, target, index );

    if( !smr_isOk( smr ) ) return( -1 );
    return( tpia_target_heated_getEnergyGrid( smr, heated, energyGrid ) );
}

tpia_target_heated* tpia_target_getHeatedTargetAtIndex_ReadIfNeeded	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index
	)

Definition at line 290 of file tpia_target.cc.

References tpia_target_heated_info_s::heatedTarget, tpia_target_s::heatedTargets, tpia_target_s::nHeatedTargets, smr_setMessageError(), and tpia_target_readHeatedTarget().

                                                                                                                                   {

    if( ( index < 0 ) || ( index >= target->nHeatedTargets ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "temperature index = %d out of range (0 <= index < %d", index, target->nHeatedTargets );
        return( NULL );
    }
    if( target->heatedTargets[index].heatedTarget == NULL ) tpia_target_readHeatedTarget( smr, target, index, 0 );
    return( target->heatedTargets[index].heatedTarget );
}

double tpia_target_getIndexChannelCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		double	T,
		xData_Int	iEg,
		double	e_in,
		int	crossSectionType
	)

Definition at line 380 of file tpia_target.cc.

References tpia_target_heated_info_s::heatedTarget, tpia_target_s::nReadHeatedTargets, tpia_target_s::readHeatedTargets, tpia_target_heated_info_s::temperature, and tpia_target_heated_getIndexChannelCrossSectionAtE().

Referenced by G4GIDI_target::sampleChannelCrossSectionAtE(), and G4GIDI_target::sumChannelCrossSectionAtE().

                           {

    int i;
    double xsec = 0., xsec1, xsec2;

    for( i = 0; i < target->nReadHeatedTargets; i++ ) if( target->readHeatedTargets[i]->temperature > T ) break;
    if( i == 0 ) {
        xsec = tpia_target_heated_getIndexChannelCrossSectionAtE( smr, target->readHeatedTargets[0]->heatedTarget, index, iEg, e_in, crossSectionType ); }
    else if( i == target->nReadHeatedTargets ) {
        xsec = tpia_target_heated_getIndexChannelCrossSectionAtE( smr, target->readHeatedTargets[i-1]->heatedTarget, index, iEg, e_in, crossSectionType ); }
    else {
        xsec1 = tpia_target_heated_getIndexChannelCrossSectionAtE(smr, target->readHeatedTargets[i-1]->heatedTarget, index, iEg, e_in, crossSectionType);
        xsec2 = tpia_target_heated_getIndexChannelCrossSectionAtE(smr, target->readHeatedTargets[i  ]->heatedTarget, index, iEg, e_in, crossSectionType);
        xsec = ( ( target->readHeatedTargets[i]->temperature - T ) * xsec1 + ( T - target->readHeatedTargets[i-1]->temperature ) * xsec2 ) / 
            ( target->readHeatedTargets[i]->temperature - target->readHeatedTargets[i-1]->temperature );
    }

    return( xsec );
}

int tpia_target_getTemperatures	(	statusMessageReporting *	smr,
		tpia_target *	target,
		double *	temperatures
	)

Definition at line 256 of file tpia_target.cc.

References tpia_target_s::heatedTargets, tpia_target_s::nHeatedTargets, and tpia_target_heated_info_s::temperature.

Referenced by G4GIDI_target::getTemperatures().

                                                                                                       {

    int i;

    if( temperatures != NULL ) for( i = 0; i < target->nHeatedTargets; i++ ) temperatures[i] = target->heatedTargets[i].temperature;
    return( target->nHeatedTargets );
}

double tpia_target_getTotalCrossSectionAtTAndE	(	statusMessageReporting *	smr,
		tpia_target *	target,
		double	T,
		xData_Int	iEg,
		double	e_in,
		int	crossSectionType
	)

Definition at line 357 of file tpia_target.cc.

References tpia_target_heated_info_s::heatedTarget, tpia_target_s::nReadHeatedTargets, tpia_target_s::readHeatedTargets, tpia_target_heated_info_s::temperature, and tpia_target_heated_getTotalCrossSectionAtE().

Referenced by G4GIDI_target::getTotalCrossSectionAtE().

                           {

    int i;
    double xsec = 0., xsec1, xsec2;

    for( i = 0; i < target->nReadHeatedTargets; i++ ) if( target->readHeatedTargets[i]->temperature > T ) break;
    if( i == 0 ) {
        xsec = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[0]->heatedTarget, iEg, e_in, crossSectionType ); }
    else if( i == target->nReadHeatedTargets ) {
        xsec = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[i-1]->heatedTarget, iEg, e_in, crossSectionType ); }
    else {
        xsec1 = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[i-1]->heatedTarget, iEg, e_in, crossSectionType );
        xsec2 = tpia_target_heated_getTotalCrossSectionAtE( smr, target->readHeatedTargets[i  ]->heatedTarget, iEg, e_in, crossSectionType );
        xsec = ( ( target->readHeatedTargets[i]->temperature - T ) * xsec1 + ( T - target->readHeatedTargets[i-1]->temperature ) * xsec2 ) / 
            ( target->readHeatedTargets[i]->temperature - target->readHeatedTargets[i-1]->temperature );
    }

    return( xsec );
}

tpia_1dData* tpia_target_getTotalCrossSectionAtTIndex	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		int	crossSectionType
	)

Definition at line 341 of file tpia_target.cc.

References smr_isOk(), smr_setMessageError(), tpia_target_heated_s::totalCrossSectionGrouped, tpia_target_heated_s::totalCrossSectionPointwise, tpia_crossSectionType_grouped, tpia_crossSectionType_pointwise, and tpia_target_getHeatedTargetAtTIndex().

                                                                                                                                           {

    tpia_target_heated *heated = tpia_target_getHeatedTargetAtTIndex( smr, target, index );

    if( !smr_isOk( smr ) ) return( NULL );
    if( crossSectionType == tpia_crossSectionType_grouped ) {
        return( &(heated->totalCrossSectionGrouped) ); }
    else if( crossSectionType == tpia_crossSectionType_pointwise ) {
        return( &(heated->totalCrossSectionPointwise) );
    }
    smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "Invalue crossSectionType = %d", crossSectionType );
    return( NULL );
}

tpia_target_heated* tpia_target_heated_create

(

statusMessageReporting *

smr

)

Definition at line 52 of file tpia_target_heated.cc.

References tpia_target_heated_initialize(), xData_free(), and xData_malloc2.

Referenced by tpia_target_heated_createRead().

                                                                             {

    tpia_target_heated *target;

    //if( ( target = xData_malloc2( smr, sizeof( tpia_target_heated ), 0, "target" ) ) == NULL ) return( NULL );
    if( ( target = (tpia_target_heated*) xData_malloc2( smr, sizeof( tpia_target_heated ), 0, "target" ) ) == NULL ) return( NULL );
    //if( tpia_target_heated_initialize( smr, target ) ) target = xData_free( smr, target );
    if( tpia_target_heated_initialize( smr, target ) ) target = (tpia_target_heated*) xData_free( smr, target );
    return( target );
}

tpia_target_heated* tpia_target_heated_createRead	(	statusMessageReporting *	smr,
		const char *	fileName,
		int	checkElememtsForAccess
	)

Definition at line 74 of file tpia_target_heated.cc.

References tpia_target_heated_create(), tpia_target_heated_read(), and xData_free().

Referenced by tpia_target_readHeatedTarget().

                                                                                                                                   {

    tpia_target_heated *target;

    if( ( target = tpia_target_heated_create( smr ) ) == NULL ) return( NULL );
    //if( tpia_target_heated_read( smr, target, fileName, checkElememtsForAccess ) != 0 ) target = xData_free( smr, target );
    if( tpia_target_heated_read( smr, target, fileName, checkElememtsForAccess ) != 0 ) target = (tpia_target_heated*) xData_free( smr, target );
    return( target );
}

tpia_target_heated* tpia_target_heated_free	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 86 of file tpia_target_heated.cc.

References tpia_target_heated_release(), and xData_free().

Referenced by tpia_target_release().

                                                                                                       {

    tpia_target_heated_release( smr, target );
    xData_free( smr, target );
    return( NULL );
}

tpia_channel* tpia_target_heated_getChannelAtIndex	(	tpia_target_heated *	target,
		int	index
	)

Definition at line 263 of file tpia_target_heated.cc.

References tpia_target_heated_s::channels.

Referenced by G4GIDI_target::init(), and tpia_target_heated_getChannelAtIndex_smr().

                                                                                            {

    tpia_channel *channel = NULL;

    if( ( index >= 0 ) && ( index < target->nChannels ) ) channel = target->channels[index];
    return( channel );
}

tpia_channel* tpia_target_heated_getChannelAtIndex_smr	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		int	index
	)

Definition at line 273 of file tpia_target_heated.cc.

References tpia_particle_s::name, tpia_target_heated_s::projectileID, smr_setMessageError(), tpia_target_heated_s::targetID, and tpia_target_heated_getChannelAtIndex().

Referenced by G4GIDI_target::getElasticFinalState(), tpia_target_heated_getIndexChannelCrossSectionAtE(), and tpia_target_heated_sampleIndexChannelProductsAtE().

                                                                                                                             {

    tpia_channel *channel = tpia_target_heated_getChannelAtIndex( target, index );

    if( channel == NULL ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "bad channel index = %d for %s + %s", index, 
            target->projectileID->name, target->targetID->name ); 
    }
    return( channel );
}

xData_Int tpia_target_heated_getEIndex	(	tpia_target_heated *	target,
		double	e_in
	)

Definition at line 305 of file tpia_target_heated.cc.

References tpia_target_heated_s::energyGrid, tpia_target_heated_s::energyGridLength, and tpia_misc_binarySearch().

Referenced by tpia_channel_getCrossSectionAtE(), and tpia_target_heated_getTotalCrossSectionAtE().

                                                                                  {

    return( tpia_misc_binarySearch( target->energyGridLength, target->energyGrid, e_in ) );
}

xData_Int tpia_target_heated_getEnergyGrid	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		double **	energyGrid
	)

Definition at line 297 of file tpia_target_heated.cc.

References tpia_target_heated_s::energyGrid, and tpia_target_heated_s::energyGridLength.

Referenced by tpia_target_getEnergyGridAtTIndex().

                                                                                                                        {

    if( energyGrid != NULL ) *energyGrid = target->energyGrid;
    return( target->energyGridLength );
}

double tpia_target_heated_getIndexChannelCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		int	index,
		xData_Int	iEg,
		double	e_in,
		int	crossSectionType
	)

Definition at line 329 of file tpia_target_heated.cc.

References tpia_channel_getCrossSectionAtE(), and tpia_target_heated_getChannelAtIndex_smr().

Referenced by tpia_target_getIndexChannelCrossSectionAtE().

                               {

    double xsec = 0.;
    tpia_channel *channel = tpia_target_heated_getChannelAtIndex_smr( smr, target, index );

    if( channel != NULL ) xsec = tpia_channel_getCrossSectionAtE( smr, channel, iEg, e_in, crossSectionType );
    return( xsec );
}

tpia_channel* tpia_target_heated_getProductionChannelAtIndex	(	tpia_target_heated *	target,
		int	index
	)

Definition at line 286 of file tpia_target_heated.cc.

References tpia_target_heated_s::productionChannels.

                                                                                                      {

    tpia_channel *channel = NULL;

    if( ( index >= 0 ) && ( index < target->nProductionChannels ) ) channel = target->productionChannels[index];
    return( channel );
}

double tpia_target_heated_getTotalCrossSectionAtE	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		xData_Int	gE,
		double	e_in,
		int	crossSectionType
	)

Definition at line 313 of file tpia_target_heated.cc.

References tpia_target_heated_s::energyGrid, tpia_target_heated_s::totalCrossSectionPointwise, tpia_crossSectionType_grouped, tpia_crossSectionType_pointwise, tpia_misc_getPointwiseCrossSectionAtE(), and tpia_target_heated_getEIndex().

Referenced by tpia_target_getTotalCrossSectionAtTAndE().

                               {

    double xsec = 0.;

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

int tpia_target_heated_initialize	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 66 of file tpia_target_heated.cc.

Referenced by tpia_target_heated_create(), and tpia_target_heated_read().

                                                                                          {

    memset( target, 0, sizeof( tpia_target_heated ) );
    return( 0 );
}

int tpia_target_heated_numberOfChannels	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 248 of file tpia_target_heated.cc.

References tpia_target_heated_s::nChannels.

Referenced by tpia_target_numberOfChannels().

                                                                                                {

    return( target->nChannels );
}

int tpia_target_heated_numberOfProductionChannels	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 256 of file tpia_target_heated.cc.

References tpia_target_heated_s::nProductionChannels.

Referenced by tpia_target_numberOfProductionChannels().

                                                                                                          {

    return( target->nProductionChannels );
}

int tpia_target_heated_read	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		const char *	fileName,
		int	checkElememtsForAccess
	)

Definition at line 126 of file tpia_target_heated.cc.

References tpia_target_heated_s::absPath, xData_element_s::attributes, tpia_target_heated_s::attributes, tpia_target_heated_s::channels, tpia_target_heated_s::contents, tpia_channel_s::crossSectionGrouped, tpia_channel_s::crossSectionPointwise, tpia_1dData_s::data, int(), xData_element_s::name, tpia_target_heated_s::nChannels, tpia_target_heated_s::nGroups, tpia_target_heated_s::nProductionChannels, tpia_target_heated_s::path, tpia_target_heated_s::productionChannels, tpia_target_heated_s::projectileID, smr_isOk(), tpia_1dData_s::start, tpia_target_heated_s::targetID, tpia_target_heated_s::totalCrossSectionGrouped, tpia_target_heated_s::totalCrossSectionPointwise, tpia_channel_createGetFromElement(), tpia_misc_pointerToAttributeIfAllOk3, tpia_misc_setMessageError_Element(), tpia_particle_getInternalID(), tpia_target_heated_initialize(), tpia_target_heated_release(), xData_addToAccessed(), xData_convertAttributeTo_xData_Int(), xData_copyAttributionList(), xData_getAttributesValue(), xData_getDocumentsElement(), xData_getFirstElement(), xData_getNextElement(), xData_getOneElementByTagName(), xData_malloc2, xData_numberOfElementsByTagName(), xData_parseFree(), xData_parseReadFile(), xDataMisc_allocateCopyString2, and xDataMisc_getAbsPath().

Referenced by tpia_target_heated_createRead().

                                                                                                                                         {
/*
*   If a target has already been read into this target, user must have called tpia_target_heated_release before calling this routine.
*   Otherwise, there will be memory leaks.
*/
    xData_document *doc = NULL;
    xData_element *element, *channelElement, *channels;
    int nChannels;
    tpia_channel *channel;
    char *name;
    xData_Int i, j;

    tpia_target_heated_initialize( smr, target );
    if( ( target->path = xDataMisc_allocateCopyString2( smr, fileName, "path" ) ) == NULL ) return( 1 );
    if( ( target->absPath = xDataMisc_getAbsPath( smr, fileName ) ) == NULL ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
    if( ( doc = xData_parseReadFile( smr, fileName, NULL, NULL ) ) == NULL ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
    element = xData_getDocumentsElement( doc );
    xData_addToAccessed( smr, element, 1 );
    if( xData_convertAttributeTo_xData_Int( smr, element, "nGroups", &i ) != 0 ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
    target->nGroups = (int) i;
    if( strcmp( element->name, "xTargetHeated" ) != 0 ) {
        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "input file's top element must be xTargetHeated and not %s", 
            element->name ); }
    else {
        xData_copyAttributionList( smr, &(target->attributes), &(element->attributes) );
        if( smr_isOk( smr ) ) target->contents = xData_getAttributesValue( &(target->attributes), "contents" );
        if( ( name = tpia_misc_pointerToAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "projectile" ) ) != NULL )
            target->projectileID = tpia_particle_getInternalID( smr, name );
        if( ( name = tpia_misc_pointerToAttributeIfAllOk3( smr, target->absPath, 1, &(target->attributes), "target" ) ) != NULL )
            target->targetID = tpia_particle_getInternalID( smr, name );
        if( smr_isOk( smr ) ) _tpia_target_heated_getEnergyGridAndAllocateTotalCrossSections( smr, target, element );
        if( smr_isOk( smr ) ) {            /* Get channels. */
            //if( ( channels = xData_getOneElementByTagName( smr, element, "channels", 1 ) ) == NULL ) 
            if( ( channels = xData_getOneElementByTagName( smr, element, (char*)"channels", 1 ) ) == NULL ) 
                return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
            xData_addToAccessed( smr, channels, 1 );
            if( ( nChannels = xData_numberOfElementsByTagName( smr, channels, "channel" ) ) > 0 ) {
                //if( ( target->channels = xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) == NULL ) 
                if( ( target->channels = (tpia_channel**) xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) == NULL ) 
                    return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) );
                for( channelElement = xData_getFirstElement( channels ); channelElement != NULL; channelElement = xData_getNextElement( channelElement ) ) {
                    if( !strcmp( channelElement->name, "channel" ) ) {
                        if( ( channel = tpia_channel_createGetFromElement( smr, target, channelElement, 1 ) ) == NULL ) break;
                        target->channels[target->nChannels] = channel;
                        target->nChannels++;
                        for( i = channel->crossSectionPointwise.start, j = 0; i < channel->crossSectionPointwise.end; i++, j++ ) 
                            target->totalCrossSectionPointwise.data[i] += channel->crossSectionPointwise.data[j];
                        for( i = channel->crossSectionGrouped.start, j = 0; i < channel->crossSectionGrouped.end; i++, j++ ) 
                            target->totalCrossSectionGrouped.data[i] += channel->crossSectionGrouped.data[j];
                    }
                }
            }
        }
        if( smr_isOk( smr ) ) {            /* Get production channels. */
            //if( ( channels = xData_getOneElementByTagName( smr, element, "productionChannels", 0 ) ) == NULL ) {
            if( ( channels = xData_getOneElementByTagName( smr, element, (char*) "productionChannels", 0 ) ) == NULL ) {
                if( !smr_isOk( smr ) ) return( _tpia_target_heated_releaseAndReturnOne( smr, doc, target ) ); }
            else {
                xData_addToAccessed( smr, channels, 1 );
                if( ( nChannels = xData_numberOfElementsByTagName( smr, channels, "channel" ) ) > 0 ) {
                    //if( ( target->productionChannels = xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) != NULL ) {
                    if( ( target->productionChannels = (tpia_channel**) xData_malloc2( smr, nChannels * sizeof( tpia_channel * ), 1, "channels" ) ) != NULL ) {
                        for( channelElement = xData_getFirstElement(channels); channelElement != NULL; channelElement = xData_getNextElement(channelElement) ) {
                            if( !strcmp( channelElement->name, "channel" ) ) {
                                channel = tpia_channel_createGetFromElement( smr, target, channelElement, 1 );
                                if( channel == NULL ) break;
                                target->productionChannels[target->nProductionChannels] = channel;
                                target->nProductionChannels++;
                            }
                        }
                    }
                }
            }
        }
    }
    if( smr_isOk( smr ) && checkElememtsForAccess ) _tpia_target_heated_checkElememtsForAccess( smr, doc );
    xData_parseFree( smr, doc );
    if( !smr_isOk( smr ) ) tpia_target_heated_release( smr, target );
    return( !smr_isOk( smr ) );
}

int tpia_target_heated_release	(	statusMessageReporting *	smr,
		tpia_target_heated *	target
	)

Definition at line 95 of file tpia_target_heated.cc.

References tpia_target_heated_s::absPath, tpia_target_heated_s::attributes, tpia_target_heated_s::channels, tpia_1dData_s::data, tpia_target_heated_s::energyGrid, tpia_target_heated_s::energyGridLength, tpia_target_heated_s::kerma, tpia_target_heated_s::nChannels, tpia_target_heated_s::nProductionChannels, tpia_target_heated_s::path, tpia_target_heated_s::productionChannels, tpia_target_heated_s::totalCrossSectionGrouped, tpia_target_heated_s::totalCrossSectionPointwise, tpia_channel_free(), xData_free(), and xData_releaseAttributionList().

Referenced by tpia_target_heated_free(), and tpia_target_heated_read().

                                                                                          {

    int i;

    //target->path = xData_free( smr, target->path );
    target->path = (char*) xData_free( smr, target->path );
    //target->absPath = xData_free( smr, target->absPath );
    target->absPath = (char*) xData_free( smr, target->absPath );
    target->energyGridLength = 0;
    //target->energyGrid = xData_free( smr, target->energyGrid );
    target->energyGrid = (double*) xData_free( smr, target->energyGrid );
    //target->kerma = xData_free( smr, target->kerma );
    target->kerma = (double*) xData_free( smr, target->kerma );
    //target->totalCrossSectionPointwise.data = xData_free( smr, target->totalCrossSectionPointwise.data );
    target->totalCrossSectionPointwise.data = (double*) xData_free( smr, target->totalCrossSectionPointwise.data );
    //target->totalCrossSectionGrouped.data = xData_free( smr, target->totalCrossSectionGrouped.data );
    target->totalCrossSectionGrouped.data = (double*) xData_free( smr, target->totalCrossSectionGrouped.data );
    xData_releaseAttributionList( smr, &(target->attributes) );
    for( i = 0; i < target->nChannels; i++ ) tpia_channel_free( smr, target->channels[i] );
    target->nChannels = 0;
    //target->channels = xData_free( smr, target->channels );
    target->channels = (tpia_channel**) xData_free( smr, target->channels );
    for( i = 0; i < target->nProductionChannels; i++ ) tpia_channel_free( smr, target->productionChannels[i] );
    target->nProductionChannels = 0;
    //target->productionChannels = xData_free( smr, target->productionChannels );
    target->productionChannels = (tpia_channel**) xData_free( smr, target->productionChannels );
    return( 0 );
}

int tpia_target_heated_sampleIndexChannelProductsAtE	(	statusMessageReporting *	smr,
		tpia_target_heated *	target,
		int	index,
		tpia_decaySamplingInfo *	decaySamplingInfo,
		int	nProductData,
		tpia_productOutgoingData *	productData
	)

Definition at line 341 of file tpia_target_heated.cc.

References tpia_channel_s::decayChannel, tpia_decayChannel_sampleProductsAtE(), and tpia_target_heated_getChannelAtIndex_smr().

Referenced by G4GIDI_target::getFinalState(), and tpia_target_sampleIndexChannelProductsAtE().

                                                                                                              {

    tpia_channel *channel = tpia_target_heated_getChannelAtIndex_smr( smr, target, index );

    if( channel == NULL ) return( -1 );
    return( tpia_decayChannel_sampleProductsAtE( smr, &(channel->decayChannel), decaySamplingInfo, nProductData, productDatas ) );
}

int tpia_target_initialize	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 61 of file tpia_target.cc.

References tpia_target_s::samplingMethods, and tpia_samplingMethods_initialize().

Referenced by tpia_target_create(), tpia_target_read(), and tpia_target_release().

                                                                               {

    memset( target, 0, sizeof( tpia_target ) );
    tpia_samplingMethods_initialize( smr, &(target->samplingMethods) );
    return( 0 );
}

int tpia_target_numberOfChannels	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 317 of file tpia_target.cc.

References tpia_target_s::baseHeatedTarget, and tpia_target_heated_numberOfChannels().

Referenced by G4GIDI_target::getChannelIDs(), G4GIDI_target::getNumberOfChannels(), and G4GIDI_target::init().

                                                                                     {

    return( tpia_target_heated_numberOfChannels( smr, target->baseHeatedTarget ) );
}

int tpia_target_numberOfProductionChannels	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 324 of file tpia_target.cc.

References tpia_target_s::baseHeatedTarget, and tpia_target_heated_numberOfProductionChannels().

Referenced by G4GIDI_target::getNumberOfProductionChannels(), and G4GIDI_target::getProductionChannelIDs().

                                                                                               {

    return( tpia_target_heated_numberOfProductionChannels( smr, target->baseHeatedTarget ) );
}

int tpia_target_read	(	statusMessageReporting *	smr,
		tpia_target *	target,
		const char *	fileName
	)

Definition at line 143 of file tpia_target.cc.

References tpia_target_s::absPath, xData_element_s::attributes, tpia_target_s::attributes, tpia_target_s::baseHeatedTarget, tpia_target_heated_info_s::contents, tpia_target_heated_info_s::heatedTarget, tpia_target_s::heatedTargets, xData_element_s::name, tpia_target_s::nHeatedTargets, tpia_target_heated_info_s::ordinal, tpia_target_heated_info_s::path, tpia_target_s::path, tpia_target_s::projectileID, tpia_target_s::readHeatedTargets, smr_isOk(), smr_setMessageError(), status, tpia_target_s::targetID, tpia_target_heated_info_s::temperature, tpia_misc_pointerToAttributeIfAllOk2, tpia_misc_setMessageError_Element(), tpia_particle_getInternalID(), tpia_target_initialize(), tpia_target_readHeatedTarget(), tpia_target_release(), xData_convertAttributeToDouble(), xData_copyAttributionList(), xData_getAttributesValueInElement(), xData_getDocumentsElement(), xData_getFirstElement(), xData_getNextElement(), xData_malloc2, xData_parseFree(), xData_parseReadFile(), xDataMisc_allocateCopyString2, and xDataMisc_getAbsPath().

Referenced by tpia_target_createRead(), and tpia_target_readFromMap().

                                                                                               {
/*
*   If a target has already been read into this target, user must have called tpia_target_release before calling this routine.
*   Otherwise, there will be memory leaks.
*/
    xData_document *doc;
    xData_element *element, *child;
    int i, iHeated, nHeated = 0, status = 1;
    double temperature;
    //fix for gcc4.6 warings 110602
    //char *pReturnValue, *name;
    char *name;
    char const *contents;

    tpia_target_initialize( smr, target );
    if( ( target->path = xDataMisc_allocateCopyString2( smr, fileName, "path" ) ) == NULL ) return( status );
    if( ( target->absPath = xDataMisc_getAbsPath( smr, fileName ) ) == NULL ) return( _tpia_target_releaseAndReturnOne( smr, target ) );
    if( ( doc = xData_parseReadFile( smr, fileName, NULL, NULL ) ) == NULL ) return( _tpia_target_releaseAndReturnOne( smr, target ) );
    element = xData_getDocumentsElement( doc );
    if( strcmp( element->name, "xTarget" ) != 0 ) {
        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "input file's top element must be xTarget and not %s", element->name ); }
    else {
        //pReturnValue = ( xData_copyAttributionList( smr, &(target->attributes), &(element->attributes) ) ) ? NULL : target->path;
        //fix for gcc4.6 warings 110602
        xData_copyAttributionList( smr, &(target->attributes),&(element->attributes) );
        name = tpia_misc_pointerToAttributeIfAllOk2( smr, element, 1, &(target->attributes), "projectile" );
        if( smr_isOk( smr ) ) target->projectileID = tpia_particle_getInternalID( smr, name );
        if( smr_isOk( smr ) && ( name = tpia_misc_pointerToAttributeIfAllOk2( smr, element, 1, &(target->attributes), "target" ) ) != NULL ) {
            if( smr_isOk( smr ) && ( target->targetID = tpia_particle_getInternalID( smr, name ) ) != NULL ) {
                status = 0;
                for( nHeated = 0, child = xData_getFirstElement( element ); child != NULL; nHeated++, child = xData_getNextElement( child ) ) {
                    if( strcmp( child->name, "target" ) != 0 ) {
                        tpia_misc_setMessageError_Element( smr, NULL, element, __FILE__, __LINE__, 1, "element can only have target sub-elements%s", 
                                element->name );
                        status = 1;
                        break;
                    }
                }
                if( status == 0 ) {
                    //if( ( target->heatedTargets = xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info ), 1, "heatedTargets" ) ) == NULL ) {
                    if( ( target->heatedTargets = (tpia_target_heated_info*) xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info ), 1, "heatedTargets" ) ) == NULL ) {
                        status = 1; }
                    else {
                        //if( ( target->readHeatedTargets = xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info * ), 1, "heatedTargets" ) ) == NULL ) 
                        if( ( target->readHeatedTargets = (tpia_target_heated_info**) xData_malloc2( smr, nHeated * sizeof( tpia_target_heated_info * ), 1, "heatedTargets" ) ) == NULL ) 
                                status = 1;
                    }
                    for( nHeated = 0, child = xData_getFirstElement( element ); ( status == 0 ) && ( child != NULL ); 
                            nHeated++, child = xData_getNextElement( child ) ) {
                        if( ( i = xData_convertAttributeToDouble( smr, child, "temperature", &temperature ) ) != 0 ) {
                            if( i > 0 ) smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "target does not have a temperature attribute" );
                            status = 1;
                            break;
                        }
                        for( iHeated = 0; iHeated < nHeated; iHeated++ ) if( target->heatedTargets[iHeated].temperature > temperature ) break;
                        if( iHeated < nHeated ) for( i = nHeated; i >= iHeated; i-- ) target->heatedTargets[i+1] = target->heatedTargets[i];
                        target->heatedTargets[iHeated].temperature = temperature;
                        target->heatedTargets[iHeated].path = NULL;
                        target->heatedTargets[iHeated].contents = NULL;
                        target->heatedTargets[iHeated].heatedTarget = NULL;
                        if( ( contents = xData_getAttributesValueInElement( child, "contents" ) ) != NULL ) {
                            if( ( target->heatedTargets[iHeated].contents = xDataMisc_allocateCopyString2( smr, contents, "contents" ) ) == NULL ) {
                                status = 1;
                                break;
                            }
                        }
                        if( ( contents = xData_getAttributesValueInElement( child, "file" ) ) == NULL ) {
                            status = 1;
                            break;
                        }
                        //if((target->heatedTargets[iHeated].path = xData_malloc2(smr, strlen( target->absPath ) + strlen( contents ) + 2, 0, "path")) == NULL) {
                        if((target->heatedTargets[iHeated].path = (char*) xData_malloc2(smr, strlen( target->absPath ) + strlen( contents ) + 2, 0, "path")) == NULL) {
                            status = 1;
                            break;
                        }
                        strcpy( target->heatedTargets[iHeated].path, target->absPath );
                        *strrchr( target->heatedTargets[iHeated].path, '/' ) = 0;
                        strcat( target->heatedTargets[iHeated].path, "/" );
                        strcat( target->heatedTargets[iHeated].path, contents );
                        target->nHeatedTargets++;
                    }
                }
            }
        }
    }
    xData_parseFree( smr, doc );
    if( status == 0 ) {
        for( i = 0; i < nHeated; i++ ) target->heatedTargets[i].ordinal = i;
        for( i = 0; i < nHeated; i++ ) if( target->heatedTargets[i].contents == NULL ) break;
        if( i == nHeated ) i = 0;                                           /* All heated targets are crossSection only. */
        if( tpia_target_readHeatedTarget( smr, target, i, 0 ) == 0 ) {
            target->baseHeatedTarget = target->heatedTargets[i].heatedTarget; }
        else {
            tpia_target_release( NULL, target );
            status = 1;
        } }
    else {
        tpia_target_release( smr, target );
    }
    return( status );
}

int tpia_target_readFromMap	(	statusMessageReporting *	smr,
		tpia_target *	target,
		tpia_map *	map,
		const char *	evaluation,
		const char *	projectileName,
		const char *	targetName
	)

Definition at line 82 of file tpia_target.cc.

References tpia_map_findTarget(), and tpia_target_read().

                                 {

    char *targetPath;

    if( ( targetPath = tpia_map_findTarget( smr, map, evaluation, projectileName, targetName ) ) == NULL ) return( 1 );
    return( tpia_target_read( smr, target, targetPath ) );
}

int tpia_target_readHeatedTarget	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		int	checkElememtsForAccess
	)

Definition at line 266 of file tpia_target.cc.

References tpia_target_heated_info_s::heatedTarget, tpia_target_s::heatedTargets, tpia_target_s::nHeatedTargets, tpia_target_s::nReadHeatedTargets, tpia_target_heated_s::ordinal, tpia_target_heated_info_s::ordinal, tpia_target_heated_info_s::path, tpia_target_s::readHeatedTargets, smr_setMessageError(), tpia_target_heated_info_s::temperature, and tpia_target_heated_createRead().

Referenced by G4GIDI_target::readTemperature(), tpia_target_getHeatedTargetAtIndex_ReadIfNeeded(), and tpia_target_read().

                                                                                                                            {

    int i;

    if( ( index < 0 ) || ( index >= target->nHeatedTargets ) ) {
        smr_setMessageError( smr, NULL, __FILE__, __LINE__, 1, "temperature index = %d out of range (0 <= index < %d", index, target->nHeatedTargets );
        return( -1 );
    }
    if( target->heatedTargets[index].heatedTarget != NULL ) return( 1 );
    target->heatedTargets[index].heatedTarget = tpia_target_heated_createRead( smr, target->heatedTargets[index].path, checkElememtsForAccess );
    if( target->heatedTargets[index].heatedTarget != NULL ) {
        target->heatedTargets[index].heatedTarget->ordinal = target->heatedTargets[index].ordinal;
        for( i = target->nReadHeatedTargets; i > 0; i-- ) {
            if( target->readHeatedTargets[i-1]->temperature < target->heatedTargets[index].temperature ) break;
            target->readHeatedTargets[i] = target->readHeatedTargets[i-1];
        }
        target->readHeatedTargets[i] = &(target->heatedTargets[i]);
        target->nReadHeatedTargets++;
    }
    return( ( target->heatedTargets[index].heatedTarget == NULL ? -1 : 0 ) );
}

int tpia_target_release	(	statusMessageReporting *	smr,
		tpia_target *	target
	)

Definition at line 117 of file tpia_target.cc.

References tpia_target_s::absPath, tpia_target_s::attributes, tpia_target_heated_info_s::contents, tpia_target_heated_info_s::heatedTarget, tpia_target_s::heatedTargets, tpia_target_s::nHeatedTargets, tpia_target_heated_info_s::path, tpia_target_s::path, tpia_target_s::readHeatedTargets, tpia_target_heated_free(), tpia_target_initialize(), xData_free(), and xData_releaseAttributionList().

Referenced by tpia_target_free(), and tpia_target_read().

                                                                            {

    int i;

    //target->path = xData_free( smr, target->path );
    target->path = (char*) xData_free( smr, target->path );
    //target->absPath = xData_free( smr, target->absPath );
    target->absPath = (char*) xData_free( smr, target->absPath );
    xData_releaseAttributionList( smr, &(target->attributes) );
    for( i = 0; i < target->nHeatedTargets; i++ ) {
        //target->heatedTargets[i].path = xData_free( smr, target->heatedTargets[i].path );
        target->heatedTargets[i].path = (char*) xData_free( smr, target->heatedTargets[i].path );
        //target->heatedTargets[i].contents = xData_free( smr, target->heatedTargets[i].contents );
        target->heatedTargets[i].contents = (char*) xData_free( smr, target->heatedTargets[i].contents );
            if( target->heatedTargets[i].heatedTarget != NULL ) tpia_target_heated_free( smr, target->heatedTargets[i].heatedTarget );
    }
    //target->heatedTargets = xData_free( smr, target->heatedTargets );
    target->heatedTargets = (tpia_target_heated_info*) xData_free( smr, target->heatedTargets );
    //target->readHeatedTargets = xData_free( smr, target->readHeatedTargets );
    target->readHeatedTargets = (tpia_target_heated_info**) xData_free( smr, target->readHeatedTargets );
    tpia_target_initialize( smr, target );
    return( 0 );
}

int tpia_target_sampleIndexChannelProductsAtE	(	statusMessageReporting *	smr,
		tpia_target *	target,
		int	index,
		double	T,
		tpia_decaySamplingInfo *	decaySamplingInfo,
		int	nProductData,
		tpia_productOutgoingData *	productData
	)

Definition at line 404 of file tpia_target.cc.

References tpia_target_s::baseHeatedTarget, and tpia_target_heated_sampleIndexChannelProductsAtE().

                                                                                                             {

    return( tpia_target_heated_sampleIndexChannelProductsAtE( smr, target->baseHeatedTarget, index, decaySamplingInfo,
        nProductData, productData ) );
}

Variable Documentation

DLL_LEND const char * tpia_productGenre_NBody_angular_energy

Definition at line 60 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char * tpia_productGenre_NBody_Legendre

Definition at line 59 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char * tpia_productGenre_NBody_pairProduction

Definition at line 62 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char * tpia_productGenre_NBody_uncorrelate_Legendre

Definition at line 61 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char * tpia_productGenre_twoBody_angular

Definition at line 57 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char * tpia_productGenre_twoBody_formFactor

Definition at line 58 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char* tpia_productGenre_unknown

Definition at line 56 of file tpia_product.cc.

Referenced by tpia_product_getFromElement().

DLL_LEND const char* tpia_samplingMethods_constant

Definition at line 45 of file tpia_samplingMethods.cc.

DLL_LEND const char * tpia_samplingMethods_linear

Definition at line 46 of file tpia_samplingMethods.cc.

Referenced by tpia_samplingMethods_initialize().