Commit c0d37340 authored by Quentin Buat's avatar Quentin Buat Committed by Jan Kretzschmar
Browse files

New JO for gamma*->tautau (tau CP developments) in R21.6

parent 553af14b
Pipeline #2986965 passed with stages
in 1 minute and 9 seconds
# JO for Pythia 8 gamma*->tau(had)tau(had)
evgenConfig.description = "Pythia 8 gamma*->tau(had)tau(had) production with A14 NNPDF23LO tune"
evgenConfig.keywords = ["SM", "drellYan", "electroweak", "2tau"]
evgenConfig.contact = ["quentin.buat@cern.ch"]
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
genSeq.Pythia8.Commands += ["WeakSingleBoson:ffbar2gmZ = on",
"WeakZ0:gmZmode = 1", # photon only
"23:onMode = off", # turn off all decays modes
"23:onIfAny = 15", # turn on the tautau decay mode
"15:offIfAny = 12 14", # turn off leptonic tau decays
"-15:offIfAny = 12 14", # turn off leptonic tau decays
"PhaseSpace:mHatMin = 60.", # lower invariant mass
"PhaseSpace:mHatMax = 7000."] # upper invariant mass
# mHatReweight hook
# -----------------
# Modifies *slope* (by mHat^(slope)) of falling gamma mass
# distribution up to mass *mHatConstMin*, after which a
# constant cross section is used. The SM gamma distribution
# is parameterised as (mHat/cme)^(p1) * (1 - mHat/cme)^(p2)
genSeq.Pythia8.UserHooks += ["mHatReweight"]
genSeq.Pythia8.Commands += [
"mHatReweight:Slope = 2.0", # Multiplier to falling gamma: mHat^(slope)
"mHatReweight:mHatConstMin = 1500.", # mHat to start constant
"mHatReweight:p1 = -4.18717",
"mHatReweight:p2 = 11.5811",
]
evgenConfig.nEventsPerJob = 10000
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