evgenLog.info('inclusive processes will be generated')
#--------------------------------------------------------------
# add some filter here
#--------------------------------------------------------------
# need more events from MG due to filter - this needs to be set before MadGraphControl_SimplifiedModelPostInclude.py is run (it's set at 2 there)
# 3 is only sufficient for large mass splittings
evt_multiplier = 3
#filter for 3lep with low dM
if '2L4andJ20or2L7' in splitConfig[-1]:
evgenLog.info('(2leptons4 and jet20) or 2leptons7 filter')
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("DileptonFilterLow")
filtSeq += MultiElecMuTauFilter("DileptonFilterHigh")
MultiElecMuTauFilter1 = filtSeq.DileptonFilterLow
MultiElecMuTauFilter1.NLeptons = 2
MultiElecMuTauFilter1.MinPt = 4000 # low pt-cut on the lepton
MultiElecMuTauFilter1.MaxEta = 2.8 # stay away from MS 2.7 just in case
MultiElecMuTauFilter1.IncludeHadTaus = 0 # don't include hadronic taus
MultiElecMuTauFilter2 = filtSeq.DileptonFilterHigh
MultiElecMuTauFilter2.NLeptons = 2
MultiElecMuTauFilter2.MinPt = 7000 #high pt cut on the lepton
MultiElecMuTauFilter2.MaxEta = 2.8 # stay away from MS 2.7 just in case
MultiElecMuTauFilter2.IncludeHadTaus = 0 # don't include hadronic taus
include('MC15JobOptions/AntiKt4TruthWZJets.py')
# need more events from MG due to filter - this needs to be set before MadGraphControl_SimplifiedModelPostInclude.py is run
# evt_multiplier = 3 is only sufficient for large mass splittings
evt_multiplier = 18 # 10 was too small for MC15.999002.MadGraphPythia8EvtGen_A14NNPDF23LO_C1N2_Stau_217p5_182p5.py
include('MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py')
if njets > 0:
genSeq.Pythia8.Commands += [
"Merging:Process = pp>{x1+,1000024}{x1-,-1000024}{n2,1000023}",
"1000024:spinType = 1", "1000023:spinType = 1"
]
#--------------------------------------------------------------
# Lepton filter
#--------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("DileptonFilter")
filtSeq += MultiElecMuTauFilter("TauFilter")
MultiElecMuTauFilter1 = filtSeq.DileptonFilter
MultiElecMuTauFilter1.NLeptons = 2
MultiElecMuTauFilter1.MinPt = 5000.
MultiElecMuTauFilter1.MaxEta = 2.7
MultiElecMuTauFilter1.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter1.IncludeHadTaus = 1 # include hadronic taus
MultiElecMuTauFilter2 = filtSeq.TauFilter
MultiElecMuTauFilter2.NLeptons = 1
MultiElecMuTauFilter2.MinPt = 1e10
MultiElecMuTauFilter2.MaxEta = 2.7
MultiElecMuTauFilter2.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter2.IncludeHadTaus = 1 # include hadronic taus
##
## Example configuration for MultiElecMuTauFilter setting up defaults
##
if not hasattr(filtSeq, "MultiElecMuTauFilter"):
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter()
pass
# Setup defaults as example
"""
MultiElecMuTauFilter = filtSeq.MultiElecMuTauFilter
MultiElecMuTauFilter.NLeptons = 3
MultiElecMuTauFilter.MaxEta = 10.0
MultiElecMuTauFilter.MinPt = 5000.0
MultiElecMuTauFilter.MinVisPtHadTau = 10000.0
MultiElecMuTauFilter.IncludeHadTaus = True
"""