Gabriele Cosmo
86 lines
2.7 KiB
Markdown
86 lines
2.7 KiB
Markdown
\page ExampleTestEm13 Example TestEm13
|
|
|
|
How to compute cross sections from the transmition coefficient.
|
|
|
|
## GEOMETRY DEFINITION
|
|
|
|
It is a single box representing a layer of finite thickness of
|
|
homogeneous material.
|
|
Two parameters define the geometry :
|
|
- the material of the box,
|
|
- the (full) size of the box.
|
|
|
|
The default geometry (1 cm of water) is constructed in
|
|
DetectorConstruction, but the above parameters can be changed
|
|
interactively via the commands defined in DetectorMessenger.
|
|
|
|
## PHYSICS LIST
|
|
|
|
The physics list contains the standard electromagnetic processes.
|
|
In order not to introduce 'artificial' constraints on the step size, the
|
|
multiple scattering is not instanciated, and all processes are
|
|
registered as discrete : there is no continuous energy loss.
|
|
|
|
## AN EVENT : THE PRIMARY GENERATOR
|
|
|
|
The primary kinematic consists of a single particle starting at the edge
|
|
of the box. The type of the particle and its energy are set in
|
|
PrimaryGeneratorAction (1 MeV gamma), and can be changed via the G4
|
|
build-in commands of G4ParticleGun class (see the macros provided with
|
|
this example).
|
|
|
|
## PHYSICS
|
|
|
|
An event is killed at the first step of the incident paticle.
|
|
Either the particle has interacted or is transmitted through the layer.
|
|
The cross section, also called absorption coefficient, is computed from
|
|
the rate of unaltered transmitted incident particles.
|
|
|
|
The result is compared with the 'input' data, i.e. with the cross
|
|
sections stored in the PhysicsTables and used by Geant4.
|
|
|
|
A set of macros defining various run conditions are provided.
|
|
The processes are actived/inactived in order to survey the processes
|
|
individually.
|
|
|
|
|
|
## VISUALIZATION
|
|
|
|
The Visualization Manager is set in the main().
|
|
The initialisation of the drawing is done via the commands
|
|
/vis/... in the macro vis.mac. To get visualisation:
|
|
```
|
|
> /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.
|
|
|
|
## HOW TO START ?
|
|
|
|
- Execute TestEm13 in 'batch' mode from macro files :
|
|
```
|
|
% ./TestEm13 compt.mac
|
|
```
|
|
|
|
- Execute TestEm13 in 'interactive mode' with visualization :
|
|
```
|
|
% ./TestEm13
|
|
Idle> control/execute vis.mac
|
|
....
|
|
Idle> type your commands
|
|
....
|
|
Idle> exit
|
|
```
|
|
|
|
Macros provided in this example:
|
|
- anni.mac: e+ (10 MeV) on 5 cm Aluminium
|
|
- compton.mac: gamma (300 keV) on 1 cm Aluminium
|
|
- conv.mac: gamma (20 MeV) on 5 mm Lead
|
|
- electron.mac: e- (100 MeV) on 1 cm Aluminium
|
|
- gamma.mac: gamma (100 keV) on 1 cm Water
|
|
- photoelec.mac: gamma (100 keV) on 1 cm Iron
|
|
- proton.mac: proton (1 GeV) on 1 cm Water
|
|
|