Warning, /geant4/examples/advanced/xray_TESdetector/README.txt is written in an unsupported language. File is not indexed.
0001 --------------------------------------------------------------------------------
0002
0003 =========================================================
0004 Geant4 - an Object-Oriented Toolkit for Simulation in HEP
0005 =========================================================
0006
0007 Xray_TESdetector
0008 ------------------
0009 P.Dondero (1), R.Stanzani (1)
0010 Dec 2022
0011
0012 1. Swhard S.r.l, Genoa (GE), Italy.
0013
0014 --------------------------------------------------------------------------------
0015
0016 Contacts: paolo.dondero@cern.ch, ronny.stanzani@cern.ch
0017
0018 --------------------------------------------------------------------------------
0019 Acknowledgements: example developed within the ESA AREMBES Project,
0020 Contract n. 4000116655/16/NL/BW. This example is a reduced mass model of the
0021 Athena X-IFU instrument based on an early configuration which is no more
0022 applicable for any evaluation. Simone Lotti provided the simplified mass
0023 model and background derived from those used in [1].
0024 --------------------------------------------------------------------------------
0025
0026 Xray_TESdetector is an example of the application of Geant4 in a space
0027 environment. It represents an x-ray detector derived from the X-IFU, the
0028 X-ray spectrometer designed and developed by the European Space Agency (ESA)
0029 for use on the ATHENA telescope.
0030 The detector is a Transition-edge sensor (TES) composed of 317 Bismuth pixels
0031 arranged in a hexagonal shape and its setup includes different layers of
0032 shielding, filters and support.
0033 The main purpose of the simulation is the estimation of the particle radiation
0034 background impacting on the detector. For execution time optimization purposes,
0035 only particle steps respecting specific conditions (e.g. hit selected
0036 volumes close to the detector) are stored on a .root file [2].
0037 An example of ROOT-based analysis of the output file is included
0038 ("./analysis/analysis.C") and can be used to obtain basic plots and histograms.
0039 Xray_TESdetector implements a physics list dedicated to space radiation interactions,
0040 developed within the ESA AREMBES Project for the ATHENA mission, called Space
0041 Physics List (SPL).
0042 Technically, this example shows how to manage a complex geometry obtained
0043 with advanced detector construction features (e.g., boolean operations,
0044 parameterisation).
0045 In addition, the example shows a way to optimize the simulation's execution time
0046 and output size by selectively saving data based on specific combined conditions
0047 (e.g. position, eventID and process name).
0048 NOTE: in a multiple-run session, the last run always overrides the root file.
0049
0050 1 - GEOMETRY DEFINITION
0051
0052 The geometry consists of a simplified version of the X-IFU detector and is composed of
0053 the following:
0054 - the TES array, the backscattering (BSC) and the
0055 Anti-coincidence detector (ACD);
0056 - the structural elements supporting the detector (e.g. the cage underneath
0057 it);
0058 - the thermal shieldings;
0059 - the structural elements of the cryostat chamber;
0060 - a hollow Aluminum sphere schematizing the satellite.
0061 Detector parameters:
0062 - Detector thickness: 3 um
0063 - Number of pixels: 317
0064 - Detector's shape: regular hexagon
0065 - Hexagon's apothem: 8.593 mm
0066 The default geometry is constructed in DetectorConstruction class.
0067 Alternatively, a GDML file is provided (xray_TESdetector.gdml).
0068 The position of each pixel is defined by a list of coordinates (x,y)
0069 contained in "pixelpos.txt".
0070
0071 2 - PARTICLE SOURCE
0072
0073 The radiation field is composed of galactic cosmic rays (GCR) protons with a
0074 flux estimated for the L1/L2 Lagrangian points, as described in [1].
0075 The energies range from 10 MeV to 100 GeV, and the particles are isotropically
0076 generated on the surface of a sphere surrounding the geometry and randomly
0077 launched toward its interior. The detector is placed in the center of the
0078 sphere and the sphere's radius is chosen to avoid intersections with geometry
0079 elements.
0080
0081 3 - PHYSICS LIST
0082
0083 This example implements a dedicated physics list called "Space Physics List",
0084 developed within the ESA AREMBES Project. This physics list has been designed
0085 focusing on the ATHENA physics processes, but contains high precision
0086 models that can be used in a more general space application.
0087 In details, this physics list provides a custom electromagnetic part combined
0088 with the QBBC hadronic physics list.
0089 In volumes near the detector, where high precision in the scattering description
0090 is needed, the use of Single Scattering (SS) model is reccomended, as shown in
0091 the "run01.mac", through the SetEmRegion command.
0092 The use of SS only in selected regions allows the simulation to reduce CPU
0093 consumption in the majority of the volumes and be very accurate near the
0094 detector.
0095 The default production cuts are selected for all volumes, i.e. 1mm.
0096
0097 4 - HOW TO RUN THE EXAMPLE
0098
0099 Compile code and execute the example in 'batch' mode from the macro file:
0100 ./XrayTESdetector run01.mac
0101 to launch it with the DetectorConstruction, or:
0102 ./XrayTESdetector run02.mac
0103 to launch it by using the provided GDML.
0104 For this example, the multi-thread (MT) capability of Geant4 is enabled by
0105 default. To specify the desired number of threads, the user can use the
0106 command "/run/numberOfThreads" in "run01.mac".
0107
0108 5 - STEPPING
0109
0110 Within the "SteppingAction" class relevant information about the particle's
0111 state are stored in Tuples [2], defined in the "HistoManager" class.
0112 The tuples contain the following information:
0113 1. event ID
0114 2. volume name
0115 3. track ID
0116 4. coordinates (x,y,z)
0117 5. angles (theta, phi)
0118 6. parent ID
0119 7. pixel number (from the TES array)
0120 8. step energy deposit
0121 9. step number
0122 10. initial kinetic energy
0123 11. kinetic energy
0124 12. particle name
0125 13. pre and post-step names
0126 14. creator process name
0127
0128 Tuples are filled with the informations listed above in two cases:
0129 - when a new particle is generated (both primaries and secondaries);
0130 - when the particles reach the volumes next to the detector and the
0131 detector itself.
0132
0133 6 - ANALYSIS
0134
0135 xray_TESdetector provides an analysis macro example (analysis.C) with several
0136 predefined histograms:
0137 - Average energy deposit per pixel (1D);
0138 - Energy deposit on the detector (2D);
0139 - Particle count per pixel (1D);
0140 - Spectra of the primaries on the detector (1D);
0141 - Total spectra on the detector (1D);
0142 - distribution of the particles on the detector (1D).
0143
0144 The first three are used to qualitatively check how the interactions are
0145 distributed on the detector pixels and what is the average deposit per pixel
0146 and particle. The 2D histogram for the Energy deposit on the detector shows
0147 the shape of the detector on the XY plane.
0148 The spectrum histograms are used to observe the following:
0149 - the initial energies of the particles (at launch or generation);
0150 - the energy deposit on the detector;
0151 - the energy of the step before the impact on the pixel.
0152 Those information are the starting point to assess the background
0153 composition and intensity on the detector, and thus optimize the
0154 detector shielding and background rejection techniques.
0155 Histograms are managed by the "analysis.C" file.
0156
0157 7 - VISUALISATION
0158
0159 The visualization manager is set via the G4VisExecutive class
0160 in the main() function in xray_TESdetector.cc.
0161 The initialisation of the drawing is done via a set of /vis/ commands
0162 in the macro vis.mac. This macro is automatically read from
0163 the main function when the example is used in interactive running mode.
0164
0165 --------------------------------------------------------------------------------
0166
0167 References
0168
0169 [1] S. Lotti, S. Molendi, C. Macculi, V. Fioretti, L. Piro et al., "Review of
0170 the Particle Background of the Athena X-IFU Instrument", The Astrophysical
0171 Journal, 2021.
0172 [2] BRUN, René, et al. "The ROOT Users Guide". CERN, http://root.cern.ch, 2003.