Display mean and deviation values on grouped boxplot in Python

Question

I want to display mean and standard deviation values above each of the boxplots in the grouped boxplot (see picture).

My code is

import pandas as pd
im