PowhegPythia8EvtGen_H2a4X_ctauY.py

Go to the documentation of this file.

#--------------------------------------------------------------
# Modify the events.lhe, since Pythia doesn't like to decay the
# SM higgs to BSM products: 25 --> 35
#--------------------------------------------------------------
 
#evgenConfig.inputfilecheck = "TXT"
#evgenConfig.inputfilecheck = "merged_lhef"
 
import os, sys, glob
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('      25     1'):
            f2.write(line.replace('      25     1','      35     1'))
        else:
            f2.write(line)
    f1.close()
    f2.close()
    os.system('mv %s %s '%(infile, infile+'.old') )
    os.system('mv %s %s '%(newfile, infile) )
 
nProcess = -1 # 0: WpH, 1: WmH, 2: ZH, 3: ggZH
ma = 0.       # 15, 55 GeV
adecay = 0    # 2, 5, 15
ctau = 0.     # 1, 10, 100 mm
 
# a->bb
if runArgs.runNumber==309851: # WpH, ma=15GeV, a->bb, ctau=1mm
    nProcess = 0
    ma = 15.
    ctau = 1.
    adecay = 5
elif runArgs.runNumber==309852: # WpH, ma=55GeV, a->bb, ctau=1mm
    nProcess = 0
    ma = 55.
    ctau = 1.
    adecay = 5
elif runArgs.runNumber==309853: # WpH, ma=15GeV, a->bb, ctau=10mm
    nProcess = 0
    ma = 15.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==309854: # WpH, ma=55GeV, a->bb, ctau=10mm
    nProcess = 0
    ma = 55.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==309855: # WpH, ma=15GeV, a->bb, ctau=100mm
    nProcess = 0
    ma = 15.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==309856: # WpH, ma=55GeV, a->bb, ctau=100mm
    nProcess = 0
    ma = 55.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==309857: # WmH, ma=15GeV, a->bb, ctau=1mm
    nProcess = 1
    ma = 15.
    ctau = 1.
    adecay = 5
elif runArgs.runNumber==309858: # WmH, ma=55GeV, a->bb, ctau=1mm
    nProcess = 1
    ma = 55.
    ctau = 1.
    adecay = 5
elif runArgs.runNumber==309859: # WmH, ma=15GeV, a->bb, ctau=10mm
    nProcess = 1
    ma = 15.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==309860: # WmH, ma=55GeV, a->bb, ctau=10mm
    nProcess = 1
    ma = 55.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==309861: # WmH, ma=15GeV, a->bb, ctau=100mm
    nProcess = 1
    ma = 15.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==309862: # WmH, ma=55GeV, a->bb, ctau=100mm
    nProcess = 1
    ma = 55.
    ctau = 100.
    adecay = 5
 
# a->tautau
elif runArgs.runNumber==309863: # WpH, ma=55GeV, a->tautau, ctau=100mm
    nProcess = 0
    ma = 55.
    ctau = 100.
    adecay = 15
elif runArgs.runNumber==309864: # WmH, ma=55GeV, a->tautau, ctau=100mm
    nProcess = 1
    ma = 55.
    ctau = 100.
    adecay = 15
 
# a->uu
elif runArgs.runNumber==309865: # WpH, ma=55GeV, a->uu, ctau=100mm
    nProcess = 0
    ma = 55.
    ctau = 100.
    adecay = 2
elif runArgs.runNumber==309866: # WmH, ma=55GeV, a->uu, ctau=100mm
    nProcess = 1
    ma = 55.
    ctau = 100.
    adecay = 2
 
# ZH, a->bb
elif runArgs.runNumber==312932: # ZH, ma=15GeV, a->bb, ctau=10mm
    nProcess = 2
    ma = 15.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312933: # ZH, ma=15GeV, a->bb, ctau=100mm
    nProcess = 2
    ma = 15.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312934: # ZH, ma=15GeV, a->bb, ctau=1000mm
    nProcess = 2
    ma = 15.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312935: # ZH, ma=25GeV, a->bb, ctau=10mm
    nProcess = 2
    ma = 25.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312936: # ZH, ma=25GeV, a->bb, ctau=100mm
    nProcess = 2
    ma = 25.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312937: # ZH, ma=25GeV, a->bb, ctau=1000mm
    nProcess = 2
    ma = 25.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312938: # ZH, ma=35GeV, a->bb, ctau=10mm
    nProcess = 2
    ma = 35.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312939: # ZH, ma=35GeV, a->bb, ctau=100mm
    nProcess = 2
    ma = 35.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312940: # ZH, ma=35GeV, a->bb, ctau=1000mm
    nProcess = 2
    ma = 35.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312941: # ZH, ma=55GeV, a->bb, ctau=10mm
    nProcess = 2
    ma = 55.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312942: # ZH, ma=55GeV, a->bb, ctau=100mm
    nProcess = 2
    ma = 35.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312943: # ZH, ma=55GeV, a->bb, ctau=1000mm
    nProcess = 2
    ma = 55.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312944: # ggZH, ma=15GeV, a->bb, ctau=10mm
    nProcess = 3
    ma = 15.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312945: # ggZH, ma=15GeV, a->bb, ctau=100mm
    nProcess = 3
    ma = 15.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312946: # ggZH, ma=15GeV, a->bb, ctau=1000mm
    nProcess = 3
    ma = 15.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312947: # ggZH, ma=25GeV, a->bb, ctau=10mm
    nProcess = 3
    ma = 25.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312948: # ggZH, ma=25GeV, a->bb, ctau=100mm
    nProcess = 3
    ma = 25.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312949: # ggZH, ma=25GeV, a->bb, ctau=1000mm
    nProcess = 3
    ma = 25.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312950: # ggZH, ma=35GeV, a->bb, ctau=10mm
    nProcess = 3
    ma = 35.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312951: # ggZH, ma=35GeV, a->bb, ctau=100mm
    nProcess = 3
    ma = 35.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312952: # ggZH, ma=35GeV, a->bb, ctau=1000mm
    nProcess = 3
    ma = 35.
    ctau = 1000.
    adecay = 5
