nam/muon-tomograpy
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

target size can also be customized

Browse Source
This commit is contained in:
Nam Tran
2021-01-18 17:03:30 -06:00
parent 636f982750
commit ebfc30a532
4 changed files with 8 additions and 48 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
macros/batch_1M.mac
View File

@@ -1,3 +1,8 @@
# customize target parameters
/muTomo/targetLocation 0 10 0 cm
/muTomo/targetSize 10 10 20 cm
/muTomo/targetMaterial G4_W
/run/initialize /run/initialize
/tracking/verbose 0 /tracking/verbose 0
/run/beamOn 1000000 /run/beamOn 1000000

23
macros/run1.mac
View File

@@ -1,23 +0,0 @@
#
# Initialize kernel
/run/initialize
#
/control/verbose 2
/run/verbose 2
/event/verbose 0
/tracking/verbose 1
#
# gamma 6 MeV to the direction (0.,0.,1.)
#
/gun/particle gamma
/gun/energy 6 MeV
#
/run/beamOn 5
#
# proton 210 MeV to the direction (0.,0.,1.)
#
/gun/particle proton
/gun/energy 210 MeV
/tracking/verbose 2
#
/run/beamOn 1

22
macros/run2.mac
View File

@@ -1,22 +0,0 @@
#/run/numberOfWorkers 4
/run/initialize
#
/control/verbose 2
/run/verbose 2
#
# gamma 6 MeV to the direction (0.,0.,1.)
# 10000 events
#
/gun/particle gamma
/gun/energy 6 MeV
#
/run/printProgress 100
/run/beamOn 1000
#
# proton 210 MeV to the direction (0.,0.,1.)
# 1000 events
#
/gun/particle proton
/gun/energy 210 MeV
#
/run/beamOn 1000

6
src/DetectorConstruction.cc
View File

@@ -23,6 +23,9 @@
DetectorConstruction::DetectorConstruction() DetectorConstruction::DetectorConstruction()
: G4VUserDetectorConstruction(), targetMaterialStr("G4_Pb"), : G4VUserDetectorConstruction(), targetMaterialStr("G4_Pb"),
fCheckOverlaps(false) { fCheckOverlaps(false) {
targetX = 10 * cm;
targetY = 10 * cm;
targetZ = 5 * cm;
fDetConMessenger = new DetectorConstructionMessenger(this); fDetConMessenger = new DetectorConstructionMessenger(this);
// DefineMaterials(); // DefineMaterials();
} }
@@ -175,9 +178,6 @@ void DetectorConstruction::ComputeGeometry() {
concreteX = 30 * cm; concreteX = 30 * cm;
concreteY = 30 * cm; concreteY = 30 * cm;
concreteZ = 30 * cm; concreteZ = 30 * cm;
targetX = 10 * cm;
targetY = 10 * cm;
targetZ = 5 * cm;
worldSizeX = 1.2 * layerSizeX; worldSizeX = 1.2 * layerSizeX;
worldSizeY = 1.2 * layerSizeY; worldSizeY = 1.2 * layerSizeY;