added core_radius set and get functions

nkelhos
1 parent 9b4007e4
Showing 6 changed files with 100 additions and 0 deletions
include/GModelSpatialRadialProfileDMBurkert.hpp
include/GModelSpatialRadialProfileDMEinasto.hpp
include/GModelSpatialRadialProfileDMZhao.hpp
pyext/GModelSpatialRadialProfileDMBurkert.i
pyext/GModelSpatialRadialProfileDMEinasto.i
pyext/GModelSpatialRadialProfileDMZhao.i
@@ -76,6 +76,8 @@ public:
     void   scale_density(const double& scale_density);
     double halo_distance(void) const ;
     void   halo_distance(const double& halo_distance);
+    double core_radius(void) const ;
+    void   core_radius(const double& core_radius);
     double mass_density( const double& radius ) const;
     double jfactor( const double& angle ) const;
  
@@ -230,4 +232,34 @@ void GModelSpatialRadialProfileDMBurkert::halo_distance(const double&amp; halo_dista
     return;
 }
  
+/***********************************************************************//**
+ * @brief Return Einasto core radius 
+ *
+ * @return core_radius (unitless).
+ *
+ * Returns the core radius in the Burkert halo density function. The DM
+ * mass density inside this radius will be the same as at this radius.
+ ***************************************************************************/
+inline
+double GModelSpatialRadialProfileDMBurkert::core_radius(void) const
+{
+    return (m_core_radius.value());
+}
+
+/***********************************************************************//**
+ * @brief Set Burkert core radius
+ *  
+ * @param[in] core radius (unitless).
+ *
+ * Sets the Burkert core radius, should be positive.  The DM mass density 
+ * inside this radius will be the same as at this radius.
+ ***************************************************************************/
+inline
+void GModelSpatialRadialProfileDMBurkert::core_radius(const double& core_radius)
+{
+    m_core_radius.value(core_radius);
+    return;
+}
+
+
 #endif /* GMODELSPATIALRADIALPROFILEDMBURKERT_HPP */
@@ -78,6 +78,8 @@ public:
     void   halo_distance(const double& halo_distance);
     double alpha(void) const;
     void   alpha(const double& alpha);
+    double core_radius(void) const ;
+    void   core_radius(const double& core_radius);
     double mass_density( const double& radius ) const ;
     double jfactor( const double& angle ) const ;
  
@@ -262,5 +264,34 @@ void GModelSpatialRadialProfileDMEinasto::alpha(const double&amp; alpha)
     return;
 }
  
+/***********************************************************************//**
+ * @brief Return Einasto core radius 
+ *
+ * @return core_radius (unitless).
+ *
+ * Returns the core radius in the Einasto halo density function.
+ ***************************************************************************/
+inline
+double GModelSpatialRadialProfileDMEinasto::core_radius(void) const
+{
+    return (m_core_radius.value());
+}
+
+/***********************************************************************//**
+ * @brief Set Einasto core radius
+ *  
+ * @param[in] core radius (unitless).
+ *
+ * Sets the Einasto core radius, should be positive.  The DM mass density
+ * will be the same inside this radius as at this radius.
+ ***************************************************************************/
+inline
+void GModelSpatialRadialProfileDMEinasto::core_radius(const double& core_radius)
+{
+    m_core_radius.value(core_radius);
+    return;
+}
+
+
 #endif /* GMODELSPATIALRADIALPROFILEDMEINASTO_HPP */
  
@@ -83,6 +83,8 @@ public:
     void   beta(const double& beta) ;
     double gamma(void) const ;
     void   gamma(const double& gamma) ;
+    double core_radius(void) const ;
+    void   core_radius(const double& core_radius) ;
     double mass_density(const double& radius) const;
     double jfactor(const double& angle) const ;
  
@@ -330,4 +332,33 @@ void GModelSpatialRadialProfileDMZhao::gamma(const double&amp; gamma)
     return;
 }
  
+/***********************************************************************//**
+ * @brief Return Zhao core radius 
+ *
+ * @return core_radius (unitless).
+ *
+ * Returns the core radius in the Zhao halo density function.
+ ***************************************************************************/
+inline
+double GModelSpatialRadialProfileDMZhao::core_radius(void) const
+{
+    return (m_core_radius.value());
+}
+
+/***********************************************************************//**
+ * @brief Set Zhao core radius
+ *  
+ * @param[in] core radius (unitless).
+ *
+ * Sets the Zhao core radius, should be positive.  The DM mass density
+ * will be the same inside this radius as at this radius.
+ ***************************************************************************/
+inline
+void GModelSpatialRadialProfileDMZhao::core_radius(const double& core_radius)
+{
+    m_core_radius.value(core_radius);
+    return;
+}
+
+
 #endif /* GMODELSPATIALRADIALPROFILEDMZHAO_HPP */
@@ -63,6 +63,8 @@ public:
     void   scale_density(const double& scale_density);
     double halo_distance(void) const ;
     void   halo_distance(const double& halo_distance);
+    double core_radius(void) const ;
+    void   core_radius(const double& core_radius);
     double mass_density( const double& radius ) const;
     double jfactor( const double& angle ) const;
 };
@@ -65,6 +65,8 @@ public:
     void   halo_distance(const double& halo_distance);
     double alpha(void) const;
     void   alpha(const double& alpha);
+    double core_radius(void) const;
+    void   core_radius(const double& core_radius);
     double mass_density( const double& radius ) const ;
     double jfactor( const double& angle ) const ;
 };
@@ -69,6 +69,8 @@ public:
     void   beta(const double& beta) ;
     double gamma(void) const ;
     void   gamma(const double& gamma) ;
+    double core_radius(void) const ;
+    void   core_radius(const double& core_radius) ;
     //double prof_val(const double& theta);
     double mass_density(const double& radius) const;
     double jfactor(const double& angle) const ;