-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathREADME
126 lines (90 loc) · 4.32 KB
/
README
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
$Id: README 93739 2015-10-30 11:04:09Z gcosmo $
-------------------------------------------------------------------
=========================================================
Geant4 - an Object-Oriented Toolkit for Simulation in HEP
=========================================================
TestEm6
-------
This example is intended to test the processes of gamma conversion
to a pair of muons and annihilation of positrons with atomic
electrons to a pair of muons.
1- GEOMETRY DEFINITION
The geometry consists of a single block of a homogenous material.
Two parameters define the geometry :
- the material of the box,
- the (full) size of the box.
The default is 500 m of iron.
In addition a transverse uniform magnetic field can be applied.
The default geometry is constructed in DetectorConstruction class,
but all of the above parameters can be changed interactively via
the commands defined in the DetectorMessenger class.
2- PHYSICS LIST
Physics Lists are based on modular design. Several modules are
instantiated:
1. Transportation
2. EM physics
3. Decays
4. StepMax - for step limitation
The electromagnetic physics is chosen from one of the Geant4 EM
physics constructors in the physics_list library.
Cross sections can be enhanced (see below).
3- AN EVENT : THE PRIMARY GENERATOR
The primary kinematic consists of a single particle which hits the
block perpendicular to the input face. The type of the particle
and its energy are set in the PrimaryGeneratorAction class, and can
changed via the G4 build-in commands of G4ParticleGun class (see
the macros provided with this example).
The default is a Gamma of 100 TeV.
In addition one can choose randomly the impact point of the incident
particle. The corresponding interactive command is built in
PrimaryGeneratorMessenger class.
A RUN is a set of events.
4- VISUALIZATION
The Visualization Manager is set in the main() (see TestEm6.cc).
The initialisation of the drawing is done via the command
> /control/execute vis.mac
The detector has a default view which is a longitudinal view of the box.
The tracks are drawn at the end of event, and erased at the end of run.
Optionally one can choose to draw all particles, only the charged ones,
or none. This command is defined in EventActionMessenger class.
5- PHYSICS DEMO
The particle's type and the physics processes which will be available
in this example are set in PhysicsList class.
In addition a build-in interactive command (/process/inactivate procname)
allows to activate/inactivate the processes one by one.
The threshold for producing secondaries can be changed.
eg: /run/particle/setCut 100 micrometer
/run/initialize
To visualize the GammaConversionToMuons :
/control/execute run01.mac
/control/execute vis.mac
/run/beamOn
To visualize the AnnihiToMuPair :
/control/execute run11.mac
/control/execute vis.mac
/run/beamOn
6- HOW TO START ?
- execute Test in 'batch' mode from macro files
% TestEm6 run01.mac
- execute Test in 'interactive mode' with visualization
% TestEm6
....
Idle> type your commands
....
Idle> exit
7- HOW TO INCREASE STATISTICS ON gamma -> mu+mu- ?
The processes of gamma -> mu+mu- and e+e- -> mu+mu-
have a low cross section but can be important
for leakage through thick absorbers and calorimeters.
Straight forward simulation will be quite time consuming.
To make the processes more visible, the cross section can be
artificially increased by some factor (here 1000)
using the commands (only effective after /run/initialize)
/testem/phys/SetGammaToMuPairFac 1000
/testem/phys/SetAnnihiToMuPairFac 1000
8- HISTOGRAMS
Testem6 produces 6 histograms which illustrate the final state of
the GammaConversionToMuons process. See their definitions in RunAction.cc
By default the histograms are saved as testem6.root
The format of the histogram file can be : root (default), xml, csv,
by selecting g4nnn.hh in RunAction.hh