bla

allowed.mac

@@ -4,7 +4,6 @@
 
 # magnetic field strength
 /VK/field/maxbfield 0.48 tesla
-/VK/detector/volumes/blinkers true
 
 # source dimensions
 /VK/source/diskShapedSource true

compile.sh

 rm -f CMakeCache.txt
 rm -rf CMakeFiles
-cmake -DCMAKE_BUILD_TYPE=DEBUG -DGeant4_DIR=/home/oliskir/src/geant-4.10/geant4.10.3-install/lib/Geant4-10.3.0 $pwd
-#cmake -DCMAKE_BUILD_TYPE=DEBUG -DGeant4_DIR=/home/oliskir/src/geant4/geant4.10.2-install/lib/Geant4-10.2.0 $pwd -DCMAKE_BUILD_TYPE=DEBUG
+#cmake -DCMAKE_BUILD_TYPE=DEBUG -DGeant4_DIR=/home/oliskir/src/geant-4.10/geant4.10.3-install/lib/Geant4-10.3.0 $pwd
+cmake -DCMAKE_BUILD_TYPE=DEBUG -DGeant4_DIR=/home/oliskir/src/geant4/geant4.10.2-install/lib/Geant4-10.2.0 $pwd -DCMAKE_BUILD_TYPE=DEBUG
 make -j4
 

electron.mac

@@ -5,9 +5,6 @@
 /VK/output/zeroSuppression true
 /VK/storeTrajectoryData true
 
-# magnetic field
-/VK/field/maxbfield 0.145 tesla
-
 # source dimensions
 /VK/source/diskShapedSource true
 /VK/source/innerDiameter 0 cm
@@ -21,6 +18,9 @@
 
 #/gps/energy 975.657 keV
 /gps/energy 1682.232 keV
-/VK/output/openFile output/electron.root                
-/run/beamOn 100000
 
+# magnetic field
+/VK/field/maxbfield 0.145 tesla
+
+/VK/output/openFile output/electron.root                
+/run/beamOn 10000

init.mac

@@ -13,6 +13,7 @@
 /VK/detector/volumes/magnet true
 /VK/detector/volumes/catcherFoil true
 /VK/detector/volumes/shield true
+/VK/detector/volumes/blinkers true
 /VK/storeTrajectoryData false
 
 # field

src/VKPhysicsList.cpp

@@ -56,8 +56,10 @@ VKPhysicsList::VKPhysicsList()
 
     // EM Physics
     RegisterPhysics(new VKEMPhysics("standard EM"));
+    // RegisterPhysics(new G4EmStandardPhysics_option4()); // not yet tested
     // RegisterPhysics(new G4EmLivermorePhysics());  // overpredicts muon energy loss by a factor of 2!
     // RegisterPhysics(new G4EmLowEPPhysics());
+    // RegisterPhysics(new G4EmPenelopePhysics()); // not yet tested
     
     // Muon Physics
     RegisterPhysics(new VKMuonPhysics("muon"));

src/VKScintillatorSensitiveDetector.cpp

@@ -59,7 +59,7 @@ void VKScintillatorSensitiveDetector::Initialize(G4HCofThisEvent* hitsCE)
     // A way to keep all the hits of this event in one place if needed
     static G4int hitsCID = -1;
     if (hitsCID<0) {
-        hitsCID = GetCollectionID(0);
+        hitsCID = GetCollectionID(0);  // "CID" = Collection ID
     }
     hitsCE -> AddHitsCollection(hitsCID, fScintCollection);
 }
@@ -89,6 +89,7 @@ G4bool VKScintillatorSensitiveDetector::ProcessHits(G4Step* aStep,G4TouchableHis
     scintHit -> SetEdep(edep);
     scintHit -> SetPos(pos);
 
+    // Add hit to collection
     fScintCollection->insert(scintHit);
 
     return true;