Sample configuration files¶
Environment settings¶
1 2 3 4 5 6 7 8 9 10 | [EnvironmentSettings]
#Must specify
home = SimTest/Home
storage = SimTest/Storage
#Optional
name2attr = {"Om":"Om0","Ol":"Ode0","w":"w0","wa":"wa","h":"h","Ob":"Ob0","si":"sigma8","ns":"ns"}
cosmo_id_digits = 3
json_tree_file = .tree.json
|
Job settings¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 | [CAMB]
#Personal settings
email = apetri@phys.columbia.edu
charge_account = TG-AST140041
#Path to executable
path_to_executable = camb
#Name of the job, output
job_name = CAMB
redirect_stdout = camb.out
redirect_stderr = camb.err
#Resources
cores_per_simulation = 1
tasks_per_node = 24
queue = development
wallclock_time = 02:00:00
#Script name
job_script_file = camb.sh
##################################################################################
[NGenIC]
#Personal settings
email = apetri@phys.columbia.edu
charge_account = TG-AST140041
#Path to executable
path_to_executable = NGenIC
#Name of the job, output
job_name = NGenIC
redirect_stdout = ngenic.out
redirect_stderr = ngenic.err
#Resources
cores_per_simulation = 16
queue = development
wallclock_time = 02:00:00
#Script name
job_script_file = ngenic.sh
##################################################################################
[Gadget2]
#Personal settings
email = apetri@phys.columbia.edu
charge_account = TG-AST140041
#Path to executable
path_to_executable = Gadget2
#Name of the job, output
job_name = Gadget2
redirect_stdout = gadget.out
redirect_stderr = gadget.err
#Resources
cores_per_simulation = 256
queue = development
wallclock_time = 02:00:00
#Script name
job_script_file = gadget.sh
##################################################################################
[LensPlanes]
#Personal settings
email = apetri@phys.columbia.edu
charge_account = TG-AST140041
#Path to executable
path_to_executable = lenstools.planes-mpi
#Name of the job, output
job_name = lenstoolsPlanes
redirect_stdout = planes.out
redirect_stderr = planes.err
#Resources
cores_per_simulation = 16
queue = development
wallclock_time = 02:00:00
#Script name
job_script_file = planes.sh
##################################################################################
[RayTracing]
#Personal settings
email = apetri@phys.columbia.edu
charge_account = TG-AST140041
#Path to executable
path_to_executable = lenstools.raytracing-mpi
#Name of the job, output
job_name = lenstoolsRay
redirect_stdout = ray.out
redirect_stderr = ray.err
#Resources
cores_per_simulation = 2
queue = development
wallclock_time = 02:00:00
#Script name
job_script_file = ray.sh
|
Cluster settings¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | [Directives]
directive_prefix = #SBATCH
job_name_switch = -J
stdout_switch = -o
stderr_switch = -e
num_cores_switch = -n
num_nodes_switch = -N
tasks_per_node_switch = None
queue_type_switch = -p
wallclock_time_switch = -t
user_email_switch = --mail-user=
user_email_type = --mail-type=all
[ClusterSpecs]
shell_prefix = #!/bin/bash
execution_preamble = None
charge_account_switch = -A
job_starter = ibrun
cores_per_node = 16
memory_per_node = 32
memory_unit = Gbyte
cores_at_execution_switch = -n
offset_switch = -o
wait_switch = wait
multiple_executables_on_node = True
|
Single redshift raytracing¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 | [PlaneSettings]
directory_name = Planes
kind = potential
override_with_local = False
format = fits
name_format = snap{0}_{1}Plane{2}_normal{3}.{4}
plane_resolution = 128
first_snapshot = 0
last_snapshot = 58
snapshots = None
snapshot_handler = Gadget2SnapshotDE
cut_points = 10.71
thickness = 3.57
length_unit = Mpc
normals = 0,1,2
##################################################################################
[MapSettings]
directory_name = Maps
override_with_local = False
format = fits
map_resolution = 128
map_angle = 3.5
angle_unit = deg
source_redshift = 2.0
#Random seed used to generate multiple map realizations
seed = 0
#Set of lens planes to be used during ray tracing
plane_set = Planes
plane_format = fits
plane_name_format = snap{0}_potentialPlane{1}_normal{2}.{3}
#N-body simulation realizations that need to be mixed
mix_nbody_realizations = 1
mix_cut_points = 0
mix_normals = 0
lens_map_realizations = 4
first_realization = 1
#Which lensing quantities do we need?
convergence = True
convergence_ks = True
shear = True
omega = True
|
Shear catalogs¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 | [PlaneSettings]
directory_name = Planes
kind = potential
override_with_local = False
format = fits
name_format = snap{0}_{1}Plane{2}_normal{3}.{4}
plane_resolution = 128
first_snapshot = 0
last_snapshot = 58
cut_points = 10.71
thickness = 3.57
length_unit = Mpc
normals = 0,1,2
##################################################################################
[CatalogSettings]
#Name of catalog batch
directory_name = Catalog
input_files = galaxy_positions.fits
total_num_galaxies = 1000
catalog_angle_unit = deg
#Use the options generated at the moment of the batch generation (advised)
override_with_local = True
#Format of the simulated catalog files
format = fits
#Random seed used to generate multiple catalog realizations
seed = 0
#Set of lens planes to be used during ray tracing
plane_set = Planes
plane_format = fits
plane_name_format = snap{0}_potentialPlane{1}_normal{2}.{3}
#N-body simulation realizations that need to be mixed
mix_nbody_realizations = 1
mix_cut_points = 0
mix_normals = 0
lens_catalog_realizations = 1
first_realization = 1
|
CMB lensing reconstructions¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | [CMBReconstruction]
input_set = kappaCMB
input_filename = WLconv*.fits
estimator = TT
unlensed_ps_filename = camb_scalCls.dat
lensed_ps_filename = camb_lensedCls.dat
ps_type = camb_dimensionless
lmax = 3500
output_type = kappa
wiener = False
noise_level_uK_arcmin = 6.0
beam_fwhm_arcmin = 1.4
save_intermediate = False
output_set = kappaCMBRec
output_fname = recTT_{0}
|
Light cone projection¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | [PlaneLightConeSettings]
directory_name = Planes
snapshot_handler = FastPMSnapshotStretchZ
override_with_local = False
zmax = 5.0
num_lenses = 40
plane_resolution = 256
normal = 2
do_lensing = True
integration_type = full
fov_deg = 1.3
fov_resolution = 128
name_format = {1}_z{2:.2f}_normal{3}.{4}
format = fits
|
Line of sight integration¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | [MapSettings]
directory_name = MapsLOS
override_with_local = False
format = fits
map_resolution = 2048
map_angle = 3.41
angle_unit = deg
source_redshift = 2.0
#Random seed used to generate multiple map realizations
seed = 0
#Set of lens planes to be used during ray tracing
plane_set = Planes
#N-body simulation realizations that need to be mixed
mix_nbody_realizations = 1
mix_cut_points = 0,1,2
mix_normals = 0,1,2
lens_map_realizations = 1024
#LOS integration type
integration_type = born
|
Integration types¶
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | """
Command line scripts that come with lenstools
"""
##############################################
#######Line of sight integration##############
##############################################
integration_types = {
"born" : "Born approximation for the convergence (linear order in the lensing potential)",
"born-rt" : "Semi-Born approximation: integrate the density on the real ray trajectory",
"postBorn2" : "Convergence at second order in the lensing potential (lens-lens + geodesic perturbation)",
"postBorn2-ll" : "Convergence due to lens-lens coupling",
"postBorn2-gp" : "Convergence due to geodesic perturbation",
"postBorn1+2" : "Convergence at Born + second post Born order in the lensing potential",
"postBorn1+2-ll" : "Convergence at Born + second post Born order (lens-lens only)",
"postBorn1+2-gp" : "Convergence at Born + second post Born order (geodesic perturbation only)",
"omega2" : "Rotation at second order in the lensing potential"
}
|