12 #if defined __cplusplus
27 return( angularEnergy );
68 char const *nativeData;
73 if( strcmp( nativeData,
"pointwise" ) == 0 ) {
75 else if( strcmp( nativeData,
"linear" ) == 0 ) {
94 double y,
norm, energyInFactor;
95 char const *energyUnit, *energyOutProbabilityUnits[2] = {
"MeV",
"1/MeV" };
117 if( ( pdfOfMuGivenE->
Ws = (
double *)
smr_malloc2( smr, V_W_XYs->
length *
sizeof(
double ), 1,
"pdfOfMuGivenE->Ws" ) ) == 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., 1
e-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;
134 for( iW = 0; iW < W_XYs->
length; iW++ ) {
135 XYs = &(W_XYs->
XYs[iW]);
144 pdfOfEpGivenEAndMu2->Ws[iW] = XYs->
value;
146 pdfOfEpGivenEAndMu2->numberOfWs++;
150 pdfOfMuGivenE->
Ws[iV] = energyInFactor * W_XYs->
value;
180 decaySamplingInfo->
frame = angularEnergy->
frame;
184 #if defined __cplusplus