From d5679171a24d1cb139239a70487db595ef026f58 Mon Sep 17 00:00:00 2001
From: oliskir <oliver.kirsebom@gmail.com>
Date: Thu, 27 Jul 2017 14:00:26 +0200
Subject: [PATCH] added Silicon as material
---
include/VKMaterials.hh | 1 +
mac/SiDet.mac | 81 ++++++++++++++++++++++++++++++++++++++++++
src/VKMaterials.cpp | 4 +++
3 files changed, 86 insertions(+)
create mode 100644 mac/SiDet.mac
diff --git a/include/VKMaterials.hh b/include/VKMaterials.hh
index eff48cb..d9c7172 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 0000000..0c7779b
--- /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 8f02e2d..984deb6 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);
--
GitLab