Geant4 Cross Reference |
1 /* 1 2 # <<BEGIN-copyright>> 3 # <<END-copyright>> 4 */ 5 #include <string.h> 6 7 #include "MCGIDI_fromTOM.h" 8 #include "MCGIDI.h" 9 #include "MCGIDI_misc.h" 10 11 #if defined __cplusplus 12 namespace GIDI { 13 using namespace GIDI; 14 #endif 15 16 static int MCGIDI_energyAngular_linear_parseFr 17 /* 18 ********************************************** 19 */ 20 MCGIDI_energyAngular *MCGIDI_energyAngular_new 21 22 MCGIDI_energyAngular *energyAngular; 23 24 if( ( energyAngular = (MCGIDI_energyAngula 25 if( MCGIDI_energyAngular_initialize( smr, 26 return( energyAngular ); 27 } 28 /* 29 ********************************************** 30 */ 31 int MCGIDI_energyAngular_initialize( statusMes 32 33 memset( energyAngular, 0, sizeof( MCGIDI_e 34 return( 0 ); 35 } 36 /* 37 ********************************************** 38 */ 39 MCGIDI_energyAngular *MCGIDI_energyAngular_fre 40 41 MCGIDI_energyAngular_release( smr, energyA 42 smr_freeMemory( (void **) &energyAngular ) 43 return( NULL ); 44 } 45 /* 46 ********************************************** 47 */ 48 int MCGIDI_energyAngular_release( statusMessag 49 50 int i; 51 52 for( i = 0; i < energyAngular->pdfOfEpGive 53 MCGIDI_sampling_pdfsOfXGivenW_release( 54 } 55 smr_freeMemory( (void **) &(energyAngular- 56 MCGIDI_sampling_pdfsOfXGivenW_release( smr 57 MCGIDI_energyAngular_initialize( smr, ener 58 59 return( 0 ); 60 } 61 /* 62 ********************************************** 63 */ 64 int MCGIDI_energyAngular_parseFromTOM( statusM 65 66 xDataTOM_element *energyAngularElement; 67 char const *nativeData; 68 69 if( ( energyAngularElement = xDataTOME_get 70 71 if( ( nativeData = xDataTOM_getAttributesV 72 if( strcmp( nativeData, "KalbachMann" ) == 73 return( MCGIDI_KalbachMann_parseFromTO 74 else if( strcmp( nativeData, "linear" ) == 75 return( MCGIDI_energyAngular_linear_pa 76 else { 77 smr_setReportError2( smr, smr_unknownI 78 goto err; 79 } 80 81 return( 0 ); 82 83 err: 84 return( 1 ); 85 } 86 /* 87 ********************************************** 88 */ 89 static int MCGIDI_energyAngular_linear_parseFr 90 91 int iV, iW; 92 double y, norm, energyInFactor, energyOutF 93 char const *energyUnit, *multiplicityProba 94 xDataTOM_element *linear; 95 ptwXY_interpolation interpolationXY, inter 96 xDataTOM_XYs *XYs; 97 xDataTOM_W_XYs *W_XYs; 98 xDataTOM_V_W_XYs *V_W_XYs; 99 MCGIDI_pdfsOfXGivenW *pdfOfEpGivenE, *pdfO 100 MCGIDI_energyAngular *energyAngular = NULL 101 ptwXYPoints *pdfXY1 = NULL, *pdfXY2 = NULL 102 nfu_status status; 103 104 if( ( linear = xDataTOME_getOneElementByNa 105 106 if( MCGIDI_fromTOM_interpolation( smr, lin 107 if( MCGIDI_fromTOM_interpolation( smr, lin 108 if( MCGIDI_fromTOM_interpolation( smr, lin 109 110 if( ( energyAngular = MCGIDI_energyAngular 111 if( ( energyAngular->frame = MCGIDI_misc_g 112 113 pdfOfEpGivenE = &(energyAngular->pdfOfEpGi 114 pdfOfEpGivenE->interpolationWY = interpola 115 pdfOfEpGivenE->interpolationXY = interpola 116 117 if( ( V_W_XYs = (xDataTOM_V_W_XYs *) xData 118 if( ( pdfOfEpGivenE->Ws = (double *) smr_m 119 if( ( pdfOfEpGivenE->dist = (MCGIDI_pdfOfX 120 if( ( pdfOfMuGivenEAndEp = (MCGIDI_pdfsOfX 121 122 energyUnit = xDataTOM_subAxes_getUnit( smr 123 if( !smr_isOk( smr ) ) goto err; 124 energyInFactor = MCGIDI_misc_getUnitConver 125 if( !smr_isOk( smr ) ) goto err; 126 127 energyUnit = xDataTOM_subAxes_getUnit( smr 128 if( !smr_isOk( smr ) ) goto err; 129 energyOutFactor = MCGIDI_misc_getUnitConve 130 if( !smr_isOk( smr ) ) goto err; 131 132 for( iV = 0; iV < V_W_XYs->length; iV++ ) 133 W_XYs = &(V_W_XYs->W_XYs[iV]); 134 pdfOfMuGivenEAndEp2 = &(pdfOfMuGivenEA 135 pdfOfMuGivenEAndEp2->interpolationWY = 136 pdfOfMuGivenEAndEp2->interpolationXY = 137 if( ( pdfXY2 = ptwXY_new( interpolatio 138 if( ( pdfOfMuGivenEAndEp2->Ws = (doubl 139 if( ( pdfOfMuGivenEAndEp2->dist = (MCG 140 for( iW = 0; iW < W_XYs->length; iW++ 141 XYs = &(W_XYs->XYs[iW]); 142 if( ( pdfXY1 = MCGIDI_misc_dataFr 143 y = ptwXY_integrateDomain( pdfXY1, 144 if( ( status = ptwXY_setValueAtX( 145 146 if( y == 0 ) { 147 if( ( status = ptwXY_add_doubl 148 } 149 pdfOfMuGivenEAndEp2->Ws[iW] = ener 150 if( MCGIDI_fromTOM_pdfOfX( smr, pd 151 pdfOfMuGivenEAndEp2->numberOfWs++; 152 153 pdfXY1 = ptwXY_free( pdfXY1 ); 154 } 155 pdfOfEpGivenE->Ws[iV] = energyInFactor 156 if( MCGIDI_fromTOM_pdfOfX( smr, pdfXY2 157 pdfOfEpGivenE->numberOfWs++; 158 159 pdfXY2 = ptwXY_free( pdfXY2 ); 160 } 161 energyAngular->pdfOfMuGivenEAndEp = pdfOfM 162 distribution->energyAngular = energyAngula 163 distribution->type = MCGIDI_distributionTy 164 165 return( 0 ); 166 167 errA: 168 smr_setReportError2( smr, smr_unknownID, 1 169 err: 170 if( pdfXY1 != NULL ) ptwXY_free( pdfXY1 ); 171 if( pdfXY2 != NULL ) ptwXY_free( pdfXY2 ); 172 if( energyAngular != NULL ) MCGIDI_energyA 173 /* ????????? Need to free pdfOfMuGivenEAndEp, 174 return( 1 ); 175 } 176 /* 177 ********************************************** 178 */ 179 int MCGIDI_energyAngular_sampleDistribution( s 180 MCGIDI_decaySamplingInfo *decaySamplin 181 182 double Ep; 183 MCGIDI_energyAngular *energyAngular = dist 184 185 MCGIDI_sampling_doubleDistribution( smr, & 186 Ep = decaySamplingInfo->mu; 187 decaySamplingInfo->mu = decaySamplingInfo- 188 decaySamplingInfo->Ep = Ep; 189 decaySamplingInfo->frame = energyAngular-> 190 191 return( 0 ); 192 } 193 194 #if defined __cplusplus 195 } 196 #endif 197 198