Difference between revisions of "SoLID full simulation and file sharing"

From Hall A Wiki
Jump to: navigation, search
(New page: = Introduction = We will run full simulation with all subsystems and generate various files at a shared central location They can used for studies like acceptance,trigger, background, GE...)
 
(Link to new wiki)
 
(26 intermediate revisions by one other user not shown)
Line 1: Line 1:
 +
This page has been moved to [https://solid.jlab.org/wiki/index.php/full_simulation_and_file_sharing https://solid.jlab.org/wiki/index.php/full_simulation_and_file_sharing] in the new [https://solid.jlab.org/wiki SoLID Wiki].  Please do not edit this page.
 +
 
= Introduction =
 
= Introduction =
  
Line 5: Line 7:
 
They can used for studies like acceptance,trigger, background, GEM etc and ensure consistent results
 
They can used for studies like acceptance,trigger, background, GEM etc and ensure consistent results
  
== Luminosity and radiation thickness ==
+
Many file size are large in order of 10-100GB, you can transfer them by using "scp jlabname@ftp.jlab.org:file" or globus here https://cc.jlab.org/filetransfers which is about twice faster
 +
 
 +
= luminosity and radiation thickness =
 +
 
 +
Z,A, and effective nuclei luminosity calculated here are used as input by eDIS generator at
 +
https://github.com/JeffersonLab/evgen_inclusive/output
 +
 
 +
PVDIS
 +
50uA, raster 0.5cm diameter
 +
LH target, 40cm, 0.071g/cm3, X=40/890.4=4.5e-2, Lumi=50e-6/1.6e-19*40*0.071*6.02e23=0.53e39/cm2/s
 +
LD target, 40cm, 0.169g/cm3, X=40/745.4=5.4e-2, Lumi=50e-6/1.6e-19*40*0.169*6.02e23=1.27e39/cm2/s
 +
Al window, 2*100um, density 2.7g/cm3, X=2*100e-4/8.897=2.25e-3, Lumi=50e-6/1.6e-19*2*100e-4*2.7*6.02e23=1e37/cm2/s
 +
 
 +
SIDIS He3
 +
15uA, raster 0.5cm diameter
 +
3he(10amg), 40cm, density=10*44.6(amg=mol/m3)*3.016(g/mol)=1.345e-3g/cm3, X=40/(67.42/1.345e-3)=0.8e-3, Lumi=15e-6/1.6e-19*40*1.345e-3*6.02e23=3e36/cm2/s, pol Lumi = 3e36/3 = 1e36cm2/s
 +
GE180 glass window, 2*120um, density 2.76g/cm3, X=2*120e-4/(19.4/2.76)=3.4e-3, Lumi=15e-6/1.6e-19*2*120e-4*2.76*6.02e23=3.74e36/cm2/s
 +
 
 +
GE180 Aluminosilicate Glass Composition,
 +
  refer to http://galileo.phys.virginia.edu/research/groups/spinphysics/glass_properties.html
 +
  Molecule Composition by weight
 +
      SiO2 60.3%
 +
      BaO 18.2%
 +
      Al2O3 14.3%
 +
      CaO 6.5%
 +
      SrO 0.25%
 +
  Z/A                   0.4829
 +
  radiation thickness   19.4246 g/cm2
 +
  radiation thickness   7.038 cm
 +
  density                2.76 g/cm3
 +
  assume Z=17 and A=35,which gives correct Z/A, then we got nuclei luminosity for each window 3.74/2/35=0.054e36/cm2/s
 +
 
 +
SIDIS proton
 +
 +
I = 100nA, raster 2.4cm diameter
 +
Target = NH3, density of NH3(den) = 0.819 g/cm^3
 +
target thickness (t) = 2.826 cm, packing fraction(PF) = 0.55
 +
NH3 target is immersed in liquid He4, density of LHe4 = 0.145 g/cm^3 at 1K
 +
Luminosity of NH3 = den * t * PF * (I/1.602e-19)*6.022e23 = 0.819*2.826*0.55*(100e-9/1.602e-19)*6.022e23 = 4.785e35/cm2/s
 +
pol Lumi = 4.785e35/17*3 = 0.844e35/cm2/s
 +
Luminosity of LHe4 = den * t * (1-PF) * (I/1.602e-19)*6.022e^23 = 0.145*2.826*0.45*(100e-9/1.602e-19)*6.022e23 = 0.69e35/cm2/s
 +
Luminosity of two LHe4 outside  =  0.145*2*0.432*(100e-9/1.602e-19)*6.022e23 = 0.47e35/cm2/s
 +
Total luminosity = (4.785+0.69+0.47) * 1e35 = 5.945e35/cm2/s
 +
Note: there are several thin Al windows before and after target, their total luminosity ~1e35/cm2/s
 +
 
 +
total length of NH3+LHe4 and two LHe4 outside, 2.826+2*0.432=3.69
 +
average density of NH3+LHe4, (0.819*0.55+0.145*0.45) = 0.5157
 +
average density of NH3+LHe4 and two LHe4 outside, (0.5157*2.826+0.145*2*0.432)/(2.826+2*0.432)=0.43
 +
Composition by weight of NH3+LHe4 and two LHe4 outside
 +
    NH3                  0.819*0.55*2.826/(2.826+2*0.432)/0.43=0.80
 +
    LHe4                  0.145*0.45*2.826/(2.826+2*0.432)/0.43=0.12
 +
    two LHe4 outside      0.145*2*0.432/(2.826+2*0.432)/0.43=0.08
 +
average Z=10*0.80+2*0.12+2*0.08=8.4
 +
average A=17*0.80+4*0.12+4*0.08=14.4
 +
Z/A=0.583
 +
assume Z=7 and A=12,which gives correct Z/A, then we got nuclei luminosity 5.945e35/12=0.495e35/cm2/s
 +
 
 +
reference:
 +
g2p thesis and target info https://hallaweb.jlab.org/wiki/index.php/G2p
 +
g2p target in simulation https://github.com/asymmetry/g2psim/blob/master/src/G2PTarget.cc
 +
 
 +
dilution factor for pure NH3 is 0.176=M_H3/(M_H3+M_N)
 +
Considering other materials, 0.176*4.785/5.945=0.142 (not including windows), 0.176*4.785/(5.945+1)=0.121 (including windows)
 +
g2p experiment has measured dilution 0.13 because particle from windows have different acceptance than target
 +
5% error of dilution factor https://userweb.jlab.org/~jones/tajima_rss/plots/2007/sep25/dfpf_oscar_report.pdf
 +
 
 +
JPsi
 +
3uA, raster 0.5cm diameter
 +
LH target, 15cm, 0.071g/cm3, X=15/890.4=1.7e-2, Lumi=3e-6/1.6e-19*15*0.071*6.02e23=1.2e37/cm2/s
 +
Al window, 2*100um, density 2.7g/cm3, X=2*100e-4/8.897=2.25e-3, Lumi=3e-6/1.6e-19*2*100e-4*2.7*6.02e23=6e35/cm2/s
  
== file sharing ==
+
= file sharing (current) =
  
The file size is large in order of 50GB, you can transfer them by using "scp jlabname@ftp.jlab.org:file" slowly or other faster ways like globus here https://cc.jlab.org/filetransfers
+
* geant4 crosssection study
 +
/mss/halla/solid/sim/solid_gemc/target_JLAB_VERSION_1.3
 +
/cache/halla/solid/sim/solid_gemc/target_JLAB_VERSION_1.3
 +
*BeamOnTargethadron*.root  (hadron only)
 +
*BeamOnTargetEM/*.root    (EM only)
 +
*BeamOnTarget/*.root      (EM+hadron)
  
* sharing location
+
* SoLID full simulation sharing location
  location on disk /cache/halla/solid/sim/solid_gemc/
+
latest as of 2020/07
 +
  /cache/halla/solid/sim/solid_gemc/PVDIS_LD2_JLAB_VERSION_1.3/pass4
 +
  /cache/halla/solid/sim/solid_gemc/SIDIS_He3_JLAB_VERSION_1.3/pass7
 +
  /cache/halla/solid/sim/solid_gemc/SIDIS_NH3_JLAB_VERSION_1.3/pass6
 +
  /cache/halla/solid/sim/solid_gemc/JPsi_LH2_JLAB_VERSION_1.3/pass5
 +
  location on disk /cache/halla/solid/sim/solid_gemc/ (The latest files are pinned on cache for max allowed 60 days. If you need older files, ask Zhiwen Zhao)
 
  location on tape /mss/halla/solid/sim/solid_gemc/
 
  location on tape /mss/halla/solid/sim/solid_gemc/
 
  dir
 
  dir
  PVDIS_LD2_JLAB_VERSION_1.2
 
 
   PVDIS_LD2_JLAB_VERSION_1.3
 
   PVDIS_LD2_JLAB_VERSION_1.3
   SIDIS_He3_JLAB_VERSION_1.2
+
   SIDIS_He3_JLAB_VERSION_1.3
 
   SIDIS_NH3_JLAB_VERSION_1.3
 
   SIDIS_NH3_JLAB_VERSION_1.3
   JPsi_LH2_JLAB_VERSION_1.2    
+
   JPsi_LH2_JLAB_VERSION_1.3    
 +
  ...
 
  subdir
 
  subdir
 
   pass1
 
   pass1
Line 26: Line 107:
 
   pass4
 
   pass4
 
   ...
 
   ...
  refer details
+
  detailed log are at
   /home/zwzhao/solid/solid_svn/study/background/log_PVDIS_LD2
+
   new and updated
   /home/zwzhao/solid/solid_svn/study/background/log_SIDIS_He3
+
  [https://github.com/JeffersonLab/solid_gemc/blob/master/analysis/log_PVDIS_LD2 log_PVDIS_LD2]
   /home/zwzhao/solid/solid_svn/study/background/log_SIDIS_NH3
+
   [https://github.com/JeffersonLab/solid_gemc/blob/master/analysis/log_SIDIS_He3 log_SIDIS_He3]
   /home/zwzhao/solid/solid_svn/study/background/log_JPsi_LH2
+
   [https://github.com/JeffersonLab/solid_gemc/blob/master/analysis/log_SIDIS_NH3 log_SIDIS_NH3]
 +
   [https://github.com/JeffersonLab/solid_gemc/blob/master/analysis/log_JPsi_LH2 log_JPsi_LH2]
  
 
* file name example
 
* file name example
Line 50: Line 132:
 
  see example here https://jlabsvn.jlab.org/svnroot/solid/study/background/background.C and https://jlabsvn.jlab.org/svnroot/solid/study/trigger/Get_PVDIS_LD2_trigger_rate.C
 
  see example here https://jlabsvn.jlab.org/svnroot/solid/study/background/background.C and https://jlabsvn.jlab.org/svnroot/solid/study/trigger/Get_PVDIS_LD2_trigger_rate.C
  
== old ==
+
= file sharing (old runs) =
  
=== run 2  ===  
+
== run 2  ==
 
It was done by Zhiwen Zhao in later 2014
 
It was done by Zhiwen Zhao in later 2014
  
Line 58: Line 140:
 
https://jlabsvn.jlab.org/svnroot/solid/study/background
 
https://jlabsvn.jlab.org/svnroot/solid/study/background
  
working dir with all output /work/halla/solid/sim/solid_gemc/*_geant4.9.6.p02 and /work/halla/solid/sim/solid_gemc/*_geant4.9.5.p01
+
location on disk /cache/halla/solid/sim/solid_gemc/  
 +
 
 +
location on tape /mss/halla/solid/sim/solid_gemc/
 +
 
 +
(ask Zhiwen Zhao if you need to use them)
  
 
table for normalization factor comparison  
 
table for normalization factor comparison  
 
[http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/normalization_factor_comparison.pdf pdf] [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/normalization_factor_comparison.pdf pptx]
 
[http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/normalization_factor_comparison.pdf pdf] [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/normalization_factor_comparison.pdf pptx]
 +
 +
file names:
 +
"EM" in file name means all real materials and beam on target which needs normalization by current/nevent (same as "EM" previously)
 +
"clean_weighted" in file name means kryptonite for all geometry and weighted generator (same as "other" previously)
 +
"dirty_normalized" in file name mean all real materials and normalized generator (same as "actual" previously)
  
 
more hit_id are added, some are changed
 
more hit_id are added, some are changed
Line 87: Line 178:
 
         10  other
 
         10  other
  
=== run 1 ===  
+
== run 1 ==  
 
It was done by Zhiwen Zhao in later 2013
 
It was done by Zhiwen Zhao in later 2013
  
Line 102: Line 193:
 
low energy neutron cross section was turned on
 
low energy neutron cross section was turned on
  
  file location for SIDIS_He3
+
(ask Zhiwen Zhao if you need to use them)
  /work/halla/solid/sim/solid_gemc/SIDIS_He3_run1
+
  file location
 +
on disk /cache/halla/solid/sim/solid_gemc/
 +
on tape  /mss/halla/solid/sim/solid_gemc/
 +
for SIDIS_He3, it's at subdir SIDIS_He3_run1/
 +
for PVDIS_LD2, it's at subdir PVDIS_LD2_run1/, main result is dir baffle_babarbafflemore1_block which has photon block before EC, while baffle_babarbafflemore1 has result without the photon block
  
 
  file names for SIDIS_He3
 
  file names for SIDIS_He3
Line 110: Line 205:
 
  "actual" in file name mean all real materials
 
  "actual" in file name mean all real materials
 
  "sum"  in file name mean sum over all real materials for gas and two windows
 
  "sum"  in file name mean sum over all real materials for gas and two windows
 
file location for PVDIS_LD2
 
/work/halla/solid/sim/solid_gemc/PVDIS_LD2_run1
 
main result is dir baffle_babarbafflemore1_block which has photon block before EC
 
while baffle_babarbafflemore1 has result without the photon block
 
  
 
  file names for PVDIS_LD2
 
  file names for PVDIS_LD2
  file with no special name is just shooting beam on target, "EM" particles dominating the results  
+
  file with no special name is just shooting beam on target, "EM" particles dominating the results, it has lead baffle and kryptonite for everything else
 
  "other" in file name means kryptonite for all geometry  
 
  "other" in file name means kryptonite for all geometry  
 
  "real" in file name mean lead baffle and kryptonite for everything else
 
  "real" in file name mean lead baffle and kryptonite for everything else

Latest revision as of 15:39, 7 December 2020

This page has been moved to https://solid.jlab.org/wiki/index.php/full_simulation_and_file_sharing in the new SoLID Wiki. Please do not edit this page.

Introduction

We will run full simulation with all subsystems and generate various files at a shared central location

They can used for studies like acceptance,trigger, background, GEM etc and ensure consistent results

Many file size are large in order of 10-100GB, you can transfer them by using "scp jlabname@ftp.jlab.org:file" or globus here https://cc.jlab.org/filetransfers which is about twice faster

luminosity and radiation thickness

Z,A, and effective nuclei luminosity calculated here are used as input by eDIS generator at https://github.com/JeffersonLab/evgen_inclusive/output

PVDIS

50uA, raster 0.5cm diameter
LH target, 40cm, 0.071g/cm3, X=40/890.4=4.5e-2, Lumi=50e-6/1.6e-19*40*0.071*6.02e23=0.53e39/cm2/s 
LD target, 40cm, 0.169g/cm3, X=40/745.4=5.4e-2, Lumi=50e-6/1.6e-19*40*0.169*6.02e23=1.27e39/cm2/s
Al window, 2*100um, density 2.7g/cm3, X=2*100e-4/8.897=2.25e-3, Lumi=50e-6/1.6e-19*2*100e-4*2.7*6.02e23=1e37/cm2/s

SIDIS He3

15uA, raster 0.5cm diameter
3he(10amg), 40cm, density=10*44.6(amg=mol/m3)*3.016(g/mol)=1.345e-3g/cm3, X=40/(67.42/1.345e-3)=0.8e-3, Lumi=15e-6/1.6e-19*40*1.345e-3*6.02e23=3e36/cm2/s, pol Lumi = 3e36/3 = 1e36cm2/s
GE180 glass window, 2*120um, density 2.76g/cm3, X=2*120e-4/(19.4/2.76)=3.4e-3, Lumi=15e-6/1.6e-19*2*120e-4*2.76*6.02e23=3.74e36/cm2/s
GE180 Aluminosilicate Glass Composition, 
  refer to http://galileo.phys.virginia.edu/research/groups/spinphysics/glass_properties.html
  Molecule 	Composition by weight
     SiO2	 	60.3%
     BaO 	 	18.2%
     Al2O3		14.3%
     CaO 	 	6.5%
     SrO 	 	0.25% 
  Z/A 	                   0.4829 
  radiation thickness 	   19.4246 g/cm2
  radiation thickness 	   7.038 cm
  density                 2.76 g/cm3
  assume Z=17 and A=35,which gives correct Z/A, then we got nuclei luminosity for each window 3.74/2/35=0.054e36/cm2/s

SIDIS proton

I = 100nA, raster 2.4cm diameter
Target = NH3, density of NH3(den) = 0.819 g/cm^3
target thickness (t) = 2.826 cm, packing fraction(PF) = 0.55
NH3 target is immersed in liquid He4, density of LHe4 = 0.145 g/cm^3 at 1K
Luminosity of NH3 = den * t * PF * (I/1.602e-19)*6.022e23 = 0.819*2.826*0.55*(100e-9/1.602e-19)*6.022e23 = 4.785e35/cm2/s
pol Lumi = 4.785e35/17*3 = 0.844e35/cm2/s
Luminosity of LHe4 = den * t * (1-PF) * (I/1.602e-19)*6.022e^23 = 0.145*2.826*0.45*(100e-9/1.602e-19)*6.022e23 = 0.69e35/cm2/s
Luminosity of two LHe4 outside   =  0.145*2*0.432*(100e-9/1.602e-19)*6.022e23 = 0.47e35/cm2/s
Total luminosity = (4.785+0.69+0.47) * 1e35 = 5.945e35/cm2/s
Note: there are several thin Al windows before and after target, their total luminosity ~1e35/cm2/s 
total length of NH3+LHe4 and two LHe4 outside, 2.826+2*0.432=3.69
average density of NH3+LHe4, (0.819*0.55+0.145*0.45) = 0.5157
average density of NH3+LHe4 and two LHe4 outside, (0.5157*2.826+0.145*2*0.432)/(2.826+2*0.432)=0.43
Composition by weight of NH3+LHe4 and two LHe4 outside
   NH3                   0.819*0.55*2.826/(2.826+2*0.432)/0.43=0.80
   LHe4                  0.145*0.45*2.826/(2.826+2*0.432)/0.43=0.12
   two LHe4 outside      0.145*2*0.432/(2.826+2*0.432)/0.43=0.08
average Z=10*0.80+2*0.12+2*0.08=8.4
average A=17*0.80+4*0.12+4*0.08=14.4
Z/A=0.583
assume Z=7 and A=12,which gives correct Z/A, then we got nuclei luminosity 5.945e35/12=0.495e35/cm2/s
reference: 
g2p thesis and target info https://hallaweb.jlab.org/wiki/index.php/G2p
g2p target in simulation https://github.com/asymmetry/g2psim/blob/master/src/G2PTarget.cc
dilution factor for pure NH3 is 0.176=M_H3/(M_H3+M_N)
Considering other materials, 0.176*4.785/5.945=0.142 (not including windows), 0.176*4.785/(5.945+1)=0.121 (including windows)
g2p experiment has measured dilution 0.13 because particle from windows have different acceptance than target
5% error of dilution factor https://userweb.jlab.org/~jones/tajima_rss/plots/2007/sep25/dfpf_oscar_report.pdf 

JPsi

3uA, raster 0.5cm diameter
LH target, 15cm, 0.071g/cm3, X=15/890.4=1.7e-2, Lumi=3e-6/1.6e-19*15*0.071*6.02e23=1.2e37/cm2/s
Al window, 2*100um, density 2.7g/cm3, X=2*100e-4/8.897=2.25e-3, Lumi=3e-6/1.6e-19*2*100e-4*2.7*6.02e23=6e35/cm2/s

file sharing (current)

  • geant4 crosssection study
/mss/halla/solid/sim/solid_gemc/target_JLAB_VERSION_1.3
/cache/halla/solid/sim/solid_gemc/target_JLAB_VERSION_1.3
*BeamOnTargethadron*.root  (hadron only)
*BeamOnTargetEM/*.root    (EM only)
*BeamOnTarget/*.root      (EM+hadron)
  • SoLID full simulation sharing location
latest as of 2020/07
  /cache/halla/solid/sim/solid_gemc/PVDIS_LD2_JLAB_VERSION_1.3/pass4
  /cache/halla/solid/sim/solid_gemc/SIDIS_He3_JLAB_VERSION_1.3/pass7
  /cache/halla/solid/sim/solid_gemc/SIDIS_NH3_JLAB_VERSION_1.3/pass6
  /cache/halla/solid/sim/solid_gemc/JPsi_LH2_JLAB_VERSION_1.3/pass5
location on disk /cache/halla/solid/sim/solid_gemc/ (The latest files are pinned on cache for max allowed 60 days. If you need older files, ask Zhiwen Zhao)
location on tape /mss/halla/solid/sim/solid_gemc/
dir
  PVDIS_LD2_JLAB_VERSION_1.3
  SIDIS_He3_JLAB_VERSION_1.3
  SIDIS_NH3_JLAB_VERSION_1.3
  JPsi_LH2_JLAB_VERSION_1.3   
  ...
subdir
  pass1
  pass2
  pass3
  pass4
  ...
detailed log are at
  new and updated
  log_PVDIS_LD2
  log_SIDIS_He3
  log_SIDIS_NH3
  log_JPsi_LH2 
  • file name example
"*BeamOnTarget_1e9_skim.root", created by shooting 1e9 electron on target, full SoLID simulation, only events has any entry in any detector is kept, physics list used is "QGSP_BERT_HP+STD+Optical", including hadron,EM and optical process
"*BeamOnTargetEM_1e9_skim.root", created by shooting 1e9 electron on target directly, only events has any entry in any detector is kept, physics list used is "STD+Optical",meaning no hardon, only EM and optical process
"*dirty_weighted_eDIS_filenum100_1e6.root,  *dirty_normalized_*_pi*HallD_filenum500_5e6.root,  *dirty_normalized_*_allnoeHallD_filenum500_5e6.root"
dirty means full SoLID simulation with physics list "QGSP_BERT_HP+STD+Optical"
weighted means the generator events distribute evenly in certain kinematic space,then have weight linked to crossection and rate
normalized means the generator events distribute according to crossection by probability and weight is a constant for each event
eDIS means electron from target by the DIS electron generator mode in "eicRate" code
pion*HallD means pion from target by the modified hallD generator
allnoeHallD means all particles from the modified hallD generator which has no electron
5e6 are number of event
filenum500 means sum of 500 small files of 1e4=5e6/500 events each
  • note about normalization factor
for BeamOnTarget skim file, (current in A)/1.6e-19/(number of event) = rate in Hz, for PVDIS_current=50uA,SIDIS_He3=15uA,SIDIS_NH3=100nA,JPsi=3uA, number of event is in file name, which is not number of tree entries in skim root file
for other file, (var8->at(0) in tree header)/filenum = rate in Hz, filenum is included in file name, this is because each small file is normalized separately
see example here https://jlabsvn.jlab.org/svnroot/solid/study/background/background.C and https://jlabsvn.jlab.org/svnroot/solid/study/trigger/Get_PVDIS_LD2_trigger_rate.C

file sharing (old runs)

run 2

It was done by Zhiwen Zhao in later 2014

The code and log files are in SVN at https://jlabsvn.jlab.org/svnroot/solid/study/background

location on disk /cache/halla/solid/sim/solid_gemc/

location on tape /mss/halla/solid/sim/solid_gemc/

(ask Zhiwen Zhao if you need to use them)

table for normalization factor comparison pdf pptx

file names:
"EM" in file name means all real materials and beam on target which needs normalization by current/nevent (same as "EM" previously)
"clean_weighted" in file name means kryptonite for all geometry and weighted generator (same as "other" previously)
"dirty_normalized" in file name mean all real materials and normalized generator (same as "actual" previously)

more hit_id are added, some are changed

 =========    hit_id and pid definition ==============
 hitid  =0 - 5 6 GEM planes, unused
         29 - 40  6 GEM plane, 1st layer (odd) and 2nd layer (even) of gas		
         20 - 25  6 GEM plane front			
         6,18    LGCC  PMT, front
         7,19    HGCC  PMT,  front
         8 - 11  FAEC front,middle,inner,rear 
         12 -15  LAEC front,middle,inner,rear 
         16,17,26 MRPC front,gas,glass
         27-28 FASPD front, inner
         41-42 LASPD front, inner
   pid =0   photon+electron+positron
        1   photon    
        2   electron + positron
        3   neutron
        4   proton
        5   pip
        6   pim
        7   Kp
        8   Km
        9   Kl
       10   other

run 1

It was done by Zhiwen Zhao in later 2013

main report for PVDIS_LD2 pptxpdf and SIDIS_He3 pptx pdf

The code are in SVN at https://jlabsvn.jlab.org/svnroot/solid/solid_gemc/analysistool/background see log files there for details

EM background is estimated from shooting beam into target with SoLID GEMC

background from eDIS, eES and hadron are using event generator at vertex as input into SoLID GEMC

low energy neutron cross section was turned on

(ask Zhiwen Zhao if you need to use them) 
file location
on disk  /cache/halla/solid/sim/solid_gemc/
on tape  /mss/halla/solid/sim/solid_gemc/
for SIDIS_He3, it's at subdir SIDIS_He3_run1/
for PVDIS_LD2, it's at subdir PVDIS_LD2_run1/, main result is dir baffle_babarbafflemore1_block which has photon block before EC, while baffle_babarbafflemore1 has result without the photon block
file names for SIDIS_He3
file with no special name is just shooting beam on target, "EM" particles dominating the results 
"other" in file name means kryptonite for all geometry 
"actual" in file name mean all real materials
"sum"  in file name mean sum over all real materials for gas and two windows
file names for PVDIS_LD2
file with no special name is just shooting beam on target, "EM" particles dominating the results, it has lead baffle and kryptonite for everything else
"other" in file name means kryptonite for all geometry 
"real" in file name mean lead baffle and kryptonite for everything else
for different detector and for different particles, we have histograms showing rate
"Eklog_R_hitid_pid" rate(kHz/mm2) at R(cm) and log10(Ek)(GeV) with bin(300, 0, 300, 200,-6,1.3)
"Eklog_R_high_hitid_pid" rate(kHz/mm2) at R(cm) and log10(Ek)(GeV) with bin(300, 0, 300, 200,-6,1.3) for Phi (0,6)deg
"Eklog_R_low_hitid_pid" rate(kHz/mm2) at R(cm) and log10(Ek)(GeV) with bin(300, 0, 300, 200,-6,1.3) for Phi (6,12)deg
"Eklog_R_Phi_hitid_pid" rate(kHz) at Phi(deg), R(cm) and log10(Ek)(GeV) with bin(48,0,12,300, 0, 300, 200,-6,1.3)
"P_R_hitid_pid" rate(kHz/mm2) at R(cm) and P(GeV) with bin(300, 0, 300, 1100,0,11)
there are many other histograms produced with similar names
Because all histograms are produced the same way. One can simply add root files together by "hadd" to look at the result in sum
 hitid =0 - 5  GEM plane 1 - 6
        6       LGCC  PMT
        18      LGCC  front
        7       HGCC  PMT         
        19      HGCC  front
        8 - 11  FAEC front,middle,inner,rear 
        12 -15  LAEC front,middle,inner,rear 
        16,17,26 MRPC front,glass,gas
        20 - 25  GEM plane 1 - 6 front
        27 - 28  SPD front,scintillator
  pid =0   photon+electron+positron
       1   photon    
       2   electron + positron
       3   neutron
       4   proton
       5   pip
       6   pim
       7   Kp
       8   Km
       9   Kl
      10   other