Difference between revisions of "Solid EC Digitization"
From Hall A Wiki
(Link to new wiki) |
|||
(One intermediate revision by one other user not shown) | |||
Line 1: | Line 1: | ||
+ | This page has been moved to [https://solid.jlab.org/wiki/index.php/EC_Digitization https://solid.jlab.org/wiki/index.php/EC_Digitization] in the new [https://solid.jlab.org/wiki SoLID Wiki]. Please do not edit this page. | ||
+ | |||
+ | 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 | ||
+ | |||
+ | Birk effect | ||
+ | https://en.wikipedia.org/wiki/Birks%27_law | ||
+ | http://hypernews.slac.stanford.edu/HyperNews/geant4/get/emprocess/1272.html?inline=1 | ||
+ | |||
From Energy deposition to number of photons are implemented in Jin Huang's standalone code and Rakitha's remoll code similarly | 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://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/CaloSimShashlik/include/Config.h | ||
Line 7: | Line 17: | ||
https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/source/2.3/hitprocess/solid_ec_ps_hitprocess.cc | https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/source/2.3/hitprocess/solid_ec_ps_hitprocess.cc | ||
add Birk effect and decay in WLS fiber for ec shower and preshower in GEMC | add Birk effect and decay in WLS fiber for ec shower and preshower in GEMC | ||
− | |||
− | |||
− | |||
− | |||
devel version for GEMC detector is here | devel version for GEMC detector is here |
Latest revision as of 15:39, 7 December 2020
This page has been moved to https://solid.jlab.org/wiki/index.php/EC_Digitization in the new SoLID Wiki. Please do not edit this page.
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
Birk effect
https://en.wikipedia.org/wiki/Birks%27_law http://hypernews.slac.stanford.edu/HyperNews/geant4/get/emprocess/1272.html?inline=1
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 https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec/source/2.3/hitprocess/solid_ec_ps_hitprocess.cc add Birk effect and decay in WLS fiber for ec shower and preshower in GEMC
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
- for shower,Edep_seg0-9,EdepB_seg0-9,Eend_seg0-9 are recorded along shower module for debugging purpose. It should give EM show profile.
- any conversion to photon is left for digitization or analysis after simulation to avoid unnecessary rerunning simulation
- refer to GEMC hit definition https://gemc.jlab.org/gemc/html/documentation/hitDefinition.html
how to test on ifarm
source /group/solid/solid_svn/set_solid 1.3 cd $HOME (or somewhere not at /work/halla/solid) svn co https://jlabsvn.jlab.org/svnroot/solid/subsystem/ec cd ec/source/2.3 scons OPT=1 -j4 cd ../../analysis ../source/2.3/solid_gemc solid_PVDIS_LD2_ec_segmented.gcard