diff --git a/mac/forbidden_dag_HO.mac b/mac/forbidden_dag_HO.mac
new file mode 100644
index 0000000000000000000000000000000000000000..95ec9366c3f52ff1bfb928da9c68e20096fd6dc4
--- /dev/null
+++ b/mac/forbidden_dag_HO.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-HO.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.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.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag.mac b/mac/forbidden_dag_HOq.mac
similarity index 84%
rename from mac/forbidden_dag.mac
rename to mac/forbidden_dag_HOq.mac
index 73b52c93cd295fd85af0820ebc9d0582f20a7766..7f4b8f50691ffcdf517581f051d179a504850a4e 100644
--- a/mac/forbidden_dag.mac
+++ b/mac/forbidden_dag_HOq.mac
@@ -41,15 +41,15 @@
 
 /VK/source/betaMinEnergy 5.00 MeV
 /VK/field/maxbfield 0.48 tesla
-/VK/output/openFile output/forbidden_cut5000_b48_shape-dag.root 
+/VK/output/openFile output/forbidden_cut5000_b48_shape-dag-HOq.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.root 
+/VK/output/openFile output/forbidden_cut5000_b50_shape-dag-HOq.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.root 
+/VK/output/openFile output/forbidden_cut5000_b56_shape-dag-HOq.root 
 /run/beamOn 1000000
diff --git a/mac/forbidden_dag_WS.mac b/mac/forbidden_dag_WS.mac
new file mode 100644
index 0000000000000000000000000000000000000000..c69266af78e2ec3f4ed32e06aacecdf5482273dc
--- /dev/null
+++ b/mac/forbidden_dag_WS.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.03267
+/VK/source/betaShapeFuncA3 -0.3674
+/VK/source/betaShapeFuncA4  0.08979
+/VK/source/betaShapeFuncA5 -0.01166
+/VK/source/betaShapeFuncA6  0.0007539
+
+# 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-WS.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-WS.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-WS.root 
+/run/beamOn 1000000
diff --git a/mac/forbidden_dag_WSq.mac b/mac/forbidden_dag_WSq.mac
new file mode 100644
index 0000000000000000000000000000000000000000..b99d655a8e916d7f0d70b73befe92145666ad96e
--- /dev/null
+++ b/mac/forbidden_dag_WSq.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.03509
+/VK/source/betaShapeFuncA3 -0.3520
+/VK/source/betaShapeFuncA4  0.08322
+/VK/source/betaShapeFuncA5 -0.01056
+/VK/source/betaShapeFuncA6  0.0006820
+
+# 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-WSq.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-WSq.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-WSq.root 
+/run/beamOn 1000000