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