ccl_redshifts.h

Go to the documentation of this file.

#ifndef __CCL_LSST_SPECS_INCLUDED__
#define __CCL_LSST_SPECS_INCLUDED__

#include <gsl/gsl_integration.h>
#include <gsl/gsl_spline.h>

#define Z_MIN_SOURCES 0.
#define Z_MAX_SOURCES 5.0

CCL_BEGIN_DECLS
typedef struct {
        double (* your_pz_func)(double, double, void *, int*); /*< Function returns the likelihood of measuring a z_ph
        * (first double) given a z_spec (second double), with a pointer to additonal arguments and a status flag.*/
        void *  your_pz_params; /*< Additional parameters to be passed into your_pz_func */
} pz_info;

typedef struct {
        double (* your_dN_func)(double, void *, int*); /*< Function returns the differential number density of galaxies wrt redshifts, 
        * with a pointer to additonal arguments and a status flag.*/
        void *  your_dN_params; /*< Additional parameters to be passed into your_dN_func */
} dNdz_info;

 typedef struct{
        double alpha;
        double beta;
        double z0;
    } smail_params;

void ccl_dNdz_tomog(double z, double bin_zmin, double bin_zmax, pz_info * photo_info,  dNdz_info * dN_info, double *tomoout, int *status);

dNdz_info* ccl_create_dNdz_info(void * params, double(*dNdz_func)(double,void*,int*));

 dNdz_info* ccl_create_Smail_dNdz_info(double alpha, double beta, double z0);


void ccl_free_dNdz_info(dNdz_info * dN_info);

pz_info* ccl_create_photoz_info(void * params, double(*pz_func)(double, double,void*,int*));

pz_info* ccl_create_gaussian_photoz_info(double sigma_z0);


void ccl_free_photoz_info(pz_info *my_photoz_info);

CCL_END_DECLS

#endif