No OneTemporary
Actions

Size

50 KB

Subscribers

None

View Options

	diff --git a/Hadronization/ColourReconnector.h b/Hadronization/ColourReconnector.h
	--- a/Hadronization/ColourReconnector.h
	+++ b/Hadronization/ColourReconnector.h
	@@ -1,265 +1,271 @@
	// -- C++ --
	//
	// ColourReconnector.h is a part of Herwig - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2017 The Herwig Collaboration
	//
	// Herwig is licenced under version 3 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	#ifndef HERWIG_ColourReconnector_H
	#define HERWIG_ColourReconnector_H

	#include <ThePEG/Interface/Interfaced.h>
	#include "CluHadConfig.h"
	#include "ColourReconnector.fh"

	namespace Herwig {


	using namespace ThePEG;

	/** \ingroup Hadronization
	* \class ColourReconnector
	* \brief Class for changing colour reconnections of partons.
	* \author Alberto Ribon, Christian Roehr
	*
	* This class does the nonperturbative colour rearrangement, after the
	* nonperturbative gluon splitting and the "normal" cluster formation.
	* It uses the list of particles in the event record, and the collections of
	* "usual" clusters which is passed to the main method. If the colour
	* reconnection is actually accepted, then the previous collections of "usual"
	* clusters is first deleted and then the new one is created.
	*
	* * @see \ref ColourReconnectorInterfaces "The interfaces"
	* defined for ColourReconnector.
	*/
	class ColourReconnector: public Interfaced {

	public:

	/**
	* Does the colour rearrangement, starting out from the list of particles in
	* the event record and the collection of "usual" clusters passed as
	* arguments. If the actual rearrangement is accepted, the initial collection of
	* clusters is overridden by the old ones.
	*/
	void rearrange(ClusterVector & clusters);

	using CluVecIt = ClusterVector::iterator;

	private:

	/** PRIVATE MEMBER FUNCTIONS */

	/**
	* @brief Calculates the sum of the squared cluster masses.
	* @arguments q, aq vectors containing the quarks and antiquarks respectively
	* @return Sum of cluster squared masses M^2_{q[i],aq[i]}.
	*/
	Energy2 _clusterMassSum(const PVector & q, const PVector & aq) const;


	/**
	* @brief Examines whether the cluster vector (under the given permutation of
	* the antiquarks) contains colour-octet clusters
	* @param cv Cluster vector
	* @param P Permutation, a vector of permutated indices from 0 to
	* cv.size()-1
	*/
	bool _containsColour8(const ClusterVector & cv, const vector<size_t> & P) const;

	/**
	* @brief A Metropolis-type algorithm which finds a local minimum in the
	* total sum of cluster masses
	* @arguments cv cluster vector
	*/
	void _doRecoStatistical(ClusterVector & cv) const;

	/**
	* @brief Plain colour reconnection as used in Herwig 2.5.0
	* @arguments cv cluster vector
	*/
	void _doRecoPlain(ClusterVector & cv) const;


	/**
	* Baryonic Colour Reconnection model
	*/
	void _doRecoBaryonic(ClusterVector & cv) const;


	void _makeBaryonicClusters(ClusterPtr &c1, ClusterPtr &c2, ClusterPtr &c3,
	ClusterPtr &newcluster1, ClusterPtr &newcluster2) const;


	/**
	* @brief Finds the cluster in cv which, if reconnected with the given
	* cluster cl, would result in the smallest sum of cluster masses.
	* If no reconnection partner can be found, a pointer to the
	* original Cluster cl is returned.
	* @arguments cv cluster vector
	* cl cluster iterator (must be from cv) which wants to have a reconnection partner
	* @return iterator to the found cluster, or the original cluster pointer if
	* no mass-reducing combination can be found
	*/


	CluVecIt _findRecoPartner(CluVecIt cl, ClusterVector & cv) const;

	CluVecIt _findPartnerRapidity(CluVecIt cl, ClusterVector & cv) const;

	CluVecIt _findPartnerBaryonic(CluVecIt cl, ClusterVector & cv,
	bool & tetraCand,
	const ClusterVector& a,
	CluVecIt &baryonic1,
	CluVecIt &baryonic2 ) const;




	/**
	* @brief Reconnects the constituents of the given clusters to the (only)
	* other possible cluster combination.
	* @return pair of pointers to the two new clusters
	*/
	pair <ClusterPtr,ClusterPtr> _reconnect(ClusterPtr &c1, ClusterPtr &c2) const;

	/**
	* Reconnection method for baryonic reconenction model
	*/
	pair <ClusterPtr,ClusterPtr> _reconnectBaryonic(ClusterPtr &c1, ClusterPtr &c2) const;


	/**
	* @brief At random, swap two antiquarks, if not excluded by the
	* constraint that there must not be any colour-octet clusters.
	* @arguments q, aq vectors containing the quarks and antiquarks respectively
	* maxtries maximal number of tries to find a non-colour-octet
	* reconfiguration
	* @return Pair of ints indicating the indices of the antiquarks to be
	* swapped. Returns (-1,-1) if no valid reconfiguration could be
	* found after maxtries trials
	*/
	pair <int,int>
	_shuffle(const PVector & q, const PVector & aq, unsigned maxtries = 10) const;


	/** DATA MEMBERS */

	/**
	* Specifies the colour reconnection algorithm to be used.
	*/
	int _algorithm = 0;

	/**
	- * The annealing factor is the ratio of two successive temperature steps:
	- * T_n = _annealingFactor * T_(n-1)
	- */
	- double _annealingFactor = 0.9;
	-
	- /**
	- * Number of temperature steps in the statistical annealing algorithm
	- */
	- unsigned int _annealingSteps = 50;
	-
	- /**
	* Do we do colour reconnections?
	*/
	int _clreco = 0;

	+
	/**
	+ * Statistical Reco:
	* Factor used to determine the initial temperature according to
	* InitialTemperature = _initTemp * median {energy changes in a few random
	* rearrangements}
	*/
	- double _initTemp = 0.1;
	+ double _initTemp = 0.01;
	+
	+ /**
	+ * Statistical Reco:
	+ * The annealing factor is the ratio of two successive temperature steps:
	+ * T_n = _annealingFactor * T_(n-1)
	+ */
	+ double _annealingFactor = 0.21;
	+
	+ /**
	+ * Statistical Reco:
	+ * Number of temperature steps in the statistical annealing algorithm
	+ */
	+ unsigned int _annealingSteps = 10;
	+
	+ /**
	+ * Statistical Reco:
	+ * The number of tries per temperature steps is the number of clusters times
	+ * this factor.
	+ */
	+ double _triesPerStepFactor = 0.66;

	/**
	* Probability that a found reconnection possibility is actually accepted.
	*/
	double _preco = 0.5;


	double _precoBaryonic = 0.5;
	-
	- /**
	- * The number of tries per temperature steps is the number of clusters times
	- * this factor.
	- */
	- double _triesPerStepFactor = 5.0;
	+
	/**
	* maximum allowed distance in the eta phi space for reconnection to occur
	*/

	/**
	* Maximium distance for reconnections
	*/
	Length _maxDistance = picometer;

	/**
	* @return true, if the two partons are splitting products of the same
	* gluon
	*/
	bool _isColour8(tcPPtr p, tcPPtr q) const;

	/**
	* Option for handling octets
	*/
	unsigned int _octetOption = 0;

	public:

	/** @name Functions used by the persistent I/O system. */
	//@{
	/**
	* Function used to write out object persistently.
	* @param os the persistent output stream written to.
	*/
	void persistentOutput(PersistentOStream & os) const;

	/**
	* Function used to read in object persistently.
	* @param is the persistent input stream read from.
	* @param version the version number of the object when written.
	*/
	void persistentInput(PersistentIStream & is, int version);
	//@}

	/**
	* Standard Init function used to initialize the interfaces.
	*/
	static void Init();

	protected:

	/** @name Clone Methods. */
	//@{
	/**
	* Make a simple clone of this object.
	* @return a pointer to the new object.
	*/
	virtual IBPtr clone() const;

	/** Make a clone of this object, possibly modifying the cloned object
	* to make it sane.
	* @return a pointer to the new object.
	*/
	virtual IBPtr fullclone() const;
	//@}


	private:

	/**
	* Private and non-existent assignment operator.
	*/
	ColourReconnector & operator=(const ColourReconnector &) = delete;


	};


	}

	#endif /* HERWIG_ColourReconnector_H */

	diff --git a/src/LHC-MB.in b/src/LHC-MB.in
	--- a/src/LHC-MB.in
	+++ b/src/LHC-MB.in
	@@ -1,98 +1,62 @@
	# -- ThePEG-repository --

	################################################################################
	# This file contains our best tune to UE data from ATLAS at 7 TeV. More recent
	# tunes and tunes for other centre-of-mass energies as well as more usage
	# instructions can be obtained from this Herwig wiki page:
	# http://projects.hepforge.org/herwig/trac/wiki/MB_UE_tunes
	# The model for soft interactions and diffractions is explained in
	# [S. Gieseke, P. Kirchgaesser, F. Loshaj, arXiv:1612.04701]
	################################################################################

	read snippets/PPCollider.in

	##################################################
	# Technical parameters for this run
	##################################################
	cd /Herwig/Generators
	##################################################
	# LHC physics parameters (override defaults here)
	##################################################
	set EventGenerator:EventHandler:LuminosityFunction:Energy 7000.0

	# Intrinsic pT tune extrapolated to LHC energy
	set /Herwig/Shower/ShowerHandler:IntrinsicPtGaussian 2.2*GeV

	# Minimum Bias
	read snippets/MB.in

	-# Read in parameters of the soft model recommended for MB/UE simulations
	-read snippets/SoftTune.in
	-
	-
	#Diffraction model
	read snippets/Diffraction.in

	-#Turn on Baryonic reconnection
	-read snippets/BaryonicReconnection.in

	-
	-# Normalization of the Min bias cross section for correct diffractive cross section
	-set /Herwig/MatrixElements/MEMinBias:csNorm 0.01
	-
	-
	-# Use LHC parametrization of the cross section
	-set /Herwig/UnderlyingEvent/MPIHandler:DLmode 2
	-
	-
	-#some preliminary parameters for the MPI model which need to be tuned
	-#TODO
	-set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 3.02
	-set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.4
	-set /Herwig/UnderlyingEvent/MPIHandler:Power 0.308
	-set /Herwig/Partons/RemnantDecayer:ladderMult 0.45
	-set /Herwig/Partons/RemnantDecayer:ladderbFactor 1.0
	-set /Herwig/UnderlyingEvent/MPIHandler:DiffractiveRatio 0.17
	-
	-#ordered pTs
	-set /Herwig/Partons/RemnantDecayer:PtDistribution 4
	-
	-# set the correct PDFs
	-read Matchbox/CT14.in
	-
	-#set /Herwig/Shower/ShowerHandler:CascadeHandler NULL # Switch off parton shower
	-#set /Herwig/Hadronization/HadronizationHandler:HadronizationHandler NULL
	-# do Soft interactions
	-
	-#set /Herwig/Hadronization/ColourReconnector:ColourReconnection No
	#################################################
	# Analyses
	##################################################

	-#Comment these lines out in order to use rivet analyses
	-cd /Herwig/Analysis
	-create ThePEG::RivetAnalysis RivetAnalysis RivetAnalysis.so
	-
	-cd /Herwig/Generators
	-insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/RivetAnalysis
	-
	-#Some example analyses
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1084540
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562
	+## Comment these lines out in order to use rivet analyses
	+# cd /Herwig/Analysis
	+# create ThePEG::RivetAnalysis RivetAnalysis RivetAnalysis.so
	+
	+# cd /Herwig/Generators
	+# insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/RivetAnalysis
	+## Some example analyses
	+# insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1084540
	+# insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562

	#set /Herwig/Analysis/Plot:EventNumber 54
	#cd /Herwig/Generators
	#insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/Plot

	#insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
	#set /Herwig/Analysis/HepMCFile:PrintEvent 1000000
	#set /Herwig/Analysis/HepMCFile:Format GenEvent
	#set /Herwig/Analysis/HepMCFile:Units GeV_mm
	#set /Herwig/Analysis/HepMCFile:Filename events.fifo

	##################################################
	# Save run for later usage with 'Herwig run'
	##################################################
	cd /Herwig/Generators
	saverun LHC-MB EventGenerator

	diff --git a/src/defaults/Cuts.in b/src/defaults/Cuts.in
	--- a/src/defaults/Cuts.in
	+++ b/src/defaults/Cuts.in
	@@ -1,118 +1,118 @@
	# -- ThePEG-repository --

	###########################################################
	# Default cuts (applied to the hard subprocess)
	#
	# Don't change values here, re-set them in your own input
	# files using these as examples.
	###########################################################

	mkdir /Herwig/Matchers
	cd /Herwig/Matchers

	create ThePEG::Matcher<Lepton> Lepton
	create ThePEG::Matcher<Neutrino> Neutrino
	create ThePEG::Matcher<ChargedLepton> ChargedLepton
	create ThePEG::Matcher<LightQuark> LightQuark
	create ThePEG::Matcher<LightParticle> LightParticle
	create ThePEG::Matcher<LightAntiQuark> LightAntiQuark
	create ThePEG::Matcher<StandardQCDParton> StandardQCDParton
	create ThePEG::Matcher<Photon> Photon
	create ThePEG::Matcher<Top> Top
	create ThePEG::Matcher<Bottom> Bottom
	create ThePEG::Matcher<WBoson> WBoson
	create ThePEG::Matcher<ZBoson> ZBoson
	create ThePEG::Matcher<HiggsBoson> HiggsBoson

	mkdir /Herwig/Cuts
	cd /Herwig/Cuts


	# cut on jet pt
	create ThePEG::SimpleKTCut JetKtCut SimpleKTCut.so
	newdef JetKtCut:Matcher /Herwig/Matchers/StandardQCDParton
	newdef JetKtCut:MinKT 20.0*GeV

	# cut on photon
	create ThePEG::SimpleKTCut PhotonKtCut SimpleKTCut.so
	newdef PhotonKtCut:Matcher /Herwig/Matchers/Photon
	newdef PhotonKtCut:MinKT 20.0*GeV
	newdef PhotonKtCut:MinEta -3.
	newdef PhotonKtCut:MaxEta 3.

	# cut on leptons
	create ThePEG::SimpleKTCut LeptonKtCut SimpleKTCut.so
	newdef LeptonKtCut:Matcher /Herwig/Matchers/Lepton
	newdef LeptonKtCut:MinKT 0.0*GeV

	# cut on charged leptons
	create ThePEG::SimpleKTCut ChargedLeptonKtCut SimpleKTCut.so
	newdef ChargedLeptonKtCut:Matcher /Herwig/Matchers/ChargedLepton
	newdef ChargedLeptonKtCut:MinKT 0.0*GeV

	# cut on top quarks
	create ThePEG::KTRapidityCut TopKtCut KTRapidityCut.so
	newdef TopKtCut:Matcher /Herwig/Matchers/Top
	newdef TopKtCut:MinKT 0.0*GeV

	# cut on W bosons
	create ThePEG::KTRapidityCut WBosonKtCut KTRapidityCut.so
	newdef WBosonKtCut:Matcher /Herwig/Matchers/WBoson
	newdef WBosonKtCut:MinKT 0.0*GeV

	# cut on Z bosons
	create ThePEG::KTRapidityCut ZBosonKtCut KTRapidityCut.so
	newdef ZBosonKtCut:Matcher /Herwig/Matchers/ZBoson
	newdef ZBosonKtCut:MinKT 0.0*GeV

	# cut on Higgs bosons
	create ThePEG::KTRapidityCut HiggsBosonKtCut KTRapidityCut.so
	newdef HiggsBosonKtCut:Matcher /Herwig/Matchers/HiggsBoson
	newdef HiggsBosonKtCut:MinKT 0.0*GeV

	# cut on Higgs bosons
	create ThePEG::KTRapidityCut LightParticleKtCut KTRapidityCut.so
	newdef LightParticleKtCut:Matcher /Herwig/Matchers/LightParticle
	newdef LightParticleKtCut:MinKT 0.0*GeV

	# create a cut on the invariant mass of lepton pairs
	create ThePEG::V2LeptonsCut MassCut V2LeptonsCut.so
	newdef MassCut:Families All
	newdef MassCut:CComb All
	newdef MassCut:MinM 20.*GeV
	newdef MassCut:MaxM 14000.*GeV

	# create a cut on Q^2 for neutral current DIS
	create ThePEG::SimpleDISCut NeutralCurrentCut SimpleDISCut.so
	newdef NeutralCurrentCut:MinQ2 20.
	newdef NeutralCurrentCut:Current Neutral

	# create a cut on Q^2 for charged current DIS
	create ThePEG::SimpleDISCut ChargedCurrentCut SimpleDISCut.so
	newdef ChargedCurrentCut:MinQ2 20.
	newdef ChargedCurrentCut:Current Charged

	# create a cut of Q^2 for charged current DIS

	# Default Cuts object
	create ThePEG::Cuts Cuts
	newdef /Herwig/Generators/EventGenerator:EventHandler:Cuts /Herwig/Cuts/Cuts
	newdef Cuts:MHatMin 20.*GeV

	# insert into hadron cuts (by defaults)
	insert Cuts:OneCuts[0] JetKtCut
	insert Cuts:OneCuts[1] PhotonKtCut
	insert Cuts:OneCuts[2] LeptonKtCut
	insert Cuts:OneCuts[3] TopKtCut
	insert Cuts:OneCuts[4] WBosonKtCut
	insert Cuts:OneCuts[5] ZBosonKtCut
	insert Cuts:OneCuts[6] HiggsBosonKtCut
	insert Cuts:OneCuts[7] ChargedLeptonKtCut
	insert Cuts:MultiCuts[0] MassCut

	# create diffrent cuts object for MinBias to avoid numerical problems
	create ThePEG::Cuts MinBiasCuts
	newdef MinBiasCuts:ScaleMin 2.0*GeV2
	-newdef MinBiasCuts:X1Min 0.055
	-newdef MinBiasCuts:X2Min 0.055
	+newdef MinBiasCuts:X1Min 0.011
	+newdef MinBiasCuts:X2Min 0.011
	newdef MinBiasCuts:MHatMin 0.0*GeV
	diff --git a/src/defaults/Hadronization.in b/src/defaults/Hadronization.in
	--- a/src/defaults/Hadronization.in
	+++ b/src/defaults/Hadronization.in
	@@ -1,95 +1,96 @@
	# -- ThePEG-repository --

	############################################################
	# Setup of default hadronization
	#
	# There are no user servicable parts inside.
	#
	# Anything that follows below should only be touched if you
	# know what you're doing.
	#############################################################

	cd /Herwig/Particles
	create ThePEG::ParticleData Cluster
	setup Cluster 81 Cluster 0.00990 0.0 0.0 0.0 0 0 0 1
	create ThePEG::ParticleData Remnant
	setup Remnant 82 Remnant 0.00990 0.0 0.0 0.0 0 0 0 1

	mkdir /Herwig/Hadronization
	cd /Herwig/Hadronization

	create Herwig::ClusterHadronizationHandler ClusterHadHandler
	create Herwig::PartonSplitter PartonSplitter
	create Herwig::ClusterFinder ClusterFinder
	create Herwig::ColourReconnector ColourReconnector
	create Herwig::ClusterFissioner ClusterFissioner
	create Herwig::LightClusterDecayer LightClusterDecayer
	create Herwig::ClusterDecayer ClusterDecayer
	create Herwig::HwppSelector HadronSelector

	newdef ClusterHadHandler:PartonSplitter PartonSplitter
	newdef ClusterHadHandler:ClusterFinder ClusterFinder
	newdef ClusterHadHandler:ColourReconnector ColourReconnector
	newdef ClusterHadHandler:ClusterFissioner ClusterFissioner
	newdef ClusterHadHandler:LightClusterDecayer LightClusterDecayer
	newdef ClusterHadHandler:ClusterDecayer ClusterDecayer
	newdef ClusterHadHandler:MinVirtuality2 0.1*GeV2
	newdef ClusterHadHandler:MaxDisplacement 1.0e-10*millimeter
	newdef ClusterHadHandler:UnderlyingEventHandler NULL

	newdef ClusterFissioner:HadronSelector HadronSelector
	newdef LightClusterDecayer:HadronSelector HadronSelector
	newdef ClusterDecayer:HadronSelector HadronSelector

	newdef ColourReconnector:ColourReconnection Yes
	-newdef ColourReconnector:ReconnectionProbability 0.652710
	-newdef ColourReconnector:Algorithm Plain
	-newdef ColourReconnector:InitialTemperature 0.01
	-newdef ColourReconnector:AnnealingFactor 0.21
	-newdef ColourReconnector:AnnealingSteps 10
	-newdef ColourReconnector:TriesPerStepFactor 0.66
	+newdef ColourReconnector:ReconnectionProbabilityBaryonic 0.57
	+newdef ColourReconnector:ReconnectionProbability 0.89
	+newdef ColourReconnector:Algorithm Baryonic
	newdef ColourReconnector:OctetTreatment All

	# Clustering parameters for light quarks
	newdef ClusterFissioner:ClMaxLight 3.649
	newdef ClusterFissioner:ClPowLight 2.780
	newdef ClusterFissioner:PSplitLight 0.899
	newdef ClusterDecayer:ClDirLight 1
	newdef ClusterDecayer:ClSmrLight 0.78

	# Clustering parameters for b-quarks
	newdef ClusterFissioner:ClMaxBottom 3.757
	newdef ClusterFissioner:ClPowBottom 0.547
	newdef ClusterFissioner:PSplitBottom 0.625
	newdef ClusterDecayer:ClDirBottom 1
	newdef ClusterDecayer:ClSmrBottom 0.078
	newdef HadronSelector:SingleHadronLimitBottom 0.000

	# Clustering parameters for c-quarks
	newdef ClusterFissioner:ClMaxCharm 3.950
	newdef ClusterFissioner:ClPowCharm 2.559
	newdef ClusterFissioner:PSplitCharm 0.994
	newdef ClusterDecayer:ClDirCharm 1
	newdef ClusterDecayer:ClSmrCharm 0.163
	newdef HadronSelector:SingleHadronLimitCharm 0.000

	# Clustering parameters for exotic quarks
	# (e.g. hadronizing Susy particles)
	newdef ClusterFissioner:ClMaxExotic 2.7*GeV
	newdef ClusterFissioner:ClPowExotic 1.46
	newdef ClusterFissioner:PSplitExotic 1.00
	newdef ClusterDecayer:ClDirExotic 1
	newdef ClusterDecayer:ClSmrExotic 0.
	newdef HadronSelector:SingleHadronLimitExotic 0.

	#
	+newdef PartonSplitter:SplitPwtSquark 0.824135
	+newdef PartonSplitter:Split uds
	+
	+#
	newdef HadronSelector:PwtDquark 1.0
	newdef HadronSelector:PwtUquark 1.0
	-newdef HadronSelector:PwtSquark 0.700
	+newdef HadronSelector:PwtSquark 0.291717
	newdef HadronSelector:PwtCquark 1.0
	newdef HadronSelector:PwtBquark 1.0
	newdef HadronSelector:PwtDIquark 0.298
	newdef HadronSelector:SngWt 0.74
	newdef HadronSelector:DecWt 0.62
	newdef HadronSelector:Mode 1
	newdef HadronSelector:BelowThreshold All
	diff --git a/src/defaults/MatrixElements.in b/src/defaults/MatrixElements.in
	--- a/src/defaults/MatrixElements.in
	+++ b/src/defaults/MatrixElements.in
	@@ -1,249 +1,250 @@
	# -- ThePEG-repository --

	##############################################################################
	# Setup of default matrix elements.
	#
	# Only one ME is activated by default, but this file lists
	# some alternatives. All available MEs can be found in the
	# 'include/Herwig/MatrixElements' subdirectory of your Herwig
	# installation.
	#
	# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
	#
	# Instead of editing this file directly, you should reset
	# the matrix elements in your own input files:
	#
	# - create your custom SubProcessHandler
	# - insert the MEs you need
	# - set your SubProcessHandler instead of the default (see HerwigDefaults.in)
	##############################################################################
	mkdir /Herwig/MatrixElements
	cd /Herwig/MatrixElements
	library HwMELepton.so
	library HwMEHadron.so
	library HwMEDIS.so

	############################################################
	# e+e- matrix elements
	############################################################
	# e+e- > q qbar
	create Herwig::MEee2gZ2qq MEee2gZ2qq
	newdef MEee2gZ2qq:MinimumFlavour 1
	newdef MEee2gZ2qq:MaximumFlavour 5
	newdef MEee2gZ2qq:AlphaQCD /Herwig/Shower/AlphaQCD
	newdef MEee2gZ2qq:AlphaQED /Herwig/Shower/AlphaQED

	# e+e- -> l+l-
	create Herwig::MEee2gZ2ll MEee2gZ2ll
	newdef MEee2gZ2ll:Allowed Charged
	set MEee2gZ2ll:AlphaQED /Herwig/Shower/AlphaQED

	# e+e- -> W+W- ZZ
	create Herwig::MEee2VV MEee2VV

	# e+e- -> ZH
	create Herwig::MEee2ZH MEee2ZH
	newdef MEee2ZH:Coupling /Herwig/Shower/AlphaQCD

	# e+e- -> e+e-H/nu_enu_ebarH
	create Herwig::MEee2HiggsVBF MEee2HiggsVBF

	############################################################
	# NLO (POWHEG e+e- matrix elements
	############################################################
	library HwPowhegMELepton.so
	create Herwig::MEee2gZ2qqPowheg PowhegMEee2gZ2qq
	newdef PowhegMEee2gZ2qq:MinimumFlavour 1
	newdef PowhegMEee2gZ2qq:MaximumFlavour 5
	newdef PowhegMEee2gZ2qq:AlphaQCD /Herwig/Shower/AlphaQCD
	newdef PowhegMEee2gZ2qq:AlphaQED /Herwig/Shower/AlphaQED
	create Herwig::MEee2gZ2llPowheg PowhegMEee2gZ2ll
	newdef PowhegMEee2gZ2ll:Allowed Charged
	set PowhegMEee2gZ2ll:AlphaQED /Herwig/Shower/AlphaQED

	############################################################
	# hadron-hadron matrix elements
	############################################################

	###################################
	# Electroweak processes
	###################################
	# q qbar -> gamma/Z -> l+l-
	create Herwig::MEqq2gZ2ff MEqq2gZ2ff
	newdef MEqq2gZ2ff:Process 3
	newdef MEqq2gZ2ff:Coupling /Herwig/Shower/AlphaQCD

	# q qbar to W -> l nu
	create Herwig::MEqq2W2ff MEqq2W2ff
	newdef MEqq2W2ff:Process 2
	newdef MEqq2W2ff:Coupling /Herwig/Shower/AlphaQCD

	# W+jet
	create Herwig::MEPP2WJet MEWJet
	newdef MEWJet:WDecay Leptons

	# Z+jet
	create Herwig::MEPP2ZJet MEZJet
	newdef MEZJet:ZDecay ChargedLeptons

	# PP->WW/WZ/ZZ
	create Herwig::MEPP2VV MEPP2VV

	# PP->WZ gamma
	create Herwig::MEPP2VGamma MEPP2VGamma

	###################################
	# Photon and jet processes
	###################################

	# qqbar/gg -> gamma gamma
	create Herwig::MEPP2GammaGamma MEGammaGamma

	# hadron-hadron to gamma+jet
	create Herwig::MEPP2GammaJet MEGammaJet

	# QCD 2-to-2
	create Herwig::MEQCD2to2 MEQCD2to2

	# MinBias
	create Herwig::MEMinBias MEMinBias
	-
	+newdef MEMinBias:csNorm 0.01
	+newdef MEMinBias:Scale 2.0
	###################################
	# Heavy Quark
	###################################

	# qqbar/gg -> t tbar
	create Herwig::MEPP2QQ MEHeavyQuark

	create Herwig::MEPP2SingleTop MESingleTopTChannel
	set MESingleTopTChannel:Process tChannel

	create Herwig::MEPP2SingleTop MESingleTopSChannel
	set MESingleTopSChannel:Process sChannel

	create Herwig::MEPP2SingleTop MESingleTopTW
	set MESingleTopTW:Process tW

	###################################
	# Higgs processes
	###################################

	# hadron-hadron to higgs
	create Herwig::MEPP2Higgs MEHiggs
	newdef MEHiggs:ShapeScheme MassGenerator
	newdef MEHiggs:Process gg
	newdef MEHiggs:Coupling /Herwig/Shower/AlphaQCD

	# hadron-hadron to higgs+jet
	create Herwig::MEPP2HiggsJet MEHiggsJet

	# PP->ZH
	create Herwig::MEPP2ZH MEPP2ZH
	newdef MEPP2ZH:Coupling /Herwig/Shower/AlphaQCD

	# PP->WH
	create Herwig::MEPP2WH MEPP2WH
	newdef MEPP2WH:Coupling /Herwig/Shower/AlphaQCD

	# PP -> Higgs via VBF
	create Herwig::MEPP2HiggsVBF MEPP2HiggsVBF
	newdef MEPP2HiggsVBF:ShowerAlphaQCD /Herwig/Shower/AlphaQCD

	# PP -> t tbar Higgs
	create Herwig::MEPP2QQHiggs MEPP2ttbarH
	newdef MEPP2ttbarH:QuarkType Top

	# PP -> b bbar Higgs
	create Herwig::MEPP2QQHiggs MEPP2bbbarH
	newdef MEPP2bbbarH:QuarkType Bottom

	##########################################################
	# Hadron-Hadron NLO matrix elements in the Powheg scheme
	##########################################################
	library HwPowhegMEHadron.so

	# q qbar -> gamma/Z -> l+l-
	create Herwig::MEqq2gZ2ffPowheg PowhegMEqq2gZ2ff
	newdef PowhegMEqq2gZ2ff:Process 3
	newdef PowhegMEqq2gZ2ff:Coupling /Herwig/Shower/AlphaQCD

	# q qbar to W -> l nu
	create Herwig::MEqq2W2ffPowheg PowhegMEqq2W2ff
	newdef PowhegMEqq2W2ff:Process 2
	newdef PowhegMEqq2W2ff:Coupling /Herwig/Shower/AlphaQCD

	# PP->ZH
	create Herwig::MEPP2ZHPowheg PowhegMEPP2ZH
	newdef PowhegMEPP2ZH:Coupling /Herwig/Shower/AlphaQCD

	# PP->WH
	create Herwig::MEPP2WHPowheg PowhegMEPP2WH
	newdef PowhegMEPP2WH:Coupling /Herwig/Shower/AlphaQCD

	# hadron-hadron to higgs
	create Herwig::MEPP2HiggsPowheg PowhegMEHiggs
	newdef PowhegMEHiggs:ShapeScheme MassGenerator
	newdef PowhegMEHiggs:Process gg
	newdef PowhegMEHiggs:Coupling /Herwig/Shower/AlphaQCD

	# PP->VV
	create Herwig::MEPP2VVPowheg PowhegMEPP2VV
	newdef PowhegMEPP2VV:Coupling /Herwig/Shower/AlphaQCD

	# PP -> Higgs via VBF
	create Herwig::MEPP2HiggsVBFPowheg PowhegMEPP2HiggsVBF
	newdef PowhegMEPP2HiggsVBF:ShowerAlphaQCD /Herwig/Shower/AlphaQCD

	# PP -> diphoton NLO
	create Herwig::MEPP2GammaGammaPowheg MEGammaGammaPowheg
	set MEGammaGammaPowheg:Process 0
	set MEGammaGammaPowheg:Contribution 1
	set MEGammaGammaPowheg:ShowerAlphaQCD /Herwig/Shower/AlphaQCD
	set MEGammaGammaPowheg:ShowerAlphaQED /Herwig/Shower/AlphaQED

	##########################################################
	# DIS matrix elements
	##########################################################

	# neutral current
	create Herwig::MENeutralCurrentDIS MEDISNC
	newdef MEDISNC:Coupling /Herwig/Shower/AlphaQCD
	newdef MEDISNC:Contribution 0
	# charged current
	create Herwig::MEChargedCurrentDIS MEDISCC
	newdef MEDISCC:Coupling /Herwig/Shower/AlphaQCD
	newdef MEDISCC:Contribution 0

	# neutral current (POWHEG)
	create Herwig::MENeutralCurrentDIS PowhegMEDISNC
	newdef PowhegMEDISNC:Coupling /Herwig/Shower/AlphaQCD
	newdef PowhegMEDISNC:Contribution 1
	# charged current (POWHEG)
	create Herwig::MEChargedCurrentDIS PowhegMEDISCC
	newdef PowhegMEDISCC:Coupling /Herwig/Shower/AlphaQCD
	newdef PowhegMEDISCC:Contribution 1

	##########################################################
	# Gamma-Gamma matrix elements
	##########################################################

	# fermion-antiferimon
	create Herwig::MEGammaGamma2ff MEgg2ff HwMEGammaGamma.so

	# W+ W-
	create Herwig::MEGammaGamma2WW MEgg2WW HwMEGammaGamma.so

	##########################################################
	# Gamma-Hadron matrix elements
	##########################################################

	# gamma parton -> 2 jets
	create Herwig::MEGammaP2Jets MEGammaP2Jets HwMEGammaHadron.so

	##########################################################
	# Set up the Subprocesses
	#
	# Generic for all colliders
	##########################################################
	create ThePEG::SubProcessHandler SubProcess
	newdef SubProcess:PartonExtractor /Herwig/Partons/PPExtractor
	diff --git a/src/defaults/Partons.in b/src/defaults/Partons.in
	--- a/src/defaults/Partons.in
	+++ b/src/defaults/Partons.in
	@@ -1,36 +1,37 @@
	# -- ThePEG-repository --

	mkdir /Herwig/Partons
	cd /Herwig/Partons

	create ThePEG::NoRemnants NoRemnants
	create ThePEG::NoPDF NoPDF
	newdef NoPDF:RemnantHandler NoRemnants

	create ThePEG::UnResolvedRemnant LeptonRemnants UnResolvedRemnant.so

	create Herwig::HwRemDecayer RemnantDecayer HwShower.so
	newdef RemnantDecayer:PomeronStructure Gluon
	-newdef RemnantDecayer:colourDisrupt 0.7983
	+newdef RemnantDecayer:colourDisrupt 0.0
	newdef RemnantDecayer:ladderPower -0.088983
	newdef RemnantDecayer:ladderNorm 1.086029
	+newdef RemnantDecayer:PtDistribution 4

	create ThePEG::SoftRemnantHandler HadronRemnants
	newdef HadronRemnants:RemnantDecayer RemnantDecayer


	read PDF.in

	# Make particle extractors
	create ThePEG::PartonExtractor EEExtractor
	newdef EEExtractor:NoPDF NoPDF
	cp EEExtractor PPExtractor
	newdef PPExtractor:FlatSHatY 1
	newdef PPExtractor:FirstPDF HardLOPDF
	newdef PPExtractor:SecondPDF HardLOPDF
	cp EEExtractor EPExtractor
	newdef EPExtractor:FlatSHatY 1

	newdef PPExtractor:FirstPDF HardLOPDF
	newdef PPExtractor:SecondPDF HardLOPDF
	newdef EPExtractor:SecondPDF HardLOPDF
	diff --git a/src/defaults/Shower.in b/src/defaults/Shower.in
	--- a/src/defaults/Shower.in
	+++ b/src/defaults/Shower.in
	@@ -1,321 +1,321 @@
	# -- ThePEG-repository --

	############################################################
	# Setup of default parton shower
	#
	# Useful switches for users are marked near the top of
	# this file.
	#
	# Don't edit this file directly, but reset the switches
	# in your own input files!
	############################################################
	library HwMPI.so
	library HwShower.so
	library HwMatching.so

	mkdir /Herwig/Shower
	cd /Herwig/Shower

	create Herwig::QTildeShowerHandler ShowerHandler
	newdef ShowerHandler:MPIHandler /Herwig/UnderlyingEvent/MPIHandler
	newdef ShowerHandler:RemDecayer /Herwig/Partons/RemnantDecayer
	# use LO PDFs for Shower, can be changed later
	newdef ShowerHandler:PDFA /Herwig/Partons/ShowerLOPDF
	newdef ShowerHandler:PDFB /Herwig/Partons/ShowerLOPDF
	newdef ShowerHandler:PDFARemnant /Herwig/Partons/RemnantPDF
	newdef ShowerHandler:PDFBRemnant /Herwig/Partons/RemnantPDF

	#####################################
	# initial setup, don't change these!
	#####################################
	create Herwig::SplittingGenerator SplittingGenerator
	create Herwig::ShowerAlphaQCD AlphaQCD
	create Herwig::ShowerAlphaQED AlphaQED
	create Herwig::PartnerFinder PartnerFinder
	newdef PartnerFinder:PartnerMethod 0
	newdef PartnerFinder:ScaleChoice 0
	create Herwig::KinematicsReconstructor KinematicsReconstructor
	newdef KinematicsReconstructor:ReconstructionOption Colour3
	newdef KinematicsReconstructor:InitialStateReconOption SofterFraction
	newdef KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
	newdef KinematicsReconstructor:FinalFinalWeight Yes

	newdef /Herwig/Partons/RemnantDecayer:AlphaS AlphaQCD
	newdef /Herwig/Partons/RemnantDecayer:AlphaEM AlphaQED

	newdef ShowerHandler:PartnerFinder PartnerFinder
	newdef ShowerHandler:KinematicsReconstructor KinematicsReconstructor
	newdef ShowerHandler:SplittingGenerator SplittingGenerator
	newdef ShowerHandler:SpinCorrelations Yes
	newdef ShowerHandler:SoftCorrelations Singular

	##################################################################
	# Intrinsic pT
	#
	# Recommended:
	# 1.9 GeV for Tevatron W/Z production.
	# 2.1 GeV for LHC W/Z production at 10 TeV
	# 2.2 GeV for LHC W/Z production at 14 TeV
	#
	# Set all parameters to 0 to disable
	##################################################################
	-newdef ShowerHandler:IntrinsicPtGaussian 1.3*GeV
	+newdef ShowerHandler:IntrinsicPtGaussian 2.2*GeV
	newdef ShowerHandler:IntrinsicPtBeta 0
	newdef ShowerHandler:IntrinsicPtGamma 0*GeV
	newdef ShowerHandler:IntrinsicPtIptmax 0*GeV
	#############################################################
	# Set up truncated shower handler.
	#############################################################

	create Herwig::PowhegShowerHandler PowhegShowerHandler
	set PowhegShowerHandler:MPIHandler /Herwig/UnderlyingEvent/MPIHandler
	set PowhegShowerHandler:RemDecayer /Herwig/Partons/RemnantDecayer
	newdef PowhegShowerHandler:PDFA /Herwig/Partons/ShowerLOPDF
	newdef PowhegShowerHandler:PDFB /Herwig/Partons/ShowerLOPDF
	newdef PowhegShowerHandler:PDFARemnant /Herwig/Partons/RemnantPDF
	newdef PowhegShowerHandler:PDFBRemnant /Herwig/Partons/RemnantPDF
	newdef PowhegShowerHandler:MPIHandler /Herwig/UnderlyingEvent/MPIHandler
	newdef PowhegShowerHandler:RemDecayer /Herwig/Partons/RemnantDecayer
	newdef PowhegShowerHandler:PDFA /Herwig/Partons/ShowerLOPDF
	newdef PowhegShowerHandler:PDFB /Herwig/Partons/ShowerLOPDF
	newdef PowhegShowerHandler:PDFARemnant /Herwig/Partons/RemnantPDF
	newdef PowhegShowerHandler:PDFBRemnant /Herwig/Partons/RemnantPDF
	newdef PowhegShowerHandler:PartnerFinder PartnerFinder
	newdef PowhegShowerHandler:KinematicsReconstructor KinematicsReconstructor
	newdef PowhegShowerHandler:SplittingGenerator SplittingGenerator
	newdef PowhegShowerHandler:Interactions QCDandQED
	newdef PowhegShowerHandler:SpinCorrelations Yes
	newdef PowhegShowerHandler:SoftCorrelations Singular
	-newdef PowhegShowerHandler:IntrinsicPtGaussian 1.3*GeV
	+newdef PowhegShowerHandler:IntrinsicPtGaussian 2.2*GeV
	newdef PowhegShowerHandler:IntrinsicPtBeta 0
	newdef PowhegShowerHandler:IntrinsicPtGamma 0*GeV
	newdef PowhegShowerHandler:IntrinsicPtIptmax 0*GeV
	newdef PowhegShowerHandler:EvolutionScheme DotProduct

	#############################################################
	# End of interesting user servicable section.
	#
	# Anything that follows below should only be touched if you
	# know what you're doing.
	#
	# Really.
	#############################################################
	#
	# a few default values
	newdef ShowerHandler:MECorrMode 1
	newdef ShowerHandler:EvolutionScheme DotProduct
	newdef AlphaQCD:ScaleFactor 1.0
	newdef AlphaQCD:NPAlphaS 2
	newdef AlphaQCD:Qmin 0.935
	newdef AlphaQCD:NumberOfLoops 2
	newdef AlphaQCD:AlphaIn 0.1186
	#
	#
	# Lets set up all the splittings
	create Herwig::HalfHalfOneSplitFn QtoQGammaSplitFn
	set QtoQGammaSplitFn:InteractionType QED
	set QtoQGammaSplitFn:ColourStructure ChargedChargedNeutral
	set QtoQGammaSplitFn:AngularOrdered Yes
	set QtoQGammaSplitFn:StrictAO Yes

	create Herwig::HalfHalfOneSplitFn QtoQGSplitFn
	newdef QtoQGSplitFn:InteractionType QCD
	newdef QtoQGSplitFn:ColourStructure TripletTripletOctet
	set QtoQGSplitFn:AngularOrdered Yes
	set QtoQGSplitFn:StrictAO Yes

	create Herwig::OneOneOneSplitFn GtoGGSplitFn
	newdef GtoGGSplitFn:InteractionType QCD
	newdef GtoGGSplitFn:ColourStructure OctetOctetOctet
	set GtoGGSplitFn:AngularOrdered Yes
	set GtoGGSplitFn:StrictAO Yes

	create Herwig::OneOneOneMassiveSplitFn WtoWGammaSplitFn
	newdef WtoWGammaSplitFn:InteractionType QED
	newdef WtoWGammaSplitFn:ColourStructure ChargedChargedNeutral
	set WtoWGammaSplitFn:AngularOrdered Yes
	set WtoWGammaSplitFn:StrictAO Yes

	create Herwig::OneHalfHalfSplitFn GtoQQbarSplitFn
	newdef GtoQQbarSplitFn:InteractionType QCD
	newdef GtoQQbarSplitFn:ColourStructure OctetTripletTriplet
	set GtoQQbarSplitFn:AngularOrdered Yes
	set GtoQQbarSplitFn:StrictAO Yes

	create Herwig::OneHalfHalfSplitFn GammatoQQbarSplitFn
	newdef GammatoQQbarSplitFn:InteractionType QED
	newdef GammatoQQbarSplitFn:ColourStructure NeutralChargedCharged
	set GammatoQQbarSplitFn:AngularOrdered Yes
	set GammatoQQbarSplitFn:StrictAO Yes

	create Herwig::HalfOneHalfSplitFn QtoGQSplitFn
	newdef QtoGQSplitFn:InteractionType QCD
	newdef QtoGQSplitFn:ColourStructure TripletOctetTriplet
	set QtoGQSplitFn:AngularOrdered Yes
	set QtoGQSplitFn:StrictAO Yes

	create Herwig::HalfOneHalfSplitFn QtoGammaQSplitFn
	newdef QtoGammaQSplitFn:InteractionType QED
	newdef QtoGammaQSplitFn:ColourStructure ChargedNeutralCharged
	set QtoGammaQSplitFn:AngularOrdered Yes
	set QtoGammaQSplitFn:StrictAO Yes
	#
	# Now the Sudakovs
	create Herwig::PTCutOff PTCutOff
	newdef PTCutOff:pTmin 0.958*GeV

	create Herwig::SudakovFormFactor SudakovCommon
	newdef SudakovCommon:Alpha AlphaQCD
	newdef SudakovCommon:Cutoff PTCutOff
	newdef SudakovCommon:PDFmax 1.0

	cp SudakovCommon QtoQGSudakov
	newdef QtoQGSudakov:SplittingFunction QtoQGSplitFn
	newdef QtoQGSudakov:PDFmax 1.9

	cp SudakovCommon QtoQGammaSudakov
	set QtoQGammaSudakov:SplittingFunction QtoQGammaSplitFn
	set QtoQGammaSudakov:Alpha AlphaQED
	set QtoQGammaSudakov:PDFmax 1.9

	cp QtoQGammaSudakov LtoLGammaSudakov
	cp PTCutOff LtoLGammaPTCutOff
	# Technical parameter to stop evolution.
	set LtoLGammaPTCutOff:pTmin 0.000001
	set LtoLGammaSudakov:Cutoff LtoLGammaPTCutOff

	cp SudakovCommon GtoGGSudakov
	newdef GtoGGSudakov:SplittingFunction GtoGGSplitFn
	newdef GtoGGSudakov:PDFmax 2.0

	cp SudakovCommon WtoWGammaSudakov
	newdef WtoWGammaSudakov:SplittingFunction WtoWGammaSplitFn
	set WtoWGammaSudakov:Alpha AlphaQED

	cp SudakovCommon GtoQQbarSudakov
	newdef GtoQQbarSudakov:SplittingFunction GtoQQbarSplitFn
	newdef GtoQQbarSudakov:PDFmax 120.0

	cp SudakovCommon GammatoQQbarSudakov
	newdef GammatoQQbarSudakov:SplittingFunction GammatoQQbarSplitFn
	set GammatoQQbarSudakov:Alpha AlphaQED
	newdef GammatoQQbarSudakov:PDFmax 120.0

	cp SudakovCommon GtobbbarSudakov
	newdef GtobbbarSudakov:SplittingFunction GtoQQbarSplitFn
	newdef GtobbbarSudakov:PDFmax 40000.0

	cp SudakovCommon GtoccbarSudakov
	newdef GtoccbarSudakov:SplittingFunction GtoQQbarSplitFn
	newdef GtoccbarSudakov:PDFmax 2000.0

	cp SudakovCommon QtoGQSudakov
	newdef QtoGQSudakov:SplittingFunction QtoGQSplitFn

	cp SudakovCommon QtoGammaQSudakov
	newdef QtoGammaQSudakov:SplittingFunction QtoGammaQSplitFn
	set QtoGammaQSudakov:Alpha AlphaQED

	cp SudakovCommon utoGuSudakov
	newdef utoGuSudakov:SplittingFunction QtoGQSplitFn
	newdef utoGuSudakov:PDFFactor OverOneMinusZ
	newdef utoGuSudakov:PDFmax 5.0

	cp SudakovCommon dtoGdSudakov
	newdef dtoGdSudakov:SplittingFunction QtoGQSplitFn
	newdef dtoGdSudakov:PDFFactor OverOneMinusZ


	#
	# Now add the final splittings
	#
	do SplittingGenerator:AddFinalSplitting u->u,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting d->d,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting s->s,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting c->c,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting b->b,g; QtoQGSudakov
	do SplittingGenerator:AddFinalSplitting t->t,g; QtoQGSudakov
	#
	do SplittingGenerator:AddFinalSplitting g->g,g; GtoGGSudakov
	#
	do SplittingGenerator:AddFinalSplitting g->u,ubar; GtoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting g->d,dbar; GtoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting g->s,sbar; GtoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting g->c,cbar; GtoccbarSudakov
	do SplittingGenerator:AddFinalSplitting g->b,bbar; GtobbbarSudakov
	do SplittingGenerator:AddFinalSplitting g->t,tbar; GtoQQbarSudakov
	#
	do SplittingGenerator:AddFinalSplitting gamma->u,ubar; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->d,dbar; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->s,sbar; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->c,cbar; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->b,bbar; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->t,tbar; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->e-,e+; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->mu-,mu+; GammatoQQbarSudakov
	do SplittingGenerator:AddFinalSplitting gamma->tau-,tau+; GammatoQQbarSudakov
	#
	do SplittingGenerator:AddFinalSplitting u->u,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddFinalSplitting d->d,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddFinalSplitting s->s,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddFinalSplitting c->c,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddFinalSplitting b->b,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddFinalSplitting t->t,gamma; QtoQGammaSudakov

	do SplittingGenerator:AddFinalSplitting e-->e-,gamma; LtoLGammaSudakov
	do SplittingGenerator:AddFinalSplitting mu-->mu-,gamma; LtoLGammaSudakov
	do SplittingGenerator:AddFinalSplitting tau-->tau-,gamma; LtoLGammaSudakov

	do SplittingGenerator:AddFinalSplitting W+->W+,gamma; WtoWGammaSudakov
	#
	# Now lets add the initial splittings. Remember the form a->b,c; means
	# that the current particle b is given and we backward branch to new
	# particle a which is initial state and new particle c which is final state
	#
	do SplittingGenerator:AddInitialSplitting u->u,g; QtoQGSudakov
	do SplittingGenerator:AddInitialSplitting d->d,g; QtoQGSudakov
	do SplittingGenerator:AddInitialSplitting s->s,g; QtoQGSudakov
	do SplittingGenerator:AddInitialSplitting c->c,g; QtoQGSudakov
	do SplittingGenerator:AddInitialSplitting b->b,g; QtoQGSudakov
	do SplittingGenerator:AddInitialSplitting u->u,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddInitialSplitting d->d,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddInitialSplitting s->s,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddInitialSplitting c->c,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddInitialSplitting b->b,gamma; QtoQGammaSudakov
	do SplittingGenerator:AddInitialSplitting t->t,gamma; QtoQGammaSudakov

	do SplittingGenerator:AddInitialSplitting g->g,g; GtoGGSudakov
	#
	do SplittingGenerator:AddInitialSplitting g->d,dbar; GtoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting g->u,ubar; GtoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting g->s,sbar; GtoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting g->c,cbar; GtoccbarSudakov
	do SplittingGenerator:AddInitialSplitting g->b,bbar; GtobbbarSudakov
	#
	do SplittingGenerator:AddInitialSplitting gamma->d,dbar; GammatoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting gamma->u,ubar; GammatoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting gamma->s,sbar; GammatoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting gamma->c,cbar; GammatoQQbarSudakov
	do SplittingGenerator:AddInitialSplitting gamma->b,bbar; GammatoQQbarSudakov
	#
	do SplittingGenerator:AddInitialSplitting d->g,d; dtoGdSudakov
	do SplittingGenerator:AddInitialSplitting u->g,u; utoGuSudakov
	do SplittingGenerator:AddInitialSplitting s->g,s; QtoGQSudakov
	do SplittingGenerator:AddInitialSplitting c->g,c; QtoGQSudakov
	do SplittingGenerator:AddInitialSplitting b->g,b; QtoGQSudakov
	do SplittingGenerator:AddInitialSplitting dbar->g,dbar; dtoGdSudakov
	do SplittingGenerator:AddInitialSplitting ubar->g,ubar; utoGuSudakov
	do SplittingGenerator:AddInitialSplitting sbar->g,sbar; QtoGQSudakov
	do SplittingGenerator:AddInitialSplitting cbar->g,cbar; QtoGQSudakov
	do SplittingGenerator:AddInitialSplitting bbar->g,bbar; QtoGQSudakov
	#
	do SplittingGenerator:AddInitialSplitting d->gamma,d; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting u->gamma,u; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting s->gamma,s; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting c->gamma,c; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting b->gamma,b; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting dbar->gamma,dbar; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting ubar->gamma,ubar; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting sbar->gamma,sbar; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting cbar->gamma,cbar; QtoGammaQSudakov
	do SplittingGenerator:AddInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov
	diff --git a/src/defaults/UnderlyingEvent.in b/src/defaults/UnderlyingEvent.in
	--- a/src/defaults/UnderlyingEvent.in
	+++ b/src/defaults/UnderlyingEvent.in
	@@ -1,75 +1,82 @@
	# -- ThePEG-repository --

	################################################
	# Set up the handler for the underlying event.
	################################################

	mkdir /Herwig/UnderlyingEvent
	cd /Herwig/UnderlyingEvent

	######## cuts #################################################
	library SimpleKTCut.so

	# cut on pt. Without a specific matcher object, it works on all
	# particles
	create ThePEG::SimpleKTCut KtCut
	newdef KtCut:MinKT 3.06

	# create the cuts object for the Underlying Event
	create ThePEG::Cuts UECuts
	# This should always be 2*MinKT!!
	newdef UECuts:MHatMin 6.12

	insert UECuts:OneCuts 0 KtCut

	######## subprocess ###########################################
	create ThePEG::SubProcessHandler FastQCD

	create Herwig::MEQCD2to2Fast MEQCD2to2Fast HwMEHadronFast.so
	insert FastQCD:MatrixElements 0 MEQCD2to2Fast
	cp /Herwig/Partons/PPExtractor /Herwig/Partons/MPIExtractor
	newdef /Herwig/Partons/MPIExtractor:FlatSHatY 1
	newdef FastQCD:PartonExtractor /Herwig/Partons/MPIExtractor

	# always use LO PDF for MPI and remnants
	newdef /Herwig/Partons/MPIExtractor:FirstPDF /Herwig/Partons/MPIPDF
	newdef /Herwig/Partons/MPIExtractor:SecondPDF /Herwig/Partons/MPIPDF

	######## MPI Handler ##########################################
	library JetCuts.so
	library HwMPI.so
	create Herwig::MPIHandler MPIHandler

	# set the subprocesses and corresponding cuts
	insert MPIHandler:SubProcessHandlers 0 FastQCD
	insert MPIHandler:Cuts 0 UECuts

	# energy extrapolation for the MinKT cut
	newdef MPIHandler:EnergyExtrapolation Power
	-newdef MPIHandler:Power 0.420445
	-newdef MPIHandler:pTmin0 3.568157
	+newdef MPIHandler:Power 0.322
	+newdef MPIHandler:pTmin0 2.55
	newdef MPIHandler:ReferenceScale 7000.*GeV
	newdef MPIHandler:EEparamA 0.585*GeV
	newdef MPIHandler:EEparamB 37.45*GeV

	+# Paramters for the ladders for soft MPIs
	+newdef /Herwig/Partons/RemnantDecayer:ladderMult 0.466
	+newdef /Herwig/Partons/RemnantDecayer:ladderbFactor 0.46
	+
	+# Ratio of diffractive events in MinBias runs
	+newdef MPIHandler:DiffractiveRatio 0.133
	+
	# The inverse hadron radius
	-newdef MPIHandler:InvRadius 1.489997
	+newdef MPIHandler:InvRadius 1.58

	# Set the details of the soft model

	# Flag to decide whether additional soft interactions (i.e. pt < ptmin)
	# should be simulated at all

	newdef MPIHandler:softInt Yes

	# Flag to decide whether to use the 2 component model. In this model,
	# the hadron radius (or better: mean parton-parton separation) can be
	# different for hard and soft scatters. The radius for the soft scatters
	# is computed during the run startup so that the elastic t-slope, B_el,
	# coincides with current measurements (for CMenergy < 1.8 TeV) or
	# Donnachie-Landshoff pomeron fit (CMenergy > 1.8 TeV). If this model is
	# disabled one has to take care that this observable is correctly
	# described by adjusting InvRadius properly.

	newdef MPIHandler:twoComp Yes
	newdef MPIHandler:DLmode 2
	diff --git a/src/defaults/quarks.in b/src/defaults/quarks.in
	--- a/src/defaults/quarks.in
	+++ b/src/defaults/quarks.in
	@@ -1,35 +1,35 @@
	# -- ThePEG-repository --

	#
	# file containing the particle data for the quarks
	#
	create ThePEG::ConstituentParticleData d
	setup d 1 d 0.005 0. 0 0 -1 3 2 1 0.325
	create ThePEG::ConstituentParticleData dbar
	setup dbar -1 dbar 0.005 0. 0 0 1 -3 2 1 0.325
	makeanti dbar d
	create ThePEG::ConstituentParticleData ubar
	setup ubar -2 ubar 0.0023 0. 0 0 -2 -3 2 1 0.325
	create ThePEG::ConstituentParticleData u
	setup u 2 u 0.0023 0. 0 0 2 3 2 1 0.325
	makeanti u ubar
	create ThePEG::ConstituentParticleData c
	setup c 4 c 1.25 0. 0 0 2 3 2 1 1.6
	create ThePEG::ConstituentParticleData cbar
	setup cbar -4 cbar 1.25 0. 0 0 -2 -3 2 1 1.6
	makeanti cbar c
	create ThePEG::ConstituentParticleData sbar
	-setup sbar -3 sbar 0.095 0. 0 0 1 -3 2 1 0.5
	+setup sbar -3 sbar 0.095 0. 0 0 1 -3 2 1 0.45
	create ThePEG::ConstituentParticleData s
	-setup s 3 s 0.095 0. 0 0 -1 3 2 1 0.5
	+setup s 3 s 0.095 0. 0 0 -1 3 2 1 0.45
	makeanti s sbar
	create ThePEG::ConstituentParticleData b
	setup b 5 b 4.2 0. 0 0 -1 3 2 1 5
	create ThePEG::ConstituentParticleData bbar
	setup bbar -5 bbar 4.2 0. 0 0 1 -3 2 1 5
	makeanti bbar b
	create ThePEG::ParticleData t
	setup t 6 t 174.2 1.4 14 0 2 3 2 0
	create ThePEG::ParticleData tbar
	setup tbar -6 tbar 174.2 1.4 14 0 -2 -3 2 0
	makeanti tbar t
	diff --git a/src/snippets/BaryonicReconnection.in b/src/snippets/BaryonicReconnection.in
	--- a/src/snippets/BaryonicReconnection.in
	+++ b/src/snippets/BaryonicReconnection.in
	@@ -1,17 +1,10 @@
	# Set strange quark mass to 0.45 in order to allow alternative gluon splitting
	set /Herwig/Particles/s:ConstituentMass 0.45*GeV
	set /Herwig/Particles/sbar:ConstituentMass 0.45*GeV

	# Use Baryonic Colour Reconnection Model
	set /Herwig/Hadronization/ColourReconnector:Algorithm Baryonic

	# Allow alternative gluon splitting
	set /Herwig/Hadronization/PartonSplitter:Split uds

	-# Parameters for the Baryonic Reconnection Model
	-set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.772606
	-set /Herwig/Hadronization/ColourReconnector:ReconnectionProbabilityBaryonic 0.477612
	-set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 3.053252
	-set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.282032
	-set /Herwig/Hadronization/HadronSelector:PwtSquark 0.291717
	-set /Herwig/Hadronization/PartonSplitter:SplitPwtSquark 0.824135
	diff --git a/src/snippets/SoftTune.in b/src/snippets/SoftTune.in
	--- a/src/snippets/SoftTune.in
	+++ b/src/snippets/SoftTune.in
	@@ -1,8 +1,16 @@
	-# Recommended set of parameters for MB/UE runs
	+# Tund for UE and MB (these are the default parameters)

	-set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.5
	-set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 3.502683
	-set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.402055
	-set /Herwig/UnderlyingEvent/MPIHandler:Power 0.416852
	-set /Herwig/Partons/RemnantDecayer:ladderPower -0.08
	-set /Herwig/Partons/RemnantDecayer:ladderNorm 0.95
	+set /Herwig/MatrixElements/MEMinBias:csNorm 0.01
	+set /Herwig/UnderlyingEvent/MPIHandler:DLmode 2
	+set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 2.55
	+set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.58
	+set /Herwig/UnderlyingEvent/MPIHandler:Power 0.322
	+set /Herwig/Partons/RemnantDecayer:ladderMult 0.466
	+set /Herwig/Partons/RemnantDecayer:ladderbFactor 0.46
	+set /Herwig/UnderlyingEvent/MPIHandler:DiffractiveRatio 0.133
	+set /Herwig/Partons/RemnantDecayer:PtDistribution 4
	+set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.89
	+set /Herwig/Hadronization/ColourReconnector:ReconnectionProbabilityBaryonic 0.57
	+set /Herwig/Hadronization/HadronSelector:PwtSquark 0.291717
	+set /Herwig/Hadronization/PartonSplitter:SplitPwtSquark 0.824135
	+

File Metadata

Mime Type: text/x-diff
Expires: Sun, Feb 23, 3:04 PM (23 m, 33 s)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4486584
Default Alt Text: (50 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions