C++ (Cpp) getMvir Examples

Example #1

File: haloes.cpp Project: surhudm/aum

/// Radii are in physical units here, so be careful
void cosmology::modelNFWhalo(double m200,double z,double &Mvir, double &Rvir, double &cvir, double &R200, double &c200)
{
    Mvir=getMvir(m200,z);
    Rvir=getRvirfromMvir(Mvir, z)/(1.+z);

    double Delta=200.;
    R200=getRDelfromMDel(m200, z, Delta)/(1.+z);

    cvir=conc(Mvir,z);

    c200=getc200(cvir,z);

}

Example #2

File: haloes.cpp Project: anupreeta27/SIMCT

/// Radii are in comoving units here
void cosmology::modelNFWhalo_com(double m200,double z,double &Mvir, double &Rvir, double &cvir)
{
    Mvir=getMvir(m200,z);

    Rvir=0.169*pow(Mvir/1.e12,1./3.);
    Rvir*=pow(Delta_crit(z)/178.0,-1.0/3.0);
    Rvir*=pow(Eofz(z),-2./3.);

    cvir=conc(Mvir,z);

    Rvir=Rvir*(1.+z);

}

Example #3

File: haloes.cpp Project: anupreeta27/SIMCT

/// Radii are in physical units here, so be careful
void cosmology::modelNFWhalo(double m200,double z,double &Mvir, double &Rvir, double &cvir, double &R200, double &c200)
{
    Mvir=getMvir(m200,z);

    Rvir=0.169*pow(Mvir/1.e12,1./3.);
    Rvir*=pow(Delta_crit(z)/178.0,-1.0/3.0);
    Rvir*=pow(Eofz(z),-2./3.);

    double Delta=200.;
    R200=pow( 1.0e12/(4./3.*Delta*3e4*0.3/(8*gee)),1./3.);
    R200*=pow(m200/1.e12,1./3.);
    R200*=pow(Omega(z)/0.3,-1.0/3.0);
    R200*=pow(Eofz(z),-2./3.);

    cvir=conc(Mvir,z);

    c200=getc200(cvir,z);

}

Example #4

File: haloes.cpp Project: surhudm/aum

/// Radii are in comoving units here
void cosmology::modelNFWhalo_com(double m200,double z,double &Mvir, double &Rvir, double &cvir, double &R200, double &c200)
{
    Mvir=getMvir(m200,z);
    Rvir=getRvirfromMvir(Mvir, z);
    
    double Delta=200.;
    R200=getRDelfromMDel(m200, z, Delta);

    /*
    R200=pow( 1.0e12/(4./3.*Delta*3e4*0.3/(8*gee)),1./3.);
    R200*=pow(m200/1.e12,1./3.);
    R200*=pow(Omega(z)/0.3,-1.0/3.0);
    R200*=pow(Eofz(z),-2./3.);
    */

    cvir=conc(Mvir,z);
    c200=getc200(cvir,z);

}

Example #5

File: cosmology.cpp Project: surhud/pevolve

/// Radii are in comoving units here
void cosmology::modelNFWhalo_com_ext(double mDel,double z,double &Mvir,
double &Rvir, double &cvir, double &RDel, double &cDel,double Del)
{
    //std::cout<<Delta_crit(z)<<std::endl;exit(0);
    Mvir=getMvir(mDel,z,Del);

    Rvir=0.169*pow(Mvir/1.e12,1./3.);
    Rvir*=pow(Delta_crit(z)/178.0,-1.0/3.0);
    Rvir*=pow(Eofz(z),-2./3.);

    double Delta=Del;
    RDel=pow( 1.0e12/(4./3.*Delta*3e4*0.3/(8*gee)),1./3.);
    RDel*=pow(mDel/1.e12,1./3.);
    RDel*=pow(Omega(z)/0.3,-1.0/3.0);
    RDel*=pow(Eofz(z),-2./3.);

    cvir=conc(Mvir,z);

    cDel=getcDel(cvir,z,Delta);

    Rvir=Rvir*(1.+z);
    RDel=RDel*(1.+z);

}