If you could do a quick review of this that would make me really happy.

The code makes three plots with information about the depth the data goes to and the number of fakes recovered.

Completeness 2D

A hexbin plot showing the faintest fake source recovered in each hex.

CompletenessHist

Histograms showing the number and fraction of fake sources inserted/recovered be magnitude bins.

CompletenessAreaDepth

A cumulative histogram showing the area covered to a given depth, where the depth is defined as being the faintest source recovered. This plot is made from the hex bin plot as a 1D summary.

I left some mostly minor code comments, except for the one about np.where(...)[0] since that might actually make a difference.

Completeness 2D:

• It's hard to see much with such a wide magnitude stretch; is this adjustable?
• Can you plot the n% completeness magnitude instead of faintest recovered?