Back to home page

EIC code displayed by LXR

 
 

    


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.