print((len(ras))) #sys.exit() p143Loc = '/astro/astronfs01/workarea/msyriac/PlanckClusters/WPR2_CMB_muK.fits' hpPlanck = hp.read_map(p143Loc) outDir = os.environ['WWW'] + "plots/kappatest/" # === COSMOLOGY === cosmologyName = 'LACosmology' # from ini file iniFile = "../SZ_filter/input/params.ini" Config = SafeConfigParser() Config.optionxform = str Config.read(iniFile) lmax = 8000 cosmoDict = dictFromSection(Config, cosmologyName) constDict = dictFromSection(Config, 'constants') cc = ClusterCosmology(cosmoDict, constDict, lmax, pickling=True) TCMB = 2.7255e6 # === QUADRATIC ESTIMATOR INITIALIZATION === from orphics.analysis import flatMaps as fmaps from alhazen.quadraticEstimator import Estimator # === RUN N SIMS === kappaStack = 0. stampStack = 0. maxN = None #10000 #None pixScale = 1.0
from szar.counts import ClusterCosmology import numpy as np import sys, os from configparser import SafeConfigParser import pickle as pickle from orphics.tools.io import dictFromSection, listFromConfig, Plotter import szar._fast as fast M = 2.e14 z = 0.1 iniFile = "input/pipeline.ini" Config = SafeConfigParser() Config.optionxform = str Config.read(iniFile) constDict = dictFromSection(Config, 'constants') fparams = {} # the for (key, val) in Config.items('params'): if ',' in val: param, step = val.split(',') fparams[key] = float(param) else: fparams[key] = float(val) clttfile = Config.get('general', 'clttfile') #fparams['H0']=100. cc = ClusterCosmology(fparams, constDict, clTTFixFile=clttfile) print((cc.Mdel_to_cdel(M, z, 500.)))
#Mexp = 14.4 #np.log10(2.e14) Mexp = np.log10(5.e14) z = 1.5 c = 1.18 clusterParams = 'LACluster' # from ini file cosmologyName = 'LACosmology' # from ini file iniFile = "../SZ_filter/input/params.ini" Config = SafeConfigParser() Config.optionxform = str Config.read(iniFile) lmax = 8000 cosmoDict = dictFromSection(Config, cosmologyName) constDict = dictFromSection(Config, 'constants') clusterDict = dictFromSection(Config, clusterParams) cc = ClusterCosmology(cosmoDict, constDict, lmax, pickling=True) theory = cc.theory beamrange = np.arange(0.5, 5.0, 0.25) # Make a CMB Noise Curve cambRoot = "data/ell28k_highacc" gradCut = 2000 halo = True #beamX = 1.0 #beamY = 1.0 noiseTX = 10.0 noisePX = 14.14
from szar.counts import ClusterCosmology from szar.szproperties import SZ_Cluster_Model from szar.foregrounds import fgNoises, f_nu from szar.ilc import ILC_simple import sys, os from configparser import SafeConfigParser import pickle as pickle from orphics.tools.io import dictFromSection, listFromConfig, Plotter from orphics.tools.cmb import noise_func iniFile = "input/pipeline.ini" Config = SafeConfigParser() Config.optionxform = str Config.read(iniFile) constDict = dictFromSection(Config, 'constants') clusterDict = dictFromSection(Config, 'cluster_params') fparams = {} # the for (key, val) in Config.items('params'): if ',' in val: param, step = val.split(',') fparams[key] = float(param) else: fparams[key] = float(val) clttfile = Config.get('general', 'clttfile') cc = ClusterCosmology(fparams, constDict, lmax=8000, pickling=True) #clTTFixFile=clttfile) fgs = fgNoises(cc.c, ksz_battaglia_test_csv="data/ksz_template_battaglia.csv",