Commit 14f66b8b authored by Ewelina Maria Lobodzinska's avatar Ewelina Maria Lobodzinska Committed by Christian Gutschow
Browse files

Add additional JOs for filters validation and harmonisation

parent fb3ecfd1
Pipeline #2960631 passed with stages
in 1 minute and 5 seconds
......@@ -2,7 +2,7 @@ from MadGraphControl.MadGraphUtils import *
nevents = runArgs.maxEvents*40. if runArgs.maxEvents>0 else 40.*evgenConfig.nEventsPerJob
gridpack_mode=False
process_def = """
import model sm
define p = g u c d s u~ c~ d~ s~
......@@ -27,7 +27,7 @@ generate(process_dir=process_dir,grid_pack=gridpack_mode,runArgs=runArgs)
arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=3,saveProcDir=False)
#### Shower
evgenConfig.description = 'MadGraph_Zee'
evgenConfig.description = 'MadGraph_Zmumu'
evgenConfig.keywords+=['Z','jets']
evgenConfig.contact = [ "zach.marshall@cern.ch","giancarlo.panizzo@cern.ch" ]
evgenConfig.generators += ["MadGraph"]
......
# based on 800889
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("Pythia8_i/Pythia8_ShowerWeights.py")
evgenConfig.description = "Dijet truth jet slice JZ2, with the A14 NNPDF23 LO tune, with c hadron filter "
evgenConfig.keywords = ["QCD", "jets", "SM"]
evgenConfig.contact = ["simone.amoroso@cern.ch","alex.emerman@cern.ch"]
evgenConfig.nEventsPerJob = 1000
genSeq.Pythia8.Commands += ["HardQCD:all = on",
"PhaseSpace:pTHatMin = 15."]
include("GeneratorFilters/FindJets.py")
CreateJets(prefiltSeq, 0.4)
if not hasattr( filtSeq, "MultiCjetFilter" ):
from GeneratorFilters.GeneratorFiltersConf import MultiCjetFilter
filtSeq += MultiCjetFilter()
print ("MultiCjetFilter added ...")
pass
filtSeq.MultiCjetFilter.NCJetsMin = 1
AddJetsFilter(filtSeq,runArgs.ecmEnergy, 0.4)
#include("GeneratorFilters/JetFilter_JZX.py")
#JZSlice(2,filtSeq)
#--------------------------------------------------------------
# Powheg Z setup starting from ATLAS defaults
#--------------------------------------------------------------
# based on 361106 (in 21.6 600701)
include('PowhegControl/PowhegControl_Z_Common.py')
PowhegConfig.decay_mode = "z > e+ e-"
# Configure Powheg setup
PowhegConfig.ptsqmin = 4.0 # needed for AZNLO tune
PowhegConfig.running_width = 1
PowhegConfig.nEvents *= 3
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering with AZNLO_CTEQ6L1 and Photos
#--------------------------------------------------------------
include('Pythia8_i/Pythia8_AZNLO_CTEQ6L1_EvtGen_Common.py')
include('Pythia8_i/Pythia8_Photospp.py')
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 Z->ee production without lepton filter and AZNLO CT10 tune'
evgenConfig.contact = ["oldrich.kepka@cern.ch"]
evgenConfig.keywords = [ 'NLO', 'SM', 'electroweak', 'Z', 'drellYan', '2electron' ]
evgenConfig.nEventsPerJob = 5000
#--------------------------------------------------------------
# FILTERS
#--------------------------------------------------------------
#include('GeneratorFilters/MultiLeptonFilter.py')
## Default cut params
#filtSeq.MultiLeptonFilter.Ptcut = 3500.
#filtSeq.MultiLeptonFilter.Etacut = 2.7
#filtSeq.MultiLeptonFilter.NLeptons = 2
include('GeneratorFilters/ChargedTrackFilter.py')
filtSeq.ChargedTracksFilter.NTracks = 21
#filtSeq.ChargedTracksFilter.NTracksMax = -1
filtSeq.ChargedTracksFilter.Ptcut = 500
filtSeq.ChargedTracksFilter.Etacut= 2.5
import MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from MadGraphControl.MadGraphUtils import *
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor=30
nevents=5000*safefactor
if runArgs.maxEvents > 0:
nevents=runArgs.maxEvents*safefactor
mode=0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters = {}
parameters['NEW'] = {'CV': '0.5000000', # CV
'C2V': '1.0000000', # C2V
'C3': '-5.0000000'} # C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters['MASS']={'25':'1.250000e+02'} #MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'scale':'125',
'dsqrt_q2fact1':'125',
'dsqrt_q2fact2':'125',
'nevents':int(nevents)}
#---------------------------------------------------------------------------------------------------
# Generating non-resonant VBF-Only HH process with MadGraph
# Parameters are set above
#---------------------------------------------------------------------------------------------------
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~
import model /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j $$ z w+ w- / a j QED=4
output -f"""
process_dir = new_process(process)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new run_card.dat from an existing one
# Using the values given in "extras" above for the selected parameters when setting up the run_card
# If not set in "extras", default values are used
#---------------------------------------------------------------------------------------------------
modify_run_card(process_dir=process_dir,runArgs=runArgs,settings=extras)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for CV, C2V, C3
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
modify_param_card(process_dir=process_dir,params=parameters)
#---------------------------------------------------------------------------------------------------
# Printing cards
#---------------------------------------------------------------------------------------------------
print_cards()
#---------------------------------------------------------------------------------------------------
# Generate events
#---------------------------------------------------------------------------------------------------
generate(process_dir=process_dir,runArgs=runArgs)
#---------------------------------------------------------------------------------------------------
# Move output files into the appropriate place, with the appropriate name
#---------------------------------------------------------------------------------------------------
arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=2,saveProcDir=True)
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8"]
evgenConfig.description = "Non-resonant LO di-Higgs production through vector-boson-fusion (VBF) which decays to bbtautau (hadhad)."
evgenConfig.keywords = ["hh","SM", "SMHiggs", "nonResonant", "VBF", "bottom"]
evgenConfig.contact = ['Tulin Varol Mete <Tulin.Varol@cern.ch>']
evgenConfig.nEventsPerJob = 10000
#evgenConfig.tune = "MMHT2014"
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("Pythia8_i/Pythia8_MadGraph.py")
genSeq.Pythia8.Commands += [ "25:oneChannel = on 0.5 100 5 -5", # bb decay
"25:addChannel = on 0.5 100 15 -15" ] # tautau decay
#--------------------------------------------------------------
# Dipole option Pythia8
#--------------------------------------------------------------
genSeq.Pythia8.Commands += [ "SpaceShower:dipoleRecoil = on" ]
#---------------------------------------------------------------------------------------------------
# Filter for bbtautau
#---------------------------------------------------------------------------------------------------
#from GeneratorFilters.GeneratorFiltersConf import ParentChildFilter
#filtSeq += ParentChildFilter("HbbFilter", PDGParent = [25], PDGChild = [5])
#filtSeq += ParentChildFilter("HTauTauFilter", PDGParent = [25], PDGChild = [15])
from GeneratorFilters.GeneratorFiltersConf import XtoVVDecayFilterExtended
filtSeq += XtoVVDecayFilterExtended("TauTauHadHadFilter")
filtSeq.TauTauHadHadFilter.PDGGrandParent = 25
filtSeq.TauTauHadHadFilter.PDGParent = 15
filtSeq.TauTauHadHadFilter.StatusParent = 2
filtSeq.TauTauHadHadFilter.PDGChild1 = [211,213,215,311,321,323,10232,10323,20213,20232,20323,30213,100213,100323,1000213]
filtSeq.TauTauHadHadFilter.PDGChild2 = [211,213,215,311,321,323,10232,10323,20213,20232,20323,30213,100213,100323,1000213]
#---------------------------------------------------------------------------------------------------
# Filter for 2 leptons (tau(had) tau(had)) with pt cuts
#---------------------------------------------------------------------------------------------------
#from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
#filtSeq += MultiElecMuTauFilter("TauPtFilter")
#filtSeq.TauPtFilter.IncludeHadTaus = True
#filtSeq.TauPtFilter.NLeptons = 2
#filtSeq.TauPtFilter.MinPt = 13000.
#filtSeq.TauPtFilter.MinVisPtHadTau = 15000.
#filtSeq.TauPtFilter.MaxEta = 3.
#---------------------------------------------------------------------------------------------------
# Leading tau filter
#---------------------------------------------------------------------------------------------------
#filtSeq += MultiElecMuTauFilter("LeadTauPtFilter")
#filtSeq.LeadTauPtFilter.IncludeHadTaus = True
#filtSeq.LeadTauPtFilter.NLeptons = 1
#filtSeq.LeadTauPtFilter.MinPt = 13000.
#filtSeq.LeadTauPtFilter.MinVisPtHadTau = 35000.
#filtSeq.LeadTauPtFilter.MaxEta = 3.
#filtSeq.Expression = "HbbFilter and HTauTauFilter and TauTauHadHadFilter and TauPtFilter and LeadTauPtFilter"
import MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from MadGraphControl.MadGraphUtils import *
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor=30
nevents=5000*safefactor
if runArgs.maxEvents > 0:
nevents=runArgs.maxEvents*safefactor
mode=0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters = {}
parameters['NEW'] = {'CV': '0.5000000', # CV
'C2V': '1.0000000', # C2V
'C3': '-5.0000000'} # C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters['MASS']={'25':'1.250000e+02'} #MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'scale':'125',
'dsqrt_q2fact1':'125',
'dsqrt_q2fact2':'125',
'nevents':int(nevents)}
#---------------------------------------------------------------------------------------------------
# Generating non-resonant VBF-Only HH process with MadGraph
# Parameters are set above
#---------------------------------------------------------------------------------------------------
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~
import model /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j $$ z w+ w- / a j QED=4
output -f"""
process_dir = new_process(process)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new run_card.dat from an existing one
# Using the values given in "extras" above for the selected parameters when setting up the run_card
# If not set in "extras", default values are used
#---------------------------------------------------------------------------------------------------
modify_run_card(process_dir=process_dir,runArgs=runArgs,settings=extras)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for CV, C2V, C3
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
modify_param_card(process_dir=process_dir,params=parameters)
#---------------------------------------------------------------------------------------------------
# Printing cards
#---------------------------------------------------------------------------------------------------
print_cards()
#---------------------------------------------------------------------------------------------------
# Generate events
#---------------------------------------------------------------------------------------------------
generate(process_dir=process_dir,runArgs=runArgs)
#---------------------------------------------------------------------------------------------------
# Move output files into the appropriate place, with the appropriate name
#---------------------------------------------------------------------------------------------------
arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=2,saveProcDir=True)
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8"]
evgenConfig.description = "Non-resonant LO di-Higgs production through vector-boson-fusion (VBF) which decays to bbtautau (hadhad)."
evgenConfig.keywords = ["hh","SM", "SMHiggs", "nonResonant", "VBF", "bottom"]
evgenConfig.contact = ['Tulin Varol Mete <Tulin.Varol@cern.ch>']
evgenConfig.nEventsPerJob = 10000
#evgenConfig.tune = "MMHT2014"
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("Pythia8_i/Pythia8_MadGraph.py")
genSeq.Pythia8.Commands += [ "25:oneChannel = on 0.5 100 5 -5", # bb decay
"25:addChannel = on 0.5 100 15 -15" ] # tautau decay
#--------------------------------------------------------------
# Dipole option Pythia8
#--------------------------------------------------------------
genSeq.Pythia8.Commands += [ "SpaceShower:dipoleRecoil = on" ]
#---------------------------------------------------------------------------------------------------
# Filter for bbtautau
#---------------------------------------------------------------------------------------------------
#from GeneratorFilters.GeneratorFiltersConf import ParentChildFilter
#filtSeq += ParentChildFilter("HbbFilter", PDGParent = [25], PDGChild = [5])
#filtSeq += ParentChildFilter("HTauTauFilter", PDGParent = [25], PDGChild = [15])
from GeneratorFilters.GeneratorFiltersConf import XtoVVDecayFilterExtended
filtSeq += XtoVVDecayFilterExtended("TauTauHadHadFilter")
filtSeq.TauTauHadHadFilter.PDGGrandParent = 25
filtSeq.TauTauHadHadFilter.PDGParent = 15
filtSeq.TauTauHadHadFilter.StatusParent = 2
filtSeq.TauTauHadHadFilter.PDGChild1 = [211,213,215,311,321,323,10232,10323,20213,20232,20323,30213,100213,100323,1000213]
filtSeq.TauTauHadHadFilter.PDGChild2 = [211,213,215,311,321,323,10232,10323,20213,20232,20323,30213,100213,100323,1000213]
#---------------------------------------------------------------------------------------------------
# Filter for 2 leptons (tau(had) tau(had)) with pt cuts
#---------------------------------------------------------------------------------------------------
#from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
#filtSeq += MultiElecMuTauFilter("TauPtFilter")
#filtSeq.TauPtFilter.IncludeHadTaus = True
#filtSeq.TauPtFilter.NLeptons = 2
#filtSeq.TauPtFilter.MinPt = 13000.
#filtSeq.TauPtFilter.MinVisPtHadTau = 15000.
#filtSeq.TauPtFilter.MaxEta = 3.
#---------------------------------------------------------------------------------------------------
# Leading tau filter
#---------------------------------------------------------------------------------------------------
#filtSeq += MultiElecMuTauFilter("LeadTauPtFilter")
#filtSeq.LeadTauPtFilter.IncludeHadTaus = True
#filtSeq.LeadTauPtFilter.NLeptons = 1
#filtSeq.LeadTauPtFilter.MinPt = 13000.
#filtSeq.LeadTauPtFilter.MinVisPtHadTau = 35000.
#filtSeq.LeadTauPtFilter.MaxEta = 3.
#filtSeq.Expression = "HbbFilter and HTauTauFilter and TauTauHadHadFilter and TauPtFilter and LeadTauPtFilter"
include("Sherpa_i/Base_Fragment.py")
include("Sherpa_i/NNPDF30NNLO.py")
evgenConfig.description = "opposite sign llvvjj-EW-noHiggs with high mjj filter"
evgenConfig.keywords = ["SM", "diboson", "2lepton", "VBS"]
evgenConfig.contact = ["atlas-generators-sherpa@cern.ch", "chris.g@cern.ch"]
evgenConfig.nEventsPerJob = 1000
genSeq.Sherpa_i.RunCard="""
(run){
EW_SCHEME=3;
ACTIVE[25]=0;
MASSIVE[5]=1;
MASSIVE[11]=1;
MASSIVE[13]=1;
MASSIVE[15]=1;
PARTICLE_CONTAINER 991[m:-1] leptons 11 -11 13 -13 15 -15;
%scales, tags for scale variations
FSCF:=1.0; RSCF:=1.0; QSCF:=1.0;
SCALES=STRICT_METS{FSCF*MU_F2}{RSCF*MU_R2}{QSCF*MU_Q2};
CORE_SCALE=VAR{Abs2(p[2]+p[3]+p[4]+p[5])};
# simplified setup as long as only 2->6 taken into account:
#SCALES=VAR{FSCF*Abs2(p[2]+p[3]+p[4]+p[5])}{RSCF*Abs2(p[2]+p[3]+p[4]+p[5])}{QSCF*Abs2(p[2]+p[3]+p[4]+p[5])};
%tags for process setup
NJET:=0; QCUT:=20.;
EXCLUSIVE_CLUSTER_MODE=1;
SOFT_SPIN_CORRELATIONS=1
% improve colour-flow treatment
CSS_CSMODE=1
% improve integration performance
PSI_ITMIN=25000;
CDXS_VSOPT=5;
INTEGRATION_ERROR 0.05;
}(run)
(processes){
# V V -> e- antinu_e e+ nu_e
Process 93 93 -> 11 -12 -11 12 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> e- antinu_e mu+ nu_mu
Process 93 93 -> 11 -12 -13 14 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> e- antinu_e tau+ nu_tau
Process 93 93 -> 11 -12 -15 16 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> mu- antinu_mu e+ nu_e
Process 93 93 -> 13 -14 -11 12 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> mu- antinu_mu mu+ nu_mu
Process 93 93 -> 13 -14 -13 14 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> mu- antinu_mu tau+ nu_tau
Process 93 93 -> 13 -14 -15 16 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> tau- antinu_tau e+ nu_e
Process 93 93 -> 15 -16 -11 12 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> tau- antinu_tau mu+ nu_mu
Process 93 93 -> 15 -16 -13 14 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> tau- antinu_tau tau+ nu_tau
Process 93 93 -> 15 -16 -15 16 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> e- e+ nu_mu antinu_mu
Process 93 93 -> 11 -11 14 -14 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> e- e+ nu_tau antinu_tau
Process 93 93 -> 11 -11 16 -16 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> mu- mu+ nu_e antinu_e
Process 93 93 -> 13 -13 12 -12 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> mu- mu+ nu_tau antinu_tau
Process 93 93 -> 13 -13 16 -16 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> tau- tau+ nu_e antinu_e
Process 93 93 -> 15 -15 12 -12 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
# V V -> tau- tau+ nu_mu antinu_mu
Process 93 93 -> 15 -15 14 -14 93 93 93{NJET};
Order (*,6); CKKW sqr(QCUT/E_CMS);
End process;
}(processes)
(selector){
Mass 11 -11 0.1 E_CMS
"PT" 991 5.0,E_CMS:5.0,E_CMS [PT_UP]
"Eta" 991 -5.0,5.0:-5.0,5.0 [PT_UP]
NJetFinder 2 15. 0. 0.4 -1;
}(selector)
"""
genSeq.Sherpa_i.Parameters += [ "EW_SCHEME=3", "GF=1.166397e-5" ]
genSeq.Sherpa_i.NCores = 16
#--------------------------------------------------------------
# VBF Mjj Interval filter
#--------------------------------------------------------------
include("GeneratorFilters/FindJets.py")
CreateJets(prefiltSeq, 0.4)
# Set up VBF filters
if not hasattr(filtSeq, "VBFForwardJetsFilter" ):
from GeneratorFilters.GeneratorFiltersConf import VBFForwardJetsFilter
VBFfilter = VBFForwardJetsFilter("VBFfilter")
filtSeq += VBFfilter
filtSeq.VBFfilter.JetMinPt=20.*GeV
filtSeq.VBFfilter.JetMaxEta=5.0
filtSeq.VBFfilter.NJets=2
filtSeq.VBFfilter.Jet1MinPt=30.*GeV
filtSeq.VBFfilter.Jet1MaxEta=5.0
filtSeq.VBFfilter.Jet2MinPt=20.*GeV
filtSeq.VBFfilter.Jet2MaxEta=5.0
filtSeq.VBFfilter.UseOppositeSignEtaJet1Jet2=False
filtSeq.VBFfilter.TruthJetContainer="AntiKt4TruthJets"
filtSeq.VBFfilter.LGMinPt=15.*GeV
filtSeq.VBFfilter.LGMaxEta=2.5
filtSeq.VBFfilter.DeltaRJLG=0.05
filtSeq.VBFfilter.RatioPtJLG=0.3
filtSeq.VBFfilter.UseLeadingJJ=True
filtSeq.VBFfilter.MassJJ = 300.*GeV
filtSeq.VBFfilter.DeltaEtaJJ = 3.0
/eos/user/c/cgutscho/mc/mc_13TeV.Sh_2211_llvvjj_os_EW6_noHiggs_highMjj.GRID.tar.gz
\ No newline at end of file
from MadGraphControl.MadGraphUtils import *
nevents = runArgs.maxEvents*40. if runArgs.maxEvents>0 else 40.*evgenConfig.nEventsPerJob
gridpack_mode=False
process_def = """
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~
generate p p > mu+ mu-
output -f"""
#Fetch default LO run_card.dat and set parameters
settings = { 'lhe_version':'3.0',
'cut_decays':'F',
'pdlabel':"'nn23lo1'",
'mmll':40,
'use_syst':"False",
'nevents':int(nevents)}
process_dir = new_process(process_def)
modify_run_card(process_dir=process_dir,runArgs=runArgs,settings=settings)
generate(process_dir=process_dir,grid_pack=gridpack_mode,runArgs=runArgs)
# saveProcDir=True for testing only
arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=3,saveProcDir=False)
#### Shower