List of Hall C Replay Variables » History » Revision 4
Revision 3 (Richard Trotta, 07/01/2020 02:47 PM) → Revision 4/8 (Richard Trotta, 07/01/2020 02:49 PM)
h1. List of Hall C Replay Variables
[[what-you-need]]
= What you need:
ROOT 5.34/35+ gcc compiler (general)
[[source-codes]]
= Source codes:
Located in *src* directory...
*TDISMC_EIC.cpp* : pion structure function with ep scattering at JLEIC
*TDISMC_EICn.cpp* : pion structure function with eD scattering at JLEIC
*TDISMC_EICK.cpp* : kaon structure function with ep scattering at JLEIC
*cteq/* : cteqpdf.h and data based call files (c++ wrapper) *cteq-tbls/*
: nucleon PDFs table *structure_functions/* : various regularization
form for pion SF/FF
[[how-to-change-inputs]]
= How to change inputs:
Located in *inputs* directory...
*kinematics.inputs* : edit this document to change simulation kinematics
(e.g. number of events, x range, Q2 range, pbeam, kbeam)
All other constants are changed in *src/TDISMC_EIC.h*
↓ Below you can see the current kinematics inputs ↓
+*In[8]:*+
[source, ipython3]
----
!more inputs/kinematics.input
----
+*Out[8]:*+
----
XMIN=0.001
XMAX=1.00
Q2MIN=1.0
Q2MAX=1100.0
NEVTS=500000
PBEAM=135.0
KBEAM=10.0
----
[[how-to-run]]
= How to run:
*./run_batch.sh * : Final states...(pi/p, pi/n, k/lambda)
↓ Below you can see an example for a pion and neutron final state
simulation ↓
+*In[7]:*+
[source, ipython3]
----
!./run_mesonMC.sh pi/n
----
+*Out[7]:*+
----
Pion with neutron final state selected
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Your kinematics: [xBj_min:xBj_max] = [ 0.001000: 1.000000]
Your kinematics: [Q2_min:Q2_max] = [ 1.000000:1100.000000]
Incident Ion Mass 0.93827 GeV
Incident Electron, Ion Momenta: 10.0000, 135.00 GeV/c | s_0 = 5400.1019 GeV^2
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Warning in <TTree::Bronch>: TLorentzVector cannot be split, resetting splitlevel to 0
Total of 57310 events out of 500000 Trials ============================] 100 %
(int) 57310
----
[[root-and-lund-outputs]]
= ROOT and LUND outputs:
In the *OUTPUTS* directory are the ROOT and LUND outputs for the
simulation for further analysis.
[[running-geant4]]
= Running GEANT4
Located in *g4e_files/* directory...
*./run_g4e.sh* : Will run the python script for the GEANT4 simulation
from the TDIS_lund.dat file *eic_g4e.py* : This python script will run
GEANT4 simulation for the lund file specified for detector='jleic' and
beamline='erhic'
_This code is maintained by Richard Trotta (trotta@cua.edu)._