diff --git a/analysis/Simulation_Output.ipynb b/analysis/Simulation_Output.ipynb
index 9e62d07..163f73b 100644
--- a/analysis/Simulation_Output.ipynb
+++ b/analysis/Simulation_Output.ipynb
@@ -13,7 +13,6 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "import uproot\n",
     "import numpy             as np\n",
     "import matplotlib.pyplot as plt\n",
     "\n",
@@ -35,7 +34,7 @@
    "source": [
     "folder    = \"/pc/choozdsk01/palomare/GEANT4/megacell_v3/\"\n",
     "sensors   = [\"Arapuca\",\"SiPM1\",\"SiPM2\"]\n",
-    "contenido = [\"megacell_v3.root\",\"megacell_v3_support_free.root\"]\n",
+    "contenido = [\"megacell_v3.root\",\"megacell_v3_rebaba_free.root\"]\n",
     "debug     = True"
    ]
   },
@@ -53,7 +52,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "### DEFINE YOUR GEOMETRY\n",
+    "### Define your geometry and sensors\n",
     "Unitary vectors:\n",
     " - x = [1,0,0]\n",
     " - y = [0,1,0]\n",
@@ -92,7 +91,7 @@
    "outputs": [],
    "source": [
     "for d,my_data in enumerate(all_data):\n",
-    "    fig_ang = plot_variable_distributions(my_data,\"Phi\",stats=(True,False),probability=True,percentile=(0,1),log=(False,False),dpi=100,save=True,debug=False)\n",
+    "    fig_ang = plot_variable_distributions(my_data,\"X\",stats=(True,False),probability=False,percentile=(0,1),log=(False,True),dpi=100,save=True,debug=False)\n",
     "    plt.show()"
    ]
   },
@@ -103,18 +102,7 @@
    "outputs": [],
    "source": [
     "for d,my_data in enumerate(all_data):\n",
-    "    fig_ang = plot_variable_distributions(my_data,\"Theta\",stats=(True,False),probability=True,percentile=(0,1),log=(False,False),dpi=100,save=True,debug=False)\n",
-    "    plt.show()"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "for d,my_data in enumerate(all_data):\n",
-    "    fig_ang = plot_variable_distributions(my_data,\"FixedIncidenceAngleDegree\",stats=(True,False),probability=True,percentile=(0,1),log=(False,True),dpi=100,save=True,debug=False)\n",
+    "    fig_ang = plot_variable_distributions(my_data,\"FixedIncidenceAngleDegree\",stats=(True,False),probability=True,percentile=(0,1),log=(False,False),dpi=100,save=True,debug=False)\n",
     "    plt.show()"
    ]
   },
@@ -135,21 +123,12 @@
    "metadata": {},
    "outputs": [],
    "source": [
+    "surface = {\"Arapuca\": [\"X\",\"Z\"], \"SiPM1\": [\"Y\",\"Z\"], \"SiPM2\": [\"Y\",\"Z\"]}\n",
     "for d,my_data in enumerate(all_data): \n",
-    "    fig_ph = plot_photon_density(my_data,sensors_info,dpi=100)\n",
+    "    fig_ph = plot_photon_density(my_data,sensors_info,surface,density=False,bins=50,dpi=150,debug=True)\n",
     "    fig_ph.set_figwidth(20)\n",
     "    fig_ph.set_figheight(5)"
    ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "for d,my_data in enumerate(all_data): \n",
-    "    fig_acum = plot_acumhits(my_data,dpi=100,semilogy=True,bins=[100,200,200])"
-   ]
   }
  ],
  "metadata": {
diff --git a/data/attempt.root b/data/attempt.root
deleted file mode 100644
index 3059d39..0000000
Binary files a/data/attempt.root and /dev/null differ
diff --git a/production/megacell_v3_rebaba_free_prod.json b/production/megacell_v3_rebaba_free_prod.json
new file mode 100644
index 0000000..0e82332
--- /dev/null
+++ b/production/megacell_v3_rebaba_free_prod.json
@@ -0,0 +1,83 @@
+{
+  "comments":"Position units given in mm",
+
+  "display": false,
+  "G4_beam_settings": 
+  [
+    "/gps/particle opticalphoton",
+    "/gps/ene/type Mono",
+    "/gps/ene/mono 9.7 eV",
+    "/gps/pos/type Plane",
+    "/gps/pos/shape Circle",
+    "/gps/pos/radius  11 mm",
+    "/gps/pos/rot1 1 0 0",
+    "/gps/pos/rot2 0 0 1 ",
+    "/gps/pos/centre 0 74.25 0 mm",
+    "/gps/direction 0 -1 0",
+    "/gps/ang/type iso",
+    "/gps/number 192000",
+    "/run/beamOn 1000",
+    "/run/initialize"
+  ],
+
+  "G4_vis_settings":
+  [
+    "/vis/viewer/set/autorefresh/true",
+    "/vis/open OGL",
+    "/vis/drawVolume",
+    "/vis/scene/add/trajectories smooth",
+    "/vis/scene/add/eventID",
+    "/vis/enable",
+    "/vis/viewer/rebuild",
+    "/vis/viewer/set/background grey",
+    "/vis/scene/add/axes",
+    "/vis/scene/endOfEventAction accumulate 10"
+  ],
+
+  "sensors_keys": 
+  [
+    "XArapuca",
+    "SiPM1",
+    "SiPM2"
+  ],
+  
+  "source": 
+  {
+    "X": 0,
+    "Y": 0,
+    "Z": 0
+  },
+  
+  "big_cajita":
+  {
+    "outter_dim": [67.0,37.425,67.0],
+    "inner_dim":  [66.7,37.125,66.7],
+    "comment":   "Rotation and position should be left as they are",
+    "rot":        [0.0, 0.0, 0.0]
+  },
+  
+  "arapuca":
+  {
+    "dim": [71.875, 1e-6, 71.875],
+    "pos": [0.0, 0.5, 0.0],
+    "comment": "Rotation should follow the same as big_cajita",
+    "rot": [0.0 , 0.0 , 0.0]
+  },
+  
+  "sipms":
+  {
+    "dim": [1e-6, 3.0, 3.0],
+    "pos1":[ 28.54, 66.58, 0],
+    "pos2":[-28.54, 66.58, 0],
+    "rot": [0.0 , 0.0 , 0.0]
+  },
+
+  "alpha_support":
+  {
+    "dim":  [12.5,14.5,0.85],
+    "angle":[45,95],
+    "comment": "Y coordinate constrained to top of big_cajita",
+    "pos":    [0.0, 0.0, 0.0],
+    "rot":  [0.0, 0.0, 0.0]
+  }
+}
\ No newline at end of file
diff --git a/scripts/setup.sh b/scripts/setup.sh
index 3da6a16..94a8caf 100755
--- a/scripts/setup.sh
+++ b/scripts/setup.sh
@@ -15,4 +15,10 @@ source /cvmfs/sft.cern.ch/lcg/views/LCG_100/x86_64-centos7-gcc8-opt/setup.sh
 source /cvmfs/sft.cern.ch/lcg/views/LCG_100/x86_64-centos7-gcc8-opt/bin/geant4.sh
 
 # Make directory for results
-mkdir -p results
\ No newline at end of file
+mkdir -p results
+
+# Check if directory data exists
+if [ ! -d "data" ]; then
+    mkdir data
+    sshfs manthey@gaeuidc1.ciemat.es:/pc/choozdsk01/palomare/GEANT4/ data
+fi
\ No newline at end of file