Updat evoid plots

notebooks/flow/void_test.ipynb

@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -286,18 +286,42 @@
    ]
   },
   {
-   "cell_type": "code",
-   "execution_count": null,
+   "cell_type": "markdown",
    "metadata": {},
-   "outputs": [],
-   "source": []
+   "source": [
+    "###  5. $V_{\\rm ext}$ along the model axis and $\\beta = 1$"
+   ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
    "metadata": {},
    "outputs": [],
-   "source": []
+   "source": [
+    "for simname in [\"IndranilVoid_exp\", \"IndranilVoid_gauss\", \"IndranilVoid_mb\"]:\n",
+    "    X = []\n",
+    "    for catalogue in [\"2MTF\", \"SFI_gals\", \"CF4_TFR_i\", \"CF4_TFR_w1\"]:\n",
+    "\n",
+    "        fname = paths.flow_validation(\n",
+    "            fdir, simname, catalogue, inference_method=\"bayes\",\n",
+    "            sample_alpha=False, sample_beta=None,\n",
+    "            no_Vext=True, zcmb_max=0.05, sample_Vmag_vax=True)\n",
+    "        \n",
+    "        X_i = samples_to_getdist(get_samples(fname, False), catalogue_to_pretty(catalogue))\n",
+    "        X.append(X_i)\n",
+    "\n",
+    "    params = [\"rLG\", \"sigma_v\", \"Vext_axis_mag\"]\n",
+    "    with plt.style.context(\"science\"):\n",
+    "        g = plots.get_subplot_plotter()\n",
+    "        g.settings.figure_legend_frame = False\n",
+    "        g.settings.alpha_filled_add = 0.75\n",
+    "\n",
+    "        g.triangle_plot(X, params=params, filled=True, legend_loc='upper right')\n",
+    "        plt.gcf().suptitle(simname_to_pretty(simname), y=1.025)\n",
+    "        plt.gcf().tight_layout()\n",
+    "        plt.gcf().show()\n",
+    "        plt.gcf().savefig(f\"../../plots/void_{simname}_Vext_along_axis_no_beta.png\", dpi=500, bbox_inches='tight')"
+   ]
   },
   {
    "cell_type": "code",