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)
1: weight by Gaussian(mu)
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: $\rho \propto r^{-1}(r+a)^{-3}$ (Hernquist 1990 model, $r_{-2}/a=1/2$)
mHubble: $\rho \propto (r^2+a^2)^{-3/2}$ (modified Hubble or non-truncated analytical King 1962 model, $r_{-2}/a=\sqrt{2}$)
isothermal: $\rho \propto (r^2+a^2)^{-1}$ (pseudo-isothermal)
Jaffe: $\rho \propto r^{-2}(r+a)^{-2}$ (Jaffe 1983)
NFW: $\rho \propto r^{-1}(r+a)^{-2}$ (Navarro, Frenk & White 1996 model, $r_{-2}/a=1$)
Plummer: $\rho \propto (r^2+a^2)^{-5/2}$ (Plummer 1911 model, $r_{-2}/a=\sqrt{2/3}$)
gPlummer: $\rho \propto r^\gamma (r^2+a^2)^{-5/2-\gamma/2}$ (generalized Plummer with free inner slope $\gamma$, $r_{-2}/a=\sqrt{(2+\gamma)/3}$)
PrugnielSimien: $\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 $b(n)$ is from analytical approximation of Ciotti & Bertin 1999, while $p(n)$ is given by Lima Neto, Gerbal & Marquez 1999)

-anismodel model [model2 … modeln]

velocity anisotropy model(s) model1 etc. of visible tracer(s); available models are:: iso: $\beta = 0$ (isotropic)
cst: $\beta$ = constant
ML: $\displaystyle \beta = {1\over 2} {r\over r+r_\beta}$ (Mamon & Łokas 2005b)
OM: $\displaystyle \beta = {r^2\over r^2+r_\beta^2}$ (Osipkov 1979; Merritt 1985)
gOM: $\displaystyle \beta = \beta_0 + (\beta_\infty-\beta_0)\,{r^2\over r^2+r_\beta^2}$ (generalized Osipkov-Merritt)

Tiret: $\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 $R_{\rm eff}$ (PrugnielSimien)
radius of logarithmic density slope $-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 $|\!$ v ${}_{\!\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: $\rho \propto (r+a)^{-1}(r^2+a^2)^{-1}$ (Burkert 1995 model, $r_{-2}/a=(1-\sqrt{26/27})^{1/3}+(1+\sqrt{26/27})^{1/3})$
Einasto: $\rho \propto \exp[-b(n)\, r^{1/n}]$ (Einasto 1965 model, $r_{-2}/a = (2n)^n$)
Hernquist: $\rho \propto r^{-1}(r+a)^{-3}$ (Hernquist 1990 model, $r_{-2}/a=1/2$)
gHernquist: $\rho \propto r^\gamma (r+a)^{-4-\gamma}$ (generalized Hernquist model with free inner slope $\gamma$, $r_{-2}/a=1+\gamma/2$)
mHubble: $\rho \propto (r^2+a^2)^{-3/2}$ (modified Hubble or non-truncated analytical King 1962 model, $r_{-2}/a=\sqrt{2}$)
isothermal: $\rho \propto (r^2+a^2)^{-1}$ (pseudo-isothermal)
Jaffe: $\rho \propto r^{-2}(r+a)^{-2}$ (Jaffe 1983)
Kazantzidis: $\rho \propto r^\gamma \exp(-r/a)$ (Kazantzidis et al. 2004 model, $r_{-2}/a=2+\gamma$)
NFW: $\rho \propto r^{-1}(r+a)^{-2}$ (Navarro, Frenk & White 1996 model, $r_{-2}/a=1$)
cNFW: $\rho \propto (r+a)^{-3}$ (cored NFW model, $r_{-2}/a=2$)
gNFW: $\rho \propto r^\gamma (r+a)^{-3-\gamma}$ (generalized NFW model with free inner slope, $r_{-2}/a=2+\gamma$)
Plummer: $\rho \propto (r^2+a^2)^{-5/2}$ (Plummer 1911 model, $r_{-2}/a=\sqrt{2/3}$)
gPlummer: $\rho \propto r^\gamma (r^2+a^2)^{-5/2-\gamma/2}$ (generalized Plummer with free inner slope $\gamma$, $r_{-2}/a=\sqrt{(2+\gamma)/3}$)

-darknormflag dark_norm_flag

dark matter normalization; available flags are [default: -1]:: -1: $\log_{10} (r_{\rm vir} \rm / kpc)$
0: $\log_{10} (M_{\rm vir} / M_\odot)$
> 0: $r_{\rm fid}$ (kpc) with norm = $M(r_{\rm fid})$

-darkscaleflag dark_scale_flag

dark matter scale flag; available values are [default: 1]:

1: scale radius

$R_{\rm eff}$ (PrugnielSimien)
scale a (isothermal or Jaffe)
radius of density logarithmic slope $-2$, $r_{-2}$ (other models, see -darkmodel)

2: concentration (virial over scale radius)

-darktotflag dark_tot_flag

dark vs. total flag; available values are [default: 1]:: 1: dark
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 $R_{\rm eff}$ (PrugnielSimien)
radius of logarithmic density slope $-2$ (other models, see -tracermodel)

-lMtr min max [min2 max2 … ]

log_10 tracer mass ($\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 ($\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: $\log_{10} (\sigma_r/\sigma_\theta)$
1: standard: $\beta = 1-\sigma_\theta^2/\sigma_r^2$
2: symetric: $\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 ($\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 $N_{\rm vir} \, r_{\rm vir}^{-2}$ v ${\!}_{\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 $r_{\rm tracer}$ (possibly using external constraints provided by -meanlrnu and -siglrnu) before fitting $p(\!$ v ${\!}_{\rm LOS}|R)$ (instead of joint fit of $(R,\!$ v ${\!}_{\rm LOS})$ [default: not set, unless meanlr $> 0$ from -meanlrtr AND sigmalr $> 0$ from -siglrtr]

-ilopflag flag

interloper flag [default: 1]:: 0: spherical halo only (no interlopers)
1: standard to infinity
2: spherical halo and interlopers in same cone treated separately

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 $z=0$, $\Omega_{\rm m}^0$ [default: 0.3]
-h h: Hubble constant at $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: $\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: $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

Usage of buildini_mamposst¶

Related Topics

This Page