Adapted from corner.py by Rodrigo Luger on 7/11/17 to produce 2D heatmaps instead of contour plots for the histograms of planet-planet occultations.