diff --git a/km3buu/cmd.py b/km3buu/cmd.py
index 5509d91d6f10cca6d69ae02382cc2d129d8353af..a32b5a2fcf8d34a3ebdf4890abde68e8b617718b 100755
--- a/km3buu/cmd.py
+++ b/km3buu/cmd.py
@@ -117,6 +117,13 @@ ARGPARSE_GENERAL_PARAMS = [{
     "action": argparse.BooleanOptionalAction,
     "help": "Do tau propagation",
     "default": False
+}, {
+    "option_strings": ["--decay"],
+    "dest": "decay",
+    "action": argparse.BooleanOptionalAction,
+    "help":
+    "Decay final state particles (according to the decays done in gSeaGen)",
+    "default": True
 }, {
     "option_strings": ["--zenith", "-z"],
     "dest": "zenith",
@@ -201,7 +208,7 @@ def main():
                                    args.target,
                                    seed=args.seed,
                                    fluxfile=fluxfile,
-                                   do_decay=False)
+                                   do_decay=args.decay)
 
     jc["neutrinoanalysis"]["outputEvents"] = True
     jc["neutrinoanalysis"]["inclusiveAnalysis"] = False