| 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 1) You must have compiled your Geant4 installation with the FULL version of the
72 CLHEP library which can handle matrix operatio 75 CLHEP library which can handle matrix operations.
73 76
74 2) The code should be compiled cmake and run w << 77 2) The variable G4ANALYSIS_USE must be set to 1.
>> 78
>> 79 3) In order to generate histograms, at least one of the AIDA implementations should be
>> 80 available.
>> 81
>> 82 4) The code should be compiled with gmake and run with :
75 83
76 ./nanobeam << 84 > $G4WORDIR/bin/$G4SYSTEM/nanobeam
77 85
78 The macro file default.mac is read by default. 86 The macro file default.mac is read by default.
79 87
80 Several macro files are provided: << 88 It can be one of the following macro files :
81 89
82 1) for the computation of intrinsic aberration 90 1) for the computation of intrinsic aberration coefficients :
83 coef-square.mac : using square magnetic field << 91 coef-square.mac : using square magnetic field model (=default.mac)
84 coef-map.mac : using 3D map magnetic field mod 92 coef-map.mac : using 3D map magnetic field model
85 coef-enge.mac : using Enge's analytical field 93 coef-enge.mac : using Enge's analytical field model
86 94
87 2) for the simulation of the beam image with a 95 2) for the simulation of the beam image with a realistic emittance :
88 image-square.mac : using square magnetic field << 96 image-square.mac : using square magnetic field model
89 image-map.mac : using 3D map magnetic field mo 97 image-map.mac : using 3D map magnetic field model
90 image-enge.mac : using Enge's analytical field 98 image-enge.mac : using Enge's analytical field model
91 99
92 3) for the simulation of grid shadow images 100 3) for the simulation of grid shadow images
93 grid-square.mac : using square magnetic field 101 grid-square.mac : using square magnetic field model
94 grid-map.mac : using 3D map magnetic field mod 102 grid-map.mac : using 3D map magnetic field model
95 grid-enge.mac : using Enge's analytical field 103 grid-enge.mac : using Enge's analytical field model
96 104
97 These macros files are stored in the ./macros 105 These macros files are stored in the ./macros directory.
98 106
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 107 ---->5. PHYSICS
107 108
108 The example runs with protons with fluctuating 109 The example runs with protons with fluctuating energies around 3 MeV.
109 Standard electromagnetic processes are activat 110 Standard electromagnetic processes are activated by default (corresponding to the
110 Physics builder G4EmStandardPhysics), includin 111 Physics builder G4EmStandardPhysics), including the G4StepLimiter process.
111 112
112 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS 113 ---->6. SIMULATION OUTPUT AND RESULT ANALYZIS
113 114
114 All results are stored in the nanobeam.root fi 115 All results are stored in the nanobeam.root file and can be displayed with the provided
115 ROOT macro file plot.C: 116 ROOT macro file plot.C:
116 * be sure to have ROOT installed on your machi 117 * be sure to have ROOT installed on your machine
117 * be sure to be in the directory where ROOT ou << 118 * be sure to be in the nanobeam directory
118 * copy plot.C into this directory << 119 * launch ROOT by typing root
119 * launch ROOT by typing root, then under your << 120 * 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 121
123 This macro file shows : 122 This macro file shows :
124 - the beam profile along the nanobeam line (on 123 - the beam profile along the nanobeam line (only for the computation of intrinsic
125 coefficients) 124 coefficients)
126 - the beam image (Y vs X) on target 125 - the beam image (Y vs X) on target
127 - the beam emittance (THETA vs X) and (PHY vs 126 - the beam emittance (THETA vs X) and (PHY vs Y) on target
128 - the grid shadow image (option) 127 - the grid shadow image (option)
129 <<
130 The output ntuples can be written as xml or cs <<
131 128
132 ---------------------------------------------- 129 ---------------------------------------------------------------------------
133 130
134 Should you have any enquiry, please do not hes 131 Should you have any enquiry, please do not hesitate to contact:
135 incerti@cenbg.in2p3.fr 132 incerti@cenbg.in2p3.fr