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 =========================================================
5 Geant4 - Nanobeam example
6 =========================================================
7
8 README file
9 ----------------------
10
11 CORRESPONDING AUTHOR
12
13 S. Incerti (a, *) et al.
14 a. Centre d'Etudes Nucleaires de Bordeaux-Gradignan
15 (CENBG), IN2P3 / CNRS / Bordeaux 1 University, 33175 Gradignan, France
16 * e-mail:incerti@cenbg.in2p3.fr
17
18 ---->1. INTRODUCTION.
19
20 The nanobeam example simulates the beam optics of the nanobeam line installed
21 on the AIFIRA electrostatic accelerator facility located at CENBG,
22 Bordeaux-Gradignan, France. For more information on this facility,
23 please visit :
24 http://www.cenbg.in2p3.fr/
25
26 The code can be used to calculate :
27 1) intrinsic aberration coefficients of the nanobeam line
28 2) beam image from a relasitic primary emittance distribution
29 3) grid shadow images
30
31 Three quadrupole field models can be used :
32 - a simple square field model
33 - a 3D mesh field model computed from OPERA3D
34 - an analytical model based on Enge's model
35
36 ---->2. GEOMETRY SET-UP.
37
38 The full magnetic configuration of the nanobeam line is simulated.
39 This configuration is made of a combination of a doublet and triplet of
40 5 Oxford Microbeams Ltd. OM50 quadrupoles.
41
42 More details on the experimental setup and its simulation with Geant4 can
43 be found in the following papers:
44
45 - A DETAILED RAY-TRACING SIMULATION OF THE HIGH RESOLUTION MICROBEAM AT THE
46 AIFIRA FACILITY
47 By F. Andersson, Ph. Barberet, S. Incerti, Ph. Moretto
48 Published in Nucl.Instrum.Meth.B266:1653-1658, 2008
49
50 - MONTE CARLO SIMULATION OF THE CENBG MICROBEAM AND NANOBEAM LINES WITH THE
51 GEANT4 TOOLKIT
52 By S. Incerti, Q. Zhang, F. Andersson, Ph. Moretto, G.W. Grime,
53 M.J. Merchant, D.T. Nguyen, C. Habchi, T. Pouthier and H. Seznec
54 Published in Nucl.Instrum.Meth.B260:20-27, 2007
55
56 - GEANT4 SIMULATION OF THE NEW CENBG MICRO AND NANO PROBES FACILITY
57 By S. Incerti, C. Habchi, Ph. Moretto, J. Olivier and H. Seznec
58 Published in Nucl.Instrum.Meth.B249:738-742, 2006
59
60 - A COMPARISON OF RAY-TRACING SOFTWARE FOR THE DESIGN OF QUADRUPOLE MICROBEAM
61 SYSTEMS
62 By S. Incerti et al.,
63 Published in Nucl.Instrum.Meth.B231:76-85, 2005
64
65 ---->3 VISUALIZATION
66
67 Visualization has not been implemented.
68
69 ---->4. HOW TO RUN THE EXAMPLE
70
71 1) You must have compiled your Geant4 installation with the FULL version of the
72 CLHEP library which can handle matrix operations.
73
74 2) The code should be compiled cmake and run with :
75
76 ./nanobeam
77
78 The macro file default.mac is read by default.
79
80 Several macro files are provided:
81
82 1) for the computation of intrinsic aberration coefficients :
83 coef-square.mac : using square magnetic field model
84 coef-map.mac : using 3D map magnetic field model
85 coef-enge.mac : using Enge's analytical field model
86
87 2) for the simulation of the beam image with a realistic emittance :
88 image-square.mac : using square magnetic field model (=default.mac)
89 image-map.mac : using 3D map magnetic field model
90 image-enge.mac : using Enge's analytical field model
91
92 3) for the simulation of grid shadow images
93 grid-square.mac : using square magnetic field model
94 grid-map.mac : using 3D map magnetic field model
95 grid-enge.mac : using Enge's analytical field model
96
97 These macros files are stored in the ./macros directory.
98
99 To run macros which include *map* in their name, copy the file OM50.grid
100 into the directory in which you run ./nanobeam.
101
102 The code can be run in MT mode, for high statistics image simulation.
103 Do not use MT for aberration coefficients calculation (32 rays only are shot).
104 The switch to MT can be made in nanobeam.cc.
105
106 ---->5. PHYSICS
107
108 The example runs with protons with fluctuating energies around 3 MeV.
109 Standard electromagnetic processes are activated by default (corresponding to the
110 Physics builder G4EmStandardPhysics), including the G4StepLimiter process.
111
112 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS
113
114 All results are stored in the nanobeam.root file and can be displayed with the provided
115 ROOT macro file plot.C:
116 * be sure to have ROOT installed on your machine
117 * be sure to be in the directory where ROOT output files are generated
118 * copy plot.C into this directory
119 * launch ROOT by typing root, then under your ROOT session, type in : .X plot.C
120 to execute the macro file
121 * or type directly: root plot.X
122
123 This macro file shows :
124 - the beam profile along the nanobeam line (only for the computation of intrinsic
125 coefficients)
126 - the beam image (Y vs X) on target
127 - the beam emittance (THETA vs X) and (PHY vs Y) on target
128 - the grid shadow image (option)
129
130 The output ntuples can be written as xml or csv files, by changing the G4AnalysisManager default file type in RunAction::BeginOfRunAction().
131
132 ---------------------------------------------------------------------------
133
134 Should you have any enquiry, please do not hesitate to contact:
135 incerti@cenbg.in2p3.fr