def METHOD_NAME(self):
    wb = Water3SiteBox(
        box=mb.Box([2.0, 2.0, 2.0]), model=water_models.WaterSPC()
    )
    assert wb.children[0].name == "WaterSPC"
    parts = [p for p in wb.children[0].particles()]
    oh_dist1 = np.linalg.norm(parts[0].pos - parts[1].pos)
    oh_dist2 = np.linalg.norm(parts[0].pos - parts[2].pos)
    assert np.isclose(oh_dist1, 0.1)
    assert np.isclose(oh_dist2, 0.1)
    v1 = np.array(parts[1].pos) - np.array(parts[0].pos)
    v2 = np.array(parts[2].pos) - np.array(parts[0].pos)
    v1_u = v1 / np.linalg.norm(v1)
    v2_u = v2 / np.linalg.norm(v2)
    angle = (
        np.arccos(np.clip(np.dot(v1_u, v2_u), -1.0, 1.0)) * 180.0 / np.pi
    )
    assert np.isclose(angle, 109.47)