Problem with exclusion regions

[matteobianucci]

Hello,

I have discovered an issue concerning the exclusion region usage with baccmod.

Here you can find two images, both are computed background models for the same observation, one is computed with an exclusion region of 0.4 degrees (top figure) and the second is computed without the exclusion region (down figure)

Here is the parameters used when computing the model :

nbins = 5.109457716

offset_axis_acc = MapAxis.from_bounds(0.0 * u.deg, 3 * u.deg, nbin=22, name='offset')
energy_axis_acc = MapAxis.from_energy_bounds(0.03 , 2, nbin=nbins*10 , per_decade=True, name='energy', unit = 'TeV')

exclude_regions = [CircleSkyRegion(center=exclusion_pos,radius=0.4 * u.deg), ]

acceptance_model_creator = Grid3DAcceptanceMapCreator(energy_axis_acc,offset_axis_acc,dynamic_energy_axis=True,
dynamic_energy_axis_target_statistics=500,
method='stack',
cos_zenith_binning_parameter_value = 3600,
initial_cos_zenith_binning = 0.01, zenith_binning_run_splitting = True,
use_mini_irf_computation = True,
mini_irf_time_resolution = 1. * u.min)
exclude_regions=exclude_regions)

acceptance_model = acceptance_model_creator.create_acceptance_map_per_observation(obs_collection,zenith_binning=True, runwise_normalisation=False)