Skip to content

Draft: Add a coalescence afterburner

Hendrik Jage requested to merge hjage_coalescence into Sim10

This MR aims to add a coalescence afterburner to Gauss. Currently the production of deuteron, triton, hypertriton, helium-3 and helium-4 is supported. It includes the following changes:

  • Add new CoalescenceProduction afterburner: Modifies HepMC::GenEvent generated by a generator in an afterburner step and combines nucleons close in phase-space into nuclei. The structure of the code is heavily based on the similar, already existing DeuteronProduction algorithm.
  • Modification of the CRMC options to allow calling the coalescence afterburner step
  • Modification of the CRMC options, to pass on full event cuts from the DecFile
  • Identify particles produced via coalescence with a new dedicated MCVertex type, which is introduced in LHCb!4838

This is still a draft, as the coalescence algorithm still needs to be reviewed internally.

Further TODOs are:

  • Figure out a way to select events with coalescence (MCVertex::Coalescence) in a full event cut of a decfile. Only selecting the pdg id of the produced ions is not enough, as ions can e.g. also be produced via spallation.
  • Add possibility to pass on afterburner options through CRMC

Merge request reports