halo.h

/*
 * halo.h
 *
 *  Created on: Feb 2, 2012
 *      Author: cgiocoli
 *
 *      This class implement the properties of an NFW-halo
 */
 
#ifndef HALO_H_
#define HALO_H_
#include <math.h>
#include <cosmo.h>
#include <utilities.h>
//TODO: Change to physical length units !!!!
class HALOCalculator{
public:
 
  HALOCalculator (COSMOLOGY *co, double mass, double redshift);
  
  double getMass(){return m;}
  double getZ(){return z;}
  
  double getRvir(int caseunit=0);
  double getR200();
  double getConcentration(int caseunit=0,double alpha0=-0.1);
  double getFormationTime(double f=0.5);
  double getFormationRedshift(double f=0.5);
  
  static double MosterStellarMassFraction(double Mtotal);
  
  void reset(double,double);
  virtual ~HALOCalculator ();
  
protected:
    COSMOLOGY *co;   // cosmological model
    double m;        // halo mass
    double z;        // redshift
    double R200;     // radius enclosing 200 times the critical density
    double Rvir;     // radius enclosing the virial overdensity
    double Omz,Omo,Oml;  // matter density parameter at redshift z and at z=0
    double sigma2M;   // the mass variance S(m)
    double deltac0;  // spherical collapse overdensity
private:
    void Set_Parameters();
};
typedef HALOCalculator *inCosmoHaloHndl;
 
#endif /* HALO_H_ */