From 291ea228336bc6191f6b5442c5b477a09bafd169 Mon Sep 17 00:00:00 2001
From: Yiming Abulaiti <yiming.abulaiti@cern.ch>
Date: Tue, 23 Apr 2024 14:46:07 +0200
Subject: [PATCH 1/2] MC20, Tau+X: VLQ decay into tautau Signal Request in R22

---
 603xxx/603472/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL200.py | 115 ++++++++++++++++++
 603xxx/603473/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL500.py | 115 ++++++++++++++++++
 603xxx/603474/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL1000.py | 115 ++++++++++++++++++
 603xxx/603475/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL2000.py | 115 ++++++++++++++++++
 603xxx/603476/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL200.py | 115 ++++++++++++++++++
 603xxx/603477/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL500.py | 115 ++++++++++++++++++
 603xxx/603478/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL1000.py | 115 ++++++++++++++++++
 603xxx/603479/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL2000.py | 115 ++++++++++++++++++
 603xxx/603480/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL200.py | 115 ++++++++++++++++++
 603xxx/603481/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL500.py | 115 ++++++++++++++++++
 603xxx/603482/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL1000.py | 115 ++++++++++++++++++
 603xxx/603483/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL2000.py | 115 ++++++++++++++++++
 603xxx/603484/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL200.py | 115 ++++++++++++++++++
 603xxx/603485/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL500.py | 115 ++++++++++++++++++
 603xxx/603486/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL1000.py | 115 ++++++++++++++++++
 603xxx/603487/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL2000.py | 115 ++++++++++++++++++
 603xxx/603488/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL200.py | 115 ++++++++++++++++++
 603xxx/603489/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL500.py | 115 ++++++++++++++++++
 603xxx/603490/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL1000.py | 115 ++++++++++++++++++
 603xxx/603491/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL2000.py | 115 ++++++++++++++++++
 603xxx/603492/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL200.py | 115 ++++++++++++++++++
 603xxx/603493/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL500.py | 115 ++++++++++++++++++
 603xxx/603494/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL1000.py | 115 ++++++++++++++++++
 603xxx/603495/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL2000.py | 115 ++++++++++++++++++
 603xxx/603496/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL200.py | 115 ++++++++++++++++++
 603xxx/603497/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL500.py | 115 ++++++++++++++++++
 603xxx/603498/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL1000.py | 115 ++++++++++++++++++
 603xxx/603499/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL2000.py | 115 ++++++++++++++++++
 603xxx/603500/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL200.py | 115 ++++++++++++++++++
 603xxx/603501/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL500.py | 115 ++++++++++++++++++
 603xxx/603502/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL1000.py | 115 ++++++++++++++++++
 603xxx/603503/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL2000.py | 115 ++++++++++++++++++
 603xxx/603504/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL200.py | 115 ++++++++++++++++++
 603xxx/603505/log.generate.short              |   6 +
 ...4NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL500.py | 115 ++++++++++++++++++
 603xxx/603506/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL1000.py | 115 ++++++++++++++++++
 603xxx/603507/log.generate.short              |   6 +
 ...NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL2000.py | 115 ++++++++++++++++++
 72 files changed, 4356 insertions(+)
 create mode 100644 603xxx/603472/log.generate.short
 create mode 100644 603xxx/603472/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL200.py
 create mode 100644 603xxx/603473/log.generate.short
 create mode 100644 603xxx/603473/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL500.py
 create mode 100644 603xxx/603474/log.generate.short
 create mode 100644 603xxx/603474/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL1000.py
 create mode 100644 603xxx/603475/log.generate.short
 create mode 100644 603xxx/603475/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL2000.py
 create mode 100644 603xxx/603476/log.generate.short
 create mode 100644 603xxx/603476/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL200.py
 create mode 100644 603xxx/603477/log.generate.short
 create mode 100644 603xxx/603477/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL500.py
 create mode 100644 603xxx/603478/log.generate.short
 create mode 100644 603xxx/603478/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL1000.py
 create mode 100644 603xxx/603479/log.generate.short
 create mode 100644 603xxx/603479/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL2000.py
 create mode 100644 603xxx/603480/log.generate.short
 create mode 100644 603xxx/603480/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL200.py
 create mode 100644 603xxx/603481/log.generate.short
 create mode 100644 603xxx/603481/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL500.py
 create mode 100644 603xxx/603482/log.generate.short
 create mode 100644 603xxx/603482/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL1000.py
 create mode 100644 603xxx/603483/log.generate.short
 create mode 100644 603xxx/603483/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL2000.py
 create mode 100644 603xxx/603484/log.generate.short
 create mode 100644 603xxx/603484/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL200.py
 create mode 100644 603xxx/603485/log.generate.short
 create mode 100644 603xxx/603485/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL500.py
 create mode 100644 603xxx/603486/log.generate.short
 create mode 100644 603xxx/603486/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL1000.py
 create mode 100644 603xxx/603487/log.generate.short
 create mode 100644 603xxx/603487/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL2000.py
 create mode 100644 603xxx/603488/log.generate.short
 create mode 100644 603xxx/603488/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL200.py
 create mode 100644 603xxx/603489/log.generate.short
 create mode 100644 603xxx/603489/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL500.py
 create mode 100644 603xxx/603490/log.generate.short
 create mode 100644 603xxx/603490/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL1000.py
 create mode 100644 603xxx/603491/log.generate.short
 create mode 100644 603xxx/603491/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL2000.py
 create mode 100644 603xxx/603492/log.generate.short
 create mode 100644 603xxx/603492/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL200.py
 create mode 100644 603xxx/603493/log.generate.short
 create mode 100644 603xxx/603493/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL500.py
 create mode 100644 603xxx/603494/log.generate.short
 create mode 100644 603xxx/603494/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL1000.py
 create mode 100644 603xxx/603495/log.generate.short
 create mode 100644 603xxx/603495/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL2000.py
 create mode 100644 603xxx/603496/log.generate.short
 create mode 100644 603xxx/603496/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL200.py
 create mode 100644 603xxx/603497/log.generate.short
 create mode 100644 603xxx/603497/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL500.py
 create mode 100644 603xxx/603498/log.generate.short
 create mode 100644 603xxx/603498/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL1000.py
 create mode 100644 603xxx/603499/log.generate.short
 create mode 100644 603xxx/603499/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL2000.py
 create mode 100644 603xxx/603500/log.generate.short
 create mode 100644 603xxx/603500/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL200.py
 create mode 100644 603xxx/603501/log.generate.short
 create mode 100644 603xxx/603501/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL500.py
 create mode 100644 603xxx/603502/log.generate.short
 create mode 100644 603xxx/603502/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL1000.py
 create mode 100644 603xxx/603503/log.generate.short
 create mode 100644 603xxx/603503/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL2000.py
 create mode 100644 603xxx/603504/log.generate.short
 create mode 100644 603xxx/603504/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL200.py
 create mode 100644 603xxx/603505/log.generate.short
 create mode 100644 603xxx/603505/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL500.py
 create mode 100644 603xxx/603506/log.generate.short
 create mode 100644 603xxx/603506/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL1000.py
 create mode 100644 603xxx/603507/log.generate.short
 create mode 100644 603xxx/603507/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL2000.py

