tt Powheg+MadSpin with special b UFO model

Make sure that all the items in the checklists are duly followed:
follow the instructions of each checklist item and report potential issues to the repository mainteners.

JIRA ticket

AGENE-2229

Checklist for approval

This checklist should be normally filled by MC contacts (but it can also be done by the person opening the merge request).
If something is unclear get in touch with the repository maintainers.

• I have followed the guidelines in McSampleRequestProcedure
• I have provided the ATLMCPROD JIRA ticket where the request has been discussed
• The request has been approved by the PA/CP group
• The pipeline has run and the status is green (still check the CI job output - if you see any errors report to the repository maintainers)
• If the pipeline is orange: check the CI job output to sort out the issues. If the issue is unclear contact the repository maintainers
• CI jobs have not been skipped (if jobs have been skipped, the requesters should provide a reason for skipping and confirmation from the PMG conveners is necessary before merging)
• Check that no log.generate.short files are included in the commit to be merged to master. If such files are present, it indicates that something went wrong in the pipeline. Check the job logs and contact the package maintainers if needed.
• Check that the title of the MR is descriptive enough (if not ask the requester to modify it)
• Check that the DSID range assigned is correct (500000-999999) and corresponds to right generator block

950xxx/950744/mc.PhPy8EG_ttbar_madspin_smspecialb_nonallhad_valid.py

+#--------------------------------------------------------------
+# EVGEN configuration
+#--------------------------------------------------------------
+evgenConfig.description = 'POWHEG+Pythia8 ttbar production with MadSpin decay, Powheg hdamp equal 1.5*top mass, A14 tune and single lepton filter.'
+evgenConfig.keywords    = [ 'SM', 'top', 'ttbar', 'lepton']
+evgenConfig.contact     = [ 'dominic.hirschbuehl@cern.ch' ]
+
+import os, sys, glob
+
+include('PowhegControl/PowhegControl_tt_Common.py')
+
+# Initial settings
+PowhegConfig.decay_mode = "t t~ > undecayed"
+
+PowhegConfig.hdamp        = 258.75
+#PowhegConfig.mu_F         = 1.0
+#PowhegConfig.mu_R         = 1.0
+#PowhegConfig.PDF          = 260000
+
+PowhegConfig.MadSpin_taus_are_leptons = True
+PowhegConfig.MadSpin_decays= ["decay t > w+ B, w+ > all all", "decay t~ > w- B~, w- > all all"]
+PowhegConfig.MadSpin_process = "generate p p > t t~ [QCD]"
+PowhegConfig.MadSpin_mode = "full"
+PowhegConfig.MadSpin_model = "/cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/sm_specialb"
+PowhegConfig.MadSpin_nFlavours = 5
+os.system('get_files -jo param_card_sm_specialb.dat')
+PowhegConfig.MadSpin_paramcard = "./param_card_sm_specialb.dat"
+
+# new line added to match MC@NLO config
+PowhegConfig.bwcutoff = 50
+
+PowhegConfig.nEvents     *= 2.2 # safety factor
+PowhegConfig.generate()
+
+for f in glob.glob("*.events"):
+    infile = f
+    f1 = open( infile )
+    newfile = infile+'.temp'
+    f2 = open(newfile,'w')
+    for line in f1:
+        if line.startswith(' 100000005'):
+            f2.write(line.replace(' 100000005','        5'))
+        elif line.startswith(' -100000005'):
+            f2.write(line.replace(' -100000005','       -5'))
+        else:
+            f2.write(line)
+
+    f1.close()
+    f2.close()
+    os.system('mv %s %s '%(infile, infile+'.old') )
+    os.system('mv %s %s '%(newfile, infile) )
+
+
+#--------------------------------------------------------------
+# Pythia8 showering
+#--------------------------------------------------------------
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+include("Pythia8_i/Pythia8_Powheg_Main31.py")
+genSeq.Pythia8.Commands += [ 'POWHEG:pTHard = 0' ]
+genSeq.Pythia8.Commands += [ 'POWHEG:NFinal = 2' ]
+genSeq.Pythia8.Commands += [ 'POWHEG:pTdef = 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' ]
+
+
+#--------------------------------------------------------------
+# Event filter
+#--------------------------------------------------------------
+include('GeneratorFilters/TTbarWToLeptonFilter.py')
+filtSeq.TTbarWToLeptonFilter.NumLeptons = -1
+filtSeq.TTbarWToLeptonFilter.Ptcut = 0.
+