Benchmark GEMC Occupancy Prediction for 11GeV Electron Beam with 0T Solenoid for LH2

From New IAC Wiki
Revision as of 14:05, 1 August 2016 by Vanwdani (talk | contribs)
Jump to navigation Jump to search

Splitting an Isotropic distribution for an 11GeV Electron bean at 0T for the Solenoid field with an LH2 target, we start with 5 files LUND_Spread_00 - LUND_Spread_04. LUND_Spread_00 - LUND_Spread_03 are equally divided with 75000 lines or 25000 Moller events in each file. LUND_Spread_04 has 37731 lines or 12577 Moller events. These LUND files are run in GEMC with:

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,LUND_Spread_LH2_00.LUND" -N=75000 eg12.gcard > tracking.dat

Using the evio file created from the GEMC run:

~/src/CLAS/evioRoot/bin/evio2root eg12.evio eg12.root 25000


LUND_Spread_01 has too many events with respect to the 75000 events GEMC can handle. Initially, LUND_Spread_01 had 25000 Moller events or 75000 lines(Header, Scattered, and Moller) in the file. Do to the scattering and secondary particles, this file is divided into two parts LUND_Spread_01_part_00 and LUND_Spread_01_part_01 of 37500 lines or 12500 Moller events. The second file is also two large and is similarly split into two parts LUND_Spread_01_part_01_00 and LUND_Spread_01_part_01_01 of 18750 lines or 6250 Moller events.


The root files can be combined with:

hadd 0T_LH2.root ./LUND_Spread_00/eg12.root LUND_Spread_01/eg12_00.evio LUND_Spread_01/eg12_01_00.evio LUND_Spread_01/eg12_01_01.evio LUND_Spread_02/eg12.root LUND_Spread_03/eg12.root LUND_Spread_04/eg12.root