.. _buildini_help: .. |br| raw:: html
Usage of buildini_mamposst -------------------------- % ``buildini_mamposst`` -*options prefix* **Options on data**: ``-data`` *file* data filename is *file* ``-dir`` *directory* data directory is *directory* (relative or absolute path) ``-chaindir`` *directory* COSMOMC chain directory is *directory* relative to ``$HOME/NOSAVE/COSMOMC/chains`` [default: ``.``] ``-weightflag`` *flag* ``1``: Gaussian weights on tracers (with standard deviation taken from that of the line-of-sight velocities of the component), else no weights [default: ``0``] ``-distflag`` *flag* flag for redshift-independent distances [default: ``0``]: ``0``: ignore mu (distance modulus) |br| ``1``: weight by Gaussian(mu) |br| ``2``: weight by density Gaussian(mu) **Options on component-dependent models and fixed parameters**: ``-n`` *num* *comp1* [*comp2* ... *compn*] use *num* tracer components with names *comp1* ... *compn* (must be placed before other tracer options) ``-CofR`` *file* [*file2* ... *filen*] use files *file* etc. (with 2 header lines) for *R* Completeness(*R*) and MAMPOSSt will spline-interpolate [default: ``none``] ``-wtofR`` *file* [*file2* ... *filen*] use files *file* etc. (with 2 header lines) for *R* weight(*R*) and MAMPOSSt will spline-interpolate [default: ``none``] ``-tracermodel`` *model* [*model2* ... *modeln*] model(s) *model* etc. of visible tracer(s); available models are: ``Hernquist``: :math:`\rho \propto r^{-1}(r+a)^{-3}` (`Hernquist 1990 `_ model, :math:`r_{-2}/a=1/2`) |br| ``mHubble``: :math:`\rho \propto (r^2+a^2)^{-3/2}` (modified Hubble or non-truncated analytical `King 1962 `_ model, :math:`r_{-2}/a=\sqrt{2}`) |br| ``isothermal``: :math:`\rho \propto (r^2+a^2)^{-1}` (pseudo-isothermal) |br| ``Jaffe``: :math:`\rho \propto r^{-2}(r+a)^{-2}` (`Jaffe 1983 `_) |br| ``NFW``: :math:`\rho \propto r^{-1}(r+a)^{-2}` (`Navarro, Frenk & White 1996 `_ model, :math:`r_{-2}/a=1`) |br| ``Plummer``: :math:`\rho \propto (r^2+a^2)^{-5/2}` (`Plummer 1911 `_ model, :math:`r_{-2}/a=\sqrt{2/3}`) |br| ``gPlummer``: :math:`\rho \propto r^\gamma (r^2+a^2)^{-5/2-\gamma/2}` (generalized Plummer with free inner slope :math:`\gamma`, :math:`r_{-2}/a=\sqrt{(2+\gamma)/3}`) |br| ``PrugnielSimien``: :math:`\displaystyle \rho \propto x^{-p(n)}\,\exp \left[-b(n)\,\left({r\over R_{\rm eff}}\right)^{1/n}\right]` (`Prugniel & Simien 1997 `_ approximation to deprojected Sersic, where :math:`b(n)` is from analytical approximation of `Ciotti & Bertin 1999 `_, while :math:`p(n)` is given by `Lima Neto, Gerbal & Marquez 1999 `_) |br| ``-anismodel`` *model* [*model2* ... *modeln*] velocity anisotropy model(s) *model1* etc. of visible tracer(s); available models are: ``iso``: :math:`\beta = 0` (isotropic) |br| ``cst``: :math:`\beta` = constant |br| ``ML``: :math:`\displaystyle \beta = {1\over 2} {r\over r+r_\beta}` (`Mamon & Łokas 2005b `_) |br| ``OM``: :math:`\displaystyle \beta = {r^2\over r^2+r_\beta^2}` (`Osipkov 1979 `_; `Merritt 1985 `_) |br| ``gOM``: :math:`\displaystyle \beta = \beta_0 + (\beta_\infty-\beta_0)\,{r^2\over r^2+r_\beta^2}` (generalized Osipkov-Merritt) |br| |br| ``Tiret``: :math:`\displaystyle \beta = \beta_0 + (\beta_\infty-\beta_0)\,{r\over r+r_\beta}` (`Tiret et al. 2007 `_) ``-meanlrtr`` *meanlr* [*meanlr2* ... *meanlrn*] externally computed log_10 tracer radius (radii) (kpc); tracer radii are defined as: * scale radius *a* (``isothermal`` or ``Jaffe``) * effective radius :math:`R_{\rm eff}` (``PrugnielSimien``) * radius of logarithmic density slope :math:`-2` (other models, see ``-tracermodel``) ``-siglrtr`` *sigmalr* [*sigmalr2* ... *sigmalrn*] uncertainty of externally computed log_10 tracer radius (radii) (kpc) ``-rfidtr`` *r* [*r2* ... *rn*] tracer fiducial radius (radii) (where tracer mass is given) [kpc, default: ``1``] ``-Rminallow`` *Rmin* [*Rmin2* ... *Rminn*] minimum allowed projected radius (radii) for fits [kpc, ``0`` for no limits, default: ``30``] ``-Rmaxallow`` *Rmax* [*Rmax2* ... *Rmaxn*] maximum allowed projected radius (radii) for fits [kpc, ``0`` for no limits, default: ``0``] ``-avzmaxallow`` *vmax* [*vmax2* ... *vmaxn*] maximum allowed :math:`|\!` *v* :math:`{}_{\!\rm LOS}|` for fits (km/s, ``0`` for no limits, default: ``0``] ``-ev`` *errv* [*errv2* ... *errvn*] velocity error (in km/s) to be automatically inserted in data file if not present [default: not set] **Options on other models and fixed parameters**: ``-darkmodel`` *model* dark matter model; available models are: ``Burkert``: :math:`\rho \propto (r+a)^{-1}(r^2+a^2)^{-1}` (`Burkert 1995 `_ model, :math:`r_{-2}/a=(1-\sqrt{26/27})^{1/3}+(1+\sqrt{26/27})^{1/3})` |br| ``Einasto``: :math:`\rho \propto \exp[-b(n)\, r^{1/n}]` (`Einasto 1965 `_ model, :math:`r_{-2}/a = (2n)^n`) |br| ``Hernquist``: :math:`\rho \propto r^{-1}(r+a)^{-3}` (`Hernquist 1990 `_ model, :math:`r_{-2}/a=1/2`) |br| ``gHernquist``: :math:`\rho \propto r^\gamma (r+a)^{-4-\gamma}` (generalized Hernquist model with free inner slope :math:`\gamma`, :math:`r_{-2}/a=1+\gamma/2`) |br| ``mHubble``: :math:`\rho \propto (r^2+a^2)^{-3/2}` (modified Hubble or non-truncated analytical `King 1962 `_ model, :math:`r_{-2}/a=\sqrt{2}`) |br| ``isothermal``: :math:`\rho \propto (r^2+a^2)^{-1}` (pseudo-isothermal) |br| ``Jaffe``: :math:`\rho \propto r^{-2}(r+a)^{-2}` (`Jaffe 1983 `_) |br| ``Kazantzidis``: :math:`\rho \propto r^\gamma \exp(-r/a)` (`Kazantzidis et al. 2004 `_ model, :math:`r_{-2}/a=2+\gamma`) |br| ``NFW``: :math:`\rho \propto r^{-1}(r+a)^{-2}` (`Navarro, Frenk & White 1996 `_ model, :math:`r_{-2}/a=1`) |br| ``cNFW``: :math:`\rho \propto (r+a)^{-3}` (cored NFW model, :math:`r_{-2}/a=2`) |br| ``gNFW``: :math:`\rho \propto r^\gamma (r+a)^{-3-\gamma}` (generalized NFW model with free inner slope, :math:`r_{-2}/a=2+\gamma`) |br| ``Plummer``: :math:`\rho \propto (r^2+a^2)^{-5/2}` (`Plummer 1911 `_ model, :math:`r_{-2}/a=\sqrt{2/3}`) |br| ``gPlummer``: :math:`\rho \propto r^\gamma (r^2+a^2)^{-5/2-\gamma/2}` (generalized Plummer with free inner slope :math:`\gamma`, :math:`r_{-2}/a=\sqrt{(2+\gamma)/3}`) |br| ``-darknormflag`` *dark_norm_flag* dark matter normalization; available flags are [default: ``-1``]: ``-1``: :math:`\log_{10} (r_{\rm vir} \rm / kpc)` |br| ``0``: :math:`\log_{10} (M_{\rm vir} / M_\odot)` |br| > ``0``: :math:`r_{\rm fid}` (kpc) with norm = :math:`M(r_{\rm fid})` ``-darkscaleflag`` *dark_scale_flag* dark matter scale flag; available values are [default: ``1``]: ``1``: scale radius * :math:`R_{\rm eff}` (``PrugnielSimien``) * scale *a* (``isothermal`` or ``Jaffe``) * radius of density logarithmic slope :math:`-2`, :math:`r_{-2}` (other models, see ``-darkmodel``) |br| ``2``: concentration (virial over scale radius) ``-darktotflag`` *dark_tot_flag* dark vs. total flag; available values are [default: ``1``]: ``1``: dark |br| ``2``: total ``-v3dmodel`` *model* model for 3D distribution of velocities [default: ``Gaussian``] ``-rmax`` *rmax* max integration radius [kpc, default: ``40000``] **Options on variable parameters:** ``-lrtr`` *min* *max* [*min2* *max2* ... ] log_10 tracer radius (kpc); tracer radii are defined as: * scale radius *a* (``isothermal`` or ``Jaffe``) * effective radius :math:`R_{\rm eff}` (``PrugnielSimien``) * radius of logarithmic density slope :math:`-2` (other models, see ``-tracermodel``) ``-lMtr`` *min* *max* [*min2* *max2* ... ] log_10 tracer mass (:math:`\rm M_\odot`) [default: ``-99`` ``-99`` [...]]; tracer masses are defined as mass within tracer radius (see ``-lrtr``) ``-lMtrtot`` *min* *max* log_10 total tracer mass (:math:`\rm M_\odot`) [default: ``-99`` ``-99``]; tracer masses are defined as mass within tracer radius (see ``-lrtr``) ``-ftr`` *min* *max* [*min2* *max2* ...] tracer mass fraction in component (if total tracer mass is set with ``-lMtrtot``) [default: ``0`` ``0`` [...]] ``-tr2`` *min* *max* [*min2* *max2* ...] 2nd tracer parameter (inner slope or Prugniel Simien index) [default: ``0`` ``0`` [...]] ``-anisflag`` *anisotropy_flag* velocity anisotropy definition flag: ``0``: logarithmic: :math:`\log_{10} (\sigma_r/\sigma_\theta)` |br| ``1``: standard: :math:`\beta = 1-\sigma_\theta^2/\sigma_r^2` |br| ``2``: symetric: :math:`\displaystyle \beta_{\rm sym} = {\sigma_r^2-\sigma_\theta^2 \over \sigma_r^2+\sigma_\theta^2}` ``-anis0`` *min* *max* [*min2* *max2* ...] inner velocity anisotropies (definition in *anisotropy_flag*) [default: ``0`` ``0`` [...]] ``-anisinf`` *min* *max* [*min2* *max2* ...] outer velocity anisotropies (definition in *anisotropy_flag*) [default: ``0`` ``0`` [...]] ``-beta0`` *min* *max* [*min2* *max2* ...] same as ``-anis0`` but force ``-anisflag 1`` [default: ``0`` ``0`` [...]] ``-betainf`` *min* *max* [*min2* *max2* ...] same as ``-anisinf`` but force ``-anisflag 1`` [default: ``0`` ``0`` [...]] ``-lranis`` *min* *max* [*min2* *max2* ...] log anisotropy radius (kpc) [default: ``0`` ``0`` [...]]; can be circumvented with ``-TALflag`` ``-norm`` *min* *max* dark or total (see ``-darktotflag``) mass normalization (see ``-darknormflag``) ``-lrdark`` *min* *max* log_10 dark or total scale radius (kpc) OR log_10 dark or total concentration (depending on *dark_scale_flag*) ``-lc`` *min* *max* same as ``-lrdark`` *min* *max*, but forces ``-darkscaleflag 2`` ``-darkpar2`` *min* *max* dark matter 2nd parameter (inner slope or Prugniel Simien index) [default: ``1`` ``1``] ``-lMBH`` *min* *max* log_10 black hole mass (:math:`\rm M_\odot`) [default: ``-99`` ``-99``] ``-lB`` *min* *max* log_10 interloper densioty in projected phase space (*B* parameter in `Mamon, Biviano & Murante 2010 `_, in virial units :math:`N_{\rm vir} \, r_{\rm vir}^{-2}` *v* :math:`{\!}_{\rm vir}^{-1}`) [default: ``-99`` ``-99``] **Options on constraints**: ``-MfL`` Mass follows light [default: not set] ``-TLM`` Tied Light Mass [default: not set] ``-TAL`` Tied Anisotropy Light [default: not set] ``-split`` solve independently for :math:`r_{\rm tracer}` (possibly using external constraints provided by ``-meanlrnu`` and ``-siglrnu``) before fitting :math:`p(\!` *v* :math:`{\!}_{\rm LOS}|R)` (instead of joint fit of :math:`(R,\!` *v* :math:`{\!}_{\rm LOS})` [default: not set, unless *meanlr* :math:`> 0` from ``-meanlrtr`` AND *sigmalr* :math:`> 0` from ``-siglrtr``] ``-ilopflag`` *flag* interloper flag [default: ``1``]: ``0``: spherical halo only (no interlopers) |br| ``1``: standard to infinity |br| ``2``: spherical halo and interlopers in same cone treated separately |br| flags ``0`` and ``2`` are only implemented for NFW tracer models **Cosmological parameters**: ``-z`` *z* redshift of object [default: ``0.0``] ``-Omegam`` *Omegam* cosmological density parameter at :math:`z=0`, :math:`\Omega_{\rm m}^0` [default: ``0.3``] ``-h`` *h* Hubble constant at :math:`z=0` (units of 100 km/s/Mpc) [default: ``0.7``] ``-Delta`` *Delta* value of virial overdensity relative to critical density of Universe at observed redshift, 0 for `Bryan & Norman 1998 `_ [default: ``200``] ``-cofM`` *a0* *a1* coefficients for concentratio-nmass relation: :math:`\log_{10} c = a_0 + a_1 \log_{10} (M/{\rm M}_\odot)` [e.g. ``2.12`` ``-0.10``, default: not set] ``-cofMdef`` assume concentration-mass relation with default coefficients [``2.12`` ``-0.10``] ``-mu0`` *mu* distance modulus of system (if known independently of redshift) **MCMC options**: ``-samples`` *N* number *N* of elements per chain before ending MCMC [default: :math:`10^4` times the number of free ajustable parameters] ``-cv`` *conv* convergence test on multiple chains: ratio of variance of means to mean of variances of chains (e.g. ``1``, ``0.1``, or ``0.01``; default: not set) ``-p`` *num* propose scale [default: ``2.4``] ``-burn_in`` *N* wait for *N* \"burn in\" elements per chain before speeding up search using covariance matrix calculations [default: ``1000``] ``-feedback`` *f* CosmoMC verbosity [default: ``0``] ``-seed`` *integer* random seed: ``0`` for seed based on clock [default: ``1234``] **Other options**: ``-v`` *vers* use version *vers* (default: not set)] ``-debug`` *level* MAMPOSSt verbosity [default: ``0``] ``-help`` display this message and exit