int main()
{
struct compoundData *cdtest;
int i;
XRayInit();
SetErrorMessages(0);
//if something goes wrong, the test will end with EXIT_FAILURE
//SetHardExit(1);
std::printf("Example of C++ program using xraylib\n");
std::printf("Density of pure Al: %f g/cm3\n", ElementDensity(13));
std::printf("Ca K-alpha Fluorescence Line Energy: %f\n",
LineEnergy(20,KA_LINE));
std::printf("Fe partial photoionization cs of L3 at 6.0 keV: %f\n",CS_Photo_Partial(26,L3_SHELL,6.0));
std::printf("Zr L1 edge energy: %f\n",EdgeEnergy(40,L1_SHELL));
std::printf("Pb Lalpha XRF production cs at 20.0 keV (jump approx): %f\n",CS_FluorLine(82,LA_LINE,20.0));
std::printf("Pb Lalpha XRF production cs at 20.0 keV (Kissel): %f\n",CS_FluorLine_Kissel(82,LA_LINE,20.0));
std::printf("Bi M1N2 radiative rate: %f\n",RadRate(83,M1N2_LINE));
std::printf("U M3O3 Fluorescence Line Energy: %f\n",LineEnergy(92,M3O3_LINE));
//parser test for Ca(HCO3)2 (calcium bicarbonate)
if ((cdtest = CompoundParser("Ca(HCO3)2")) == NULL)
return 1;
std::printf("Ca(HCO3)2 contains %g atoms, %i elements and has a molar mass of %g g/mol\n", cdtest->nAtomsAll, cdtest->nElements, cdtest->molarMass);
for (i = 0 ; i < cdtest->nElements ; i++)
std::printf("Element %i: %f %% and %g atoms\n", cdtest->Elements[i], cdtest->massFractions[i]*100.0, cdtest->nAtoms[i]);
FreeCompoundData(cdtest);
//parser test for SiO2 (quartz)
if ((cdtest = CompoundParser("SiO2")) == NULL)
return 1;
std::printf("SiO2 contains %g atoms, %i elements and has a molar mass of %g g/mol\n", cdtest->nAtomsAll, cdtest->nElements, cdtest->molarMass);
for (i = 0 ; i < cdtest->nElements ; i++)
std::printf("Element %i: %f %% and %g atoms\n", cdtest->Elements[i], cdtest->massFractions[i]*100.0, cdtest->nAtoms[i]);
FreeCompoundData(cdtest);
std::printf("Ca(HCO3)2 Rayleigh cs at 10.0 keV: %f\n",CS_Rayl_CP("Ca(HCO3)2",10.0f) );
std::printf("CS2 Refractive Index at 10.0 keV : %f - %f i\n",Refractive_Index_Re("CS2",10.0f,1.261f),Refractive_Index_Im("CS2",10.0f,1.261f));
std::printf("C16H14O3 Refractive Index at 1 keV : %f - %f i\n",Refractive_Index_Re("C16H14O3",1.0f,1.2f),Refractive_Index_Im("C16H14O3",1.0f,1.2f));
std::printf("SiO2 Refractive Index at 5 keV : %f - %f i\n",Refractive_Index_Re("SiO2",5.0f,2.65f),Refractive_Index_Im("SiO2",5.0f,2.65f));
std::printf("Compton profile for Fe at pz = 1.1 : %f\n",ComptonProfile(26,1.1f));
std::printf("M5 Compton profile for Fe at pz = 1.1 : %f\n",ComptonProfile_Partial(26,M5_SHELL,1.1f));
std::printf("K atomic level width for Fe: %f\n", AtomicLevelWidth(26,K_SHELL));
std::printf("M1->M5 Coster-Kronig transition probability for Au : %f\n",CosKronTransProb(79,FM15_TRANS));
std::printf("L1->L3 Coster-Kronig transition probability for Fe : %f\n",CosKronTransProb(26,FL13_TRANS));
std::printf("Au Ma1 XRF production cs at 10.0 keV (Kissel): %f\n", CS_FluorLine_Kissel(79,MA1_LINE,10.0f));
std::printf("Au Mb XRF production cs at 10.0 keV (Kissel): %f\n", CS_FluorLine_Kissel(79,MB_LINE,10.0f));
std::printf("Au Mg XRF production cs at 10.0 keV (Kissel): %f\n", CS_FluorLine_Kissel(79,MG_LINE,10.0f));
std::printf("Bi L2-M5M5 Auger non-radiative rate: %f\n",AugerRate(86,L2_M5M5_AUGER));
std::printf("Bi L3 Auger yield: %f\n", AugerYield(86, L3_SHELL));
std::printf("Sr anomalous scattering factor Fi at 10.0 keV: %f\n", Fi(38, 10.0));
std::printf("Sr anomalous scattering factor Fii at 10.0 keV: %f\n", Fii(38, 10.0));
char *symbol = AtomicNumberToSymbol(26);
std::printf("Symbol of element 26 is: %s\n",symbol);
xrlFree(symbol);
std::printf("Number of element Fe is: %i\n",SymbolToAtomicNumber("Fe"));
std::printf("Pb Malpha XRF production cs at 20.0 keV with cascade effect: %f\n",CS_FluorLine_Kissel(82,MA1_LINE,20.0));
std::printf("Pb Malpha XRF production cs at 20.0 keV with radiative cascade effect: %f\n",CS_FluorLine_Kissel_Radiative_Cascade(82,MA1_LINE,20.0));
std::printf("Pb Malpha XRF production cs at 20.0 keV with non-radiative cascade effect: %f\n",CS_FluorLine_Kissel_Nonradiative_Cascade(82,MA1_LINE,20.0));
std::printf("Pb Malpha XRF production cs at 20.0 keV without cascade effect: %f\n",CS_FluorLine_Kissel_no_Cascade(82,MA1_LINE,20.0));
std::printf("Al mass energy-absorption cs at 20.0 keV: %f\n", CS_Energy(13, 20.0));
std::printf("Pb mass energy-absorption cs at 40.0 keV: %f\n", CS_Energy(82, 40.0));
std::printf("CdTe mass energy-absorption cs at 40.0 keV: %f\n", CS_Energy_CP("CdTe", 40.0));
/* Si Crystal structure */
Crystal_Struct* cryst = Crystal_GetCrystal("Si", NULL);
if (cryst == NULL) return 1;
std::printf ("Si unit cell dimensions are %f %f %f\n", cryst->a, cryst->b, cryst->c);
std::printf ("Si unit cell angles are %f %f %f\n", cryst->alpha, cryst->beta, cryst->gamma);
std::printf ("Si unit cell volume is %f\n", cryst->volume);
std::printf ("Si atoms at:\n");
std::printf (" Z fraction X Y Z\n");
Crystal_Atom* atom;
for (i = 0; i < cryst->n_atom; i++) {
atom = &cryst->atom[i];
std::printf (" %3i %f %f %f %f\n", atom->Zatom, atom->fraction, atom->x, atom->y, atom->z);
}
/* Si diffraction parameters */
std::printf ("\nSi111 at 8 KeV. Incidence at the Bragg angle:\n");
double energy = 8;
double debye_temp_factor = 1.0;
double rel_angle = 1.0;
double bragg = Bragg_angle (cryst, energy, 1, 1, 1);
std::printf (" Bragg angle: Rad: %f Deg: %f\n", bragg, bragg*180/PI);
double q = Q_scattering_amplitude (cryst, energy, 1, 1, 1, rel_angle);
std::printf (" Q Scattering amplitude: %f\n", q);
double f0, fp, fpp;
Atomic_Factors (14, energy, q, debye_temp_factor, &f0, &fp, &fpp);
std::printf (" Atomic factors (Z = 14) f0, fp, fpp: %f, %f, i*%f\n", f0, fp, fpp);
xrlComplex FH, F0;
FH = Crystal_F_H_StructureFactor (cryst, energy, 1, 1, 1, debye_temp_factor, rel_angle);
std::printf (" FH(1,1,1) structure factor: (%f, %f)\n", FH.re, FH.im);
F0 = Crystal_F_H_StructureFactor (cryst, energy, 0, 0, 0, debye_temp_factor, rel_angle);
std::printf (" F0=FH(0,0,0) structure factor: (%f, %f)\n", F0.re, F0.im);
/* Diamond diffraction parameters */
cryst = Crystal_GetCrystal("Diamond", NULL);
std::printf ("\nDiamond 111 at 8 KeV. Incidence at the Bragg angle:\n");
bragg = Bragg_angle (cryst, energy, 1, 1, 1);
std::printf (" Bragg angle: Rad: %f Deg: %f\n", bragg, bragg*180/PI);
q = Q_scattering_amplitude (cryst, energy, 1, 1, 1, rel_angle);
std::printf (" Q Scattering amplitude: %f\n", q);
Atomic_Factors (6, energy, q, debye_temp_factor, &f0, &fp, &fpp);
std::printf (" Atomic factors (Z = 6) f0, fp, fpp: %f, %f, i*%f\n", f0, fp, fpp);
FH = Crystal_F_H_StructureFactor (cryst, energy, 1, 1, 1, debye_temp_factor, rel_angle);
std::printf (" FH(1,1,1) structure factor: (%f, %f)\n", FH.re, FH.im);
F0 = Crystal_F_H_StructureFactor (cryst, energy, 0, 0, 0, debye_temp_factor, rel_angle);
std::printf (" F0=FH(0,0,0) structure factor: (%f, %f)\n", F0.re, F0.im);
xrlComplex FHbar = Crystal_F_H_StructureFactor (cryst, energy, -1, -1, -1, debye_temp_factor, rel_angle);
double dw = 1e10 * 2 * (R_E / cryst->volume) * (KEV2ANGST * KEV2ANGST/ (energy * energy)) *
std::sqrt(c_abs(c_mul(FH, FHbar))) / PI / std::sin(2*bragg);
std::printf (" Darwin width: %f micro-radians\n", 1e6*dw);
/* Alpha Quartz diffraction parameters */
cryst = Crystal_GetCrystal("AlphaQuartz", NULL);
std::printf ("\nAlpha Quartz 020 at 8 KeV. Incidence at the Bragg angle:\n");
bragg = Bragg_angle (cryst, energy, 0, 2, 0);
std::printf (" Bragg angle: Rad: %f Deg: %f\n", bragg, bragg*180/PI);
q = Q_scattering_amplitude (cryst, energy, 0, 2, 0, rel_angle);
std::printf (" Q Scattering amplitude: %f\n", q);
Atomic_Factors (8, energy, q, debye_temp_factor, &f0, &fp, &fpp);
std::printf (" Atomic factors (Z = 8) f0, fp, fpp: %f, %f, i*%f\n", f0, fp, fpp);
FH = Crystal_F_H_StructureFactor (cryst, energy, 0, 2, 0, debye_temp_factor, rel_angle);
std::printf (" FH(0,2,0) structure factor: (%f, %f)\n", FH.re, FH.im);
F0 = Crystal_F_H_StructureFactor (cryst, energy, 0, 0, 0, debye_temp_factor, rel_angle);
std::printf (" F0=FH(0,0,0) structure factor: (%f, %f)\n", F0.re, F0.im);
/* Muscovite diffraction parameters */
cryst = Crystal_GetCrystal("Muscovite", NULL);
std::printf ("\nMuscovite 331 at 8 KeV. Incidence at the Bragg angle:\n");
bragg = Bragg_angle (cryst, energy, 3, 3, 1);
std::printf (" Bragg angle: Rad: %f Deg: %f\n", bragg, bragg*180/PI);
q = Q_scattering_amplitude (cryst, energy, 3, 3, 1, rel_angle);
std::printf (" Q Scattering amplitude: %f\n", q);
Atomic_Factors (19, energy, q, debye_temp_factor, &f0, &fp, &fpp);
std::printf (" Atomic factors (Z = 19) f0, fp, fpp: %f, %f, i*%f\n", f0, fp, fpp);
FH = Crystal_F_H_StructureFactor (cryst, energy, 3, 3, 1, debye_temp_factor, rel_angle);
std::printf (" FH(3,3,1) structure factor: (%f, %f)\n", FH.re, FH.im);
F0 = Crystal_F_H_StructureFactor (cryst, energy, 0, 0, 0, debye_temp_factor, rel_angle);
std::printf (" F0=FH(0,0,0) structure factor: (%f, %f)\n", F0.re, F0.im);
char **crystals;
crystals = Crystal_GetCrystalsList(NULL, NULL);
std::printf ("List of available crystals:\n");
for (i = 0 ; crystals[i] != NULL ; i++) {
std::printf (" Crystal %i: %s\n", i, crystals[i]);
xrlFree(crystals[i]);
}
xrlFree(crystals);
std::printf ("\n");
/* compoundDataNIST tests */
struct compoundDataNIST *cdn;
cdn = GetCompoundDataNISTByName("Uranium Monocarbide");
std::printf ("Uranium Monocarbide\n");
std::printf (" Name: %s\n", cdn->name);
std::printf (" Density: %lf g/cm3\n", cdn->density);
for (i = 0 ; i < cdn->nElements ; i++) {
std::printf(" Element %i: %lf %%\n",cdn->Elements[i],cdn->massFractions[i]*100.0);
}
FreeCompoundDataNIST(cdn);
cdn = NULL;
cdn = GetCompoundDataNISTByIndex(NIST_COMPOUND_BRAIN_ICRP);
std::printf ("NIST_COMPOUND_BRAIN_ICRP\n");
std::printf (" Name: %s\n", cdn->name);
std::printf (" Density: %lf g/cm3\n", cdn->density);
for (i = 0 ; i < cdn->nElements ; i++) {
std::printf(" Element %i: %lf %%\n",cdn->Elements[i],cdn->massFractions[i]*100.0);
}
FreeCompoundDataNIST(cdn);
cdn = NULL;
char **nistCompounds = GetCompoundDataNISTList(NULL);
std::printf ("List of available NIST compounds:\n");
for (i = 0 ; nistCompounds[i] != NULL ; i++) {
std::printf (" Compound %i: %s\n", i, nistCompounds[i]);
xrlFree(nistCompounds[i]);
}
xrlFree(nistCompounds);
std::printf ("\n");
/* radioNuclideData tests */
struct radioNuclideData *rnd;
rnd = GetRadioNuclideDataByName("109Cd");
std::printf ("109Cd\n");
std::printf (" Name: %s\n", rnd->name);
std::printf (" Z: %i\n", rnd->Z);
std::printf (" A: %i\n", rnd->A);
std::printf (" N: %i\n", rnd->N);
std::printf (" Z_xray: %i\n", rnd->Z_xray);
std::printf (" X-rays:\n");
for (i = 0 ; i < rnd->nXrays ; i++)
std::printf (" %f keV -> %f\n", LineEnergy(rnd->Z_xray, rnd->XrayLines[i]), rnd->XrayIntensities[i]);
std::printf (" Gamma rays:\n");
for (i = 0 ; i < rnd->nGammas ; i++)
std::printf (" %f keV -> %f\n", rnd->GammaEnergies[i], rnd->GammaIntensities[i]);
FreeRadioNuclideData(rnd);
rnd = GetRadioNuclideDataByIndex(RADIO_NUCLIDE_125I);
std::printf ("RADIO_NUCLIDE_125I\n");
std::printf (" Name: %s\n", rnd->name);
std::printf (" Z: %i\n", rnd->Z);
std::printf (" A: %i\n", rnd->A);
std::printf (" N: %i\n", rnd->N);
std::printf (" Z_xray: %i\n", rnd->Z_xray);
std::printf (" X-rays:\n");
for (i = 0 ; i < rnd->nXrays ; i++)
std::printf (" %f keV -> %f\n", LineEnergy(rnd->Z_xray, rnd->XrayLines[i]), rnd->XrayIntensities[i]);
std::printf (" Gamma rays:\n");
for (i = 0 ; i < rnd->nGammas ; i++)
std::printf (" %f keV -> %f\n", rnd->GammaEnergies[i], rnd->GammaIntensities[i]);
FreeRadioNuclideData(rnd);
char **radioNuclides;
radioNuclides = GetRadioNuclideDataList(NULL);
std::printf ("List of available radionuclides:\n");
for (i = 0 ; radioNuclides[i] != NULL ; i++) {
std::printf (" Radionuclide %i: %s\n", i, radioNuclides[i]);
xrlFree(radioNuclides[i]);
}
xrlFree(radioNuclides);
std::printf ("\n--------------------------- END OF XRLEXAMPLE6 -------------------------------\n");
return 0;
}
int main(int argc, char **argv) {
Glib::set_application_name("app3");
#if defined(G_OS_WIN32)
setlocale(LC_ALL,"English_United States");
//g_setenv("LANG","en_US",TRUE);
#else
Glib::setenv("LANG","en_US",TRUE);
#endif
setbuf(stdout,NULL);
if (bam_xmlLoadCatalog() == 0) {
return 1;
}
if (xmi_xmlLoadCatalog() == 0) {
return 1;
}
static struct xmi_main_options options = xmi_get_default_main_options();
static std::string spe_file_noconv;
static std::string spe_file_conv;
static std::string csv_file_noconv;
static std::string csv_file_conv;
static std::string svg_file_noconv;
static std::string svg_file_conv;
static std::string htm_file_noconv;
static std::string htm_file_conv;
static std::string custom_detector_response;
static std::string hdf5_file;
static std::string xmimsim_hdf5_solid_angles;
static std::string xmimsim_hdf5_escape_ratios;
static double conv_threshold(BAM_QUANT_CONV_THRESHOLD_DEFAULT*100.0);
bool version = false;
std::vector<std::string> argv_files;
Glib::OptionContext option_context("inputfile XMSO-file");
option_context.set_help_enabled();
Glib::OptionGroup option_group("XMI-MSIM", "XMI-MSIM options", "Show all XMI-MSIM options");
//option_group.add_entry(App3OptionEntry( "enable-M-lines", 0, "Enable M lines (default)"), (bool&) options.use_M_lines);
option_group.add_entry(App3OptionEntry( "disable-M-lines", 0, "Disable M lines", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_M_lines);
//option_group.add_entry(App3OptionEntry( "enable-auger-cascade", 0, "Enable Auger (non radiative) cascade effects (default)"), (bool&) options.use_cascade_auger);
option_group.add_entry(App3OptionEntry( "disable-auger-cascade", 0, "Disable Auger (non radiative) cascade effects", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_cascade_auger);
//option_group.add_entry(App3OptionEntry( "enable-radiative-cascade", 0, "Enable radiative cascade effects (default)"), (bool&) options.use_cascade_radiative);
option_group.add_entry(App3OptionEntry( "disable-radiative-cascade", 0, "Disable radiative cascade effects", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_cascade_radiative);
//option_group.add_entry(App3OptionEntry( "enable-variance-reduction", 0, "Enable variance reduction (default)"), (bool&) options.use_variance_reduction);
option_group.add_entry(App3OptionEntry( "disable-variance-reduction", 0, "Disable variance reduction", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_variance_reduction);
option_group.add_entry_filename(App3OptionEntry( "with-hdf5-data", 0, "Select a HDF5 data file (advanced usage)", "file", Glib::OptionEntry::FLAG_HIDDEN), hdf5_file);
option_group.add_entry_filename(App3OptionEntry( "with-solid-angles-data", 0, "Select a HDF5 solid angles file (advanced usage)", "file", Glib::OptionEntry::FLAG_HIDDEN), xmimsim_hdf5_solid_angles);
option_group.add_entry_filename(App3OptionEntry( "with-escape-ratios-data", 0, "Select a HDF5 escape ratios file (advanced usage)", "file", Glib::OptionEntry::FLAG_HIDDEN), xmimsim_hdf5_escape_ratios);
option_group.add_entry_filename(App3OptionEntry( "spe-file", 0, "Write detector convoluted spectra to file", "file"), spe_file_conv);
option_group.add_entry_filename(App3OptionEntry( "spe-file-unconvoluted", 0, "Write detector unconvoluted spectra to file", "file"), spe_file_noconv);
option_group.add_entry_filename(App3OptionEntry( "csv-file", 0, "Write detector convoluted spectra to CSV file", "file"), csv_file_conv);
option_group.add_entry_filename(App3OptionEntry( "csv-file-unconvoluted", 0, "Write detector unconvoluted spectra to CSV file", "file"), csv_file_noconv);
option_group.add_entry_filename(App3OptionEntry( "svg-file", 0, "Write detector convoluted spectra to SVG file", "file"), svg_file_conv);
option_group.add_entry_filename(App3OptionEntry( "svg-file-unconvoluted", 0, "Write detector unconvoluted spectra to SVG file", "file"), svg_file_noconv);
option_group.add_entry_filename(App3OptionEntry( "htm-file", 0, "Write detector convoluted spectra to HTML file", "file"), htm_file_conv);
option_group.add_entry_filename(App3OptionEntry( "htm-file-unconvoluted", 0, "Write detector unconvoluted spectra to HTML file", "file"), htm_file_noconv);
option_group.add_entry(App3OptionEntry( "enable-pile-up", 0, "Enable pile-up"), (bool&) options.use_sum_peaks);
//option_group.add_entry(App3OptionEntry( "disable-pile-up", 0, "Disable pile-up (default)", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_sum_peaks);
//option_group.add_entry(App3OptionEntry( "enable-escape-peaks", 0, "Enable escape peaks (default)"), (bool&) options.use_escape_peaks);
option_group.add_entry(App3OptionEntry( "disable-escape-peaks", 0, "Disable escape peaks", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_escape_peaks);
option_group.add_entry(App3OptionEntry( "enable-poisson", 0, "Generate Poisson noise in the spectra"), (bool&) options.use_poisson);
//option_group.add_entry(App3OptionEntry( "disable-poisson", 0, "Disable the generating of spectral Poisson noise (default)", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_poisson);
option_group.add_entry(App3OptionEntry( "enable-opencl", 0, "Enable OpenCL"), (bool&) options.use_opencl);
//option_group.add_entry(App3OptionEntry( "disable-opencl", 0, "Disable OpenCL (default)", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_opencl);
option_group.add_entry(App3OptionEntry( "enable-advanced-compton", 0, "Enable advanced yet slower Compton simulation"), (bool&) options.use_advanced_compton);
//option_group.add_entry(App3OptionEntry( "disable-advanced-compton", 0, "Disable advanced yet slower Compton simulation (default)", "", Glib::OptionEntry::FLAG_REVERSE), (bool&) options.use_advanced_compton);
option_group.add_entry_filename(App3OptionEntry( "custom-detector-response", 0, "Use the supplied library for the detector response routine", "file"), custom_detector_response);
option_group.add_entry(App3OptionEntry( "set-threads", 0, "Set the number of threads (default=max)", "N"), options.omp_num_threads);
option_group.add_entry(App3OptionEntry( "verbose", 'v', "Verbose mode"), (bool&) options.verbose);
option_group.add_entry(App3OptionEntry( "very-verbose", 'V', "Even more verbose mode"), (bool&) options.extra_verbose);
option_group.add_entry(App3OptionEntry( "version", 0, "Display version information"), version);
Glib::OptionGroup option_group2("BAM-quant options", "BAM-quant options");
option_group2.add_entry(App3OptionEntry( "convergence-threshold", 0, "Convergence threshold for each of the relative X-ray intensities (default = 1%)", "VAL"), conv_threshold);
option_group2.add_entry_filename(App3OptionEntry( G_OPTION_REMAINING), argv_files);
option_context.set_main_group(option_group2);
option_context.add_group(option_group);
option_context.set_ignore_unknown_options();
try {
if (option_context.parse(argc, argv) == false) {
std::cerr << option_context.get_help(true) << std::endl;
return 1;
}
}
catch (Glib::OptionError &e) {
std::cerr << "Exception caught: " << e.what() << std::endl;
std::cerr << option_context.get_help(true) << std::endl;
return 1;
}
if (conv_threshold < 0) {
std::cerr << "convergence-threshold must be greater than 0!" << std::endl;
return 1;
}
//convert to fraction
conv_threshold /= 100.0;
if (version) {
std::cout << xmi_version_string() << std::endl;
return 0;
}
if (argv_files.size() != 2) {
std::cerr << option_context.get_help(true) << std::endl;
return 1;
}
if (options.extra_verbose) {
options.verbose = 1;
std::cout << "Option M-lines: " << options.use_M_lines << std::endl;
std::cout << "Option non-radiative cascade: " << options.use_cascade_auger << std::endl;
std::cout << "Option radiative cascade: " << options.use_cascade_radiative << std::endl;
std::cout << "Option variance reduction: " << options.use_variance_reduction << std::endl;
std::cout << "Option pile-up: " << options.use_sum_peaks << std::endl;
std::cout << "Option Poisson noise: " << options.use_poisson << std::endl;
std::cout << "Option escape peaks: " << options.use_escape_peaks << std::endl;
std::cout << "Option OpenCL: " << options.use_opencl << std::endl;
std::cout << "Option number of threads: " << options.omp_num_threads << std::endl;
std::cout << "Option advanced Compton: " << options.use_advanced_compton << std::endl;
}
try {
BAM::Job::XMSI::RandomNumberAcquisitionStart();
//disable xraylib's error messages
SetErrorMessages(0);
//start by reading in the inputfile
BAM::File::RXI::Common *common = BAM::File::RXI::Parse(argv_files[0]);
BAM::Job::Quant job(common, argv_files[1], options, conv_threshold);
BAM::Job::XMSI::RandomNumberAcquisitionStop();
}
catch (BAM::Exception &e) {
std::cerr << "Fatal BAM exception: " << e.what() << std::endl;
return 1;
}
catch (std::exception &e) {
std::cerr << "Fatal non-BAM exception: " << e.what() << std::endl;
return 1;
}
catch (Glib::Exception &e) {
std::cerr << "Fatal Glib exception: " << e.what() << std::endl;
return 1;
}
catch (...) {
std::cerr << "Unknown exception occurred" << std::endl;
return 1;
}
return 0;
}
