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