handling Multidark simulations¶

This class handles the MultiDark data

class MultiDark.MultiDark¶

The class MultiDark is a wrapper to handle Multidark simulations results / outputs.

class MultiDark.MultiDarkSimulation(Lbox=<Quantity 2500.0 Mpc>, wdir='/data2/DATA/eBOSS/Multidark-lightcones/', boxDir='MD_2.5Gpc', snl=array([], dtype=float64), zsl=None, zArray=array([ 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1., 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2., 2.1, 2.2, 2.3]), Hbox=<Quantity 67.77 km / (Mpc s)>, Melement=23593750000.0)[source]¶

Loads the environement proper to the Multidark simulations. This is the fixed framework of the simulation.

Parameters:

Parameters:	Lbox – length of the box in Mpc/h wdir – Path to the multidark lightcone directory boxDir – box directory name snl – list of snapshots available zsl – list of redshift corresponding to the snapshots zArray – redshift array to be considered to interpolate the redshift – distance conversion Hbox – Hubble constant at redshift 0 of the box Melement – Mass of the resolution element in solar masses. columnDict – dictionnary to convert column name into the index to find it in the snapshots

Lbox – length of the box in Mpc/h
wdir – Path to the multidark lightcone directory
boxDir – box directory name
snl – list of snapshots available
zsl – list of redshift corresponding to the snapshots
zArray – redshift array to be considered to interpolate the redshift – distance conversion
Hbox – Hubble constant at redshift 0 of the box
Melement – Mass of the resolution element in solar masses.
columnDict – dictionnary to convert column name into the index to find it in the snapshots

combinesDoubleDistributionFunction(ii, nameA, nameB, binsA, binsB, type='Central')[source]¶: Coombines the outputs of computeDoubleDistributionFunction. :param ii: index of the snapshot :param name: name of the quantity studies :param bins: bins the histogram was done with :param type: “Central” or “Satellite”

