PHONON_sim/gamma_source_scint.mac at main · jzettle/PHONON_sim · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
# Macro file for example B2
#
# Can be run in batch, without graphic
# or interactively: Idle> /control/execute run1.mac
#
# Change the default number of workers (in multi-threading mode)
#/run/numberOfThreads 4
#
/PHONON/det/setGDMLFile gdml/detector_shield/ScintillatorShield.gdml
/PHONON/det/setType scintillator
/PHONON/run/useHDF5 false
#/PHONON/run/setOutputFile shielded_testfiletypesetting
#/PHONON/run/setOutputFile shielded_testfiletypesetting_higherEprimary
/PHONON/run/setOutputFile gamma_spectrum.root
#/PHONON/run/setOutputFile scrap.root

# Initialize kernel
/run/initialize
#
# Default kinematics:
# proton 3 GeV in direction (0.,0.,1.)
# 1 event with tracking/verbose
#
#/tracking/verbose 1
#/run/beamOn 1
#
# 1 event with printing hits
#
#/tracking/verbose 0
#/hits/verbose 2
#/run/beamOn 1
#
# set target and chamber material
#
#/B2/det/setTargetMaterial G4_WATER
#/B2/det/setChamberMaterial G4_Ar
#/run/beamOn 3
#
# set a magnetic field and max allowed step length
# 3 event with printing hits
#
#/globalField/verbose 1
#/globalField/setValue 0.2 0 0 tesla
#/B2/det/stepMax 1.0 mm
#/gun/energy 0.3 GeV
#/run/beamOn 3
#
# muon 300 MeV in the direction (0.,0.,1.)
# 3 event with detailed printing
#

/g4cmp/producePhonons 0
/g4cmp/produceCharges 0
/g4cmp/maxLukePhonons 0
#/g4cmp/sampleLuke 0.00001
/g4cmp/phononBounces 1
/g4cmp/minEPhonons 0.001 eV
#/g4cmp/verbose 1

/process/inactivate phononScattering
/process/inactivate phononDownconversion

/PHONONEvent/offset 0

/gps/particle gamma
#/gps/ene/type Mono
#/gps/energy 8.0 MeV
#/gps/energy 0.03 eV
#/gps/particle phononL
#/gps/pos/type Point
#/gps/ang/type iso
#/gps/number 1
/gps/ene/type Lin
/gps/ene/min 0.0 MeV
/gps/ene/max 10.0 MeV
/gps/ene/gradient 0.0
/gps/ene/intercept 0.001
/gps/pos/type Plane
/gps/pos/shape Square
/gps/pos/centre -0.51 0.0 0.0 m
/gps/pos/halfx 0.5 m
/gps/pos/halfy 0.5 m
/gps/pos/rot1 0 0 1
/gps/pos/rot2 0 1 0
/gps/ang/type cos
#/gps/pos/confine NeutronPlaneVol

/run/beamOn 10000000
#/run/beamOn 1000