hadronic/Hadr04/README.md

0001 \page ExampleHadr04 Example Hadr04
0002
0003    This example is focused on neutronHP physics, especially neutron transport,
0004    including thermal scattering.
0005    See A.R. Garcia, E. Mendoza, D. Cano-Ott presentation at G4 Hadronic group
0006    meeting (04/2013) and note on G4NeutronHP package
0007
0008 ## MATERIALS AND GEOMETRY DEFINITION
0009
0010    It is a single box representing a 'quasi infinite' homogeneous medium.
0011    Two parameters define the geometry :
0012         - the material of the box,
0013         - the (full) size of the box.
0014
0015    The default geometry (1 m3 of pressurized water) is built in
0016    DetectorConstruction, but the above parameters can be changed interactively
0017    via commands defined in DetectorMessenger.
0018
0019    A function, and its associated UI command, allows to build a material
0020    directly from a single isotope.
0021
0022    To be identified by the ThermalScattering module, the elements composing a
0023    material must have a specific name (see G4NeutronHPThermalScatteringNames.cc)
0024    Examples of such materials are build in DetectorConstruction.
0025
0026 ## PHYSICS LIST
0027
0028    Only processes of neutronHP package are registered : neutronElastic
0029    (including thermalScattering), neutronInelastic, nCapture; nFission.
0030    See class NeutronHPphysics. No other hadronic nor electromagnetic processes.
0031
0032    A command allows to select or not ThermalScattering model.
0033
0034    Several hadronic physics options are controlled by environment variables.
0035    To select them, see Hadr04.cc
0036
0037    NB. class NeutronHPphysics can be reused with other physicsConstructors,
0038    as neutron processes are deleted before to be re-created.
0039
0040 ## AN EVENT : THE PRIMARY GENERATOR
0041
0042    The primary kinematic is a single particle randomly shooted at the
0043    centre of the box. The type of the particle and its energy are set in
0044    PrimaryGeneratorAction (neutron 2 MeV), and can be changed via the G4
0045    build-in commands of ParticleGun class (see the macros provided with
0046    this example).
0047
0048 ## PHYSICS
0049
0050    All secondaries are killed in StackingAction. Therefore an event consists of
0051    the transport of the primary neutron. Then one survey the thermal and non
0052    thermal part of this parcours.
0053
0054
0055 ## HISTOGRAMS
0056
0057    The test contains 7 built-in 1D histograms, which are managed by
0058    G4AnalysisManager and its Messenger. The histos can be individually
0059    activated with the command :
0060    /analysis/h1/set id nbBins  valMin valMax unit
0061    where unit is the desired unit for the histo (MeV or keV, etc..)
0062    (see the macros xxxx.mac).
0063    ```
0064      1   "incident neutron: nb of collisions above 1 eV"
0065      2   "incident neutron: total track length above 1 eV"
0066      3   "incident neutron: time of flight above 1 eV"
0067      4   "incident neutron: nb of collisions below 1 eV"
0068      5   "incident neutron: total track length below 1*eV"
0069      6   "incident neutron: time of flight below 1 eV"
0070      7   "incident neutron: energy distribution below 1*eV"
0071    ```
0072    The histograms are managed by the HistoManager class and its Messenger.
0073    The histos can be individually activated with the command :
0074    ```
0075    /analysis/h1/set id nbBins  valMin valMax unit
0076    ```
0077    where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..)
0078
0079    One can control the name of the histograms file with the command:
0080    ```
0081    /analysis/setFileName  name  (default Hadr04)
0082    ```
0083    It is possible to choose the format of the histogram file : root (default),
0084    xml, csv, by using namespace in HistoManager.hh
0085
0086    It is also possible to print selected histograms on an ascii file:
0087    ```
0088    /analysis/h1/setAscii id
0089 ```
0090    All selected histos will be written on a file name.ascii (default Hadr04)
0091
0092 ## VISUALIZATION
0093
0094    The Visualization Manager is set in the main().
0095    The initialisation of the drawing is done via the commands
0096    /vis/... in the macro vis.mac. To get visualisation:
0097 ```
0098 > /control/execute vis.mac
0099 ```
0100
0101    The detector has a default view which is a longitudinal view of the box.
0102    The tracks are drawn at the end of event, and erased at the end of run.
0103
0104 ## HOW TO START ?
0105
0106    Execute Hadr04 in 'batch' mode from macro files :
0107 ```
0108 % ./Hadr04  run01.mac
0109 ```
0110
0111    Execute Hadr04 in 'interactive mode' with visualization :
0112 ```
0113 % ./Hadr04
0114 Idle> control/execute vis.mac
0115  ....
0116  Idle> type your commands
0117  ....
0118  Idle> exit
0119 ```
0120
0121  Macros provided in this example:
0122   - hadr04.in: macro used in Geant4 testing
0123   - graphite.mac: neutron (2 MeV) in graphite
0124   - run01.mac: neutron (2 MeV) in Water_ts
0125
0126  Macros to be run interactively:
0127   - debug.mac: neutron (2 MeV) in Water_ts
0128   - vis.mac: To activate visualization