From 7d1eaa7b729207ee73c7d53627f95752851a7047 Mon Sep 17 00:00:00 2001
From: cetin <cetin@physi.uni-heidelberg.de>
Date: Mon, 2 Oct 2023 16:55:09 +0200
Subject: [PATCH] rework mon

---
 B_rework.ipynb | 93 +++++++++++++++-----------------------------------
 1 file changed, 27 insertions(+), 66 deletions(-)

diff --git a/B_rework.ipynb b/B_rework.ipynb
index fe31740..29a4921 100644
--- a/B_rework.ipynb
+++ b/B_rework.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 21,
+   "execution_count": 1,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -22,7 +22,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 22,
+   "execution_count": 2,
    "metadata": {},
    "outputs": [
     {
@@ -31,7 +31,7 @@
        "9056"
       ]
      },
-     "execution_count": 22,
+     "execution_count": 2,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -50,7 +50,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 23,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [
     {
@@ -59,7 +59,7 @@
        "0.8606728758791105"
       ]
      },
-     "execution_count": 23,
+     "execution_count": 3,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -75,7 +75,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 24,
+   "execution_count": 19,
    "metadata": {},
    "outputs": [
     {
@@ -85,58 +85,19 @@
       "sample size:  32\n",
       "eff (cutoff =  0 ) =  0.96875\n",
       "sample size:  32\n",
-      "eff (cutoff =  50 ) =  0.96875\n",
-      "sample size:  32\n",
       "eff (cutoff =  100 ) =  0.96875\n",
-      "sample size:  43\n",
-      "eff (cutoff =  150 ) =  0.9767441860465116\n",
       "sample size:  65\n",
       "eff (cutoff =  200 ) =  0.9692307692307692\n",
-      "sample size:  97\n",
-      "eff (cutoff =  250 ) =  0.9587628865979382\n",
       "sample size:  129\n",
       "eff (cutoff =  300 ) =  0.9457364341085271\n",
-      "sample size:  150\n",
-      "eff (cutoff =  350 ) =  0.9533333333333334\n",
       "sample size:  169\n",
       "eff (cutoff =  400 ) =  0.9408284023668639\n",
-      "sample size:  197\n",
-      "eff (cutoff =  450 ) =  0.9390862944162437\n",
       "sample size:  227\n",
       "eff (cutoff =  500 ) =  0.920704845814978\n",
-      "sample size:  257\n",
-      "eff (cutoff =  550 ) =  0.9260700389105059\n",
-      "sample size:  297\n",
-      "eff (cutoff =  600 ) =  0.9326599326599326\n",
-      "sample size:  334\n",
-      "eff (cutoff =  650 ) =  0.9281437125748503\n",
-      "sample size:  366\n",
-      "eff (cutoff =  700 ) =  0.9289617486338798\n",
-      "sample size:  400\n",
-      "eff (cutoff =  750 ) =  0.925\n",
-      "sample size:  436\n",
-      "eff (cutoff =  800 ) =  0.9151376146788991\n",
-      "sample size:  468\n",
-      "eff (cutoff =  850 ) =  0.9102564102564102\n",
-      "sample size:  500\n",
-      "eff (cutoff =  900 ) =  0.912\n",
-      "sample size:  533\n",
-      "eff (cutoff =  950 ) =  0.9136960600375235\n",
-      "sample size:  562\n",
-      "eff (cutoff =  1000 ) =  0.9163701067615658\n",
-      "\n",
-      "sample size:  150\n"
+      "cutoff energy = 350MeV\n",
+      "sample size: 150\n",
+      "eff =  0.9533333333333334\n"
      ]
-    },
-    {
-     "data": {
-      "text/plain": [
-       "0.9533333333333334"
-      ]
-     },
-     "execution_count": 24,
-     "metadata": {},
-     "output_type": "execute_result"
     }
    ],
    "source": [
@@ -158,7 +119,7 @@
     "\n",
     "\n",
     "\n",
-    "for cutoff_energy in range(0,1050,50):\n",
+    "for cutoff_energy in range(0,550,100):\n",
     "\tnobrem_f = found[ak.all(found[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n",
     "\tnobrem_l = lost[ak.all(lost[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n",
     "\tprint(\"sample size: \",ak.num(nobrem_f,axis=0)+ak.num(nobrem_l,axis=0))\n",
@@ -175,8 +136,8 @@
     "nobrem_found = found[ak.all(found[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n",
     "nobrem_lost = lost[ak.all(lost[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n",
     "\n",
-    "print(\"\\nsample size: \",ak.num(nobrem_found,axis=0)+ak.num(nobrem_lost,axis=0))\n",
-    "t_eff(nobrem_found, nobrem_lost)"
+    "print(\"cutoff energy = 350MeV\\nsample size:\",ak.num(nobrem_found,axis=0)+ak.num(nobrem_lost,axis=0))\n",
+    "print(\"eff = \",t_eff(nobrem_found, nobrem_lost))"
    ]
   },
   {
@@ -188,7 +149,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 25,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
@@ -197,7 +158,7 @@
        "0.8593328191284226"
       ]
      },
-     "execution_count": 25,
+     "execution_count": 5,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -222,7 +183,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 26,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [
     {
@@ -243,7 +204,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 27,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [
     {
@@ -278,7 +239,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 28,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
@@ -315,7 +276,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 29,
+   "execution_count": 9,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -344,7 +305,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 30,
+   "execution_count": 10,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -390,7 +351,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 31,
+   "execution_count": 11,
    "metadata": {},
    "outputs": [
     {
@@ -414,7 +375,7 @@
        "<Array [101, 101, 101, 101, 101, ..., 101, 101, 101, 0] type='11 * float32'>"
       ]
      },
-     "execution_count": 31,
+     "execution_count": 11,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -426,7 +387,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 32,
+   "execution_count": 12,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -485,7 +446,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 33,
+   "execution_count": 13,
    "metadata": {},
    "outputs": [
     {
@@ -509,7 +470,7 @@
        "<Array [101, 101, 101, 101, 101, ..., 101, 101, 101, 0] type='11 * float64'>"
       ]
      },
-     "execution_count": 33,
+     "execution_count": 13,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -530,7 +491,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 34,
+   "execution_count": 14,
    "metadata": {},
    "outputs": [
     {
@@ -610,7 +571,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 35,
+   "execution_count": 15,
    "metadata": {},
    "outputs": [
     {
@@ -661,7 +622,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 36,
+   "execution_count": 16,
    "metadata": {},
    "outputs": [
     {