fsky = Config.getfloat(expName, 'fsky')
try:
v3mode = Config.getint(expName, 'V3mode')
except:
v3mode = -1
fparams = {}
for (key, val) in Config.items('params'):
if ',' in val:
param, step = val.split(',')
fparams[key] = float(param)
else:
fparams[key] = float(val)
cc = ClusterCosmology(fparams, constDict, clTTFixFile=clttfile)
HMF = Halo_MF(cc, mexp_edges, z_edges)
SZCluster = SZ_Cluster_Model(cc,
clusterDict,
rms_noises=noise,
fwhms=beam,
freqs=freq,
lknee=lknee,
alpha=alpha,
v3mode=v3mode,
fsky=fsky)
HMF.updateSigN(SZCluster)
HMF.updatePfunc(SZCluster)
np.save(bigDataDir + "Sel_func" + saveId, HMF.Pfunc)
cc = ClusterCosmology(fparams,constDict,clTTFixFile=clttfile)
HMF = Halo_MF(cc,mgrid,zgrid)
qs = listFromConfig(Config,'general','qbins')
qspacing = Config.get('general','qbins_spacing')
if qspacing=="log":
qbin_edges = np.logspace(np.log10(qs[0]),np.log10(qs[1]),int(qs[2])+1)
elif qspacing=="linear":
qbin_edges = np.linspace(qs[0],qs[1],int(qs[2])+1)
else:
raise ValueError
SZProp = SZ_Cluster_Model(cc,clusterDict,rms_noises = noise,fwhms=beam,freqs=freq,lknee=lknee,alpha=alpha)
Mwl = 10**HMF.Mexp
z_arr = HMF.zarr
M_arr = np.outer(HMF.M,np.ones([len(z_arr)]))
if HMF.sigN is None: HMF.updateSigN(SZProp)
sigN = HMF.sigN
q_arr = old_div((qbin_edges[1:]+qbin_edges[:-1]),2.)
blah = SZProp.P_of_qn_corr(SZProp.lnY,M_arr,z_arr,sigN,q_arr,Mwl)#,lndM)
#dN_dmqz_corr = HMF.N_of_mqz_SZ_corr(lndM,qbin_edges,SZProp)
#dN_dmqz = HMF.N_of_mqz_SZ(lndM,qbin_edges,SZProp)
