Line data    Source code

       1              : #include "crpropa/massDistribution/Cordes.h"
       2              : #include "crpropa/Common.h"
       3              : 
       4              : #include <sstream>
       5              : 
       6              : namespace crpropa {
       7              : 
       8            4 : double Cordes::getHIIDensity(const Vector3d &position) const {
       9              :         double n=0;  // in m^-3
      10              : 
      11            4 :         double z=position.z;
      12            4 :         double R = sqrt(position.x*position.x+position.y*position.y);  //radius in galactic disk
      13              : 
      14            4 :         n += 0.025/ccm*exp(-fabs(z)/(1*kpc))*exp(-pow_integer<2>(R/(20*kpc)));  // galactocentric component
      15            4 :         n += 0.2/ccm*exp(-fabs(z)/(0.15*kpc))*exp(-pow_integer<2>((R-4*kpc)/(2*kpc)));  // anular component
      16              : 
      17            4 :         return n;
      18              : }
      19              : 
      20            2 : double Cordes::getDensity(const Vector3d &position) const {
      21            2 :         return Cordes::getHIIDensity(position);
      22              : }
      23              : 
      24            1 : double Cordes::getNucleonDensity(const Vector3d &position) const {
      25            1 :         return getHIIDensity(position);
      26              : }
      27              : 
      28            1 : bool Cordes::getIsForHI() {
      29            1 :         return isforHI;
      30              : }
      31              : 
      32            1 : bool Cordes::getIsForHII() {
      33            1 :         return isforHII;
      34              : }
      35              : 
      36            1 : bool Cordes::getIsForH2() {
      37            1 :         return isforH2;
      38              : }
      39              : 
      40            0 : std::string Cordes::getDescription() {
      41            0 :         std::stringstream s;
      42            0 :         s << "Density Cordes: includes only an HII component";
      43            0 :         return s.str();
      44            0 : }
      45              : 
      46              : }  // namespace crpropa