|
|
(29 intermediate revisions by 2 users not shown) |
Line 1: |
Line 1: |
− | = some estimation =
| + | This page has been moved to [https://solid.jlab.org/wiki/index.php/Background https://solid.jlab.org/wiki/index.php/Background] in the new [https://solid.jlab.org/wiki SoLID Wiki]. Please do not edit this page. |
| | | |
− | PVDIS
| |
− | 50uA,
| |
− | 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 3he
| |
− | 15uA
| |
− | 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
| |
− | 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,
| + | = rate = |
− | Molecule Composition by weight
| + | |
− | SiO 260.3%
| + | |
− | BaO 18.2%
| + | |
− | Al2O3 14.3%
| + | |
− | CaO 6.5%
| + | |
− | SrO 0.25%
| + | |
− | Z/A 0.4829
| + | |
− | radiation thickness 19.4246 g/cm2
| + | |
− | denisty 2.76 g/cm3
| + | |
− | radition thickness 7.038 cm
| + | |
− | refer to http://galileo.phys.virginia.edu/research/groups/spinphysics/glass_properties.html
| + | |
| | | |
− | SIDIS proton
| + | The result is at /work/halla/solid/evgen |
− | I = 100nA, Target = NH3, density of NH3(den) = 0.917 g/cm^3, target thickness (t) = 3 cm, packing fraction(PF) = 0.55.
| + | |
− | | + | |
− | NH3 target is immersed in liquid He4, density of LHe4 = 0.145 g/cm^3
| + | |
− | | + | |
− | Luminosity of NH3 = den * t * PF * (I/1.602*10^-19) * 6.022*10^23 = 0.917 * 3 * 0.55 * (100*10^-9/1.602*10^-19) * 6.022*10^23 = 5.68*10^35
| + | |
− | | + | |
− | Luminosity of He4 = den * t * PF * (I/1.602*10^-19) * 6.022*10^23 = 0.145 * 3 * 0.45 * (100*10^-9/1.602*10^-19) * 6.022*10^23 = 0.73*10^35
| + | |
− | | + | |
− | Total luminosity = 5.68+0.73 * 10^35 = 6.41 * 10^35/cm2/s
| + | |
− |
| + | |
− | | + | |
− | | + | |
− | JPsi
| + | |
− | 3uA
| + | |
− | 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
| + | |
− | | + | |
− | = rate =
| + | |
| | | |
| == SIDIS e and hadron rate == | | == SIDIS e and hadron rate == |
Line 51: |
Line 13: |
| | | |
| The result is here | | The result is here |
− | http://hallaweb.jlab.org/12GeV/SoLID/download/sim/rate/SIDIS_He3/double/
| |
| | | |
| == single particle rate == | | == single particle rate == |
− |
| |
− | individual particle rate and pi/e ratio at vertex are available here
| |
| | | |
| various particle from vertex is estimated by "eicRate" in our event generator repository. | | various particle from vertex is estimated by "eicRate" in our event generator repository. |
Line 61: |
Line 20: |
| https://jlabsvn.jlab.org/svnroot/solid/evgen/eicRate_20101102 | | https://jlabsvn.jlab.org/svnroot/solid/evgen/eicRate_20101102 |
| https://jlabsvn.jlab.org/svnroot/solid/evgen/eicRate_20101102/output | | https://jlabsvn.jlab.org/svnroot/solid/evgen/eicRate_20101102/output |
− |
| |
− | The result is here
| |
− | http://hallaweb.jlab.org/12GeV/SoLID/download/sim/rate
| |
− | *root is eicRate output tree, *.lund can be used as input for SoLID GEMC
| |
− |
| |
− | note: in eicRate, eDIS and eES crosssection is from CTEQ6 fit, all DIS hadrons crossection is from Wiser fit particles from resonance region is still missing as 2013/03
| |
| | | |
| == rate study == | | == rate study == |
Line 76: |
Line 29: |
| [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/wiser_table.pdf pdf] study by Nguyen Ton and Xiaochao Zheng (outdated) | | [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/talk/wiser_table.pdf pdf] study by Nguyen Ton and Xiaochao Zheng (outdated) |
| | | |
− | = background =
| |
− |
| |
− | == intro ==
| |
− | EM background is estimated from shooting beam into target with SoLID GEMC
| |
− |
| |
− | background from eDIS, eES and hadron are using the rate from code "eicRate" at vertex as input into SoLID GEMC
| |
− |
| |
− | low energy neutron crsssection was turned
| |
− |
| |
− | == run 1 ==
| |
− | It was done by Zhiwen Zhao in later 2013
| |
− |
| |
− | The code is in SVN at
| |
− | https://jlabsvn.jlab.org/svnroot/solid/solid_gemc/analysistool/background
| |
− |
| |
− | read the log files "log_PVDIS_LD2" and "log_SIDIS_He3" in SVN
| |
− |
| |
− | The output files are here
| |
− | http://hallaweb.jlab.org/12GeV/SoLID/download/sim/background
| |
− |
| |
− | working dir with all files on ifarm
| |
− | /home/zwzhao/work_halla/solid/solid_svn/solid/solid_gemc/analysistool/background
| |
− |
| |
− | 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]
| |
− |
| |
− | 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
| |
− |
| |
− | == run 2 ==
| |
− | It was done by Zhiwen Zhao in later 2014
| |
− |
| |
− | The code is in SVN at
| |
− | https://jlabsvn.jlab.org/svnroot/solid/study/background
| |
− |
| |
− | read the log files "log_PVDIS_LD2" and "log_SIDIS_He3" in SVN
| |
− |
| |
− | working dir with all files on ifarm
| |
− | /home/zwzhao/work_halla/solid/solid_svn/solid/study/background
| |
− |
| |
− | table for normalization factor comparison
| |
− | [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/background/normalization_factor_comparison.pdf pdf] [http://hallaweb.jlab.org/12GeV/SoLID/download/sim/background/normalization_factor_comparison.pdf pptx]
| |
| | | |
− | more hit_id are added, some are changed
| |
− | ========= hit_id and pid defination ==============
| |
− | 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
| |
| | | |
| = other staff = | | = other staff = |
Rates and flux from GEMC for PVDIS follow. Validation against GEANT3 can be found here: Compare_to_geant3_result.
Rates on GEMs for 50uA on LH2. Solid lines are unbaffled, dashed are with lead baffles:
Flux in EM calorimeter for 50uA on LH2. Solid lines are unbaffled, dashed are with lead baffles. The true rate should be taken from the black lines.