def extract_yields_stats(yields):
    """Extract coverage, mean and std of yields."""
    coverage_mask = (yields != 0.0)
    return coverage_mask.mean(), yields[coverage_mask].mean(), yields[coverage_mask].std()