diff --git a/mac/forbidden_dag_HO-May2019.mac b/mac/forbidden_dag_HO-May2019.mac
new file mode 100644
index 0000000000000000000000000000000000000000..74505eb40c4b949b78ccdced946cd56f0352ec93
--- /dev/null
+++ b/mac/forbidden_dag_HO-May2019.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2 -0.02049
+/VK/source/betaShapeFuncA3 -0.4022
+/VK/source/betaShapeFuncA4  0.1017
+/VK/source/betaShapeFuncA5 -0.01317
+/VK/source/betaShapeFuncA6  0.0008100
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-HO-May2019.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-HO-May2019.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-HO-May2019.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_final-1.mac b/mac/forbidden_dag_final-1.mac
new file mode 100644
index 0000000000000000000000000000000000000000..11a3393c12bbd352fef035917e76c02bb2b18082
--- /dev/null
+++ b/mac/forbidden_dag_final-1.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2  0.02666
+/VK/source/betaShapeFuncA3 -0.1903
+/VK/source/betaShapeFuncA4  0.02885
+/VK/source/betaShapeFuncA5 -0.004433
+/VK/source/betaShapeFuncA6  0.001540
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-CVC.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_final-2.mac b/mac/forbidden_dag_final-2.mac
new file mode 100644
index 0000000000000000000000000000000000000000..32884f89ab7a80e9479656469068965224c0fe2c
--- /dev/null
+++ b/mac/forbidden_dag_final-2.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2  0.03275
+/VK/source/betaShapeFuncA3 -0.2594
+/VK/source/betaShapeFuncA4  0.04412
+/VK/source/betaShapeFuncA5 -0.004527
+/VK/source/betaShapeFuncA6  0.001200
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-CVC.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_final-3.mac b/mac/forbidden_dag_final-3.mac
new file mode 100644
index 0000000000000000000000000000000000000000..b521b03dc0ab8b314f7101187609cf9f6a7b0b31
--- /dev/null
+++ b/mac/forbidden_dag_final-3.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2 -0.02049
+/VK/source/betaShapeFuncA3 -0.4022
+/VK/source/betaShapeFuncA4  0.1017
+/VK/source/betaShapeFuncA5 -0.01317
+/VK/source/betaShapeFuncA6  0.0008100
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-CVC.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_final-4.mac b/mac/forbidden_dag_final-4.mac
new file mode 100644
index 0000000000000000000000000000000000000000..26225849e7672500c486c01ac3ec4a34cd006bd4
--- /dev/null
+++ b/mac/forbidden_dag_final-4.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2 -0.02581
+/VK/source/betaShapeFuncA3 -0.3686
+/VK/source/betaShapeFuncA4  0.08695
+/VK/source/betaShapeFuncA5 -0.01066
+/VK/source/betaShapeFuncA6  0.0006468
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-CVC.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_final-5.mac b/mac/forbidden_dag_final-5.mac
new file mode 100644
index 0000000000000000000000000000000000000000..d0931bbc18d9dbac03d5db2456780d3746745b97
--- /dev/null
+++ b/mac/forbidden_dag_final-5.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2 -0.02572
+/VK/source/betaShapeFuncA3 -0.3697
+/VK/source/betaShapeFuncA4  0.08748
+/VK/source/betaShapeFuncA5 -0.01075
+/VK/source/betaShapeFuncA6  0.0006531
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-CVC.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_final-6.mac b/mac/forbidden_dag_final-6.mac
new file mode 100644
index 0000000000000000000000000000000000000000..3b28ea0d39fd9a90d01fcc2d2a4e652637ac7492
--- /dev/null
+++ b/mac/forbidden_dag_final-6.mac
@@ -0,0 +1,55 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta  90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1  1.00
+/VK/source/betaShapeFuncA2 -0.02876
+/VK/source/betaShapeFuncA3 -0.3504
+/VK/source/betaShapeFuncA4  0.07935
+/VK/source/betaShapeFuncA5 -0.009422
+/VK/source/betaShapeFuncA6  0.0005710
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-CVC.root 
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-CVC.root 
+/run/beamOn 1000000