rad_lhe_to_slic_job.py

Go to the documentation of this file.

"""!
@file rad_lhe_to_slic_job.py
 
Radiative (rad) events from a pre-generated LHE file through SLIC.
 
"""
from hpsmc.generators import StdHepConverter
from hpsmc.tools import Unzip, AddMotherFullTruth, BeamCoords, SLIC
 
job.description = 'rad from LHE to slic'
 
 
inputs = list(job.input_files.values())
 
nevents = job.params.get('nevents', 10000)
 
 
converter = job.params.get('lhe_converter', 'lhe_rad')
 
 
unzip = Unzip(inputs=inputs, outputs=['rad.lhe'])
 
 
cnv = StdHepConverter(name=converter, inputs=inputs)
 
 
mom = AddMotherFullTruth(inputs=[cnv.output_files()[0], unzip.output_files()[0]],
                         outputs=['rad_mom.stdhep'])
 
 
rot = BeamCoords(inputs=mom.output_files(), outputs=['rad_mom_rot.stdhep'])
 
 
slic = SLIC(inputs=rot.output_files(), outputs=['rad_mom_rot.slcio'],
            nevents=nevents + 1)
 
 
job.add([unzip, cnv, mom, rot, slic])