elif runArgs.runNumber==312953: # ggZH, ma=55GeV, a->bb, ctau=10mm
    nProcess = 3
    ma = 55.
    ctau = 10.
    adecay = 5
elif runArgs.runNumber==312954: # ggZH, ma=55GeV, a->bb, ctau=100mm
    nProcess = 3
    ma = 55.
    ctau = 100.
    adecay = 5
elif runArgs.runNumber==312955: # ggZH, ma=55GeV, a->bb, ctau=1000mm
    nProcess = 3
    ma = 55.
    ctau = 1000.
    adecay = 5
 
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include('Pythia8_i/Pythia8_AZNLO_CTEQ6L1_EvtGen_Common.py')
#--------------------------------------------------------------
# Pythia8 main31 matching
#--------------------------------------------------------------
 
if nProcess==0:
    genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 3']
elif nProcess==1:
    genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 3']
elif nProcess==2:
    genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 3']
elif nProcess==3:
    genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 2']
 
#--------------------------------------------------------------
# Higgs->bbar at Pythia8
#--------------------------------------------------------------
genSeq.Pythia8.Commands += [
                            'Higgs:useBSM = on',
 
                            '35:m0 = 125',
                            '35:mWidth = 0.00407',
                            '35:doForceWidth = on',
                            '35:onMode = off',
                            '35:onIfMatch = 36 36', # h->aa
 
                            '36:onMode = off', # decay of the a
                            '36:onIfAny = %d' % adecay, # decay a->XX
                            '36:m0 = %.1f' % ma, #scalar mass
                            '36:mMin = 0',
                            '36:tau0 = %.1f' % ctau, #nominal proper lifetime (in mm/c)
                            ]
 
genSeq.Pythia8.Commands = [i for i in genSeq.Pythia8.Commands if (("limitTau0" not in i) and ("tau0Max" not in i))]
genSeq.Pythia8.Commands += [
                            'ParticleDecays:tau0Max = 100000.0',
                            'ParticleDecays:limitTau0 = off'
                           ]
 
testSeq.TestHepMC.MaxTransVtxDisp = 100000000 #in mm
testSeq.TestHepMC.MaxVtxDisp = 100000000 #in mm
 
 
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
if nProcess==0:
    if adecay==5:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+W+jet->l+vbbbarbbbar production"
        evgenConfig.process = "WpH, H->2a->4b, W->lv"
    elif adecay==15:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+W+jet->l+vtau+tau-tau+tau- production"
        evgenConfig.process = "WpH, H->2a->4tau, W->lv"
    elif adecay==2:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+W+jet->l+vuubaruubar production"
        evgenConfig.process = "WpH, H->2a->4u, W->lv"
elif nProcess==1:
    if adecay==5:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+W+jet->l-vbbbarbbbar production"
        evgenConfig.process = "WmH, H->2a->4b, W->lv"
    elif adecay==15:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+W+jet->l-vtau+tau-tau+tau- production"
        evgenConfig.process = "WmH, H->2a->4tau, W->lv"
    elif adecay==2:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+W+jet->l-vuubaruubar production"
        evgenConfig.process = "WmH, H->2a->4u, W->lv"
elif nProcess==2:
    if adecay==5:
        evgenConfig.description = "POWHEG+MiNLO+Pythia8 H+Z+jet->l+l-bbbarbbbar production"
        evgenConfig.process     = "ZH, H->aa->4b, Z->ll"
elif nProcess==3:
    if adecay==5:
        evgenConfig.description = "POWHEG+Pythia8 gg->H+Z->l+l-bbbarbbbar production"
        evgenConfig.process     = "gg->ZH, H->aa->4b, Z->ll"
evgenConfig.keywords    = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
evgenConfig.contact     = [ 'manfredi.ronzani@cern.ch' ]
evgenConfig.minevents = 50