diff --git a/603xxx/603472/log.generate.short b/603xxx/603472/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603472/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603472/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL200.py b/603xxx/603472/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603472/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603473/log.generate.short b/603xxx/603473/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603473/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603473/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL500.py b/603xxx/603473/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603473/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603474/log.generate.short b/603xxx/603474/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603474/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603474/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL1000.py b/603xxx/603474/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603474/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603475/log.generate.short b/603xxx/603475/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603475/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603475/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL2000.py b/603xxx/603475/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603475/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_0_33R1_0_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603476/log.generate.short b/603xxx/603476/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603476/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603476/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL200.py b/603xxx/603476/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603476/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603477/log.generate.short b/603xxx/603477/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603477/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603477/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL500.py b/603xxx/603477/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603477/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603478/log.generate.short b/603xxx/603478/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603478/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603478/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL1000.py b/603xxx/603478/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603478/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603479/log.generate.short b/603xxx/603479/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603479/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603479/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL2000.py b/603xxx/603479/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603479/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_0_33R1_0_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603480/log.generate.short b/603xxx/603480/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603480/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603480/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL200.py b/603xxx/603480/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603480/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603481/log.generate.short b/603xxx/603481/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603481/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603481/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL500.py b/603xxx/603481/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603481/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603482/log.generate.short b/603xxx/603482/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603482/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603482/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL1000.py b/603xxx/603482/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603482/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603483/log.generate.short b/603xxx/603483/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603483/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603483/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL2000.py b/603xxx/603483/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603483/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_0_33R1_0_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603484/log.generate.short b/603xxx/603484/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603484/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603484/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL200.py b/603xxx/603484/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603484/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603485/log.generate.short b/603xxx/603485/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603485/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603485/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL500.py b/603xxx/603485/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603485/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603486/log.generate.short b/603xxx/603486/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603486/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603486/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL1000.py b/603xxx/603486/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603486/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603487/log.generate.short b/603xxx/603487/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603487/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603487/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL2000.py b/603xxx/603487/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603487/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L1_5_33R1_5_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603488/log.generate.short b/603xxx/603488/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603488/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603488/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL200.py b/603xxx/603488/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603488/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603489/log.generate.short b/603xxx/603489/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603489/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603489/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL500.py b/603xxx/603489/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603489/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603490/log.generate.short b/603xxx/603490/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603490/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603490/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL1000.py b/603xxx/603490/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603490/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603491/log.generate.short b/603xxx/603491/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603491/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603491/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL2000.py b/603xxx/603491/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603491/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L1_5_33R1_5_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603492/log.generate.short b/603xxx/603492/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603492/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603492/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL200.py b/603xxx/603492/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603492/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603493/log.generate.short b/603xxx/603493/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603493/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603493/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL500.py b/603xxx/603493/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603493/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603494/log.generate.short b/603xxx/603494/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603494/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603494/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL1000.py b/603xxx/603494/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603494/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603495/log.generate.short b/603xxx/603495/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603495/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603495/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL2000.py b/603xxx/603495/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603495/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L1_5_33R1_5_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603496/log.generate.short b/603xxx/603496/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603496/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603496/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL200.py b/603xxx/603496/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603496/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603497/log.generate.short b/603xxx/603497/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603497/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603497/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL500.py b/603xxx/603497/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603497/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603498/log.generate.short b/603xxx/603498/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603498/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603498/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL1000.py b/603xxx/603498/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603498/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603499/log.generate.short b/603xxx/603499/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603499/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603499/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL2000.py b/603xxx/603499/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603499/mc.PhPy8EG_A14NNPDF23_DYV_MU1000_33L2_0_33R2_0_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603500/log.generate.short b/603xxx/603500/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603500/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603500/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL200.py b/603xxx/603500/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603500/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603501/log.generate.short b/603xxx/603501/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603501/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603501/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL500.py b/603xxx/603501/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603501/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603502/log.generate.short b/603xxx/603502/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603502/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603502/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL1000.py b/603xxx/603502/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603502/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603503/log.generate.short b/603xxx/603503/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603503/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603503/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL2000.py b/603xxx/603503/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603503/mc.PhPy8EG_A14NNPDF23_DYV_MU2000_33L2_0_33R2_0_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603504/log.generate.short b/603xxx/603504/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603504/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603504/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL200.py b/603xxx/603504/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL200.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603504/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL200.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603505/log.generate.short b/603xxx/603505/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603505/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603505/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL500.py b/603xxx/603505/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL500.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603505/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL500.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603506/log.generate.short b/603xxx/603506/log.generate.short
new file mode 100644
index 0000000000..875c004959
--- /dev/null
+++ b/603xxx/603506/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603506/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL1000.py b/603xxx/603506/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL1000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603506/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL1000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
diff --git a/603xxx/603507/log.generate.short b/603xxx/603507/log.generate.short
new file mode 100644
index 0000000000..98b705f86b
--- /dev/null
+++ b/603xxx/603507/log.generate.short
@@ -0,0 +1,6 @@
+- estimated CPU for CI job = 0.00 hrs
+- using release =  AthGeneration-23.6.16
+- ecmEnergy =  13000.0
+- randomSeed =  42
+- EVNT to EVNT =  False
+- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603507/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL2000.py b/603xxx/603507/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL2000.py
new file mode 100644
index 0000000000..7571637117
--- /dev/null
+++ b/603xxx/603507/mc.PhPy8EG_A14NNPDF23_DYV_MU3000_33L2_0_33R2_0_SL2000.py
@@ -0,0 +1,115 @@
+# Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
+
+#--------------------------------------------------------------
+# This is an example joboption to generate events with Powheg
+# using ATLAS' interface. Users should optimise and carefully
+# validate the settings before making an official sample request.
+#--------------------------------------------------------------
+import os
+import re
+from MCJobOptionUtils.JOsupport import get_physics_short
+
+#get JO name
+job_option_name = get_physics_short()
+print(job_option_name)
+
+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = "POWHEG+Pythia8 Drell-Yan Vector LeptoQuark production with A14 NNPDF2.3 tune."
+evgenConfig.keywords = ["BSM", "leptoquark"]
+evgenConfig.contact = ["tpelzer@cern.ch"]
+
+# --------------------------------------------------------------
+# Load ATLAS defaults for the Powheg Drell-Yan Vector LeptoQuark process
+# --------------------------------------------------------------
+include("PowhegControl/PowhegControl_DY_VLQ_Common.py")
+
+# --------------------------------------------------------------
+# Relevant parameters for this process
+# --------------------------------------------------------------
+# 4321 Model Parameters
+PowhegConfig.g4 = 1.414214 # SU(4) coupling strength g4
+
+# Get betaL33 from jo file
+matches_L33 = re.search("33L([0-9]_[0-9]+).*", job_option_name)
+if matches_L33 is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    L33=0.0
+else:
+    L33= float(matches_L33.group(1).replace("_", "."))
+PowhegConfig.betaL3x3 = L33 # Relative coupling strength to left-handed fermions (3rd generation)
+
+# Get betaR33 from jo file
+matches_R33 = re.search("33R([0-9]_[0-9]+).*", job_option_name)
+if matches_R33 is None:
+    print("Cannot find betaR 33 string in job option name: {:s}.".format(job_option_name))
+    R33=0.0
+else:
+    R33= float(matches_R33.group(1).replace("_", "."))
+PowhegConfig.betaR3x3 = R33 # Relative coupling strength to right-handed fermions (3rd generation)
+
+# Get MU from jo file
+matches_MU = re.search("MU([0-9]+).*", job_option_name)
+if matches_MU is None:
+    print("Cannot find MU string in job option name: {:s}.".format(job_option_name))
+    MU=0.0
+else:
+    MU= float(matches_MU.group(1).replace("_", "."))
+PowhegConfig.MU1 = MU # Mass of vector leptoquark U1
+PowhegConfig.MGp = MU # Mass of the coloron Gp, take MGp=MU now, can also take MGp=2.5MU
+
+# General Leptoquark (LQ) Parameters
+PowhegConfig.SM = 0 # Include SM contribution
+PowhegConfig.LQ = 1 # Include basic LQ contributions
+PowhegConfig.LQ_Int = 1 # Include the interference between the SM and the LQ contributions
+PowhegConfig.bornonly = 0 # NOT Include NLO
+PowhegConfig.mass_t = 172.5 # top-quark (running) mass
+
+# Get mass_low from jo file
+matches_slice = re.search("SL([0-9]+).*", job_option_name)
+if matches_slice is None:
+    print("Cannot find betaL 33 string in job option name: {:s}.".format(job_option_name))
+    slice_low=0.0
+else:
+    slice_low= float(matches_slice.group(1))
+PowhegConfig.mass_low = slice_low # lower limit for dilepton mass
+if(slice_low==0.):
+    PowhegConfig.mass_low = -1 # upper limit for dilepton mass
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+elif(slice_low==200.):
+    PowhegConfig.mass_high = 500 # upper limit for dilepton mass
+elif(slice_low==500.):
+    PowhegConfig.mass_high = 1000 # upper limit for dilepton mass
+elif(slice_low==2000.):
+    PowhegConfig.mass_high = -1 # upper limit for dilepton mass
+else:
+    PowhegConfig.mass_high = slice_low + 1000. # upper limit for dilepton mass
+#PowhegConfig.runningscale = 1
+#PowhegConfig.new_damp = 1
+#PowhegConfig.hnew_damp = 0.5
+#PowhegConfig.hdamp = 1.0
+PowhegConfig.decay_mode = "ta+ ta-"
+
+# --------------------------------------------------------------
+# Integration settings
+# --------------------------------------------------------------
+PowhegConfig.ncall1 		= 10000
+PowhegConfig.ncall2 		= 10000
+PowhegConfig.nubound 		= 10000
+PowhegConfig.smartsig       = 0
+PowhegConfig.bornzerodamp   = 0
+PowhegConfig.hfact          = -1
+
+# --------------------------------------------------------------
+# Generate events
+# --------------------------------------------------------------
+PowhegConfig.generate()
+
+#--------------------------------------------------------------
+# Pythia8 showering with the A14 NNPDF2.3 tune, main31 routine
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+# Setting the appropriate number of final state particles for the main31 routine
+genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
-- 
GitLab


