def METHOD_NAME(field, data):
    x_H = 1.0 - data[(ptype, "He_metallicity")] - data[(ptype, "metallicity")]
    return (
        x_H
        * data[(ptype, "density")]
        * data[(ptype, "NeutralHydrogenAbundance")]
    )