SoLID full simulation and file sharing

From Hall A Wiki
Revision as of 18:04, 8 November 2017 by Zwzhao (Talk | contribs) (New page: = 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, GE...)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

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

Luminosity and radiation thickness

file sharing

The file size is large in order of 50GB, you can transfer them by using "scp jlabname@ftp.jlab.org:file" slowly or other faster ways like globus here https://cc.jlab.org/filetransfers

  • sharing location
location on disk /cache/halla/solid/sim/solid_gemc/
location on tape /mss/halla/solid/sim/solid_gemc/
dir
  PVDIS_LD2_JLAB_VERSION_1.2
  PVDIS_LD2_JLAB_VERSION_1.3
  SIDIS_He3_JLAB_VERSION_1.2
  SIDIS_NH3_JLAB_VERSION_1.3
  JPsi_LH2_JLAB_VERSION_1.2   
subdir
  pass1
  pass2
  pass3
  pass4
  ...
refer details
  /home/zwzhao/solid/solid_svn/study/background/log_PVDIS_LD2
  /home/zwzhao/solid/solid_svn/study/background/log_SIDIS_He3
  /home/zwzhao/solid/solid_svn/study/background/log_SIDIS_NH3
  /home/zwzhao/solid/solid_svn/study/background/log_JPsi_LH2
  • 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

old

run 2

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

working dir with all output /work/halla/solid/sim/solid_gemc/*_geant4.9.6.p02 and /work/halla/solid/sim/solid_gemc/*_geant4.9.5.p01

table for normalization factor comparison pdf pptx

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

It was done by Zhiwen Zhao in later 2013

main report for PVDIS_LD2 pptxpdf and SIDIS_He3 pptx 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

file location for SIDIS_He3
/work/halla/solid/sim/solid_gemc/SIDIS_He3_run1
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 location for PVDIS_LD2
/work/halla/solid/sim/solid_gemc/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 PVDIS_LD2
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 
"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