ctobservation.hpp
6.94 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
/***************************************************************************
* ctobservation - Base class for observation tools *
* ----------------------------------------------------------------------- *
* copyright (C) 2016-2019 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 ctobservation.hpp
* @brief Observation tool base class interface definition
* @author Juergen Knoedlseder
*/
#ifndef CTOBSERVATION_HPP
#define CTOBSERVATION_HPP
/* __ Includes ___________________________________________________________ */
#include "ctool.hpp"
/* __Definitions _________________________________________________________ */
/***********************************************************************//**
* @class ctobservation
*
* @brief Base class for observation tools
*
* This is the baseclass for observation tools. Observation tools are ctools
* that hold an observation container.
***************************************************************************/
class ctobservation : public ctool {
public:
// Constructors and destructors
ctobservation(const std::string& name,
const std::string& version);
ctobservation(const std::string& name,
const std::string& version,
const GApplicationPars& pars);
ctobservation(const std::string& name,
const std::string& version,
int argc,
char* argv[]);
ctobservation(const std::string& name,
const std::string& version,
const GObservations& obs);
ctobservation(const ctobservation& app);
virtual ~ctobservation(void);
// Operators
ctobservation& operator=(const ctobservation& app);
// Pure virtual methods
virtual void clear(void) = 0;
virtual void run(void) = 0;
virtual void save(void) = 0;
// Methods
void obs(const GObservations& obs);
const GObservations& obs(void) const;
#ifndef SWIG
protected:
#endif
// Protected methods
GCTAObservation* first_unbinned_observation(void);
GCTAObservation* next_unbinned_observation(void);
const GCTAObservation* first_unbinned_observation(void) const;
const GCTAObservation* next_unbinned_observation(void) const;
void read_ogip_keywords(GFitsHDU* hdu) const;
void write_ogip_keywords(GFitsHDU* hdu) const;
void set_obs_statistic(const std::string& statistic);
void set_obs_bounds();
void save_events_fits(void);
void save_events_xml(void);
protected:
// Protected methods
void init_members(void);
void copy_members(const ctobservation& app);
void free_members(void);
// Protected members
GObservations m_obs; //!< Observation container
private:
// Private members
mutable std::string m_ogip_telescope; //!< Name of telescope
mutable GTime m_ogip_tstart; //!< Start time for OGIP keywords
mutable GTime m_ogip_tstop; //!< Stop time for OGIP keywords
mutable double m_ogip_telapse; //!< Elapsed time
mutable double m_ogip_exposure; //!< Exposure time
mutable double m_ogip_ontime; //!< Ontime for OGIP keywords
mutable double m_ogip_livetime; //!< Livetime for OGIP keywords
mutable int m_index_unbinned; //!< Current index of unbinned observation
};
/***********************************************************************//**
* @brief Set observation container
*
* @param[in] obs Observation container.
*
* Set observation container.
***************************************************************************/
inline
void ctobservation::obs(const GObservations& obs)
{
m_obs = obs;
return;
}
/***********************************************************************//**
* @brief Return observation container
*
* @return Reference to observation container.
*
* Returns a reference to the observation container.
***************************************************************************/
inline
const GObservations& ctobservation::obs(void) const
{
return m_obs;
}
/***********************************************************************//**
* @brief Return first unbinned CTA observation
*
* @return Pointer to first unbinned CTA observation
*
* Returns a pointer to the first unbinned CTA observation in the container.
* If no CTA observation exists a NULL pointer is returned.
*
* The method calls next_unbinned_observation(). See the method for details.
***************************************************************************/
inline
GCTAObservation* ctobservation::first_unbinned_observation(void)
{
return (const_cast<GCTAObservation*>
(static_cast<const ctobservation&>
(*this).first_unbinned_observation()));
}
/***********************************************************************//**
* @brief Return next unbinned CTA observation
*
* @return Pointer to next unbinned CTA observation
*
* Returns a pointer to the next unbinned CTA observation in the container.
* If no CTA observation exists any more a NULL pointer is returned.
*
* The method writes for each encountered observation a level 3 header into
* the logger. It will also signal when an observation was skipped because
* it either was not a CTA observation or not an unbinned observation.
***************************************************************************/
inline
GCTAObservation* ctobservation::next_unbinned_observation(void)
{
return (const_cast<GCTAObservation*>
(static_cast<const ctobservation&>
(*this).next_unbinned_observation()));
}
#endif /* CTOBSERVATION_HPP */