fix savefile argument in scalar_total_xsec and fix u = 0 argument for example 2

Author: mhostert

examples/Example_1_total_xsecs.ipynb

@@ -100,14 +100,14 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "maj_test   = dn.model.ThreePortalModel(name='maj_test', m4 = 0.1, Umu4=np.sqrt(2.2e-7), epsilon2=4.6e-4, alphaD=0.4, UD4=1, HNLtype='majorana', mzprime=1.25)\n",
-    "dir_test   = dn.model.ThreePortalModel(name='dirac_test', m4 = 0.1, Umu4=np.sqrt(2.2e-7), epsilon2=4.6e-4, alphaD=0.4, UD4=1, HNLtype='dirac', mzprime=1.25)\n",
+    "maj_test   = dn.model.ThreePortalModel(name='maj_test', m4 = 0.2, Umu4=np.sqrt(2.2e-7), epsilon2=4.6e-4, alphaD=0.4, UD4=1, HNLtype='majorana', mzprime=1.25)\n",
+    "dir_test   = dn.model.ThreePortalModel(name='dirac_test', m4 = 0.2, Umu4=np.sqrt(2.2e-7), epsilon2=4.6e-4, alphaD=0.4, UD4=1, HNLtype='dirac', mzprime=1.25)\n",
     "\n",
     "calc_maj_coh   = dn.UpscatteringProcess(TheoryModel = maj_test, scattering_regime = 'coherent', **common_kwargs) \n",
     "calc_dir_coh   = dn.UpscatteringProcess(TheoryModel = dir_test, scattering_regime = 'coherent', **common_kwargs) \n",
     "\n",
-    "print(f'maj: tot_sigma = {calc_maj_coh.total_xsec(1.0, NEVAL=1000, NINT=10)}')\n",
-    "print(f'dir: tot_sigma = {calc_dir_coh.total_xsec(1.0, NEVAL=1000, NINT=10)}')\n"
+    "print(f'maj: tot_sigma = {calc_maj_coh.total_xsec(Enu=1.0, NEVAL=1000, NINT=10)}')\n",
+    "print(f'dir: tot_sigma = {calc_dir_coh.total_xsec(Enu=1.0, NEVAL=1000, NINT=10)}')\n"
    ]
   },
   {
@@ -121,8 +121,8 @@
     "\n",
     "fig, ax = dn.plot_tools.std_fig(figsize=(6,3))\n",
     "\n",
-    "ax.plot(q2,calc_maj_coh.diff_xsec_Q2(1,q2),label='majorana', lw=2)\n",
-    "ax.plot(q2,calc_dir_coh.diff_xsec_Q2(1,q2),ls='--',label='dirac', lw=2)\n",
+    "ax.plot(q2,calc_maj_coh.diff_xsec_Q2(1,q2, diagrams=['total']),label='majorana', lw=2)\n",
+    "ax.plot(q2,calc_dir_coh.diff_xsec_Q2(1,q2, diagrams=['total']),ls='--',label='dirac', lw=2)\n",
     "ax.set_yscale(\"log\")\n",
     "ax.set_xscale(\"log\")\n",
     "ax.set_xlabel(r\"$Q^2$/GeV$^2$\")\n",