lenstool package
Submodules
lenstool.constant module
lenstool.dimension module
lenstool.lenstool module
- class lenstool.lenstool.Lenstool(infile=None)[source]
Bases:
objectA lenstool instance initialized with a paramater file, and providing some high-level functions of the C version.
The simplest way to initialize a Lenstool object is to load a standard lenstool parameter file
>>> lenstool = lenstool.Lenstool('test_galaxies.par')
Make sure to have all the file dependencies mentioned in the parameter file available in your current directory.
Note that a Lenstool instance is not thread-safe. The global variables initialized in memory by the underlying C library will be shared between multiple Lenstool instances.
- Parameters:
infile: The path to the parameter file
- property C
Return the C structure
- property CL
Return the CL structure
- property F
Return the F structure
- property G
Return the G structure
- property I
Return the I structure
- property M
Return the M global variable, with the runmode content
- property S
Return the S structure
- property arclet
Return the arclet array global variable
- create_field_wcs(nx, ny)[source]
Return a astropy.wcs.WCS object based on field boundary and reference coordinates in F and M global variables
- create_wcs(nx, ny, xmin, xmax, ymin, ymax, ra, dec)[source]
Return a astropy.wcs.WCS object based on boundary, size, and reference center
- criticnew(zs=None, xmin=None, ymin=None, xmax=None, ymax=None)[source]
Compute critical and caustic lines and return two arrays of structure.biline for the tangential and radial critic and caustic lines. These can be parsed using lenstool.pcl.parse_cline() function.
If zs is None, self.S.zs is used. If one of xmin, xmax, ymin, ymax is None, corresponding field limit is taken from self.F
Modified internal variables: CL.cz[0], CL.nplan, F.xmin, F.xmax, F.ymin, F.ymax, tangent, radial
Return
tangent, radial : arrays of structure.biline, each for critic and caustic lines
- e_pixel_source(np)[source]
Compute and return the image of the arcs and arclets in the source plane
- get_chains(return_ndarray=True)[source]
- Return the content of the bayes.dat. If return_ndarray is True, the
returned numpy.ndarray contains a copy of the internal memory content.
Parameters
- return_ndarray Return a numpy.ndarray if True, otherwise return the
content in the internal ctype square_double_t type.
- property image
Return the image array global variable
- init_cosmoratio(cp0=None)[source]
Initialize the internal cosmological distance look-up-table
Parameters
- cp0: list of tuples
Each list element should contain a tuple (z1, z2, angular_diameter_distance(z1, z2))
- property lens
Return the lens array
- property multi
Return the multi array global variable
- property narclet
Return the narclet global variable
- property nrline
Return the nrline global variable
- property ntline
Return the ntline global variable
- o_chi_lhood0()[source]
Return the chi2 and likelihood corresponding to the model and the constraints previously stored in memory.
Returns
error, chi, lhood0 The error value (0 is good, 1 is bad model), the chi2 and the normalization likelihood term
- property ps
Return the ps structure
- readConstraints()[source]
Read in memory the constraints files (multiple images, arclets, etc.) mentioned in the input parameter file.
- rescaleCube(cube, return_rescaled=True)[source]
Rescale the cube values to their corresponding parameter values in Lenstool.
- Parameters
- cube[nparams] The array containing the parameters for all the potentials. Cube values must
follow uniform distributions in the range [0,1].
return_rescaled If True, cube values are replaced in-place by their rescaled values.
- Returns
Return 1 if the rescaled parameters are meaningfull in the parameter space, 0 otherwise.
- setBayesModel(method=-4)[source]
Select a model from the bayes.dat samples read with readBayesModels() function.
- Parameters:
method Specify the row from the bayes.dat from which a chires.dat file should be calculated. Special values are -4 to use the maximum likelihood row (default).
- set_cosmology(H0, omegaM, omegaX, wX=-1, wa=0, kcourb=0, model=1)[source]
Initialize internal cosmology variables
Parameters
- model
Several models are implemented (see cosmoratio.c). model=1 is Lambda CDM model.
- set_field(field, reset_grid=True)[source]
Set the field section (internal F global variable)
Parameter
- fieldlist or scalar
If a list, set the xmin, xmax, ymin, ymax value. If a scalar, set the dmax value centered on the lens
- property source
Return the source array global variable
- write_bayes_header()[source]
Write the header of the bayes.dat file. Open the burnin.dat and bayes.dat files on the disk.
Return the number of free parameters/dimensions
- lenstool.lenstool.sdarray_ndarray(array, nrows, ncols)[source]
Convert Lenstool square_double array to numpy.ndarray. The memory content of the input sd_array is freed upon the call to this function.
Parameters
array Object of ctypes type square_double_t
nrows, ncols Size of the first and second dimensions
Returns
Array of type numpy.ndarray with shape (nrows, ncols)
lenstool.optimize module
lenstool.pcl module
lenstool.plot module
lenstool.potentials module
- class lenstool.potentials.dpie(x, y, emass, theta, z, sigma, rc=None, rcut=None, rckpc=None, rcutkpc=None)[source]
Bases:
pot