Commit eeac213e authored by mcjo's avatar mcjo
Browse files

add 410644-5, 410436, 410642-3, 345929


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@812961 4525493e-7705-40b1-a816-d608a930855b
parent 130e0c1f
2017-02-09 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-09-06
* add 410644-5, 410436, 410642-3, 345929
2017-02-07 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-09-05
* add 410633-7, update 309708, 395244-58, replace 345873
......
from MadGraphControl.MadGraphUtils import *
import fileinput
import shutil
import subprocess
import os
import sys
# General settings
minevents=50000
nevents=55000
#mode=0
# MG merging settings
maxjetflavor=4
ickkw=0
dyn_scale = '10'
### DSID lists (extensions can include systematics samples)
tZ =[510049]
systScalefactUp=[510153]
systScalefactDn=[510154]
DSID = runArgs.runNumber
if DSID in tZ + systScalefactUp + systScalefactDn:
mgproc = """generate p p > t b~ j l+ l- $$ W+ W- [QCD] \n"""
mgproc += """add process p p > t~ b j l+ l- $$ W+ W- [QCD]"""
name='tllq'
keyword=['SM','singleTop','tZ','lepton']
mode = 1
cluster_type= 'condor'
cluster_queue= 'None'
else:
raise RuntimeError("runNumber %i not recognised in these jobOptions."%runArgs.runNumber)
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_'+str(name)
get_param_card = subprocess.Popen(['get_files','-data','aMcAtNlo_param_card_tt.dat'])
get_param_card.wait()
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
"""+mgproc+"""
output -f
""")
fcard.close()
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
# Decay with MadSpin
madspin_card_loc='madspin_card.dat'
mscard = open(madspin_card_loc,'w')
mscard.write("""#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#set seed %i
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
decay t > w+ b, w+ > l+ vl
decay t~ > w- b~, w- > l- vl~
decay z > l+ l-
launch """)
mscard.close()
pdflabel="lhapdf"
lhaid=260400
#Fetch default LO run_card.dat and set parameters
extras = { 'lhe_version' : '3.0',
'pdlabel' : "'"+pdflabel+"'",
'lhaid' : lhaid,
'maxjetflavor' : maxjetflavor,
'ickkw' : 0,
'ptj' : 1,
'drll' : 0.0,
'ptl' : 1,
'etal' : 5,
'etaj' : 5,
'bwcutoff' : 50,
'dynamical_scale_choice' : dyn_scale,
'reweight_scale': 'True',
'reweight_PDF': 'True',
'rw_rscale': "1.0, 0.5, 2.0",
'rw_fscale': "1.0, 0.5, 2.0",
'PDF_set_min': 260401,
'PDF_set_max':260500,
'store_rwgt_info':True
}
scalefact = 1.0
if runArgs.runNumber in systScalefactUp:
name=name+'_facsc2_rensc2'
scalefact = 2.0
elif runArgs.runNumber in systScalefactDn:
name=name+'_facsc0p5_rensc0p5'
scalefact = 0.5
process_dir = new_process()
# Cook the setscales file for the user defined dynamical scale
fileN = process_dir+'/SubProcesses/setscales.f'
mark = ' elseif(dynamical_scale_choice.eq.10) then'
rmLines = ['ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc',
'cc USER-DEFINED SCALE: ENTER YOUR CODE HERE cc',
'cc to use this code you must set cc',
'cc dynamical_scale_choice = 10 cc',
'cc in the run_card (run_card.dat) cc',
'ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc',
'write(*,*) "User-defined scale not set"',
'stop 1',
'temp_scale_id=\'User-defined dynamical scale\' ! use a meaningful string',
'tmp = 0d0']
dyn_scale_fact = 1.0
flag=0
for line in fileinput.input(fileN, inplace=1):
toKeep = True
for rmLine in rmLines:
if line.find(rmLine) >= 0:
toKeep = False
break
if toKeep:
print line,
if line.startswith(mark) and flag==0:
flag +=1
print """
c Q^2= mb^2 + 0.5*(pt^2+ptbar^2)
c rscale=0d0 !factorization scale**2 for pdf1
c !factorization scale**2 for pdf2
c xmtc=dot(P(0,6),P(0,6))
c rscale = 4d0 * sqrt(pt(P(0,6))**2+xmtc)
do i=3,nexternal
xm2=dot(pp(0,i),pp(0,i))
if ( xm2 < 30 .and. xm2 > 10 ) then
tmp = 4d0 * dsqrt(pt(pp(0,i))**2+xm2)
c write(*,*) i, pt(P(0,i))**2, xmtc, sqrt(pt(P(0,i))**2+xmtc), rscale
c write(*,*) i, xmtc, pt(P(0,i)), rscale
endif
enddo
c it can be 4 or 0.25
tmp = %f*tmp""" % dyn_scale_fact
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,xqcut=0.,
extras=extras,scalefact=scalefact)
print_cards()
generate(run_card_loc='run_card.dat',
param_card_loc='aMcAtNlo_param_card_tt.dat',
madspin_card_loc=madspin_card_loc,
mode=mode,
cluster_type=cluster_type,
cluster_queue=cluster_queue,
proc_dir=process_dir)
arrange_output(proc_dir=process_dir,outputDS=stringy+'._00001.events.tar.gz')
## Shower
evgenConfig.description = 'MadGraph_'+str(name)
evgenConfig.keywords+=keyword
evgenConfig.inputfilecheck = stringy
evgenConfig.minevents = minevents
runArgs.inputGeneratorFile=stringy+'._00001.events.tar.gz'
#--------------------------------------------------------------
# POWHEG+Herwig7 gg->H+Z->vvbbar production
#--------------------------------------------------------------
#--------------------------------------------------------------
# Herwig7 showering
#--------------------------------------------------------------
include("MC15JobOptions/Herwig7_701_H7UE_MMHT2014lo68cl_NNPDF3ME_LHEF_EvtGen_Common.py")
include("MC15JobOptions/Herwig7_701_StripWeights.py")
from Herwig7_i import config as hw
genSeq.Herwig7.Commands += hw.powhegbox_cmds().splitlines()
## only consider H->bb decays
genSeq.Herwig7.Commands += [
'## force H->bb decays',
'do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;',
'do /Herwig/Particles/h0:PrintDecayModes' # print out decays modes and branching ratios to the terminal/log.generate
]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = "POWHEG+Herwig7 gg->H+Z->vvbbar production"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs", "mH125" , "ZHiggs" ]
evgenConfig.contact = [ 'carlo.enrico.pandini@cern.ch', 'paolo.francavilla@cern.ch' ]
evgenConfig.generators += [ 'Powheg', 'Herwig7' ]
evgenConfig.inputfilecheck = "TXT"
evgenConfig.process = "gg->ZH, H->bb, Z->vv"
evgenConfig.description = 'MG+Pythia8 NLO single top + Z production with Pythia A14 tune'
evgenConfig.keywords = [ 'SM','singleTop','tZ','lepton','NLO']
evgenConfig.contact = ['dylan.frizzell@cern.ch' ]
#Common MG_Control file
include("MC15JobOptions/MadGraphControl_tllq_NLO_Pythia8_A14.py")
#Below is showering
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
include("MC15JobOptions/Pythia8_ShowerWeights.py")
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8+EvtGen single-top-quark t-channel (2->3) production (top),MadSpin, A14 tune, ME NNPDF3.04f NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'lepton']
evgenConfig.contact = [ 'dominic.hirschbuehl@cern.ch']
evgenConfig.minevents = 500
include('PowhegControl/PowhegControl_t_tch_4FS_Common.py')
PowhegConfig.decay_mode = "t > undecayed"
PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R = [1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5] # List of renormalisation scales
PowhegConfig.PDF = [260400, 25410, 13191, 92000] # NNPDF30, MMHT, CT14, PDF4LHC - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260401, 260501)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(92001 , 92031 )) # Include the PDF4LHC error set
# Information on how to run with multiple weights: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/PowhegForATLAS#Running_with_multiple_scale_PDF
# PDFs - you can see a listing here: https://lhapdf.hepforge.org/pdfsets.html; picked these three as they are the inputs to the PDF4LHC2015 prescription (http://arxiv.org/pdf/1510.03865v2.pdf).
# Define a weight group configuration for scale variations with different PDFs
# Nominal mu_F = mu_R = 1.0 is not required as this is captured by the PDF variation above
PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25410] )
# Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13191] )
# Scale variations, PDF4LHC15_nlo_30
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15', parameter_values=[ 2.0, 1.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15', parameter_values=[ 0.5, 1.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15', parameter_values=[ 1.0, 2.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15', parameter_values=[ 1.0, 0.5, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15', parameter_values=[ 0.5, 0.5, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15', parameter_values=[ 2.0, 2.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15', parameter_values=[ 0.5, 2.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15', parameter_values=[ 2.0, 0.5, 92000] )
PowhegConfig.nEvents *= 1.1 # Safety factor
PowhegConfig.MadSpin_decays= ["decay t > w+ b, w+ > l+ vl"]
PowhegConfig.MadSpin_process = "generate p p > t b~ j $$ w+ w- [QCD]"
# Parameters are hardcoded to avoid madspin picking up it's own defaults.
PowhegConfig.alphaem_inv = 1.323489e+02
PowhegConfig.G_F = 1.166370e-05
PowhegConfig.alphaqcd = 1.184000e-01
PowhegConfig.mass_W = 80.399
PowhegConfig.mass_Z = 9.118760e+01
PowhegConfig.mass_H = 125.0
# use diagonal CKM matrix
PowhegConfig.CKM_Vud = 1.0
PowhegConfig.CKM_Vus = 0.000001
PowhegConfig.CKM_Vub = 0.000001
PowhegConfig.CKM_Vcd = 0.000001
PowhegConfig.CKM_Vcs = 1.0
PowhegConfig.CKM_Vcb = 0.000001
PowhegConfig.CKM_Vtd = 0.000001
PowhegConfig.CKM_Vts = 0.000001
PowhegConfig.CKM_Vtb = 1.0
# Initial settings
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering with the A14 NNPDF2.3 tune
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8+EvtGen single-top-quark t-channel (2->3) production (top),MadSpin, A14 tune, ME NNPDF3.04f NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'lepton']
evgenConfig.contact = [ 'dominic.hirschbuehl@cern.ch']
evgenConfig.minevents = 500
include('PowhegControl/PowhegControl_t_tch_4FS_Common.py')
PowhegConfig.decay_mode = "t~ > undecayed"
PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R = [1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5] # List of renormalisation scales
PowhegConfig.PDF = [260400, 25410, 13191, 92000] # NNPDF30, MMHT, CT14, PDF4LHC - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260401, 260501)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(92001 , 92031 )) # Include the PDF4LHC error set
# Information on how to run with multiple weights: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/PowhegForATLAS#Running_with_multiple_scale_PDF
# PDFs - you can see a listing here: https://lhapdf.hepforge.org/pdfsets.html; picked these three as they are the inputs to the PDF4LHC2015 prescription (http://arxiv.org/pdf/1510.03865v2.pdf).
# Define a weight group configuration for scale variations with different PDFs
# Nominal mu_F = mu_R = 1.0 is not required as this is captured by the PDF variation above
PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25410] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25410] )
# Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13191] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13191] )
# Scale variations, PDF4LHC15_nlo_30
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15', parameter_values=[ 2.0, 1.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15', parameter_values=[ 0.5, 1.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15', parameter_values=[ 1.0, 2.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15', parameter_values=[ 1.0, 0.5, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15', parameter_values=[ 0.5, 0.5, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15', parameter_values=[ 2.0, 2.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15', parameter_values=[ 0.5, 2.0, 92000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15', parameter_values=[ 2.0, 0.5, 92000] )
PowhegConfig.nEvents *= 1.1 # Safety factor
PowhegConfig.MadSpin_decays= ["decay t~ > w- b~, w- > l- vl~"]
PowhegConfig.MadSpin_process = "generate p p > t~ b j $$ w+ w- [QCD]"
# Parameters are hardcoded to avoid madspin picking up it's own defaults.
PowhegConfig.alphaem_inv = 1.323489e+02
PowhegConfig.G_F = 1.166370e-05
PowhegConfig.alphaqcd = 1.184000e-01
PowhegConfig.mass_W = 80.399
PowhegConfig.mass_Z = 9.118760e+01
PowhegConfig.mass_H = 125.0
# use diagonal CKM matrix
PowhegConfig.CKM_Vud = 1.0
PowhegConfig.CKM_Vus = 0.000001
PowhegConfig.CKM_Vub = 0.000001
PowhegConfig.CKM_Vcd = 0.000001
PowhegConfig.CKM_Vcs = 1.0
PowhegConfig.CKM_Vcb = 0.000001
PowhegConfig.CKM_Vtd = 0.000001
PowhegConfig.CKM_Vts = 0.000001
PowhegConfig.CKM_Vtb = 1.0
# Initial settings
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering with the A14 NNPDF2.3 tune
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8+EvtGen single-top-quark s-channel production (top), inclusive, A14 tune, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'top', 'singleTop', 'sChannel' ]
evgenConfig.contact = [ 'cescobar@cern.ch']
evgenConfig.generators += [ 'Powheg' ]
evgenConfig.minevents = 1000
#--------------------------------------------------------------
# Powheg single-top s-channel setup
#--------------------------------------------------------------
include('PowhegControl/PowhegControl_t_sch_Common.py')
PowhegConfig.decay_mode = 't > b l+ vl'
PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R = [1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5] # List of renormalisation scales
PowhegConfig.PDF = [260000, 25200, 13165, 90900] # NNPDF30, MMHT, CT14, PDF4LHC - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260001, 260101)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(90901 , 90931 )) # Include the PDF4LHC error set
# Information on how to run with multiple weights: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/PowhegForATLAS#Running_with_multiple_scale_PDF
# PDFs - you can see a listing here: https://lhapdf.hepforge.org/pdfsets.html; picked these three as they are the inputs to the PDF4LHC2015 prescription (http://arxiv.org/pdf/1510.03865v2.pdf).
# Define a weight group configuration for scale variations with different PDFs
# Nominal mu_F = mu_R = 1.0 is not required as this is captured by the PDF variation above
PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25200] )
# Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13165] )
# Scale variations, PDF4LHC15_nlo_30
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15', parameter_values=[ 2.0, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15', parameter_values=[ 0.5, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15', parameter_values=[ 1.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15', parameter_values=[ 1.0, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15', parameter_values=[ 0.5, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15', parameter_values=[ 2.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15', parameter_values=[ 0.5, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15', parameter_values=[ 2.0, 0.5, 90900] )
PowhegConfig.nEvents *= 1.1 # Safety factor
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering with the A14 NNPDF2.3 tune
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8+EvtGen single-top-quark s-channel production (anti-top), inclusive, A14 tune, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'top', 'singleTop', 'sChannel' ]
evgenConfig.contact = [ 'cescobar@cern.ch']
evgenConfig.generators += [ 'Powheg' ]
evgenConfig.minevents = 1000
#--------------------------------------------------------------
# Powheg single-top s-channel setup
#--------------------------------------------------------------
include('PowhegControl/PowhegControl_t_sch_Common.py')
PowhegConfig.decay_mode = 't~ > b~ l- vl~'
PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R = [1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5] # List of renormalisation scales
PowhegConfig.PDF = [260000, 25200, 13165, 90900] # NNPDF30, MMHT, CT14, PDF4LHC - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260001, 260101)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(90901 , 90931 )) # Include the PDF4LHC error set
# Information on how to run with multiple weights: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/PowhegForATLAS#Running_with_multiple_scale_PDF
# PDFs - you can see a listing here: https://lhapdf.hepforge.org/pdfsets.html; picked these three as they are the inputs to the PDF4LHC2015 prescription (http://arxiv.org/pdf/1510.03865v2.pdf).