combinesSingleDistributionFunction(ii, name='Vpeak', bins=array([ 1.00000000e+00, 1.02329299e+00, 1.04712855e+00, 1.07151931e+00, 1.09647820e+00, 1.12201845e+00, 1.14815362e+00, 1.17489755e+00, 1.20226443e+00, 1.23026877e+00, 1.25892541e+00, 1.28824955e+00, 1.31825674e+00, 1.34896288e+00, 1.38038426e+00, 1.41253754e+00, 1.44543977e+00, 1.47910839e+00, 1.51356125e+00, 1.54881662e+00, 1.58489319e+00, 1.62181010e+00, 1.65958691e+00, 1.69824365e+00, 1.73780083e+00, 1.77827941e+00, 1.81970086e+00, 1.86208714e+00, 1.90546072e+00, 1.94984460e+00, 1.99526231e+00, 2.04173794e+00, 2.08929613e+00, 2.13796209e+00, 2.18776162e+00, 2.23872114e+00, 2.29086765e+00, 2.34422882e+00, 2.39883292e+00, 2.45470892e+00, 2.51188643e+00, 2.57039578e+00, 2.63026799e+00, 2.69153480e+00, 2.75422870e+00, 2.81838293e+00, 2.88403150e+00, 2.95120923e+00, 3.01995172e+00, 3.09029543e+00, 3.16227766e+00, 3.23593657e+00, 3.31131121e+00, 3.38844156e+00, 3.46736850e+00, 3.54813389e+00, 3.63078055e+00, 3.71535229e+00, 3.80189396e+00, 3.89045145e+00, 3.98107171e+00, 4.07380278e+00, 4.16869383e+00, 4.26579519e+00, 4.36515832e+00, 4.46683592e+00, 4.57088190e+00, 4.67735141e+00, 4.78630092e+00, 4.89778819e+00, 5.01187234e+00, 5.12861384e+00, 5.24807460e+00, 5.37031796e+00, 5.49540874e+00, 5.62341325e+00, 5.75439937e+00, 5.88843655e+00, 6.02559586e+00, 6.16595002e+00, 6.30957344e+00, 6.45654229e+00, 6.60693448e+00, 6.76082975e+00, 6.91830971e+00, 7.07945784e+00, 7.24435960e+00, 7.41310241e+00, 7.58577575e+00, 7.76247117e+00, 7.94328235e+00, 8.12830516e+00, 8.31763771e+00, 8.51138038e+00, 8.70963590e+00, 8.91250938e+00, 9.12010839e+00, 9.33254301e+00, 9.54992586e+00, 9.77237221e+00, 1.00000000e+01, 1.02329299e+01, 1.04712855e+01, 1.07151931e+01, 1.09647820e+01, 1.12201845e+01, 1.14815362e+01, 1.17489755e+01, 1.20226443e+01, 1.23026877e+01, 1.25892541e+01, 1.28824955e+01, 1.31825674e+01, 1.34896288e+01, 1.38038426e+01, 1.41253754e+01, 1.44543977e+01, 1.47910839e+01, 1.51356125e+01, 1.54881662e+01, 1.58489319e+01, 1.62181010e+01, 1.65958691e+01, 1.69824365e+01, 1.73780083e+01, 1.77827941e+01, 1.81970086e+01, 1.86208714e+01, 1.90546072e+01, 1.94984460e+01, 1.99526231e+01, 2.04173794e+01, 2.08929613e+01, 2.13796209e+01, 2.18776162e+01, 2.23872114e+01, 2.29086765e+01, 2.34422882e+01, 2.39883292e+01, 2.45470892e+01, 2.51188643e+01, 2.57039578e+01, 2.63026799e+01, 2.69153480e+01, 2.75422870e+01, 2.81838293e+01, 2.88403150e+01, 2.95120923e+01, 3.01995172e+01, 3.09029543e+01, 3.16227766e+01, 3.23593657e+01, 3.31131121e+01, 3.38844156e+01, 3.46736850e+01, 3.54813389e+01, 3.63078055e+01, 3.71535229e+01, 3.80189396e+01, 3.89045145e+01, 3.98107171e+01, 4.07380278e+01, 4.16869383e+01, 4.26579519e+01, 4.36515832e+01, 4.46683592e+01, 4.57088190e+01, 4.67735141e+01, 4.78630092e+01, 4.89778819e+01, 5.01187234e+01, 5.12861384e+01, 5.24807460e+01, 5.37031796e+01, 5.49540874e+01, 5.62341325e+01, 5.75439937e+01, 5.88843655e+01, 6.02559586e+01, 6.16595002e+01, 6.30957344e+01, 6.45654229e+01, 6.60693448e+01, 6.76082975e+01, 6.91830971e+01, 7.07945784e+01, 7.24435960e+01, 7.41310241e+01, 7.58577575e+01, 7.76247117e+01, 7.94328235e+01, 8.12830516e+01, 8.31763771e+01, 8.51138038e+01, 8.70963590e+01, 8.91250938e+01, 9.12010839e+01, 9.33254301e+01, 9.54992586e+01, 9.77237221e+01, 1.00000000e+02, 1.02329299e+02, 1.04712855e+02, 1.07151931e+02, 1.09647820e+02, 1.12201845e+02, 1.14815362e+02, 1.17489755e+02, 1.20226443e+02, 1.23026877e+02, 1.25892541e+02, 1.28824955e+02, 1.31825674e+02, 1.34896288e+02, 1.38038426e+02, 1.41253754e+02, 1.44543977e+02, 1.47910839e+02, 1.51356125e+02, 1.54881662e+02, 1.58489319e+02, 1.62181010e+02, 1.65958691e+02, 1.69824365e+02, 1.73780083e+02, 1.77827941e+02, 1.81970086e+02, 1.86208714e+02, 1.90546072e+02, 1.94984460e+02, 1.99526231e+02, 2.04173794e+02, 2.08929613e+02, 2.13796209e+02, 2.18776162e+02, 2.23872114e+02, 2.29086765e+02, 2.34422882e+02, 2.39883292e+02, 2.45470892e+02, 2.51188643e+02, 2.57039578e+02, 2.63026799e+02, 2.69153480e+02, 2.75422870e+02, 2.81838293e+02, 2.88403150e+02, 2.95120923e+02, 3.01995172e+02, 3.09029543e+02, 3.16227766e+02, 3.23593657e+02, 3.31131121e+02, 3.38844156e+02, 3.46736850e+02, 3.54813389e+02, 3.63078055e+02, 3.71535229e+02, 3.80189396e+02, 3.89045145e+02, 3.98107171e+02, 4.07380278e+02, 4.16869383e+02, 4.26579519e+02, 4.36515832e+02, 4.46683592e+02, 4.57088190e+02, 4.67735141e+02, 4.78630092e+02, 4.89778819e+02, 5.01187234e+02, 5.12861384e+02, 5.24807460e+02, 5.37031796e+02, 5.49540874e+02, 5.62341325e+02, 5.75439937e+02, 5.88843655e+02, 6.02559586e+02, 6.16595002e+02, 6.30957344e+02, 6.45654229e+02, 6.60693448e+02, 6.76082975e+02, 6.91830971e+02, 7.07945784e+02, 7.24435960e+02, 7.41310241e+02, 7.58577575e+02, 7.76247117e+02, 7.94328235e+02, 8.12830516e+02, 8.31763771e+02, 8.51138038e+02, 8.70963590e+02, 8.91250938e+02, 9.12010839e+02, 9.33254301e+02, 9.54992586e+02, 9.77237221e+02, 1.00000000e+03, 1.02329299e+03, 1.04712855e+03, 1.07151931e+03, 1.09647820e+03, 1.12201845e+03, 1.14815362e+03, 1.17489755e+03, 1.20226443e+03, 1.23026877e+03, 1.25892541e+03, 1.28824955e+03, 1.31825674e+03, 1.34896288e+03, 1.38038426e+03, 1.41253754e+03, 1.44543977e+03, 1.47910839e+03, 1.51356125e+03, 1.54881662e+03, 1.58489319e+03, 1.62181010e+03, 1.65958691e+03, 1.69824365e+03, 1.73780083e+03, 1.77827941e+03, 1.81970086e+03, 1.86208714e+03, 1.90546072e+03, 1.94984460e+03, 1.99526231e+03, 2.04173794e+03, 2.08929613e+03, 2.13796209e+03, 2.18776162e+03, 2.23872114e+03, 2.29086765e+03, 2.34422882e+03, 2.39883292e+03, 2.45470892e+03, 2.51188643e+03, 2.57039578e+03, 2.63026799e+03, 2.69153480e+03, 2.75422870e+03, 2.81838293e+03, 2.88403150e+03, 2.95120923e+03, 3.01995172e+03, 3.09029543e+03]), type='Central')[source]¶: Coombines the outputs of computeSingleDistributionFunction. :param ii: index of the snapshot :param name: name of the quantity studies :param bins: bins the histogram was done with :param type: “Central” or “Satellite”

