After the comments from the referee the algorithm for AgeWizard has been refined somewhat:
The individual star age probability distribution functions are now such that the probability that a star has an age older than the current age of the Universe is 0
The total probability density function is now calculated from the individual star PDFs, not from their unnormalised distributions. The original algorithm effectively leads to a higher weighting of stars found in densely populated regions of the HRD and CMDs. As a result, stars in the giant branch carried a significantly lower weight than main sequence stars. This is problematic because these stars carry the most reliable information about age since their evolutionary phase is much shorter lived than the Main Sequence. Consequently, normalising the distributions over time before combining them is a first step towards making the age determination more accurate. Further steps can be taken by the human and these will be described in soon-to-come tutorials and publications. Feel free to contact hfstevance@gmail.com if you need more detail about this.
After the comments from the referee the algorithm for AgeWizard has been refined somewhat: