NEWS
6.6 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
New Features and Important Changes in GammaLib 0.11.0
Author(s) : Juergen Knoedlseder
28 October 2015
1. Introduction
---------------
This document briefly summarizes module by module the most important changes
made since the last release of GammaLib.
2. Public interface modifications
---------------------------------
The GOptimizerLM logger constructor now takes a pointer to a logger instead of
a reference. This allows setting of a NULL pointer, disabling the logging of the
optimizer.
The GSkymap class has been renamed to GSkyMap to be consistent with the other
classes in the sky module.
The small angle approximation has been removed from GModelSpatialRadialShell.
3. Configuration
-----------------
None
4. Application module
---------------------
GammaLib applications now only look for the application parameter file in the
users pfiles directory if the PFILES environment variable has been not set. This
assure that the PFILES environment variable will always take precedence over the
users pfiles directory, and if the paramter file is missing in the directory
pointed by the PFILES environment variable, it will be copied over from the
syspfiles directory.
The GApplication class now has a need_help() method that returns true is the
--help option is specified on the command line. This allows the client to display
any help information.
5. Base module
--------------
None
6. FITS I/O Interface
---------------------
None
7. Linear algebra module
------------------------
None
8. Model module
---------------
Implement a rejection method for Monte Carlo simulations of photons in
GModelSpatialDiffuseMap and GModelSpatialDiffuseCube. This assures that the
simulated photon distribution follows the intensity distribution of the skymap
as obtained using the sky direction access operator.
The algorithm used to draw random numbers in GModelSpatialDiffuseConst::mc()
was wrong (sky directions were not distributed isotropically on the sky). This
has now been fixed (#1480).
The small angle approximation has been removed from GModelSpatialRadialShell
as this complexifies the class and the approximation is not really needed.
A GModelSpatialDiffuseMap::set_mc_cone method has been introduced to handle
correctly Monte Carlo simulations for unnormalized diffuse maps. This method
now computes the normalization for a map within the acceptance cone which is
then returned by the GModelSpatialDiffuseMap::mc_norm method.
The GModelSpatial::norm method has been renamed to GModelSpatial::mc_norm so that
the usage of this method is now explicitly given by the method name.
save() and write() methods have been added to GModelSpatialDiffuseCube that allow
the saving of a diffuse cube spatial model to disk or into a FITS objects.
The computations in GModelSpatialDiffuseCube were not precise enough to allow
an accurate modelling of the simulated events. The eval() and eval_gradients()
methods now use a log-log interpolation. Furthermore, the map is selected
in mc() randomly based on the interpolation weights of th neighbouring maps.
And some bugs related to the map indexing and maximum map value computation
have been corrected.
9. Numerics module
------------------
None
10. Observation module
----------------------
Integration in npred() is now done using a fixed number of Romberg iterations
instead of using the integration precision. This will assure that the log-
likelihood function is smooth.
The number of integration iterations for the energy dispersion was changed from
4 to 5 to make it compatible with the default integration order.
An interval constructor and the set_lin() and set_log() methods have been added
to the GEnergies class to allow quick creation of a linear or logarithmic series
of energies.
11. Optimization module
-----------------------
The gradient computation has been corrected. Instead of multiplying the gradient
factor with the scale factor we have to divide the gradient factor by the scale
factor.
An logger() method has been added to GOptimizerLM to set the logger for an
optimizer (before, a logger could only be set during construction).
12. Sky module
--------------
Added the flux() method that provides estimates for the flux value of a pixel. The
flux is defined as the intensity (computed by bilinear interpolation over the pixel)
multiplied by the solid angle of the pixel.
The pixel containment check has been corrected. Before, the containment check returned
true even if the pixel was by one beyond the skymap limits. On of the unwanted side
effects in the sky map access operator was that pixels outside the valid sky map
array were accessed, leading to unpredictable results at the edge of the images.
This quite severe bug was now corrected.
The GSkymap::operator/=(GSkymap) operator now sets all pixels to zero for which the
nominator is zero. No exception is thrown any longer.
Add +, -, * and / binary operators to the GSkymap class.
13. Support module
------------------
None
14. XML module
--------------
None
15. VO module
-------------
None
16. Xspec module
----------------
None
17. Test module
---------------
None
18. Multi-wavelength interface
------------------------------
None
19. Fermi-LAT interface
----------------------
None
20. CTA interface
-----------------
All integrations are now done using a fixed number of Romberg iterations
instead of using the integration precision. This will assure that the log-
likelihood function is smooth.
The number of integration iterations for the energy dispersion was changed from
4 to 5 to make it compatible with the default integration order.
An is_valid() method has been added to GCTARoi that returns true if the ROI
has a positive radius. DS keywords are only written for a valid ROI into the
event list file.
The OBJECT attribute of an event file is now stored in to m_object member of
GCTAObservation, and object() methods have been added to set and get the
attribute value.
The deadtime correction has been removed from GCTAModelCubeBackground as the
background cube rates are already given per livetime. We formerly did the
correction twice (#1503).
The GCTAModelAeffBackground class has been added to generate a background model
from the effective area information (#1525). The usage of this model is similar
to that of the IRF background model. The model type in the XML file is
"CTAAeffBackground".
The precision for the IRF computation of radial shell models has been increased
to reduce a bias in the fitting of the model for unbinned analysis (#1521). There
is still a slight bias for the shell parameters in an unbinned analysis when the
width of the shell becomes comparable to the angular resolution.
21. COMPTEL interface
---------------------
None