def METHOD_NAME():
    model = SkyModel.create("pl", "gauss", name="test")
    model.spectral_model.amplitude.min = 1e-15
    model.spectral_model.amplitude.max = 1e-10
    estimator = FluxEstimator(source="test", selection_optional=[], reoptimize=True)
    scale_model = estimator.get_scale_model(Models([model]))
    assert_allclose(scale_model.norm.min, 1e-3)
    assert_allclose(scale_model.norm.max, 1e2)
    assert scale_model.norm.interp == "log"