| 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