WG4 Modelling of Physics Processes and Software Tools

Contact Conveners: You can contact the conveners by composing their email addresses as firstname.lastname@cern.ch

• Ozkan Sahin
• Rob Veenhof
• Piet Verwilligen

Main Research Topics

Disclaimer: here below you can find a selected list of papers and talks. This list is by definition incomplete and items are not ordered neither chronologically neither according to importance.

Integration, Maintenance & Development of Simulation Tools

Simulation of MPGD detectors (Examples)

Simulation of charging-up & Space Charge

Other Items

Lectures & Training

Open Issues

We are looking for interested and motivated people that can improve the simulation tools. There are several items where we have ideas on how to improve the simulation codes, but where we lack manpower:

• Discrepancy in Gain between data and simulation in GEM-structures: Simulation underestimates with a factor 2(? check) the gain in a single GEM-foil, and this propagates also to multi-GEM structures. More information can be found [[][here (insert link]]. A possible candidate for this under-estimation of the gain is the fact that the Electric field changes substantially within the volume of a mesh element. Garfield(++) however takes the central value and asks Magboltz to calculate the Mean Free Path before collision. If the electric field varies strongly, the mean free path is over-estimated and the multiplication factor (e-/cm) is under-estimated. A possible solution continue
• Reading COMSOL fieldmap in Garfield++: Garfield++ has only a function to read a COMSOL fieldmap made of first-order elements (tetraeders with 4 nodes), while more accurate results could be obtained with second order elements (tetraeders with 10 nodes), because these second order elements allow a more precise meshing of volumes with convex or concave structures. Instead the fortran version of Garfield is able to read 2nd order element fieldmaps.
• Update Magboltz in Garfield++: Garfield++ is interfaced with Magboltz 10.1 (%REDcheck), while the latest version is 11.6 (Nov 2018).
• Update HEED in Garfield++: the version of HEED actually implemented in Garfield++ does not simulate the Coulomb scattering when passing through a material. Furthermore the specific primary ionisation is calculated based upon the initial kinetic energy of the particle and is not updated if the particle looses energy inside the detector. While this works fine for a Minimum Ionizing Particle passing through a thin absorber, this does not work for layered detectors where some material needs to be perforated or for the simulation of sub-relativistic electrons that loose energy at much higher rate than a MIP and might get stopped in the detector.
• ...