Page MenuHomeHEPForge

No OneTemporary

diff --git a/src/LHC-TTBA.in b/src/LHC-TTBA.in
new file mode 100644
--- /dev/null
+++ b/src/LHC-TTBA.in
@@ -0,0 +1,142 @@
+ ##################################################
+# Example generator for the Leptoquark model
+# in hadron collisions
+# The best way to use this is to make your own
+# copy of this file and edit that as you require.
+#
+# The first section loads the model file which
+# does not contain anything that users need to touch.
+#
+# The second section contains the user settings.
+###################################################
+read TTBA.model
+
+#Set up semi-leptonic (e or mu) top/hadronic topbar
+
+set /Herwig/Particles/t:Synchronized Not_synchronized
+set /Herwig/Particles/tbar:Synchronized Not_synchronized
+
+set /Herwig/Particles/t/t->b,bbar,c;:OnOff Off
+set /Herwig/Particles/t/t->b,c,dbar;:OnOff Off
+set /Herwig/Particles/t/t->b,c,sbar;:OnOff Off
+set /Herwig/Particles/t/t->b,sbar,u;:OnOff Off
+set /Herwig/Particles/t/t->b,u,dbar;:OnOff Off
+set /Herwig/Particles/t/t->nu_e,e+,b;:OnOff On
+set /Herwig/Particles/t/t->nu_mu,mu+,b;:OnOff On
+set /Herwig/Particles/t/t->nu_tau,tau+,b;:OnOff Off
+set /Herwig/Particles/tbar/tbar->b,bbar,cbar;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,cbar,d;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,cbar,s;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,s,ubar;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,ubar,d;:OnOff On
+set /Herwig/Particles/tbar/tbar->nu_ebar,e-,bbar;:OnOff Off
+set /Herwig/Particles/tbar/tbar->nu_mubar,mu-,bbar;:OnOff Off
+set /Herwig/Particles/tbar/tbar->nu_taubar,tau-,bbar;:OnOff Off
+
+cd /Herwig/Generators
+
+##################################################
+# LHC physics parameters (override defaults here)
+##################################################
+
+# Intrinsic pT tune extrapolated to LHC energy
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
+
+##################################################
+# Energy-dependent MPI parameters
+##################################################
+
+################################################################################
+# NOTE: Set the correct parameters for the collision energy of your choice
+# here.
+################################################################################
+
+#######################
+## sqrt(s) = 900 GeV ##
+#######################
+#set LHCGenerator:EventHandler:LuminosityFunction:Energy 900.0
+#set /Herwig/UnderlyingEvent/KtCut:MinKT 1.86
+#set /Herwig/UnderlyingEvent/UECuts:MHatMin 3.72
+
+########################
+## sqrt(s) = 2760 GeV ##
+########################
+#set LHCGenerator:EventHandler:LuminosityFunction:Energy 2760.0
+#set /Herwig/UnderlyingEvent/KtCut:MinKT 2.62
+#set /Herwig/UnderlyingEvent/UECuts:MHatMin 5.24
+
+########################
+## sqrt(s) = 7000 GeV ##
+########################
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 7000.0
+set /Herwig/UnderlyingEvent/KtCut:MinKT 3.06
+set /Herwig/UnderlyingEvent/UECuts:MHatMin 6.12
+
+#########################
+## sqrt(s) = 14000 GeV ##
+#########################
+#set LHCGenerator:EventHandler:LuminosityFunction:Energy 14000.0
+#set /Herwig/UnderlyingEvent/KtCut:MinKT 3.53
+#set /Herwig/UnderlyingEvent/UECuts:MHatMin 7.06
+
+
+##################################################
+#
+# This section contains the user defined settings
+#
+##################################################
+# Example hard process: Incoming proton, outgoing leptoquarks
+cd /Herwig/NewPhysics
+
+set HPConstructor:Processes Exclusive
+
+insert HPConstructor:Incoming 0 /Herwig/Particles/u
+insert HPConstructor:Incoming 1 /Herwig/Particles/ubar
+insert HPConstructor:Incoming 2 /Herwig/Particles/d
+insert HPConstructor:Incoming 3 /Herwig/Particles/dbar
+insert HPConstructor:Incoming 4 /Herwig/Particles/s
+insert HPConstructor:Incoming 5 /Herwig/Particles/sbar
+insert HPConstructor:Incoming 6 /Herwig/Particles/c
+insert HPConstructor:Incoming 7 /Herwig/Particles/cbar
+insert HPConstructor:Incoming 8 /Herwig/Particles/b
+insert HPConstructor:Incoming 9 /Herwig/Particles/bbar
+insert HPConstructor:Incoming 10 /Herwig/Particles/g
+
+# Comment out the following line if Axigluon model is selected
+insert /Herwig/NewPhysics/HPConstructor:Excluded 0 /Herwig/Particles/Ag
+
+insert HPConstructor:Outgoing 0 /Herwig/Particles/t
+insert HPConstructor:Outgoing 1 /Herwig/Particles/tbar
+
+#This is available for comparison to the LO Hw++ result
+#cd /Herwig/MatrixElements
+#insert SimpleQCD:MatrixElements[0] MEHeavyQuark
+
+cd /Herwig/Shower
+set SplittingGenerator:ISR Yes
+set SplittingGenerator:FSR Yes
+
+cd /Herwig/EventHandlers
+#set LHCHandler:CascadeHandler NULL
+#set LHCHandler:HadronizationHandler NULL
+#set LHCHandler:DecayHandler NULL
+
+# Other parameters for run
+cd /Herwig/Generators
+
+# Uncomment to turn off MPI
+#set LHCGenerator:EventHandler:CascadeHandler:MPIHandler NULL
+
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:PrintEvent 100
+set LHCGenerator:MaxErrors 10000
+
+#HepMC analysis
+#insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+#set /Herwig/Analysis/HepMCFile:PrintEvent 10000
+#set /Herwig/Analysis/HepMCFile:Format GenEvent
+#set /Herwig/Analysis/HepMCFile:Units GeV_mm
+
+saverun LHC-TTBA LHCGenerator
diff --git a/src/TTBA.model b/src/TTBA.model
--- a/src/TTBA.model
+++ b/src/TTBA.model
@@ -1,159 +1,159 @@
##################################################
# Common setup for Leptoquark models
#
-# See LHC-LQ.in or ILC-TTBA.in for example usage
+# See LHC-TTBA.in or ILC-TTBA.in for example usage
#
# This file does not contain anything that
# users need to touch. User settings are in
# ???-TTBA.in
#
###################################################
#
# Particle Data object for the new resonances
#
###################################################
cd /Herwig/Particles
# W prime model
create /ThePEG/ParticleData Wp
setup Wp 34 Wp 400.0 0.0 0.0 0.0 3 0 3 0
create /ThePEG/ParticleData Wpbar
setup Wpbar -34 Wpbar 400.0 0.0 0.0 0.0 -3 0 3 0
makeanti Wp Wpbar
# Abelian Z prime model
create /ThePEG/ParticleData Zp
setup Zp 32 Zp 700.0 0.0 0.0 0.0 0 0 3 0
# Axigluon model
create /ThePEG/ParticleData Ag
setup Ag 63 Ag 2000.0 0.0 0 0.0 0 8 3 0
# Non-Abelian SU(2)_X model
create /ThePEG/ParticleData Vz
setup Vz 70 Vz 80.0 0.0 0.0 0.0 0 0 3 0
create /ThePEG/ParticleData Yp
setup Yp 71 Yp 200.0 0.0 0.0 0.0 0 0 3 0
create /ThePEG/ParticleData Ym
setup Ym -71 Ym 200.0 0.0 0.0 0.0 0 0 3 0
makeanti Yp Ym
###################################################
#
# Main directory and model object
#
###################################################
mkdir /Herwig/NewPhysics/TTbA
cd /Herwig/NewPhysics/TTbA
create Herwig::TTbAModel Model HwTTbAModel.so
# SM couplings
set Model:QCD/RunningAlphaS /Herwig/AlphaS
set Model:EW/RunningAlphaEM /Herwig/AlphaEM
set Model:EW/CKM /Herwig/CKM
set Model:RunningMass /Herwig/RunningMass
# Select the active model (0: W', 1: Z', 2: Axigluon, 3: non-Abelian SU(2)_X)
-set Model:modelselect 3
+set Model:modelselect 0
# W prime couplings
set Model:WPTDLCoupling 0.0
set Model:WPTDRCoupling 2.0
# Z prime couplings
set Model:ZPTULCoupling 0.0
set Model:ZPTURCoupling 1.8
set Model:ZPUULCoupling 0.0
set Model:ZPUURCoupling 0.1
set Model:ZPCCLCoupling 0.0
set Model:ZPCCRCoupling 0.1
# Axigluon couplings
set Model:AGQQLCoupling -3.0
set Model:AGQQRCoupling 3.0
set Model:AGTTLCoupling 3.0
set Model:AGTTRCoupling -3.0
-# SU(2)_X parameters
+# SU(2)_X parameters (see: arXiv:1103.4835)
set Model:SU2Xcostheta 0.95
set Model:SU2Xalpha 0.060
###################################################
#
# Vertices
#
###################################################
# create model vertices
mkdir /Herwig/Vertices/TTbA
cd /Herwig/Vertices/TTbA
library HwTTbAModel.so
create Herwig::TTbAModelWPTDVertex TTbA_WPTDVertex
create Herwig::TTbAModelZPQQVertex TTbA_ZPQQVertex
create Herwig::TTbAModelAGQQVertex TTbA_AGQQVertex
create Herwig::TTbAModelSU2XVertex TTbA_SU2XVertex
cd /Herwig/NewPhysics/TTbA
# SM vertices
set Model:Vertex/FFZ /Herwig/Vertices/FFZVertex
set Model:Vertex/FFW /Herwig/Vertices/FFWVertex
set Model:Vertex/FFH /Herwig/Vertices/FFHVertex
set Model:Vertex/FFG /Herwig/Vertices/FFGVertex
set Model:Vertex/FFP /Herwig/Vertices/FFPVertex
set Model:Vertex/GGG /Herwig/Vertices/GGGVertex
set Model:Vertex/GGGG /Herwig/Vertices/GGGGVertex
set Model:Vertex/WWH /Herwig/Vertices/WWHVertex
set Model:Vertex/WWW /Herwig/Vertices/WWWVertex
set Model:Vertex/WWWW /Herwig/Vertices/WWWWVertex
set Model:Vertex/HGG /Herwig/Vertices/HGGVertex
set Model:Vertex/HHH /Herwig/Vertices/HHHVertex
set Model:Vertex/WWHH /Herwig/Vertices/WWHHVertex
set Model:Vertex/HHH /Herwig/Vertices/HHHVertex
set Model:Vertex/HPP /Herwig/Vertices/HPPVertex
# model vertices
set Model:Vertex/WPTD /Herwig/Vertices/TTbA/TTbA_WPTDVertex
set Model:Vertex/ZPQQ /Herwig/Vertices/TTbA/TTbA_ZPQQVertex
set Model:Vertex/AGQQ /Herwig/Vertices/TTbA/TTbA_AGQQVertex
set Model:Vertex/SU2X /Herwig/Vertices/TTbA/TTbA_SU2XVertex
###################################################
#
# Set up spin correlation Decayers
#
###################################################
cd /Herwig/NewPhysics
set TwoBodyDC:CreateDecayModes Yes
set ThreeBodyDC:CreateDecayModes No
insert NewModel:DecayParticles 0 /Herwig/Particles/Wp
insert NewModel:DecayParticles 1 /Herwig/Particles/Zp
insert NewModel:DecayParticles 2 /Herwig/Particles/Ag
insert NewModel:DecayParticles 3 /Herwig/Particles/Yp
insert NewModel:DecayParticles 4 /Herwig/Particles/Ym
insert NewModel:DecayParticles 5 /Herwig/Particles/Vz
###################################################
# Set up the model framework
###################################################
set TTbA/Model:ModelGenerator NewModel
###################################################
#
# Choose Model over SM
#
###################################################
cd /Herwig/Generators
set LEPGenerator:StandardModelParameters /Herwig/NewPhysics/TTbA/Model
set LHCGenerator:StandardModelParameters /Herwig/NewPhysics/TTbA/Model
diff --git a/src/TVT-TTBA.in b/src/TVT-TTBA.in
--- a/src/TVT-TTBA.in
+++ b/src/TVT-TTBA.in
@@ -1,135 +1,101 @@
-##################################################
+ ##################################################
# Example generator for the Leptoquark model
# in hadron collisions
# The best way to use this is to make your own
# copy of this file and edit that as you require.
#
# The first section loads the model file which
# does not contain anything that users need to touch.
#
# The second section contains the user settings.
###################################################
read TTBA.model
+#Set up semi-leptonic (e or mu) top/hadronic topbar
+
set /Herwig/Particles/t:Synchronized Not_synchronized
set /Herwig/Particles/tbar:Synchronized Not_synchronized
-
-#set /Herwig/Particles/t/t->b,bbar,c;:OnOff On
-#set /Herwig/Particles/t/t->b,c,dbar;:OnOff On
-#set /Herwig/Particles/t/t->b,c,sbar;:OnOff On
-#set /Herwig/Particles/t/t->b,sbar,u;:OnOff On
-#set /Herwig/Particles/t/t->b,u,dbar;:OnOff On
-#set /Herwig/Particles/t/t->nu_e,e+,b;:OnOff Off
-#set /Herwig/Particles/t/t->nu_mu,mu+,b;:OnOff Off
-#set /Herwig/Particles/t/t->nu_tau,tau+,b;:OnOff Off
-#set /Herwig/Particles/tbar/tbar->b,bbar,cbar;:OnOff Off
-#set /Herwig/Particles/tbar/tbar->bbar,cbar,d;:OnOff Off
-#set /Herwig/Particles/tbar/tbar->bbar,cbar,s;:OnOff Off
-#set /Herwig/Particles/tbar/tbar->bbar,s,ubar;:OnOff Off
-#set /Herwig/Particles/tbar/tbar->bbar,ubar,d;:OnOff Off
-#set /Herwig/Particles/tbar/tbar->nu_ebar,e-,bbar;:OnOff On
-#set /Herwig/Particles/tbar/tbar->nu_mubar,mu-,bbar;:OnOff On
-#set /Herwig/Particles/tbar/tbar->nu_taubar,tau-,bbar;:OnOff Off
-
set /Herwig/Particles/t/t->b,bbar,c;:OnOff Off
set /Herwig/Particles/t/t->b,c,dbar;:OnOff Off
set /Herwig/Particles/t/t->b,c,sbar;:OnOff Off
set /Herwig/Particles/t/t->b,sbar,u;:OnOff Off
set /Herwig/Particles/t/t->b,u,dbar;:OnOff Off
set /Herwig/Particles/t/t->nu_e,e+,b;:OnOff On
set /Herwig/Particles/t/t->nu_mu,mu+,b;:OnOff On
set /Herwig/Particles/t/t->nu_tau,tau+,b;:OnOff Off
set /Herwig/Particles/tbar/tbar->b,bbar,cbar;:OnOff On
set /Herwig/Particles/tbar/tbar->bbar,cbar,d;:OnOff On
set /Herwig/Particles/tbar/tbar->bbar,cbar,s;:OnOff On
set /Herwig/Particles/tbar/tbar->bbar,s,ubar;:OnOff On
set /Herwig/Particles/tbar/tbar->bbar,ubar,d;:OnOff On
set /Herwig/Particles/tbar/tbar->nu_ebar,e-,bbar;:OnOff Off
set /Herwig/Particles/tbar/tbar->nu_mubar,mu-,bbar;:OnOff Off
set /Herwig/Particles/tbar/tbar->nu_taubar,tau-,bbar;:OnOff Off
-
cd /Herwig/Generators
set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+#set /Herwig/UnderlyingEvent/KtCut:MinKT 2.55
+#set /Herwig/UnderlyingEvent/UECuts:MHatMin 5.1
+
set LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
##################################################
#
# This section contains the user defined settings
#
##################################################
# Example hard process: Incoming proton, outgoing leptoquarks
cd /Herwig/NewPhysics
set HPConstructor:Processes Exclusive
-
insert HPConstructor:Incoming 0 /Herwig/Particles/u
insert HPConstructor:Incoming 1 /Herwig/Particles/ubar
insert HPConstructor:Incoming 2 /Herwig/Particles/d
insert HPConstructor:Incoming 3 /Herwig/Particles/dbar
insert HPConstructor:Incoming 4 /Herwig/Particles/s
insert HPConstructor:Incoming 5 /Herwig/Particles/sbar
insert HPConstructor:Incoming 6 /Herwig/Particles/c
insert HPConstructor:Incoming 7 /Herwig/Particles/cbar
-insert HPConstructor:Incoming 8 /Herwig/Particles/g
+insert HPConstructor:Incoming 8 /Herwig/Particles/b
+insert HPConstructor:Incoming 9 /Herwig/Particles/bbar
+insert HPConstructor:Incoming 10 /Herwig/Particles/g
# Comment out the following line if Axigluon model is selected
insert /Herwig/NewPhysics/HPConstructor:Excluded 0 /Herwig/Particles/Ag
-
insert HPConstructor:Outgoing 0 /Herwig/Particles/t
insert HPConstructor:Outgoing 1 /Herwig/Particles/tbar
-
-# Intrinsic pT tune extrapolated to LHC energy
-set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
-
-
-cd /Herwig/Analysis
-set Basics:CheckQuark false
-
-
-
-# Cuts
-cd /Herwig/Cuts
-set JetKtCut:MinKT 0.0*GeV
-set QCDCuts:MHatMax 14000.0*GeV
-set EECuts:MHatMin 0.0*GeV
-set PhotonKtCut:MinKT 0.0*GeV
-set PhotonKtCut:MinEta -10.
-set PhotonKtCut:MaxEta 10.
-set MassCut:MinM 0.*GeV
+#This is available for comparison to the LO Hw++ result
+#cd /Herwig/MatrixElements
+#insert SimpleQCD:MatrixElements[0] MEHeavyQuark
cd /Herwig/Shower
set SplittingGenerator:ISR Yes
set SplittingGenerator:FSR Yes
cd /Herwig/EventHandlers
#set LHCHandler:CascadeHandler NULL
#set LHCHandler:HadronizationHandler NULL
#set LHCHandler:DecayHandler NULL
# Other parameters for run
cd /Herwig/Generators
-# Turn off MPI
+# Uncomment to turn off MPI
#set LHCGenerator:EventHandler:CascadeHandler:MPIHandler NULL
-
-
-
set LHCGenerator:DebugLevel 1
set LHCGenerator:NumberOfEvents 10000000
set LHCGenerator:RandomNumberGenerator:Seed 31122001
-set LHCGenerator:PrintEvent 1000
+set LHCGenerator:PrintEvent 100
set LHCGenerator:MaxErrors 10000
-
-insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
-set /Herwig/Analysis/HepMCFile:PrintEvent 10000
-set /Herwig/Analysis/HepMCFile:Format GenEvent
-set /Herwig/Analysis/HepMCFile:Units GeV_mm
+#insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+#set /Herwig/Analysis/HepMCFile:PrintEvent 10000
+#set /Herwig/Analysis/HepMCFile:Format GenEvent
+#set /Herwig/Analysis/HepMCFile:Units GeV_mm
saverun TVT-TTBA LHCGenerator

File Metadata

Mime Type
text/x-diff
Expires
Sat, May 3, 6:26 AM (1 d, 14 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4983022
Default Alt Text
(15 KB)

Event Timeline