def METHOD_NAME(self):
    """Check freud against data generated by MSM reference code.
    https://github.com/ArvDee/Minkowski-structure-metrics-calculator
    """
    rvd_data = _get_reference_data("RvD_MSM_avq4avq6avw4avw6")
    gsd_frame = _read_gsd_snapshot("Test_Configuration.gsd")
    freud_data = _compute_comparison_data(gsd_frame, average=True, weighted=True)
    names = ["ave. q'4", "ave. q'6", "ave. w'4", "ave. w'6"]
    for i, name in enumerate(names):
        # If most of the particles agree, this test is probably correct.
        # Close investigation showed that numerical issues (truncated
        # particle positions?) and bonds with facet weights of the order of
        # 1e-9 caused some deviations in the neighbor counts, which in turn
        # affected the averaged ql values. This affected 37 out of 3288
        # particles in the sample data for one of the cases checked.
        close_values = np.isclose(rvd_data[:, i], freud_data[:, i], atol=5e-5)
        assert np.sum(close_values) > (0.985 * len(freud_data)), f"{name} failed"