Commit 4f188fc3 authored by Nishu Nishu's avatar Nishu Nishu Committed by Joany Manjarres
Browse files

Excited e- and mu-lepton signal MC (DSIDs 801003-801114)

parent afb46583
Pipeline #3001482 passed with stages
in 1 minute and 7 seconds
evgenConfig.description = "Excited Electron production: eStarE -> eqQE, meStar=100 Lambda=10000, Pythia8 with the A14 NNPDF23LO tune"
evgenConfig.process = "excitedElectron + electron -> electron + electron + quark + quark"
evgenConfig.keywords = ["BSM","exotic","excitedElectron","electron","quark"]
evgenConfig.contact = ["Alexei Myagkov <alexei.miagkov@cern.ch>"]
evgenConfig.nEventsPerJob = 10000
# Excited Lepton ID
leptID = 4000011
# Excited lepton Mass (in GeV)
M_ExLep = 100
# Excited lepton Width (in GeV)
W_ExLep = 2.591575E-05
# Mass Scale parameter (Lambda, in GeV)
M_Lam = 10000
# Coupling constants
f = 1.0
fPrime = 1.0
# Branching ratios
br_GI_gamma=7.038577E-01
br_GI_Z=1.340343E-02
br_GI_W=2.799340E-01
br_CI_ddbar=3.834193E-04
br_CI_uubar=3.834541E-04
br_CI_ssbar=3.834502E-04
br_CI_ccbar=3.832110E-04
br_CI_bbbar=3.772494E-04
br_CI_ttbar=0.000000E+00
br_CI_eebar=2.556206E-04
br_CI_nuenuebar=1.278103E-04
br_CI_mumubar=1.278142E-04
br_CI_numunumubar=1.278103E-04
br_CI_tautaubar=1.271659E-04
br_CI_nutaunutaubar=1.278103E-04
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
if leptID > 0 :
SMleptID=leptID-4000000
if SMleptID% 2 == 1 :
WleptID=-24
SMPairleptID = SMleptID+1
else :
SMPairleptID = SMleptID-1
WleptID=24
else :
SMleptID=leptID+4000000
if SMleptID% 2 == 1 :
WleptID=24
SMPairleptID = SMleptID-1
else :
WleptID=-24
SMPairleptID = SMleptID+1
if leptID == 4000011 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2eStare = on" ]
elif leptID == 4000013 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2muStarmu = on"]
elif leptID == 4000015 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2tauStartau = on"]
elif leptID == 4000012 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nueStarnue = on" ]
elif leptID == 4000014 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2numuStarnumu = on"]
elif leptID == 4000016 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nutauStarnutau = on"]
else :
from AthenaCommon.Logging import logging
logging.error("Can not generate this process for particle "+str(leptID))
genSeq.Pythia8.Commands += [
str(leptID)+":m0 = "+str(M_ExLep),
str(leptID)+":mWidth = "+str(W_ExLep),
str(leptID)+":doForceWidth = on",
"ExcitedFermion:Lambda = "+str(M_Lam),
"ExcitedFermion:coupF = "+str(f),
"ExcitedFermion:coupFprime = "+str(fPrime),
"ExcitedFermion:coupFcol = "+str(f),
str(leptID)+":onMode = off",
str(leptID)+':oneChannel = 0 '+str(br_GI_gamma)+' 102 22 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_Z)+' 102 23 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_W)+' 102 '+str(WleptID)+' '+str(SMPairleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ddbar)+' 102 1 -1 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_uubar)+' 102 2 -2 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ssbar)+' 102 3 -3 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ccbar)+' 102 4 -4 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_bbbar)+' 102 5 -5 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ttbar)+' 102 6 -6 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_eebar)+' 102 11 -11 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nuenuebar)+' 102 12 -12 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_mumubar)+' 102 13 -13 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_numunumubar)+' 102 14 -14 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_tautaubar)+' 102 15 -15 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nutaunutaubar)+' 102 16 -16 '+str(SMleptID),
]
evgenConfig.description = "Excited Electron production: eStarE -> eqQE, meStar=200 Lambda=10000, Pythia8 with the A14 NNPDF23LO tune"
evgenConfig.process = "excitedElectron + electron -> electron + electron + quark + quark"
evgenConfig.keywords = ["BSM","exotic","excitedElectron","electron","quark"]
evgenConfig.contact = ["Alexei Myagkov <alexei.miagkov@cern.ch>"]
evgenConfig.nEventsPerJob = 10000
# Excited Lepton ID
leptID = 4000011
# Excited lepton Mass (in GeV)
M_ExLep = 200
# Excited lepton Width (in GeV)
W_ExLep = 4.624195E-04
# Mass Scale parameter (Lambda, in GeV)
M_Lam = 10000
# Coupling constants
f = 1.0
fPrime = 1.0
# Branching ratios
br_GI_gamma=3.155793E-01
br_GI_Z=1.036699E-01
br_GI_W=5.757110E-01
br_CI_ddbar=6.877305E-04
br_CI_uubar=6.876656E-04
br_CI_ssbar=6.876873E-04
br_CI_ccbar=6.875575E-04
br_CI_bbbar=6.848543E-04
br_CI_ttbar=0.000000E+00
br_CI_eebar=4.584798E-04
br_CI_nuenuebar=2.292291E-04
br_CI_mumubar=2.292291E-04
br_CI_numunumubar=2.292291E-04
br_CI_tautaubar=2.289263E-04
br_CI_nutaunutaubar=2.292291E-04
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
if leptID > 0 :
SMleptID=leptID-4000000
if SMleptID% 2 == 1 :
WleptID=-24
SMPairleptID = SMleptID+1
else :
SMPairleptID = SMleptID-1
WleptID=24
else :
SMleptID=leptID+4000000
if SMleptID% 2 == 1 :
WleptID=24
SMPairleptID = SMleptID-1
else :
WleptID=-24
SMPairleptID = SMleptID+1
if leptID == 4000011 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2eStare = on" ]
elif leptID == 4000013 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2muStarmu = on"]
elif leptID == 4000015 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2tauStartau = on"]
elif leptID == 4000012 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nueStarnue = on" ]
elif leptID == 4000014 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2numuStarnumu = on"]
elif leptID == 4000016 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nutauStarnutau = on"]
else :
from AthenaCommon.Logging import logging
logging.error("Can not generate this process for particle "+str(leptID))
genSeq.Pythia8.Commands += [
str(leptID)+":m0 = "+str(M_ExLep),
str(leptID)+":mWidth = "+str(W_ExLep),
str(leptID)+":doForceWidth = on",
"ExcitedFermion:Lambda = "+str(M_Lam),
"ExcitedFermion:coupF = "+str(f),
"ExcitedFermion:coupFprime = "+str(fPrime),
"ExcitedFermion:coupFcol = "+str(f),
str(leptID)+":onMode = off",
str(leptID)+':oneChannel = 0 '+str(br_GI_gamma)+' 102 22 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_Z)+' 102 23 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_W)+' 102 '+str(WleptID)+' '+str(SMPairleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ddbar)+' 102 1 -1 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_uubar)+' 102 2 -2 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ssbar)+' 102 3 -3 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ccbar)+' 102 4 -4 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_bbbar)+' 102 5 -5 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ttbar)+' 102 6 -6 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_eebar)+' 102 11 -11 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nuenuebar)+' 102 12 -12 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_mumubar)+' 102 13 -13 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_numunumubar)+' 102 14 -14 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_tautaubar)+' 102 15 -15 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nutaunutaubar)+' 102 16 -16 '+str(SMleptID),
]
evgenConfig.description = "Excited Electron production: eStarE -> eqQE, meStar=300 Lambda=10000, Pythia8 with the A14 NNPDF23LO tune"
evgenConfig.process = "excitedElectron + electron -> electron + electron + quark + quark"
evgenConfig.keywords = ["BSM","exotic","excitedElectron","electron","quark"]
evgenConfig.contact = ["Alexei Myagkov <alexei.miagkov@cern.ch>"]
evgenConfig.nEventsPerJob = 10000
# Excited Lepton ID
leptID = 4000011
# Excited lepton Mass (in GeV)
M_ExLep = 300
# Excited lepton Width (in GeV)
W_ExLep = 1.766914E-03
# Mass Scale parameter (Lambda, in GeV)
M_Lam = 10000
# Coupling constants
f = 1.0
fPrime = 1.0
# Branching ratios
br_GI_gamma=2.787346E-01
br_GI_Z=1.139331E-01
br_GI_W=5.973126E-01
br_CI_ddbar=1.366677E-03
br_CI_uubar=1.366620E-03
br_CI_ssbar=1.366563E-03
br_CI_ccbar=1.366620E-03
br_CI_bbbar=1.364186E-03
br_CI_ttbar=0.000000E+00
br_CI_eebar=9.111932E-04
br_CI_nuenuebar=4.555909E-04
br_CI_mumubar=4.555457E-04
br_CI_numunumubar=4.555909E-04
br_CI_tautaubar=4.553136E-04
br_CI_nutaunutaubar=4.555909E-04
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
if leptID > 0 :
SMleptID=leptID-4000000
if SMleptID% 2 == 1 :
WleptID=-24
SMPairleptID = SMleptID+1
else :
SMPairleptID = SMleptID-1
WleptID=24
else :
SMleptID=leptID+4000000
if SMleptID% 2 == 1 :
WleptID=24
SMPairleptID = SMleptID-1
else :
WleptID=-24
SMPairleptID = SMleptID+1
if leptID == 4000011 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2eStare = on" ]
elif leptID == 4000013 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2muStarmu = on"]
elif leptID == 4000015 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2tauStartau = on"]
elif leptID == 4000012 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nueStarnue = on" ]
elif leptID == 4000014 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2numuStarnumu = on"]
elif leptID == 4000016 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nutauStarnutau = on"]
else :
from AthenaCommon.Logging import logging
logging.error("Can not generate this process for particle "+str(leptID))
genSeq.Pythia8.Commands += [
str(leptID)+":m0 = "+str(M_ExLep),
str(leptID)+":mWidth = "+str(W_ExLep),
str(leptID)+":doForceWidth = on",
"ExcitedFermion:Lambda = "+str(M_Lam),
"ExcitedFermion:coupF = "+str(f),
"ExcitedFermion:coupFprime = "+str(fPrime),
"ExcitedFermion:coupFcol = "+str(f),
str(leptID)+":onMode = off",
str(leptID)+':oneChannel = 0 '+str(br_GI_gamma)+' 102 22 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_Z)+' 102 23 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_W)+' 102 '+str(WleptID)+' '+str(SMPairleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ddbar)+' 102 1 -1 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_uubar)+' 102 2 -2 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ssbar)+' 102 3 -3 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ccbar)+' 102 4 -4 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_bbbar)+' 102 5 -5 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ttbar)+' 102 6 -6 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_eebar)+' 102 11 -11 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nuenuebar)+' 102 12 -12 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_mumubar)+' 102 13 -13 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_numunumubar)+' 102 14 -14 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_tautaubar)+' 102 15 -15 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nutaunutaubar)+' 102 16 -16 '+str(SMleptID),
]
evgenConfig.description = "Excited Electron production: eStarE -> eqQE, meStar=400 Lambda=10000, Pythia8 with the A14 NNPDF23LO tune"
evgenConfig.process = "excitedElectron + electron -> electron + electron + quark + quark"
evgenConfig.keywords = ["BSM","exotic","excitedElectron","electron","quark"]
evgenConfig.contact = ["Alexei Myagkov <alexei.miagkov@cern.ch>"]
evgenConfig.nEventsPerJob = 10000
# Excited Lepton ID
leptID = 4000011
# Excited lepton Mass (in GeV)
M_ExLep = 400
# Excited lepton Width (in GeV)
W_ExLep = 4.385962E-03
# Mass Scale parameter (Lambda, in GeV)
M_Lam = 10000
# Coupling constants
f = 1.0
fPrime = 1.0
# Branching ratios
br_GI_gamma=2.661674E-01
br_GI_Z=1.164100E-01
br_GI_W=6.003928E-01
br_CI_ddbar=2.320130E-03
br_CI_uubar=2.320130E-03
br_CI_ssbar=2.320130E-03
br_CI_ccbar=2.320130E-03
br_CI_bbbar=2.317622E-03
br_CI_ttbar=1.807289E-05
br_CI_eebar=1.546821E-03
br_CI_nuenuebar=7.734221E-04
br_CI_mumubar=7.733765E-04
br_CI_numunumubar=7.734221E-04
br_CI_tautaubar=7.732169E-04
br_CI_nutaunutaubar=7.734221E-04
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
if leptID > 0 :
SMleptID=leptID-4000000
if SMleptID% 2 == 1 :
WleptID=-24
SMPairleptID = SMleptID+1
else :
SMPairleptID = SMleptID-1
WleptID=24
else :
SMleptID=leptID+4000000
if SMleptID% 2 == 1 :
WleptID=24
SMPairleptID = SMleptID-1
else :
WleptID=-24
SMPairleptID = SMleptID+1
if leptID == 4000011 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2eStare = on" ]
elif leptID == 4000013 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2muStarmu = on"]
elif leptID == 4000015 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2tauStartau = on"]
elif leptID == 4000012 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nueStarnue = on" ]
elif leptID == 4000014 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2numuStarnumu = on"]
elif leptID == 4000016 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nutauStarnutau = on"]
else :
from AthenaCommon.Logging import logging
logging.error("Can not generate this process for particle "+str(leptID))
genSeq.Pythia8.Commands += [
str(leptID)+":m0 = "+str(M_ExLep),
str(leptID)+":mWidth = "+str(W_ExLep),
str(leptID)+":doForceWidth = on",
"ExcitedFermion:Lambda = "+str(M_Lam),
"ExcitedFermion:coupF = "+str(f),
"ExcitedFermion:coupFprime = "+str(fPrime),
"ExcitedFermion:coupFcol = "+str(f),
str(leptID)+":onMode = off",
str(leptID)+':oneChannel = 0 '+str(br_GI_gamma)+' 102 22 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_Z)+' 102 23 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_W)+' 102 '+str(WleptID)+' '+str(SMPairleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ddbar)+' 102 1 -1 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_uubar)+' 102 2 -2 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ssbar)+' 102 3 -3 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ccbar)+' 102 4 -4 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_bbbar)+' 102 5 -5 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ttbar)+' 102 6 -6 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_eebar)+' 102 11 -11 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nuenuebar)+' 102 12 -12 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_mumubar)+' 102 13 -13 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_numunumubar)+' 102 14 -14 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_tautaubar)+' 102 15 -15 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nutaunutaubar)+' 102 16 -16 '+str(SMleptID),
]
evgenConfig.description = "Excited Electron production: eStarE -> eqQE, meStar=500 Lambda=10000, Pythia8 with the A14 NNPDF23LO tune"
evgenConfig.process = "excitedElectron + electron -> electron + electron + quark + quark"
evgenConfig.keywords = ["BSM","exotic","excitedElectron","electron","quark"]
evgenConfig.contact = ["Alexei Myagkov <alexei.miagkov@cern.ch>"]
evgenConfig.nEventsPerJob = 10000
# Excited Lepton ID
leptID = 4000011
# Excited lepton Mass (in GeV)
M_ExLep = 500
# Excited lepton Width (in GeV)
W_ExLep = 8.800819E-03
# Mass Scale parameter (Lambda, in GeV)
M_Lam = 10000
# Coupling constants
f = 1.0
fPrime = 1.0
# Branching ratios
br_GI_gamma=2.590782E-01
br_GI_Z=1.167505E-01
br_GI_W=5.979557E-01
br_CI_ddbar=3.528990E-03
br_CI_uubar=3.528762E-03
br_CI_ssbar=3.528876E-03
br_CI_ccbar=3.528649E-03
br_CI_bbbar=3.526717E-03
br_CI_ttbar=3.399684E-04
br_CI_eebar=2.352395E-03
br_CI_nuenuebar=1.176254E-03
br_CI_mumubar=1.176368E-03
br_CI_numunumubar=1.176254E-03
br_CI_tautaubar=1.176141E-03
br_CI_nutaunutaubar=1.176254E-03
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
if leptID > 0 :
SMleptID=leptID-4000000
if SMleptID% 2 == 1 :
WleptID=-24
SMPairleptID = SMleptID+1
else :
SMPairleptID = SMleptID-1
WleptID=24
else :
SMleptID=leptID+4000000
if SMleptID% 2 == 1 :
WleptID=24
SMPairleptID = SMleptID-1
else :
WleptID=-24
SMPairleptID = SMleptID+1
if leptID == 4000011 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2eStare = on" ]
elif leptID == 4000013 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2muStarmu = on"]
elif leptID == 4000015 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2tauStartau = on"]
elif leptID == 4000012 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nueStarnue = on" ]
elif leptID == 4000014 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2numuStarnumu = on"]
elif leptID == 4000016 :
genSeq.Pythia8.Commands += ["ExcitedFermion:qqbar2nutauStarnutau = on"]
else :
from AthenaCommon.Logging import logging
logging.error("Can not generate this process for particle "+str(leptID))
genSeq.Pythia8.Commands += [
str(leptID)+":m0 = "+str(M_ExLep),
str(leptID)+":mWidth = "+str(W_ExLep),
str(leptID)+":doForceWidth = on",
"ExcitedFermion:Lambda = "+str(M_Lam),
"ExcitedFermion:coupF = "+str(f),
"ExcitedFermion:coupFprime = "+str(fPrime),
"ExcitedFermion:coupFcol = "+str(f),
str(leptID)+":onMode = off",
str(leptID)+':oneChannel = 0 '+str(br_GI_gamma)+' 102 22 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_Z)+' 102 23 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_GI_W)+' 102 '+str(WleptID)+' '+str(SMPairleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ddbar)+' 102 1 -1 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_uubar)+' 102 2 -2 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ssbar)+' 102 3 -3 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ccbar)+' 102 4 -4 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_bbbar)+' 102 5 -5 '+str(SMleptID),
str(leptID)+':addChannel = 1 '+str(br_CI_ttbar)+' 102 6 -6 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_eebar)+' 102 11 -11 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nuenuebar)+' 102 12 -12 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_mumubar)+' 102 13 -13 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_numunumubar)+' 102 14 -14 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_tautaubar)+' 102 15 -15 '+str(SMleptID),
str(leptID)+':addChannel = 0 '+str(br_CI_nutaunutaubar)+' 102 16 -16 '+str(SMleptID),
]
evgenConfig.description = "Excited Electron production: eStarE -> eqQE, meStar=600 Lambda=10000, Pythia8 with the A14 NNPDF23LO tune"
evgenConfig.process = "excitedElectron + electron -> electron + electron + quark + quark"
evgenConfig.keywords = ["BSM","exotic","excitedElectron","electron","quark"]
evgenConfig.contact = ["Alexei Myagkov <alexei.miagkov@cern.ch>"]
evgenConfig.nEventsPerJob = 10000
# Excited Lepton ID
leptID = 4000011
# Excited lepton Mass (in GeV)
M_ExLep = 600
# Excited lepton Width (in GeV)
W_ExLep = 1.553405E-02
# Mass Scale parameter (Lambda, in GeV)
M_Lam = 10000
# Coupling constants
f = 1.0
fPrime = 1.0
# Branching ratios
br_GI_gamma=2.536364E-01
br_GI_Z=1.161320E-01
br_GI_W=5.926079E-01
br_CI_ddbar=4.974878E-03
br_CI_uubar=4.975779E-03
br_CI_ssbar=4.974749E-03
br_CI_ccbar=4.975135E-03
br_CI_bbbar=4.972625E-03
br_CI_ttbar=1.142394E-03
br_CI_eebar=3.316778E-03
br_CI_nuenuebar=1.658293E-03
br_CI_mumubar=1.658357E-03
br_CI_numunumubar=1.658293E-03
br_CI_tautaubar=1.658164E-03