/***********************************************************************//**
* @brief Return number of bins in mean PSF
*
* @return Number of bins in mean PSF
***************************************************************************/
inline
int GLATMeanPsf::size(void) const
{
// Compute size
int size = m_energy.size()*m_offset.size();
// Return size
return size;
}
/***********************************************************************//**
* @brief Return maximum delta value in radians
*
* @return Maximum delta value (radians).
***************************************************************************/
inline
double GCTACubePsf::delta_max(void) const
{
// Get maximum delta value
double delta_max = (m_deltas.size() > 0) ? m_deltas[m_deltas.size()-1] : 0.0;
// Return
return (delta_max * gammalib::deg2rad);
}
/***********************************************************************//**
* @brief Update 1D cache
*
* @param[in] arg Argument.
*
* Updates the 1D interpolation cache. The interpolation cache is composed
* of two indices and weights that define 2 data values of the 2D table
* that are used for linear interpolation.
*
* @todo Write down formula
*
* @todo Makes GNodeArray::set_value method const and use mutable members
***************************************************************************/
void GCTAResponseTable::update(const double& arg) const
{
// Get pointer to node array (circumvent const correctness)
GNodeArray* nodes = const_cast<GNodeArray*>(&(m_axis_nodes[0]));
// Set value for node array
nodes->set_value(arg);
// Set indices and weighting factors for interpolation
m_inx_left = nodes->inx_left();
m_inx_right = nodes->inx_right();
m_wgt_left = nodes->wgt_left();
m_wgt_right = nodes->wgt_right();
// Return
return;
}
/***********************************************************************//**
* @brief Return number of node energies in response
*
* @return Number of node energies.
***************************************************************************/
inline
int GCTABackgroundPerfTable::size(void) const
{
return (m_logE.size());
}
/***********************************************************************//**
* @brief Return number of energy bins
*
* @return Number of energy bins.
***************************************************************************/
inline
int GLATMeanPsf::nenergies(void) const
{
return m_energy.size();
}
/***********************************************************************//**
* @brief Return number of offset bins
*
* @return Number of offset bins.
***************************************************************************/
inline
int GLATMeanPsf::noffsets(void) const
{
return m_offset.size();
}
/***********************************************************************//**
* @brief Return map offset
*
* @return Map offset.
***************************************************************************/
inline
int GCTACubePsf::offset(const int& idelta, const int& iebin) const
{
return (idelta + iebin*m_deltas.size());
}
