guilhem_lavaux/cosmotool
0
0
Fork 0
Code Issues Pull requests Projects Releases 1 Packages Wiki Activity Actions

Added help strings and ksz support

Browse source
This commit is contained in:
Guilhem Lavaux 2014-06-18 16:13:16 +02:00
parent ff01447ab6
commit 247d8512e4
7 changed files with 333 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

9
python_sample/build_2lpt_skymap.py
View file

@ -18,6 +18,7 @@ parser.add_argument('--maxval', type=float, default=60)
parser.add_argument('--depth_min', type=float, default=10)
parser.add_argument('--depth_max', type=float, default=60)
parser.add_argument('--iid', type=int, default=0)
parser.add_argument('--proj_cat', type=bool, default=False)
args = parser.parse_args()
INDATA="/nethome/lavaux/Copy/PlusSimulation/BORG/Input_Data/2m++.npy"
@ -61,8 +62,12 @@ for i in xrange(args.start,args.end,args.step):
plt.clf()
d = np.load(args.base_cic % i)
proj = build_sky_proj(1+d, dmin=args.depth_min,dmax=args.depth_max,iid=args.iid)
proj /= (args.depth_max-args.depth_min)
hp.mollview(proj, fig=1, coord='CG', cmap=plt.cm.coolwarm, title='Sample %d' % i, min=args.minval, max=args.maxval)
hp.projscatter(b[idx], l[idx], lw=0, color='g', s=2.0, alpha=0.7)
hp.write_map("skymaps/proj_map_%d.fits" % i, proj)
hp.mollview(proj, fig=1, coord='CG', cmap=plt.cm.copper, title='Sample %d' % i, min=args.minval, max=args.maxval)
if args.proj_cat:
hp.projscatter(b[idx], l[idx], lw=0, color=[0.1,0.8,0.8], s=2.0, alpha=0.7)
ff.savefig(args.base_fig % i)