Difference between revisions of "Solid Background"

From Hall A Wiki
Jump to: navigation, search
(→‎Luminosity and radiation thickness)
(→‎Luminosity and radiation thickness)
Line 27: Line 27:
  
 
SIDIS proton
 
SIDIS proton
−
  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.
+
  I = 100nA, Target = NH3, density of NH3(den) = 0.917 g/cm^3 (from hallB note and thesis like https://books.google.com/books?id=Jn-SXhcvHekC&pg=PA113&lpg=PA113&dq=0.917+ammonia+nh3+density&source=bl&ots=PPi0pdCS6B&sig=QHt7LlIWmS5PTPnxNG4wT-cfJik&hl=en&sa=X&ved=0ahUKEwjH0KWN4aTSAhUJOCYKHfWrCugQ6AEIGjAA#v=onepage&q=0.917%20ammonia%20nh3%20density&f=false)
 +
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
 
  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
 
  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

Revision as of 12:40, 23 February 2017

Luminosity and radiation thickness

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, 
  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, Target = NH3, density of NH3(den) = 0.917 g/cm^3 (from hallB note and thesis like https://books.google.com/books?id=Jn-SXhcvHekC&pg=PA113&lpg=PA113&dq=0.917+ammonia+nh3+density&source=bl&ots=PPi0pdCS6B&sig=QHt7LlIWmS5PTPnxNG4wT-cfJik&hl=en&sa=X&ved=0ahUKEwjH0KWN4aTSAhUJOCYKHfWrCugQ6AEIGjAA#v=onepage&q=0.917%20ammonia%20nh3%20density&f=false)
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 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 = (5.36+0.69+0.47) * 1e35 = 6.52e35/cm2/s
Note: there are several thin Al windows, but their luminosities are small

total length of NH3+LHe4 and two LHe4 outside, 2.826+2*0.432=3.69
average density of NH3+LHe4, (0.917*0.55+0.145*0.45) = 0.5696
average density of NH3+LHe4 and two LHe4 outside, (0.5696*2.826+0.145*2*0.432)/(2.826+2*0.432)=0.47
Composition by weight of NH3+LHe4 and two LHe4 outside
   NH3                   0.917*0.55*2.826/(2.826+2*0.432)/0.47=0.82
   LHe4                  0.145*0.45*2.826/(2.826+2*0.432)/0.47=0.11
   two LHe4 outside      0.145*2*0.432/(2.826+2*0.432)/0.470=0.07
average Z=10*0.82+2*0.11+2*0.07=8.56
average A=17*0.82+4*0.11+4*0.07=14.66
Z/A=0.584
assume Z=73 and A=125,which gives correct Z/A, then we got nuclei luminosity 6.52e35/125=0.0522e35/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

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

The result is here http://hallaweb.jlab.org/12GeV/SoLID/download/sim/rate/SIDIS_He3/double/

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. 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

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

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)

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 pptxpdf and SIDIS_He3 pptx 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 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

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

Retrieved from "https://hallaweb.jlab.org/wiki/index.php?title=Solid_Background&oldid=32358"