From 6d09e24e191b4ca5738784d63839fc226d39984a Mon Sep 17 00:00:00 2001
From: mcgensvc <mcgensvc@cern.ch>
Date: Tue, 23 Apr 2024 14:58:51 +0200
Subject: [PATCH 2/2] Commit: 8b392e71 pipeline succeeded. Removing
 log.generate files

---
 603xxx/603472/log.generate.short | 6 ------
 603xxx/603473/log.generate.short | 6 ------
 603xxx/603474/log.generate.short | 6 ------
 603xxx/603475/log.generate.short | 6 ------
 603xxx/603476/log.generate.short | 6 ------
 603xxx/603477/log.generate.short | 6 ------
 603xxx/603478/log.generate.short | 6 ------
 603xxx/603479/log.generate.short | 6 ------
 603xxx/603480/log.generate.short | 6 ------
 603xxx/603481/log.generate.short | 6 ------
 603xxx/603482/log.generate.short | 6 ------
 603xxx/603483/log.generate.short | 6 ------
 603xxx/603484/log.generate.short | 6 ------
 603xxx/603485/log.generate.short | 6 ------
 603xxx/603486/log.generate.short | 6 ------
 603xxx/603487/log.generate.short | 6 ------
 603xxx/603488/log.generate.short | 6 ------
 603xxx/603489/log.generate.short | 6 ------
 603xxx/603490/log.generate.short | 6 ------
 603xxx/603491/log.generate.short | 6 ------
 603xxx/603492/log.generate.short | 6 ------
 603xxx/603493/log.generate.short | 6 ------
 603xxx/603494/log.generate.short | 6 ------
 603xxx/603495/log.generate.short | 6 ------
 603xxx/603496/log.generate.short | 6 ------
 603xxx/603497/log.generate.short | 6 ------
 603xxx/603498/log.generate.short | 6 ------
 603xxx/603499/log.generate.short | 6 ------
 603xxx/603500/log.generate.short | 6 ------
 603xxx/603501/log.generate.short | 6 ------
 603xxx/603502/log.generate.short | 6 ------
 603xxx/603503/log.generate.short | 6 ------
 603xxx/603504/log.generate.short | 6 ------
 603xxx/603505/log.generate.short | 6 ------
 603xxx/603506/log.generate.short | 6 ------
 603xxx/603507/log.generate.short | 6 ------
 36 files changed, 216 deletions(-)
 delete mode 100644 603xxx/603472/log.generate.short
 delete mode 100644 603xxx/603473/log.generate.short
 delete mode 100644 603xxx/603474/log.generate.short
 delete mode 100644 603xxx/603475/log.generate.short
 delete mode 100644 603xxx/603476/log.generate.short
 delete mode 100644 603xxx/603477/log.generate.short
 delete mode 100644 603xxx/603478/log.generate.short
 delete mode 100644 603xxx/603479/log.generate.short
 delete mode 100644 603xxx/603480/log.generate.short
 delete mode 100644 603xxx/603481/log.generate.short
 delete mode 100644 603xxx/603482/log.generate.short
 delete mode 100644 603xxx/603483/log.generate.short
 delete mode 100644 603xxx/603484/log.generate.short
 delete mode 100644 603xxx/603485/log.generate.short
 delete mode 100644 603xxx/603486/log.generate.short
 delete mode 100644 603xxx/603487/log.generate.short
 delete mode 100644 603xxx/603488/log.generate.short
 delete mode 100644 603xxx/603489/log.generate.short
 delete mode 100644 603xxx/603490/log.generate.short
 delete mode 100644 603xxx/603491/log.generate.short
 delete mode 100644 603xxx/603492/log.generate.short
 delete mode 100644 603xxx/603493/log.generate.short
 delete mode 100644 603xxx/603494/log.generate.short
 delete mode 100644 603xxx/603495/log.generate.short
 delete mode 100644 603xxx/603496/log.generate.short
 delete mode 100644 603xxx/603497/log.generate.short
 delete mode 100644 603xxx/603498/log.generate.short
 delete mode 100644 603xxx/603499/log.generate.short
 delete mode 100644 603xxx/603500/log.generate.short
 delete mode 100644 603xxx/603501/log.generate.short
 delete mode 100644 603xxx/603502/log.generate.short
 delete mode 100644 603xxx/603503/log.generate.short
 delete mode 100644 603xxx/603504/log.generate.short
 delete mode 100644 603xxx/603505/log.generate.short
 delete mode 100644 603xxx/603506/log.generate.short
 delete mode 100644 603xxx/603507/log.generate.short

