Revision as of 00:29, 8 November 2013

some estimation

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.7e36/cm2/s

SIDIS proton

100nA
NH3

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

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 textbook formula, all DIS hadrons crossection is from Wiser fit. particles from resonance region is still missing as 2013/03

background

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

background from eDIS, eES and DIS hadron are using the rate at vertex as input into SoLID GEMC

The code is in SVN at https://jlabsvn.jlab.org/svnroot/solid/solid_gemc/analysistool/background

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

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, 50,-6,1.3)
"Eklog_R_high_hitid_pid" rate(kHz/mm2) at R(cm) and log10(Ek)(GeV) with bin(300, 0, 300, 50,-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, 50,-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, 50,-6,1.3)
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 MRPC front,within
        20 - 25  GEM plane 1 - 6 front
  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:

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.

SIDIS Geant4 modeled background from Target

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

See here