Commit 00a693b6 authored by Ewelina Maria Lobodzinska's avatar Ewelina Maria Lobodzinska
Browse files

Signal samples for yy-WW-semileptonic analysis (Run-2 and HL-LHC -- new JO's + EVNT + FastSim))

parent 62700cf9
Pipeline #3125786 passed with stages
in 1 minute and 49 seconds
- estimated CPU for CI job = 0.01 hrs
- using release = AthGeneration-21.6.78
- ecmEnergy = 13000.0
- LHEonly = False
import MadGraphControl.MadGraphUtils
from MadGraphControl.MadGraphUtils import *
MadGraphControl.MadGraphUtils.MADGRAPH_PDFSETTING={
'central_pdf':13300,
'pdf_variations': [13300,26400,26500,26100,26200],#None, #[26400,26500,26100,26200],
'alternative_pdfs':[13300,26500,26400,26100,26200],# None,# [26400,26500,26100,26200],
'scale_variations':[0.5,1.,2.]
}
# ----------------------------------------------
# Some global production settings
# ----------------------------------------------
# Make some excess events to allow for Pythia8 failures
nevents=2*runArgs.maxEvents if runArgs.maxEvents>0 else 5500
#mode=0
process = """
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define vll = ve vm vt ve~ vm~ vt~
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define l+ = e+ mu+ ta+
define l = e+ mu+ ta+ e- mu- ta-
define l- = e- mu- ta-
generate a a > w+ w- > l vll j j
output -f"""
process_dir = new_process(process)
#Fetch default LO run_card.dat and set parameters
settings = {'lhe_version':'3.0',
'lpp1' : '2',
'lpp2' : '2',
'fixed_fac_scale1' : 'F',
'fixed_fac_scale2' : 'F',
'ptl' : 5.,
'etal' : 4.5,
'drll': 0.2,
'drjj': 0.2,
'nevents' : int(nevents) }
modify_run_card(process_dir=process_dir,runArgs=runArgs,settings=settings)
generate(process_dir=process_dir,runArgs=runArgs,grid_pack=False)
arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=3,saveProcDir=True)
evgenConfig.generators = ["MadGraph"]
runName='PROC_sm_0/run_01'
############################
# Shower JOs will go here
#runName=
include("Pythia8_i/Py8_NNPDF23_NNLO_as118_QED_EE_Common.py")
include("Pythia8_i/Pythia8_ShowerWeights.py")
include("Pythia8_i//Pythia8_MadGraph.py")
include('GeneratorFilters/MultiLeptonFilter.py')
MultiLeptonFilter = filtSeq.MultiLeptonFilter
MultiLeptonFilter.Ptcut = 3500.
MultiLeptonFilter.Etacut = 4.5
MultiLeptonFilter.NLeptons = 1
#--------------------------------------------------------------
# Evgen
#--------------------------------------------------------------
evgenConfig.description = 'MadGraphPy8EG_yyWW'
evgenConfig.keywords+=["SM", "Z"]
evgenConfig.contact = ["Kristin Lohwasser <kristin@cern.ch"]
# JobOption fragment for Herwig >=7.0.4
# Common job option for gamma gamma processes in Herwig7
# Contact: Oldrich Kepka
assert hasattr(genSeq, "Herwig7")
Herwig7Config.add_commands(Herwig7Config.pdf_gammagamma_cmds())
cmds = """
set /Herwig/Partons/PPExtractor:FirstPDF /Herwig/Partons/BudnevPDF
set /Herwig/Partons/PPExtractor:SecondPDF /Herwig/Partons/BudnevPDF
set /Herwig/Partons/BudnevPDF:Q2Max 4.
create ThePEG::O1AlphaS /Herwig/AlphaQCD_O1 O1AlphaS.so
set /Herwig/Generators/EventGenerator:StandardModelParameters:QCD/RunningAlphaS /Herwig/AlphaQCD_O1
# Cuts
cd /Herwig/Cuts
set Cuts:ScaleMin 0.0*GeV2
set Cuts:X1Min 0
set Cuts:X2Min 0
set Cuts:X1Max 1.
set Cuts:X2Max 1.
erase Cuts:MultiCuts 0
"""
Herwig7Config.add_commands(cmds)
del cmds
# To avoid warning from displaced vertices, bugfix needed in herwig++
testSeq.TestHepMC.MaxTransVtxDisp = 1000000
testSeq.TestHepMC.MaxVtxDisp = 1000000000
- estimated CPU for CI job = 0.00 hrs
- using release = AthGeneration-21.6.78
- ecmEnergy = 13000.0
- LHEonly = False
evgenConfig.description = "gammagamma->WW with Budnev parameterization, leptonic decays, central lepton filter pt>15 GeV, "
evgenConfig.keywords = ["2lepton", "diboson", "SM", "WW", "electroweak", "exclusive"]
evgenConfig.contact = ["Kristin Lohwasser <kristin@cern.ch"]
evgenConfig.generators = ["Herwig7", "EvtGen"]
evgenConfig.tune = "H7.2-Default"
## Initialise Herwig7 for run with built-in/old-style matrix elements
include("Herwig7_i/Herwig7_BuiltinME.py")
Herwig7Config.pdf_gammagamma_cmds()
include("./Herwig7_QED_EvtGen_Common.py")
commands = """
#do /Herwig/Particles/W-:PrintDecayModes
#do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+; W+->nu_tau,tau+;
#do /Herwig/Particles/W-:SelectDecayModes W-->nu_ebar,e-; W-->nu_mubar,mu-; W-->nu_taubar,tau-;
# Selected the hard process
cd /Herwig/MatrixElements
insert SubProcess:MatrixElements 0 MEgg2WW
"""
Herwig7Config.add_commands(commands)
del commands
Herwig7Config.run()
include('GeneratorFilters/MultiLeptonFilter.py')
MultiLeptonFilter = filtSeq.MultiLeptonFilter
MultiLeptonFilter.Ptcut = 3500.
MultiLeptonFilter.Etacut = 4.5
MultiLeptonFilter.NLeptons = 1
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment