Solid EC Digitization
From Hall A Wiki
From Energy deposition to number of photons are implemented in Jin Huang's standalone code and Rakitha's remoll code similarly
https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/CaloSimShashlik/include/Config.h https://github.com/JeffersonLab/remoll_solid/blob/master/src/remollCalDetector.cc
devel version for GEMC EC hit processing is here
https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/source/2.3/hitprocess/solid_ec_hitprocess.cc add Birk effect and decay in WLS fiber for EC in GEMC
refer to example from CLAS12 code
https://github.com/gemc/source/blob/master/hitprocess/clas12/cnd_hitprocess.cc https://github.com/gemc/source/blob/master/hitprocess/bdx/cormo_hitprocess.cc
devel version for GEMC detector is here
https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/ec_segmented_solid_gemc/ note the bank definition have to match hit processing, the time window in hit definition can affect number of hits output and running speed
- EC detector in GEMC is implemented this way:
- the parameter file define all module id and location, and each's module's characteristics. then entire EC is built from the parameter file by perl scripts
- define one preshower module, then make copy to produce all the rest preshower modules
- define one shower module, then make copy to produce all the rest modules
- for one shower module, define full shower as scintillation with one id as mother, then define one lead layer,one gap front and one gap back, then make copy of lead and gap layer all as daughter.
- EC hit processing is implemented this way:
- in output root file, each entry in the solid_ec vector tree is one hit
- each hit is a sum of all steps from all tracks in all scintillation layers of a module within one time windows (100ns) because all scintillation layers for one module has one id
- for each hit,
- totEdep sum all steps deposited energy,
- every step Edep is corrected by Birk effect to EdepB, then sum into totEdepB
- every step EdepB is corrected by light decay in shower WLS fiber to Eend, then sum into totEend which give the final result
- Edep_seg0-9,EdepB_seg0-9,Eend_seg0-9 are recorded along shower module for debugging purpose. It should give EM show profile.
- any convertiion photon is left for digitization or analysis after simulation to avoid unnecessary rerun simulation
- currently preshower share same hit processing with shower, so the Birk effect is ok, but light decay is wrong, we need to separate its hit processing from shower
refer to GEMC hit definition https://gemc.jlab.org/gemc/html/documentation/hitDefinition.html