Forgive me if that is a silly question, but I was not able to dig its answer in a simple way from the documentation and I am not quite conversant with Python...

I have a data file with 3 columns of best fit parameters (say, [math]) from a simulation I ran in a completely independent way from COSMOMC (or CAMB). Thus my data file does not have the usual format arising from COSMOMC output. Nevertheless, I would still like, if possible, to use GetDist (via its python wrapper) to:

(1) create a kernel density estimate, from the corresponding scatterplot of $(a, b, c)$ triples, for the (normalized) probability density function [math].

(2) from item (1), determine the mode (maximum) of the estimated pdf

(3) from item (1), determine the contour lines for, e.g., the 68%, 95% and 99% highest confidence/credibility regions in the distinct 2d parameter planes

(4) still from item (1), determine the 1d marginal pdf's [math] and [math]

Finally, I would like to plot all this information in a single "triangle_plot" (à la https://getdist.readthedocs.io/en/latest/plot_gallery.html) but,

**additionally**, with:

- 2d scatterplots of the data
- 2d modes

Perhaps all this is achievable, if I manage to get the identity of each of the subplots constituting the triangle plot.

Thanks in advance!