Line data    Source code

       1              : #include "crpropa/massDistribution/ConstantDensity.h"
       2              : 
       3              : #include "kiss/logger.h"
       4              : 
       5              : #include <sstream>
       6              : 
       7              : namespace crpropa{
       8              : 
       9            3 : ConstantDensity::ConstantDensity(double HI, double HII, double H2) {
      10              :         // set all types active which are not equal 0 and change number density
      11            3 :         if(HI!=0)
      12            3 :                 setHI(true, HI);
      13            3 :         if(HII!=0)
      14            3 :                 setHII(true, HII);
      15            3 :         if(H2!=0)
      16            3 :                 setH2(true, H2);
      17            3 : }
      18              : 
      19            4 : double ConstantDensity::getDensity(const Vector3d &position) const {
      20              :         double n = 0;
      21              : 
      22            4 :         if(isHI)
      23            3 :                 n += HIdensitynumber;
      24            4 :         if(isHII)
      25            3 :                 n += HIIdensitynumber;
      26            4 :         if(isH2)
      27            3 :                 n += H2densitynumber;
      28              : 
      29              :         // check if all densities are deactivated and raise warning if so
      30            4 :         if((isHI || isHII || isH2) == false){
      31            2 :                 KISS_LOG_WARNING
      32            1 :                         << "\nCalled getNucleonDensity on fully deactivated ConstantDensity "
      33            1 :                         << "gas density model. In this case the density is allways set to 0. \n";
      34              :         }
      35              : 
      36            4 :         return n;
      37              : }
      38              : 
      39            4 : double ConstantDensity::getNucleonDensity(const Vector3d &position) const {
      40              :         double n = 0;
      41              : 
      42            4 :         if(isHI)
      43            3 :                 n += HIdensitynumber;
      44            4 :         if(isHII)
      45            3 :                 n += HIIdensitynumber;
      46            4 :         if(isH2)
      47            3 :                 n += 2*H2densitynumber;
      48              : 
      49              :         // check if all densities are deactivated and raise warning if so
      50            4 :         if((isHI || isHII || isH2) == false){
      51            2 :                 KISS_LOG_WARNING
      52            1 :                         << "\nCalled getNucleonDensity on fully deactivated ConstantDensity "
      53            1 :                         << "gas density model. In this case the density is allways set to 0. \n";
      54              :         }
      55            4 :         return n;
      56              : }
      57              : 
      58            4 : double ConstantDensity::getHIDensity(const Vector3d &position) const {
      59            4 :         return HIdensitynumber;
      60              : }
      61              : 
      62            4 : double ConstantDensity::getHIIDensity(const Vector3d &position) const{
      63            4 :         return HIIdensitynumber;
      64              : }
      65              : 
      66            4 : double ConstantDensity::getH2Density(const Vector3d &position) const{
      67            4 :         return H2densitynumber;
      68              : }
      69              : 
      70            2 : bool ConstantDensity::getIsForHI() {
      71            2 :         return isHI;
      72              : }
      73              : 
      74            2 : bool ConstantDensity::getIsForHII() {
      75            2 :         return isHII;
      76              : }
      77              : 
      78            2 : bool ConstantDensity::getIsForH2() {
      79            2 :         return isH2;
      80              : }
      81              : 
      82            5 : void ConstantDensity::setHI(bool activate, double densitynumber) {
      83            5 :         isHI = activate;
      84            5 :         HIdensitynumber = densitynumber;
      85            5 : }
      86              : 
      87            1 : void ConstantDensity::setHI(bool activate) {
      88            1 :         setHI(activate, HIdensitynumber);
      89            1 : }
      90              : 
      91            1 : void ConstantDensity::setHI(double densitynumber) {
      92            1 :         setHI(isHI, densitynumber);
      93            1 : }
      94              : 
      95            5 : void ConstantDensity::setHII(bool activate, double densitynumber) {
      96            5 :         isHII = activate;
      97            5 :         HIIdensitynumber = densitynumber;
      98            5 : }
      99              : 
     100            1 : void ConstantDensity::setHII(bool activate) {
     101            1 :         setHII(activate, HIIdensitynumber);
     102            1 : }
     103              : 
     104            1 : void ConstantDensity::setHII(double densitynumber) {
     105            1 :         setHII(isHII, densitynumber);
     106            1 : }
     107              : 
     108            5 : void ConstantDensity::setH2(bool activate, double densitynumber) {
     109            5 :         isH2 = activate;
     110            5 :         H2densitynumber = densitynumber;
     111            5 : }
     112              : 
     113            1 : void ConstantDensity::setH2(bool activate) {
     114            1 :         setH2(activate, H2densitynumber);
     115            1 : }
     116              : 
     117            1 : void ConstantDensity::setH2(double densitynumber) {
     118            1 :         setH2(isH2, densitynumber);
     119            1 : }
     120              : 
     121            0 : std::string ConstantDensity::getDescription() {
     122            0 :         std::stringstream s;
     123            0 :         s << "ConstantDensity:\n";
     124            0 :         s<< "HI component is ";
     125            0 :         if(!isHI)
     126            0 :                 s<< "not ";
     127            0 :         s<< "active and has a density of " << HIdensitynumber/ccm << " cm^-3" << "\nHII component is ";
     128            0 :         if(!isHII)
     129            0 :                 s<< "not ";
     130            0 :         s<<"active and has a density of " << HIIdensitynumber/ccm<<" cm^-3" <<  "\nH2 component is ";
     131            0 :         if(!isH2)
     132            0 :                 s<<"not ";
     133            0 :         s<<"active and has a density of " << H2densitynumber/ccm << " cm^-3";
     134            0 :         return s.str();
     135            0 : }
     136              : 
     137              : }  // namespace crpropa