Difference between revisions of "Solid Background"

From Hall A Wiki
Jump to: navigation, search
(Luminosity and radiation thickness)
(Link to new wiki)
 
(18 intermediate revisions by one other user not shown)
Line 1: Line 1:
= Luminosity and radiation thickness =
+
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 He3
 
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, 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,
+
= rate =
  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
+
  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 (from hallB note and thesis), 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
+
Luminosity of NH3 = den * t * PF * (I/1.602e-19)*6.022e23 = 0.917*2.826*0.55*(100e-9/1.602e-19)*6.022e23 = 5.36e35/cm2/s
+
pol Lumi = 5.36e35/17*3 = 0.946e35/cm2/s
+
Luminosity of He4 = 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 He4 outside  =  0.145*2*0.432*(100e-9/1.602e-19)*6.022e23 = 0.47e35/cm2/s
+
Luminosity of Al foil window = 2.7*2*0.001778*(100e-9/1.602e-19)*6.022e23 = 0.036e35/cm2/s
+
Total luminosity = (5.36+0.69+0.47+0.036) * 1e35 = 6.56e35/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 48: 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 58: 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 73: 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 =

Latest revision as of 15:39, 7 December 2020

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


rate

The result is at /work/halla/solid/evgen

SIDIS e and hadron rate

generated by "collider", the code and script for analyzing result are at https://jlabsvn.jlab.org/svnroot/solid/evgen/collider

The result is here

single particle rate

various particle from vertex is estimated by "eicRate" in our event generator repository. The code and script for analyzing result are at https://jlabsvn.jlab.org/svnroot/solid/evgen/eicRate_20101102 https://jlabsvn.jlab.org/svnroot/solid/evgen/eicRate_20101102/output

rate study

pdf Comparison with 6GeV transversity data by Xin Qian

pptx pdf event generator rate comparison, study by Zhiwen Zhao and Zhihong Ye

pdf study by Nguyen Ton and Xiaochao Zheng (outdated)


other staff

PVDIS Geant4 modeled background from Target

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:

Geant4 truegemrate.png

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.

Emflux.png

SIDIS Geant4 modeled background from Target

Sidisgem.png

<math>\pi</math>/e Ratio

See here