.. _physics-surveyGeometry: Survey Geometry =============== Class providing galaxy survey geometries---the selection function of an observed galaxy sample, including the survey solid angle, the minimum and maximum detection distances as a function of galaxy mass, luminosity, or star formation rate, and the 3-D window functions needed for clustering analyses. These quantities are used in the on-the-fly output analysis classes to compute galaxy stellar mass functions, luminosity functions, and two-point correlation functions by applying the appropriate :math:`V_\mathrm{max}` weighting and survey volume corrections. **Default implementation:** ``surveyGeometryLiWhite2009SDSS`` Methods ------- ``fieldCount`` → ``integer`` Returns the number of distinct fields included in the survey. ``distanceMinimum`` → ``double precision`` Returns the minimum distance (in Mpc) at which a galaxy of the specified ``mass`` (in :math:`\mathrm{M}_\odot`) would be included in the survey. * ``double precision, intent(in ), optional :: mass , magnitudeAbsolute, luminosity, starFormationRate`` * ``integer , intent(in ), optional :: field`` ``distanceMaximum`` → ``double precision`` Returns the maximum distance (in Mpc) at which a galaxy of the specified ``mass`` (in :math:`\mathrm{M}_\odot`) could be detected. * ``double precision, intent(in ), optional :: mass , magnitudeAbsolute, luminosity, starFormationRate`` * ``integer , intent(in ), optional :: field`` ``solidAngle`` → ``double precision`` Return the solid angle (in steradians) of the survey. * ``integer, intent(in ), optional :: field`` ``volumeMaximum`` → ``double precision`` Returns the maximum volume (in Mpc\ :math:`^3`) at which a galaxy of the specified ``mass`` (in :math:`\mathrm{M}_\odot`) could be detected. * ``double precision, intent(in ) :: mass`` * ``integer , intent(in ), optional :: field`` ``windowFunctionAvailable`` → ``logical`` Returns true if survey 3-D window functions are available. ``angularPowerAvailable`` → ``logical`` Returns true if angular power spectrum of survey window function is available. ``windowFunctions`` → ``void`` Returns the window functions on a grid of the specified size (``gridCount`` cells in each dimension) for galaxies of the specified ``mass1`` and ``mass2`` (in :math:`\mathrm{M}_\odot`). The ``boxLength`` should be set to an appropriate value to fully enclose (with sufficient buffering to allow for Fourier transformation) the two window functions. * ``double precision , intent(in ) :: mass1 , mass2`` * ``integer , intent(in ) :: gridCount`` * ``double precision , intent( out) :: boxLength`` * ``complex (c_double_complex), intent( out), dimension(gridCount,gridCount,gridCount) :: windowFunction1, windowFunction2`` ``angularPower`` → ``double precision`` Return :math:`C^{ij}_\ell`, where :math:`(2\ell+1) C^{ij}_\ell = \sum_{m=-\ell}^{+\ell} \Psi^i_{\ell m} \Psi^{j*}_{\ell m}`, and :math:`\Psi^i_{\ell m}` are the coefficients of the spherical harmonic expansion of the :math:`i^\mathrm{th}` field. * ``integer, intent(in ) :: i,j,l`` ``angularPowerMaximumDegree`` → ``integer`` Return the maximum degree, :math:`\ell_\mathrm{max}`, for which the angular power is available. ``pointIncluded`` → ``logical`` Return true if the given Cartesian point lies within the survey bounds for the given mass limit. * ``double precision, intent(in ), dimension(3) :: point`` * ``double precision, intent(in ) :: mass`` .. _physics-surveyGeometryBaldry2012GAMA: ``surveyGeometryBaldry2012GAMA`` -------------------------------- A survey geometry class that describes the survey geometry of :cite:t:`baldry_galaxy_2012`. For the angular mask we use the specifications of the G09, G12, and G15 fields given by :cite:t:`driver_galaxy_2011` to construct mangle polygon files. To determine the depth as a function of stellar mass, we make use of the publicly available tabulated mass function, :math:`\phi`, and number of galaxies per bin, :math:`N`. The effective volume of each bin is found as :math:`V_i = N_i/\phi_i\Delta\log_{10}M_\star`, where :math:`\Delta\log_{10}M_\star` is the width of the bin. The GAMA survey consists of three fields, each of the same solid angle, but with differing depths. We assume that the relative depths in terms of stellar mass scale with the depth in terms of flux. Given this assumption, these volumes are converted to maximum distances in each field using the solid angle quoted above. The resulting mass vs. distance relation in each field is fit with a :math:`1^\mathrm{st}`-order polynomial in log-log space over the range where the maximum volume is limited by the survey depth and not by the imposed :math:`z=0.06` upper limit to redshift. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum distance at which such a galaxy would be included in the sample. Points indicate results from GAMA, while the line shows a polynomial fit: .. math:: :label: eq-BaldryDepthPolynomial \log_{10} \left[ {D_\mathrm{max}(M_\star) \over \hbox{Mpc}}\right] = \left\{ \begin{array}{ll} -0.521 + 0.319m & \hbox{fields G09/G15} \\ -0.361 + 0.319m & \hbox{field G12} \end{array} \right. where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. .. figure:: /_figures/BaldryGAMAMassDistanceRelation.png :name: fig-BaldryGAMADepthFit The maximum distance at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`baldry_galaxy_2012`. Points show the results obtained from data provided by Baldry, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-BaldryDepthPolynomial`). Note that above :math:`10^9\mathrm{M}_\odot` the distance is limited by the imposed upper limit of :math:`z=0.06` in the GAMA sample---the polynomial fit does not consider these points. .. _physics-surveyGeometryBernardi2013SDSS: ``surveyGeometryBernardi2013SDSS`` ---------------------------------- A survey geometry class that describes the survey geometry of :cite:t:`bernardi_massive_2013`. For the angular mask, we make use of the :term:`mangle` polygon file provided by the :term:`mangle` project\footnoteSpecifically, `https://zenodo.org/records/10998446/files/sdss_dr72safe0_res6d.pol.gz `_. The solid angle of this mask, computed using the :term:`mangle` ``harmonize`` command is 2.232262776405 sr. To determine the depth as a function of stellar mass, we make use of results provided by M. Bernardi (private communication), giving the mean maximum volume, :math:`V_\mathrm{max}`, as a function of stellar mass for galaxies in this sample. These maximum volumes are converted to maximum distances using the solid angle quoted above. The results mass vs. distance relation is fit with a :math:`5^\mathrm{th}`-order polynomial. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum distance at which such a galaxy would be included in the sample. Points indicate results from Bernardi, while the line shows a polynomial fit: .. math:: :label: eq-BernardiDepthPolynomial \log_{10} \left[ {D_\mathrm{max}(M_\star) \over \hbox{Mpc}}\right] = 1282.11+m (-626.644+m (122.091+m (-11.8431+m (0.572399+m (-0.0110301))))) where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. .. figure:: /_figures/BernardiSDSSMassLuminosityRelation.png :name: fig-BernardiSDSSDepthFit The maximum distance at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`bernardi_massive_2013`. Points show the results obtained from data provided by Bernardi, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-BernardiDepthPolynomial`). .. _physics-surveyGeometryCaputi2011UKIDSSUDS: ``surveyGeometryCaputi2011UKIDSSUDS`` ------------------------------------- A survey geometry class which implements the UKIDSS UDS survey used by :cite:t:`caputi_stellar_2011`. The survey window function is determined from the set of galaxy positions provided by Caputi (private communication), by finding a suitable bounding box and then cutting out empty regions (corresponding to regions that were removed around bright stars). A set of random points are then found within this mask and are used to find the Fourier transform of the survey volume. To estimate the depth of the :cite:t:`caputi_stellar_2011` sample as a function of galaxy stellar mass we make use of semi-analytic models in the Millennium Database. Specifically, we use the :term:`SAM` of :cite:t:`guo_dwarf_2011` and :cite:t:`henriques_confronting_2012` specifically the ``Guo2010a..MR`` and ``Henriques2012a.wmap1.BC03_001`` tables in the Millennium Database. For each snapshot in the database, we extract the stellar masses and observed-frame IRAC 4.5\ :math:`\mu`\ m apparent magnitudes (including dust extinction), and determine the median apparent magnitude as a function of stellar mass. Using the limiting apparent magnitude of the :cite:t:`caputi_stellar_2011` sample, :math:`i_{4.5}=24`, we infer the corresponding absolute magnitude at each redshift and, using our derived apparent magnitude--stellar mass relation, infer the corresponding stellar mass. The end result of this procedure is the limiting stellar mass as a function of redshift, accounting for k-corrections, evolution, and the effects of dust. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum redshift at which such a galaxy would be included in the sample. Points indicate measurements from the :term:`SAM`, while the line shows a polynomial fit: .. math:: :label: eq-UKIDSSUDSDepthPolynomial z(M_\star) = -56.247 + 5.881 m, where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. .. figure:: /_figures/UKIDSSUDSMassLuminosityRelation.png :name: fig-UKIDSSUDSMassRedshift The maximum redshift at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`caputi_stellar_2011`. Points show the results obtained using the :cite:t:`henriques_confronting_2012` model from the Millennium Database, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-UKIDSSUDSDepthPolynomial`). **Parameters** * ``[redshiftBin]`` — The redshift bin (0, 1, or 2) of the :cite:t:`caputi_stellar_2011` to use. .. _physics-surveyGeometryCombined: ``surveyGeometryCombined`` -------------------------- Implements a survey geometry which combines multiple other surveys. .. _physics-surveyGeometryDavidzon2013VIPERS: ``surveyGeometryDavidzon2013VIPERS`` ------------------------------------ A survey geometry class that describes the survey geometry of :cite:t:`davidzon_vimos_2013`. For the angular mask, we make use of ``mangle`` polygon files provided by I. Davidzon (private communication) corresponding to the VIPERS fields. The solid angle of each mask is computed using the ``mangle`` ``harmonize`` command. To determine the depth as a function of stellar mass, we make use of the tabulated mass function, :math:`\phi`, and number of galaxies per bin, :math:`N`, supplied by I. Davidzon (private communication). The effective volume of each bin is found as :math:`V_i = N_i/f_\mathrm{complete}\phi_i\Delta\log_{10}M_\star`, where :math:`\Delta\log_{10}M_\star` is the width of the bin, and :math:`f_\mathrm{complete}` is the completeness of the survey, estimated to be approximately 40% :cite:p:`guzzo_vimos_2013`. These volumes are converted to maximum distances in each field using the survey solid angle. The resulting mass vs. distance relation in each field is fit with a :math:`1^\mathrm{st}`-order polynomial in log-log space over the range where the maximum volume is limited by the survey depth and not by the imposed upper limit to redshift. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum distance at which such a galaxy would be included in the sample. Points indicate results from VIPERS, while the lines show polynomial fits: .. math:: :label: eq-DavidzonDepthPolynomial \log_{10} \left[ {D_\mathrm{max}(M_\star) \over \hbox{Mpc}}\right] = \left\{ \begin{array}{ll} 3.207 + 0.0124m & 0.5 < z < 0.6 \\ 3.148 + 0.0268m & 0.6 < z < 0.8 \\ 3.207 + 0.0273m & 0.8 < z < 1.0 \end{array} \right. where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. .. figure:: /_figures/DavidzonVIPERSMassDistanceRelation.png :name: fig-Davidzon2013DepthFit The maximum distance at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`davidzon_vimos_2013`. Points show the results obtained from data provided by Davidzon, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-DavidzonDepthPolynomial`). Note that at high masses the distance is limited by the imposed upper limit---the polynomial fit does not consider these points. **Parameters** * ``[redshiftBin]`` — The redshift bin (0, 1, 2) of the :cite:t:`davidzon_vimos_2013` mass function to use. .. _physics-surveyGeometryFullSky: ``surveyGeometryFullSky`` ------------------------- Implements a survey geometry covering the full sky over a specified distance or redshift range. The survey volume is bounded by ``[redshiftMinimum]`` and ``[redshiftMaximum]``, with all sky positions included, making this suitable for theoretical or volume-limited analyses without angular masking. **Parameters** * ``[redshiftMinimum]`` (default ``0.0d0``) — The minimum redshift of the full-sky survey volume; sources below this redshift are excluded from the survey sample. * ``[redshiftMaximum]`` (default ``huge(1.0d0)``) — The maximum redshift of the full-sky survey volume; sources above this redshift are excluded from the survey sample. .. _physics-surveyGeometryGunawardhana2013SDSS: ``surveyGeometryGunawardhana2013SDSS`` -------------------------------------- Implements the geometry of the SDSS survey of :cite:t:`gunawardhana_galaxy_2013`. .. _physics-surveyGeometryHearin2014SDSS: ``surveyGeometryHearin2014SDSS`` -------------------------------- Implements the survey geometry of the SDSS sample used by :cite:t:`hearin_dark_2013`. .. _physics-surveyGeometryKelvin2014GAMAnear: ``surveyGeometryKelvin2014GAMAnear`` ------------------------------------ Implements the geometry of the GAMAnear survey of :cite:t:`kelvin_galaxy_2014-1`. .. _physics-surveyGeometryLiWhite2009SDSS: ``surveyGeometryLiWhite2009SDSS`` --------------------------------- A survey geometry class that describes the survey geometry of :cite:t:`li_distribution_2009`. For the angular mask, we make use of the catalog of random points within the survey footprint provided by the NYU-VAGC\footnoteSpecifically, `https://zenodo.org/records/10257229/files/lss_random-0.dr72.dat `_ (which is a copy of the dataset originally found at the, now defunct, URL ``http://sdss.physics.nyu.edu/lss/dr72/random/lss_random-0.dr72.dat``). (:cite:author:`blanton_new_2005` :cite:year:`blanton_new_2005`; see also :cite:t:`adelman-mccarthy_sixth_2008,padmanabhan_improved_2008`). :cite:t:`li_distribution_2009` consider only the main, contiguous region and so we keep only those points which satisfy RA\ :math:`>100^\circ`, RA\ :math:`<300^\circ`, and RA\ :math:`<247^\circ` or :math:`\delta< 51^\circ`. When the survey window function is needed, these points are used to determine which elements of a 3D grid fall within the window function. To estimate the depth of the :cite:t:`li_distribution_2009` sample as a function of galaxy stellar mass we make use of semi-analytic models in the Millennium Database. Specifically, we use the :term:`SAM` of :cite:author:`de_lucia_hierarchical_2007` (:cite:year:`de_lucia_hierarchical_2007`; specifically the ``millimil..DeLucia2006a`` and ``millimil..DeLucia2006a_sdss2mass`` tables in the Millennium Database). For each snapshot in the database, we extract the stellar masses and observed-frame SDSS r-band absolute magnitudes (including dust extinction), and determine the median absolute magnitude as a function of stellar mass. Using the limiting apparent magnitude of the :cite:t:`li_distribution_2009` sample, :math:`r=17.6`, we infer the corresponding absolute magnitude at each redshift and, using our derived absolute magnitude--stellar mass relation, infer the corresponding stellar mass. The end result of this procedure is the limiting stellar mass as a function of redshift, accounting for k-corrections, evolution, and the effects of dust. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum redshift at which such a galaxy would be included in the sample. Points indicate measurements from the :term:`SAM`, while the line shows a polynomial fit: .. math:: :label: eq-DepthPolynomial z(M_\star) & = -5.950 + 2.638 m - 0.4211 m^2 \nonumber \\ & + 2.852\times 10^{-2} m^3 - 6.783 \times 10^{-4} m^4, where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. We adopt a solid angle of :math:`2.1901993` sr :cite:p:`percival_shape_2007` for the sample. .. figure:: /_figures/SDSSMassLuminosityRelation.png :name: fig-SDSSDepthFit The maximum redshift at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`li_distribution_2009`. Points show the results obtained using the :cite:t:`de_lucia_hierarchical_2007` model from the Millennium Database, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-DepthPolynomial`). **(Default implementation)** **Parameters** * ``[redshiftMinimum]`` (default ``0.0d0``) — The minimum redshift of the :cite:t:`li_distribution_2009` survey volume; sources below this redshift are excluded from the survey sample. * ``[redshiftMaximum]`` (default ``huge(1.0d0)``) — The maximum redshift of the :cite:t:`li_distribution_2009` survey volume; sources above this redshift are excluded from the survey sample. .. _physics-surveyGeometryLocalGroupClassical: ``surveyGeometryLocalGroupClassical`` ------------------------------------- Implements a survey geometry corresponding to the detectability of classical Local Group satellite galaxies, defined by :term:`mangle` polygon sky coverage. The maximum detection distance is set by ``[distanceMaximumSurvey]``, and a stellar mass threshold ``[massThreshold]`` separates classical from ultra-faint satellites. **Parameters** * ``[distanceMaximumSurvey]`` (default ``300.0d-3``) — The maximum distance for the sample of classical Local Group galaxies. * ``[massThreshold]`` (default ``1.0d5``) — The minimum stellar mass for a classical Local Group dwarf galaxy. .. _physics-surveyGeometryLocalGroupDES: ``surveyGeometryLocalGroupDES`` ------------------------------- Implements the angular footprint and maximum detection distance of the Dark Energy Survey (DES) as applied to Local Group dwarf galaxy searches. The survey volume is defined by :term:`mangle` polygon files, with the maximum survey distance set by ``[distanceMaximumSurvey]``. **Parameters** * ``[distanceMaximumSurvey]`` (default ``300.0d-3``) — The maximum distance at which galaxies are to be included in the survey. .. _physics-surveyGeometryLocalGroupSDSS: ``surveyGeometryLocalGroupSDSS`` -------------------------------- Implements the angular footprint of the SDSS survey adapted for Local Group dwarf galaxy detectability, with the maximum survey depth set by ``[distanceMaximumSurvey]``. This geometry is used to model the detection volume of faint dwarf galaxies in the nearby Universe. **Parameters** * ``[distanceMaximumSurvey]`` (default ``300.0d-3``) — The maximum distance at which galaxies are to be included in the survey. .. _physics-surveyGeometryMangle: ``surveyGeometryMangle`` ------------------------ Implements an abstract survey geometry using :term:`mangle` polygons. **Methods** * ``mangleDirectory`` — Return the directory containing :term:`mangle` files for this survey geometry. * ``mangleFiles`` — Return array of :term:`mangle` filenames for this survey geometry. * ``initialize`` — Initialize an instance of the :term:`mangle` survey geometry class. .. _physics-surveyGeometryMartin2010ALFALFA: ``surveyGeometryMartin2010ALFALFA`` ----------------------------------- A survey geometry class that describes the survey geometry of :cite:t:`martin_arecibo_2010`. For the angular mask we use the three disjoint regions defined by 07\ :math:`^\mathrm{h}`\ 30\ :math:`^\mathrm{m}` :math:`<` R.A. :math:`<` 16\ :math:`^\mathrm{h}`\ 30\ :math:`^\mathrm{m}`, +04\ :math:`^\circ` :math:`<` decl. :math:`<` +16\ :math:`^\circ`, and +24\ :math:`^\circ` :math:`<` decl. :math:`<` +28\ :math:`^\circ` and 22\ :math:`^\mathrm{h}` :math:`<` R.A. :math:`<` 03\ :math:`^\mathrm{h}`, +14\ :math:`^\circ` :math:`<` decl. :math:`<` +16\ :math:`^\circ`, and +24\ :math:`^\circ` :math:`<` decl. :math:`<` +32\ :math:`^\circ` corresponding to the sample of :cite:t:`martin_arecibo_2010`. When the survey window function is needed we generate randomly distributed points within this angular mask and out to the survey depth. These points are used to determine which elements of a 3D grid fall within the window function. To estimate the depth of the :cite:t:`martin_arecibo_2010` sample as a function of galaxy HI mass we first infer the median line width corresponding to that mass. To do so, we have fit the median line width-mass relation from the :math:`\alpha.40` sample with power-law function as shown in Fig. :numref:`{number} `. We find that the median line width can be approximated by .. math:: :label: eq-ALFALFALineWidthMassRelation \log_{10} (W_\mathrm{50}/\hbox{km s}^{-1}) = c_0 + c_1 \log_10(M_\mathrm{HI}/\mathrm{M}_\odot), with :math:`c_0=-0.770` and :math:`c_1=0.315`. Given the line width, the corresponding integrated flux limit, :math:`S_\mathrm{int}`, for a signal-to-noise of :math:`6.5` is inferred using equation (A1) of :cite:t:`haynes_arecibo_2011`. Finally, this integrated flux limit is converted to maximum distance at which the source could be detected using the expression given in the text of section 2.2 of :cite:t:`martin_arecibo_2010`: .. math:: M_\mathrm{HI} = 2.356\times10^5 \left({D\over \hbox{Mpc}}\right)^2 \left({S_\mathrm{int}\over\hbox{Jy km s}^{-1}}\right). .. figure:: /_figures/alfalfaHILineWidthMassRelation.png :name: fig-ALFALFALineWidthMassRelation HI line width vs. HI mass as measured from the :math:`\alpha.40` survey of :cite:t:`martin_arecibo_2010`. Red points with error bars show individual measurements, while the larger circles indicate the running median of these data. The green line is a power-law fit to the running median as described in eqn. (:eq:`eq-ALFALFALineWidthMassRelation`). .. _physics-surveyGeometryMonteroDorta2009SDSS: ``surveyGeometryMonteroDorta2009SDSS`` -------------------------------------- Implements the geometry of the SDSS survey as used by :cite:t:`montero-dorta_sdss_2009` for their luminosity function measurements. The effective survey volume is defined by redshift limits ``[redshiftMinimum]`` and ``[redshiftMaximum]`` and apparent magnitude limits in a specified photometric ``[band]``. **Parameters** * ``[band]`` — The band for which the survey geometry should be computed. .. _physics-surveyGeometryMoustakas2013PRIMUS: ``surveyGeometryMoustakas2013PRIMUS`` ------------------------------------- A survey geometry class that describes the survey geometry of :cite:t:`moustakas_primus:_2013`. For the angular mask, we make use of :term:`mangle` polygon files provided by J. Moustakas (private communication) corresponding to the give PRIMUS fields. The solid angle of each mask is computed using the :term:`mangle` ``harmonize`` command. To determine the depth as a function of stellar mass, we make use of completeness limits for "All" galaxies given in Table 2 of :cite:t:`moustakas_primus:_2013`. These are fit, for each field with a second order polynomial to give the limiting redshift as a function of stellar mass. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum redshift at which such a galaxy would be included in the sample. Points indicate results from :cite:t:`moustakas_primus:_2013`, while the line shows a polynomial fits: .. math:: :label: eq-MoustakasDepthPolynomial z_\mathrm{max}(M_\star) = +3.51+m(-0.941+m(+0.0651)) & \hbox{COSMOS} \\ z_\mathrm{max}(M_\star) = +2.46+m(-0.730+m(+0.0542)) & \hbox{XMM-SXDS} \\ z_\mathrm{max}(M_\star) = -3.60+m(+0.500+m(-0.0078)) & \hbox{XMM-CFHTLS} \\ z_\mathrm{max}(M_\star) = +5.87+m(-1.528+m(+0.0982)) & \hbox{CDFS} \\ z_\mathrm{max}(M_\star) = +6.87+m(-1.656+m(+0.1003)) & \hbox{ELAIS-S1} where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. .. figure:: /_figures/MoustakasPRIMUSMassRedshiftRelation.png :name: fig-MoustakasPRIMUSDepthFit The maximum distance at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`moustakas_primus:_2013`. Points show the results obtained from completeness limit data taken from Table 2 of :cite:t:`moustakas_primus:_2013`, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-MoustakasDepthPolynomial`). **Parameters** * ``[redshiftBin]`` — The redshift bin (0, 1, 2, 3, 4, 5, or 5) of the :cite:t:`moustakas_primus:_2013` mass function to use. .. _physics-surveyGeometryMuzzin2013ULTRAVISTA: ``surveyGeometryMuzzin2013ULTRAVISTA`` -------------------------------------- A survey geometry class that describes the survey geometry of :cite:t:`muzzin_evolution_2013`. For the angular mask, we generate a :term:`mangle` polygon file, by first defining a rectangle encompassing the bounds of the ULTRAVISTA field (:math:`149.373^\circ < \alpha < 150.779^\circ` and :math:`1.604^\circ < \delta < 2.81^\circ`). From this rectangle, we then remove circles of radii :math:`75^{\prime\prime}` around bright stars (i.e. those bright than 10\ :math:`^\mathrm{th}` and :math:`8^\mathrm{th}` magnitudes in the USNO and 2MASS star lists respectively) and radii :math:`30^{\prime\prime}` around medium stars (i.e. those bright than :math:`13^\mathrm{th}` and :math:`10.5^\mathrm{th}` magnitudes in the USNO and 2MASS star lists respectively). Finally, we mask regions of one detector for which 75% of pixels are dead by clipping pixels with weights below :math:`0.02` in the K\ :math:`_\mathrm{s}`-band weight map. These choices match those made in the ULTRAVISTA survey (A. Muzzin, private communication). The solid angle of each mask is computed using the :term:`mangle` ``harmonize`` command. To determine the depth as a function of stellar mass, we simply fit the `tabulated relations `_ provided by the ULTRAVISTA survey: .. math:: :label: eq-MuzzinDepthPolynomial z_\mathrm{max}(M_\star) = {-8364.45 + m (4331.82 + m (-896.596 + m (92.6999 + m (-4.78750 + m (0.0988215))))) \over 1 - \exp[(m-11.24)/0.02] } where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. .. figure:: /_figures/MuzzinULTRAVISTAMassRedshiftRelation.png :name: fig-MuzzinULTRAVISTADepthFit The maximum distance at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`muzzin_evolution_2013`. The dotted line shows the results obtained from the ULTRAVISTA survey :cite:p:`muzzin_evolution_2013`, while the solid line shows the polynomial fit to these results (given in eqn. :eq:`eq-MuzzinDepthPolynomial`). **Parameters** * ``[redshiftBin]`` — The redshift bin (0, 1, 2, 3, 4, 5, or 6) of the :cite:t:`muzzin_evolution_2013` mass function to use. .. _physics-surveyGeometryRandomPoints: ``surveyGeometryRandomPoints`` ------------------------------ Implements survey geometries defined by random points. **Methods** * ``randomsInitialize`` — Initialize arrays of random points to define the survey angular geometry. .. _physics-surveyGeometryTomczak2014ZFOURGE: ``surveyGeometryTomczak2014ZFOURGE`` ------------------------------------ A survey geometry class that describes the survey geometry of :cite:t:`tomczak_galaxy_2014`. For the angular mask, we make use of ``mangle`` polygon files constructed by hand using vertices matched approximately to the distribution of galaxies in the survey (positions of which were provided by R. Quadri; private communication). The solid angle of each mask is computed using the ``mangle`` ``harmonize`` command. To determine the depth as a function of stellar mass, we make use of the tabulated mass completeness limits as a function of redshift for ZFOURGE and NMBS fields provided by R. Quadri (private communication). These are fit with fourth-order polynomials. Figure :numref:`{number} ` shows the resulting relation between stellar mass and the maximum redshift at which such a galaxy would be included in the sample. Dotted lines indicate the tabulated result from ZFOURGE, while the lines show polynomial fits: .. math:: :label: eq-TomczakDepthPolynomial z_\mathrm{max}(M_\star) = \left\{ \begin{array}{ll} -114.66+m*(45.901+m*(-6.1617+m*(0.27822))) & \hbox{ZFOURGE fields} \\ -58.483+m*(20.250+m*(-2.3563+m*(0.092705))) & \hbox{NMBS fields} \end{array} \right. where :math:`m= \log_{10}(M_\star/\mathrm{M}_\odot)`. We use this polynomial fit to determine the depth of the sample as a function of stellar mass. .. figure:: /_figures/TomczakZFOURGEMassRedshiftRelation.png :name: fig-Tomczak2014DepthFit The maximum redshift at which a galaxy of given stellar mass can be detected in the sample of :cite:t:`tomczak_galaxy_2014`. Points show the results obtained from data provided by Davidzon, while the lines shows a polynomial fit to these results (given in eqn. :eq:`eq-TomczakDepthPolynomial`). **Parameters** * ``[redshiftBin]`` — The redshift bin (0, 1, 2, 3, 4, 5, 6, or 7) of the :cite:t:`tomczak_galaxy_2014` mass function to use.