diff --git a/examples/cherenkov_muon.png b/examples/cherenkov_muon.png
index 4903664b57e87a5e3070f6f22d56ef7ed4208614..ba04e55084a7079bae64845b4d8fe4c8f2f43a96 100644
Binary files a/examples/cherenkov_muon.png and b/examples/cherenkov_muon.png differ
diff --git a/examples/cherenkov_muon.py b/examples/cherenkov_muon.py
index d323cc3a698c9ccec6da02562fc95ada70e4317c..7c1637825e7c7536184459336335309728888226 100755
--- a/examples/cherenkov_muon.py
+++ b/examples/cherenkov_muon.py
@@ -11,7 +11,7 @@ Usage:
 
 Options:
     -f FILE_DESCRIPTOR  Descriptor of the PDFs (e.g. $JPP_DIR/data/J%p.dat).
-    --energy=<GeV>      Muon energy [default: 1.0].
+    --energy=<GeV>      Muon energy [default: 1e3].
     --distance=<m>      Distance of approach [default: 50].
     --tts=<ns>          PMT time smearing [default: 0.0].
     --xmin=<ns>         X-axis minimum [default: -20].
@@ -62,6 +62,9 @@ def main():
         ax.plot(xs, list(map(pdf, xs)), label=direction)
         ax.set_yscale("log")
 
+    plt.title(
+        "PDF of Cherenkov light from muon with E = {} GeV and R = {} m".format(
+            args['--energy'], args['--distance']))
     plt.ylim([float(args['--ymin']), float(args['--ymax'])])
     plt.ylabel(r"$\frac{dP}{dt}$ / (npe/ns)")
     plt.xlabel(r"$\Delta$t / ns")