GLAMER - II. Multiple-plane gravitational lensing

We present an extension to multiple planes of the gravitational lensing code GLAMER. The method entails projecting the mass in the observed light-cone on to a discrete number of lens planes and inverse ray-shooting from the image to the source plane. The mass on each plane can be represented as haloes, simulation particles, a projected mass map extracted form a numerical simulation or any combination of these. The image finding is done in a source-oriented fashion, where only regions of interest are iteratively refined on an initially coarse image plane grid. The calculations are performed in parallel on shared memory machines. The code is able to handle different types of analytic haloes (NFW, NSIE, power law, etc.), haloes extracted from numerical simulations and clusters constructed from semi-analytic models (MOKA). Likewise, there are several different options for modelling the source(s) which can be distributed throughout the light-cone. The distribution of matter in the light-cone can be either taken from a pre-existing N-body numerical simulations, from halo catalogues, or are generated from an analytic mass function. We present several tests of the code and demonstrate some of its applications such as generating mock images of galaxy and galaxy cluster lenses.