Exemplo n.º 1
0
"""
run this code to create the directories
"""

import os
import sys
sys.path.append(os.path.abspath('../'))
from shear_KL_source import params
params.load('base_params.dat')

for param_file in ('base_params.dat',
                   'params_f0.5.dat',
                   'params_sig0.1.dat',
                   'params_sig0.1_f0.5.dat',
                   'params_sig0.1_f0.35.dat'):
    params.load(param_file)
    for p in ['scratch_dir', 'shear_in_dir', 'shear_recons_dir', 'Map_dir',
              'kappa_dir', 'mask_outdir', 'condorlog']:
        F = params[p]
        if not os.path.exists(F):
            print "mkdir %s" % F
            os.system('mkdir %s' % F)
Exemplo n.º 2
0
    if logplot:
        pylab.semilogx(xvals,SN,label='Signal+Noise')
        pylab.semilogx(xvals,S,label='Signal')
        pylab.semilogx(xvals,N,'--k',label='Noise')
        ax.xaxis.set_major_formatter(FuncFormatter(lambda x,*args: '%i'%x))
    else:
        pylab.plot(xvals,SN,label='Signal+Noise')
        pylab.plot(xvals,S,label='Signal')
        pylab.plot(xvals,N,'--k',label='Noise')
        pylab.legend(loc=4)
    
    pylab.xlim(1,4096)
    pylab.grid(True,c='#AAAAAA')
    pylab.xlabel('mode number')
    pylab.ylabel('cumulative value')

if __name__ == '__main__':
    import os
    import sys
    sys.path.append(os.path.abspath('../'))
    from shear_KL_source import params
    params.load('../run/base_params.dat')

    from shear_KL_source.DES_tile.tools import get_basis_filename
    
    plot_evals( get_basis_filename() )
    pylab.savefig('figs/fig03_eigenvalues.eps')
    pylab.savefig('figs/fig03_eigenvalues.pdf')
    pylab.show()
Exemplo n.º 3
0
def plot_reconstruction_2():
    nmodes=900
    alpha=0.15
    RAmin=11.5
    DECmin=36
    xlim=(15,45)
    ylim=(20,50)
    
    RAlim = (RAmin,RAmin+1)
    DEClim = (DECmin,DECmin+1)

    RAside = (RAlim[1]-RAlim[0])
    DECside = (DEClim[1]-DEClim[0])
    
    params.load('../run/base_params.dat')
    params.load('../run/params_sig0.1.dat')
    gamma_out_masked_1, noise = get_reconstruction(nmodes,
                                                 alpha,
                                                 True,
                                                 RAmin,
                                                 DECmin)
    gamma_out_unmasked_1, noise = get_reconstruction(nmodes,
                                                   alpha,
                                                   False,
                                                   RAmin,
                                                   DECmin)
    mask_1 = get_mask(False,RAlim,DEClim)
    
    params.load('../run/params_sig0.1_f0.35.dat')
    gamma_out_masked_2, noise = get_reconstruction(nmodes,
                                                 alpha,
                                                 True,
                                                 RAmin,
                                                 DECmin)
    mask_2 = get_mask(False,RAlim,DEClim)
    
    params.load('../run/params_sig0.1_f0.5.dat')
    gamma_out_masked_3, noise = get_reconstruction(nmodes,
                                                 alpha,
                                                 True,
                                                 RAmin,
                                                 DECmin)
    mask_3 = get_mask(False,RAlim,DEClim)

    offset = (RAmin,DECmin)
    dtheta = params.dtheta/60.
    
    xlim = (offset[0] + xlim[0]/60.,
            offset[0] + xlim[1]/60.)
    ylim = (offset[1] + ylim[0]/60.,
            offset[1] + ylim[1]/60.)

    extent=(offset[0],offset[0]+RAside,
            offset[1],offset[1]+DECside)
    
    #--------------------------------------------------
    pylab.figure(figsize=(8,8))
    ax = pylab.subplot(221)
    pylab.imshow(mask_1.T*0,
                 origin='lower',
                 extent=extent,
                 interpolation='nearest',
                 cmap=pylab.cm.binary)
    whiskerplot(gamma_out_unmasked_1,dtheta,dtheta,offset=offset)
    pylab.xlim(xlim)
    pylab.ylim(ylim)
    #pylab.title(r'$\mathdefault{unmasked\ KL}$')
    #pylab.xlabel('')
    pylab.text(0.97,0.97,'no mask',
               transform = ax.transAxes,
               va = 'top',
               ha = 'right',
               bbox=dict(facecolor='w',edgecolor='w') )
    ax.xaxis.set_major_locator(MultipleLocator(0.2))
    ax.yaxis.set_major_locator(MultipleLocator(0.2))

    #--------------------------------------------------
    ax = pylab.subplot(222)
    pylab.imshow(mask_1.T,
                 origin='lower',
                 extent=extent,
                 interpolation='nearest',
                 cmap=GreyWhite)
    pylab.clim(0,1)
    whiskerplot(gamma_out_masked_1,dtheta,dtheta,offset=offset)
    pylab.xlim(xlim)
    pylab.ylim(ylim)
    #pylab.title(r'$\mathdefault{masked\ KL}$')
    #pylab.xlabel('')
    pylab.ylabel('')
    pylab.text(0.97,0.97,'20% mask',
               transform = ax.transAxes,
               va = 'top',
               ha = 'right',
               bbox=dict(facecolor='w',edgecolor='w') )
    ax.xaxis.set_major_locator(MultipleLocator(0.2))
    ax.yaxis.set_major_locator(MultipleLocator(0.2))

    #--------------------------------------------------
    ax = pylab.subplot(223)
    pylab.imshow(mask_2.T,
                 origin='lower',
                 extent=extent,
                 interpolation='nearest',
                 cmap=GreyWhite)
    pylab.clim(0,1)
    whiskerplot(gamma_out_masked_2,dtheta,dtheta,offset=offset)
    pylab.xlim(xlim)
    pylab.ylim(ylim)
    #pylab.title(r'$\mathdefault{masked\ KL}$')
    pylab.text(0.97,0.97,'35% mask',
               transform = ax.transAxes,
               va = 'top',
               ha = 'right',
               bbox=dict(facecolor='w',edgecolor='w') )
    ax.xaxis.set_major_locator(MultipleLocator(0.2))
    ax.yaxis.set_major_locator(MultipleLocator(0.2))

    #--------------------------------------------------
    ax = pylab.subplot(224)
    pylab.imshow(mask_3.T,
                 origin='lower',
                 extent=extent,
                 interpolation='nearest',
                 cmap=GreyWhite)
    pylab.clim(0,1)
    whiskerplot(gamma_out_masked_3,dtheta,dtheta,offset=offset)
    pylab.xlim(xlim)
    pylab.ylim(ylim)
    #pylab.title(r'$\mathdefault{masked\ KL}$')
    pylab.ylabel('')
    pylab.text(0.97,0.97,'50% mask',
               transform = ax.transAxes,
               va = 'top',
               ha = 'right',
               bbox=dict(facecolor='w',edgecolor='w') )
    ax.xaxis.set_major_locator(MultipleLocator(0.2))
    ax.yaxis.set_major_locator(MultipleLocator(0.2))
    
    #--------------------------------------------------

    pylab.figtext(0.5,0.97,r'$\mathdefault{KL\ Reconstructions}$',
                  ha='center',va='top',
                  fontsize=16)
    pylab.figtext(0.5,0.94,r'$\mathdefault{n=900\ modes,\ \ n_{gal}=100/arcmin^2}$',
                  ha='center',va='top',
                  fontsize=14)
