GammaLib

Martina Cardillo reported a problem during the event simulations using some custom IRFs. For some reason, the simulated energies are outside the energy range of interest. Here the error message that occurs with the attached test script:

ValueError: *** ERROR in GCTABackground3D::mc(GEnergy&, GTime&, GRan&): Invalid value. Event energy 395.910126393249 TeV is outside the energy range [630.957365036011 GeV, 350.392200724172 TeV] covered by the background response table. Please restrict the energy range of the simulation to the validity range of the background response table.

It turned out that the number of internal nodes were not correctly computed in GCTABackground3D::init_mc_cache. Instead of

double logE_bin = (m_logE_max - m_logE_min) / (nbins+1);

the line should be

double logE_bin = (m_logE_max - m_logE_min) / (nbins-1);