Adapt joint analysis tutorial
"\n",
" # optionally define a safe energy threshold\n",
" emin = None\n",
- " dataset.mask_safe = dataset.counts.geom.energy_mask(emin=emin)\n",
+ " data = dataset.counts.geom.energy_mask(emin=emin)\n",
+ " dataset.mask_safe = Map.from_geom(geom=dataset.counts.geom, data=data)\n",
" datasets.append(dataset)"
]
},
"\n",
"for dataset in datasets:\n",
" residuals = dataset.residuals()\n",
- " coords = residuals.geom.get_coord()\n",
- "\n",
- " residuals_stacked.fill_by_coord(coords, residuals.data)"
+ " residuals_stacked.stack(residuals)"
]
},
{
"metadata": {},
"outputs": [],
"source": [
- "residuals_stacked.sum_over_axes().smooth(\"0.1 deg\").plot(\n",
+ "residuals_stacked.sum_over_axes().smooth(\"0.08 deg\").plot(\n",
" vmin=-1, vmax=1, cmap=\"coolwarm\", add_cbar=True, stretch=\"linear\"\n",
");"
]