Exemplo n.º 4
0
"""
run this code to compute the aperture mass map from the input shear
"""

import os
import sys
sys.path.append(os.path.abspath('../'))
from shear_KL_source import params
params.load('base_params.dat')

from shear_KL_source.DES_tile.mass_maps import compute_Map_input

for usemask in ['y','n']:
    for normed in [True,False]:
        if usemask=='n' and normed: continue
        for add_signal in (True, False):
            for add_noise in (True, False):
                if add_signal == False and add_noise == False:
                    continue
                compute_Map_input(add_signal = add_signal,
                                  add_noise = add_noise,
                                  usemask = usemask,
                                  normed = normed)
Exemplo n.º 5
0
"""
run this code to create the cfg file for the condor runs
"""

import os
import sys
sys.path.append(os.path.abspath('../'))
from shear_KL_source import params
params.load('base_params.dat')

from shear_KL_source.DES_tile.setup_utils import create_cfg_file

create_cfg_file( cfg_file = 'DES_tile.cfg',
                 alphas = (0.15,),
                 NMODES = (900,),
                 use_noise = (True,),
                 use_mask = (True,False),
                 noise_only = False,
                 compute_shear_noise = False,
                 compute_Map = True,
                 compute_Map_noise = True,
                 )

create_cfg_file( cfg_file = 'DES_tile_noiseonly.cfg',
                 alphas = (0.15,),
                 NMODES = (900,),
                 use_noise = (True,),
                 use_mask = (True,False),
                 noise_only = True,
                 compute_shear_noise = False,
                 compute_Map = True,
Exemplo n.º 6
0
"""
run this to create the mask tiles for the condor runs

This creates masks with different mask fractions
"""

import os
import sys
sys.path.append(os.path.abspath('../'))
from shear_KL_source import params

from shear_KL_source.DES_tile.setup_utils import create_mask_tiles

import numpy
params.load('base_params.dat')

params.load('params_f0.5.dat')
numpy.random.seed(3)
create_mask_tiles( params.mask_outdir,
                   fmask = params.fmask )

params.load('params_f0.35.dat')
numpy.random.seed(3)
create_mask_tiles( params.mask_outdir,
                   fmask = params.fmask )