Commit 226861dc authored by mcjo's avatar mcjo
Browse files

Add 406000-5


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@666348 4525493e-7705-40b1-a816-d608a930855b
parent 301fa76f
2015-05-11 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* Add 406000-5
2015-05-11 James Monk <jmonk@cern.ch>
* Add DSID 429020 for Epos min bias
......@@ -16,7 +20,7 @@
* Add 341173-341176
* Add 423300-423303
* Add common/Filters/JetFilter.py
2015-05-07 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* Add 361250-3, 300998-9, BSignalFilter.py
......
# This comes after all Simplified Model setup files
from MadGraphControl.MadGraphUtils import SUSY_SM_Generation
# Set maximum number of events if the event multiplier has been modified
if evt_multiplier>0:
if runArgs.maxEvents>0:
nevts=runArgs.maxEvents*evt_multiplier
else:
nevts=5000*evt_multiplier
if njets<0:
evgenLog.fatal('njets is not set')
raise RuntimeError('njet is not set')
# Set up for grid pack running
gridpackDirName=None
if hasattr(runArgs, "inputGenConfFile"):
gridpackDirName='madevent/'
if 'writeGridpack' not in dir():
writeGridpack=False
if not SLHAonly:
[qcut,outputDS] = SUSY_SM_Generation(runArgs=runArgs,process=process,gentype=gentype,decaytype=decaytype,masses=masses,nevts=nevts,syst_mod=syst_mod, writeGridpack=writeGridpack, gridpackDirName=gridpackDirName, keepOutput=keepOutput, pdlabel="'nn23lo1'")
else:
[qcut,outputDS] = SUSY_SM_Generation(runArgs=runArgs,process=process,gentype=gentype,decaytype=decaytype,masses=masses,nevts=10000,syst_mod=None, SLHAonly=True, keepOutput=keepOutput)
if qcut<0 or outputDS is None or ''==outputDS:
evgenLog.warning('Looks like something went wrong with the MadGraph generation - bailing out!')
raise RuntimeError('Error in MadGraph generation')
import os
if 'ATHENA_PROC_NUMBER' in os.environ:
evgenLog.info('Noticed that you have run with an athena MP-like whole-node setup. Will re-configure now to make sure that the remainder of the job runs serially.')
njobs = os.environ.pop('ATHENA_PROC_NUMBER')
# Try to modify the opts underfoot
if not hasattr(opts,'nprocs'): mglog.warning('Did not see option!')
else: opts.nprocs = 0
print opts
runArgs.qcut = qcut
runArgs.inputGeneratorFile = outputDS
if 'syst_mod' in dir():
runArgs.syst_mod = syst_mod
runArgs.decaytype = decaytype
runArgs.gentype = gentype
# Pythia8 setup
genSeq.Pythia8.Commands += ["Init:showAllParticleData = on",
"Next:numberShowLHA = 10",
"Next:numberShowEvent = 10",
]
if njets>0:
genSeq.Pythia8.Commands += ["Merging:mayRemoveDecayProducts = on",
"Merging:nJetMax = "+str(njets),
"Merging:doKTMerging = on",
"Merging:TMS = "+str(qcut),
"Merging:ktType = 1",
"Merging:Dparameter = 0.4",
"Merging:nQuarksMerge = 4"]
if hasattr(runArgs,'syst_mod') and runArgs.syst_mod is not None:
evgenLog.error('No known equivalent for systematics in Pythia8!')
# Configuration for EvgenJobTransforms
#--------------------------------------------------------------
evgenConfig.generators += ["EvtGen"]
if not hasattr(runArgs,'inputGeneratorFile'):
print 'ERROR: something wasnt write in file name transfer from the fragment.'
runArgs.inputGeneratorFile='madgraph.*._events.tar.gz'
evgenConfig.description = 'SUSY Simplified Model generation with MadGraph/Pythia8 in MC15'
evgenConfig.keywords += ["SUSY"]
evgenConfig.inputfilecheck = runArgs.inputGeneratorFile.split('._0')[0]
# Generator transform pre-include
# Gets us ready for on-the-fly SUSY SM generation
# Simple variable setups
njets = -1
masses = {}
# Event multipliers for getting more events out of madgraph to feed through athena (esp. for filters)
evt_multiplier = 2.0
# Useful definitions
squarks = []
squarksl = []
for anum in [1,2,3,4]:
squarks += [str(1000000+anum),str(-1000000-anum),str(2000000+anum),str(-2000000-anum)]
squarksl += [str(1000000+anum),str(-1000000-anum)]
dict_index_syst = {0:'scalefactup',
1:'scalefactdown',
2:'alpsfactup',
3:'alpsfactdown',
4:'moreFSR',
5:'lessFSR',
6:'qup',
7:'qdown'}
# Basic settings for production and filters
SLHAonly = False
syst_mod=None
if 'EventMultiplier' in dir(): evt_multiplier=EventMultiplier
if hasattr(runArgs,'EventMultiplier'): evt_multiplier=runArgs.EventMultiplier
nevts=5000*evt_multiplier
evt_multiplier=-1
# in case someone needs to be able to keep the output directory for testing
keepOutput = False
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000005'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
if masses['1000022']<0.5: masses['1000022']=0.5
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > b1 b1~ $ go susylq susylq~ b2 t1 t2 b2~ t1~ t2~ @1
add process p p > b1 b1~ j $ go susylq susylq~ b2 t1 t2 b2~ t1~ t2~ @2
add process p p > b1 b1~ j j $ go susylq susylq~ b2 t1 t2 b2~ t1~ t2~ @2
'''
njets = 2
evgenLog.info('Registered generation of sbottom pair production, sbottom to b+LSP; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000005']))
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'sbottom']
evgenConfig.description = 'sbottom direct pair production, sb->b+LSP in simplified model'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
genSeq.Pythia8.Commands += ["Merging:Process = pp>{b1,1000005}{b1~,-1000005}"]
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000024'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000023'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
mslep = (masses['1000024']+masses['1000022'])*0.5
masses['1000011'] = mslep
masses['1000012'] = mslep
masses['1000013'] = mslep
masses['1000014'] = mslep
masses['1000015'] = mslep
masses['1000016'] = mslep
if masses['1000022']<0.5: masses['1000022']=0.5
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > x1+ n2 $ susystrong @1
add process p p > x1- n2 $ susystrong @1
add process p p > x1+ n2 j $ susystrong @2
add process p p > x1- n2 j $ susystrong @2
add process p p > x1+ n2 j j $ susystrong @3
add process p p > x1- n2 j j $ susystrong @3
'''
njets = 2
evgenLog.info('Registered generation of ~chi1+/- ~chi20 production, decay via slepton; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000024']) + ' ' + str(masses['1000022']))
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['gaugino', 'chargino', 'neutralino', 'slepton']
evgenConfig.description = '~chi1+/- ~chi20 production, decay via slepton in simplified model'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
genSeq.Pythia8.Commands += [ "Merging:Process = pp>{x1+,1000024}{x1-,-1000024}{n2,1000023}",
"1000024:spinType = 1",
"1000023:spinType = 1" ]
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000024'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000023'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
if masses['1000022']<0.5: masses['1000022']=0.5
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > x1+ n2 $ susystrong @1
add process p p > x1- n2 $ susystrong @1
add process p p > x1+ n2 j $ susystrong @2
add process p p > x1- n2 j $ susystrong @2
add process p p > x1+ n2 j j $ susystrong @3
add process p p > x1- n2 j j $ susystrong @3
'''
njets = 2
evgenLog.info('Registered generation of ~chi1+/- ~chi20 production, decay via WZ; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000024']) + ' ' + str(masses['1000022']))
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['gaugino', 'chargino', 'neutralino']
evgenConfig.description = '~chi1+/- ~chi20 production, decay via WZ in simplified model'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if len(runArgs.jobConfig[0].split('_'))>6 and '3lep' in runArgs.jobConfig[0].split('_')[6].split('.'):
evgenLog.info('3lepton processes will be only generated')
genSeq.Pythia8.Commands += [
"23:onMode = off",
"23:onIfAny = 11 13 15", # Z->ll
"24:onMode = off",
"24:onIfAny = 11 12 13 14 15 16" # W->lnu
]
else:
evgenLog.info('inclusive processes will be generated')
if njets>0:
genSeq.Pythia8.Commands += [ "Merging:Process = pp>{x1+,1000024}{x1-,-1000024}{n2,1000023}",
"1000024:spinType = 1",
"1000023:spinType = 1" ]
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000021'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
if masses['1000022']<0.5: masses['1000022']=0.5
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > go go $ susysq susysq~ @1
add process p p > go go j $ susysq susysq~ @2
add process p p > go go j j $ susysq susysq~ @3
'''
njets = 2
evgenLog.info('Registered generation of gluino grid '+str(runArgs.runNumber))
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel','gluino']
evgenConfig.description = 'gluino production in simplified model'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
genSeq.Pythia8.Commands += ["Merging:Process = pp>{go,1000021}{go,1000021}"]
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
for q in squarks: masses[q] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
if masses['1000022']<0.5: masses['1000022']=0.5
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > susylq susylq~ $ go susyweak @1
add process p p > susylq susylq~ j $ go susyweak @2
add process p p > susylq susylq~ j j $ go susyweak @3
'''
njets = 2
evgenLog.info('Registered generation of squark grid '+str(runArgs.runNumber))
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'squark']
evgenConfig.description = 'squark production in simplified model'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
genSeq.Pythia8.Commands += ["Merging:Process = pp>{ul,1000002}{ul~,-1000002}{ur,2000002}{ur~,-2000002}{dl,1000001}{dl~,-1000001}{dr,2000001}{dr~,-2000001}{sl,1000003}{sl~,-1000003}{sr,2000003}{sr~,-2000003}{cl,1000004}{cl~,-1000004}{cr,2000004}{cr~,-2000004}"]
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000006'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
if masses['1000022']<0.5: masses['1000022']=0.5
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > t1 t1~ $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
add process p p > t1 t1~ j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @2
add process p p > t1 t1~ j j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @2
'''
njets = 2
evgenLog.info('Registered generation of stop pair production, stop to t+LSP; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000006']))
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'stop']
evgenConfig.description = 'stop direct pair production, st->t+LSP in simplified model'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
genSeq.Pythia8.Commands += ["Merging:Process = pp>{t1,1000006}{t1~,-1000006}"]
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_SS_direct.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_TT_directTT.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_BB_direct.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_C1N2_WZ.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_C1N2_Slep.py' )
......@@ -816,7 +816,6 @@ MC15.301939.Pythia8EvtGen_A14NNPDF23LO_Bstar3000.py
MC15.301940.Pythia8EvtGen_A14NNPDF23LO_Bstar4000.py
MC15.301941.Pythia8EvtGen_A14NNPDF23LO_Bstar5000.py
DSID302xxx:
DSID303xxx:
......@@ -1425,6 +1424,12 @@ DSID403xxx:
DSID404xxx:
DSID406xxx:
MC15.406000.MadGraphPythia8EvtGen_A14NNPDF23LO_GG_direct_1500_0.py
MC15.406001.MadGraphPythia8EvtGen_A14NNPDF23LO_SS_direct_1050_0.py
MC15.406002.MadGraphPythia8EvtGen_A14NNPDF23LO_TT_directTT_800_1.py
MC15.406003.MadGraphPythia8EvtGen_A14NNPDF23LO_BB_direct_800_1.py
MC15.406004.MadGraphPythia8EvtGen_A14NNPDF23LO_C1N2_WZ_400_100_3lep.py
MC15.406005.MadGraphPythia8EvtGen_A14NNPDF23LO_C1N2_Slep_600_300.py
DSID407xxx:
MC15.407013.Sherpa_CT10_2DP50_DirectPhotonFilter.py
......@@ -2085,3 +2090,4 @@ MC15.429008.Pythia8_AU2MSTW2008LO_zprime1000_tt.py
MC15.429009.Pythia8EvtGen_AU2MSTW2008LO_zprime1000_tt.py
MC15.429010.Pythia8_A14NNPDF23LO_jetjet_JZ4.py
MC15.429011.PowhegPythia8_AZNLOCTEQ6L1_Ztautau.py
MC15.429020.Epos_minbias_inelastic.py
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment