def METHOD_NAME():
    energy = MapAxis.from_energy_bounds(0.1, 10, 3.0, unit="TeV", name="energy_true")
    template = WcsNDMap.create(npix=10, axes=[energy], unit="cm-2 s-1 sr-1 TeV-1")
    spatial = TemplateSpatialModel(template, normalize=False)
    spectral = PowerLawNormSpectralModel()
    model = SkyModel(spectral_model=spectral, spatial_model=spatial, name="test")
    estimator = FluxEstimator(source="test", selection_optional=[], reoptimize=True)
    with pytest.raises(ValueError, match="`NormSpectralModel` are not supported"):
        estimator.get_scale_model(Models([model]))