7.1. Galactic (Milky Way): \(\tt -g\)

Run for a Galactic halo + list of haloes (optional)

Note

All default parameter files from the following run examples are generated with the -D flag.

$ clumpy -gX -D

7.1.1. \(\tt -g0\): mass profile (1D)

Calculates \(M_{\rm Gal.}(r)\) in one dimension:

$ clumpy -g0 -i clumpy_params_g0.txt
_images/g0D.png

Fig. 7.1 Mass profile of the Galactic DM halo.

Note that a triaxial halo can be studied as well.


7.1.2. \(\tt -g1\): density profiles (1D)

Calculates \(\rho_{\rm Gal.}(r)\) (smooth, sub-continuum and total) of a spherically symmetric halo in one dimension:

$ clumpy -g1 -i clumpy_params_g1.txt
_images/g1D.png

Fig. 7.2 Density profile of the Galactic DM halo. This is Fig. 5.1 already created in the Quick start tutorial.


7.1.3. \(\tt -g2\): J-factor profiles (1D)

Calculates \(J(\alpha_{\rm int};\,l,b)\) towards a user-defined direction \((l,b)\) in the Galactic halo. Note that a triaxial halo can be studied as well.

$ clumpy -g2 -i clumpy_params_g2.txt
_images/g2D-1.png

Fig. 7.3 J-factor (of smooth, sub-continuum and total components) towards the Galactic anti-center as a function of the integration radius \(\alpha_{\rm int}\) of the search cone \(\Delta\Omega\).

_images/g2D-2.png

Fig. 7.4 Flux (here: integrated in energy) towards the Galactic anti-center as a function of the integration radius \(\alpha_{\rm int}\).


7.1.4. \(\tt -g3\): \({\rm d}J/{\rm d}\Omega\) (1D)

Calculates \({\rm d}J/{\rm d}\Omega(\theta)\) as a function of the distance \(\theta\) from the centre of a spherically symmetric Galactic halo.

$ clumpy -g3 -i clumpy_params_g3.txt
_images/g3D-1.png

Fig. 7.5 \({\rm d}J/{\rm d}\Omega\) (of smooth, sub-continuum and total components) over the distance \(\theta\) from the Galactic centre.

Like in the -g2 module, also the corresponding \(\gamma\)-ray and neutrino intensity can be calculated if enabled in the parameter file:

_images/g3D-3.png

Fig. 7.6 \({\rm d}\Phi/{\rm d}\Omega\) in \(\gamma\)-rays and neutrinos over the distance \(\theta\) from the Galactic centre.


7.1.5. \(\tt -g4\): \(J_{\alpha_{\rm int}}(\theta)\) (1D)

Returns the \(J(\theta)\) and optionally, the corresponding fluxes, for a fixed, user-defined integration angle \(\alpha_{\rm int}\) as a function of the distance \(\theta\) from the Galactic centre. Any J-factors \(J_{\rm list}(\alpha_{\rm int})\) from haloes defined in the file gLIST_HALOES are overplotted as well:

$ clumpy -g4 -i clumpy_params_g4.txt
_images/g4D-1.png

Fig. 7.7 Galactic \(J(\theta;\,\alpha_{\rm int})\) (of smooth, sub-continuum and total components) over the distance \(\theta\) from the Galactic centre.

_images/g4D-2.png

Fig. 7.8 Corresponding Galactic DM fluxes \(\Phi\) in \(\gamma\)-rays and neutrinos over the distance \(\theta\) from the Galactic centre.


7.1.6. \(\tt -g5\): \(\langle J_{\rm pixel}\rangle\) and \(\langle {\rm d}J/{\rm d}\Omega\rangle\) (2D)

$ clumpy -g5 -i clumpy_params_g5.txt
_images/g5D.png

Fig. 7.9 \(4^{\circ}\times 4^{\circ}\) skymap of \(J(\rm pixel)\) towards the Galactic anti-center without resolving individual Galactic subhaloes.

Corresponding flux- and intensity maps for all -g5 to -g8 are calculated with the gSIM_IS_WRITE_FLUXMAPS = True enabled.


7.1.7. \(\tt -g6\): g5 + list haloes (2D)

Same as -g5, but allows to draw individual resolved haloes from a list defined in gLIST_HALOES into the skymap.

$ clumpy -g6 -i clumpy_params_g6.txt
_images/g6D.png

Fig. 7.10 \(90^{\circ}\times 45^{\circ}\) skymap of \(J(\rm pixel)\) in cartesian (left) and Hammer-Aitoff (right) projection towards the Galactic centre without resolving Galactic subhaloes, but including a user-defined halo at \((l,b) =(20^{\circ},20^{\circ})\). Note the user-drawn halo at \((l,b) =(20^{\circ},20^{\circ})\) faintly visible in the right panel.

Note

Only Galactic objects with redshift \(z\ll 1\) can be added in these maps in the current release.


7.1.8. \(\tt -g7\): \(J_{\rm pixel}\) and \({\rm d}J/{\rm d}\Omega\) + drawn subhaloes (2D)

Same as -g5, but allows to draw individual randomly drawn resolved Galactic subhaloes from a list defined in gLIST_HALOES into the skymap.

$ clumpy -g7 -i clumpy_params_g7.txt
_images/g7D.png

Fig. 7.11 \(4^{\circ}\times 4^{\circ}\) skymap of \(J(\rm pixel)\) in cartesian projection towards the Galactic anti-center with resolved and unresolved Galactic subhaloes. The left panel shows the components (emission from smooth Galactic halo, unresolved halos, cross-product, and resolved halos), the right panel the sum.

In this option run, for the selected gSIM_USER_RSE threshold, 54 subhaloes are resolved. If more than 20 subhaloes are resolved in a run, the following interactive ROOT population study plots are generated:

_images/g7D-popstudy.png

Fig. 7.12 Plots to study the proerties of the resolved subhalo population in the default -g7 run.


7.1.9. \(\tt -g8\): g7 + list haloes (2D)

The -g8 module does the same as -g7, however, allows to additionally draw user-defined haloes from a list, stored in a file called by the gLIST_HALOES parameter:

$ clumpy -g8 -i clumpy_params_g8.txt
_images/g8D.png

Fig. 7.13 \(4^{\circ}\times 4^{\circ}\) skymap of \(J(\rm pixel)\) in cartesian projection towards the Galactic anti-centre with resolved and unresolved Galactic subhaloes plus the example halo rs10_gamma10 read from the definition file $CLUMPY/data/list_generic.txt.

Note

You need the CLUMPY environmental variable pointing to your installation directory in order to find list_generic.txt.