Commit fbe383c5 authored by Valentina Fioretti's avatar Valentina Fioretti
Browse files

Edit examples.md

parent 6a2043a1
Loading
Loading
Loading
Loading
+168 −1
Original line number Diff line number Diff line
@@ -113,7 +113,6 @@ Volume selection:
`IOXYZ.volumeselection.0.end = 1000000000`



Number of rows for each file:

`IOXYZ.NRows  = 5000000`
@@ -128,6 +127,8 @@ Exclude the world (NB: the world contains the initial energy and momentum):

----> Physics list:

Name of the physics class:

`PHYS.VERSION = OPTPhys`

Production cut applied to all volumes in the simulation:
@@ -218,6 +219,172 @@ and Single Scattering physics lists. Several custom parameters,
including the rotation of the target, can be adapted at runtime
using the configuration file.

The example selects the Geant4 reference single scattering electromagnetic physics list (emstandardSS_QBBC) but the [AREMBESPhys](https://www.ict.inaf.it/gitlab/icsc_g4_hpc/BoGEMMS-HPC/-/blob/main/phys/AREMBESPhys.cc?ref_type=heads) general physics list allows the user to select other options, for which custom user limits could be required.

To run the example: 

`> bogemms -c ProtonScattering.conf -m sp_300keV.mac`

###### Configuration file ProtonScattering.conf:

Selection of output files:

`IO.FILETYPE.FITS = true`

`IO.FILETYPE.FITS.COMPRESS = true`

`IO.FILE.XYZ = true`


IO XYZ file configuration:

`IOXYZ.EVT_ID = true`

`IOXYZ.TRK_ID = true`

`IOXYZ.PARENT_TRK_ID = true`

`IOXYZ.VOLUME_ID = true`

`IOXYZ.VOLUME_NAME = true`

`IOXYZ.MOTHER_ID = true`

`IOXYZ.E_DEP = true`

`IOXYZ.X_ENT = true`

`IOXYZ.Y_ENT = true`

`IOXYZ.Z_ENT = true`

`IOXYZ.X_EXIT = true`

`IOXYZ.Y_EXIT = true`

`IOXYZ.Z_EXIT = true`

`IOXYZ.E_KIN_ENT = true`

`IOXYZ.E_KIN_EXIT = true`

`IOXYZ.MDX_ENT = true`

`IOXYZ.MDY_ENT = true`

`IOXYZ.MDZ_ENT = true`

`IOXYZ.MDX_EXIT = true`

`IOXYZ.MDY_EXIT = true`

`IOXYZ.MDZ_EXIT = true`

`IOXYZ.GTIME_ENT = true`

`IOXYZ.GTIME_EXIT = true`

`IOXYZ.PARTICLE_ID = true`

`IOXYZ.PARTICLE_NAME = false`

`IOXYZ.PROCESS_ID = true`

`IOXYZ.PROCESS_NAME = false`

Volume selection:

`IOXYZ.volumeselection.number = 1`

`IOXYZ.volumeselection.0.start = 0`

`IOXYZ.volumeselection.0.end = 10000000`


Number of rows for each file:

`IOXYZ.NRows  = 5000000`

Write only volume with energy deposit (this do not exclude the world):

`IOXYZ.writeonlyifenergydeposit = 0`

Exclude the world (NB: the world contains the initial energy and momentum):

`IOXYZ.removeWorldEvent = 0`

Volume ID handling (# 1 = summing mother and grandmother ID in case of replicas):

`GEOM.VOLID.TYPE = 1`

Multithreading:

`MT.NUM.THREADS = 1`

----> Physics list:

Name of the physics class:

`PHYS.VERSION = AREMBESPhys`

Options:
`AREMBESPhys.PHYSLIST = emstandardSS_QBBC`

`AREMBESPhys.DEFAULT.CUT = 0.001`


Activating SS EM parameters:

`AREMBESPhys.SSPHYS.ACTIVATE = 1`


User limit (activating the custom user limit to the target, to be used in multiple scattering is chosen instead of SS):

`AREMBESPhys.STEPLIMITER.ACTIVATE = 0`

Region cuts applied to the target:

`REGION.TARGET.CUT = 0.000001`

`REGION.TARGET.PROTON.CUT = 0.000001`

Maximum step limiter (if step limiter is activated):

`GEOM.MAXSTEP = 0.00000001`

----> Geometry configuration:

Setting the world dimension [mm]:

`WORLD.BOX.SIDE = 5000`

Name of the geometry class:

`GEOM.VERSION = ProtonScattering`

Experiment type 1 (target and detector)

`GEOM.EXPERIMENT = 1`

`GEOM.MAXSTEP = 0.1`

# selecting the chamber material
GEOM.VACUUM = 1

# collimator distance (from right side to target center) in mm
GEOM.COLLIMATOR.DISTANCE = 87.

##### TARGET
# target type: 1 = eRosita
GEOM.TARGET = 1

GEOM.INCIDENT.ANGLE = 0.64

# detector distance (distance from the target center to the detector surface) in mm
GEOM.PANEL.DISTANCE = 933.



</details>