fee_gen_to_recon_job.py

Go to the documentation of this file.

"""!
Simulation of beam, beam signals in detector (using SLIC), and readout.
The simulation is followed by reconstruction of the events.
"""
from hpsmc.generators import EGS5
from hpsmc.tools import BeamCoords, RandomSample, SLIC, ExtractEventsWithHitAtHodoEcal, JobManager
 
job.description = 'beam from generation to slic'
 
 
inputs = list(job.input_files.values())
 
if 'nevents' in job.params:
    nevents = job.params['nevents']
else:
    nevents = 250000
 
 
egs5 = EGS5(name="beam_v7_%s" % job.params['run_params'])
 
 
rot = BeamCoords()
 
 
sample = RandomSample()
 
 
slic = SLIC(nevents=nevents + 1)
 
 
space_events = ExtractEventsWithHitAtHodoEcal(event_interval=250, num_hodo_hits=0)
 
 
readout = JobManager(steering='readout')
 
 
recon = JobManager(steering='recon')
 
 
job.add([egs5, rot, sample, slic, space_events, readout, recon])