def METHOD_NAME():
    import naima
    ECPL = naima.models.ExponentialCutoffPowerLaw(
        1e36 * u.Unit("1/eV"), 1 * u.TeV, 2.1, 13 * u.TeV
    )
    IC = naima.models.InverseCompton(ECPL, seed_photon_fields=["CMB"])
    naima_model = NaimaSpectralModel(IC)
    model = SkyModel(spectral_model=naima_model, name="test")
    estimator = FluxEstimator(source="test", selection_optional=[], reoptimize=True)
    scale_model = estimator.get_scale_model(Models([model]))
    assert_allclose(scale_model.norm.min, np.nan)
    assert_allclose(scale_model.norm.max, np.nan)