8 #include "MCGIDI_fromTOM.h"
9 #include "MCGIDI_misc.h"
10 #include "MCGIDI_private.h"
12 #if defined __cplusplus
23 MCGIDI_angularEnergy *angularEnergy;
25 if( ( angularEnergy = (MCGIDI_angularEnergy *) smr_malloc2( smr,
sizeof( MCGIDI_angularEnergy ), 0,
"angularEnergy" ) ) == NULL )
return( NULL );
27 return( angularEnergy );
34 memset( angularEnergy, 0,
sizeof( MCGIDI_angularEnergy ) );
53 for( i = 0; i < angularEnergy->pdfOfMuGivenE.numberOfWs; i++ ) {
67 xDataTOM_element *angularEnergyElement, *pointwise = NULL;
68 char const *nativeData;
73 if( strcmp( nativeData,
"pointwise" ) == 0 ) {
75 else if( strcmp( nativeData,
"linear" ) == 0 ) {
78 smr_setReportError2( smr, smr_unknownID, 1,
"angularEnergy nativeData = '%s' not supported", nativeData );
94 double y, norm, energyInFactor;
95 char const *energyUnit, *energyOutProbabilityUnits[2] = {
"MeV",
"1/MeV" };
96 MCGIDI_angularEnergy *angularEnergy = NULL;
97 ptwXY_interpolation interpolationXY, interpolationWY, interpolationVY;
99 xDataTOM_W_XYs *W_XYs;
100 xDataTOM_V_W_XYs *V_W_XYs;
101 MCGIDI_pdfsOfXGivenW *pdfOfMuGivenE, *pdfOfEpGivenEAndMu = NULL, *pdfOfEpGivenEAndMu2 = NULL;
102 ptwXYPoints *pdfXY1 = NULL, *pdfXY2 = NULL;
112 pdfOfMuGivenE = &(angularEnergy->pdfOfMuGivenE);
113 pdfOfMuGivenE->interpolationWY = interpolationVY;
114 pdfOfMuGivenE->interpolationXY = interpolationWY;
116 if( ( V_W_XYs = (xDataTOM_V_W_XYs *)
xDataTOME_getXDataIfID( smr, pointwise,
"V_W_XYs" ) ) == NULL )
goto err;
117 if( ( pdfOfMuGivenE->Ws = (
double *) smr_malloc2( smr, V_W_XYs->length *
sizeof(
double ), 1,
"pdfOfMuGivenE->Ws" ) ) == NULL )
goto err;
118 if( ( pdfOfMuGivenE->dist = (MCGIDI_pdfOfX *) smr_malloc2( smr, V_W_XYs->length *
sizeof( MCGIDI_pdfOfX ), 0,
"pdfOfMuGivenE->dist" ) ) == NULL )
goto err;
119 if( ( pdfOfEpGivenEAndMu = (MCGIDI_pdfsOfXGivenW *) smr_malloc2( smr, V_W_XYs->length *
sizeof( MCGIDI_pdfsOfXGivenW ), 1,
"pdfOfEpGivenEAndMu" ) ) == NULL )
goto err;
126 for( iV = 0; iV < V_W_XYs->length; iV++ ) {
127 W_XYs = &(V_W_XYs->W_XYs[iV]);
128 pdfOfEpGivenEAndMu2 = &(pdfOfEpGivenEAndMu[iV]);
129 pdfOfEpGivenEAndMu2->interpolationWY = interpolationWY;
130 pdfOfEpGivenEAndMu2->interpolationXY = interpolationXY;
131 if( ( pdfXY2 =
ptwXY_new( interpolationWY, NULL, 2., 1e-6, W_XYs->length, 10, &status, 0 ) ) == NULL )
goto errA;
132 if( ( pdfOfEpGivenEAndMu2->Ws = (
double *) smr_malloc2( smr, W_XYs->length *
sizeof(
double ), 1,
"pdfOfEpGivenEAndMu2->Ws" ) ) == NULL )
goto err;
133 if( ( pdfOfEpGivenEAndMu2->dist = (MCGIDI_pdfOfX *) smr_malloc2( smr, W_XYs->length *
sizeof( MCGIDI_pdfOfX ), 0,
"pdfOfEpGivenEAndMu2->dist" ) ) == NULL )
goto err;
134 for( iW = 0; iW < W_XYs->length; iW++ ) {
135 XYs = &(W_XYs->XYs[iW]);
138 if( ( status =
ptwXY_setValueAtX( pdfXY2, XYs->value, y ) ) != nfu_Okay )
goto errA;
139 if( status != nfu_Okay )
goto errA;
144 pdfOfEpGivenEAndMu2->Ws[iW] = XYs->value;
146 pdfOfEpGivenEAndMu2->numberOfWs++;
150 pdfOfMuGivenE->Ws[iV] = energyInFactor * W_XYs->value;
152 pdfOfMuGivenE->numberOfWs++;
157 angularEnergy->pdfOfEpGivenEAndMu = pdfOfEpGivenEAndMu;
158 distribution->angularEnergy = angularEnergy;
159 distribution->type = MCGIDI_distributionType_angularEnergy_e;
164 smr_setReportError2( smr, smr_unknownID, 1,
"ptwXY_integrateDomain err = %d: %s\n", status,
nfu_statusMessage( status ) );
176 MCGIDI_decaySamplingInfo *decaySamplingInfo ) {
180 decaySamplingInfo->frame = angularEnergy->frame;
184 #if defined __cplusplus
int MCGIDI_angularEnergy_parseFromTOM(statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution)
xDataTOM_element * xDataTOME_getOneElementByName(statusMessageReporting *smr, xDataTOM_element *element, char const *name, int required)
int MCGIDI_angularEnergy_initialize(statusMessageReporting *, MCGIDI_angularEnergy *angularEnergy)
static int MCGIDI_angularEnergy_parsePointwiseFromTOM(statusMessageReporting *smr, xDataTOM_element *pointwise, MCGIDI_distribution *distribution)
MCGIDI_angularEnergy * MCGIDI_angularEnergy_new(statusMessageReporting *smr)
nfu_status ptwXY_add_double(ptwXYPoints *ptwXY, double value)
ptwXYPoints * ptwXY_free(ptwXYPoints *ptwXY)
nfu_status ptwXY_setValueAtX(ptwXYPoints *ptwXY, double x, double y)
int MCGIDI_sampling_pdfsOfXGivenW_release(statusMessageReporting *smr, MCGIDI_pdfsOfXGivenW *dists)
ptwXYPoints * MCGIDI_misc_dataFromXYs2ptwXYPointsInUnitsOf(statusMessageReporting *smr, xDataTOM_XYs *XYs, ptwXY_interpolation interpolation, char const *toUnits[2])
char const * xDataTOM_subAxes_getUnit(statusMessageReporting *smr, xDataTOM_subAxes *subAxes, int index)
int MCGIDI_angularEnergy_sampleDistribution(statusMessageReporting *smr, MCGIDI_angularEnergy *angularEnergy, MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo)
int MCGIDI_sampling_doubleDistribution(statusMessageReporting *smr, MCGIDI_pdfsOfXGivenW *pdfOfWGivenV, MCGIDI_pdfsOfXGivenW *pdfOfXGivenVAndW, MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo)
int MCGIDI_fromTOM_pdfOfX(statusMessageReporting *smr, ptwXYPoints *pdfXY, MCGIDI_pdfOfX *dist, double *norm)
int smr_isOk(statusMessageReporting *smr)
void * xDataTOME_getXDataIfID(statusMessageReporting *smr, xDataTOM_element *TE, char const *ID)
void * smr_freeMemory(void **p)
double MCGIDI_misc_getUnitConversionFactor(statusMessageReporting *smr, char const *fromUnit, char const *toUnit)
char const * xDataTOM_getAttributesValueInElement(xDataTOM_element *element, char const *name)
int MCGIDI_angularEnergy_release(statusMessageReporting *smr, MCGIDI_angularEnergy *angularEnergy)
MCGIDI_angularEnergy * MCGIDI_angularEnergy_free(statusMessageReporting *smr, MCGIDI_angularEnergy *angularEnergy)
const char * nfu_statusMessage(nfu_status status)
enum xDataTOM_frame MCGIDI_misc_getProductFrame(statusMessageReporting *smr, xDataTOM_element *frameElement)
ptwXYPoints * ptwXY_new(ptwXY_interpolation interpolation, ptwXY_interpolationOtherInfo const *interpolationOtherInfo, double biSectionMax, double accuracy, int64_t primarySize, int64_t secondarySize, nfu_status *status, int userFlag)
double ptwXY_integrateDomain(ptwXYPoints *ptwXY, nfu_status *status)
int MCGIDI_fromTOM_interpolation(statusMessageReporting *smr, xDataTOM_element *element, int index, ptwXY_interpolation *interpolation)