Add CTA background cube skymap grid constructor

Jürgen Knödlseder
1 parent 571d0bc5
Showing 6 changed files with 68 additions and 7 deletions
inst/cta/include/GCTACubeBackground.hpp
inst/cta/pyext/GCTACubeBackground.i
inst/cta/src/GCTACubeBackground.cpp
inst/cta/src/GCTACubeEdisp.cpp
inst/cta/src/GCTACubeExposure.cpp
inst/cta/src/GCTACubePsf.cpp
@@ -55,6 +55,15 @@ public:
     explicit GCTACubeBackground(const GFilename& filename);
     explicit GCTACubeBackground(const GCTAEventCube& cube);
     GCTACubeBackground(const GCTACubeBackground& bgd);
+    GCTACubeBackground(const std::string&   wcs,
+                       const std::string&   coords,
+                       const double&        x,
+                       const double&        y,
+                       const double&        dx,
+                       const double&        dy,
+                       const int&           nx,
+                       const int&           ny,
+                       const GEnergies&     energies);
     virtual ~GCTACubeBackground(void);
     // Operators
@@ -42,6 +42,15 @@ public:
     explicit GCTACubeBackground(const GFilename& filename);
     explicit GCTACubeBackground(const GCTAEventCube& cube);
     GCTACubeBackground(const GCTACubeBackground& bgd);
+    GCTACubeBackground(const std::string&   wcs,
+                       const std::string&   coords,
+                       const double&        x,
+                       const double&        y,
+                       const double&        dx,
+                       const double&        dy,
+                       const int&           nx,
+                       const int&           ny,
+                       const GEnergies&     energies);
     virtual ~GCTACubeBackground(void);
     // Operators
@@ -158,6 +158,49 @@ GCTACubeBackground::GCTACubeBackground(const GCTACubeBackground&amp; bgd)
 /***********************************************************************//**
+ * @brief Background cube constructor
+ *
+ * @param[in] wcs      World Coordinate System.
+ * @param[in] coords   Coordinate System (CEL or GAL).
+ * @param[in] x        X coordinate of sky map centre (deg).
+ * @param[in] y        Y coordinate of sky map centre (deg).
+ * @param[in] dx       Pixel size in x direction at centre (deg/pixel).
+ * @param[in] dy       Pixel size in y direction at centre (deg/pixel).
+ * @param[in] nx       Number of pixels in x direction.
+ * @param[in] ny       Number of pixels in y direction.
+ * @param[in] energies Energies.
+ *
+ * Constructs a background cube by specifying the sky map grid and the
+ * energies.
+ ***************************************************************************/
+GCTACubeBackground::GCTACubeBackground(const std::string&   wcs,
+                                       const std::string&   coords,
+                                       const double&        x,
+                                       const double&        y,
+                                       const double&        dx,
+                                       const double&        dy,
+                                       const int&           nx,
+                                       const int&           ny,
+                                       const GEnergies&     energies)
+{
+    // Initialise class members
+    init_members();
+
+    // Store energies
+    m_energies = energies;
+
+    // Set GNodeArray used for interpolation
+    set_eng_axis();
+
+    // Create sky map
+    m_cube = GSkyMap(wcs, coords, x, y, dx, dy, nx, ny, m_energies.size());
+
+    // Return
+    return;
+}
+
+
+/***********************************************************************//**
  * @brief Destructor
  ***************************************************************************/
 GCTACubeBackground::~GCTACubeBackground(void)
@@ -202,7 +202,7 @@ GCTACubeEdisp::GCTACubeEdisp(const GCTAEventCube&amp; cube, const double&amp; mmax,
 /***********************************************************************//**
- * @brief Mean PSF cube constructor
+ * @brief Energy dispersion cube constructor
  *
  * @param[in] wcs      World Coordinate System.
  * @param[in] coords   Coordinate System (CEL or GAL).
@@ -219,7 +219,8 @@ GCTACubeEdisp::GCTACubeEdisp(const GCTAEventCube&amp; cube, const double&amp; mmax,
  * @exception GException::invalid_argument
  *            Maximum energy migration or number of migration bins invalid.
  *
- * Construct an energy dispersion cube with all elements set to zero.
+ * Constructs an energy dispersion cube by specifying the sky map grid and
+ * the energies.
  ***************************************************************************/
 GCTACubeEdisp::GCTACubeEdisp(const std::string&   wcs,
                              const std::string&   coords,
@@ -161,8 +161,8 @@ GCTACubeExposure::GCTACubeExposure(const GCTAEventCube&amp; cube)
  * @param[in] ny       Number of pixels in y direction.
  * @param[in] energies Energies.
  *
- * Constructs an exposure cube by specifying the sky map grid and the energy
- * boundaries.
+ * Constructs an exposure cube by specifying the sky map grid and the
+ * energies.
  ***************************************************************************/
 GCTACubeExposure::GCTACubeExposure(const std::string&   wcs,
                                    const std::string&   coords,
@@ -181,7 +181,7 @@ GCTACubePsf::GCTACubePsf(const GCTAEventCube&amp; cube, const double&amp; dmax,
 /***********************************************************************//**
- * @brief Mean PSF cube constructor
+ * @brief PSF cube constructor
  *
  * @param[in] wcs      World Coordinate System.
  * @param[in] coords   Coordinate System (CEL or GAL).
@@ -195,8 +195,7 @@ GCTACubePsf::GCTACubePsf(const GCTAEventCube&amp; cube, const double&amp; dmax,
  * @param[in] dmax     Maximum delta (deg).
  * @param[in] ndbins   Number of delta bins.
  *
- * Constructs a mean PSF cube by computing the mean PSF from all CTA
- * observations found in the observation container.
+ * Constructs a PSF cube by specifying the sky map grid and the energies.
  ***************************************************************************/
 GCTACubePsf::GCTACubePsf(const std::string&   wcs,
                          const std::string&   coords,