initial commit of script to plot GENIE atmospheric Monte Carlo information

1 files changed, 21 insertions, 0 deletions

diff --git a/convert-genie-to-gst b/convert-genie-to-gst
new file mode 100755
index 0000000..2661e7d
--- /dev/null
+++ b/convert-genie-to-gst
@@ -0,0 +1,21 @@
+#!/usr/bin/env python
+from subprocess import check_call
+from os.path import split, splitext, join
+import os
+
+if __name__ == '__main__':
+    import argparse
+
+    parser = argparse.ArgumentParser("script to convert full GENIE root files to the reduced gst ROOT format")
+    parser.add_argument("filenames", nargs="+", help="GENIE root files")
+    parser.add_argument("--dest", required=True, help="destination directory")
+    args = parser.parse_args()
+
+    for filename in args.filenames:
+        head, tail = split(filename)
+        root, ext = splitext(tail)
+        output = join(args.dest, tail) + ".ntuple.root"
+        cmd = ["gntpc","-f","gst","-i",filename,"-o",output]
+        print(" ".join(cmd))
+        with open(os.devnull,"w") as devnull:
+            check_call(cmd, stdout=devnull, stderr=devnull)