def METHOD_NAME(pig,Lbox, hh):
    """Get a conventionally unitted halo mass function in a resolved region.
    Lbox is box size in Mpc (not Mpc/h!)"""
    #Change units to M_sun
    fofmasses = pig['FOFGroups/Mass'][:]*10**10/hh
    #Find minimum halo mass
    rsl = 2*min(fofmasses[fofmasses>0])
    #Find the mass function
    smf = massfunc(fofmasses[fofmasses>rsl],Lbox)
    return smf