Solid EC Digitization

From Hall A Wiki
Jump to: navigation, search

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