def METHOD_NAME(self):
    """Are proatom densities imported correctly?"""
    self.parse()
    self.analysis = Hirshfeld(self.data, self.volume, os.path.dirname(os.path.realpath(__file__)))
    refH_den = [
        2.66407645e-01,
        2.66407645e-01,
        2.66407643e-01,
        2.66407612e-01,
        2.66407322e-01,
    ]  # Hydrogen first five densities
    refH_r = [
        1.17745807e-07,
        4.05209491e-06,
        3.21078677e-05,
        1.39448474e-04,
        4.35643929e-04,
    ]  # Hydrogen first five radii
    refO_den = [
        2.98258510e02,
        2.98258510e02,
        2.98258509e02,
        2.98258487e02,
        2.98258290e02,
    ]  # Oxygen first five densities
    refO_r = [
        5.70916728e-09,
        1.97130512e-07,
        1.56506399e-06,
        6.80667366e-06,
        2.12872046e-05,
    ]  # Oxygen first five radii
    assert_allclose(self.analysis.proatom_density[0][0:5], refO_den, rtol=1e-3)
    assert_allclose(self.analysis.proatom_density[1][0:5], refH_den, rtol=1e-3)
    assert_allclose(self.analysis.proatom_density[2][0:5], refH_den, rtol=1e-3)