View source for Full simulation and file sharing

= 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 
 I=50uA, raster 0.5cm diameter, target center Z=10cm relative to SoLID solenoid center z=0cm
 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*120um, density 2.7g/cm3, X=2*120e-4/8.897=2.7e-3, Lumi=50e-6/1.6e-19*2*120e-4*2.7*6.02e23=1.2e37/cm2/s
 material after vertex and before GEM, target side wall Al 400um, scattering chamber window Al 200um, air from target to GEM 100cm, X=(400e-4+200e-4)/8.897+(100/3.039E+04)=1e-2

SIDIS He3
 [https://logbooks.jlab.org/entry/3852850 Jlab HallC 12GeV He3 target ladder]
 I=15uA, raster 0.5cm diameter, target center Z=-350cm relative to SoLID solenoid center z=0cm
 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

{| Border = "1"
|-
! He3 target at 15uA
! Be up
! N2 up
! glass up
! He3
! glass down
! N2 down
! Be down
|-
! z (cm)
! -375
! -372.5
! -370
! -350
! -330
! -327.5
! -325
|-
! length (cm)
! 0.01*2.54
! 5
! 0.012
! 40
! 0.012
! 5
! 0.02*2.54
|-
! Rad L (1e-3)
! 0.7
! 0.15
! 1.7
! 0.8
! 1.7
! 0.15
! 1.4
|-
! lumi (1e36/cm2/s)
! 2.65
! 0.64
! 1.87
! 3
! 1.87
! 0.64
! 5.3
|-
|}

(HallC A1n He3 2020 run used two Al windows at up/downstream 0.0002inch/0.001inch to cover Be windows. it created white powders on target cell we won't use them in future)

 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

 material after vertex and before GEM, target cell side wall GE180 0.1cm, air from target to GEM 200cm, X=0.1/7.038+(200/3.039E+04)=2e-2, (window collimator not included)

SIDIS NH3
 
 I = 100nA, raster 2.4cm diameter,target center Z=-350cm relative to SoLID solenoid center z=0cm

 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 [https://nvlpubs.nist.gov/nistpubs/Legacy/TN/nbstechnicalnote1334.pdf]
 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, ~0.1cm lumi=100e-9/1.6e-19*0.1*2.7*6.02e23=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

 X0_NH3=40.87g/cm2, X0_LHe=94.32g/cm2,target X0=(0.819*0.55+0.145*0.45)/(0.819*0.55/40.87+0.145*0.45/94.32)=44g/cm2 (refer to https://cds.cern.ch/record/1279627/files/PH-EP-Tech-Note-2010-013.pdf)
 target X=2.826/(44/(0.819*0.55+0.145*0.45))=2e-2 
 two LHe4 outside, X=2*0.432/(94.32/0.145)=1.3e-3
 material before target, Al windows 0.03cm, X=0.03/8.897=3.4e-3
 material after vertex and before GEM, Al windows 0.07cm, air from target to GEM 200cm, X=0.07/8.897+200/3.039E+04=1.44e-2

 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
 I=3uA, raster 0.5cm diameter,target center Z=-315cm relative to SoLID solenoid center z=0cm
 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
 LD target, 15cm, 0.169g/cm3, X=15/745.4=2.0e-2, Lumi=1.25e-6/1.6e-19*15*0.169*6.02e23=1.2e37/cm2/s
 Al window, 2*120um, density 2.7g/cm3, X=2*120e-4/8.897=2.7e-3, Lumi=3e-6/1.6e-19*2*120e-4*2.7*6.02e23=7.3e35/cm2/s
 material after vertex and before GEM, target side wall Al 178um, air from target to GEM 200cm, X=178e-4/8.897+200/3.039E+04=0.85e-2

= file sharing =

location

work dir /work/halla/solid/sim/solid_gemc/ 

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

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

temp output work dir /volatile/halla/solid/sim/solid_gemc/ 

== run 3  (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)

* PVDIS
<pre>
PVDIS_LD2_JLAB_VERSION_1.2/pass1,in 2016/01, new GEM, no LGC, no EC
PVDIS_LD2_JLAB_VERSION_1.2/pass2,in 2016/02, has LGC, has EC as whole, with mother particles
PVDIS_LD2_JLAB_VERSION_1.3/pass3,in 2016/03, add LGC virtual plane,add ec_photon_block,add neutron shielding
PVDIS_LD2_JLAB_VERSION_1.3/pass4,2017/11, svn 1193,in 2016/12, use CLEO2 baffle to replace babarmore1 baffle, use ec_segmented to replace ec as a whole, change scattering chamber front and back windows from 200um Al to vacuum, GEM size match 21-36 deg
PVDIS_LD2_JLAB_VERSION_1.3/pass5,2017/12, same like pass4, but on centos7.2 instead of centos 6.5
PVDIS_LD2_JLAB_VERSION_1.3/pass6,2020 spring, same like pass5, using solid_gemc_1.0.0 with container

PVDIS_JLAB_VERSION_2.5/pass3, 2022 fall  , solid_gemc commit1765a83_20221123, solenoid_v4,ec with rod,lgc optimized for small angle, ec small overlap fixed
</pre>

* SIDIS He3
<pre>
SIDIS_He3_JLAB_VERSION_1.2/pass1, new GEM, no LGC, no EC, in (2016/01)
SIDIS_He3_JLAB_VERSION_1.2/pass2, optimize hgc,has LGC, has EC as whole, (2016/02)
SIDIS_He3_JLAB_VERSION_1.3/pass3, turn on mother particle, add LGC virtual plane, use TungstenPowder for He3 collimator which has 60% density of Tungsten, use external materiel for H3 as SL_target_He3_He3_10amg instead of internal He3_10amg and the result should be same (2016/07)
SIDIS_He3_JLAB_VERSION_1.3/pass4, ,change SPD scintilator material from PVT alike to PS, 
SIDIS_He3_JLAB_VERSION_1.3/pass5, fix LASPD overlap with GEM, change LASPD from 5mm to 20mm, FASPD from 3mm to 5mm, use SPD with scintilator material PVT alike, use segmented EC with scintilator material PS
SIDIS_He3_JLAB_VERSION_1.3/pass6,svn rev1153,change LGC mirror,HGC,GEM sector numbering to make its 1st sector center at 96deg, LGC PMT 1st is center at 97.7deg, fix LGC H12700 QE near 200nm from 0.03 to 0.3,GEM has 7 planes and grouped
SIDIS_He3_JLAB_VERSION_1.3/pass7,svn rev1169,GEM change back to 6 planes without grouping, used for 2019 precdr
SIDIS_He3_JLAB_VERSION_1.3/pass8,2019 spring, same like pass7, study HGC window thickness
SIDIS_He3_JLAB_VERSION_1.3/pass9,2020 spring, same like pass7, using solid_gemc_1.0.0 with container
SIDIS_He3_JLAB_VERSION_1.3/pass10,2020 fall, same like pass7

SIDIS_He3_JLAB_VERSION_2.5/pass1, 2022 spring, solid_gemc commitf42e8a8_20220331, longer endcap with solenoid_v3,ec with rod,hgc Al window,SIDIS_He3 new collimator
SIDIS_He3_JLAB_VERSION_2.5/pass2, 2022 spring, solid_gemc commitf9df5d0_20220811, longer endcap with solenoid_v4, fix SIDIS_He3 new collimator overlap bug
SIDIS_He3_JLAB_VERSION_2.5/pass3, 2022 fall  , solid_gemc commit1765a83_20221123, lgc optimized for small angle, ec small overlap fixed

</pre>

* SIDIS NH3
<pre>
SIDIS_NH3_JLAB_VERSION_1.3/pass1,2017/11,svn rev1235, NH3 target updated, detector similar to SIDIS_He3_JLAB_VERSION_1.3/pass7, wrong target field map "solenoid_ptarget.dat" 
SIDIS_NH3_JLAB_VERSION_1.3/pass2,2017/11,svn rev1235, NH3 target updated, detector similar to SIDIS_He3_JLAB_VERSION_1.3/pass7, g2p target field map "g2p_ptarget.dat"
SIDIS_NH3_JLAB_VERSION_1.3/pass3,2017/11, svn rev1235, NH3 target updated, detector similar to SIDIS_He3_JLAB_VERSION_1.3/pass7, oxford target field map "oxford_ptarget.dat" 
SIDIS_NH3_JLAB_VERSION_1.3/pass4,2020 spring,github 3fbcee8, after gemc cylindrical-x transverse field bug fix 
SIDIS_NH3_JLAB_VERSION_1.3/pass5,2020 spring,github 43940e5, geometry same as pass4, test NH3 target blocks with Tungstengpowder near z=-320cm 
SIDIS_NH3_JLAB_VERSION_1.3/pass6,2020 spring,github 3fa5cd2, change NH3 scattering chamber size from 500cm diameter to 680cm diameter and window size from 25 to 28 deg 
SIDIS_NH3_JLAB_VERSION_1.3/pass7,2020 spring,github 3fa5cd2, geometry same as pass6, test NH3 target blocks with Tungstengpowder near z=-300cm
</pre>

*JPsi
<pre>
JPsi_LH2_JLAB_VERSION_1.2/pass1,in 2016/01, new GEM, no LGC, no EC,
JPsi_LH2_JLAB_VERSION_1.2/pass3,has LGC,in 2016/03, has EC as whole, with mother particles,add LGC virtual plane
JPsi_LH2_JLAB_VERSION_1.2/pass4,in 2016/05,add another virtual plane for muon large angle
JPsi_LH2_JLAB_VERSION_1.2/pass5,svn rev1171, in 2016/11, similar to SIDIS_He3_JLAB_VERSION_1.3/pass7
</pre>

* version used for preCDR study
   /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

* 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

== run 2 (old) ==
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


(ask Zhiwen Zhao if you need to use them)

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]

 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 (old) == 
It was done by Zhiwen Zhao in later 2013

main report for PVDIS_LD2 [http://hallaweb.jlab.org/12GeV/SoLID/download/baffle/talk/solid_baffle_zwzhao_20130910.pptx pptx][http://hallaweb.jlab.org/12GeV/SoLID/download/baffle/talk/solid_baffle_zwzhao_20130910.pdf pdf] and SIDIS_He3 [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/solid_SIDIS_background_zwzhao_20131118.pptx pptx] [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/solid_SIDIS_background_zwzhao_20131118.pdf 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