Introduction

Welcome to manticore! Manticore solves for the temperature and column density of molecular gas clouds given multi-wavelength infrared maps of the source. The data for each wavelength must be provided in a separate FITS file containing a single HDU (header-data unit), which must include at least the following keywords:

Keyword	Contents
BAND	band center (microns)
BUNIT	map data units (MJy/sr)
OBJECT	target object name
TELESCOP	facility name

These keywords should naturally be supplemented by coordinate system specifications CRVALn, CDELTn, etc., but these are merely copied to the output and not used during processing. Note also the following assumptions on input maps:

Unit conversion is not supported; in particular, the maps must be in units of MJy/sr.
All maps are assumed to be coincident on the sky and must have the same pixel dimensions.

Manticore requires detailed knowledge of the detector bandpass response to process data, encapsulated in the Detector class. The following detectors are supported:

Herschel PACS (70, 100 and 160 micron bands)
Herschel SPIRE (250, 350 and 500 micron bands)

Adding support for a new detector requires creation of a new concrete class and updating logic in solve.cc to associate it with an input map as appropriate.

Calculations support the use of arbitrary dust opacity models, as encapsulated by the Dust class. Tabular models must be compiled into the program by adding the table to Dust.cc and updating dust::modelMap accordingly. Dust models should be specified for a gas-to-dust ratio of 100.

Manticore Usage

Manticore is a simple command line utility. To see the full list of program options, run manticore -vh (or --verbose --help). At the time of writing, the output was:

Manticore version 1.2, released Aug 18 2018.

Usage:  manticore [OPTIONS] BAND1.fits BAND2.fits BAND3.fits [...]

Options:

  -2, --two-param
        Fit two parameter dust emission model (default).
        (Single-temperature cloud model.)

  -3, --three-param
        Fit three parameter dust emission model (default: two-param).
        (Two-temperature cloud model with prescribed halo temperature.)

  -4, --four-param
        Fit four parameter dust emission model (default: two-param).
        (Two-temperature cloud model, all parameters free.)

  -A, --accuracy=RELERR
        Maximum relative error for least-squares fit (default: 1e-2).

  -a, --abs-diff
        Output absolute differences in delta maps (default: sigma-relative).

  -B, --no-band-output
        Do not append input band data to output FITS file.

  -C, --check
        Generate comparative PACS color correction tables and exit.

  -c, --chi2-max=CHI2
        Maximum chi^2/DOF for meaningful 2-param result (default: 1.0).

  -D, --dust-model=NAME
        Dust opacity model (default: OH5).
        (Available models: OH5)

  -e, --error-map=FILE
        Input FITS error map (BAND and TELESCOP keys must match parent image).

  -E, --error-pct=PCT
        Percent flux error when error map unavailable (default: 3.0).

  -F, --force
        Force overwrite of existing files.

  -g, --gas-to-dust=RATIO
        Gas-to-dust ratio (default: 100.0).

  -h, --help
        Display command line help summary and exit.

  -l, --log-file=FILE
        Redirect log output to FILE (default is stderr).

  -M, --model-map
        Replace input data with model observation grid.

  -m, --max-iterations=NMAX
        Maximum least-squares iterations per pixel (default: 25).

  -n, --num-bands=NMIN
        Minimum number of bands required for least-squares fit (default: 3).

  -o, --output-file=FILE
        Output result to FITS file FILE (default: mtcore.fits).

  -p, --point-source
        Use point source detector response (default is extended sources).

  -q, --quiet
        Quiet mode; decrement logging level (opposite of -v).
        (Specify multiple times for added effect.)

  -s, --single-temp=FILE
        FITS output file from the single-temperature (2-param) fit.

  -T, --temp0=T0
        Initial guess for (halo) temperature in Kelvin (default: 16.5).

  -t, --threads=NTHREADS
        Use NTHREADS threads for processing (default: 1).

  -v, --verbose
        Verbose mode; increment logging level (opposite of -q).
        (Specify multiple times for added effect.)

  -V, --version
        Display program version and exit.

Compilation Requirements

GNU Science Library, available in the base repository on RHEL 6/7. Version 1.15 is known to work. Version 2.2 and later are not supported at this time as the non-linear fitting API was completely rewritten and is not compatible.
cfitsio, C libraries for manipulating FITS files. Available in the EPEL repository if not installed.
CCfits, a C++ wrapper library for the cfitsio library.
MathJax is required for the Doxygen documentation. If the MathJax/ directory (or link) is present, the generated Doxyfile will use it; otherwise, a web link will be used (hence requiring live web access to see equations in the documentation).

To compile, just type make and make docs. The manticore executable will be in the top-level directory, documentation in html/.