forked from ActarSimGroup/Actarsim
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathkine_batch.mac
279 lines (272 loc) · 9.06 KB
/
kine_batch.mac
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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
################################################################
#*-- AUTHOR : Hector Alvarez-Pol
#*-- Date: 05/2005
#*-- Last Update: 15/05/08
#*-- Copyright: GENP (Univ. Santiago de Compostela)
# --------------------------------------------------------------
# Comments:
# - 15/05/08 Multidetector geometry
# - 05/05/06 Updating to new ActarSim (geant4.8) code
# - 22/11/05 Updated including CINE options
# - 31/05/05 Macro for ACTAR simulation
#
################################################################
# Macro file for testing online jobs
################################################################
# verbosity levels and saveHistory
/control/verbose 1
/control/saveHistory
/run/verbose 0
/event/verbose 0
/tracking/verbose 0
#
# Setting the Physics modules; valid values are here listed:
# em: standard, lowenergy, penelope, (choose one from this three)
# common: decay,
# hadronic: elastic, binary, binary_ion, gamma_nuc,
# ion low-energy: ion-LowE, ion-LowE-ziegler1977, ion-LowE-ziegler1985,
# ion-LowE-ziegler2000, ion-standard
#
#/ActarSim/phys/addPhysics local
#/ActarSim/phys/addPhysics standard
#/ActarSim/phys/addPhysics standardGS
#/ActarSim/phys/addPhysics standardSS
#/ActarSim/phys/addPhysics standardWVI
#/ActarSim/phys/addPhysics decay
#/ActarSim/phys/addPhysics elastic
#/ActarSim/phys/addPhysics binary
#/ActarSim/phys/addPhysics binary_ion
#/ActarSim/phys/addPhysics gamma_nuc
#/ActarSim/phys/addPhysics lowenergy
#/ActarSim/phys/addPhysics ion-LowE
#/ActarSim/phys/addPhysics ion-LowE-ziegler1977
#/ActarSim/phys/addPhysics ion-LowE-ziegler1985
#/ActarSim/phys/addPhysics ion-LowE-ziegler2000
#/ActarSim/phys/addPhysics ion-standard
#/ActarSim/phys/addPhysics penelope
#/ActarSim/phys/addPhysics emstandard_opt0
/ActarSim/phys/addPhysics emstandard_opt3
#/ActarSim/phys/addPhysics emstandard_opt4
/ActarSim/phys/addPhysics ionGasModels
#/process/eLoss/useCutAsFinalRange
#/process/eLoss/fluct false
#/process/eLoss/subsec
#/process/eLoss/minsubsec
#/process/eLoss/StepFunction 0.2 1. mm
#/process/eLoss/minKinEnergy 100. MeV
#/process/eLoss/maxKinEnergy 300. MeV
#/process/eLoss/verbose 1
#/process/eLoss/linLossLimit 0.5
#
# Cuts for the particles (incomplete list, see README)
#
# /ActarSim/phys/setGCut 0.1
# /ActarSim/phys/setECut 0.1
# /ActarSim/phys/setPCut 0.1
# /ActarSim/phys/setCuts 1. mm
/ActarSim/phys/verbose 0
#
# Initialization is moved here from the main allowing PhysicsList
#
/run/initialize
#
# DETECTOR CHARACTERIZATION
#
# GENERAL COMMANDS
#
# Control of the material outside the chamber (by default Air)
#/ActarSim/det/setMediumMat Galactic
#/ActarSim/det/setMediumMat Air
#/ActarSim/det/setMediumMat Water
#Electric and Magnetic fields
#/ActarSim/det/setEleField 0 5e-3 0
#/ActarSim/det/setMagField 0 0 0 T
#
#
# GAS material:
#
/ActarSim/det/gas/setGasMat CH4_STP
#
#
/ActarSim/det/MaikoGeoIncludedFlag off
#
# GAS DETECTOR
#
/ActarSim/det/gasGeoIncludedFlag on
# if box
/ActarSim/det/gas/setDetectorGeometry box
/ActarSim/det/gas/luciteBoxIncluded off
#
#/ActarSim/det/gas/setXLengthGasBox 32. mm
#/ActarSim/det/gas/setYLengthGasBox 80. mm
#/ActarSim/det/gas/setZLengthGasBox 64. mm
#
/ActarSim/det/gas/setXLengthGasBox 37. mm
/ActarSim/det/gas/setYLengthGasBox 80.9 mm
/ActarSim/det/gas/setZLengthGasBox 69. mm
#
# Beam shield? tube or off
/ActarSim/det/gas/setBeamShield off
#
# SILICON DETECTOR
#
/ActarSim/det/silGeoIncludedFlag on
#Options for Silicon and scintillator coverage:
# 6 bits to indicate which sci wall is present (1) or absent (0)
# order is:
# bit1 (lsb) beam output wall 1
# bit2 lower (gravity based) wall 2
# bit3 upper (gravity based) wall 4
# bit4 left (from beam point of view) wall 8
# bit5 right (from beam point of view) wall 16
# bit6 (msb) beam entrance wall 32
# examples: 63 full coverage; 3 only output and bottom walls ...
/ActarSim/det/sil/sideCoverage 25
/ActarSim/det/sil/xBoxHalfLength 95. mm
/ActarSim/det/sil/yBoxHalfLength 105. mm
/ActarSim/det/sil/zBoxHalfLength 120. mm
/ActarSim/det/sil/print
#
# SCINTILLATOR DETECTOR
#
/ActarSim/det/sciGeoIncludedFlag off
# see above explanation in the equivalent command for the Silicons
/ActarSim/det/sci/sideCoverage 25
/ActarSim/det/sci/xBoxHalfLength 95. mm
/ActarSim/det/sci/yBoxHalfLength 105. mm
/ActarSim/det/sci/zBoxHalfLength 120. mm
/ActarSim/det/sci/print
#
#Control of the output on the ROOT file
#all the tracks are stored (note: huge space comsumption)
#Note: it should come before the update!!!
/ActarSim/analControl/storeTracks off
/ActarSim/analControl/storeTrackHistos on
/ActarSim/analControl/storeEvents on
/ActarSim/analControl/storeHistograms on
/ActarSim/analControl/storeSimpleTracks on
#/ActarSim/analControl/setMinStrideLength 1.
#/ActarSim/phys/stepMax 1000000.
#
# Update is mandatory after any material,field or detector change
#
/ActarSim/det/update
#/ActarSim/det/print
#
# Control of the primary events
#For all cases the possibility to have realistic beam distribution
/ActarSim/gun/beamInteraction on
/ActarSim/gun/realisticBeam off
#/ActarSim/gun/beamRadiusAtEntrance 1.5 mm
#/ActarSim/gun/emittance 200.0
#
# Realistic Event-Generator on
/ActarSim/gun/reactionFromEvGen off
#
# Reaction from Cine and Event-Generator
/ActarSim/gun/reactionFromCine off
/ActarSim/gun/reactionFromCrossSection off
#
# A) Track a particle or set of particles defined from the Particles list
#
#/ActarSim/gun/List
#/ActarSim/gun/particle proton
#if you want to use an ion, write "ion" in the previous command
#and set the Z, A and charge state in the next...
#/ActarSim/gun/ion 3 11 3
#/ActarSim/gun/energy 11. MeV
#/ActarSim/gun/position 0 0 -56 mm
#/ActarSim/gun/direction 0 0 1
#/ActarSim/gun/time 0
#/ActarSim/gun/polarization 0
#/ActarSim/gun/number 1
/ActarSim/gun/randomVertexZPosition off
#/ActarSim/gun/randomVertexZRange 0 70 mm
/ActarSim/gun/vertexZPosition 69 mm
#
# B) Track a predefined reaction from a file:
#
/ActarSim/gun/reactionFromFile off
#
#C) Track a reacion calculated from CINE (program from W. Mittig)
#
/ActarSim/gun/reactionFromCine off
#
#D) Track a reacion calculated from KINE
#8He(12C,16C*)4He
#IMPORTANT: SINCE WE DO NOT HAVE MESSENGER FOR DECAY JUST REMEMBER TO CHANGE EXCITATION ENERGY OF MOTHER NUCLEI IN PRIMARYGENERATOR (AND ALSO GAS). FOR KINEMATICS CHECK THE FLAG AND THE MASSES OF DECAY PRODUCTS
# Record masses here 11Be:11.021657749 D:2.01410177785 10Li:10.035481259 3He:3.01602931914 9Li:9.026789 t:3.01604928 p:1.007825 n:1.008665 4He:4.00260325 12Be:12.026921 10Be:10.013533
# For Two body reaction+decay we chose CASE * (with Kine). For two body reaction we chose CASE 4 with Kine. We have to modify this at PrimaryGenerationAction until we create a messenger.
# Added realistic excitation energy for scattered particle. Modify excitation energy in PrimaryGenerationAction
#/ActarSim/gun/reactionFromDecay off
/ActarSim/gun/reactionFromKine on
#
#/ActarSim/gun/Kine/incidentIon 1 2 1 0.0 2.0141
#/ActarSim/gun/Kine/incidentIon 2 4 2 0.0 4.001506179125
#/ActarSim/gun/Kine/incidentIon 2 4 2 0.0 4.00260325
/ActarSim/gun/Kine/incidentIon 6 13 6 0.0 13.00335483778
#/ActarSim/gun/Kine/incidentIon 50 132 50 0.0 131.9178
/ActarSim/gun/Kine/targetIon 1 1 1 0.0 1.00782503207
#/ActarSim/gun/Kine/targetIon 1 2 1 0.0 2.0141
#/ActarSim/gun/Kine/scatteredIon 1 2 1 0.0 2.0141
#/ActarSim/gun/Kine/scatteredIon 2 4 2 0.0 4.001506179125
#/ActarSim/gun/Kine/scatteredIon 2 4 2 0.0 4.00260325
/ActarSim/gun/Kine/scatteredIon 6 13 6 0.0 13.00335483778
#/ActarSim/gun/Kine/scatteredIon 50 133 50 0.854 132.9238
/ActarSim/gun/Kine/recoilIon 1 1 1 0.0 1.00782503207
/ActarSim/gun/Kine/labEnergy 650 MeV
/ActarSim/gun/Kine/randomThetaCM on
/ActarSim/gun/Kine/randomThetaRange 0.0 180.0
/ActarSim/gun/Kine/randomPhiAngle off
#/ActarSim/gun/Kine/userThetaCM 90.0 deg
#/ActarSim/gun/Kine/userPhiAngle 0. deg
#/ActarSim/gun/Kine/vertexPosition 0. 0. 7 cm
#
# VISUALIZATION
#
# Draw the whole geometry tree with details as function of verbosity
#/vis/ASCIITree/verbose 10
#/vis/drawTree
# visualization
#/vis/scene/create
#/vis/open OGLIX 600x600-0+0
#/vis/open OGLSQt 600x600-0+0
#/vis/open VRML2FILE
#/vis/viewer/set/autoRefresh 0
#/vis/viewer/flush
# set camera
#/vis/viewer/reset
#/vis/viewer/set/hiddenEdge 1
#/vis/viewer/set/lightsThetaPhi 120 40
#/vis/viewer/set/viewpointThetaPhi 115. 145.
#/vis/viewer/set/viewpointThetaPhi 90. 90.
#/vis/viewer/zoom 1.0
#/vis/viewer/set/background 1 1 1 1
#/vis/viewer/flush
#
# drawing style
#/vis/viewer/set/style surface
#/vis/viewer/set/style wireframe
#/vis/viewer/flush
#
#/vis/modeling/trajectories/create/drawByParticleID
#/vis/modeling/trajectories/drawByParticleID-0/set gamma red
#/vis/modeling/trajectories/drawByParticleID-0/setRGBA e- 1 0 1 1
#/vis/modeling/trajectories/drawByParticleID-0/setRGBA genericIon 1 0 1 1
#
# drawing the tracks
#/tracking/storeTrajectory 1
#/vis/scene/endOfEventAction accumulate
#/vis/viewer/set/autoRefresh 1
#
# create an empty scene and add the detector geometry to it
#/vis/drawVolume
#/vis/scene/add/axes 0 0 0 0.1 m
#/vis/scene/add/trajectories
#/vis/scene/add/hits
#/ActarSim/event/drawTracks all
/ActarSim/event/printModulo 1000
#
# RUN: number of events
#/run/beamOn 100000
/run/beamOn 10000