Add howto entry
@@ -86,6 +86,22 @@ extraction can then be performed as usual with the
`~gammapy.estimators.LightCurveEstimator`. This is demonstrated in the `Light curve -
Flare <tutorials/time/light_curve_flare.html>`__ tutorial.
+Choose units for plotting
++++++++++++++++++++++++++
+Units for plotting are handled with a combination of `matplotlib` and `astropy.units`.
+For most plotting methods Gammapy forwards additional keywords to the corresponding
+matplotlib plot method, including the `xunits` and `yunits` keywords, which allows
+you to define the x and y axis units using `astropy.units`. Here is a minimal example:
+
+.. code::
+
+ from gammapy.estimators import FluxPoints
+ from astropy import units as u
+
+ filename = "$GAMMAPY_DATA/hawc_crab/HAWC19_flux_points.fits"
+ fp = FluxPoints.read(filename)
+ fp.plot(sed_type="e2dnde", xunits=u.erg, yunits=u.Unit("erg cm-2 s-1"))
+
Compute source significance
+++++++++++++++++++++++++++