diff --git a/603xxx/603472/log.generate.short b/603xxx/603472/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603472/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603473/log.generate.short b/603xxx/603473/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603473/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603474/log.generate.short b/603xxx/603474/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603474/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603475/log.generate.short b/603xxx/603475/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603475/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603476/log.generate.short b/603xxx/603476/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603476/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603477/log.generate.short b/603xxx/603477/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603477/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603478/log.generate.short b/603xxx/603478/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603478/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603479/log.generate.short b/603xxx/603479/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603479/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603480/log.generate.short b/603xxx/603480/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603480/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603481/log.generate.short b/603xxx/603481/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603481/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603482/log.generate.short b/603xxx/603482/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603482/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603483/log.generate.short b/603xxx/603483/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603483/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603484/log.generate.short b/603xxx/603484/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603484/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603485/log.generate.short b/603xxx/603485/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603485/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603486/log.generate.short b/603xxx/603486/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603486/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603487/log.generate.short b/603xxx/603487/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603487/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603488/log.generate.short b/603xxx/603488/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603488/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603489/log.generate.short b/603xxx/603489/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603489/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603490/log.generate.short b/603xxx/603490/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603490/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603491/log.generate.short b/603xxx/603491/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603491/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603492/log.generate.short b/603xxx/603492/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603492/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603493/log.generate.short b/603xxx/603493/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603493/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603494/log.generate.short b/603xxx/603494/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603494/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603495/log.generate.short b/603xxx/603495/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603495/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603496/log.generate.short b/603xxx/603496/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603496/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603497/log.generate.short b/603xxx/603497/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603497/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603498/log.generate.short b/603xxx/603498/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603498/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603499/log.generate.short b/603xxx/603499/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603499/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603500/log.generate.short b/603xxx/603500/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603500/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603501/log.generate.short b/603xxx/603501/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603501/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603502/log.generate.short b/603xxx/603502/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603502/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603503/log.generate.short b/603xxx/603503/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603503/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603504/log.generate.short b/603xxx/603504/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603504/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603505/log.generate.short b/603xxx/603505/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603505/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603506/log.generate.short b/603xxx/603506/log.generate.short
deleted file mode 100644
index 875c004959..0000000000
--- a/603xxx/603506/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
diff --git a/603xxx/603507/log.generate.short b/603xxx/603507/log.generate.short
deleted file mode 100644
index 98b705f86b..0000000000
--- a/603xxx/603507/log.generate.short
+++ /dev/null
@@ -1,6 +0,0 @@
-- estimated CPU for CI job = 0.00 hrs
-- using release =  AthGeneration-23.6.16
-- ecmEnergy =  13000.0
-- randomSeed =  42
-- EVNT to EVNT =  False
-- LHEonly =  False
\ No newline at end of file
-- 
GitLab