ccl_cls.h

Go to the documentation of this file.

#ifndef __CCL_CLS_H_INCLUDED__
#define __CCL_CLS_H_INCLUDED__

typedef enum ccl_tracer_t
{
  ccl_number_counts_tracer = 1,
  ccl_weak_lensing_tracer  = 2,
  ccl_cmb_lensing_tracer   = 3,
} ccl_tracer_t;

typedef enum ccl_tracer_func_t
{
  ccl_trf_nz = 201, //Redshift distribution
  ccl_trf_bz = 202, //Clustering bias
  ccl_trf_sz = 203, //Magnification bias
  ccl_trf_rf = 204, //Aligned fraction
  ccl_trf_ba = 205, //Alignment bias
  ccl_trf_wL = 206, //Weak lensing window function
  ccl_trf_wM = 207, //Magnification window function
} ccl_tracer_func_t;

CCL_BEGIN_DECLS

typedef struct {
  int tracer_type; //Type (see above)
  double prefac_lensing; //3*O_M*H_0^2/2
  double chimax; //Limits in chi where we care about this tracer
  double chimin;
  double zmin; //Limits in chi where we care about this tracer
  double zmax;
  double chi_source; //Comoving distance to the source (for CMB lensing)
  int has_rsd;
  int has_magnification;
  int has_intrinsic_alignment;
  SplPar *spl_nz; //Spline for normalized N(z)
  SplPar *spl_bz; //Spline for linear bias
  SplPar *spl_sz; //Spline for magnification bias
  SplPar *spl_rf; //Spline for red fraction
  SplPar *spl_ba; //Spline for alignment bias
  SplPar *spl_wL; //Spline for lensing kernel
  SplPar *spl_wM; //Spline for magnification
} CCL_ClTracer;


CCL_ClTracer *ccl_cl_tracer(ccl_cosmology *cosmo,int tracer_type,
                int has_rsd,int has_magnification,int has_intrinsic_alignment,
                int nz_n,double *z_n,double *n,
                int nz_b,double *z_b,double *b,
                int nz_s,double *z_s,double *s,
                int nz_ba,double *z_ba,double *ba,
                int nz_rf,double *z_rf,double *rf,
                double z_source,int * status);

CCL_ClTracer *ccl_cl_tracer_number_counts(ccl_cosmology *cosmo,
                          int has_rsd,int has_magnification,
                          int nz_n,double *z_n,double *n,
                          int nz_b,double *z_b,double *b,
                          int nz_s,double *z_s,double *s, int * status);


CCL_ClTracer *ccl_cl_tracer_number_counts_simple(ccl_cosmology *cosmo,
                             int nz_n,double *z_n,double *n,
                             int nz_b,double *z_b,double *b, int * status);

CCL_ClTracer *ccl_cl_tracer_lensing(ccl_cosmology *cosmo,
                    int has_alignment,
                    int nz_n,double *z_n,double *n,
                    int nz_ba,double *z_ba,double *ba,
                    int nz_rf,double *z_rf,double *rf, int * status);

CCL_ClTracer *ccl_cl_tracer_lensing_simple(ccl_cosmology *cosmo,
                           int nz_n,double *z_n,double *n, int * status);


CCL_ClTracer *ccl_cl_tracer_cmblens(ccl_cosmology *cosmo,double z_source,int *status);

void ccl_cl_tracer_free(CCL_ClTracer *clt);

double ccl_get_tracer_fa(ccl_cosmology *cosmo,CCL_ClTracer *clt,double a,int func_code,int *status);

int ccl_get_tracer_fas(ccl_cosmology *cosmo,CCL_ClTracer *clt,int na,double *a,double *fa,
               int func_code,int *status);

//Workspace for C_ell computations
typedef struct {
  int lmax; //*Maximum multipole
  int l_limber; //*All power spectra for l>l_limber will be computed using Limber's approximation
  double l_logstep; //*Logarithmic step factor used at low l
  int l_linstep; //*Linear step used at high l
  int n_ls; //Number of multipoles that result from the previous combination of parameters
  int *l_arr; //*Array of multipole values resulting from the previous parameters
} CCL_ClWorkspace;

//CCL_ClWorkspace constructor
CCL_ClWorkspace *ccl_cl_workspace_new(int lmax,int l_limber,
                      double l_logstep,int l_linstep,
                      int *status);
//CCL_ClWorkspace simplified constructor
CCL_ClWorkspace *ccl_cl_workspace_new_limber(int lmax, double l_logstep,int l_linstep,int *status);
//CCL_ClWorkspace destructor
void ccl_cl_workspace_free(CCL_ClWorkspace *w);

void ccl_angular_cls(ccl_cosmology *cosmo,CCL_ClWorkspace *w,
             CCL_ClTracer *clt1,CCL_ClTracer *clt2,
             int nl_out,int *l,double *cl,int *status);

CCL_END_DECLS


#endif