diff --git a/clustering/emcal_electrons.sh b/clustering/emcal_electrons.sh
index d5bc5abaae442c582a00cdd1f66fdd179eff35f7..8302547fad67de8c77721406cc4ffd2a1e2c9b03 100644
--- a/clustering/emcal_electrons.sh
+++ b/clustering/emcal_electrons.sh
@@ -10,8 +10,8 @@ ls -lrth
 root -b -q "datasets/emcal_electrons.cxx(1e4, \"emcal_uniform_electrons.hepmc\")"
 cd topside && ls -l
 npsim --runType batch --numberOfEvents 1000 --compactFile topside.xml --inputFiles ../emcal_uniform_electrons.hepmc --outputFile  ../sim_output/sim_emcal_electrons.root
-cd ..
 pwd
 /usr/local/bin/xenv --xml /usr/local/Juggler.xenv /usr/local/scripts/gaudirun.py ./options/example_crystal.py
-root -b -q ./scripts/makeplot.C 
+cd ..
+root -b -q "./scripts/makeplot.C(\"sim_output/sim_emcal_electrons_output.root\")"
 
diff --git a/scripts/makeplot.C b/scripts/makeplot.C
index c49cbaa26b1abdd0f02eb3006975c213223d8bc3..638aa5c6c48438e53377fbe7956704f32d462709 100644
--- a/scripts/makeplot.C
+++ b/scripts/makeplot.C
@@ -2,7 +2,7 @@
 // Read reconstruction ROOT output file
 // Plot variables
 ////////////////////////////////////////
-int makeplot(const char* input_fname = "sim_output/rec_electron_0GeVto30GeV_100k_output.root")
+int makeplot(const char* input_fname = "../sim_output/sim_emcal_electrons_output.root")
 {
   // Setting figures
   gROOT->SetStyle("Plain");