| Geant4 Cross Reference |
1 ---------------------------------------------- 1 -------------------------------------------------------------------
>> 2 $Id: README,v 1.9 2010-10-06 12:16:59 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, 06/10/2010
>> 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:
44 47
45 - A DETAILED RAY-TRACING SIMULATION OF THE HIG 48 - A DETAILED RAY-TRACING SIMULATION OF THE HIGH RESOLUTION MICROBEAM AT THE
46 AIFIRA FACILITY 49 AIFIRA FACILITY
47 By F. Andersson, Ph. Barberet, S. Incerti, Ph. << 50 F. Andersson, Ph. Barberet, S. Incerti, Ph. Moretto (CENBG, Gradignan) . Dec 2007.
48 Published in Nucl.Instrum.Meth.B266:1653-1658, 51 Published in Nucl.Instrum.Meth.B266:1653-1658, 2008
49 52
50 - MONTE CARLO SIMULATION OF THE CENBG MICROBEA 53 - MONTE CARLO SIMULATION OF THE CENBG MICROBEAM AND NANOBEAM LINES WITH THE
51 GEANT4 TOOLKIT 54 GEANT4 TOOLKIT
52 By S. Incerti, Q. Zhang, F. Andersson, Ph. Mor 55 By S. Incerti, Q. Zhang, F. Andersson, Ph. Moretto, G.W. Grime,
53 M.J. Merchant, D.T. Nguyen, C. Habchi, T. Pout 56 M.J. Merchant, D.T. Nguyen, C. Habchi, T. Pouthier and H. Seznec
54 Published in Nucl.Instrum.Meth.B260:20-27, 200 57 Published in Nucl.Instrum.Meth.B260:20-27, 2007
55 58
56 - GEANT4 SIMULATION OF THE NEW CENBG MICRO AND 59 - GEANT4 SIMULATION OF THE NEW CENBG MICRO AND NANO PROBES FACILITY
57 By S. Incerti, C. Habchi, Ph. Moretto, J. Oliv << 60 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 61 Published in Nucl.Instrum.Meth.B249:738-742, 2006
59 62
60 - A COMPARISON OF RAY-TRACING SOFTWARE FOR THE 63 - A COMPARISON OF RAY-TRACING SOFTWARE FOR THE DESIGN OF QUADRUPOLE MICROBEAM
61 SYSTEMS 64 SYSTEMS
62 By S. Incerti et al., 65 By S. Incerti et al.,
63 Published in Nucl.Instrum.Meth.B231:76-85, 200 66 Published in Nucl.Instrum.Meth.B231:76-85, 2005
64 67
65 ---->3 VISUALIZATION 68 ---->3 VISUALIZATION
66 69
67 Visualization has not been implemented. 70 Visualization has not been implemented.
68 71
69 ---->4. HOW TO RUN THE EXAMPLE 72 ---->4. HOW TO RUN THE EXAMPLE
70 73
71 1) You must have compiled your Geant4 installa << 74 The variable G4ANALYSIS_USE must be set to 1.
72 CLHEP library which can handle matrix operatio <<
73 75
74 2) The code should be compiled cmake and run w << 76 In order to generate histograms, at least one of the AIDA implementations should be
>> 77 available.
>> 78
>> 79 The code should be compiled with gmake and run with :
75 80
76 ./nanobeam << 81 > $G4WORDIR/bin/$G4SYSTEM/Nanobeam
77 82
78 The macro file default.mac is read by default. 83 The macro file default.mac is read by default.
79 84
80 Several macro files are provided: << 85 It can be one of the following macro files :
81 86
82 1) for the computation of intrinsic aberration 87 1) for the computation of intrinsic aberration coefficients :
83 coef-square.mac : using square magnetic field << 88 coef-square.mac : using square magnetic field model (=default.mac)
84 coef-map.mac : using 3D map magnetic field mod 89 coef-map.mac : using 3D map magnetic field model
85 coef-enge.mac : using Enge's analytical field 90 coef-enge.mac : using Enge's analytical field model
86 91
87 2) for the simulation of the beam image with a 92 2) for the simulation of the beam image with a realistic emittance :
88 image-square.mac : using square magnetic field << 93 image-square.mac : using square magnetic field model
89 image-map.mac : using 3D map magnetic field mo 94 image-map.mac : using 3D map magnetic field model
90 image-enge.mac : using Enge's analytical field 95 image-enge.mac : using Enge's analytical field model
91 96
92 3) for the simulation of grid shadow images 97 3) for the simulation of grid shadow images
93 grid-square.mac : using square magnetic field 98 grid-square.mac : using square magnetic field model
94 grid-map.mac : using 3D map magnetic field mod 99 grid-map.mac : using 3D map magnetic field model
95 grid-enge.mac : using Enge's analytical field 100 grid-enge.mac : using Enge's analytical field model
96 101
97 These macros files are stored in the ./macros 102 These macros files are stored in the ./macros directory.
98 103
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 104 ---->5. PHYSICS
107 105
108 The example runs with protons with fluctuating 106 The example runs with protons with fluctuating energies around 3 MeV.
109 Standard electromagnetic processes are activat 107 Standard electromagnetic processes are activated by default (corresponding to the
110 Physics builder G4EmStandardPhysics), includin 108 Physics builder G4EmStandardPhysics), including the G4StepLimiter process.
111 109
112 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS 110 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS
113 111
114 All results are stored in the nanobeam.root fi 112 All results are stored in the nanobeam.root file and can be displayed with the provided
115 ROOT macro file plot.C: 113 ROOT macro file plot.C:
116 * be sure to have ROOT installed on your machi 114 * be sure to have ROOT installed on your machine
117 * be sure to be in the directory where ROOT ou << 115 * be sure to be in the nanobeam directory
118 * copy plot.C into this directory << 116 * launch ROOT by typing root
119 * launch ROOT by typing root, then under your << 117 * under your ROOT session, type in : .X plot.C to execute the macro file
120 to execute the macro file <<
121 * or type directly: root plot.X <<
122 118
123 This macro file shows : 119 This macro file shows :
124 - the beam profile along the nanobeam line (on 120 - the beam profile along the nanobeam line (only for the computation of intrinsic
125 coefficients) 121 coefficients)
126 - the beam image (Y vs X) on target 122 - the beam image (Y vs X) on target
127 - the beam emittance (THETA vs X) and (PHY vs 123 - the beam emittance (THETA vs X) and (PHY vs Y) on target
128 - the grid shadow image (option) 124 - the grid shadow image (option)
129 <<
130 The output ntuples can be written as xml or cs <<
131 125
132 ---------------------------------------------- 126 ---------------------------------------------------------------------------
133 127
134 Should you have any enquiry, please do not hes 128 Should you have any enquiry, please do not hesitate to contact:
135 incerti@cenbg.in2p3.fr 129 incerti@cenbg.in2p3.fr