Commit 932dea87 authored by Franziska Iltzsche's avatar Franziska Iltzsche
Browse files

add 309810 - 309812


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@813581 4525493e-7705-40b1-a816-d608a930855b
parent 4b400ed0
2018-03-29 franziska iltzsche <franziska.speiser@cern.ch>
* add 309810 - 309812
* add JO and csv entry for 415300
* add topcontrol/MadGraph_param_card_ttFCNC_NLO_FixedFCNCBR.dat
* update common/Powheg/PowhegControl_ttFCNC_NLO.py
......
from MadGraphControl.MadGraphUtils import *
import re
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
import model VLQ_UFO_for_interference
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define bb = b b~
define WW = w+ w-
define tt = t t~
define ferm = ve vm vt ve~ vm~ vt~ mu- ta- e- mu+ ta+ e+
define VLQ = y y~
generate p p > j VLQ bb /h Z, (VLQ > WW bb)
output -f
""")
fcard.flush()
process_dir = new_process()
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
beamEnergy = 6500.
mass=int(re.findall(r'\d+',re.findall(r'\d+LH',runArgs.jobConfig[0])[0])[0])
print mass
# allow usage of all PDF sets
os.environ['LHAPATH']=os.environ["LHAPDF_DATA_PATH"]=(os.environ['LHAPATH'].split(':')[0])+":/cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current/"
# especially 263400 = NNPDF30_lo_as_0130_nf_4
extras = { 'lhe_version' : '2.0',
'cut_decays' : 'F',
'pdlabel' : "'lhapdf'",
'lhaid' : '263400',
'bwcutoff' : '10000'
}
safe_factor = 1.1
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir), run_card_new='run_card.dat', nevts=runArgs.maxEvents * safe_factor, rand_seed=runArgs.randomSeed, beamEnergy=beamEnergy, xqcut=0., extras=extras)
if not os.access(process_dir+'/Cards/param_card.dat',os.R_OK):
print 'ERROR: Could not get param card'
#elif os.access('param_card.dat',os.R_OK):
# print 'ERROR: Old run card in the current directory. Dont want to clobber it. Please move it first.'
else:
oldcard = open(process_dir+'/Cards/param_card.dat','r')
newcard = open('param_card.dat','w')
for line in oldcard:
if ' # MX ' in line:
newcard.write(' 6000005 %i # MX \n'%(mass))
elif ' # MTP ' in line:
newcard.write(' 6000006 %i # MTP \n'%(mass))
elif ' # MBP ' in line:
newcard.write(' 6000007 %i # MBP \n'%(mass))
elif ' # MY ' in line:
newcard.write(' 6000008 %i # MY \n'%(mass))
elif ' # KX ' in line:
newcard.write(' 1 0.4 # KX \n')
elif ' # KT ' in line:
newcard.write(' 2 0.4 # KT \n')
elif ' # KB ' in line:
newcard.write(' 3 0.4 # KB \n')
elif ' # KY ' in line:
newcard.write(' 4 0.4 # KY \n')
elif ' # WX ' in line:
newcard.write('DECAY 6000005 auto # WX \n')
elif ' # WTP ' in line:
newcard.write('DECAY 6000006 auto # WTP \n')
elif ' # WBP ' in line:
newcard.write('DECAY 6000007 auto # WBP \n')
elif ' # WY ' in line:
newcard.write('DECAY 6000008 auto # WY \n')
elif ' # zetaXuL ' in line:
newcard.write(' 1 0.000000e-01 # zetaXuL \n')
elif ' # zetaXcL ' in line:
newcard.write(' 2 0.000000e-01 # zetaXcL \n')
elif ' # zetaXtL ' in line:
newcard.write(' 3 1.000000e+00 # zetaXtL \n')
elif ' # zetaTuL ' in line:
newcard.write(' 4 0.000000e-01 # zetaTuL \n')
elif ' # zetaTcL ' in line:
newcard.write(' 5 0.000000e-01 # zetaTcL \n')
elif ' # zetaTtL ' in line:
newcard.write(' 6 1.000000e+00 # zetaTtL \n')
elif ' # zetaBdL ' in line:
newcard.write(' 7 0.000000e-01 # zetaBdL \n')
elif ' # zetaBsL ' in line:
newcard.write(' 8 0.000000e-01 # zetaBsL \n')
elif ' # zetaBbL ' in line:
newcard.write(' 9 1.000000e+00 # zetaBdL \n')
elif ' # zetaYdL ' in line:
newcard.write(' 10 0.000000e+00# zetaYdL \n')
elif ' # zetaYsL ' in line:
newcard.write(' 11 0.000000e+00 # zetaYsL \n')
elif ' # zetaYbL ' in line:
newcard.write(' 12 1.000000e+00 # zetaYbL \n')
elif ' # zetaXuR ' in line:
newcard.write(' 13 0.000000e+00 # zetaXuR \n')
elif ' # zetaXcR ' in line:
newcard.write(' 14 0.000000e+00 # zetaXcR \n')
elif ' # zetaXtR ' in line:
newcard.write(' 15 0.000000e+00 # zetaXtR \n')
elif ' # zetaTuR ' in line:
newcard.write(' 16 0.000000e+00 # zetaTuR \n')
elif ' # zetaTcR ' in line:
newcard.write(' 17 0.000000e+00 # zetaTcR \n')
elif ' # zetaTtR ' in line:
newcard.write(' 18 0.000000e+00 # zetaTtR \n')
elif ' # zetaBdR ' in line:
newcard.write(' 19 0.000000e+00 # zetaBdR \n')
elif ' # zetaBsR ' in line:
newcard.write(' 20 0.000000e+00 # zetaBsR \n')
elif ' # zetaBbR ' in line:
newcard.write(' 21 0.000000e+00 # zetaBdR \n')
elif ' # zetaYdR ' in line:
newcard.write(' 22 0.000000e+00 # zetaYdR \n')
elif ' # zetaYsR ' in line:
newcard.write(' 23 0.000000e+00 # zetaYsR \n')
elif ' # zetaYbR ' in line:
newcard.write(' 24 0.000000e+00 # zetaYbR \n')
# set some sm parameters to atlas defaults
elif ' # MB ' in line:
newcard.write(' 5 4.950000e+00 # MB \n')
elif ' # MT ' in line:
newcard.write(' 6 1.725000e+02 # MT \n')
elif ' # WT ' in line:
newcard.write('DECAY 6 1.350000e+00 # WT \n')
elif ' # WZ ' in line:
newcard.write('DECAY 23 2.495200e+00 # WZ \n')
elif ' # WW ' in line:
newcard.write('DECAY 24 2.085000e+00 # WW \n')
else:
newcard.write(line)
oldcard.close()
newcard.close()
runName='run_01'
generate(run_card_loc='./run_card.dat',param_card_loc='./param_card.dat',run_name=runName,proc_dir=process_dir)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
evgenConfig.description = "MadGraph+Pythia8 production JO with NNPDF30LN and A15NNPDF23LO for VLQ single Y to Wb with Y produced via W"
evgenConfig.keywords = ["BSM", "BSMtop", "exotic"]
evgenConfig.process = "Y_Wb"
evgenConfig.contact = ['fschenck@cern.ch']
evgenConfig.minevents = 5000
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
from MadGraphControl.MadGraphUtils import *
import re
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
import model VLQ_UFO_for_interference
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define bb = b b~
define WW = w+ w-
define tt = t t~
define ferm = ve vm vt ve~ vm~ vt~ mu- ta- e- mu+ ta+ e+
define VLQ = y y~
generate p p > j WW bb bb / p tp tp~ bp bp~ x x~ z h a YWB^2==2 TWB^2==2
output -f
""")
fcard.flush()
process_dir = new_process()
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
beamEnergy = 6500.
mass=int(re.findall(r'\d+',re.findall(r'\d+LH',runArgs.jobConfig[0])[0])[0])
print mass
# allow usage of all PDF sets
os.environ['LHAPATH']=os.environ["LHAPDF_DATA_PATH"]=(os.environ['LHAPATH'].split(':')[0])+":/cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current/"
# especially 263400 = NNPDF30_lo_as_0130_nf_4
extras = { 'lhe_version' : '2.0',
'cut_decays' : 'F',
'pdlabel' : "'lhapdf'",
'lhaid' : '263400',
'bwcutoff' : '10000'
}
safe_factor = 1.1
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir), run_card_new='run_card.dat', nevts=runArgs.maxEvents * safe_factor, rand_seed=runArgs.randomSeed, beamEnergy=beamEnergy, xqcut=0., extras=extras)
if not os.access(process_dir+'/Cards/param_card.dat',os.R_OK):
print 'ERROR: Could not get param card'
#elif os.access('param_card.dat',os.R_OK):
# print 'ERROR: Old run card in the current directory. Dont want to clobber it. Please move it first.'
else:
oldcard = open(process_dir+'/Cards/param_card.dat','r')
newcard = open('param_card.dat','w')
for line in oldcard:
if ' # MX ' in line:
newcard.write(' 6000005 %i # MX \n'%(mass))
elif ' # MTP ' in line:
newcard.write(' 6000006 %i # MTP \n'%(mass))
elif ' # MBP ' in line:
newcard.write(' 6000007 %i # MBP \n'%(mass))
elif ' # MY ' in line:
newcard.write(' 6000008 %i # MY \n'%(mass))
elif ' # KX ' in line:
newcard.write(' 1 0.4 # KX \n')
elif ' # KT ' in line:
newcard.write(' 2 0.4 # KT \n')
elif ' # KB ' in line:
newcard.write(' 3 0.4 # KB \n')
elif ' # KY ' in line:
newcard.write(' 4 0.4 # KY \n')
elif ' # WX ' in line:
newcard.write('DECAY 6000005 auto # WX \n')
elif ' # WTP ' in line:
newcard.write('DECAY 6000006 auto # WTP \n')
elif ' # WBP ' in line:
newcard.write('DECAY 6000007 auto # WBP \n')
elif ' # WY ' in line:
newcard.write('DECAY 6000008 auto # WY \n')
elif ' # zetaXuL ' in line:
newcard.write(' 1 0.000000e-01 # zetaXuL \n')
elif ' # zetaXcL ' in line:
newcard.write(' 2 0.000000e-01 # zetaXcL \n')
elif ' # zetaXtL ' in line:
newcard.write(' 3 1.000000e+00 # zetaXtL \n')
elif ' # zetaTuL ' in line:
newcard.write(' 4 0.000000e-01 # zetaTuL \n')
elif ' # zetaTcL ' in line:
newcard.write(' 5 0.000000e-01 # zetaTcL \n')
elif ' # zetaTtL ' in line:
newcard.write(' 6 1.000000e+00 # zetaTtL \n')
elif ' # zetaBdL ' in line:
newcard.write(' 7 0.000000e-01 # zetaBdL \n')
elif ' # zetaBsL ' in line:
newcard.write(' 8 0.000000e-01 # zetaBsL \n')
elif ' # zetaBbL ' in line:
newcard.write(' 9 1.000000e+00 # zetaBdL \n')
elif ' # zetaYdL ' in line:
newcard.write(' 10 0.000000e+00# zetaYdL \n')
elif ' # zetaYsL ' in line:
newcard.write(' 11 0.000000e+00 # zetaYsL \n')
elif ' # zetaYbL ' in line:
newcard.write(' 12 1.000000e+00 # zetaYbL \n')
elif ' # zetaXuR ' in line:
newcard.write(' 13 0.000000e+00 # zetaXuR \n')
elif ' # zetaXcR ' in line:
newcard.write(' 14 0.000000e+00 # zetaXcR \n')
elif ' # zetaXtR ' in line:
newcard.write(' 15 0.000000e+00 # zetaXtR \n')
elif ' # zetaTuR ' in line:
newcard.write(' 16 0.000000e+00 # zetaTuR \n')
elif ' # zetaTcR ' in line:
newcard.write(' 17 0.000000e+00 # zetaTcR \n')
elif ' # zetaTtR ' in line:
newcard.write(' 18 0.000000e+00 # zetaTtR \n')
elif ' # zetaBdR ' in line:
newcard.write(' 19 0.000000e+00 # zetaBdR \n')
elif ' # zetaBsR ' in line:
newcard.write(' 20 0.000000e+00 # zetaBsR \n')
elif ' # zetaBbR ' in line:
newcard.write(' 21 0.000000e+00 # zetaBdR \n')
elif ' # zetaYdR ' in line:
newcard.write(' 22 0.000000e+00 # zetaYdR \n')
elif ' # zetaYsR ' in line:
newcard.write(' 23 0.000000e+00 # zetaYsR \n')
elif ' # zetaYbR ' in line:
newcard.write(' 24 0.000000e+00 # zetaYbR \n')
# set some sm parameters to atlas defaults
elif ' # MB ' in line:
newcard.write(' 5 4.950000e+00 # MB \n')
elif ' # MT ' in line:
newcard.write(' 6 1.725000e+02 # MT \n')
elif ' # WT ' in line:
newcard.write('DECAY 6 1.350000e+00 # WT \n')
elif ' # WZ ' in line:
newcard.write('DECAY 23 2.495200e+00 # WZ \n')
elif ' # WW ' in line:
newcard.write('DECAY 24 2.085000e+00 # WW \n')
else:
newcard.write(line)
oldcard.close()
newcard.close()
runName='run_01'
generate(run_card_loc='./run_card.dat',param_card_loc='./param_card.dat',run_name=runName,proc_dir=process_dir)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
evgenConfig.description = "MadGraph+Pythia8 production JO with NNPDF30LN and A15NNPDF23LO for VLQ single Y to Wb with Y produced via W"
evgenConfig.keywords = ["BSM", "BSMtop", "exotic"]
evgenConfig.process = "Y_Wb"
evgenConfig.contact = ['fschenck@cern.ch']
evgenConfig.minevents = 5000
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
from MadGraphControl.MadGraphUtils import *
import re
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
import model VLQ_UFO_for_interference
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define bb = b b~
define WW = w+ w-
define tt = t t~
define ferm = ve vm vt ve~ vm~ vt~ mu- ta- e- mu+ ta+ e+
define VLQ = y y~
generate p p > j WW bb bb / p tp tp~ bp bp~ x x~ tp y y~ z h a TWB^2==4
output -f
""")
fcard.flush()
process_dir = new_process()
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
beamEnergy = 6500.
mass=int(re.findall(r'\d+',re.findall(r'\d+LH',runArgs.jobConfig[0])[0])[0])
print mass
# allow usage of all PDF sets
os.environ['LHAPATH']=os.environ["LHAPDF_DATA_PATH"]=(os.environ['LHAPATH'].split(':')[0])+":/cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current/"
# especially 263400 = NNPDF30_lo_as_0130_nf_4
extras = { 'lhe_version' : '2.0',
'cut_decays' : 'F',
'pdlabel' : "'lhapdf'",
'lhaid' : '263400',
'bwcutoff' : '10000'
}
safe_factor = 5.0
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir), run_card_new='run_card.dat', nevts=runArgs.maxEvents * safe_factor, rand_seed=runArgs.randomSeed, beamEnergy=beamEnergy, xqcut=0., extras=extras)
if not os.access(process_dir+'/Cards/param_card.dat',os.R_OK):
print 'ERROR: Could not get param card'
#elif os.access('param_card.dat',os.R_OK):
# print 'ERROR: Old run card in the current directory. Dont want to clobber it. Please move it first.'
else:
oldcard = open(process_dir+'/Cards/param_card.dat','r')
newcard = open('param_card.dat','w')
for line in oldcard:
if ' # MX ' in line:
newcard.write(' 6000005 %i # MX \n'%(mass))
elif ' # MTP ' in line:
newcard.write(' 6000006 %i # MTP \n'%(mass))
elif ' # MBP ' in line:
newcard.write(' 6000007 %i # MBP \n'%(mass))
elif ' # MY ' in line:
newcard.write(' 6000008 %i # MY \n'%(mass))
elif ' # KX ' in line:
newcard.write(' 1 0.4 # KX \n')
elif ' # KT ' in line:
newcard.write(' 2 0.4 # KT \n')
elif ' # KB ' in line:
newcard.write(' 3 0.4 # KB \n')
elif ' # KY ' in line:
newcard.write(' 4 0.4 # KY \n')
elif ' # WX ' in line:
newcard.write('DECAY 6000005 auto # WX \n')
elif ' # WTP ' in line:
newcard.write('DECAY 6000006 auto # WTP \n')
elif ' # WBP ' in line:
newcard.write('DECAY 6000007 auto # WBP \n')
elif ' # WY ' in line:
newcard.write('DECAY 6000008 auto # WY \n')
elif ' # zetaXuL ' in line:
newcard.write(' 1 0.000000e-01 # zetaXuL \n')
elif ' # zetaXcL ' in line:
newcard.write(' 2 0.000000e-01 # zetaXcL \n')
elif ' # zetaXtL ' in line:
newcard.write(' 3 1.000000e+00 # zetaXtL \n')
elif ' # zetaTuL ' in line:
newcard.write(' 4 0.000000e-01 # zetaTuL \n')
elif ' # zetaTcL ' in line:
newcard.write(' 5 0.000000e-01 # zetaTcL \n')
elif ' # zetaTtL ' in line:
newcard.write(' 6 1.000000e+00 # zetaTtL \n')
elif ' # zetaBdL ' in line:
newcard.write(' 7 0.000000e-01 # zetaBdL \n')
elif ' # zetaBsL ' in line:
newcard.write(' 8 0.000000e-01 # zetaBsL \n')
elif ' # zetaBbL ' in line:
newcard.write(' 9 1.000000e+00 # zetaBdL \n')
elif ' # zetaYdL ' in line:
newcard.write(' 10 0.000000e+00# zetaYdL \n')
elif ' # zetaYsL ' in line:
newcard.write(' 11 0.000000e+00 # zetaYsL \n')
elif ' # zetaYbL ' in line:
newcard.write(' 12 1.000000e+00 # zetaYbL \n')
elif ' # zetaXuR ' in line:
newcard.write(' 13 0.000000e+00 # zetaXuR \n')
elif ' # zetaXcR ' in line:
newcard.write(' 14 0.000000e+00 # zetaXcR \n')
elif ' # zetaXtR ' in line:
newcard.write(' 15 0.000000e+00 # zetaXtR \n')
elif ' # zetaTuR ' in line:
newcard.write(' 16 0.000000e+00 # zetaTuR \n')
elif ' # zetaTcR ' in line:
newcard.write(' 17 0.000000e+00 # zetaTcR \n')
elif ' # zetaTtR ' in line:
newcard.write(' 18 0.000000e+00 # zetaTtR \n')
elif ' # zetaBdR ' in line:
newcard.write(' 19 0.000000e+00 # zetaBdR \n')
elif ' # zetaBsR ' in line:
newcard.write(' 20 0.000000e+00 # zetaBsR \n')
elif ' # zetaBbR ' in line:
newcard.write(' 21 0.000000e+00 # zetaBdR \n')
elif ' # zetaYdR ' in line:
newcard.write(' 22 0.000000e+00 # zetaYdR \n')
elif ' # zetaYsR ' in line:
newcard.write(' 23 0.000000e+00 # zetaYsR \n')
elif ' # zetaYbR ' in line:
newcard.write(' 24 0.000000e+00 # zetaYbR \n')
# set some sm parameters to atlas defaults
elif ' # MB ' in line:
newcard.write(' 5 4.950000e+00 # MB \n')
elif ' # MT ' in line:
newcard.write(' 6 1.725000e+02 # MT \n')
elif ' # WT ' in line:
newcard.write('DECAY 6 1.350000e+00 # WT \n')
elif ' # WZ ' in line:
newcard.write('DECAY 23 2.495200e+00 # WZ \n')
elif ' # WW ' in line:
newcard.write('DECAY 24 2.085000e+00 # WW \n')
else:
newcard.write(line)
oldcard.close()
newcard.close()
runName='run_01'
generate(run_card_loc='./run_card.dat',param_card_loc='./param_card.dat',run_name=runName,proc_dir=process_dir)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
# Configure the HT filter
include('MC15JobOptions/HTFilter.py')
filtSeq.HTFilter.MinJetPt = 35.*GeV # Min pT to consider jet in HT
filtSeq.HTFilter.MaxJetEta = 2.5 # Max eta to consider jet in HT
filtSeq.HTFilter.MinHT = 250.*GeV # Min HT to keep event
filtSeq.HTFilter.MaxHT = 9000.*GeV # Max HT to keep event
filtSeq.HTFilter.UseLeptonsFromWZTau = True # Include e/mu from the MC event in the HT
filtSeq.HTFilter.MinLeptonPt = 25.*GeV # Min pT to consider muon in HT
filtSeq.HTFilter.MaxLeptonEta = 2.5 # Max eta to consider muon in HT
evgenConfig.description = "MadGraph+Pythia8 production JO with NNPDF30LN and A15NNPDF23LO for VLQ single Y to Wb with Y produced via W"
evgenConfig.keywords = ["BSM", "BSMtop", "exotic"]
evgenConfig.process = "Y_Wb"
evgenConfig.contact = ['fschenck@cern.ch']
evgenConfig.minevents = 5000
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
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