Alexander Ziegler / gammalib

…s only used for Monte Carlo simulation and has a quite specific behaviour.

…y intergation is done over all energy ranges in an energy boundary object (see #1450).

Updated at the same time the Doxygen information.

…ll GModelData derivatives

The CTA response was refactored. GCTAResponse was renamed GCTAResponseIrf, and GCTAResponse became an abstract base class for all possible CTA response functions. GCTAResponseIrf derives from GCTAResponse. Also GCTAResponseCube derives from GCTAResponse.

To make this work, the GObservation::response() method now returns a pointer instead of a reference. The interface definition is such that the pointer never can be NULL. The code has been adapted to comply to this new interface.

Upon reading an XML file, GCTAObservation now determines the appropriate reponse type and calls its reader. If no response is specified in the XML file, GCTAResponseIrf will be used by default.

The norm() method has been added to the abstract GModelSpatial interface that returns the model normalization within a circular region. This is used in the GModelSky::mc method to correctly reflect the normalization specified in the XML file. The method is also used to flag sources outside the simulation cone (applies so for the point sources, radial source and elliptical sources).

Minor adjustments were needed to make the code compile.

The energy integration has been implemented in the GModelSky::integrate_energy method. For that purpose, the GModelSky::edisp_kern class has been added to GModelSky that provides the kernel for the integration.

The GResponse::src_ebounds method has been added to determine the true energy boundaries for a given observed event energy. The method is virtual in GResponse and always throws a "feature not implemented" exception. The method has then also been implemented in GCTAResponse to provide the true energy boundaries.

This change was prompted by the addition of the ctobssim option
to apply energy dispersion. Even if an instrument has an energy
dispersion, it is useful to be able to apply it or not.

So the meaning now is roughly this:
use_edisp = (has_edisp && apply_edisp)
although the apply_edisp flag has so far been only added to
GCTAResponse, use_edisp = False.

GResponse::has_tdisp() was also renamed to GResponse::use_tdisp()
for consistency.

The attribute for identifiers is "id". This has been corrected in GCTAModelBackground. Further, writing of the attribute has been added to GModelSky. Attributes are also only written if the identifier string is not empty.

The coding and design conventions (CDC) document was also updated.

The abstract pointing class has been removed as so far no abstract pointing interface is needed.

…ling to the simulation cone.

This assures that the simulated spectrum corresponds indeed to the selected zone in the diffuse map cube, and also considerably speeds up computations for large maps.

The mc() method now explicitely support Monte Carlo simulation for diffuse map cubes. The method detects whether a spatial component is of type GModelSpatialDiffuseCube. In case, it takes the diffuse map cube spectrum as the spectral component, implemented by a GModelSpectralNodes object. Any existing spectral component is multiplied to the nodes model by evaluation the spectral function at the node energies.

Note that the method still has a limitation in that the simulation cone is not taken into account when extracting the diffuse map cube spectrum.

Add GChatter argument to print() method.

The print() method now takes an argument of type GChatter to control the chattiness of the output in a homogenous manner.
Pass NORMAL chatter level in __str__ method.
Add chatter argument to print() methods.
Use GLog::chatter() method to control the chattiness of the log stream operator.
Add chatter() methods and provide member access through const references.
Add chatter() methods, provide member access through const references, and declare one liners as inline methods.

The chatter() methods control the new m_chatter attribute that is used to control the chattiness of the object's print() methods.
name() and version() methods return const references.
logFileOpen() method now sets the chattiness from the applications's chatter parameter, and name() and version() methods now return a const reference.