added method to extract backscal as GNdarray

L. Tibaldo
1 parent 53177ac4
Showing 3 changed files with 21 additions and 0 deletions
include/GPha.hpp
pyext/GPha.i
src/xspec/GPha.cpp
@@ -101,6 +101,7 @@ public:
     const double&    areascal(const int& index) const;
     void             backscal(const int& index, const double& backscal);
     const double&    backscal(const int& index) const;
+    GNdarray         backscal_spectrum(void) const;
     const double&    underflow(void) const;
     const double&    overflow(void) const;
     const double&    outflow(void) const;
@@ -62,6 +62,7 @@ public:
     const double&    areascal(const int& index) const;
     void             backscal(const int& index, const double& backscal);
     const double&    backscal(const int& index) const;
+    GNdarray         backscal_spectrum(void) const;
     const double&    underflow(void) const;
     const double&    overflow(void) const;
     const double&    outflow(void) const;
@@ -601,6 +601,25 @@ GNdarray GPha::counts_spectrum(void) const
     return counts;
 }
  
+/***********************************************************************//**
+ * @brief Number of counts in spectrum as GNdArray
+ *
+ * Returns the number of counts in the spectrum as a GNdArray.
+ ***************************************************************************/
+GNdarray GPha::backscal_spectrum(void) const
+{
+    // Initialise array
+    int size = m_backscal.size();
+    GNdarray alpha = GNdarray(size);
+    
+    // Compute content
+    for (int i = 0; i < size; ++i) {
+      alpha(i) = m_backscal[i];
+    }
+
+    // Return counts
+    return alpha;
+}
  
 /***********************************************************************//**
  * @brief Fill spectrum with a value.