def METHOD_NAME(self):
    r"""Compare the SM $m_{\ell\ell}$ spectrum from flavio to madgraph prediction in 10 bins from 200 GeV to 2 TeV
        MG events simulated with NNPDF30_nnlo_as_0118 PDFs
        2% acceptable error"""
    spectrum_mg = np.asarray([2.66250518e-02, 8.76624628e-03, 2.92981721e-03, 9.64842852e-04, 3.12854454e-04, 9.77018348e-05, 2.84022912e-05, 8.08026977e-06, 2.02974822e-06])
    # Number of SM events generated in MG [60169., 25284., 10791., 25514., 52059., 20750., 42577., 45628., 14659.]
    s = 13e3**2
    logbins = np.linspace(np.log10(200),np.log10(1800),10)
    bins = 10**logbins
    nbins = len(bins)-1
    for i in range(nbins):
        center = 0.5*(bins[i]+bins[i+1])
        width = bins[i+1]-bins[i]
        spectrum = ppll.sigma_qqll_int(s, bins[i], bins[i+1], 'mu', center**2, 0, par2, 0, newphys=False)*GeVtopb/width
        err = (spectrum-spectrum_mg[i])/spectrum_mg[i]
        self.assertAlmostEqual(err,0,delta=0.02,msg=f'error in bin {i}: {err}')