Project

General

Profile

Actions

List of Hall C Replay Variables » History » Revision 4

« Previous | Revision 4/8 (diff) | Next »
Richard Trotta, 07/01/2020 02:49 PM


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 ↓

In8:
[source, ipython3]
----
!more inputs/kinematics.input
----

Out8:
----
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 ↓

In7:
[source, ipython3]
----
!./run_mesonMC.sh pi/n
----

Out7:
----

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

Updated by Richard Trotta almost 4 years ago · 4 revisions