Adapt analysis 1 tutorial
"# We now fix the energy axis for the counts map\n",
"config.datasets.geom.axes.energy.min = \"1 TeV\"\n",
"config.datasets.geom.axes.energy.max = \"10 TeV\"\n",
- "config.datasets.geom.axes.energy.nbins = 4\n",
+ "config.datasets.geom.axes.energy.nbins = 10\n",
"\n",
"# We now fix the energy axis for the IRF maps (exposure, etc)\n",
"config.datasets.geom.axes.energy_true.min = \"0.5 TeV\"\n",
"config.datasets.geom.axes.energy_true.max = \"20 TeV\"\n",
- "config.datasets.geom.axes.energy.nbins = 10"
+ "config.datasets.geom.axes.energy_true.nbins = 20"
]
},
{
"source": [
"config.fit.fit_range.min = 1 * u.TeV\n",
"config.fit.fit_range.max = 10 * u.TeV\n",
- "config.flux_points.energy = {\"min\": \"1 TeV\", \"max\": \"10 TeV\", \"nbins\": 3}\n",
+ "config.flux_points.energy = {\"min\": \"1 TeV\", \"max\": \"10 TeV\", \"nbins\": 4}\n",
"config.excess_map.correlation_radius = 0.1 * u.deg"
]
},
"metadata": {},
"outputs": [],
"source": [
+ "%%time\n",
"analysis.run_fit()"
]
},
"outputs": [],
"source": [
"analysis.get_excess_map()\n",
- "analysis.excess_map[\"sqrt_ts\"].plot(add_cbar=True);"
+ "analysis.excess_map[\"sqrt_ts\"].plot(\n",
+ " add_cbar=True, cmap=\"RdBu\", vmin=-5, vmax=5\n",
+ ");"
]
},
{
],
"metadata": {
"kernelspec": {
- "display_name": "Python 3",
+ "display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
- "version": "3.7.7"
+ "version": "3.9.12"
},
"widgets": {
"application/vnd.jupyter.widget-state+json": {
