MKONE (Mock Light Cone Utilities)

Generate mock light cones with filamentary strucure or Gaussian blobs added

Note the structure does not intend to be realistic and the mock catalogs can be useful for teaching or testing correlation functions but not to do real science (except for routine uniform_sky())

Main Routines

mcone_filam Generate a mock light cone of random 3D points emulating filamentary structure
mcone_gaussblobs Generate a mock light cone of random 3D points with added gaussian blobs or "clusters"
filament_box Generate random 3D points emulating filamentary structure over a box
fill_cone Fill the volume of an observation light cone with small cubes of 3D points
uniform_sky Generate uniform random points (ra,dec) in a given area of sky

Example of use

# Create a 60x60 deg light cone with z=[0.01,0.15], filled with filaments
import mkone as mk
ralim  = [0,60]
declim = [0,60]
zlim   = [0.01,0.15]
kone = mk.mcone_filam(ralim,declim,zlim,npts=80000,nvoids=2000,nstep=100,b=150.)
# Send to Topcat (start Topcat first)
mk.send(kone,'kone')

# Create a 40x20 deg light cone with z=[0.01,0.15], filled with gaussian "clusters"
import mkone as mk
ralim   = [10,50]
declim  = [20,40]
zlim    = [0.01,0.15]
cradlim = [0.5,20]
kone = mk.mcone_gaussblobs(ralim,declim,zlim,n=20000,ufrac=0.7,ncen=500,
                           cradlim=cradlim,oformat='array')
# Send to Topcat (start Topcat first)
mk.send(kone,'kone',cols=['ra','dec','z','comd','px','py','pz'])

Dependencies

astropy : to provide table output and cosmological functions
matplotlib,mplot3d : to provide basic graphics
scipy.spatial : to provide cKDTree for fast searching nearest voids
sampc : to exchange data over SAMP with some VO apps (e.g. Topcat)

To Do

Add seed random as input parameter where needed

Credits