RealisticCalo* is a re-orgainsation of the calorimeter digitisation and reconstruction code, which was previously in ILDCaloDigi.
It's also designed to work only with DD4hep geometry information.
The task has been split into several tasks:
1) hit digitisation
conversion of geant4 energy deposit to some calibrated electronics signal.
e.g. for silicon, converted to the MIP scale, for scintillator to the SiPM pixel scale.
effects such as miscalibration, dead cells, pe statistics are applied at this stage
2) hit reconstruction
the conversion of these hits to the EM energy scale.
This includes unfolding of SiPM response, and correction for the calorimeter sampling ratio.
** note that the calibration factors given here are for converting from the MIP scale to the total deposited energy scale.
** this is different to LDC/ILDCaloDigi, where the conversion was from the energy deposited in the active material to total deposited energy.
3) gap correction
If hits are situated directly across a dead area of the detector, their energies and the distance between them is used to guess how much energy was lost.
"Gap correction hits" are now written into a new calorimeterhit collection (previously, the energy of existing hits at wafer boundaries was increased).
For now, a BruteForceEcalGapFiller has been written, which uses the distance between hits to decide if they are situated across a dead region.
This is close to what was done previously in e.g. ILDCaloDigi.
In principle, it should be possible (and probably more efficient) to use cell indices to decide whether a hit is at the edge of a wafer.
----------------------
One instance of each processor should be used for each collection of hits (or set of collections with identical properties).
There are virtual base classes RealisticCaloDigi and RealisticCaloReco.
Derived classes for a particular technology inherit from these: RealisticCaloDigiSilicon, RealisticCaloDigiScinPpd and RealisticCaloRecoSilicon, RealisticCaloRecoSilicon.
example of steering:
EcalBarrelCollection
EcalBarrelCollectionDigi
EcalBarrelCollectionDigiRelation
0.5
MIP
1
1
-10
100
1.7e-4
0.1
2000
layer
3.6
EcalBarrelCollectionDigi
EcalBarrelCollectionDigiRelation
EcalBarrelCollectionReco
EcalBarrelCollectionRecoRelation
20 100
6.58e-3 13.16e-3
1
layer
HcalBarrelRegCollection
HcalBarrelRegCollectionDigi
HcalBarrelRegCollectionDigiRelation
0.5
MIP
1
1
-10
100
5.3e-4
0.1
500
layer
15
1600
0
0
HcalBarrelRegCollectionDigi
HcalBarrelRegCollectionDigiRelation
HcalBarrelRegCollectionReco
HcalBarrelRegCollectionRecoRelation
100
3.0e-2
15
1600
layer
EcalBarrelCollectionReco
EcalBarrelCollectionGapHits
layer
module
stave