Commit ecaf98d9 authored by Cyril Becot's avatar Cyril Becot
Browse files

ATLMCPROD-8044

parent 49c496a4
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 1.5*top mass, A14 tune, at least one lepton + at least one soft muon filter, ME NNPDF30 NLO, A14 NNPDF23 LO, PS evolution 2nd order'
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'simone.amoroso@cern.ch','marco.vanadia@cern.ch']
evgenConfig.minevents = 1000
#include('PowhegControl/PowhegControl_tt_Common.py')
## Initial settings
##PowhegConfig.topdecaymode = 22222 # Inclusive
#if hasattr(PowhegConfig, "topdecaymode"):
## Use PowhegControl-00-02-XY (and earlier) syntax
#PowhegConfig.topdecaymode = 22222 # inclusive top decays
#else:
## Use PowhegControl-00-03-XY (and later) syntax
#PowhegConfig.decay_mode = "t t~ > all"
#PowhegConfig.quark_mass = 172.5
#PowhegConfig.width_t = 1.320 #from the script get_t_width.py
#PowhegConfig.hdamp = 258.75 # 1.5 * mtop
#PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0]
#PowhegConfig.mu_R = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0]
#PowhegConfig.PDF = [260000, 25200, 13165, 90900, 265000, 266000, 303400]
#PowhegConfig.PDF.extend(range(260001, 260101))
#PowhegConfig.PDF.extend(range(90901 , 90931 )) # Include the PDF4LHC error set
## Scale variations, PDF4LHC15_nlo_30
#PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15', parameter_values=[ 2.0, 1.0, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15', parameter_values=[ 0.5, 1.0, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15', parameter_values=[ 1.0, 2.0, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15', parameter_values=[ 1.0, 0.5, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15', parameter_values=[ 0.5, 0.5, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15', parameter_values=[ 2.0, 2.0, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15', parameter_values=[ 0.5, 2.0, 90900] )
#PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15', parameter_values=[ 2.0, 0.5, 90900] )
#PowhegConfig.nEvents *= 15. # compensate filter efficiency
#PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
include("MC15JobOptions/Pythia8_ShowerWeights.py")
genSeq.Pythia8.UserModes += [ 'Main31:pTHard = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 2' ]
genSeq.Pythia8.UserModes += [ 'Main31:pTdef = 2' ]
genSeq.Pythia8.UserModes += [ 'Main31:veto = 1' ]
genSeq.Pythia8.UserModes += [ 'Main31:vetoCount = 3' ]
genSeq.Pythia8.UserModes += [ 'Main31:pTemt = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:emitted = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:MPIveto = 0' ]
#rB value fitted at LEP=1.05
genSeq.Pythia8.Commands += [ 'StringZ:rFactB = 1.05' ]
genSeq.Pythia8.Commands += [ 'TimeShower:alphaSorder=2' ]
#--------------------------------------------------------------
# Event filter
#--------------------------------------------------------------
include('MC15JobOptions/TTbarWToLeptonFilter.py')
filtSeq.TTbarWToLeptonFilter.NumLeptons = -1
filtSeq.TTbarWToLeptonFilter.Ptcut = 20000.
from GeneratorFilters.GeneratorFiltersConf import TTbarWToLeptonFilter
filtSeq += TTbarWToLeptonFilter("onetWlep")
filtSeq.onetWlep.NumLeptons = 1
filtSeq.onetWlep.Ptcut = 0.
filtSeq += TTbarWToLeptonFilter("twotWlep")
filtSeq.twotWlep.NumLeptons = 2
filtSeq.twotWlep.Ptcut = 0.
from GeneratorFilters.GeneratorFiltersConf import ParentChildFilter
filtSeq += ParentChildFilter("Wmu")
filtSeq.Wmu.PDGParent = [24]
filtSeq.Wmu.PDGChild = [13]
filtSeq += ParentChildFilter("Wetau")
filtSeq.Wetau.PDGParent = [24]
filtSeq.Wetau.PDGChild = [11,15]
from GeneratorFilters.GeneratorFiltersConf import MultiMuonFilter
filtSeq += MultiMuonFilter("OneMuonFilter")
filtSeq.OneMuonFilter.Ptcut = 3250.
filtSeq.OneMuonFilter.Etacut = 2.8
filtSeq.OneMuonFilter.NMuons = 1
filtSeq += MultiMuonFilter("TwoMuonsFilter")
filtSeq.TwoMuonsFilter.Ptcut = 3250.
filtSeq.TwoMuonsFilter.Etacut = 2.8
filtSeq.TwoMuonsFilter.NMuons = 2
filtSeq += MultiMuonFilter("ThreeMuonsFilter")
filtSeq.ThreeMuonsFilter.Ptcut = 3250.
filtSeq.ThreeMuonsFilter.Etacut = 2.8
filtSeq.ThreeMuonsFilter.NMuons = 3
from GeneratorFilters.GeneratorFiltersConf import ElectronFilter
filtSeq += ElectronFilter("ElectronFilter")
filtSeq.ElectronFilter.Ptcut = 20000.
filtSeq.ElectronFilter.Etacut = 2.8
from GeneratorFilters.GeneratorFiltersConf import MuonFilter
filtSeq += MuonFilter("MuonFilter")
filtSeq.MuonFilter.Ptcut = 20000.
filtSeq.MuonFilter.Etacut = 2.8
filtSeq.Expression="TTbarWToLeptonFilter and (ElectronFilter or MuonFilter) and ((onetWlep and ((Wmu and TwoMuonsFilter) or (Wetau and OneMuonFilter))) or (twotWlep and ( ((not Wmu) and OneMuonFilter) or (Wmu and Wetau and TwoMuonsFilter) or (Wmu and (not Wetau) and ThreeMuonsFilter) )) )"
Supports Markdown
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