added Silicon as material

include/VKMaterials.hh

@@ -85,6 +85,7 @@ class VKMaterials {
     G4Element *fPb;
     
     G4Material *fLabVacuum;
+    G4Material *fSilicon;
     G4Material *fCopper;
     G4Material *fIron;
     G4Material *fBrass;

mac/SiDet.mac

+
+/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

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);