11 #if defined __cplusplus
26 return( energyAngular );
67 char const *nativeData;
72 if( strcmp( nativeData,
"KalbachMann" ) == 0 ) {
74 else if( strcmp( nativeData,
"linear" ) == 0 ) {
92 double y,
norm, energyInFactor, energyOutFactor;
93 char const *energyUnit, *multiplicityProbabilityUnits[2] = {
"",
"1/MeV" };
118 if( ( pdfOfEpGivenE->
Ws = (
double *)
smr_malloc2( smr, V_W_XYs->
length *
sizeof(
double ), 1,
"pdfOfEpGivenE->Ws" ) ) == NULL )
goto err;
132 for( iV = 0; iV < V_W_XYs->
length; iV++ ) {
133 W_XYs = &(V_W_XYs->
W_XYs[iV]);
134 pdfOfMuGivenEAndEp2 = &(pdfOfMuGivenEAndEp[iV]);
135 pdfOfMuGivenEAndEp2->interpolationWY = interpolationWY;
136 pdfOfMuGivenEAndEp2->interpolationXY = interpolationXY;
137 if( ( pdfXY2 =
ptwXY_new( interpolationWY, NULL, 2., 1
e-6, W_XYs->
length, 10, &status, 0 ) ) == NULL )
goto errA;
138 if( ( pdfOfMuGivenEAndEp2->Ws = (
double *)
smr_malloc2( smr, W_XYs->
length *
sizeof(
double ), 1,
"pdfOfMuGivenEAndEp2->Ws" ) ) == NULL )
goto err;
140 for( iW = 0; iW < W_XYs->
length; iW++ ) {
141 XYs = &(W_XYs->
XYs[iW]);
150 pdfOfMuGivenEAndEp2->Ws[iW] = energyOutFactor * XYs->
value;
152 pdfOfMuGivenEAndEp2->numberOfWs++;
156 pdfOfEpGivenE->
Ws[iV] = energyInFactor * W_XYs->
value;
187 Ep = decaySamplingInfo->
mu;
188 decaySamplingInfo->
mu = decaySamplingInfo->
Ep;
189 decaySamplingInfo->
Ep = Ep;
190 decaySamplingInfo->
frame = energyAngular->
frame;
195 #if defined __cplusplus