compute2PCF(catalogList, vmin=65, rmax=200, dlogBin=0.05, Nmax=4000000.0, dr=1.0, name='')[source]¶: Extracts the 2PCF out of a catalog of halos :param catalog: where the catalog is :param vmin: minimum circular velocity. :param dlogBin: bin width. :param rmax: maximum distance

compute2PCF_MASS(catalogList, rmax=200, dr=0.1, vmin=9, dlogBin=0.05, Nmax=2000000.0, name='')[source]¶: Extracts the 2PCF out of a catalog of halos :param catalog: where the catalog is :param vmin: minimum circular velocity. :param dlogBin: bin width. :param rmax: maximum distance

computeDoubleDistributionFunction(ii, nameA, nameB, binsA, binsB, Mfactor=100.0)[source]¶: Extracts the distributions of two quantity and their correlation ‘name’ out of all snapshots of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param name: name of the quantity of interest, mass, velocity. :param index: of the quantity of interest in the snapshots. :param bins: binning scheme to compute the historgram. :param Mfactor: only halos with Mvir > Mfact* Melement are used.

computeSingleDistributionFunction(ii, name, bins, Mfactor=10.0)[source]¶: Extracts the distribution of quantity ‘name’ out of all snapshots of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param name: name of the quantity of interest, mass, velocity. :param index: of the quantity of interest in the snapshots. :param bins: binning scheme to compute the historgram. :param Mfactor: only halos with Mvir > Mfact* Melement are used.

