From db522996c8b761eb42285fcadec008c4e4f6acaa Mon Sep 17 00:00:00 2001
From: Yang Liu <yang.l@cern.ch>
Date: Fri, 20 Dec 2024 15:09:06 +0800
Subject: [PATCH 1/2] Generation of aMC@NLO samples for the ttZ with dilep
decay (MC20/23)
---
525xxx/525727/log.generate.short | 7 +
...8EG_NNPDF30NLO_A14N23LO_ttZnunu_ttdilep.py | 1 +
525xxx/525727/ttZ_NLO_ttdilep.py | 226 ++++++++++++++++++
3 files changed, 234 insertions(+)
create mode 100644 525xxx/525727/log.generate.short
create mode 100644 525xxx/525727/mc.aMCPy8EG_NNPDF30NLO_A14N23LO_ttZnunu_ttdilep.py
create mode 100644 525xxx/525727/ttZ_NLO_ttdilep.py
diff --git a/525xxx/525727/log.generate.short b/525xxx/525727/log.generate.short
new file mode 100644
index 0000000000..72bc845fe9
--- /dev/null
+++ b/525xxx/525727/log.generate.short
@@ -0,0 +1,7 @@
+- estimated CPU for CI job = 0.01 hrs
+- using release = AthGeneration-23.6.23
+- ecmEnergy = 13600.0
+- randomSeed = 1
+- EVNT to EVNT = False
+- LHEonly = False
+- platform = Linux-5.14.0-503.11.1.el9_5.x86_64-x86_64-with-glibc2.34
\ No newline at end of file
diff --git a/525xxx/525727/mc.aMCPy8EG_NNPDF30NLO_A14N23LO_ttZnunu_ttdilep.py b/525xxx/525727/mc.aMCPy8EG_NNPDF30NLO_A14N23LO_ttZnunu_ttdilep.py
new file mode 100644
index 0000000000..19a452ed04
--- /dev/null
+++ b/525xxx/525727/mc.aMCPy8EG_NNPDF30NLO_A14N23LO_ttZnunu_ttdilep.py
@@ -0,0 +1 @@
+include('../100001/ttZ_NLO_ttdilep.py')
\ No newline at end of file
diff --git a/525xxx/525727/ttZ_NLO_ttdilep.py b/525xxx/525727/ttZ_NLO_ttdilep.py
new file mode 100644
index 0000000000..5878ecb9bd
--- /dev/null
+++ b/525xxx/525727/ttZ_NLO_ttdilep.py
@@ -0,0 +1,226 @@
+# Original JO 522023/ttZ_NLO_run3.py
+# In this CO file, the dilep decay of the two tops have been added. No other changes have been made wrt the original JO.
+
+from MadGraphControl.MadGraphUtils import *
+from MadGraphControl.MadGraphUtilsHelpers import get_physics_short
+from MadGraphControl.MadGraphParamHelpers import set_top_params
+import MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
+import os
+
+# General settings
+nevents = int(1.1*runArgs.maxEvents) # safety factor
+mllcut = 5 # mininum m(ll) mass (to cut off DY)
+zdecay = 'decay z > all all' # placeholder
+ptgmincut = 20 # photons with pT>20 GeV
+lhe_version = 3 # LHE version to save reweighted scales+PDF
+gridpack_dir = 'madevent/' # to generate from gridpack
+gridpack_mode = False # we have gridpacks for everything so far
+
+beamEnergy = -999 # safety check
+if hasattr(runArgs,'ecmEnergy'):
+ beamEnergy = runArgs.ecmEnergy / 2.
+else:
+ raise RuntimeError("No center of mass energy found.")
+
+# Parse the name of the JO to extract information about the process
+jo_name = get_physics_short()
+ttee = 'ttee' in jo_name
+ttmumu = 'ttmumu' in jo_name
+tttautau = 'tttautau' in jo_name
+ttZnunu = 'ttZnunu' in jo_name
+ttZqq = 'ttZqq' in jo_name
+lowmass = 'lowmass' in jo_name
+a14Var1Up = 'A14Var1Up' in jo_name
+a14Var1Down = 'A14Var1Down' in jo_name
+a14Var2Up = 'A14Var2Up' in jo_name
+a14Var2Down = 'A14Var2Down' in jo_name
+a14Var3aUp = 'A14Var3aUp' in jo_name
+a14Var3aDown = 'A14Var3aDown' in jo_name
+a14Var3bUp = 'A14Var3bUp' in jo_name
+a14Var3bDown = 'A14Var3bDown' in jo_name
+a14Var3cUp = 'A14Var3cUp' in jo_name
+a14Var3cDown = 'A14Var3cDown' in jo_name
+pythia = 'Py8' in jo_name
+herwig = 'H7' in jo_name
+
+# Set MadGraph generate commands
+if ttee:
+ mgproc = '''generate p p > t t~ e+ e- [QCD]'''
+ name = 'ttee_NLO_ttdilep'
+elif ttmumu:
+ mgproc = '''generate p p > t t~ mu+ mu- [QCD]'''
+ name = 'ttmumu_NLO_ttdilep'
+elif tttautau:
+ mgproc = '''generate p p > t t~ ta+ ta- $$ h [QCD]'''
+ name = 'tttautau_NLO_ttdilep'
+elif ttZnunu:
+ mgproc = '''generate p p > t t~ z [QCD]'''
+ name = 'ttZnunu_NLO_ttdilep'
+ zdecay = 'decay z > vl vl~'
+elif ttZqq:
+ mgproc = '''generate p p > t t~ z [QCD]'''
+ name = 'ttZqq_NLO_ttdilep'
+ zdecay = 'decay z > q q~'
+else:
+ raise RuntimeError("If your JO doesn't concern ttee, ttmumu, tttautau, ttZnunu or ttZqq, you've imported the wrong control file, buddy!")
+
+# Set the parton shower
+if pythia:
+ psname = 'PYTHIA8'
+elif herwig:
+ psname = 'HERWIGPP'
+else:
+ raise RuntimeError("Can't tell whether the parton shower is supposed to be Pythia8 or Herwig7, stopping this madness now!")
+
+# Write the proc_card
+process = '''
+import model loop_sm-no_b_mass
+define p = g u c d s b u~ c~ d~ s~ b~
+define j = g u c d s b u~ c~ d~ s~ b~
+define w = w+ w-
+define l+ = e+ mu+ ta+
+define l- = e- mu- ta-
+'''+mgproc+'''
+output -f
+'''
+process_dir = new_process(process)
+
+# Low mass cuts: 1-5 GeV
+if lowmass:
+ mllcut = 1
+ get_cuts_file = subprocess.Popen(['get_files','-data','ttZ_lowmass_cuts.f'])
+ get_cuts_file.wait()
+ import os, shutil
+ if not os.path.exists('ttZ_lowmass_cuts.f'):
+ raise RuntimeError("Cannot find ttZ_lowmass_cuts.f")
+ shutil.copy('ttZ_lowmass_cuts.f', os.path.join(process_dir, 'SubProcesses/cuts.f'))
+
+# Write the run_card
+settings = {
+ 'nevents': nevents,
+ 'maxjetflavor': 5,
+ 'parton_shower': psname,
+ 'ptgmin': ptgmincut,
+ 'mll_sf': mllcut,
+}
+
+modify_run_card(process_dir=process_dir,
+ runArgs=runArgs,
+ settings=settings)
+
+# Modify the param card (important! this corrects the BR(W->had) issue in previous ttZ samples
+set_top_params(process_dir,mTop=172.5,FourFS=False)
+
+# Write the MadSpin card
+madspin_card = process_dir+'/Cards/madspin_card.dat'
+if os.access(madspin_card,os.R_OK):
+ os.unlink(madspin_card)
+mscard = open(madspin_card,'w')
+mscard.write("""#************************************************************
+#* MadSpin *
+#* *
+#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
+#* *
+#* Part of the MadGraph5_aMC@NLO Framework: *
+#* The MadGraph5_aMC@NLO Development Team - Find us at *
+#* https://server06.fynu.ucl.ac.be/projects/madgraph *
+#* *
+#************************************************************
+#Some options (uncomment to apply)
+#
+set seed %i
+set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
+set BW_cut 15 # cut on how far the particle can be off-shell
+set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
+# specify the decay for the final state particles
+define q = u d s c b
+define q~ = u~ d~ s~ c~ b~
+decay t > w+ b, w+ > l+ vl
+decay t~ > w- b~, w- > l- vl~
+decay w+ > l+ vl
+decay w- > l- vl~
+%s
+# running the actual code
+launch"""%(runArgs.randomSeed,zdecay))
+mscard.close()
+
+# Print and generate
+print_cards()
+generate(process_dir=process_dir,
+ grid_pack=gridpack_mode,
+ gridpack_compile=False, # set to True to generate the gridpack, to False to run from it
+ required_accuracy=0.001, # gridpack accuracy, ignored if running OTF
+ runArgs=runArgs)
+outputDS = arrange_output(process_dir=process_dir,
+ runArgs=runArgs,
+ saveProcDir=False,
+ lhe_version=lhe_version)
+
+# Metadata
+evgenConfig.description = 'aMcAtNlo_'+str(name)
+evgenConfig.keywords += ['SM','ttZ']
+evgenConfig.contact = ['baptiste.ravina@cern.ch']
+evgenConfig.generators = ['aMcAtNlo','EvtGen']
+if pythia:
+ evgenConfig.generators += ['Pythia8']
+if herwig:
+ evgenConfig.generators += ['Herwig7']
+if a14Var1Up:
+ evgenConfig.description += '_A14Var1Up'
+if a14Var1Down:
+ evgenConfig.description += '_A14Var1Down'
+if a14Var2Up:
+ evgenConfig.description += '_A14Var2Up'
+if a14Var2Down:
+ evgenConfig.description += '_A14Var2Down'
+if a14Var3aUp:
+ evgenConfig.description += '_A14Var3aUp'
+if a14Var3aDown:
+ evgenConfig.description += '_A14Var3aDown'
+if a14Var3bUp:
+ evgenConfig.description += '_A14Var3bUp'
+if a14Var3bDown:
+ evgenConfig.description += '_A14Var3bDown'
+if a14Var3cUp:
+ evgenConfig.description += '_A14Var3cUp'
+if a14Var3cDown:
+ evgenConfig.description += '_A14Var3cDown'
+
+# go back to single-thread
+check_reset_proc_number(opts)
+# shower
+if pythia:
+ if a14Var1Up:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var1Up_EvtGen_Common.py")
+ elif a14Var1Down:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var1Down_EvtGen_Common.py")
+ elif a14Var2Up:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var2Up_EvtGen_Common.py")
+ elif a14Var2Down:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var2Down_EvtGen_Common.py")
+ elif a14Var3aUp:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var3aUp_EvtGen_Common.py")
+ elif a14Var3aDown:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var3aDown_EvtGen_Common.py")
+ elif a14Var3bUp:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var3bUp_EvtGen_Common.py")
+ elif a14Var3bDown:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var3bDown_EvtGen_Common.py")
+ elif a14Var3cUp:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var3cUp_EvtGen_Common.py")
+ elif a14Var3cDown:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_Var3cDown_EvtGen_Common.py")
+ else:
+ include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+ include("Pythia8_i/Pythia8_aMcAtNlo_decayMEC.py")
+elif herwig:
+ evgenConfig.tune = "H7.2-Default"
+ include("Herwig7_i/Herwig7_LHEF.py")
+ include("Herwig7_i/Herwig7_EvtGen.py")
+ runArgs.inputGeneratorFile = outputDS+'.events'
+ Herwig7Config.me_pdf_commands(order="NLO",
+ name="NNPDF30_nlo_as_0118")
+ Herwig7Config.tune_commands()
+ Herwig7Config.lhef_mg5amc_commands(lhe_filename=runArgs.inputGeneratorFile,
+ me_pdf_order="NLO")
+ Herwig7Config.run()
--
GitLab
From 268dbc6e92d9cf608a369b1bbcdbabda87e6d54b Mon Sep 17 00:00:00 2001
From: mcgensvc <mcgensvc@cern.ch>
Date: Fri, 20 Dec 2024 09:12:23 +0100
Subject: [PATCH 2/2] Commit: 3b002c17 pipeline succeeded. Removing
log.generate files
---
525xxx/525727/log.generate.short | 7 -------
1 file changed, 7 deletions(-)
delete mode 100644 525xxx/525727/log.generate.short
diff --git a/525xxx/525727/log.generate.short b/525xxx/525727/log.generate.short
deleted file mode 100644
index 72bc845fe9..0000000000
--- a/525xxx/525727/log.generate.short
+++ /dev/null
@@ -1,7 +0,0 @@
-- estimated CPU for CI job = 0.01 hrs
-- using release = AthGeneration-23.6.23
-- ecmEnergy = 13600.0
-- randomSeed = 1
-- EVNT to EVNT = False
-- LHEonly = False
-- platform = Linux-5.14.0-503.11.1.el9_5.x86_64-x86_64-with-glibc2.34
\ No newline at end of file
--
GitLab