diff --git a/include/VKMaterials.hh b/include/VKMaterials.hh
index eff48cbe7976487db464bfd5ba093f6eca6dc6b4..d9c717238375291e8a76878e49c7dbe263dbb4da 100644
--- a/include/VKMaterials.hh
+++ b/include/VKMaterials.hh
@@ -85,6 +85,7 @@ class VKMaterials {
     G4Element *fPb;
     
     G4Material *fLabVacuum;
+    G4Material *fSilicon;
     G4Material *fCopper;
     G4Material *fIron;
     G4Material *fBrass;
diff --git a/mac/SiDet.mac b/mac/SiDet.mac
new file mode 100644
index 0000000000000000000000000000000000000000..0c7779bd38bc1aa11d6e2310217c5d34262eee5c
--- /dev/null
+++ b/mac/SiDet.mac
@@ -0,0 +1,81 @@
+
+/control/execute mac/init.mac
+
+# setup
+/VK/detector/volumes/betaDetector true
+/VK/detector/volumes/magnet false
+/VK/detector/volumes/catcherFoil false
+/VK/detector/volumes/shield false
+/VK/detector/volumes/blinkers false
+/VK/storeTrajectoryData false
+
+# field
+/VK/field/enable false
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 6 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 0 -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/betaEndPoint 5.40 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 5.40 MeV
+
+# output 
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+# run simulations
+
+###/VK/field/maxbfield 0.10 tesla
+###/VK/output/openFile output/f20scan/080617/b10.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.15 tesla
+###/VK/output/openFile output/f20scan/080617/b15.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.20 tesla
+###/VK/output/openFile output/f20scan/080617/b20.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.25 tesla
+###/VK/output/openFile output/f20scan/080617/b25.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.30 tesla
+###/VK/output/openFile output/f20scan/080617/b30.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.35 tesla
+###/VK/output/openFile output/f20scan/080617/b35.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.40 tesla
+###/VK/output/openFile output/f20scan/080617/b40.root                
+###/run/beamOn 100000
+
+###/VK/field/maxbfield 0.44 tesla
+###/VK/output/openFile output/f20scan/080617/b44.root                
+###/run/beamOn 100000
+
+/VK/field/maxbfield 0.46 tesla
+/VK/output/openFile output/f20scan/080617/b46.root                
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/f20scan/080617/b48.root                
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/f20scan/080617/b50.root                
+/run/beamOn 1000000
diff --git a/src/VKMaterials.cpp b/src/VKMaterials.cpp
index 8f02e2de1128d5000ce8b9993e8863487d1a0e8c..984deb6c9c24ebdb66c0477f9863fb2233ad3280 100644
--- a/src/VKMaterials.cpp
+++ b/src/VKMaterials.cpp
@@ -121,6 +121,10 @@ void VKMaterials::CreateMaterials(){
     fLabVacuum->AddMaterial(fOxygen,   fractionalMass=0.2315);
     fLabVacuum->AddMaterial(fArgon,    fractionalMass=0.0128);
 
+    // Copper
+    fSilicon = new G4Material(name="SILICON", density=2.33*g/cm3, components=1);
+    fSilicon -> AddElement(fSi, 1);
+
     // Copper
     fCopper = new G4Material(name="COPPER", density=8.96*g/cm3, components=1);
     fCopper -> AddElement(fCu, 1);