GObservation.i
5.85 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
/***************************************************************************
* GObservation.i - Abstract observation base class *
* ----------------------------------------------------------------------- *
* copyright (C) 2008-2021 by Juergen Knoedlseder *
* ----------------------------------------------------------------------- *
* *
* This program is free software: you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or *
* (at your option) any later version. *
* *
* This program is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
* GNU General Public License for more details. *
* *
* You should have received a copy of the GNU General Public License *
* along with this program. If not, see <http://www.gnu.org/licenses/>. *
* *
***************************************************************************/
/**
* @file GObservation.i
* @brief Abstract observation base class interface definition
* @author Juergen Knoedlseder
*/
%{
/* Put headers and other declarations here that are needed for compilation */
#include "GObservation.hpp"
#include "GEvents.hpp"
#include "GEventList.hpp"
#include "GEventCube.hpp"
#include "GTools.hpp"
%}
/* __ Typemaps ___________________________________________________________ */
/*
* Convert double output argument from and into a Python float object. This
* typemap is needed for the likelihood() method.
*/
%typemap(in) double *npred (double temp) {
temp = PyFloat_AsDouble($input);
$1 = &temp;
}
%typemap(argout) double *npred {
$result = PyFloat_FromDouble(*$1);
}
/***********************************************************************//**
* @class GObservation
*
* @brief Abstract observation base class
***************************************************************************/
class GObservation : public GBase {
public:
// Constructors and destructors
GObservation(void);
GObservation(const GObservation& obs);
virtual ~GObservation(void);
// Pure virtual methods
virtual void clear(void) = 0;
virtual GObservation* clone(void) const = 0;
virtual std::string classname(void) const = 0;
virtual void response(const GResponse& rsp) = 0;
virtual const GResponse* response(void) const = 0;
virtual std::string instrument(void) const = 0;
virtual double ontime(void) const = 0;
virtual double livetime(void) const = 0;
virtual double deadc(const GTime& time) const = 0;
virtual void read(const GXmlElement& xml) = 0;
virtual void write(GXmlElement& xml) const = 0;
// Virtual methods
virtual GEvents* events(void);
virtual void events(const GEvents& events);
virtual double likelihood(const GModels& models,
GVector* gradients,
GMatrixSparse* curvature,
double* npred) const;
virtual double model(const GModels& models,
const GEvent& event,
GVector* gradients = NULL) const;
virtual GVector model(const GModels& models,
GMatrixSparse* gradients = NULL) const;
virtual int nobserved(void) const;
virtual double npred(const GModels& models,
GVector* gradients = NULL) const;
virtual double npred(const GModel& model) const;
virtual double model_grad(const GModel& model,
const GModelPar& par,
const GEvent& event) const;
virtual double npred_grad(const GModel& model,
const GModelPar& par) const;
virtual void remove_response_cache(const std::string& name);
// Implemented methods
bool has_events(void) const;
bool has_gradient(const GModel& model,
const GModelPar& par) const;
void name(const std::string& name);
void id(const std::string& id);
void statistic(const std::string& statistic);
const std::string& name(void) const;
const std::string& id(void) const;
const std::string& statistic(void) const;
void computed_gradient(const GModel& model,
const GModelPar& par) const;
};
/***********************************************************************//**
* @brief GObservation class extension
***************************************************************************/
%extend GObservation {
%pythoncode {
def __getstate__(self):
if self.has_events():
events = self.events()
else:
events = None
args = (self.name(), self.id(), self.statistic(), events)
return args
def __setstate__(self, state):
self.__init__()
self.name(state[0])
self.id(state[1])
self.statistic(state[2])
if state[3] != None:
self.events(state[3])
}
};