diff --git a/950xxx/950787/mc.PhPy8EG_VINCIA_tt_hdamp258p75_nonallhad_valid.py b/950xxx/950787/mc.PhPy8EG_VINCIA_tt_hdamp258p75_nonallhad_valid.py new file mode 100644 index 0000000000000000000000000000000000000000..87afcd11917bdba0c599d71f51bd738b3a07dba6 --- /dev/null +++ b/950xxx/950787/mc.PhPy8EG_VINCIA_tt_hdamp258p75_nonallhad_valid.py @@ -0,0 +1,45 @@ +evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 1.5*top mass, A14 tune, at least one lepton filter, ME NNPDF30 NLO, VINCIA shower, DSID 410450 LHE files with Shower Weights added ' +evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton'] +evgenConfig.contact = [ 'dominic.hirschbuehl@cern.ch','m.fenton@cern.ch'] +evgenConfig.nEventsPerJob = 10000 +evgenConfig.inputFilesPerJob = 2 + +#-------------------------------------------------------------- +# Pythia8 showering +#-------------------------------------------------------------- +include('Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py') +include("Pythia8_i/Pythia8_Powheg_Main31.py") +include("Pythia8_i/Pythia8_ShowerWeights.py") + + +genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 1' ] +genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 1' ] + +genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ] +genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ] +genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ] +genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ] +genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ] +genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ] + + +genSeq.Pythia8.Commands += ['PartonShowers:model = 2'] +genSeq.Pythia8.Commands += ['Vincia:pTmaxMatch = 2'] # recommended in https://arxiv.org/pdf/2106.10987.pdf +genSeq.Pythia8.Commands += ['Vincia:ewMode = 0'] +genSeq.Pythia8.Commands += ['Vincia:interleaveResDec= on'] +genSeq.Pythia8.Commands += ['Vincia:helicityShower = off'] + +genSeq.Pythia8.Commands += ['BeamRemnants:primordialKThard = 0.4'] +genSeq.Pythia8.Commands += ['ColourReconnection:range = 1.75'] + +genSeq.Pythia8.Commands += ['MultipartonInteractions:alphaSvalue = 0.119'] + +genSeq.Pythia8.Commands += ['MultipartonInteractions:pT0Ref = 2.24'] +genSeq.Pythia8.Commands += ['SigmaProcess:alphaSvalue = 0.119'] + +#-------------------------------------------------------------- +# Event filter +#-------------------------------------------------------------- +include('GeneratorFilters/TTbarWToLeptonFilter.py') +filtSeq.TTbarWToLeptonFilter.NumLeptons = -1 +filtSeq.TTbarWToLeptonFilter.Ptcut = 0.