def METHOD_NAME():
    # Tests the exponential cutoff
    ks = np.geomspace(1E-2, 15., 128)
    t = ccl.nl_pt.PTNumberCountsTracer(b1=1.0)
    ptc1 = ccl.nl_pt.LagrangianPTCalculator(cosmo=COSMO)
    pk2d1 = ptc1.get_biased_pk2d(t, tracer2=t)
    pk1 = pk2d1(ks, 1.0, cosmo=COSMO)
    ptc2 = ccl.nl_pt.LagrangianPTCalculator(k_cutoff=10.,
                                            n_exp_cutoff=2.,
                                            cosmo=COSMO)
    pk2d2 = ptc2.get_biased_pk2d(t, tracer2=t)
    pk2 = pk2d2(ks, 1.0, cosmo=COSMO)
    expcut = np.exp(-(ks/10.)**2)
    assert np.allclose(pk1*expcut, pk2, atol=0, rtol=1E-3)