Commit 553af14b authored by Jan Kretzschmar's avatar Jan Kretzschmar Committed by Christian Gutschow
Browse files

Powheg jj LHE optimised v2

parent 657079bb
Pipeline #2978453 passed with stages
in 1 minute and 3 seconds
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = "POWHEG+Pythia8 dijet production with A14 NNPDF2.3 tune. v2 with tighter bornktmin cut and extreme folding"
evgenConfig.keywords = ["SM", "QCD", "jets", "2jet"]
evgenConfig.contact = ["amoroso@cern.ch", "jan.kretzschmar@cern.ch"]
evgenConfig.nEventsPerJob = 100
evgenConfig.generators = ["Powheg"]
# --------------------------------------------------------------
# Load ATLAS defaults for the Powheg jj process
# --------------------------------------------------------------
include("PowhegControl/PowhegControl_jj_Common.py")
# raise bornktmin to kill badly behaved events with low-pT and high/negative weight
PowhegConfig.bornktmin =10
PowhegConfig.bornsuppfact =60
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 = [14000, 304400, 25300, 14200, 14100, 14300, 14400, 25100, 91500, 42560, 317500, 303600, 319300, 319500, 322500, 322700, 322900, 323100, 323300, 323500, 323700, 323900, 27400]
PowhegConfig.rwl_group_events = 10000
PowhegConfig.ncall1 = 160000
PowhegConfig.ncall2 = 160000
PowhegConfig.nubound = 9600000
PowhegConfig.itmx1 = 20
PowhegConfig.itmx2 = 20
### Extreme fold parameter to kill all negative weights
PowhegConfig.foldcsi = 10
PowhegConfig.foldphi = 10
PowhegConfig.foldy = 10
PowhegConfig.nEvents = runArgs.maxEvents if runArgs.maxEvents>0 else evgenConfig.nEventsPerJob
# PDF variations with nominal scale variation:
# CT18NNLO,
# NNPDF31_nnlo_as_0118_hessian,
# MMHT2014nnlo68cl,
# CT18ANNLO,
# CT18ZNNLO,
# CT18XNNLO,
# CT18NLO,
# MMHT2014nlo68cl,
# PDF4LHC15_nnlo_mc,
# ABMP16_5_nnlo,
# NNPDF31_nnlo_as_0118_nojets,
# NNPDF31_nnlo_as_0118
# NNPDF31_nnlo_as_0116
# NNPDF31_nnlo_as_0120
# NNPDF31_nnlo_as_0108
# NNPDF31_nnlo_as_0110
# NNPDF31_nnlo_as_0112
# NNPDF31_nnlo_as_0114
# NNPDF31_nnlo_as_0117
# NNPDF31_nnlo_as_0119
# NNPDF31_nnlo_as_0122
# NNPDF31_nnlo_as_0124
# MSHT20nnlo_as118
PowhegConfig.PDF.extend(range(14001, 14059)) # Include the CT18NNLO error set
PowhegConfig.PDF.extend(range(304401, 304501)) # Include the NNPDF31_nnlo_as_0118_hessian error set
#PowhegConfig.PDF.extend(range(25301, 25351)) # Include the MMHT2014nnlo68cl error set
PowhegConfig.PDF.extend(range(27401, 27465)) # Include the MSHT20nnlo_as118 error set
# --------------------------------------------------------------
# Generate events
# --------------------------------------------------------------
PowhegConfig.generate()
../../950xxx/950330/mc_13TeV.Ph_jj_JZ1_cutfold.GRID.tar.gz
\ No newline at end of file
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = "POWHEG+Pythia8 dijet production with A14 NNPDF2.3 tune. v2 with tighter bornktmin cut"
evgenConfig.keywords = ["SM", "QCD", "jets", "2jet"]
evgenConfig.contact = ["amoroso@cern.ch", "jan.kretzschmar@cern.ch"]
evgenConfig.nEventsPerJob = 2000
evgenConfig.generators = ["Powheg"]
# --------------------------------------------------------------
# Load ATLAS defaults for the Powheg jj process
# --------------------------------------------------------------
include("PowhegControl/PowhegControl_jj_Common.py")
# raise bornktmin to kill badly behaved events with low-pT and high/negative weight
PowhegConfig.bornktmin =30
PowhegConfig.bornsuppfact =160
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 = [14000, 304400, 25300, 14200, 14100, 14300, 14400, 25100, 91500, 42560, 317500, 303600, 319300, 319500, 322500, 322700, 322900, 323100, 323300, 323500, 323700, 323900, 27400]
PowhegConfig.rwl_group_events = 10000
PowhegConfig.ncall1 = 160000
PowhegConfig.ncall2 = 160000
PowhegConfig.nubound = 9600000
PowhegConfig.itmx1 = 20
PowhegConfig.itmx2 = 20
### Fold parameter reducing the negative eventweight fraction
PowhegConfig.foldcsi = 2
PowhegConfig.foldphi = 2
PowhegConfig.foldy = 5
PowhegConfig.nEvents = runArgs.maxEvents if runArgs.maxEvents>0 else evgenConfig.nEventsPerJob
# PDF variations with nominal scale variation:
# CT18NNLO,
# NNPDF31_nnlo_as_0118_hessian,
# MMHT2014nnlo68cl,
# CT18ANNLO,
# CT18ZNNLO,
# CT18XNNLO,
# CT18NLO,
# MMHT2014nlo68cl,
# PDF4LHC15_nnlo_mc,
# ABMP16_5_nnlo,
# NNPDF31_nnlo_as_0118_nojets,
# NNPDF31_nnlo_as_0118
# NNPDF31_nnlo_as_0116
# NNPDF31_nnlo_as_0120
# NNPDF31_nnlo_as_0108
# NNPDF31_nnlo_as_0110
# NNPDF31_nnlo_as_0112
# NNPDF31_nnlo_as_0114
# NNPDF31_nnlo_as_0117
# NNPDF31_nnlo_as_0119
# NNPDF31_nnlo_as_0122
# NNPDF31_nnlo_as_0124
# MSHT20nnlo_as118
PowhegConfig.PDF.extend(range(14001, 14059)) # Include the CT18NNLO error set
PowhegConfig.PDF.extend(range(304401, 304501)) # Include the NNPDF31_nnlo_as_0118_hessian error set
#PowhegConfig.PDF.extend(range(25301, 25351)) # Include the MMHT2014nnlo68cl error set
PowhegConfig.PDF.extend(range(27401, 27465)) # Include the MSHT20nnlo_as118 error set
# --------------------------------------------------------------
# Generate events
# --------------------------------------------------------------
PowhegConfig.generate()
../../950xxx/950331/mc_13TeV.Ph_jj_JZ2_cut.GRID.tar.gz
\ No newline at end of file
Markdown is supported
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