def METHOD_NAME(self):
    """Check tetragonal lattice."""
    a = 2
    b = 4
    c = 8
    xen = range(0, self.box_l, a)
    yen = range(0, self.box_l, b)
    zen = range(0, self.box_l, c)
    for i, j, k in itertools.product(xen, yen, zen):
        self.system.part.add(type=self.part_ty, pos=(i, j, k))
    wavevectors, intensities = self.system.analysis.structure_factor(
        sf_types=[self.part_ty], sf_order=self.sf_order)
    intensities = np.around(intensities, 8)
    # no reflection conditions on (h,k,l)
    peaks_ref = self.generate_peaks(a, b, c, lambda h, k, l: True)
    peaks = self.peak_orders(wavevectors[np.nonzero(intensities)])
    np.testing.assert_array_equal(peaks, peaks_ref[:len(peaks)])