Commit a5c40282 authored by Stephen Farry's avatar Stephen Farry
Browse files

dec files for ttbar, single top and Wbb powheg production

parent 3ed9f336
Pipeline #1795596 failed with stages
in 2 minutes and 5 seconds
# EventType: 41900073
#
# Descriptor: pp => (W- => ...) (t => b ...) ...
#
# NickName: ST_wtch_DR_t_1l17GeV
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/OneLepFromTop
#
# InsertPythonCode:
# from Configurables import Generation, Special, PowhegProduction
# Generation().addTool( Special )
# Generation().Special.addTool( PowhegProduction )
# Generation().Special.PowhegProduction.Commands += [ "ttype 1" ]
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "OneLepFromTop" )
# tracksInAcc = Generation().OneLepFromTop
# tracksInAcc.Code = " ( ( count ( isGoodBfromT ) > 1 ) & ( count ( isGoodLepton ) > 0 ) ) "
# tracksInAcc.Preambulo += [
# "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodBfromT = ( ( 6 == GABSID ) & GINTREE( GBEAUTY ) )"
# , "isGoodLepton = ( GLEPTON & GCHARGED & ( GTHETA < 400.0*mrad ) & ( GPT > 17*GeV ) )"
# ]
# EndInsertPythonCode
#
#
# ExtraOptions: Top_ST_wtch_DR
#
# Documentation: Single top Production through WT Channel. 1 lepton within LHCb acceptance with pT greater than 17GeV.
# EndDocumentation
#
# PhysicsWG: EW
# Tested: Yes
# Responsible: Stephen Farry
# Email: stephen.farry@cern.ch
# Date: 20200707
#
End
# EventType: 41900074
#
# Descriptor: pp => (W+ => ...) (t~ => b~ ...) ...
#
# NickName: ST_wtch_DR_tbar_1l17GeV
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/OneLepFromTop
#
# InsertPythonCode:
# from Configurables import Generation, Special, PowhegProduction
# Generation().addTool( Special )
# Generation().Special.addTool( PowhegProduction )
# Generation().Special.PowhegProduction.Commands += [ "ttype -1" ]
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "OneLepFromTop" )
# tracksInAcc = Generation().OneLepFromTop
# tracksInAcc.Code = " ( ( count ( isGoodBfromT ) > 1 ) & ( count ( isGoodLepton ) > 0 ) ) "
# tracksInAcc.Preambulo += [
# "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodBfromT = ( ( 6 == GABSID ) & GINTREE( GBEAUTY ) )"
# , "isGoodLepton = ( GLEPTON & GCHARGED & ( GTHETA < 400.0*mrad ) & ( GPT > 17*GeV ) )"
# ]
# EndInsertPythonCode
#
#
# ExtraOptions: Top_ST_wtch_DR
#
# Documentation: Single anti-top Production through WT Channel. 1 lepton within LHCb acceptance with pT greater than 17GeV.
# EndDocumentation
#
# PhysicsWG: EW
# Tested: Yes
# Responsible: Stephen Farry
# Email: stephen.farry@cern.ch
# Date: 20200707
#
End
# EventType: 42100072
#
# Descriptor: pp => (W- => mu- nu_mu) b b
#
# NickName: W_mumbb_powheg=mu17
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/OneLep
#
# InsertPythonCode:
# from Configurables import Generation, Special, PowhegProduction
# Generation().addTool( Special )
# Generation().Special.addTool( PowhegProduction )
# Generation().Special.PowhegProduction.Commands += [ "idvecbos -24", "vdecaymode 2" ]
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "OneLep" )
# tracksInAcc = Generation().OneLep
# tracksInAcc.Code = " ( ( count ( isGoodLepton ) > 0 ) ) "
# tracksInAcc.Preambulo += [
# "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodLepton = ( ( 13 == GABSID ) & GCHARGED & ( GTHETA < 400.0*mrad ) & ( GPT > 17*GeV ) )"
# ]
# EndInsertPythonCode
#
#
# ExtraOptions: Wbb_powheg
#
# Documentation: Wbb production with negative W decaying to muon. 1 muon within LHCb acceptance with pT greater than 17GeV.
# EndDocumentation
#
# PhysicsWG: EW
# Tested: Yes
# Responsible: Stephen Farry
# Email: stephen.farry@cern.ch
# Date: 20200707
#
End
# EventType: 42100071
#
# Descriptor: pp => (W+ => mu+ nu_mu) b b
#
# NickName: W_mupbb_powheg=mu17
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/OneLep
#
# InsertPythonCode:
# from Configurables import Generation, Special, PowhegProduction
# Generation().addTool( Special )
# Generation().Special.addTool( PowhegProduction )
# Generation().Special.PowhegProduction.Commands += [ "idvecbos 24", "vdecaymode 2" ]
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "OneLep" )
# tracksInAcc = Generation().OneLep
# tracksInAcc.Code = " ( ( count ( isGoodLepton ) > 0 ) ) "
# tracksInAcc.Preambulo += [
# "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodLepton = ( ( 13 == GABSID ) & GCHARGED & ( GTHETA < 400.0*mrad ) & ( GPT > 17*GeV ) )"
# ]
# EndInsertPythonCode
#
#
# ExtraOptions: Wbb_powheg
#
# Documentation: Wbb production with positive W decaying to muon. 1 muon within LHCb acceptance with pT greater than 17GeV.
#
# PhysicsWG: EW
# Tested: Yes
# Responsible: Stephen Farry
# Email: stephen.farry@cern.ch
# Date: 20200707
#
End
# EventType: 41900071
#
# Descriptor: pp => (t => b ...) (t~ => b~ ...) ...
#
# NickName: ttbar_powheg_1l17GeV
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/OneLepFromTop
#
# InsertPythonCode:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "OneLepFromTop" )
# tracksInAcc = Generation().OneLepFromTop
# tracksInAcc.Code = " ( ( count ( isGoodBfromT ) > 1 ) & ( count ( isGoodLepton ) > 0 ) ) "
# tracksInAcc.Preambulo += [
# "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodBfromT = ( ( 6 == GABSID ) & GINTREE( GBEAUTY ) )"
# , "isGoodLepton = ( GLEPTON & GCHARGED & ( GTHETA < 400.0*mrad ) & ( GPT > 17*GeV ) )"
# ]
# EndInsertPythonCode
#
#
# ExtraOptions: Top_ttbar_powheg
#
# Documentation: 1 lepton within LHCb acceptance with pT greater than 17GeV. Powheg performs inclusive top decays.
# EndDocumentation
#
# PhysicsWG: EW
# Tested: Yes
# Responsible: Stephen Farry
# Email: stephen.farry@cern.ch
# Date: 20200707
#
End
# EventType: 41900072
#
# Descriptor: pp => (t => b ...) (t~ => b~ ...) ...
#
# NickName: ttbar_powheg_2l17GeV
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/OneLepFromTop
#
# InsertPythonCode:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "OneLepFromTop" )
# tracksInAcc = Generation().OneLepFromTop
# tracksInAcc.Code = " ( ( count ( isGoodBfromT ) > 1 ) & ( count ( isGoodLepton ) > 2 ) ) "
# tracksInAcc.Preambulo += [
# "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodBfromT = ( ( 6 == GABSID ) & GINTREE( GBEAUTY ) )"
# , "isGoodLepton = ( GLEPTON & GCHARGED & ( GTHETA < 400.0*mrad ) & ( GPT > 17*GeV ) )"
# ]
# EndInsertPythonCode
#
#
# ExtraOptions: Top_ttbar_powheg_dilepton
#
# Documentation: 2 lepton within LHCb acceptance with pT greater than 17GeV. Dilepton decay forced by Powheg
# EndDocumentation
#
# PhysicsWG: EW
# Tested: Yes
# Responsible: Stephen Farry
# Email: stephen.farry@cern.ch
# Date: 20200707
#
End
# Powheg options for t tbar
from Configurables import Generation
from Gaudi.Configuration import *
Generation().PileUpTool = "FixedLuminosityForRareProcess"
importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
from Configurables import Special, PowhegProduction
Generation().addTool( Special )
Generation().Special.addTool( PowhegProduction )
# Powheg options.
Generation().Special.addTool(PowhegProduction)
Generation().Special.PowhegProduction.Process = "ST_wtch_DR"
Generation().Special.PowhegProduction.Commands += [
"lhans1 10770", # Change the first proton PDF.
"lhans2 10770", # Change the second proton PDF.
"wdecaymode 11100",
"min_w_mass 10d0",
"max_w_mass 1000d0",
"withnegweights 0",
'tdec/elbranching 0.108',
"hdamp 150d0", #default taken from testrun-lhc in Wbb package
"topdecaymode 11100",
"topmass 172.5",
"bmass 4.75",
"wmass 80.398",
"sthw2 0.2226459",
"alphaem_inv 132.3384323",
"charmthr 1.5",
"bottomthr 4.75",
"smartsig 1",
"withdamp 1"
]
# Powheg options for t tbar
from Configurables import Generation
from Gaudi.Configuration import *
Generation().PileUpTool = "FixedLuminosityForRareProcess"
importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
from Configurables import Special, PowhegProduction
Generation().addTool( Special )
Generation().Special.addTool( PowhegProduction )
# Powheg options.
Generation().Special.addTool(PowhegProduction)
Generation().Special.PowhegProduction.Process = "hvq"
Generation().Special.PowhegProduction.Commands += [
"lhans1 10770", # Change the first proton PDF.
"lhans2 10770", # Change the second proton PDF.
"qmass 172.5",
"topdecaymode 22222", # let it decay to all modes
"topmass 172.5",
"hdamp 258.75",
"tdec/wmass 80.4",
"tdec/wwidth 2.141",
"tdec/bmass 5",
"tdec/twidth 1.31",
"tdec/elbranching 0.108",
"tdec/emass 0.00051",
"tdec/mumass 0.1057",
"tdec/taumass 1.777",
"tdec/dmass 0.100",
"tdec/umass 0.100",
"tdec/smass 0.200",
"tdec/cmass 1.5",
"tdec/sin2cabibbo 0.051",
]
\ No newline at end of file
# Powheg options for t tbar
from Configurables import Generation
from Gaudi.Configuration import *
Generation().PileUpTool = "FixedLuminosityForRareProcess"
importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
from Configurables import Special, PowhegProduction
Generation().addTool( Special )
Generation().Special.addTool( PowhegProduction )
# Powheg options.
Generation().Special.addTool(PowhegProduction)
Generation().Special.PowhegProduction.Process = "hvq"
Generation().Special.PowhegProduction.Commands += [
"lhans1 10770", # Change the first proton PDF.
"lhans2 10770", # Change the second proton PDF.
"qmass 172.5",
"topdecaymode 22200", # let it decay to all leptonic modes
"topmass 172.5",
"hdamp 258.75",
"tdec/wmass 80.4",
"tdec/wwidth 2.141",
"tdec/bmass 5",
"tdec/twidth 1.31",
"tdec/elbranching 0.108",
"tdec/emass 0.00051",
"tdec/mumass 0.1057",
"tdec/taumass 1.777",
"tdec/dmass 0.100",
"tdec/umass 0.100",
"tdec/smass 0.200",
"tdec/cmass 1.5",
"tdec/sin2cabibbo 0.051",
]
# Powheg options for t tbar
from Configurables import Generation
from Gaudi.Configuration import *
Generation().PileUpTool = "FixedLuminosityForRareProcess"
importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
from Configurables import Special, PowhegProduction
Generation().addTool( Special )
Generation().Special.addTool( PowhegProduction )
# Powheg options.
Generation().Special.addTool(PowhegProduction)
Generation().Special.PowhegProduction.Process = "Wbb_dec"
Generation().Special.PowhegProduction.Commands += [
"lhans1 10770", # Change the first proton PDF.
"lhans2 10770", # Change the second proton PDF.
"idvecbos 24",
"vdecaymode 2",
"min_w_mass 10d0",
"max_w_mass 1000d0",
"withnegweights 0",
"hdamp 150d0" #default taken from testrun-lhc in Wbb package
]
\ No newline at end of file
Markdown is supported
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