Adding computation of HDP

@charlie.hoy as discussed, for production analyses it would be great if the computation of the HDP interval to check the level of exclusion of equal mass could be implemented directly in pesummary rather than be done in an additional post-processing step.

Example code from the GW170817 analysis can be found in this notebook

I am happy to help with the review.

Many thanks, Patricia