def multiscreen_halo( specfile, params, amin=AMIN, alpha=ALPHA, **kwargs ):
x1, x2, logNH1, logNH2, amax, p = params
s1 = MH.simulate_screen( specfile, xg=x1, NH=np.power(10.0,logNH1), \
a0=AMIN, a1=amax, p=p, alpha=alpha, **kwargs )
s2 = MH.simulate_screen( specfile, xg=x2, NH=np.power(10.0,logNH2), \
a0=AMIN, a1=amax, p=p, alpha=alpha, **kwargs )
return sum_interp( s1, s2 )
def uniscreen(specfile, params, alpha=ALPHA, **kwargs):
logNHu, logNHs, a_u, a_s, p_u, p_s, x_s = params
nhu = np.power(10.0, logNHu)
UU = MH.simulate_uniform(specfile, NH=nhu, a0=AMIN, a1=a_u, p=p_u, alpha=alpha, **kwargs)
nhs = np.power(10.0, logNHs)
SS = MH.simulate_screen(specfile, xg=x_s, NH=nhs, a0=AMIN, a1=a_s, p=p_s, alpha=alpha, **kwargs)
return sum_interp(UU, SS)
def uniscreen( specfile, params, alpha=ALPHA, **kwargs ):
logNHu, logNHs, a_u, a_s, p_u, p_s, x_s = params
nhu = np.power( 10.0, logNHu )
UU = MH.simulate_uniform( specfile, NH=nhu, \
a0=AMIN, a1=a_u, p=p_u, alpha=alpha, **kwargs )
nhs = np.power( 10.0, logNHs )
SS = MH.simulate_screen( specfile, xg=x_s, NH=nhs, \
a0=AMIN, a1=a_s, p=p_s, alpha=alpha, **kwargs )
return sum_interp( UU, SS )
def screen_halo(filename, params, alpha=ALPHA, **kwargs):
xg, logNH, amax, p = params
return MH.simulate_screen(filename, xg=xg, NH=np.power(10.0, logNH), a0=AMIN, a1=amax, p=p, alpha=alpha, **kwargs)
def screen_halo( filename, params, alpha=ALPHA, **kwargs ):
xg, logNH, amax, p = params
return MH.simulate_screen( filename, \
xg=xg, NH=np.power(10.0,logNH), a0=AMIN, a1=amax, p=p, alpha=alpha, **kwargs)
