No OneTemporary
Actions

Size

56 KB

Subscribers

None

View Options

	diff --git a/src/LH.model b/src/LH.model
	--- a/src/LH.model
	+++ b/src/LH.model
	@@ -1,161 +1,174 @@
	# -- ThePEG-repository --

	##################################################
	# Common setup for Little Higgs models
	#
	#
	# This file does not contain anything that
	# users need to touch. User settings are in
	# ???-LH.in
	#
	###################################################
	#
	# Create particle content
	#
	###################################################
	cd /Herwig/Particles

	create ThePEG::ParticleData A_H
	setup A_H 32 A_H 500.0 0. 0. 0.0 0 0 3 1
	create ThePEG::ParticleData Z_H
	setup Z_H 33 Z_H 900.0 0.0 0.0 0.0 0 0 3 1
	create ThePEG::ParticleData W_H+
	setup W_H+ 34 W_H+ 500.0 0. 0. 0.0 3 0 3 1
	create ThePEG::ParticleData W_H-
	setup W_H- -34 W_H- 500.0 0. 0. 0.0 -3 0 3 1
	makeanti W_H+ W_H-
	create ThePEG::ParticleData T
	setup T 8 T 1000.0 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData Tbar
	setup Tbar -8 Tbar 1000.0 0.0 0.0 0.0 -2 -3 2 1
	makeanti T Tbar
	create ThePEG::ParticleData Phi0
	setup Phi0 35 Phi0 1000.0 0. 0. 0.0 0 0 1 1
	create ThePEG::ParticleData PhiP
	setup PhiP 36 PhiP 1000.0 0. 0. 0.0 0 0 1 1
	create ThePEG::ParticleData Phi+
	setup Phi+ 37 Phi+ 1000.0 0. 0. 0.0 3 0 1 1
	create ThePEG::ParticleData Phi-
	setup Phi- -37 Phi- 1000.0 0. 0. 0.0 -3 0 1 1
	makeanti Phi+ Phi-
	create ThePEG::ParticleData Phi++
	setup Phi++ 38 Phi++ 1000.0 0. 0. 0.0 +6 0 1 1
	create ThePEG::ParticleData Phi--
	setup Phi-- -38 Phi-- 1000.0 0. 0. 0.0 -6 0 1 1
	makeanti Phi++ Phi--

	###################################################
	#
	# Main directory and model object
	#
	###################################################
	mkdir /Herwig/NewPhysics/LH
	cd /Herwig/NewPhysics/LH

	create Herwig::LHModel Model HwLHModel.so

	# SM couplings
	set Model:QCD/RunningAlphaS /Herwig/Couplings/NLOAlphaS
	set Model:EW/RunningAlphaEM /Herwig/Couplings/AlphaEM
	set Model:EW/CKM /Herwig/CKM
	set Model:RunningMass /Herwig/RunningMass
	###################################################
	#
	# Vertices
	#
	###################################################
	mkdir /Herwig/Vertices/LH
	cd /Herwig/Vertices/LH
	create Herwig::LHFFZVertex FFZVertex
	create Herwig::LHFFPVertex FFPVertex
	create Herwig::LHFFGVertex FFGVertex
	create Herwig::LHFFWVertex FFWVertex
	create Herwig::LHFFHVertex FFHVertex
	create Herwig::LHWWWVertex WWWVertex
	create Herwig::LHWWHVertex WWHVertex
	create Herwig::LHWWWWVertex WWWWVertex
	create Herwig::LHWHHVertex WHHVertex
	create Herwig::LHWWHHVertex WWHHVertex

	# LHWHHVertex.h LHWWHH



	cd /Herwig/NewPhysics/LH
	# SM vertices
	set Model:Vertex/GGG /Herwig/Vertices/GGGVertex
	set Model:Vertex/GGGG /Herwig/Vertices/GGGGVertex
	set Model:Vertex/HGG /Herwig/Vertices/HGGVertex
	set Model:Vertex/HPP /Herwig/Vertices/HPPVertex
	set Model:Vertex/HHH /Herwig/Vertices/HHHVertex
	# LH vertices
	set Model:Vertex/FFZ /Herwig/Vertices/LH/FFZVertex
	set Model:Vertex/FFP /Herwig/Vertices/LH/FFPVertex
	set Model:Vertex/FFG /Herwig/Vertices/LH/FFGVertex
	set Model:Vertex/FFW /Herwig/Vertices/LH/FFWVertex
	set Model:Vertex/FFH /Herwig/Vertices/LH/FFHVertex
	set Model:Vertex/WWW /Herwig/Vertices/LH/WWWVertex
	set Model:Vertex/WWH /Herwig/Vertices/LH/WWHVertex
	set Model:Vertex/WWWW /Herwig/Vertices/LH/WWWWVertex
	set Model:Vertex/WHH /Herwig/Vertices/LH/WHHVertex
	set Model:Vertex/WWHH /Herwig/Vertices/LH/WWHHVertex

	# Switch off the SM Higgs decays
	cd /Herwig/Particles

	set h0:Width_generator NULL
	set h0:Mass_generator NULL

	set h0/h0->b,bbar;:Active No
	set h0/h0->b,bbar;:BranchingRatio 0.0
	set h0/h0->W+,W-;:Active No
	set h0/h0->W+,W-;:BranchingRatio 0.0
	set h0/h0->tau-,tau+;:Active No
	set h0/h0->tau-,tau+;:BranchingRatio 0.0
	set h0/h0->g,g;:Active No
	set h0/h0->g,g;:BranchingRatio 0.0
	set h0/h0->c,cbar;:Active No
	set h0/h0->c,cbar;:BranchingRatio 0.0
	set h0/h0->Z0,Z0;:Active No
	set h0/h0->Z0,Z0;:BranchingRatio 0.0
	set h0/h0->gamma,gamma;:Active No
	set h0/h0->gamma,gamma;:BranchingRatio 0.0
	set h0/h0->Z0,gamma;:Active No
	set h0/h0->Z0,gamma;:BranchingRatio 0.0
	set h0/h0->mu-,mu+;:Active No
	set h0/h0->mu-,mu+;:BranchingRatio 0.0
	set h0/h0->t,tbar;:Active No
	set h0/h0->t,tbar;:BranchingRatio 0.0

	cd /Herwig/NewPhysics

	# Set the particles to be decayed
	insert NewModel:DecayParticles 0 /Herwig/Particles/A_H
	insert NewModel:DecayParticles 0 /Herwig/Particles/Z_H
	insert NewModel:DecayParticles 0 /Herwig/Particles/W_H+
	insert NewModel:DecayParticles 0 /Herwig/Particles/T
	insert NewModel:DecayParticles 0 /Herwig/Particles/Phi0
	insert NewModel:DecayParticles 0 /Herwig/Particles/PhiP
	insert NewModel:DecayParticles 0 /Herwig/Particles/Phi+
	insert NewModel:DecayParticles 0 /Herwig/Particles/Phi++
	insert NewModel:DecayParticles 0 /Herwig/Particles/h0

	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 32
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 33
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 34
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 8
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 35
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 36
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 37
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 38

	# Set up the model framework
	set LH/Model:ModelGenerator NewModel

	# Shower stuff for coloured particles
	cd /Herwig/Shower
	+
	+create Herwig::ZeroZeroOneSplitFn PhitoPhiGammaSplitFn
	+set PhitoPhiGammaSplitFn:InteractionType QED
	+set PhitoPhiGammaSplitFn:ColourStructure ChargedChargedNeutral
	+set PhitoPhiGammaSplitFn:AngularOrdered Yes
	+
	+cp SudakovCommon PhitoPhiGammaSudakov
	+set PhitoPhiGammaSudakov:SplittingFunction PhitoPhiGammaSplitFn
	+
	do SplittingGenerator:AddFinalSplitting T->T,g; QtoQGSudakov

	+do SplittingGenerator:AddFinalSplitting T->T,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting W_H+->W_H+,gamma; WtoWGammaSudakov
	+do SplittingGenerator:AddFinalSplitting Phi+->Phi+,gamma; PhitoPhiGammaSudakov
	+do SplittingGenerator:AddFinalSplitting Phi++->Phi++,gamma; PhitoPhiGammaSudakov
	cd /Herwig/Generators
	# Choose Little Higgs over SM
	set EventGenerator:StandardModelParameters /Herwig/NewPhysics/LH/Model
	diff --git a/src/LHTP.model b/src/LHTP.model
	--- a/src/LHTP.model
	+++ b/src/LHTP.model
	@@ -1,273 +1,298 @@
	# -- ThePEG-repository --

	##################################################
	# Common setup for Little Higgs models with T-parity
	#
	#
	# This file does not contain anything that
	# users need to touch. User settings are in
	# ???-LHTP.in
	#
	###################################################
	#
	# Create particle content
	#
	###################################################
	cd /Herwig/Particles
	#
	# using PDG number for excited gauge bosons A_H=Z' Z_H=Z'' W_H=W'
	#
	create ThePEG::ParticleData A_H
	setup A_H 32 A_H 1000. 0.0 0.0 0.0 0 0 3 1
	create ThePEG::ParticleData Z_H
	setup Z_H 33 Z_H 1000. 0.0 0.0 0.0 0 0 3 1
	create ThePEG::ParticleData W_H+
	setup W_H+ 34 W_H+ 1000.0 0.0 0.0 0.0 3 0 3 1
	create ThePEG::ParticleData W_H-
	setup W_H- -34 W_H- 1000.0 0.0 0.0 0.0 -3 0 3 1
	makeanti W_H+ W_H-
	#
	# extra top quark using 4th generation number
	#
	create ThePEG::ParticleData TpT
	setup TpT 8 TpT 1000.0 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData TpTbar
	setup TpTbar -8 TpTbar 1000.0 0.0 0.0 0.0 -2 -3 2 1
	makeanti TpT TpTbar
	#
	# triplet Higgs bosons (Phi0,PhiP,Phi+/- using numbers
	# for doublet + Phi++/-- 38)
	create ThePEG::ParticleData Phi0
	setup Phi0 35 Phi0 1000.0 0. 0. 0.0 0 0 1 1
	create ThePEG::ParticleData PhiP
	setup PhiP 36 PhiP 1000.0 0. 0. 0.0 0 0 1 1
	create ThePEG::ParticleData Phi+
	setup Phi+ 37 Phi+ 1000.0 0. 0. 0.0 3 0 1 1
	create ThePEG::ParticleData Phi-
	setup Phi- -37 Phi- 1000.0 0. 0. 0.0 -3 0 1 1
	makeanti Phi+ Phi-
	create ThePEG::ParticleData Phi++
	setup Phi++ 38 Phi++ 1000.0 0. 0. 0.0 +6 0 1 1
	create ThePEG::ParticleData Phi--
	setup Phi-- -38 Phi-- 1000.0 0. 0. 0.0 -6 0 1 1
	makeanti Phi++ Phi--
	#
	# T-odd particles using numbers for excited particles
	#
	create ThePEG::ParticleData dmT
	setup dmT 4000001 dmT 1000.0 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData dmTbar
	setup dmTbar -4000001 dmTbar 1000.0 0.0 0.0 0.0 1 -3 2 1
	makeanti dmT dmTbar
	create ThePEG::ParticleData umT
	setup umT 4000002 umT 1000.0 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData umTbar
	setup umTbar -4000002 umTbar 1000.0 0.0 0.0 0.0 -2 -3 2 1
	makeanti umT umTbar
	create ThePEG::ParticleData smT
	setup smT 4000003 smT 1000.0 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData smTbar
	setup smTbar -4000003 smTbar 1000.0 0.0 0.0 0.0 1 -3 2 1
	makeanti smT smTbar
	create ThePEG::ParticleData cmT
	setup cmT 4000004 cmT 1000.0 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData cmTbar
	setup cmTbar -4000004 cmTbar 1000.0 0.0 0.0 0.0 -2 -3 2 1
	makeanti cmT cmTbar
	create ThePEG::ParticleData bmT
	setup bmT 4000005 bmT 1000.0 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData bmTbar
	setup bmTbar -4000005 bmTbar 1000.0 0.0 0.0 0.0 1 -3 2 1
	makeanti bmT bmTbar
	create ThePEG::ParticleData tmT
	setup tmT 4000006 tmT 1000.0 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData tmTbar
	setup tmTbar -4000006 tmTbar 1000.0 0.0 0.0 0.0 -2 -3 2 1
	makeanti tmT tmTbar
	create ThePEG::ParticleData TmT
	setup TmT 4000008 TmT 1000.0 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData TmTbar
	setup TmTbar -4000008 TmTbar 1000.0 0.0 0.0 0.0 -2 -3 2 1
	makeanti TmT TmTbar
	create ThePEG::ParticleData emT-
	setup emT- 4000011 emT- 1000. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData emT+
	setup emT+ -4000011 emT+ 1000. 0.0 0.0 0.0 3 0 2 1
	makeanti emT- emT+
	create ThePEG::ParticleData mumT-
	setup mumT- 4000013 mumT- 1000. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData mumT+
	setup mumT+ -4000013 mumT+ 1000. 0. 0.0 0.0 3 0 2 1
	makeanti mumT+ mumT-
	create ThePEG::ParticleData taumT+
	setup taumT+ -4000015 taumT+ 1000. 0.0 0.0 0.0 3 0 2 1
	create ThePEG::ParticleData taumT-
	setup taumT- 4000015 taumT- 1000. 0.0 0.0 0.0 -3 0 2 1
	makeanti taumT- taumT+
	create ThePEG::ParticleData nu_emT
	setup nu_emT 4000012 nu_emT 1000. 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData nu_emTbar
	setup nu_emTbar -4000012 nu_emTbar 1000 0.0 0.0 0.0 0 0 2 1
	makeanti nu_emTbar nu_emT
	create ThePEG::ParticleData nu_mumT
	setup nu_mumT 4000014 nu_mumT 1000. 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData nu_mumTbar
	setup nu_mumTbar -4000014 nu_mumTbar 0.0 0.0 0.0 0 0 0 2 1
	makeanti nu_mumTbar nu_mumT
	create ThePEG::ParticleData nu_taumT
	setup nu_taumT 4000016 nu_taumT 1000. 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData nu_taumTbar
	setup nu_taumTbar -4000016 nu_taumTbar 1000. 0.0 0.0 0.0 0 0 2 1
	makeanti nu_taumTbar nu_taumT
	###################################################
	#
	# Main directory and model object
	#
	###################################################

	mkdir /Herwig/NewPhysics/LHTP
	cd /Herwig/NewPhysics/LHTP

	create Herwig::LHTPModel Model HwLHTPModel.so

	# SM couplings
	set Model:QCD/RunningAlphaS /Herwig/Couplings/NLOAlphaS
	set Model:EW/RunningAlphaEM /Herwig/Couplings/AlphaEM
	set Model:EW/CKM /Herwig/CKM
	set Model:RunningMass /Herwig/RunningMass
	###################################################
	#
	# Vertices
	#
	###################################################
	mkdir /Herwig/Vertices/LHTP
	cd /Herwig/Vertices/LHTP
	create Herwig::LHTPFFGVertex FFGVertex
	create Herwig::LHTPFFPVertex FFPVertex
	create Herwig::LHTPFFWVertex FFWVertex
	create Herwig::LHTPFFZVertex FFZVertex
	create Herwig::LHTPWWHVertex WWHVertex
	create Herwig::LHTPWWWVertex WWWVertex
	create Herwig::LHTPHHHVertex HHHVertex
	create Herwig::LHTPWHHVertex WHHVertex
	create Herwig::LHTPFFHVertex FFHVertex

	cd /Herwig/NewPhysics/LHTP
	set Model:Vertex/FFG /Herwig/Vertices/LHTP/FFGVertex
	set Model:Vertex/FFZ /Herwig/Vertices/LHTP/FFZVertex
	set Model:Vertex/FFW /Herwig/Vertices/LHTP/FFWVertex
	set Model:Vertex/FFP /Herwig/Vertices/LHTP/FFPVertex
	set Model:Vertex/WWH /Herwig/Vertices/LHTP/WWHVertex
	set Model:Vertex/WWW /Herwig/Vertices/LHTP/WWWVertex
	set Model:Vertex/HHH /Herwig/Vertices/LHTP/HHHVertex
	set Model:Vertex/WHH /Herwig/Vertices/LHTP/WHHVertex
	set Model:Vertex/FFH /Herwig/Vertices/LHTP/FFHVertex

	# Not implemented LHTP yet
	set Model:Vertex/WWWW /Herwig/Vertices/WWWWVertex
	set Model:Vertex/HGG /Herwig/Vertices/HGGVertex
	set Model:Vertex/HPP /Herwig/Vertices/HPPVertex

	# No LHTP equivalent
	set Model:Vertex/GGG /Herwig/Vertices/GGGVertex
	set Model:Vertex/GGGG /Herwig/Vertices/GGGGVertex

	# Switch off the SM Higgs decays
	cd /Herwig/Particles

	set h0:Width_generator NULL
	set h0:Mass_generator NULL

	set h0/h0->b,bbar;:Active No
	set h0/h0->b,bbar;:BranchingRatio 0.0
	set h0/h0->W+,W-;:Active No
	set h0/h0->W+,W-;:BranchingRatio 0.0
	set h0/h0->tau-,tau+;:Active No
	set h0/h0->tau-,tau+;:BranchingRatio 0.0
	set h0/h0->g,g;:Active No
	set h0/h0->g,g;:BranchingRatio 0.0
	set h0/h0->c,cbar;:Active No
	set h0/h0->c,cbar;:BranchingRatio 0.0
	set h0/h0->Z0,Z0;:Active No
	set h0/h0->Z0,Z0;:BranchingRatio 0.0
	set h0/h0->gamma,gamma;:Active No
	set h0/h0->gamma,gamma;:BranchingRatio 0.0
	set h0/h0->Z0,gamma;:Active No
	set h0/h0->Z0,gamma;:BranchingRatio 0.0
	set h0/h0->mu-,mu+;:Active No
	set h0/h0->mu-,mu+;:BranchingRatio 0.0
	set h0/h0->t,tbar;:Active No
	set h0/h0->t,tbar;:BranchingRatio 0.0

	cd /Herwig/NewPhysics
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/HPPVertex
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/HGGVertex

	# Set up the model framework
	set LHTP/Model:ModelGenerator NewModel

	# --- Perturbative Decays ---
	insert NewModel:DecayParticles 0 /Herwig/Particles/h0
	insert NewModel:DecayParticles 0 /Herwig/Particles/A_H
	insert NewModel:DecayParticles 0 /Herwig/Particles/Z_H
	insert NewModel:DecayParticles 0 /Herwig/Particles/W_H+
	insert NewModel:DecayParticles 0 /Herwig/Particles/TpT
	insert NewModel:DecayParticles 0 /Herwig/Particles/Phi0
	insert NewModel:DecayParticles 0 /Herwig/Particles/PhiP
	insert NewModel:DecayParticles 0 /Herwig/Particles/Phi+
	insert NewModel:DecayParticles 0 /Herwig/Particles/Phi++
	insert NewModel:DecayParticles 0 /Herwig/Particles/dmT
	insert NewModel:DecayParticles 0 /Herwig/Particles/umT
	insert NewModel:DecayParticles 0 /Herwig/Particles/smT
	insert NewModel:DecayParticles 0 /Herwig/Particles/cmT
	insert NewModel:DecayParticles 0 /Herwig/Particles/bmT
	insert NewModel:DecayParticles 0 /Herwig/Particles/tmT
	insert NewModel:DecayParticles 0 /Herwig/Particles/TmT
	insert NewModel:DecayParticles 0 /Herwig/Particles/emT-
	insert NewModel:DecayParticles 0 /Herwig/Particles/mumT-
	insert NewModel:DecayParticles 0 /Herwig/Particles/taumT-
	insert NewModel:DecayParticles 0 /Herwig/Particles/nu_emT
	insert NewModel:DecayParticles 0 /Herwig/Particles/nu_mumT
	insert NewModel:DecayParticles 0 /Herwig/Particles/nu_taumT
	###################################################
	#
	# Shower stuff for coloured particles
	#
	###################################################
	cd /Herwig/Shower

	+create Herwig::ZeroZeroOneSplitFn PhitoPhiGammaSplitFn
	+set PhitoPhiGammaSplitFn:InteractionType QED
	+set PhitoPhiGammaSplitFn:ColourStructure ChargedChargedNeutral
	+set PhitoPhiGammaSplitFn:AngularOrdered Yes
	+
	+cp SudakovCommon PhitoPhiGammaSudakov
	+set PhitoPhiGammaSudakov:SplittingFunction PhitoPhiGammaSplitFn
	+
	do SplittingGenerator:AddFinalSplitting dmT->dmT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting umT->umT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting smT->smT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting cmT->cmT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting bmT->bmT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting tmT->tmT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting TmT->TmT,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting TpT->TpT,g; QtoQGSudakov

	+do SplittingGenerator:AddFinalSplitting dmT->dmT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting umT->umT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting smT->smT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting cmT->cmT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting bmT->bmT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting tmT->tmT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting TmT->TmT,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting TpT->TpT,gamma; QtoQGammaSudakov
	+
	+do SplittingGenerator:AddFinalSplitting emT-->emT-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting mumT-->mumT-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting taumT-->taumT-,gamma; LtoLGammaSudakov
	+
	+do SplittingGenerator:AddFinalSplitting W_H+->W_H+,gamma; WtoWGammaSudakov
	+do SplittingGenerator:AddFinalSplitting Phi+->Phi+,gamma; PhitoPhiGammaSudakov
	+do SplittingGenerator:AddFinalSplitting Phi++->Phi++,gamma; PhitoPhiGammaSudakov
	+
	cd /Herwig/Generators
	# Choose Little Higgs over SM
	set EventGenerator:StandardModelParameters /Herwig/NewPhysics/LHTP/Model

	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 32
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 33
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 34
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 8
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 35
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 36
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 37
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 38
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000001
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000002
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000003
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000004
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000005
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000006
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000008
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000011
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000012
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000013
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000014
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000015
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 4000016

	diff --git a/src/MSSM.model b/src/MSSM.model
	--- a/src/MSSM.model
	+++ b/src/MSSM.model
	@@ -1,472 +1,466 @@
	# -- ThePEG-repository --

	##################################################
	# Common setup for the MSSM
	#
	# See LHC-MSSM.in or ILC-MSSM.in for example usage
	#
	# This file does not contain anything that
	# users need to touch. User settings are in
	# ???-MSSM.in
	#
	###################################################
	#
	# Create particle content
	#
	###################################################
	cd /Herwig/Particles
	create ThePEG::ParticleData H0
	setup H0 35 H0 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData A0
	setup A0 36 A0 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData H+
	setup H+ 37 H+ 500.0 0.0 0.0 0.0 3 0 1 1
	create ThePEG::ParticleData H-
	setup H- -37 H- 500.0 0.0 0.0 0.0 -3 0 1 1
	makeanti H+ H-
	create ThePEG::ParticleData ~d_L
	setup ~d_L 1000001 ~d_L 500.0 0.0 0.0 0.0 -1 3 1 1
	create ThePEG::ParticleData ~d_Lbar
	setup ~d_Lbar -1000001 ~d_Lbar 500.0 0.0 0.0 0.0 1 -3 1 1
	makeanti ~d_L ~d_Lbar
	create ThePEG::ParticleData ~u_L
	setup ~u_L 1000002 ~u_L 500.0 0.0 0.0 0.0 2 3 1 0
	create ThePEG::ParticleData ~u_Lbar
	setup ~u_Lbar -1000002 ~u_Lbar 500.0 0.0 0.0 0.0 -2 -3 1 1
	makeanti ~u_L ~u_Lbar
	create ThePEG::ParticleData ~s_L
	setup ~s_L 1000003 ~s_L 500.0 0.0 0.0 0.0 -1 3 1 1
	create ThePEG::ParticleData ~s_Lbar
	setup ~s_Lbar -1000003 ~s_Lbar 500.0 0.0 0.0 0.0 1 -3 1 1
	makeanti ~s_L ~s_Lbar
	create ThePEG::ParticleData ~c_L
	setup ~c_L 1000004 ~c_L 500.0 0.0 0.0 0.0 2 3 1 1
	create ThePEG::ParticleData ~c_Lbar
	setup ~c_Lbar -1000004 ~c_Lbar 500.0 0.0 0.0 0.0 -2 -3 1 1
	makeanti ~c_L ~c_Lbar
	create ThePEG::ParticleData ~b_1
	setup ~b_1 1000005 ~b_1 500.0 0.0 0.0 0.0 -1 3 1 1
	create ThePEG::ParticleData ~b_1bar
	setup ~b_1bar -1000005 ~b_1bar 500.0 0.0 0.0 0.0 1 -3 1 1
	makeanti ~b_1 ~b_1bar
	create ThePEG::ParticleData ~t_1
	setup ~t_1 1000006 ~t_1 500.0 0.0 0.0 0.0 2 3 1 1
	create ThePEG::ParticleData ~t_1bar
	setup ~t_1bar -1000006 ~t_1bar 500.0 0.0 0.0 0.0 -2 -3 1 1
	makeanti ~t_1 ~t_1bar
	create ThePEG::ParticleData ~e_L-
	setup ~e_L- 1000011 ~e_L- 500.0 0.0 0.0 0.0 -3 0 1 1
	create ThePEG::ParticleData ~e_L+
	setup ~e_L+ -1000011 ~e_L+ 500.0 0.0 0.0 0.0 3 0 1 1
	makeanti ~e_L- ~e_L+
	create ThePEG::ParticleData ~nu_eL
	setup ~nu_eL 1000012 ~nu_eL 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData ~nu_eLbar
	setup ~nu_eLbar -1000012 ~nu_eLbar 500.0 0.0 0.0 0.0 0 0 1 1
	makeanti ~nu_eL ~nu_eLbar
	create ThePEG::ParticleData ~mu_L-
	setup ~mu_L- 1000013 ~mu_L- 500.0 0.0 0.0 0.0 -3 0 1 1
	create ThePEG::ParticleData ~mu_L+
	setup ~mu_L+ -1000013 ~mu_L+ 500.0 0.0 0.0 0.0 3 0 1 1
	makeanti ~mu_L- ~mu_L+
	create ThePEG::ParticleData ~nu_muL
	setup ~nu_muL 1000014 ~nu_muL 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData ~nu_muLbar
	setup ~nu_muLbar -1000014 ~nu_muLbar 500.0 0.0 0.0 0.0 0 0 1 1
	makeanti ~nu_muL ~nu_muLbar
	create ThePEG::ParticleData ~tau_1-
	setup ~tau_1- 1000015 ~tau_1- 500.0 0.0 0.0 0.0 -3 0 1 1
	create ThePEG::ParticleData ~tau_1+
	setup ~tau_1+ -1000015 ~tau_1+ 500.0 0.0 0.0 0.0 3 0 1 1
	makeanti ~tau_1- ~tau_1+
	create ThePEG::ParticleData ~nu_tauL
	setup ~nu_tauL 1000016 ~nu_tauL 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData ~nu_tauLbar
	setup ~nu_tauLbar -1000016 ~nu_tauLbar 500.0 0.0 0.0 0.0 0 0 1 1
	makeanti ~nu_tauL ~nu_tauLbar
	create ThePEG::ParticleData ~g
	setup ~g 1000021 ~g 500.0 0.0 0.0 0.0 0 8 2 1
	create ThePEG::ParticleData ~chi_10
	setup ~chi_10 1000022 ~chi_10 500.0 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData ~chi_20
	setup ~chi_20 1000023 ~chi_20 500.0 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData ~chi_1+
	setup ~chi_1+ 1000024 ~chi_1+ 500.0 0.0 0.0 0.0 3 0 2 1
	create ThePEG::ParticleData ~chi_1-
	setup ~chi_1- -1000024 ~chi_1- 500.0 0.0 0.0 0.0 -3 0 2 1
	makeanti ~chi_1+ ~chi_1-
	create ThePEG::ParticleData ~chi_30
	setup ~chi_30 1000025 ~chi_30 500.0 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData ~chi_40
	setup ~chi_40 1000035 ~chi_40 500.0 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData ~chi_2+
	setup ~chi_2+ 1000037 ~chi_2+ 500.0 0.0 0.0 0.0 3 0 2 1
	create ThePEG::ParticleData ~chi_2-
	setup ~chi_2- -1000037 ~chi_2- 500.0 0.0 0.0 0.0 -3 0 2 1
	makeanti ~chi_2+ ~chi_2-
	create ThePEG::ParticleData ~gravitino
	setup ~gravitino 1000039 ~gravitino 500.0 0.0 0.000010 0.0 0 0 4 1
	create ThePEG::ParticleData ~d_R
	setup ~d_R 2000001 ~d_R 500.0 0.0 0.0 0.0 -1 3 1 1
	create ThePEG::ParticleData ~d_Rbar
	setup ~d_Rbar -2000001 ~d_Rbar 500.0 0.0 0.0 0.0 1 -3 1 1
	makeanti ~d_R ~d_Rbar
	create ThePEG::ParticleData ~u_R
	setup ~u_R 2000002 ~u_R 500.0 0.0 0.0 0.0 2 3 1 1
	create ThePEG::ParticleData ~u_Rbar
	setup ~u_Rbar -2000002 ~u_Rbar 500.0 0.0 0.0 0.0 -2 -3 1 1
	makeanti ~u_R ~u_Rbar
	create ThePEG::ParticleData ~s_R
	setup ~s_R 2000003 ~s_R 500.0 0.0 0.0 0.0 -1 3 1 1
	create ThePEG::ParticleData ~s_Rbar
	setup ~s_Rbar -2000003 ~s_Rbar 500.0 0.0 0.0 0.0 1 -3 1 0
	makeanti ~s_R ~s_Rbar
	create ThePEG::ParticleData ~c_R
	setup ~c_R 2000004 ~c_R 500.0 0.0 0.0 0.0 2 3 1 1
	create ThePEG::ParticleData ~c_Rbar
	setup ~c_Rbar -2000004 ~c_Rbar 500.0 0.0 0.0 0.0 -2 -3 1 1
	makeanti ~c_R ~c_Rbar
	create ThePEG::ParticleData ~b_2
	setup ~b_2 2000005 ~b_2 500.0 0.0 0.0 0.0 -1 3 1 1
	create ThePEG::ParticleData ~b_2bar
	setup ~b_2bar -2000005 ~b_2bar 500.0 0.0 0.0 0.0 1 -3 1 1
	makeanti ~b_2 ~b_2bar
	create ThePEG::ParticleData ~t_2
	setup ~t_2 2000006 ~t_2 500.0 0.0 0.0 0.0 2 3 1 1
	create ThePEG::ParticleData ~t_2bar
	setup ~t_2bar -2000006 ~t_2bar 500.0 0.0 0.0 0.0 -2 -3 1 1
	makeanti ~t_2 ~t_2bar
	create ThePEG::ParticleData ~e_R-
	setup ~e_R- 2000011 ~e_R- 500.0 0.0 0.0 0.0 -3 0 1 1
	create ThePEG::ParticleData ~e_R+
	setup ~e_R+ -2000011 ~e_R+ 500.0 0.0 0.0 0.0 3 0 1 1
	makeanti ~e_R- ~e_R+
	create ThePEG::ParticleData ~mu_R-
	setup ~mu_R- 2000013 ~mu_R- 500.0 0.0 0.0 0.0 -3 0 1 1
	create ThePEG::ParticleData ~mu_R+
	setup ~mu_R+ -2000013 ~mu_R+ 500.0 0.0 0.0 0.0 3 0 1 1
	makeanti ~mu_R- ~mu_R+
	create ThePEG::ParticleData ~tau_2-
	setup ~tau_2- 2000015 ~tau_2- 500.0 0.0 0.0 0.0 -3 0 1 1
	create ThePEG::ParticleData ~tau_2+
	setup ~tau_2+ -2000015 ~tau_2+ 500.0 0.0 0.0 0.0 3 0 1 1
	makeanti ~tau_2- ~tau_2+
	###################################################
	#
	# Main directory and model object
	#
	###################################################
	mkdir /Herwig/NewPhysics/MSSM
	cd /Herwig/NewPhysics/MSSM
	create Herwig::MSSM Model HwSusy.so
	# SM couplings
	set Model:QCD/RunningAlphaS /Herwig/Couplings/NLOAlphaS
	set Model:EW/RunningAlphaEM /Herwig/Couplings/AlphaEM
	set Model:EW/CKM /Herwig/CKM
	set Model:RunningMass /Herwig/RunningMass
	###################################################
	#
	# Vertices
	#
	###################################################
	# create MSSM vertices
	mkdir /Herwig/Vertices/MSSM
	cd /Herwig/Vertices/MSSM
	# MSSM
	library HwSusy.so
	create Herwig::SSWSSVertex MSSM_WSS
	create Herwig::SSWWSSVertex MSSM_WWSS
	create Herwig::SSNFSVertex MSSM_NFS
	create Herwig::SSGFSVertex MSSM_GFS
	create Herwig::SSHSFSFVertex MSSM_HSS
	create Herwig::SSCFSVertex MSSM_CFS
	create Herwig::SSGSSVertex MSSM_GSS
	create Herwig::SSGNGVertex MSSM_GNG
	create Herwig::SSGGSQSQVertex MSSM_GGSS
	create Herwig::SSWGSSVertex MSSM_WGSS
	create Herwig::SSGSGSGVertex MSSM_GGOGO
	create Herwig::SSNNZVertex MSSM_NNZ
	create Herwig::SSNNPVertex MSSM_NNP
	create Herwig::SSCCZVertex MSSM_CCZ
	create Herwig::SSCNWVertex MSSM_CNW
	create Herwig::SSFFHVertex MSSM_FFH
	create Herwig::SSGOGOHVertex MSSM_GOGOH
	create Herwig::SSWWHVertex MSSM_WWH
	create Herwig::SSWHHVertex MSSM_WHH
	create Herwig::SSHHHVertex MSSM_HHH
	create Herwig::SSHGGVertex MSSM_HGG
	create Herwig::SSHPPVertex MSSM_HPP
	create Herwig::SSWWHHVertex MSSM_WWHH
	create Herwig::SSNCTVertex MSSM_NCT
	create Herwig::SSGVNHVertex MSSM_GVNH
	create Herwig::SSGVFSVertex MSSM_GVFS
	create Herwig::SSGVNVVertex MSSM_GVNV
	cd /Herwig/NewPhysics/MSSM
	# SM vertices
	set Model:Vertex/FFZ /Herwig/Vertices/FFZVertex
	set Model:Vertex/FFW /Herwig/Vertices/FFWVertex
	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/WWW /Herwig/Vertices/WWWVertex
	set Model:Vertex/WWWW /Herwig/Vertices/WWWWVertex
	# MSSM feynman rules
	set Model:Vertex/WSFSF /Herwig/Vertices/MSSM/MSSM_WSS
	set Model:Vertex/WWSFSF /Herwig/Vertices/MSSM/MSSM_WWSS
	set Model:Vertex/NFSF /Herwig/Vertices/MSSM/MSSM_NFS
	set Model:Vertex/GFSF /Herwig/Vertices/MSSM/MSSM_GFS
	set Model:Vertex/HSFSF /Herwig/Vertices/MSSM/MSSM_HSS
	set Model:Vertex/CFSF /Herwig/Vertices/MSSM/MSSM_CFS
	set Model:Vertex/GSFSF /Herwig/Vertices/MSSM/MSSM_GSS
	set Model:Vertex/GGSQSQ /Herwig/Vertices/MSSM/MSSM_GGSS
	set Model:Vertex/WGSQSQ /Herwig/Vertices/MSSM/MSSM_WGSS
	set Model:Vertex/GSGSG /Herwig/Vertices/MSSM/MSSM_GGOGO
	set Model:Vertex/GNG /Herwig/Vertices/MSSM/MSSM_GNG
	set Model:Vertex/NNZ /Herwig/Vertices/MSSM/MSSM_NNZ
	set Model:Vertex/NNP /Herwig/Vertices/MSSM/MSSM_NNP
	set Model:Vertex/CCZ /Herwig/Vertices/MSSM/MSSM_CCZ
	set Model:Vertex/CNW /Herwig/Vertices/MSSM/MSSM_CNW
	set Model:Vertex/FFH /Herwig/Vertices/MSSM/MSSM_FFH
	set Model:Vertex/GOGOH /Herwig/Vertices/MSSM/MSSM_GOGOH
	set Model:Vertex/WWH /Herwig/Vertices/MSSM/MSSM_WWH
	set Model:Vertex/SSWHH /Herwig/Vertices/MSSM/MSSM_WHH
	set Model:Vertex/HHH /Herwig/Vertices/MSSM/MSSM_HHH
	set Model:Vertex/HGG /Herwig/Vertices/MSSM/MSSM_HGG
	set Model:Vertex/HPP /Herwig/Vertices/MSSM/MSSM_HPP
	set Model:Vertex/WWHH /Herwig/Vertices/MSSM/MSSM_WWHH
	set Model:Vertex/NCT /Herwig/Vertices/MSSM/MSSM_NCT
	set Model:Vertex/GVNH /Herwig/Vertices/MSSM/MSSM_GVNH
	set Model:Vertex/GVFS /Herwig/Vertices/MSSM/MSSM_GVFS
	set Model:Vertex/GVNV /Herwig/Vertices/MSSM/MSSM_GVNV
	# Switch off the SM Higgs decays
	cd /Herwig/Particles

	set h0:Width_generator NULL
	set h0:Mass_generator NULL

	set h0/h0->b,bbar;:Active No
	set h0/h0->b,bbar;:BranchingRatio 0.0
	set h0/h0->W+,W-;:Active No
	set h0/h0->W+,W-;:BranchingRatio 0.0
	set h0/h0->tau-,tau+;:Active No
	set h0/h0->tau-,tau+;:BranchingRatio 0.0
	set h0/h0->g,g;:Active No
	set h0/h0->g,g;:BranchingRatio 0.0
	set h0/h0->c,cbar;:Active No
	set h0/h0->c,cbar;:BranchingRatio 0.0
	set h0/h0->Z0,Z0;:Active No
	set h0/h0->Z0,Z0;:BranchingRatio 0.0
	set h0/h0->gamma,gamma;:Active No
	set h0/h0->gamma,gamma;:BranchingRatio 0.0
	set h0/h0->Z0,gamma;:Active No
	set h0/h0->Z0,gamma;:BranchingRatio 0.0
	set h0/h0->mu-,mu+;:Active No
	set h0/h0->mu-,mu+;:BranchingRatio 0.0
	set h0/h0->t,tbar;:Active No
	set h0/h0->t,tbar;:BranchingRatio 0.0
	###################################################
	#
	# Set up spin correlation Decayers
	#
	###################################################
	cd /Herwig/NewPhysics

	set TwoBodyDC:CreateDecayModes No
	set ThreeBodyDC:CreateDecayModes No
	set WeakDecayConstructor:CreateDecayModes No
	set FourBodyDC:CreateDecayModes No
	# switch on decays via weak current
	insert DecayConstructor:NBodyDecayConstructors 0 WeakDecayConstructor
	# exclude gravitino and effective vertices from three body decays
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVNH
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVFS
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVNV
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_HGG
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_HPP
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_NNP
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GNG
	insert ThreeBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_NCT
	# and four body decays
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVNH
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVFS
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVNV
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_HGG
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_HPP
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_NNP
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GNG
	insert FourBodyDC:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_NCT
	# and hard processes
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVNH
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVFS
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GVNV
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_HGG
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_HPP
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_NNP
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_GNG
	insert HPConstructor:ExcludedVertices 0 /Herwig/Vertices/MSSM/MSSM_NCT
	# exclude intermediate higgs from 4 body decays
	insert FourBodyDC:ExcludedIntermediates 0 /Herwig/Particles/h0
	insert FourBodyDC:ExcludedIntermediates 0 /Herwig/Particles/H0
	insert FourBodyDC:ExcludedIntermediates 0 /Herwig/Particles/A0
	insert FourBodyDC:ExcludedIntermediates 0 /Herwig/Particles/H+

	# which particles get the off-shell treatment
	set NewModel:WhichOffshell All
	# particles for which decays are included
	insert NewModel:DecayParticles 0 /Herwig/Particles/~d_L
	insert NewModel:DecayParticles 0 /Herwig/Particles/~u_L
	insert NewModel:DecayParticles 0 /Herwig/Particles/~s_L
	insert NewModel:DecayParticles 0 /Herwig/Particles/~c_L
	insert NewModel:DecayParticles 0 /Herwig/Particles/~b_1
	insert NewModel:DecayParticles 0 /Herwig/Particles/~t_1
	insert NewModel:DecayParticles 0 /Herwig/Particles/~d_R
	insert NewModel:DecayParticles 0 /Herwig/Particles/~u_R
	insert NewModel:DecayParticles 0 /Herwig/Particles/~s_R
	insert NewModel:DecayParticles 0 /Herwig/Particles/~c_R
	insert NewModel:DecayParticles 0 /Herwig/Particles/~b_2
	insert NewModel:DecayParticles 0 /Herwig/Particles/~t_2
	insert NewModel:DecayParticles 0 /Herwig/Particles/~e_L-
	insert NewModel:DecayParticles 0 /Herwig/Particles/~nu_eL
	insert NewModel:DecayParticles 0 /Herwig/Particles/~mu_L-
	insert NewModel:DecayParticles 0 /Herwig/Particles/~nu_muL
	insert NewModel:DecayParticles 0 /Herwig/Particles/~tau_1-
	insert NewModel:DecayParticles 0 /Herwig/Particles/~nu_tauL
	insert NewModel:DecayParticles 0 /Herwig/Particles/~e_R-
	insert NewModel:DecayParticles 0 /Herwig/Particles/~mu_R-
	insert NewModel:DecayParticles 0 /Herwig/Particles/~tau_2-
	insert NewModel:DecayParticles 0 /Herwig/Particles/~g
	insert NewModel:DecayParticles 0 /Herwig/Particles/~chi_10
	insert NewModel:DecayParticles 0 /Herwig/Particles/~chi_20
	insert NewModel:DecayParticles 0 /Herwig/Particles/~chi_30
	insert NewModel:DecayParticles 0 /Herwig/Particles/~chi_40
	insert NewModel:DecayParticles 0 /Herwig/Particles/~chi_1+
	insert NewModel:DecayParticles 0 /Herwig/Particles/~chi_2+
	insert NewModel:DecayParticles 0 /Herwig/Particles/h0
	insert NewModel:DecayParticles 0 /Herwig/Particles/H0
	insert NewModel:DecayParticles 0 /Herwig/Particles/A0
	insert NewModel:DecayParticles 0 /Herwig/Particles/H+
	###################################################
	#
	# Shower stuff for coloured sparticles
	#
	###################################################
	cd /Herwig/Shower

	create Herwig::ZeroZeroOneSplitFn PhitoPhiGSplitFn
	set PhitoPhiGSplitFn:InteractionType QCD
	set PhitoPhiGSplitFn:ColourStructure TripletTripletOctet

	create Herwig::ZeroZeroOneSplitFn PhitoPhiGammaSplitFn
	set PhitoPhiGammaSplitFn:InteractionType QED
	set PhitoPhiGammaSplitFn:ColourStructure ChargedChargedNeutral
	set PhitoPhiGammaSplitFn:AngularOrdered Yes

	-create Herwig::HalfHalfOneSplitFn CharginotoCharginoGammaSplitFn
	-set CharginotoCharginoGammaSplitFn:InteractionType QED
	-set CharginotoCharginoGammaSplitFn:ColourStructure ChargedChargedNeutral
	-set CharginotoCharginoGammaSplitFn:AngularOrdered Yes

	create Herwig::HalfHalfOneSplitFn GluinotoGluinoGSplitFn
	set GluinotoGluinoGSplitFn:InteractionType QCD
	set GluinotoGluinoGSplitFn:ColourStructure OctetOctetOctet

	cp SudakovCommon PhitoPhiGSudakov
	set PhitoPhiGSudakov:SplittingFunction PhitoPhiGSplitFn

	cp SudakovCommon PhitoPhiGammaSudakov
	set PhitoPhiGammaSudakov:SplittingFunction PhitoPhiGammaSplitFn

	cp SudakovCommon GluinotoGluinoGSudakov
	set GluinotoGluinoGSudakov:SplittingFunction GluinotoGluinoGSplitFn

	-cp SudakovCommon CharginotoCharginoGammaSudakov
	-set CharginotoCharginoGammaSudakov:SplittingFunction CharginotoCharginoGammaSplitFn

	do SplittingGenerator:AddFinalSplitting ~u_L->~u_L,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~d_L->~d_L,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~s_L->~s_L,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~c_L->~c_L,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~b_1->~b_1,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~t_1->~t_1,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~u_R->~u_R,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~d_R->~d_R,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~s_R->~s_R,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~c_R->~c_R,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~b_2->~b_2,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~t_2->~t_2,g; PhitoPhiGSudakov
	do SplittingGenerator:AddFinalSplitting ~g->~g,g; GluinotoGluinoGSudakov

	do SplittingGenerator:AddFinalSplitting ~u_L->~u_L,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~d_L->~d_L,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~s_L->~s_L,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~c_L->~c_L,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~b_1->~b_1,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~t_1->~t_1,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~u_R->~u_R,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~d_R->~d_R,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~s_R->~s_R,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~c_R->~c_R,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~b_2->~b_2,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~t_2->~t_2,gamma; PhitoPhiGammaSudakov

	do SplittingGenerator:AddFinalSplitting ~e_L-->~e_L-,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~mu_L-->~mu_L-,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~tau_1-->~tau_1-,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~e_R-->~e_R-,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~mu_R-->~mu_R-,gamma; PhitoPhiGammaSudakov
	do SplittingGenerator:AddFinalSplitting ~tau_2-->~tau_2-,gamma; PhitoPhiGammaSudakov

	do SplittingGenerator:AddFinalSplitting H+->H+,gamma; PhitoPhiGammaSudakov
	-do SplittingGenerator:AddFinalSplitting ~chi_1+->~chi_1+,gamma; CharginotoCharginoGammaSudakov
	-do SplittingGenerator:AddFinalSplitting ~chi_2+->~chi_2+,gamma; CharginotoCharginoGammaSudakov
	+do SplittingGenerator:AddFinalSplitting ~chi_1+->~chi_1+,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting ~chi_2+->~chi_2+,gamma; LtoLGammaSudakov
	###################################################
	# Set up the bsm framework
	###################################################
	cd /Herwig/NewPhysics
	set MSSM/Model:ModelGenerator NewModel
	###################################################
	#
	# Choose MSSM over SM
	#
	###################################################
	cd /Herwig/Generators
	set EventGenerator:StandardModelParameters /Herwig/NewPhysics/MSSM/Model

	###################################################
	#
	# Get Shower to decay SUSY particles
	#
	###################################################
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000001 # SUSY_d_L
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000002 # SUSY_u_L
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000003 # SUSY_s_L
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000004 # SUSY_c_L
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000005 # SUSY_b_1
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000006 # SUSY_t_1
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000011 # SUSY_e_Lminus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000012 # SUSY_nu_eL
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000013 # SUSY_mu_Lminus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000014 # SUSY_nu_muL
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000015 # SUSY_tau_1minus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000016 # SUSY_nu_tauL
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000021 # SUSY_g
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000022 # SUSY_chi_10
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000023 # SUSY_chi_20
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000024 # SUSY_chi_1plus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000025 # SUSY_chi_30
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000035 # SUSY_chi_40
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000037 # SUSY_chi_2plus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 1000039 # SUSY_gravitino
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000001 # SUSY_d_R
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000002 # SUSY_u_R
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000003 # SUSY_s_R
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000004 # SUSY_c_R
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000005 # SUSY_b_2
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000006 # SUSY_t_2
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000011 # SUSY_e_Rminus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000013 # SUSY_mu_Rminus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 2000015 # SUSY_tau_2minus
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 35 # H0
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 36 # A0
	insert /Herwig/Shower/ShowerHandler:DecayInShower 0 37 # H+
	diff --git a/src/MUED.model b/src/MUED.model
	--- a/src/MUED.model
	+++ b/src/MUED.model
	@@ -1,335 +1,369 @@
	# -- ThePEG-repository --

	##################################################
	# Common setup for Universal extra dimensions
	#
	# See LHC-MUED.in or ILC-MUED.in for example usage
	#
	# This file does not contain anything that
	# users need to touch. User settings are in
	# ???-MUED.in
	#
	###################################################
	#
	# Create particle content
	#
	###################################################
	# The UED level-1 particles
	cd /Herwig/Particles
	create ThePEG::ParticleData KK1_d_L
	setup KK1_d_L 5100001 KK1_d_L 500. 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData KK1_d_Lbar
	setup KK1_d_Lbar -5100001 KK1_d_Lbar 500. 0.0 0.0 0.0 1 -3 2 1
	makeanti KK1_d_L KK1_d_Lbar
	create ThePEG::ParticleData KK1_u_L
	setup KK1_u_L 5100002 KK1_u_L 500. 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData KK1_u_Lbar
	setup KK1_u_Lbar -5100002 KK1_u_Lbar 500. 0.0 0.0 0.0 -2 -3 2 1
	makeanti KK1_u_L KK1_u_Lbar
	create ThePEG::ParticleData KK1_s_L
	setup KK1_s_L 5100003 KK1_s_L 500. 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData KK1_s_Lbar
	setup KK1_s_Lbar -5100003 KK1_s_Lbar 500. 0.0 0.0 0.0 1 -3 2 1
	makeanti KK1_s_L KK1_s_Lbar
	create ThePEG::ParticleData KK1_c_L
	setup KK1_c_L 5100004 KK1_c_L 500. 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData KK1_c_Lbar
	setup KK1_c_Lbar -5100004 KK1_c_Lbar 500. 0.0 0.0 0.0 -2 -3 2 1
	makeanti KK1_c_L KK1_c_Lbar
	create ThePEG::ParticleData KK1_b_1
	setup KK1_b_1 5100005 KK1_b_1 500. 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData KK1_b_1bar
	setup KK1_b_1bar -5100005 KK1_b_1bar 500. 0.0 0.0 0.0 1 -3 2 1
	makeanti KK1_b_1 KK1_b_1bar
	create ThePEG::ParticleData KK1_t_1
	setup KK1_t_1 5100006 KK1_t_1 500. 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData KK1_t_1bar
	setup KK1_t_1bar -5100006 KK1_t_1bar 500. 0.0 0.0 0.0 -2 -3 2 1
	makeanti KK1_t_1 KK1_t_1bar
	# Right
	create ThePEG::ParticleData KK1_d_R
	setup KK1_d_R 6100001 KK1_d_R 500. 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData KK1_d_Rbar
	setup KK1_d_Rbar -6100001 KK1_d_Rbar 500. 0.0 0.0 0.0 1 -3 2 1
	makeanti KK1_d_R KK1_d_Rbar
	create ThePEG::ParticleData KK1_u_R
	setup KK1_u_R 6100002 KK1_u_R 500. 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData KK1_u_Rbar
	setup KK1_u_Rbar -6100002 KK1_u_Rbar 500. 0.0 0.0 0.0 -2 -3 2 1
	makeanti KK1_u_R KK1_u_Rbar
	create ThePEG::ParticleData KK1_s_R
	setup KK1_s_R 6100003 KK1_s_R 500. 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData KK1_s_Rbar
	setup KK1_s_Rbar -6100003 KK1_s_Rbar 500. 0.0 0.0 0.0 1 -3 2 1
	makeanti KK1_s_R KK1_s_Rbar
	create ThePEG::ParticleData KK1_c_R
	setup KK1_c_R 6100004 KK1_c_R 500. 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData KK1_c_Rbar
	setup KK1_c_Rbar -6100004 KK1_c_Rbar 500. 0.0 0.0 0.0 -2 -3 2 1
	makeanti KK1_c_R KK1_c_Rbar
	create ThePEG::ParticleData KK1_b_2
	setup KK1_b_2 6100005 KK1_b_2 500. 0.0 0.0 0.0 -1 3 2 1
	create ThePEG::ParticleData KK1_b_2bar
	setup KK1_b_2bar -6100005 KK1_b_2bar 500. 0.0 0.0 0.0 1 -3 2 1
	makeanti KK1_b_2 KK1_b_2bar
	create ThePEG::ParticleData KK1_t_2
	setup KK1_t_2 6100006 KK1_t_2 500. 0.0 0.0 0.0 2 3 2 1
	create ThePEG::ParticleData KK1_t_2bar
	setup KK1_t_2bar -6100006 KK1_t_2bar 500. 0.0 0.0 0.0 -2 -3 2 1
	makeanti KK1_t_2 KK1_t_2bar
	create ThePEG::ParticleData KK1_e_L-
	setup KK1_e_L- 5100011 KK1_e_L- 500. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData KK1_e_L+
	setup KK1_e_L+ -5100011 KK1_e_L+ 500. 0.0 0.0 0.0 3 0 2 1
	makeanti KK1_e_L- KK1_e_L+
	create ThePEG::ParticleData KK1_nu_eL
	setup KK1_nu_eL 5100012 KK1_nu_eL 500. 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData KK1_nu_eLbar
	setup KK1_nu_eLbar -5100012 KK1_nu_eLbar 500. 0.0 0.0 0.0 0 0 2 1
	makeanti KK1_nu_eL KK1_nu_eLbar
	create ThePEG::ParticleData KK1_mu_L-
	setup KK1_mu_L- 5100013 KK1_mu_L- 500. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData KK1_mu_L+
	setup KK1_mu_L+ -5100013 KK1_mu_L+ 500. 0.0 0.0 0.0 3 0 2 1
	makeanti KK1_mu_L- KK1_mu_L+
	create ThePEG::ParticleData KK1_nu_muL
	setup KK1_nu_muL 5100014 KK1_nu_muL 500. 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData KK1_nu_muLbar
	setup KK1_nu_muLbar -5100014 KK1_nu_muLbar 500. 0.0 0.0 0.0 0 0 2 1
	makeanti KK1_nu_muL KK1_nu_muLbar
	create ThePEG::ParticleData KK1_tau_1-
	setup KK1_tau_1- 5100015 KK1_tau_1- 500. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData KK1_tau_1+
	setup KK1_tau_1+ -5100015 KK1_tau_1+ 500. 0.0 0.0 0.0 3 0 2 1
	makeanti KK1_tau_1- KK1_tau_1+
	create ThePEG::ParticleData KK1_nu_tauL
	setup KK1_nu_tauL 5100016 KK1_nu_tauL 500. 0.0 0.0 0.0 0 0 2 1
	create ThePEG::ParticleData KK1_nu_tauLbar
	setup KK1_nu_tauLbar -5100016 KK1_nu_tauLbar 500. 0.0 0.0 0.0 0 0 2 1
	makeanti KK1_nu_tauL KK1_nu_tauLbar
	create ThePEG::ParticleData KK1_e_R-
	setup KK1_e_R- 6100011 KK1_e_R- 500. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData KK1_e_R+
	setup KK1_e_R+ -6100011 KK1_e_R+ 500. 0.0 0.0 0.0 3 0 2 1
	makeanti KK1_e_R- KK1_e_R+
	create ThePEG::ParticleData KK1_mu_R-
	setup KK1_mu_R- 6100013 KK1_mu_R- 500. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData KK1_mu_R+
	setup KK1_mu_R+ -6100013 KK1_mu_R+ 500. 0.0 0.0 0.0 3 0 2 1
	makeanti KK1_mu_R- KK1_mu_R+
	create ThePEG::ParticleData KK1_tau_2-
	setup KK1_tau_2- 6100015 KK1_tau_2- 500. 0.0 0.0 0.0 -3 0 2 1
	create ThePEG::ParticleData KK1_tau_2+
	setup KK1_tau_2+ -6100015 KK1_tau_2+ 500. 0.0 0.0 0.0 3 0 2 1
	makeanti KK1_tau_2- KK1_tau_2+
	create ThePEG::ParticleData KK1_g
	setup KK1_g 5100021 KK1_g 500. 0.0 0.0 0.0 0 8 3 1
	create ThePEG::ParticleData KK1_gamma
	setup KK1_gamma 5100022 KK1_gamma 500. 0.0 0.0 0.0 0 0 3 1
	create ThePEG::ParticleData KK1_Z0
	setup KK1_Z0 5100023 KK1_Z0 500. 0.0 0.0 0.0 0 0 3 1
	create ThePEG::ParticleData KK1_W+
	setup KK1_W+ 5100024 KK1_W+ 500. 0.0 0.0 0.0 3 0 3 1
	create ThePEG::ParticleData KK1_W-
	setup KK1_W- -5100024 KK1_W- 500. 0.0 0.0 0.0 -3 0 3 1
	makeanti KK1_W+ KK1_W-
	create ThePEG::ParticleData KK1_h0
	setup KK1_h0 5100025 KK1_h0 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData KK1_A0
	setup KK1_A0 5100036 KK1_A0 500.0 0.0 0.0 0.0 0 0 1 1
	create ThePEG::ParticleData KK1_H+
	setup KK1_H+ 5100037 KK1_H+ 500.0 0.0 0.0 0.0 3 0 1 1
	create ThePEG::ParticleData KK1_H-
	setup KK1_H- -5100037 KK1_H- 500.0 0.0 0.0 0.0 -3 0 1 1
	makeanti KK1_H+ KK1_H-
	###################################################
	#
	# Main directory and model object
	#
	###################################################
	mkdir /Herwig/NewPhysics/MUED
	cd /Herwig/NewPhysics/MUED
	create Herwig::UEDBase Model HwUED.so
	# SM couplings
	set Model:QCD/RunningAlphaS /Herwig/Couplings/NLOAlphaS
	set Model:EW/RunningAlphaEM /Herwig/Couplings/AlphaEM
	set Model:EW/CKM /Herwig/CKM
	set Model:RunningMass /Herwig/RunningMass
	###################################################
	#
	# Vertices
	#
	###################################################
	mkdir /Herwig/Vertices/MUED
	cd /Herwig/Vertices/MUED
	# create MUED vertices
	library HwUED.so
	create Herwig::UEDF1F1Z0Vertex UED_F1F1Z
	create Herwig::UEDF1F1P0Vertex UED_F1F1P
	create Herwig::UEDF1F1W0Vertex UED_F1F1W
	create Herwig::UEDF1F1G0Vertex UED_F1F1G0
	create Herwig::UEDF1F0G1Vertex UED_F1F0G1
	create Herwig::UEDF1F0W1Vertex UED_F1F0W1
	create Herwig::UEDF1F0H1Vertex UED_F1F0H1
	create Herwig::UEDG1G1G0Vertex UED_G1G1G0
	create Herwig::UEDG0G0G1G1Vertex UED_G0G0G1G1
	create Herwig::UEDP0H1H1Vertex UED_P0H1H1
	create Herwig::UEDZ0H1H1Vertex UED_Z0H1H1
	create Herwig::UEDW0A1H1Vertex UED_W0A1H1
	create Herwig::UEDZ0A1h1Vertex UED_Z0A1h1
	create Herwig::UEDW0W1W1Vertex UED_W0W1W1
	# SM vertices
	cd /Herwig/NewPhysics/MUED
	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/HPP /Herwig/Vertices/HPPVertex
	# MUED Feynman rules
	set Model:Vertex/F1F1Z /Herwig/Vertices/MUED/UED_F1F1Z
	set Model:Vertex/F1F1P /Herwig/Vertices/MUED/UED_F1F1P
	set Model:Vertex/F1F1W /Herwig/Vertices/MUED/UED_F1F1W
	set Model:Vertex/F1F1G0 /Herwig/Vertices/MUED/UED_F1F1G0
	set Model:Vertex/F1F0G1 /Herwig/Vertices/MUED/UED_F1F0G1
	set Model:Vertex/F1F0W1 /Herwig/Vertices/MUED/UED_F1F0W1
	set Model:Vertex/F1F0H1 /Herwig/Vertices/MUED/UED_F1F0H1
	set Model:Vertex/G1G1G0 /Herwig/Vertices/MUED/UED_G1G1G0
	set Model:Vertex/G0G0G1G1 /Herwig/Vertices/MUED/UED_G0G0G1G1
	set Model:Vertex/P0H1H1 /Herwig/Vertices/MUED/UED_P0H1H1
	set Model:Vertex/Z0H1H1 /Herwig/Vertices/MUED/UED_Z0H1H1
	set Model:Vertex/W0A1H1 /Herwig/Vertices/MUED/UED_W0A1H1
	set Model:Vertex/Z0A1h1 /Herwig/Vertices/MUED/UED_Z0A1h1
	set Model:Vertex/W0W1W1 /Herwig/Vertices/MUED/UED_W0W1W1
	###################################################
	#
	# Set up spin correlation Decayers and remove SM
	# Higgs related things
	#
	###################################################
	cd /Herwig/Particles

	set h0:Width_generator NULL
	set h0:Mass_generator NULL

	set h0/h0->b,bbar;:Active No
	set h0/h0->b,bbar;:BranchingRatio 0.0
	set h0/h0->W+,W-;:Active No
	set h0/h0->W+,W-;:BranchingRatio 0.0
	set h0/h0->tau-,tau+;:Active No
	set h0/h0->tau-,tau+;:BranchingRatio 0.0
	set h0/h0->g,g;:Active No
	set h0/h0->g,g;:BranchingRatio 0.0
	set h0/h0->c,cbar;:Active No
	set h0/h0->c,cbar;:BranchingRatio 0.0
	set h0/h0->Z0,Z0;:Active No
	set h0/h0->Z0,Z0;:BranchingRatio 0.0
	set h0/h0->gamma,gamma;:Active No
	set h0/h0->gamma,gamma;:BranchingRatio 0.0
	set h0/h0->Z0,gamma;:Active No
	set h0/h0->Z0,gamma;:BranchingRatio 0.0
	set h0/h0->mu-,mu+;:Active No
	set h0/h0->mu-,mu+;:BranchingRatio 0.0
	set h0/h0->t,tbar;:Active No
	set h0/h0->t,tbar;:BranchingRatio 0.0

	cd /Herwig/NewPhysics

	set TwoBodyDC:CreateDecayModes Yes
	set ThreeBodyDC:CreateDecayModes Yes

	# which particles get the off-shell treatment
	set NewModel:WhichOffshell All
	# particles for which decays are included
	insert NewModel:DecayParticles 0 /Herwig/Particles/KK1_d_L
	insert NewModel:DecayParticles 1 /Herwig/Particles/KK1_u_L
	insert NewModel:DecayParticles 2 /Herwig/Particles/KK1_s_L
	insert NewModel:DecayParticles 3 /Herwig/Particles/KK1_c_L
	insert NewModel:DecayParticles 4 /Herwig/Particles/KK1_b_1
	insert NewModel:DecayParticles 5 /Herwig/Particles/KK1_t_1
	insert NewModel:DecayParticles 6 /Herwig/Particles/KK1_d_R
	insert NewModel:DecayParticles 7 /Herwig/Particles/KK1_u_R
	insert NewModel:DecayParticles 8 /Herwig/Particles/KK1_s_R
	insert NewModel:DecayParticles 9 /Herwig/Particles/KK1_c_R
	insert NewModel:DecayParticles 10 /Herwig/Particles/KK1_b_2
	insert NewModel:DecayParticles 11 /Herwig/Particles/KK1_t_2
	insert NewModel:DecayParticles 12 /Herwig/Particles/KK1_e_L-
	insert NewModel:DecayParticles 13 /Herwig/Particles/KK1_nu_eL
	insert NewModel:DecayParticles 14 /Herwig/Particles/KK1_mu_L-
	insert NewModel:DecayParticles 15 /Herwig/Particles/KK1_nu_muL
	insert NewModel:DecayParticles 16 /Herwig/Particles/KK1_tau_1-
	insert NewModel:DecayParticles 17 /Herwig/Particles/KK1_nu_tauL
	insert NewModel:DecayParticles 18 /Herwig/Particles/KK1_e_R-
	insert NewModel:DecayParticles 19 /Herwig/Particles/KK1_mu_R-
	insert NewModel:DecayParticles 20 /Herwig/Particles/KK1_tau_2-
	insert NewModel:DecayParticles 21 /Herwig/Particles/KK1_g
	insert NewModel:DecayParticles 22 /Herwig/Particles/KK1_Z0
	insert NewModel:DecayParticles 23 /Herwig/Particles/KK1_W+
	insert NewModel:DecayParticles 24 /Herwig/Particles/h0
	insert NewModel:DecayParticles 25 /Herwig/Particles/KK1_h0
	insert NewModel:DecayParticles 26 /Herwig/Particles/KK1_A0
	insert NewModel:DecayParticles 27 /Herwig/Particles/KK1_H+
	###################################################
	#
	# Shower stuff for coloured sparticles
	#
	###################################################
	cd /Herwig/Shower
	insert ShowerHandler:DecayInShower 0 5100001
	insert ShowerHandler:DecayInShower 1 5100002
	insert ShowerHandler:DecayInShower 2 5100003
	insert ShowerHandler:DecayInShower 3 5100004
	insert ShowerHandler:DecayInShower 4 5100005
	insert ShowerHandler:DecayInShower 5 5100006
	insert ShowerHandler:DecayInShower 6 5100011
	insert ShowerHandler:DecayInShower 7 5100012
	insert ShowerHandler:DecayInShower 8 5100013
	insert ShowerHandler:DecayInShower 9 5100014
	insert ShowerHandler:DecayInShower 10 5100015
	insert ShowerHandler:DecayInShower 11 5100016
	insert ShowerHandler:DecayInShower 12 6100001
	insert ShowerHandler:DecayInShower 13 6100002
	insert ShowerHandler:DecayInShower 14 6100003
	insert ShowerHandler:DecayInShower 15 6100004
	insert ShowerHandler:DecayInShower 16 6100005
	insert ShowerHandler:DecayInShower 17 6100006
	insert ShowerHandler:DecayInShower 18 6100011
	insert ShowerHandler:DecayInShower 19 6100013
	insert ShowerHandler:DecayInShower 20 6100015
	insert ShowerHandler:DecayInShower 21 5100021
	insert ShowerHandler:DecayInShower 22 5100022
	insert ShowerHandler:DecayInShower 23 5100023
	insert ShowerHandler:DecayInShower 24 5100024
	insert ShowerHandler:DecayInShower 25 5100025
	insert ShowerHandler:DecayInShower 26 5100036
	insert ShowerHandler:DecayInShower 27 5100037
	+
	+create Herwig::ZeroZeroOneSplitFn PhitoPhiGammaSplitFn
	+set PhitoPhiGammaSplitFn:InteractionType QED
	+set PhitoPhiGammaSplitFn:ColourStructure ChargedChargedNeutral
	+set PhitoPhiGammaSplitFn:AngularOrdered Yes
	+
	+cp SudakovCommon PhitoPhiGammaSudakov
	+set PhitoPhiGammaSudakov:SplittingFunction PhitoPhiGammaSplitFn
	+
	do SplittingGenerator:AddFinalSplitting KK1_u_L->KK1_u_L,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_d_L->KK1_d_L,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_s_L->KK1_s_L,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_c_L->KK1_c_L,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_b_1->KK1_b_1,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_t_1->KK1_t_1,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_u_R->KK1_u_R,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_d_R->KK1_d_R,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_s_R->KK1_s_R,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_c_R->KK1_c_R,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_b_2->KK1_b_2,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting KK1_t_2->KK1_t_2,g; QtoQGSudakov
	+
	+do SplittingGenerator:AddFinalSplitting KK1_g->KK1_g,g; GtoGGSudakov
	+
	+do SplittingGenerator:AddFinalSplitting KK1_u_L->KK1_u_L,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_d_L->KK1_d_L,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_s_L->KK1_s_L,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_c_L->KK1_c_L,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_b_1->KK1_b_1,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_t_1->KK1_t_1,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_u_R->KK1_u_R,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_d_R->KK1_d_R,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_s_R->KK1_s_R,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_c_R->KK1_c_R,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_b_2->KK1_b_2,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_t_2->KK1_t_2,gamma; QtoQGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_e_L-->KK1_e_L-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_mu_L-->KK1_mu_L-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_tau_1-->KK1_tau_1-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_e_R-->KK1_e_R-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_mu_R-->KK1_mu_R-,gamma; LtoLGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_tau_2-->KK1_tau_2-,gamma; LtoLGammaSudakov
	+
	+do SplittingGenerator:AddFinalSplitting KK1_W+->KK1_W+,gamma; WtoWGammaSudakov
	+do SplittingGenerator:AddFinalSplitting KK1_H+->KK1_H+,gamma; PhitoPhiGammaSudakov
	+
	###################################################
	# Set up the model framework
	###################################################
	cd /Herwig/NewPhysics
	set MUED/Model:ModelGenerator NewModel
	###################################################
	#
	# Choose MSSM over SM
	#
	###################################################
	cd /Herwig/Generators
	set EventGenerator:StandardModelParameters /Herwig/NewPhysics/MUED/Model

File Metadata

Mime Type: text/x-diff
Expires: Sun, Feb 23, 2:47 PM (1 d, 7 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4486703
Default Alt Text: (56 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions