Skip to content

Output Format and Options

Jump to bottom
Pearl Sandick edited this page Oct 4, 2019 · 1 revision

Output Formatting

Output may print either to the screen or to a file in the Output directory. If this directory does not exist, it will be created with the first run of MADHAT. Output will include input parameters for the run, and some or all of the following quantities (depending on the Run Option and input provided):

  • PhiPP [cm^3 s^-1 GeV^-2] - "particle physics factor"
    • prints if J-factors are provided for dwarfs
  • +/-dPhiPP [cm^3 s^-1 GeV^-2] - uncertainty due to varying J-factor
    • prints if J-factor uncertainties are provided for dwarfs
  • sigv [cm^3 s^-1] - annihilation cross section
    • prints if dark matter model parameters (mass and integrated spectrum) are provided
  • +/-dsigv [cm^3 s^-1] - uncertainty due to varying J-factor
    • prints if dark matter model parameters and J-factor uncertainties are provided

For detailed parameter definitions, see https://arxiv.org/1802.03826. Note: The cross section limit calculation assumes that the dark matter particles are their own antiparticles. If the dark matter is Dirac rather than Majorana, multiply the cross section by two. If the particle has some asymmetry between the particle and antiparticle, then the prefactor needs to be calculated.

Output Options

There are two optional output flags located in the file flags.dat, both of which default to "false".

  1. bool "fileout", on line 9 of flags.dat: If this is set to "true" then all run options print to a file [ ]_[beta].out in the Output directory rather than to the screen.

  2. bool "printbeta", on line 10 of flags.dat: If this is set to "true" then for running with two or zero arguments (run options B and C), the output will be a table of Nbound(beta) for Nbound increasing as an integer up to the value [beta] passed as input. Note that this output option is not defined for running with three inputs (Run Option A).

    • Each line of the output will either have the format "beta, Nbound(beta), PhiPP, +dPhiPP, -dPhiPP, sigv, +dsigv, -dsigv" or "beta, Nbound(beta), PhiPP, +dPhiPP, -dPhiPP."
    • To output all values of (beta, Nbound), use [beta]=0.999 as input.

Clone this wiki locally