seutils.debug()
cmssw = qondor.init_cmssw()
mz = int(qondor.get_var('mz'))
max_events = 10000 # runSVJ wants maxEvents to be passed, just set it large enough
physics = {
'year': 2018,
'mz': mz,
'mdark': 20.,
'rinv': 0.3,
'boost': 0.,
}
part = qondor.get_proc_id() + 1
for i_rootfile, rootfile in enumerate(qondor.get_item()):
# Ensure unique part name per rootfile across jobs; Not sure if necessary for SIM>MINIAOD
svjpart = 1000 * part + i_rootfile
qondor.logger.info('Job %s: Processing %s', qondor.get_proc_id(), rootfile)
extended_physics = dict(physics, part=svjpart, max_events=max_events)
expected_infile = osp.join(
cmssw.cmssw_src, 'SVJ/Production/test',
svjqondor.formatted_filename('step1_LHE-GEN', **extended_physics))
seutils.cp(rootfile, expected_infile)
svjqondor.run_step_cmd(cmssw,
inpre='step1_LHE-GEN',
outpre='step1_SIM',
#$ pip module-install svjqondor
#$ file cmssw_tarball root://cmseos.fnal.gov//store/user/klijnsma/semivis/cmssw_tarballs/CMSSW_10_2_18_htfilterprintout.tar.gz
#$ items 150 250 450 650
import qondor, svjqondor, os.path as osp
cmssw = qondor.init_cmssw()
physics = {
'year': 2018,
'mz': float(qondor.get_item()),
'mdark': 20.,
'rinv': 0.3,
'boost': 1000.,
}
svjqondor.download_mg_tarball(dst=osp.join(cmssw.cmssw_src,
'SVJ/Production/test'),
**physics)
# Aim for 25 events times a factor
max_events = {
250: 2000 * 2,
450: 150 * 10,
50: 15000, # aim for 10 instead of 25
650: 85 * 10
}[int(qondor.get_item())]
expected_outfile = svjqondor.run_step_cmd(cmssw,
inpre='step0_GRIDPACK',
outpre='step1_LHE-GEN',
import qondor, seutils
cmssw = qondor.init_cmssw()
cmssw.run_commands([
'cd TreeMaker/Production/test', 'cmsRun runMakeTreeFromMiniAOD_cfg.py'
' numevents=-1'
' outfile=outfile'
' scenario=Autumn18sig'
' lostlepton=0'
' doZinv=0'
' systematics=0'
' deepAK8=0'
' deepDoubleB=0'
' doPDFs=0'
' nestedVectors=False'
' debugjets=1'
' splitLevel=99' + (' dump=1' if args.dump else '') +
(' inputFiles=' + ','.join(qondor.get_item()))
])
expected_outfile = osp.join(cmssw.cmssw_src, 'TreeMaker/Production/test',
'outfile_RA2AnalysisTree.root')
seutils.cp(
expected_outfile,
'root://cmseos.fnal.gov//store/user/lpcsusyhad/SVJ2017/boosted/treemaker/ht1000_nonested_{date}_mz{mz}/{part}.root'
.format(date=qondor.get_submission_time().strftime('%b%d'),
mz=qondor.get_var('mz'),
part=qondor.get_proc_id()))
import os.path as osp, argparse, re
import qondor, seutils, svjqondor
cmssw = qondor.init_cmssw()
mz = int(qondor.get_var('mz'))
max_events = 10000 # runSVJ wants maxEvents to be passed, just set it large enough
physics = {
'year': 2018,
'mz': mz,
'mdark': 20.,
'rinv': 0.3,
'boost': 1000.,
}
rootfiles = qondor.get_item()
qondor.logger.info('Job %s: Processing %s', qondor.get_proc_id(), rootfiles)
for rootfile in rootfiles:
part = int(osp.basename(rootfile).replace(
'.root', '')) # Get part from the gen step
extended_physics = dict(physics, part=part, max_events=max_events)
expected_infile = osp.join(
cmssw.cmssw_src, 'SVJ/Production/test',
svjqondor.formatted_filename('step1_LHE-GEN', **extended_physics))
if not osp.isfile(expected_infile): seutils.cp(rootfile, expected_infile)
svjqondor.run_step_cmd(cmssw,
inpre='step1_LHE-GEN',
outpre='step1_SIM',
**extended_physics)