| Geant4 Cross Reference |
1 ---------------------------------------------- 1 -------------------------------------------------------------------
>> 2 $Id: README,v 1.8 2008/06/02 10:00:40 sincerti Exp $
2 ---------------------------------------------- 3 -------------------------------------------------------------------
3 4
4 ========================================= 5 =========================================================
5 Geant4 - Nanobeam example 6 Geant4 - Nanobeam example
6 ========================================= 7 =========================================================
7 8
8 README file 9 README file
9 -------------------- 10 ----------------------
10 11
11 CORRESPONDING AUTHO 12 CORRESPONDING AUTHOR
12 13
13 S. Incerti (a, *) et al. 14 S. Incerti (a, *) et al.
14 a. Centre d'Etudes Nucleaires de Bordeaux-Grad 15 a. Centre d'Etudes Nucleaires de Bordeaux-Gradignan
15 (CENBG), IN2P3 / CNRS / Bordeaux 1 University, 16 (CENBG), IN2P3 / CNRS / Bordeaux 1 University, 33175 Gradignan, France
16 * e-mail:incerti@cenbg.in2p3.fr 17 * e-mail:incerti@cenbg.in2p3.fr
17 18
>> 19 Last modified by S. Incerti, 10/04/2008
>> 20
18 ---->1. INTRODUCTION. 21 ---->1. INTRODUCTION.
19 22
20 The nanobeam example simulates the beam optics 23 The nanobeam example simulates the beam optics of the nanobeam line installed
21 on the AIFIRA electrostatic accelerator facili 24 on the AIFIRA electrostatic accelerator facility located at CENBG,
22 Bordeaux-Gradignan, France. For more informati 25 Bordeaux-Gradignan, France. For more information on this facility,
23 please visit : 26 please visit :
24 http://www.cenbg.in2p3.fr/ 27 http://www.cenbg.in2p3.fr/
25 28
26 The code can be used to calculate : 29 The code can be used to calculate :
27 1) intrinsic aberration coefficients of the na 30 1) intrinsic aberration coefficients of the nanobeam line
28 2) beam image from a relasitic primary emittan 31 2) beam image from a relasitic primary emittance distribution
29 3) grid shadow images 32 3) grid shadow images
30 33
31 Three quadrupole field models can be used : 34 Three quadrupole field models can be used :
32 - a simple square field model 35 - a simple square field model
33 - a 3D mesh field model computed from OPERA3D 36 - a 3D mesh field model computed from OPERA3D
34 - an analytical model based on Enge's model 37 - an analytical model based on Enge's model
35 38
36 ---->2. GEOMETRY SET-UP. 39 ---->2. GEOMETRY SET-UP.
37 40
38 The full magnetic configuration of the nanobea 41 The full magnetic configuration of the nanobeam line is simulated.
39 This configuration is made of a combination of 42 This configuration is made of a combination of a doublet and triplet of
40 5 Oxford Microbeams Ltd. OM50 quadrupoles. 43 5 Oxford Microbeams Ltd. OM50 quadrupoles.
41 44
42 More details on the experimental setup and its 45 More details on the experimental setup and its simulation with Geant4 can
43 be found in the following papers: << 46 be found in the following papers, which may be found on the SLAC-SPIRES
>> 47 online database (http://www.slac.stanford.edu/spires/) :
44 48
45 - A DETAILED RAY-TRACING SIMULATION OF THE HIG 49 - A DETAILED RAY-TRACING SIMULATION OF THE HIGH RESOLUTION MICROBEAM AT THE
46 AIFIRA FACILITY 50 AIFIRA FACILITY
47 By F. Andersson, Ph. Barberet, S. Incerti, Ph. << 51 F. Andersson, Ph. Barberet, S. Incerti, Ph. Moretto (CENBG, Gradignan) . Dec 2007.
48 Published in Nucl.Instrum.Meth.B266:1653-1658, << 52 In press in Nucl.Instrum.Meth.B266:1653-1658, 2008
49 53
50 - MONTE CARLO SIMULATION OF THE CENBG MICROBEA 54 - MONTE CARLO SIMULATION OF THE CENBG MICROBEAM AND NANOBEAM LINES WITH THE
51 GEANT4 TOOLKIT 55 GEANT4 TOOLKIT
52 By S. Incerti, Q. Zhang, F. Andersson, Ph. Mor 56 By S. Incerti, Q. Zhang, F. Andersson, Ph. Moretto, G.W. Grime,
53 M.J. Merchant, D.T. Nguyen, C. Habchi, T. Pout 57 M.J. Merchant, D.T. Nguyen, C. Habchi, T. Pouthier and H. Seznec
54 Published in Nucl.Instrum.Meth.B260:20-27, 200 << 58 In press in Nucl.Instrum.Meth.B260:20-27, 2007
55 59
56 - GEANT4 SIMULATION OF THE NEW CENBG MICRO AND 60 - GEANT4 SIMULATION OF THE NEW CENBG MICRO AND NANO PROBES FACILITY
57 By S. Incerti, C. Habchi, Ph. Moretto, J. Oliv << 61 By S. Incerti, C. Habchi, Ph. Moretto, J. Olivier and H. Seznec. May 2006. 5pp.
58 Published in Nucl.Instrum.Meth.B249:738-742, 2 62 Published in Nucl.Instrum.Meth.B249:738-742, 2006
59 63
60 - A COMPARISON OF RAY-TRACING SOFTWARE FOR THE 64 - A COMPARISON OF RAY-TRACING SOFTWARE FOR THE DESIGN OF QUADRUPOLE MICROBEAM
61 SYSTEMS 65 SYSTEMS
62 By S. Incerti et al., 66 By S. Incerti et al.,
63 Published in Nucl.Instrum.Meth.B231:76-85, 200 67 Published in Nucl.Instrum.Meth.B231:76-85, 2005
64 68
65 ---->3 VISUALIZATION << 69 ---->3. SET-UP
>> 70
>> 71 - a standard Geant4 example GNUmakefile is provided
>> 72
>> 73 setup with:
>> 74 compiler = gcc-3.4.6
>> 75 G4SYSTEM = linux-g++
>> 76
>> 77 The following section gives the necessary environment variables.
>> 78
>> 79 ------->>3.1 ENVIRONMENT VARIABLES
>> 80
>> 81 All variables are defined with their default value.
>> 82
>> 83 - G4SYSTEM = Linux-g++
>> 84
>> 85 - G4INSTALL points to the installation directory of GEANT4;
>> 86
>> 87 - G4LIB point to the compiled libraries of GEANT4;
>> 88
>> 89 - G4WORKDIR points to the work directory;
>> 90
>> 91 - CLHEP_BASE_DIR points to the installation directory of CHLEP;
>> 92
>> 93 - G4LEDATA points to the low energy electromagnetic libraries;
>> 94
>> 95 - LD_LIBRARY_PATH = $CLHEP_BASE_DIR/lib
>> 96
>> 97 - G4LEVELGAMMADATA points to the photoevaporation library;
>> 98
>> 99 - G4NEUTRONHPDATA points to the neutron data files;
>> 100
>> 101 - G4RADIOACTIVEDATA points to the libraries for radio-active decay
>> 102 hadronic processes;
>> 103
>> 104 - G4ABLADATA points to the libraries fo ablation processes;
>> 105
>> 106 However, the
>> 107 $G4LEDATA,
>> 108 $G4LEVELGAMMADATA,
>> 109 $G4NEUTRONHPDATA,
>> 110 $G4RADIOACTIVEDATA
>> 111 and $G4ABLADATA
>> 112 variables do not need to be defined for this example.
>> 113
>> 114 Once these variables have been set, simply type gmake to compile the Nanobeam
>> 115 example.
>> 116
>> 117 ------->>3.2 VISUALIZATION
66 118
67 Visualization has not been implemented. 119 Visualization has not been implemented.
>> 120 All results are stored in text files and can be displayed with the provided
>> 121 ROOT macro file plot.C.
68 122
69 ---->4. HOW TO RUN THE EXAMPLE << 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
70 129
71 1) You must have compiled your Geant4 installa << 130 ---->4. HOW TO RUN THE EXAMPLE
72 CLHEP library which can handle matrix operatio <<
73 131
74 2) The code should be compiled cmake and run w << 132 The code can be run with :
75 133
76 ./nanobeam << 134 > $G4WORDIR/bin/$G4SYSTEM/Nanobeam
77 135
78 The macro file default.mac is read by default. 136 The macro file default.mac is read by default.
79 137
80 Several macro files are provided: << 138 It can be one of the following macro files :
81 139
82 1) for the computation of intrinsic aberration 140 1) for the computation of intrinsic aberration coefficients :
83 coef-square.mac : using square magnetic field << 141 coef-square.mac : using square magnetic field model (=default.mac)
84 coef-map.mac : using 3D map magnetic field mod 142 coef-map.mac : using 3D map magnetic field model
85 coef-enge.mac : using Enge's analytical field 143 coef-enge.mac : using Enge's analytical field model
86 144
87 2) for the simulation of the beam image with a 145 2) for the simulation of the beam image with a realistic emittance :
88 image-square.mac : using square magnetic field << 146 image-square.mac : using square magnetic field model
89 image-map.mac : using 3D map magnetic field mo 147 image-map.mac : using 3D map magnetic field model
90 image-enge.mac : using Enge's analytical field 148 image-enge.mac : using Enge's analytical field model
91 149
92 3) for the simulation of grid shadow images 150 3) for the simulation of grid shadow images
93 grid-square.mac : using square magnetic field 151 grid-square.mac : using square magnetic field model
94 grid-map.mac : using 3D map magnetic field mod 152 grid-map.mac : using 3D map magnetic field model
95 grid-enge.mac : using Enge's analytical field 153 grid-enge.mac : using Enge's analytical field model
96 154
97 These macros files are stored in the ./macros 155 These macros files are stored in the ./macros directory.
98 156
99 To run macros which include *map* in their nam <<
100 into the directory in which you run ./nanobeam <<
101 <<
102 The code can be run in MT mode, for high stati <<
103 Do not use MT for aberration coefficients calc <<
104 The switch to MT can be made in nanobeam.cc. <<
105 <<
106 ---->5. PHYSICS 157 ---->5. PHYSICS
107 158
108 The example runs with protons with fluctuating << 159 The example runs with protons with fluctuationg energies around 3 MeV.
109 Standard electromagnetic processes are activat << 160 Standard electromagnetic processes are activated by default.
110 Physics builder G4EmStandardPhysics), includin <<
111 161
112 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS 162 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS
113 163
114 All results are stored in the nanobeam.root fi << 164 The output results consists in several .txt files.
115 ROOT macro file plot.C: <<
116 * be sure to have ROOT installed on your machi <<
117 * be sure to be in the directory where ROOT ou <<
118 * copy plot.C into this directory <<
119 * launch ROOT by typing root, then under your <<
120 to execute the macro file <<
121 * or type directly: root plot.X <<
122 165
123 This macro file shows : << 166 These files can be easily analyzed using for example the provided ROOT macro
124 - the beam profile along the nanobeam line (on << 167 file plot.C; to do so :
125 coefficients) << 168 * be sure to have ROOT installed on your machine
126 - the beam image (Y vs X) on target << 169 * be sure to be in the nanobeam directory
127 - the beam emittance (THETA vs X) and (PHY vs << 170 * launch ROOT by typing root
128 - the grid shadow image (option) << 171 * under your ROOT session, type in : .X plot.C to execute the macro file
129 <<
130 The output ntuples can be written as xml or cs <<
131 172
132 ---------------------------------------------- 173 ---------------------------------------------------------------------------
133 174
134 Should you have any enquiry, please do not hes 175 Should you have any enquiry, please do not hesitate to contact:
135 incerti@cenbg.in2p3.fr 176 incerti@cenbg.in2p3.fr