diff --git a/results/.gitignore b/results/.gitignore
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/results/zdc_neutrons_reader.pdf b/results/zdc_neutrons_reader.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..a4588e6e3d7f805bf6ae5da718615352665eab4f
Binary files /dev/null and b/results/zdc_neutrons_reader.pdf differ
diff --git a/results/zdc_neutrons_reader.png b/results/zdc_neutrons_reader.png
new file mode 100644
index 0000000000000000000000000000000000000000..7e7c13fb695d72feda97178d2158a78d90470531
Binary files /dev/null and b/results/zdc_neutrons_reader.png differ
diff --git a/neutrons_zdc.cxx b/zdc_neutrons.cxx
similarity index 100%
rename from neutrons_zdc.cxx
rename to zdc_neutrons.cxx
diff --git a/zdc_neutrons_reader.cxx b/zdc_neutrons_reader.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..d6e8c5ee33ce55a8a9b7c5aa9d906c13a7de987b
--- /dev/null
+++ b/zdc_neutrons_reader.cxx
@@ -0,0 +1,44 @@
+#include "HepMC3/GenEvent.h"
+#include "HepMC3/ReaderAscii.h"
+#include "HepMC3/WriterAscii.h"
+#include "HepMC3/Print.h"
+
+#include "TH1F.h"
+#include <iostream>
+using namespace HepMC3;
+
+void zdc_neutrons_reader(){
+
+  //-------------------------------------
+  ReaderAscii hepmc_input("data/neutrons_zdc.hepmc");
+  int        events_parsed = 0;
+  GenEvent   evt(Units::GEV, Units::MM);
+
+  TH1F* h_neutron_energy = new TH1F("n energy","; E [GeV]",100,0,200);
+
+  while(!hepmc_input.failed()) {
+    // Read event from input file
+    hepmc_input.read_event(evt);
+
+    // If reading failed - exit loop
+    if( hepmc_input.failed() ) break;
+
+    for(const auto& v : evt.vertices() ) {
+      for(const auto& p : v->particles_out() ) {
+        if(p->pid() == 2112) {
+          h_neutron_energy->Fill(p->momentum().e());
+        }
+      }
+    }
+    evt.clear();
+    events_parsed++;
+  }
+  std::cout << "Events parsed and written: " << events_parsed << std::endl;
+
+  TCanvas* c = new TCanvas();
+  h_neutron_energy->Draw();
+  c->SaveAs("results/zdc_neutrons_reader.png");
+  c->SaveAs("results/zdc_neutrons_reader.pdf");
+}
+
+