computeSingleDistributionFunctionJKresampling(fileList, rootname, name, bins, Ljk=100.0, overlap=1.0)[source]¶: Extracts the distribution of quantity ‘name’ out of all snapshots of the Multidark simulation. Resamples the box in smaller boxes of length Ljk in Mpc/h :param ii: index of the snapshot in the list self.snl :param name: name of the quantity of interest, mass, velocity. :param index: of the quantity of interest in the snapshots. :param bins: binning scheme to compute the historgram. :param Ljk: length of the resampled box :param overlap: allowed overlap between resampled realizations : 1 = no overlap 2 : 50% overlap ...

computeSingleDistributionFunctionV2(fileList, rootname, name, bins)[source]¶: Extracts the distribution of quantity ‘name’ out of all snapshots of the Multidark simulation. Resamples the box in smaller boxes of length Ljk in Mpc/h :param ii: index of the snapshot in the list self.snl :param name: name of the quantity of interest, mass, velocity. :param index: of the quantity of interest in the snapshots. :param bins: binning scheme to compute the historgram. :param Ljk: length of the resampled box :param overlap: allowed overlap between resampled realizations : 1 = no overlap 2 : 50% overlap ...

cornerLCpositionCatalog(ii, DMIN=0.0, DMAX=1000.0, vmin=190, vmax=100000, NperBatch=10000000)[source]¶: Extracts the positions and velocity out of a snapshot of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param DMIN:maximum distance for the pointto be included. :param DMAX:maximum distance for the pointto be included. :param vmin: name of the quantity of interest, mass, velocity. :param vmax: of the quantity of interest in the snapshots. :param NperBatch: number of line per fits file, default: 1000000

writePositionCatalogPM(ii, vmin=30.0, mmin=100000000, NperBatch=20000000)[source]¶: Extracts the positions and velocity out of a snapshot of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param vmin: name of the quantity of interest, mass, velocity. :param vmax: of the quantity of interest in the snapshots. :param NperBatch: number of line per fits file, default: 1000000

writePositionCatalogVmax(ii, vmin=190, vmax=10000, NperBatch=10000000)[source]¶: Extracts the positions and velocity out of a snapshot of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param vmin: name of the quantity of interest, mass, velocity. :param vmax: of the quantity of interest in the snapshots. :param NperBatch: number of line per fits file, default: 1000000

writePositionCatalogVmaxM200c(ii, vmin=190, vmax=10000, NperBatch=10000000)[source]¶: Extracts the positions and velocity out of a snapshot of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param vmin: name of the quantity of interest, mass, velocity. :param vmax: of the quantity of interest in the snapshots. :param NperBatch: number of line per fits file, default: 1000000

writeSAMcatalog(ii, mmin=100000000, NperBatch=2000000)[source]¶: Extracts the positions and mass out of a snapshot of the Multidark simulation. :param ii: index of the snapshot in the list self.snl :param vmin: name of the quantity of interest, mass, velocity. :param vmax: of the quantity of interest in the snapshots. :param NperBatch: number of line per fits file, default: 1000000