Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/lend/include/MCGIDI.h

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 /*
  2 # <<BEGIN-copyright>>
  3 # <<END-copyright>>
  4 */
  5 #ifndef MCGIDI_h_included
  6 #define MCGIDI_h_included
  7 
  8 #define MCGIDI_VERSION_MAJOR 1
  9 #define MCGIDI_VERSION_MINOR 0
 10 #define MCGIDI_VERSION_PATCHLEVEL 0
 11 
 12 #include <GIDI_settings.hh>
 13 #include <map>
 14 #include <vector>
 15 
 16 #include <statusMessageReporting.h>
 17 #include <ptwXY.h>
 18 #include <xDataTOM.h>
 19 
 20 #include "MCGIDI_mass.h"
 21 #include "MCGIDI_map.h"
 22 
 23 /* Disable Effective C++ warnings in GIDI code. */
 24 #if __INTEL_COMPILER > 1399
 25 #pragma warning( disable:2021 )
 26 #pragma warning( disable:593 )
 27 #pragma warning( disable:111 )
 28 #elif __INTEL_COMPILER > 1199
 29 #pragma warning( disable:2304 )
 30 #endif
 31 
 32 #if defined __cplusplus
 33     extern "C" {
 34     namespace GIDI {
 35 #endif
 36 
 37 typedef struct MCGIDI_GammaBranching_s MCGIDI_GammaBranching;
 38 typedef struct MCGIDI_POP_s MCGIDI_POP;
 39 typedef struct MCGIDI_POPs_s MCGIDI_POPs;
 40 typedef struct MCGIDI_particle_s MCGIDI_particle;
 41 typedef struct MCGIDI_target_s MCGIDI_target;
 42 typedef struct MCGIDI_target_heated_info_s MCGIDI_target_heated_info;
 43 typedef struct MCGIDI_target_heated_sorted_s MCGIDI_target_heated_sorted;
 44 typedef struct MCGIDI_target_heated_s MCGIDI_target_heated;
 45 typedef struct MCGIDI_reaction_s MCGIDI_reaction;
 46 typedef struct MCGIDI_outputChannel_s MCGIDI_outputChannel;
 47 typedef struct MCGIDI_product_s MCGIDI_product;
 48 typedef struct MCGIDI_distribution_s MCGIDI_distribution;
 49 typedef struct MCGIDI_KalbachMann_s MCGIDI_KalbachMann;
 50 typedef struct MCGIDI_KalbachMann_ras_s MCGIDI_KalbachMann_ras;
 51 typedef struct MCGIDI_pdfOfX_s MCGIDI_pdfOfX;
 52 typedef struct MCGIDI_pdfsOfXGivenW_s MCGIDI_pdfsOfXGivenW;
 53 typedef struct MCGIDI_pdfsOfXGivenW_sampled_s MCGIDI_pdfsOfXGivenW_sampled;
 54 typedef struct MCGIDI_angular_s MCGIDI_angular;
 55 typedef struct MCGIDI_energyWeightedFunctional_s MCGIDI_energyWeightedFunctional;
 56 typedef struct MCGIDI_energyWeightedFunctionals_s MCGIDI_energyWeightedFunctionals;
 57 typedef struct MCGIDI_energyNBodyPhaseSpace_s MCGIDI_energyNBodyPhaseSpace;
 58 typedef struct MCGIDI_energy_s MCGIDI_energy;
 59 typedef struct MCGIDI_energyAngular_s MCGIDI_energyAngular;
 60 typedef struct MCGIDI_angularEnergy_s MCGIDI_angularEnergy;
 61 
 62 typedef struct MCGIDI_decaySamplingInfo_s MCGIDI_decaySamplingInfo;
 63 typedef struct MCGIDI_productsInfo_s MCGIDI_productsInfo;
 64 typedef struct MCGIDI_productInfo_s MCGIDI_productInfo;
 65 typedef struct MCGIDI_sampledProductsData_s MCGIDI_sampledProductsData;
 66 typedef struct MCGIDI_sampledProductsDatas_s MCGIDI_sampledProductsDatas;
 67 
 68 #if defined __cplusplus
 69     }
 70     }
 71 #endif
 72 
 73 enum MCGIDI_quantityLookupMode {      
 74     MCGIDI_quantityLookupMode_pointwise     /**< Pointwise data are used to determine a quantity's value an energy E. */,
 75     MCGIDI_quantityLookupMode_grouped       /**< Grouped data are used to determine a quantity's value an energy E. */
 76 };
 77 
 78 class MCGIDI_quantitiesLookupModes {
 79 
 80     private:
 81         int mProjectilesPOPID;
 82         double mProjectileEnergy;
 83         int mGroupIndex;
 84         double mProjectileEnergyForGroupIndex;
 85         double mTemperature;
 86         enum MCGIDI_quantityLookupMode mCrossSectionMode;
 87         enum MCGIDI_quantityLookupMode mMultiplicityMode;
 88 
 89     public:
 90         MCGIDI_quantitiesLookupModes( int projectilesPOPID );
 91         ~MCGIDI_quantitiesLookupModes( );
 92 
 93         inline double getProjectileEnergy( void ) const { return( mProjectileEnergy ); }
 94         void setProjectileEnergy( double e_in ) { mProjectileEnergy = e_in; }
 95 
 96         inline int getGroupIndex( void ) const { return( mGroupIndex ); }
 97         int setGroupIndex( GIDI_settings const &settings, bool encloseOutOfRange );
 98 
 99         inline double getTemperature( void ) const { return( mTemperature ); }
100         void setTemperature( double temperature ) { mTemperature = temperature; }
101 
102         enum MCGIDI_quantityLookupMode getMode( std::string const &quantity ) const;
103         enum MCGIDI_quantityLookupMode getCrossSectionMode( void ) const { return( mCrossSectionMode ); };
104         std::vector<std::string> getListOfLookupQuanities( ) const;
105         void setMode( std::string const &quantity, enum MCGIDI_quantityLookupMode mode );
106         void setCrossSectionMode( enum MCGIDI_quantityLookupMode mode ) { mCrossSectionMode = mode; };
107         void setModeAll( enum MCGIDI_quantityLookupMode mode );
108 };
109 
110 typedef struct MCGIDI_samplingMultiplicityBias_s MCGIDI_samplingMultiplicityBias;
111 
112 struct MCGIDI_samplingMultiplicityBias_s {
113     int PoPID;
114     double multiplicityFactor;
115 };
116 
117 class MCGIDI_samplingMethods {
118 
119     public:
120         MCGIDI_samplingMethods( );
121         ~MCGIDI_samplingMethods( );
122 };
123 
124 class MCGIDI_samplingSettings {
125 
126     private:                    // This is user input.
127         enum GIDI::xDataTOM_frame mWantFrame;
128         bool mWantVelocities;
129         double (*mRng)( void * );
130         void *mRngState;
131         std::vector<struct MCGIDI_samplingMultiplicityBias_s> mSamplingMultiplicityBiases;
132 
133     public:                     // Temporary variables used in MCGIDI sampling routines.
134         enum GIDI::xDataTOM_frame mGotFrame;
135         GIDI::MCGIDI_POP *mPoP;
136         double mMu;
137         double mEp;
138 
139     public:
140         MCGIDI_samplingSettings( enum GIDI::xDataTOM_frame frame, bool wantVelocities, double (*rng)( void * ), void *rngState );
141         ~MCGIDI_samplingSettings( );
142 
143         inline double getProductMultiplicityBias( int PoPID ) const {
144                 for( int i1 = 0; i1 < (int) mSamplingMultiplicityBiases.size( ); ++i1 ) {
145                     if( PoPID == mSamplingMultiplicityBiases[i1].PoPID ) return( mSamplingMultiplicityBiases[i1].multiplicityFactor );
146                 }
147                 return( 1. ); }
148         int setProductMultiplicityBias( GIDI::statusMessageReporting *smr, int PoPID, double fractor );
149 };
150 
151 #if defined __cplusplus
152     extern "C" {
153     namespace GIDI {
154 #endif
155 
156 enum MCGIDI_transportability {                          /**< This enum is used to give the transportability status for a particle in a reaction or target. */
157     MCGIDI_transportability_unknown,                    /**< Particle is not a product of this reaction or target. */
158     MCGIDI_transportability_none,                       /**< Particle is a product but has not distribution data. */
159     MCGIDI_transportability_partial,                    /**< Particle is a product and has some distribution data. */
160     MCGIDI_transportability_full };                     /**< Particle is a product and all needed distribution data. */
161 
162 #if defined __cplusplus
163     }
164     }
165 #endif
166 
167 typedef std::map<int, enum GIDI::MCGIDI_transportability> transportabilitiesMap;
168 
169 #if defined __cplusplus
170     extern "C" {
171     namespace GIDI {
172 #endif
173 
174 #define MCGIDI_crossSectionType_grouped 1
175 #define MCGIDI_crossSectionType_pointwise 2
176 
177 #define MCGIDI_nullReaction -10001
178 
179 #define MCGIDI_speedOfLight_cm_sec 2.99792458e10
180 #define MCGIDI_AMU2MeV 931.494028
181 
182 enum MCGIDI_reactionType { 
183     MCGIDI_reactionType_unknown_e,                      /* This should never happen. */
184     MCGIDI_reactionType_null_e,                         /* Only occurs when sampling with from grouped cross sections and the projectile is below threshold. */
185     MCGIDI_reactionType_elastic_e,                      /* A nuclear elastic reaction. */
186     MCGIDI_reactionType_scattering_e,                   /* A nuclear reaction where the projectile and target are products as well as gammas, 
187                                                             excluding reactions that are MCGIDI_reactionType_elastic_e and 
188                                                             MCGIDI_reactionType_nuclearLevelTransition_e. */
189     MCGIDI_reactionType_nuclearIsomerTransmutation_e,   /* A nuclear that changes N or Z and is not one of the others.*/
190     MCGIDI_reactionType_nuclearLevelTransition_e,       /* Reaction in which the residual is the same isotope as the target but in a
191                                                             different nuclear level. Mainly for meta-stables. */
192     MCGIDI_reactionType_capture_e,                      /* A nuclear capture reaction. */
193     MCGIDI_reactionType_fission_e,                      /* A nuclear fission reaction. */
194     MCGIDI_reactionType_sumOfRemainingOutputChannels_e, /* ENDF MT 5 reactions. */
195     MCGIDI_reactionType_atomic_e                        
196 };
197 
198 enum MCGIDI_channelGenre { MCGIDI_channelGenre_undefined_e, MCGIDI_channelGenre_twoBody_e, MCGIDI_channelGenre_uncorrelated_e, 
199     MCGIDI_channelGenre_sumOfRemaining_e, MCGIDI_channelGenre_twoBodyDecay_e, MCGIDI_channelGenre_uncorrelatedDecay_e };
200 
201 enum MCGIDI_productMultiplicityType { MCGIDI_productMultiplicityType_invalid_e, MCGIDI_productMultiplicityType_unknown_e, MCGIDI_productMultiplicityType_integer_e,
202     MCGIDI_productMultiplicityType_energyDependent_e, MCGIDI_productMultiplicityType_gammaBranching_e, MCGIDI_productMultiplicityType_mixed_e };
203 
204 enum MCGIDI_distributionType { MCGIDI_distributionType_none_e, MCGIDI_distributionType_unknown_e, MCGIDI_distributionType_angular_e, 
205     MCGIDI_distributionType_KalbachMann_e, MCGIDI_distributionType_uncorrelated_e, MCGIDI_distributionType_energyAngular_e, 
206     MCGIDI_distributionType_angularEnergy_e };
207 
208 enum MCGIDI_angularType { MCGIDI_angularType_isotropic, MCGIDI_angularType_recoil, MCGIDI_angularType_linear };
209 
210 enum MCGIDI_energyType { MCGIDI_energyType_unknown, MCGIDI_energyType_primaryGamma, MCGIDI_energyType_discreteGamma, 
211     MCGIDI_energyType_linear, MCGIDI_energyType_generalEvaporation, MCGIDI_energyType_simpleMaxwellianFission, MCGIDI_energyType_evaporation, 
212         MCGIDI_energyType_Watt, MCGIDI_energyType_MadlandNix, MCGIDI_energyType_NBodyPhaseSpace, MCGIDI_energyType_weightedFunctional };
213 
214 extern const char *MCGIDI_productGenre_unknown, *MCGIDI_productGenre_twoBody_angular, *MCGIDI_productGenre_twoBody_formFactor,
215     *MCGIDI_productGenre_NBody_angular_energy, *MCGIDI_productGenre_NBody_pairProduction;
216 
217 #define MCGIDI_particleLevel_continuum -1
218 #define MCGIDI_particleLevel_sum -2
219 
220 struct MCGIDI_GammaBranching_s {
221     MCGIDI_POP *finalLevel;
222     double probability;
223 };
224 
225 struct MCGIDI_POP_s {
226     MCGIDI_POP *next;
227     MCGIDI_POP *parent;
228     char *name;
229     int globalPoPsIndex;        /* Index of particle in the PoPs library if particle can be return to packages using */
230     int Z, A, level, m;         /* this library. Otherwise, -1. */
231     double mass_MeV;
232     double level_MeV;
233     int numberOfGammaBranchs;
234     MCGIDI_GammaBranching *gammas;
235 };
236 
237 struct MCGIDI_POPs_s {
238     int numberOfPOPs, size, increment;
239     MCGIDI_POP *first, *last, **sorted;
240 };
241 
242 struct MCGIDI_particle_s {
243     MCGIDI_particle *prior;
244     MCGIDI_particle *next;
245     int ordinal;
246     int Z, A, m;
247     double mass_MeV;
248     char *name;
249 };
250 
251 struct MCGIDI_decaySamplingInfo_s {
252     enum xDataTOM_frame frame;                  /* The frame the product data are in. */
253     int isVelocity;                             /* See struct MCGIDI_sampledProductsData_s for meaning. This is user input. */
254     double (*rng)( void * );                    /* User supplied rng. */
255     void *rngState;                             /* User supplied rng state. */
256     MCGIDI_POP *pop;                            /* pop for the sampled product. */
257     double mu;                                  /* mu = cos( theta ) for the sampled product. Frame is given by frame member. */
258     double Ep;                                  /* Energy of the product. Frame is given by frame member. */
259 };
260 
261 struct MCGIDI_productInfo_s {
262     int globalPoPsIndex;
263     enum MCGIDI_productMultiplicityType productMultiplicityType;
264     int multiplicity;
265     int transportable;
266 };
267 
268 struct MCGIDI_productsInfo_s {
269     int numberOfProducts;
270     int numberOfAllocatedProducts;
271     MCGIDI_productInfo *productInfo;
272 };
273 
274 struct MCGIDI_sampledProductsData_s {
275     int isVelocity;             /* If true, px_vx, py_vy and pz_vz are velocities otherwise momenta. */
276     MCGIDI_POP *pop;
277     double kineticEnergy;
278     double px_vx;
279     double py_vy;
280     double pz_vz;
281     int delayedNeutronIndex;
282     double delayedNeutronRate;
283     double birthTimeSec;        /* Some products, like delayed fission neutrons, are to appear (be born) later. */
284 };
285 
286 struct MCGIDI_sampledProductsDatas_s {
287     int numberOfProducts;
288     int numberAllocated;
289     int incrementSize;
290     MCGIDI_sampledProductsData *products;
291 };
292 
293 struct MCGIDI_pdfOfX_s {
294     int numberOfXs;
295     double *Xs;
296     double *pdf;
297     double *cdf;
298 };
299 
300 struct MCGIDI_pdfsOfXGivenW_s {
301     int numberOfWs;
302     ptwXY_interpolation interpolationWY, interpolationXY;
303     double *Ws;
304     MCGIDI_pdfOfX *dist;
305 };
306 
307 struct MCGIDI_pdfsOfXGivenW_sampled_s {
308     statusMessageReporting *smr;
309     ptwXY_interpolation interpolationWY, interpolationXY;
310     int iW, iX1, iX2;
311     double x, w, frac;
312 };
313 
314 struct MCGIDI_angular_s {
315     enum xDataTOM_frame frame;
316     enum MCGIDI_angularType type;
317     MCGIDI_angular *recoilProduct;
318     MCGIDI_pdfsOfXGivenW dists;
319     double projectileMass_MeV, targetMass_MeV, productMass_MeV, residualMass_MeV;
320 };
321 
322 struct MCGIDI_energyWeightedFunctional_s {
323     ptwXYPoints *weight;
324     MCGIDI_energy *energy;
325 };
326 
327 struct MCGIDI_energyWeightedFunctionals_s {
328     int numberOfWeights;
329     MCGIDI_energyWeightedFunctional weightedFunctional[4];      /* ??????????? Hardwired for no good reason. Will handle up to a (z,4n) reaction. */
330 };
331 
332 struct MCGIDI_energyNBodyPhaseSpace_s {
333     int numberOfProducts;
334     double mass, massFactor, e_inCOMFactor, Q_MeV;
335 };
336 
337 struct MCGIDI_energy_s {
338     enum xDataTOM_frame frame;
339     enum MCGIDI_energyType type;
340     double gammaEnergy_MeV;
341     double primaryGammaMassFactor;
342     double e_inCOMFactor;
343     MCGIDI_pdfsOfXGivenW dists;
344     double U;
345     ptwXYPoints *theta, *Watt_a, *Watt_b;
346     ptwXY_interpolation gInterpolation;
347     MCGIDI_pdfOfX g;
348     MCGIDI_energyWeightedFunctionals weightedFunctionals;
349     MCGIDI_energyNBodyPhaseSpace NBodyPhaseSpace;
350 };
351 
352 struct MCGIDI_energyAngular_s {
353     enum xDataTOM_frame frame;
354     MCGIDI_pdfsOfXGivenW pdfOfEpGivenE;
355     MCGIDI_pdfsOfXGivenW *pdfOfMuGivenEAndEp;   /* The number of MCGIDI_pdfsOfXGivenW allocated is given by pdfOfEpGivenE.numberOfWs. */
356 };
357 
358 struct MCGIDI_angularEnergy_s {
359     enum xDataTOM_frame frame;
360     MCGIDI_pdfsOfXGivenW pdfOfMuGivenE;
361     MCGIDI_pdfsOfXGivenW *pdfOfEpGivenEAndMu;   /* The number of MCGIDI_pdfsOfXGivenW allocated is given by pdfOfMuGivenE.numberOfWs. */
362 };
363 
364 struct MCGIDI_KalbachMann_ras_s {
365     double *rs;
366     double *as;
367 };
368 
369 struct MCGIDI_KalbachMann_s {
370     enum xDataTOM_frame frame;
371     double energyToMeVFactor, massFactor, Sa, Sb, Ma, mb;           /* Needed if a(E,E') is caluclated from the formula. */
372     MCGIDI_pdfsOfXGivenW dists;                             /* Sa currently not used. */
373     MCGIDI_KalbachMann_ras *ras;
374 };
375 
376 struct MCGIDI_distribution_s {
377     MCGIDI_product *product;
378     enum MCGIDI_distributionType type;
379     MCGIDI_angular *angular;                /* All distribution forms must have a frame member. */
380     MCGIDI_energy *energy;
381     MCGIDI_energyAngular *energyAngular;
382     MCGIDI_angularEnergy *angularEnergy;
383     MCGIDI_KalbachMann *KalbachMann;
384 };
385 
386 struct MCGIDI_outputChannel_s {
387     enum MCGIDI_channelGenre genre;
388     MCGIDI_reaction *reaction;              /* This is only used for output channels. */
389     MCGIDI_product *parent;                 /* This is only used for decay channels. */
390     int QIsFloat;
391     double Q;
392     int numberOfProducts;
393     MCGIDI_product *products;
394 };
395 
396 struct MCGIDI_product_s {
397     MCGIDI_POP *pop;
398     char *label;
399     MCGIDI_outputChannel *outputChannel;
400     int multiplicity;                                       /* If 0, the multiplicity is either 'energyDependent' or 'partialProduction'. */
401     int delayedNeutronIndex;
402     double delayedNeutronRate;
403     ptwXYPoints *multiplicityVsEnergy;
404     ptwXYPoints *norms;
405     int numberOfPiecewiseMultiplicities;
406     ptwXYPoints **piecewiseMultiplicities;
407     MCGIDI_distribution distribution;
408     MCGIDI_outputChannel decayChannel;
409 };
410 
411 struct MCGIDI_reaction_s {
412     MCGIDI_target_heated *target;
413     int ENDF_MT, ENDL_C, ENDL_S;
414     enum MCGIDI_reactionType reactionType;
415     char const *outputChannelStr;
416     xDataTOM_attributionList attributes;            /* Do not free, owned by attributes. */
417     int domainValuesPresent;                        /* True if cross section data defined so EMin and EMax are value. */
418     int thresholdGroupIndex;                        /* For grouped data, the group index where threshold starts. */
419     double thresholdGroupDomain;                    /* This is groupEnergy[thresholdGroupIndex+1] - EMin. */
420     double thresholdGroupedDeltaCrossSection;       /* The adjusted group cross section in group thresholdGroupIndex. */
421     double EMin, EMax, finalQ;                      /* BRB, EMin is used as threshold. However, some reactions, especially charged particle */
422     ptwXYPoints *crossSection;                      /* have effective thresholds much higher than EMin, may need to handle these differently??????? */
423     ptwXPoints *crossSectionGrouped;
424     MCGIDI_outputChannel outputChannel;
425     MCGIDI_productsInfo productsInfo;               /* See MCGIDI_reaction_ParseDetermineReactionProducts for description. */
426     transportabilitiesMap *transportabilities;
427 };
428 
429 struct MCGIDI_target_heated_s {
430     int ordinal;
431     char *path;            /* Partial path of input file. */
432     char *absPath;         /* Full absolute path of input file. */
433     MCGIDI_POPs pops;
434     MCGIDI_POP *projectilePOP;
435     MCGIDI_POP *targetPOP;
436     xDataTOM_attributionList attributes;
437     char *contents;
438     double temperature_MeV;
439     double EMin, EMax;
440     ptwXYPoints *crossSection;
441     ptwXPoints *crossSectionGrouped;
442     ptwXPoints *crossSectionGroupedForSampling;
443     int numberOfReactions;
444     MCGIDI_reaction *reactions;
445     transportabilitiesMap *transportabilities;
446 };
447 
448 struct MCGIDI_target_heated_info_s {
449     int ordinal;
450     double temperature;
451     char *path;                 /* Full path of input file. */
452     char *contents;
453     MCGIDI_target_heated *heatedTarget;
454 };
455 
456 struct MCGIDI_target_s {
457     char *path;                 /* Full path of input file. */
458     char *absPath;              /* Full absolute path of input file. */
459     MCGIDI_POP *projectilePOP;
460     MCGIDI_POP *targetPOP;
461     xDataTOM_attributionList attributes;
462     int nHeatedTargets, nReadHeatedTargets;
463     MCGIDI_target_heated *baseHeatedTarget;   /* The lowest temperature whose contents is "all" data, (e.g, not just "crossSection"). */
464     MCGIDI_target_heated_info *heatedTargets;         /* List of heated targets in order by temperature. */
465     MCGIDI_target_heated_info **readHeatedTargets;    /* List of "read in" heated targets in order by temperature. */
466 };
467 
468 char const *MCGIDI_version( void );
469 int MCGIDI_versionMajor( void );
470 int MCGIDI_versionMinor( void );
471 int MCGIDI_versionPatchLevel( void );
472 
473 /*
474 * Routines in MCGIDI_target.c
475 */
476 MCGIDI_target *MCGIDI_target_new( statusMessageReporting *smr );
477 int MCGIDI_target_initialize( statusMessageReporting *smr, MCGIDI_target *target );
478 MCGIDI_target *MCGIDI_target_newRead( statusMessageReporting *smr, const char *fileName );
479 int MCGIDI_target_readFromMapViaPoPIDs( statusMessageReporting *smr, MCGIDI_target *target, MCGIDI_map *map, const char *evaluation,
480         int projectile_PoPID, int target_PoPID );
481 int MCGIDI_target_readFromMap( statusMessageReporting *smr, MCGIDI_target *target, MCGIDI_map *map, const char *evaluation, const char *projectileName, 
482     const char *targetName );
483 MCGIDI_target *MCGIDI_target_newReadFromMapViaPoPIDs( statusMessageReporting *smr, MCGIDI_map *map, const char *evaluation,
484         int projectile_PoPID, int target_PoPID );
485 MCGIDI_target *MCGIDI_target_newReadFromMap( statusMessageReporting *smr, MCGIDI_map *map, const char *evaluation, const char *projectileName, 
486     const char *targetName );
487 MCGIDI_target *MCGIDI_target_free( statusMessageReporting *smr, MCGIDI_target *target );
488 int MCGIDI_target_release( statusMessageReporting *smr, MCGIDI_target *target );
489 int MCGIDI_target_read( statusMessageReporting *smr, MCGIDI_target *target, const char *fileName );
490 char const *MCGIDI_target_getAttributesValue( statusMessageReporting *smr, MCGIDI_target *target, char const *name );
491 int MCGIDI_target_getTemperatures( statusMessageReporting *smr, MCGIDI_target *target, double *temperatures );
492 int MCGIDI_target_readHeatedTarget( statusMessageReporting *smr, MCGIDI_target *target, int index );
493 MCGIDI_target_heated *MCGIDI_target_getHeatedTargetAtIndex_ReadIfNeeded( statusMessageReporting *smr, MCGIDI_target *target, int index );
494 MCGIDI_target_heated *MCGIDI_target_getHeatedTargetAtTIndex( statusMessageReporting *smr, MCGIDI_target *target, int index );
495 
496 int MCGIDI_target_numberOfReactions( statusMessageReporting *smr, MCGIDI_target *target );
497 enum MCGIDI_reactionType MCGIDI_target_getReactionTypeAtIndex( statusMessageReporting *smr, MCGIDI_target *target, int index );
498 MCGIDI_reaction *MCGIDI_target_getReactionAtIndex( MCGIDI_target *target, int index );
499 MCGIDI_reaction *MCGIDI_target_getReactionAtIndex_smr( statusMessageReporting *smr, MCGIDI_target *target, int index );
500 int MCGIDI_target_numberOfProductionReactions( statusMessageReporting *smr, MCGIDI_target *target );
501 
502 transportabilitiesMap const *MCGIDI_target_getUniqueProducts( statusMessageReporting *smr, MCGIDI_target *target );
503 int MCGIDI_target_recast( statusMessageReporting *smr, MCGIDI_target *target, GIDI_settings &settings );
504 
505 int MCGIDI_target_getDomain( statusMessageReporting *smr, MCGIDI_target *target, double *EMin, double *EMax );
506 double MCGIDI_target_getTotalCrossSectionAtTAndE( statusMessageReporting *smr, MCGIDI_target *target, MCGIDI_quantitiesLookupModes &modes,
507         bool sampling );
508 double MCGIDI_target_getIndexReactionCrossSectionAtE( statusMessageReporting *smr, MCGIDI_target *target, int index, MCGIDI_quantitiesLookupModes &modes,
509         bool sampling );
510 int MCGIDI_target_sampleReaction( statusMessageReporting *smr, MCGIDI_target *target, MCGIDI_quantitiesLookupModes &modes, double totalXSec, 
511         double (*userrng)( void * ), void *rngState );
512 int MCGIDI_target_sampleNullReactionProductsAtE( statusMessageReporting *smr, MCGIDI_target *target,
513     MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productDatas );
514 int MCGIDI_target_sampleIndexReactionProductsAtE( statusMessageReporting *smr, MCGIDI_target *target, int index, 
515         MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productData );
516 double MCGIDI_target_getIndexReactionFinalQ( statusMessageReporting *smr, MCGIDI_target *target, int index, MCGIDI_quantitiesLookupModes &modes );
517 
518 /*
519 * Routines in MCGIDI_target_heated.c
520 */
521 MCGIDI_target_heated *MCGIDI_target_heated_new( statusMessageReporting *smr );
522 int MCGIDI_target_heated_initialize( statusMessageReporting *smr, MCGIDI_target_heated *target );
523 MCGIDI_target_heated *MCGIDI_target_heated_newRead( statusMessageReporting *smr, const char *fileName );
524 MCGIDI_target_heated *MCGIDI_target_heated_free( statusMessageReporting *smr, MCGIDI_target_heated *target );
525 int MCGIDI_target_heated_release( statusMessageReporting *smr, MCGIDI_target_heated *target );
526 int MCGIDI_target_heated_read( statusMessageReporting *smr, MCGIDI_target_heated *target, const char *fileName );
527 int MCGIDI_target_heated_numberOfReactions( statusMessageReporting *smr, MCGIDI_target_heated *target );
528 int MCGIDI_target_heated_numberOfProductionReactions( statusMessageReporting *smr, MCGIDI_target_heated *target );
529 MCGIDI_reaction *MCGIDI_target_heated_getReactionAtIndex( MCGIDI_target_heated *target, int index );
530 MCGIDI_reaction *MCGIDI_target_heated_getReactionAtIndex_smr( statusMessageReporting *smr, MCGIDI_target_heated *target, int index );
531 #if 0
532 MCGIDI_reaction *MCGIDI_target_heated_getProductionReactionAtIndex( MCGIDI_target_heated *target, int index );
533 #endif
534 MCGIDI_POP *MCGIDI_target_heated_getPOPForProjectile( statusMessageReporting *smr, MCGIDI_target_heated *target );
535 MCGIDI_POP *MCGIDI_target_heated_getPOPForTarget( statusMessageReporting *smr, MCGIDI_target_heated *target );
536 double MCGIDI_target_heated_getProjectileMass_MeV( statusMessageReporting *smr, MCGIDI_target_heated *target );
537 double MCGIDI_target_heated_getTargetMass_MeV( statusMessageReporting *smr, MCGIDI_target_heated *target );
538 int MCGIDI_target_heated_getEnergyGrid( statusMessageReporting *smr, MCGIDI_target_heated *target, double **energyGrid );
539 double MCGIDI_target_heated_getTotalCrossSectionAtE( statusMessageReporting *smr, MCGIDI_target_heated *target, MCGIDI_quantitiesLookupModes &modes,
540         bool sampling );
541 double MCGIDI_target_heated_getIndexReactionCrossSectionAtE( statusMessageReporting *smr, MCGIDI_target_heated *target, int index, 
542         MCGIDI_quantitiesLookupModes &modes, bool sampling );
543 int MCGIDI_target_heated_sampleIndexReactionProductsAtE( statusMessageReporting *smr, MCGIDI_target_heated *target, int index, 
544         MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productData );
545 double MCGIDI_target_heated_getReactionsThreshold( statusMessageReporting *smr, MCGIDI_target_heated *target, int index );
546 int MCGIDI_target_heated_getReactionsDomain( statusMessageReporting *smr, MCGIDI_target_heated *target, int index, double *EMin, double *EMax );
547 double MCGIDI_target_heated_getIndexReactionFinalQ( statusMessageReporting *smr, MCGIDI_target_heated *target, int index, 
548         MCGIDI_quantitiesLookupModes &modes );
549 
550 transportabilitiesMap const *MCGIDI_target_heated_getUniqueProducts( statusMessageReporting *smr, MCGIDI_target_heated *target );
551 int MCGIDI_target_heated_recast( statusMessageReporting *smr, MCGIDI_target_heated *target, GIDI_settings &settings );
552 
553 /*
554 * Routines in MCGIDI_reaction.c
555 */
556 MCGIDI_reaction *MCGIDI_reaction_new( statusMessageReporting *smr );
557 int MCGIDI_reaction_initialize( statusMessageReporting *smr, MCGIDI_reaction *reaction );
558 MCGIDI_reaction *MCGIDI_reaction_free( statusMessageReporting *smr, MCGIDI_reaction *reaction );
559 int MCGIDI_reaction_release( statusMessageReporting *smr, MCGIDI_reaction *reaction );
560 int MCGIDI_reaction_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_target_heated *target, 
561     MCGIDI_POPs *pops, MCGIDI_reaction *reaction );
562 enum MCGIDI_reactionType MCGIDI_reaction_getReactionType( statusMessageReporting *smr, MCGIDI_reaction *reaction );
563 MCGIDI_target_heated *MCGIDI_reaction_getTargetHeated( statusMessageReporting *smr, MCGIDI_reaction *reaction );
564 double MCGIDI_reaction_getProjectileMass_MeV( statusMessageReporting *smr, MCGIDI_reaction *reaction );
565 double MCGIDI_reaction_getTargetMass_MeV( statusMessageReporting *smr, MCGIDI_reaction *reaction );
566 int MCGIDI_reaction_getDomain( statusMessageReporting *smr, MCGIDI_reaction *reaction, double *EMin, double *EMax );
567 int MCGIDI_reaction_fixDomains( statusMessageReporting *smr, MCGIDI_reaction *reaction, double EMin, double EMax, nfu_status *status );
568 double MCGIDI_reaction_getCrossSectionAtE( statusMessageReporting *smr, MCGIDI_reaction *reaction, MCGIDI_quantitiesLookupModes &modes, bool sampling );
569 double MCGIDI_reaction_getFinalQ( statusMessageReporting *smr, MCGIDI_reaction *reaction, MCGIDI_quantitiesLookupModes &modes );
570 int MCGIDI_reaction_getENDF_MTNumber( MCGIDI_reaction *reaction );
571 int MCGIDI_reaction_getENDL_CSNumbers( MCGIDI_reaction *reaction, int *S );
572 int MCGIDI_reaction_recast( statusMessageReporting *smr, MCGIDI_reaction *reaction, GIDI_settings &settings, 
573     GIDI_settings_particle const *projectileSettings, double temperature_MeV, ptwXPoints *totalGroupedCrossSection );
574 
575 MCGIDI_productsInfo *MCGIDI_reaction_getProductsInfo( MCGIDI_reaction *reaction );
576 int MCGIDI_productsInfo_getNumberOfUniqueProducts( MCGIDI_productsInfo *productsInfo );
577 int MCGIDI_productsInfo_getPoPsIndexAtIndex( MCGIDI_productsInfo *productsInfo, int index );
578 enum MCGIDI_productMultiplicityType MCGIDI_productsInfo_getMultiplicityTypeAtIndex( MCGIDI_productsInfo *productsInfo, int index );
579 int MCGIDI_productsInfo_getIntegerMultiplicityAtIndex( MCGIDI_productsInfo *productsInfo, int index );
580 int MCGIDI_productsInfo_getTransportableAtIndex( MCGIDI_productsInfo *productsInfo, int index );
581 
582 /*
583 * Routines in MCGIDI_pop.c
584 */
585 MCGIDI_POPs *MCGIDI_POPs_new( statusMessageReporting *smr, int size );
586 int MCGIDI_POPs_initial( statusMessageReporting *smr, MCGIDI_POPs *pops, int size );
587 void *MCGIDI_POPs_free( MCGIDI_POPs *pops );
588 int MCGIDI_POPs_release( MCGIDI_POPs *pops );
589 MCGIDI_POP *MCGIDI_POPs_addParticleIfNeeded( statusMessageReporting *smr, MCGIDI_POPs *pops, char const *name, double mass_MeV, 
590     double level_MeV, MCGIDI_POP *parent, int globalParticle );
591 int MCGIDI_POPs_findParticleIndex( MCGIDI_POPs *pops, char const *name );
592 MCGIDI_POP *MCGIDI_POPs_findParticle( MCGIDI_POPs *pops, char const *name );
593 void MCGIDI_POPs_writeSortedList( MCGIDI_POPs *pops, FILE *f );
594 void MCGIDI_POPs_printSortedList( MCGIDI_POPs *pops );
595 
596 MCGIDI_POP *MCGIDI_POP_new( statusMessageReporting *smr, char const *name, double mass_MeV, double level_MeV, MCGIDI_POP *parent );
597 MCGIDI_POP *MCGIDI_POP_free( MCGIDI_POP *pop );
598 MCGIDI_POP *MCGIDI_POP_release( MCGIDI_POP *pop );
599 double MCGIDI_POP_getMass_MeV( MCGIDI_POP *pop );
600 
601 /*
602 * Routines in MCGIDI_particle.c
603 */
604 MCGIDI_particle *MCGIDI_particle_new( statusMessageReporting *smr );
605 int MCGIDI_particle_initialize( statusMessageReporting *smr, MCGIDI_particle *particle );
606 MCGIDI_particle *MCGIDI_particle_free( statusMessageReporting *smr, MCGIDI_particle *particle );
607 int MCGIDI_particle_release( statusMessageReporting *smr, MCGIDI_particle *particle );
608 int MCGIDI_particle_freeInternalList( statusMessageReporting *smr );
609 MCGIDI_particle *MCGIDI_particle_getInternalID( statusMessageReporting *smr, const char * const name, MCGIDI_POPs *pops );
610 int MCGIDI_particle_printInternalSortedList( statusMessageReporting *smr );
611 
612 /*
613 * Routines in MCGIDI_outputChannel.c
614 */
615 MCGIDI_outputChannel *MCGIDI_outputChannel_new( statusMessageReporting *smr );
616 int MCGIDI_outputChannel_initialize( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel );
617 MCGIDI_outputChannel *MCGIDI_outputChannel_free( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel );
618 int MCGIDI_outputChannel_release( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel );
619 int MCGIDI_outputChannel_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_POPs *pops, MCGIDI_outputChannel *outputChannel,
620     MCGIDI_reaction *reaction, MCGIDI_product *parent );
621 int MCGIDI_outputChannel_numberOfProducts( MCGIDI_outputChannel *outputChannel );
622 MCGIDI_product *MCGIDI_outputChannel_getProductAtIndex( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel, int i );
623 int MCGIDI_outputChannel_getDomain( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel, double *EMin, double *EMax );
624 MCGIDI_target_heated *MCGIDI_outputChannel_getTargetHeated( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel );
625 double MCGIDI_outputChannel_getProjectileMass_MeV( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel );
626 double MCGIDI_outputChannel_getTargetMass_MeV( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel );
627 double MCGIDI_outputChannel_getQ_MeV( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel, double e_in );
628 double MCGIDI_outputChannel_getFinalQ( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel, double e_in );
629 int MCGIDI_outputChannel_sampleProductsAtE( statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel, MCGIDI_quantitiesLookupModes &modes,
630     MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productDatas, double *masses );
631 
632 /*
633 * Routines in MCGIDI_product.c
634 */
635 MCGIDI_product *MCGIDI_product_new( statusMessageReporting *smr );
636 int MCGIDI_product_initialize( statusMessageReporting *smr, MCGIDI_product *product );
637 MCGIDI_product *MCGIDI_product_free( statusMessageReporting *smr, MCGIDI_product *product );
638 int MCGIDI_product_release( statusMessageReporting *smr, MCGIDI_product *product );
639 int MCGIDI_product_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_outputChannel *outputChannel,
640         MCGIDI_POPs *pops, MCGIDI_product *product, int *delayedNeutronIndex );
641 int MCGIDI_product_getDomain( statusMessageReporting *smr, MCGIDI_product *product, double *EMin, double *EMax );
642 int MCGIDI_product_setTwoBodyMasses( statusMessageReporting *smr, MCGIDI_product *product, double projectileMass_MeV, double targetMass_MeV,
643     double productMass_MeV, double residualMass_MeV );
644 double MCGIDI_product_getMass_MeV( statusMessageReporting *smr, MCGIDI_product *product );
645 MCGIDI_target_heated *MCGIDI_product_getTargetHeated( statusMessageReporting *smr, MCGIDI_product *product );
646 double MCGIDI_product_getProjectileMass_MeV( statusMessageReporting *smr, MCGIDI_product *product );
647 double MCGIDI_product_getTargetMass_MeV( statusMessageReporting *smr, MCGIDI_product *product );
648 int MCGIDI_product_sampleMultiplicity( statusMessageReporting *smr, MCGIDI_product *product, double e_in, double r );
649 int MCGIDI_product_sampleMu( statusMessageReporting *smr, MCGIDI_product *product, MCGIDI_quantitiesLookupModes &modes,
650     MCGIDI_decaySamplingInfo *decaySamplingInfo );
651 
652 int MCGIDI_sampledProducts_initialize( statusMessageReporting *smr, MCGIDI_sampledProductsDatas *sampledProductsDatas, int incrementSize );
653 int MCGIDI_sampledProducts_release( statusMessageReporting *smr, MCGIDI_sampledProductsDatas *sampledProductsDatas );
654 int MCGIDI_sampledProducts_remalloc( statusMessageReporting *smr, MCGIDI_sampledProductsDatas *sampledProductsDatas );
655 int MCGIDI_sampledProducts_addProduct( statusMessageReporting *smr, MCGIDI_sampledProductsDatas *sampledProductsDatas, 
656     MCGIDI_sampledProductsData *sampledProductsData );
657 int MCGIDI_sampledProducts_number( MCGIDI_sampledProductsDatas *sampledProductsDatas );
658 MCGIDI_sampledProductsData *MCGIDI_sampledProducts_getProductAtIndex( MCGIDI_sampledProductsDatas *sampledProductsDatas, int index );
659 
660 /*
661 * Routines in MCGIDI_distribution.c
662 */
663 MCGIDI_distribution *MCGIDI_distribution_new( statusMessageReporting *smr );
664 int MCGIDI_distribution_initialize( statusMessageReporting *smr, MCGIDI_distribution *distribution );
665 MCGIDI_distribution *MCGIDI_distribution_free( statusMessageReporting *smr, MCGIDI_distribution *distribution );
666 int MCGIDI_distribution_release( statusMessageReporting *smr, MCGIDI_distribution *distribution );
667 int MCGIDI_distribution_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_product *product, MCGIDI_POPs *pops, ptwXYPoints *norms );
668 
669 /*
670 * Routines in MCGIDI_angular.c
671 */
672 MCGIDI_angular *MCGIDI_angular_new( statusMessageReporting *smr );
673 int MCGIDI_angular_initialize( statusMessageReporting *smr, MCGIDI_angular *angular );
674 MCGIDI_angular *MCGIDI_angular_free( statusMessageReporting *smr, MCGIDI_angular *angular );
675 int MCGIDI_angular_release( statusMessageReporting *smr, MCGIDI_angular *angular );
676 int MCGIDI_angular_setTwoBodyMasses( statusMessageReporting *smr, MCGIDI_angular *angular, double projectileMass_MeV, double targetMass_MeV,
677     double productMass_MeV, double residualMass_MeV );
678 int MCGIDI_angular_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution, ptwXYPoints *norms );
679 int MCGIDI_angular_sampleMu( statusMessageReporting *smr, MCGIDI_angular *angular, MCGIDI_quantitiesLookupModes &modes,
680     MCGIDI_decaySamplingInfo *decaySamplingInfo );
681 
682 /*
683 * Routines in MCGIDI_energy.c
684 */
685 MCGIDI_energy *MCGIDI_energy_new( statusMessageReporting *smr );
686 int MCGIDI_energy_initialize( statusMessageReporting *smr, MCGIDI_energy *energy );
687 MCGIDI_energy *MCGIDI_energy_free( statusMessageReporting *smr, MCGIDI_energy *energy );
688 int MCGIDI_energy_release( statusMessageReporting *smr, MCGIDI_energy *energy );
689 int MCGIDI_energy_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution, ptwXYPoints *norms,
690     enum MCGIDI_energyType energyType, double gammaEnergy_MeV );
691 int MCGIDI_energy_sampleEnergy( statusMessageReporting *smr, MCGIDI_energy *energy, MCGIDI_quantitiesLookupModes &modes, 
692     MCGIDI_decaySamplingInfo *decaySamplingInfo );
693 
694 /*
695 * Routines in MCGIDI_energyAngular.c
696 */
697 int MCGIDI_energyAngular_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution );
698 MCGIDI_energyAngular *MCGIDI_energyAngular_new( statusMessageReporting *smr );
699 int MCGIDI_energyAngular_initialize( statusMessageReporting *smr, MCGIDI_energyAngular *energyAngular );
700 MCGIDI_energyAngular *MCGIDI_energyAngular_free( statusMessageReporting *smr, MCGIDI_energyAngular *energyAngular );
701 int MCGIDI_energyAngular_release( statusMessageReporting *smr, MCGIDI_energyAngular *energyAngular );
702 int MCGIDI_energyAngular_sampleDistribution( statusMessageReporting *smr, MCGIDI_distribution *distribution, MCGIDI_quantitiesLookupModes &modes,
703     MCGIDI_decaySamplingInfo *decaySamplingInfo );
704 
705 /*
706 * Routines in MCGIDI_angularEnergy.c
707 */
708 MCGIDI_angularEnergy *MCGIDI_angularEnergy_new( statusMessageReporting *smr );
709 int MCGIDI_angularEnergy_initialize( statusMessageReporting *smr, MCGIDI_angularEnergy *energyAngular );
710 MCGIDI_angularEnergy *MCGIDI_angularEnergy_free( statusMessageReporting *smr, MCGIDI_angularEnergy *energyAngular );
711 int MCGIDI_angularEnergy_release( statusMessageReporting *smr, MCGIDI_angularEnergy *energyAngular );
712 int MCGIDI_angularEnergy_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution );
713 int MCGIDI_angularEnergy_sampleDistribution( statusMessageReporting *smr, MCGIDI_angularEnergy *angularEnergy, MCGIDI_quantitiesLookupModes &modes,
714     MCGIDI_decaySamplingInfo *decaySamplingInfo );
715 
716 /*
717 * Routines in MCGIDI_KalbachMann.c
718 */
719 MCGIDI_KalbachMann *MCGIDI_KalbachMann_new( statusMessageReporting *smr, ptwXY_interpolation interpolationWY, ptwXY_interpolation interpolationXY );
720 int MCGIDI_KalbachMann_initialize( statusMessageReporting *smr, MCGIDI_KalbachMann *KalbachMann, ptwXY_interpolation interpolationWY, ptwXY_interpolation interpolationXY );
721 MCGIDI_KalbachMann *MCGIDI_KalbachMann_free( statusMessageReporting *smr, MCGIDI_KalbachMann *KalbachMann );
722 int MCGIDI_KalbachMann_release( statusMessageReporting *smr, MCGIDI_KalbachMann *KalbachMann );
723 int MCGIDI_KalbachMann_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution );
724 int MCGIDI_KalbachMann_sampleEp( statusMessageReporting *smr, MCGIDI_KalbachMann *KalbachMann, MCGIDI_quantitiesLookupModes &modes, 
725     MCGIDI_decaySamplingInfo *decaySamplingInfo );
726 
727 /*
728 * Routines in MCGIDI_uncorrelated.c
729 */
730 int MCGIDI_uncorrelated_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution, ptwXYPoints *norms,
731     enum MCGIDI_energyType energyType, double gammaEnergy_MeV );
732 int MCGIDI_uncorrelated_sampleDistribution( statusMessageReporting *smr, MCGIDI_distribution *distribution, MCGIDI_quantitiesLookupModes &modes,
733     MCGIDI_decaySamplingInfo *decaySamplingInfo );
734 
735 /*
736 * Routines in MCGIDI_LLNLAngular_angularEnergy.c
737 */
738 int MCGIDI_LLNLAngular_angularEnergy_parseFromTOM( statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_distribution *distribution );
739 
740 /*
741 * Routines in MCGIDI_kinetics.c
742 */
743 int MCGIDI_kinetics_2BodyReaction( statusMessageReporting *smr, MCGIDI_angular *angular, double K, double mu, double phi,
744         MCGIDI_sampledProductsData *outgoingData );
745 int MCGIDI_kinetics_COMKineticEnergy2LabEnergyAndMomentum( statusMessageReporting *smr, double beta, double e_kinetic_com, double mu, double phi,
746         double m3cc, double m4cc, MCGIDI_sampledProductsData *outgoingData );
747 int MCGIDI_kinetics_COM2Lab( statusMessageReporting *smr, MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, double masses[3] );
748 
749 /*
750 * Routines in MCGIDI_sampling.c
751 */
752 int MCGIDI_sampling_pdfsOfXGivenW_initialize( statusMessageReporting *smr, MCGIDI_pdfsOfXGivenW *dists );
753 int MCGIDI_sampling_pdfsOfXGivenW_release( statusMessageReporting *smr, MCGIDI_pdfsOfXGivenW *dists );
754 int MCGIDI_sampling_pdfsOfX_release( statusMessageReporting *smr, MCGIDI_pdfOfX *dist );
755 int MCGIDI_sampling_sampleX_from_pdfsOfXGivenW( MCGIDI_pdfsOfXGivenW *dists, MCGIDI_pdfsOfXGivenW_sampled *sampled, double r );
756 int MCGIDI_sampling_sampleX_from_pdfOfX( MCGIDI_pdfOfX *dist, MCGIDI_pdfsOfXGivenW_sampled *sampled, double r );
757 int MCGIDI_sampling_doubleDistribution( statusMessageReporting *smr, MCGIDI_pdfsOfXGivenW *pdfOfWGivenV, MCGIDI_pdfsOfXGivenW *pdfOfXGivenVAndW,  
758         MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo );
759 int MCGIDI_sampling_interpolationValues( statusMessageReporting *smr, ptwXY_interpolation interpolation, double *ws, double y1, double y2, double *y );
760 double MCGIDI_sampling_ptwXY_getValueAtX( ptwXYPoints *ptwXY, double x1 );
761 
762 /*
763 * Routines in MCGIDI_misc.c
764 */
765 int MCGIDI_misc_NumberOfZSymbols( void );
766 const char *MCGIDI_misc_ZToSymbol( int iZ );
767 int MCGIDI_misc_symbolToZ( const char *Z );
768 int MCGIDI_miscNameToZAm( statusMessageReporting *smr, const char *name, int *Z, int *A, int *m, int *level );
769 xDataTOM_Int MCGIDI_misc_binarySearch( xDataTOM_Int n, double *ds, double d );
770 int MCGIDI_misc_PQUStringToDouble( statusMessageReporting *smr, char const *str, char const *unit, double conversion, double *value );
771 int MCGIDI_misc_PQUStringToDoubleInUnitOf( statusMessageReporting *smr, char const *str, char const *toUnit, double *value );
772 void MCGIDI_misc_updateTransportabilitiesMap( transportabilitiesMap *transportabilities, int PoPID, enum MCGIDI_transportability transportability );
773 void MCGIDI_misc_updateTransportabilitiesMap2( transportabilitiesMap *transportabilities, int PoPID, int transportable );
774 
775 #if defined __cplusplus
776     }
777     }
778 #endif
779 
780 #endif          /* End of MCGIDI_h_included. */
781