Version:
[ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]
1 -------------------------------------------------------------------
2
3 =========================================================
4 Geant4 - an Object-Oriented Toolkit for Simulation in HEP
5 =========================================================
6
7 Example ch2
8 -----------
9 A. Sytov, G. Paternò
10 INFN Ferrara Division, sytov@fe.infn.it, paterno@fe.infn.it
11
12 INTRODUCTION
13 Example ch2 is an enhanced version of ch1, providing the user with the full functionality of
14 both the G4ChannelingFastSimModel and G4BaierKatkov, with parameters set up via a macro,
15 in order to simulate the physics of channeling and channeling radiation/coherent bremsstrahlung.
16
17 The example can be exploited for a wide range of cases to study coherent effects in
18 a straight, bent or periodically bent crystal (crystalline undulator). Channeling
19 physics in ch2 is active for protons, ions, muons, pions, electrons and their antiparticles.
20 Any other charged particle can also be activated.
21
22 DESCRIPTION
23
24 The setup of the example ch2 is identical to ch1. As ch1, this example includes a bent crystal
25 and a detector positioned behind it. Like ch1, it is based on the experiments on
26 channeling [1] and channeling radiation [2] in a bent crystal, carried out at
27 Mainz Mikrotron MAMI with 855 MeV electrons. The experimental validation of
28 G4ChannelingFastSimModel is described in [3].
29
30 However, since ch2 parameters are fully set up in the macro run.mac, this example
31 is quite flexible and can be easily adapted for entirely different cases.
32
33 A description of all the available options is provided in run.mac.
34 It includes crystal and detector geometry, activation flags for
35 G4ChannelingFastSimModel and G4BaierKatkov and various options.
36
37 The example also provides detailed descriptions of various options for
38 G4ChannelingFastSimModel and G4BaierKatkov, which can adjust model parameters
39 depending on the specific case (see ConstructSDandField in DetectorConstruction).
40
41 The front surface of the crystal is placed at z=0 (with z as the beam direction),
42 while the front position of the detector can be set up via run.mac.
43
44 The output is recorded into the file results.root as a set of root ntuples.
45 These ntuples include:
46 - crystal: particles recorded at the crystal entrance,
47 - detector: all particles (except photons) recorded at the detector entrance,
48 - detector_photons: photons recorded at the detector entrance.
49
50 The format of every ntuple includes the following 10 variables (columns):
51
52 "eventID", "volume", "x", "y", "angle_x", "angle_y", "Ekin", "particle", "particleID", "parentID"
53
54 The variables represent:
55 - the event number within the run (column 0),
56 - the volume, either the crystal or the detector (column 1),
57 - the coordinate (x,y) and the angles (x'=dx/dz, y'=dy/dz) of the impinging particles (columns 2-6),
58 - the kinetic energy of the particle (column 7),
59 - the particle name (column 8),
60 - the particle ID (column 9),
61 - the parent ID of the particle (column 10).
62
63 To visualize these data, one should either open results.root using root TBrowser or use the python script analysis_ch2.py.
64
65 REFERENCES
66 [1] A. Mazzolari et al. Phys. Rev. Lett. 112, 135503 (2014).
67 [2] L. Bandiera et al. Phys. Rev. Lett. 115, 025504 (2015).
68 [3] A. Sytov et al. Journal of the Korean Physical Society 83, 132–139 (2023).
