lesgourg / class_public

Public repository of the Cosmic Linear Anisotropy Solving System (master for the most recent version of the standard code; GW_CLASS to include Cosmic Gravitational Wave Background anisotropies; classnet branch for acceleration with neutral networks; ExoCLASS branch for exotic energy injection; class_matter branch for FFTlog)
230 stars 285 forks source link

Selection_bias as a function #359

Open FigoHerranz opened 4 years ago

FigoHerranz commented 4 years ago

Hi,

I am currently trying to implement a model where the bias is a function of redshift with some parameters, let say for simplicity b(z)=b_0*f(z) where I only need to have f(z) valuated in theselection_mean redshifts, no need to have it as a continuous function . I know the parameter of CLASS selection_bias can be used for different bias in the different redshift but I would like to have b_0 as a free parameter read from the .ini file and then for future applications, for example with MontePython as a parameter. Any idea on how to implement this b(z)=b_0*f(z) in the CLASS selection_bias ?

Thanks everyone in this forum, I hope we can find an answer!

ShazAlvi commented 4 years ago

Radical idea. How about you use the selection_bias as the free parameter and edit the relevant pieces of the code in transfer.c to accommodate the redshift dependent function. You can code to use the bin_means read from the ini file. But, I warn, this is a radical idea completely off the top of my head and I assume no liability for it. Try at your own risk and don't try to find me :P