15 #if defined __cplusplus
103 double crossSectionInit[4] = { 0., 0., 0., 0., };
106 int subtag1_Notice = 0;
116 if( strcmp( version,
"GND 1.3" ) != 0 ) {
121 if( strcmp( element->
name,
"reactionSuite" ) != 0 ) {
132 while( isspace( *e1 ) ) ++e1;
152 if( strcmp( child->
name,
"particles" ) == 0 )
continue;
153 if( strcmp( child->
name,
"styles" ) == 0 )
continue;
154 if( strcmp( child->
name,
"documentations" ) == 0 )
continue;
155 if( strcmp( child->
name,
"resonances" ) == 0 )
continue;
156 if( strcmp( child->
name,
"summedReaction" ) == 0 )
continue;
157 if( strcmp( child->
name,
"fissionComponent" ) == 0 )
continue;
158 if( strcmp( child->
name,
"reaction" ) == 0 ) {
164 if( ir == 0 ) { target->
EMin = EMin; target->
EMax = EMax; }
165 if( EMin < target->EMin ) target->
EMin = EMin;
166 if( EMax > target->
EMax ) target->
EMax = EMax;
167 for( transportabilitiesMap::const_iterator iter = reaction->
transportabilities->begin( );
172 else if( strcmp( child->
name,
"production" ) == 0 ) {
174 else if( strcmp( child->
name,
"aliases" ) == 0 ) {
176 else if( strcmp( child->
name,
"partialGammaProduction" ) == 0 ) {
177 if( subtag1_Notice == 0 )
printf(
"Unsupported reactionSuite sub-tag = '%s'\n", child->
name );
180 printf(
"Unsupported reactionSuite sub-tag = '%s'\n", child->
name );
183 crossSectionInit[0] = target->
EMin;
184 crossSectionInit[2] = target->
EMax;
221 if( strcmp( child->
name,
"particle" ) ) {
263 int globalParticle = 1;
264 char const *
name = NULL, *mass = NULL;
281 if( strcmp( child->
name,
"level" ) ) {
299 int globalParticle = 0;
300 char const *
name, *level, *aliasValue;
301 double level_MeV = 0.;
319 if( strcmp( aliasValue, name ) == 0 ) globalParticle = 1;
337 char const *finalLevelString;
341 if( strcmp( child->
name,
"gamma" ) ) {
347 if( gammaCounts > 0 ) {
403 if( ( index >= 0 ) && ( index < target->numberOfReactions ) )
return( &(target->
reactions[index]) );
413 if( reaction == NULL ) {
424 MCGIDI_channel *channel = NULL;
426 if( ( index >= 0 ) && ( index < target->nProductionReactions ) ) channel = target->productionReactions[index];
469 if( e_in < target->EMin ) e_in = target->
EMin;
470 if( e_in > target->
EMax ) e_in = target->
EMax;
482 if( xSecP != NULL ) {
514 if( reaction == NULL )
return( -1 );
524 if( reaction == NULL )
return( -1 );
525 return( reaction->
EMin );
534 if( reaction == NULL )
return( -1 );
535 *EMin = reaction->
EMin;
536 *EMax = reaction->
EMax;
547 if( reaction == NULL )
return( 0. );
567 if( projectileSettings == NULL ) {
576 if( ( totalGroupedCrossSection =
ptwX_createLine( numberOfGroups, numberOfGroups, 0, 0, &status_nf ) ) == NULL ) {
604 totalGroupedCrossSection = NULL;
614 #if defined __cplusplus
double MCGIDI_target_heated_getReactionsThreshold(statusMessageReporting *smr, MCGIDI_target_heated *target, int index)
MCGIDI_reaction * MCGIDI_target_heated_getReactionAtIndex_smr(statusMessageReporting *smr, MCGIDI_target_heated *target, int index)
xDataTOM_element * xDataTOME_getNextElement(xDataTOM_element *element)
int MCGIDI_target_heated_read(statusMessageReporting *smr, MCGIDI_target_heated *target, const char *fileName)
ptwXPoints * ptwX_free(ptwXPoints *ptwX)
static int MCGIDI_target_heated_parseParticleGammas(statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, char const *name)
MCGIDI_target_heated * MCGIDI_target_heated_newRead(statusMessageReporting *smr, const char *fileName)
char * xDataTOMMisc_getAbsPath(statusMessageReporting *smr, const char *fileName)
ptwXPoints * crossSectionGroupedForSampling
GIDI_settings_particle const * getParticle(int PoPId) const
static int MCGIDI_target_heated_parseReaction(statusMessageReporting *smr, xDataTOM_element *child, MCGIDI_target_heated *target, MCGIDI_POPs *pops, MCGIDI_reaction *reaction)
transportabilitiesMap * transportabilities
double MCGIDI_misc_getUnitConversionFactor(statusMessageReporting *smr, char const *fromUnit, char const *toUnit)
#define smr_allocateCopyString2(smr, s, forItem)
#define smr_setReportError2(smr, libraryID, code, fmt,...)
int getNumberOfGroups(void) const
MCGIDI_POP * MCGIDI_POPs_addParticleIfNeeded(statusMessageReporting *smr, MCGIDI_POPs *pops, char const *name, double mass_MeV, double level_MeV, MCGIDI_POP *parent, int globalParticle)
void MCGIDI_misc_updateTransportabilitiesMap(transportabilitiesMap *transportabilities, int PoPID, enum MCGIDI_transportability transportability)
MCGIDI_POP * MCGIDI_target_heated_getPOPForProjectile(statusMessageReporting *smr, MCGIDI_target_heated *target)
MCGIDI_reaction * reactions
double MCGIDI_target_heated_getTargetMass_MeV(statusMessageReporting *smr, MCGIDI_target_heated *target)
xDataTOM_TOM * xDataTOM_importFile(statusMessageReporting *smr, const char *fileName)
int xDataTOM_numberOfElementsByName(statusMessageReporting *smr, xDataTOM_element *element, char const *name)
ptwXYPoints * crossSection
double MCGIDI_target_heated_getIndexReactionCrossSectionAtE(statusMessageReporting *smr, MCGIDI_target_heated *target, int index, MCGIDI_quantitiesLookupModes &modes, bool sampling)
static int MCGIDI_target_heated_parsePOPs(statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, xDataTOM_element *particleAliases)
int64_t ptwX_length(ptwXPoints *ptwX)
int MCGIDI_target_heated_numberOfReactions(statusMessageReporting *smr, MCGIDI_target_heated *target)
int MCGIDI_POPs_initial(statusMessageReporting *smr, MCGIDI_POPs *pops, int size)
int MCGIDI_reaction_getDomain(statusMessageReporting *smr, MCGIDI_reaction *reaction, double *EMin, double *EMax)
xDataTOM_element * xDataTOME_getFirstElement(xDataTOM_element *element)
double ptwX_getPointAtIndex_Unsafely(ptwXPoints *ptwX, int64_t index)
int MCGIDI_misc_PQUStringToDoubleInUnitOf(statusMessageReporting *smr, char const *str, char const *toUnit, double *value)
MCGIDI_target_heated * MCGIDI_target_heated_free(statusMessageReporting *smr, MCGIDI_target_heated *target)
double * ptwX_getPointAtIndex(ptwXPoints *ptwX, int64_t index)
int MCGIDI_target_heated_numberOfProductionReactions(statusMessageReporting *smr, MCGIDI_target_heated *target)
int MCGIDI_target_heated_sampleIndexReactionProductsAtE(statusMessageReporting *smr, MCGIDI_target_heated *target, int index, MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productData)
void * xDataTOM_freeTOM(statusMessageReporting *smr, xDataTOM_TOM **TOM)
ptwXYPoints * crossSection
transportabilitiesMap const * MCGIDI_target_heated_getUniqueProducts(statusMessageReporting *smr, MCGIDI_target_heated *target)
int MCGIDI_POPs_release(MCGIDI_POPs *pops)
int getGroupIndex(void) const
nfu_status ptwX_setPointAtIndex(ptwXPoints *ptwX, int64_t index, double x)
int MCGIDI_target_heated_getReactionsDomain(statusMessageReporting *smr, MCGIDI_target_heated *target, int index, double *EMin, double *EMax)
#define smr_malloc2(smr, size, zero, forItem)
xDataTOM_element * xDataTOME_getOneElementByName(statusMessageReporting *smr, xDataTOM_element *element, char const *name, int required)
MCGIDI_POP * projectilePOP
int MCGIDI_target_heated_recast(statusMessageReporting *smr, MCGIDI_target_heated *target, GIDI_settings &settings)
char const * xDataTOM_getAttributesValueInElement(xDataTOM_element *element, char const *name)
int MCGIDI_target_heated_release(statusMessageReporting *smr, MCGIDI_target_heated *target)
int MCGIDI_reaction_recast(statusMessageReporting *smr, MCGIDI_reaction *reaction, GIDI_settings &settings, GIDI_settings_particle const *projectileSettings, double temperature_MeV, ptwXPoints *totalGroupedCrossSection)
double MCGIDI_target_heated_getIndexReactionFinalQ(statusMessageReporting *smr, MCGIDI_target_heated *target, int index, MCGIDI_quantitiesLookupModes &modes)
MCGIDI_POP * MCGIDI_target_heated_getPOPForTarget(statusMessageReporting *smr, MCGIDI_target_heated *target)
void xDataTOMAL_release(xDataTOM_attributionList *attributes)
double thresholdGroupedDeltaCrossSection
MCGIDI_outputChannel outputChannel
double MCGIDI_reaction_getCrossSectionAtE(statusMessageReporting *smr, MCGIDI_reaction *reaction, MCGIDI_quantitiesLookupModes &modes, bool sampling)
int smr_isOk(statusMessageReporting *smr)
bool isEnergyMode_grouped(void) const
enum nfu_status_e nfu_status
printf("%d Experimental points found\n", nlines)
nfu_status ptwXY_getValueAtX(ptwXYPoints *ptwXY, double x, double *y)
transportabilitiesMap * transportabilities
enum MCGIDI_quantityLookupMode getCrossSectionMode(void) const
#define smr_setReportError2p(smr, libraryID, code, fmt)
xDataTOM_attributionList attributes
int MCGIDI_target_heated_initialize(statusMessageReporting *smr, MCGIDI_target_heated *target)
static int MCGIDI_target_heated_parseParticle(statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, xDataTOM_element *particleAliases)
int MCGIDI_misc_PQUStringToDouble(statusMessageReporting *smr, char const *str, char const *unit, double conversion, double *value)
void * smr_freeMemory(void **p)
double MCGIDI_target_heated_getProjectileMass_MeV(statusMessageReporting *smr, MCGIDI_target_heated *target)
MCGIDI_POP * MCGIDI_POPs_findParticle(MCGIDI_POPs *pops, char const *name)
int xDataTOMAL_copyAttributionList(statusMessageReporting *smr, xDataTOM_attributionList *desc, xDataTOM_attributionList *src)
xDataTOM_element * xDataTOM_getDocumentsElement(xDataTOM_TOM *TOM)
ptwXPoints * crossSectionGrouped
ptwXYPoints * ptwXY_free(ptwXYPoints *ptwXY)
xDataTOM_attributionList attributes
ptwXYPoints * ptwXY_create(ptwXY_interpolation interpolation, ptwXY_interpolationOtherInfo const *interpolationOtherInfo, double biSectionMax, double accuracy, int64_t primarySize, int64_t secondarySize, int64_t length, double const *xy, nfu_status *status, int userFlag)
double getProjectileEnergy(void) const
double MCGIDI_POP_getMass_MeV(MCGIDI_POP *pop)
int MCGIDI_outputChannel_sampleProductsAtE(statusMessageReporting *smr, MCGIDI_outputChannel *outputChannel, MCGIDI_quantitiesLookupModes &modes, MCGIDI_decaySamplingInfo *decaySamplingInfo, MCGIDI_sampledProductsDatas *productDatas, double *masses)
double MCGIDI_target_heated_getTotalCrossSectionAtE(statusMessageReporting *smr, MCGIDI_target_heated *target, MCGIDI_quantitiesLookupModes &modes, bool sampling)
#define MCGIDI_misc_pointerToTOMAttributeIfAllOk3(smr, path, required, attributes, name)
ptwXPoints * ptwX_clone(ptwXPoints *ptwX, nfu_status *status)
static int MCGIDI_target_heated_parseParticleLevel(statusMessageReporting *smr, MCGIDI_target_heated *target, xDataTOM_element *element, MCGIDI_POP *parent, double mass_MeV, xDataTOM_element *particleAliases)
ptwXYPoints * ptwXY_add_ptwXY(ptwXYPoints *ptwXY1, ptwXYPoints *ptwXY2, nfu_status *status)
int xDataTOME_convertAttributeToDouble(statusMessageReporting *smr, xDataTOM_element *element, char const *name, double *d)
std::map< int, enum GIDI::MCGIDI_transportability > transportabilitiesMap
const char * nfu_statusMessage(nfu_status status)
MCGIDI_target_heated * MCGIDI_target_heated_new(statusMessageReporting *smr)
int MCGIDI_reaction_release(statusMessageReporting *smr, MCGIDI_reaction *reaction)
int MCGIDI_reaction_fixDomains(statusMessageReporting *smr, MCGIDI_reaction *reaction, double EMin, double EMax, nfu_status *status)
ptwXPoints * ptwX_createLine(int64_t size, int64_t length, double slope, double offset, nfu_status *status)
MCGIDI_reaction * MCGIDI_target_heated_getReactionAtIndex(MCGIDI_target_heated *target, int index)
int MCGIDI_reaction_parseFromTOM(statusMessageReporting *smr, xDataTOM_element *element, MCGIDI_target_heated *target, MCGIDI_POPs *pops, MCGIDI_reaction *reaction)
double MCGIDI_reaction_getFinalQ(statusMessageReporting *smr, MCGIDI_reaction *reaction, MCGIDI_quantitiesLookupModes &modes)
MCGIDI_GammaBranching * gammas