No OneTemporary
Actions

Size

55 KB

Subscribers

None

View Options

	diff --git a/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.cc b/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.cc
	--- a/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.cc
	+++ b/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.cc
	@@ -1,110 +1,176 @@
	// -- C++ --
	//
	// MatchboxTopIndividualMTScale.cc is a part of Herwig - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2012 The Herwig Collaboration
	//
	// Herwig is licenced under version 2 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	//
	// This is the implementation of the non-inlined, non-templated member
	// functions of the MatchboxTopIndividualMTScale class.
	//

	#include "MatchboxTopIndividualMTScale.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	+#include "ThePEG/Interface/Switch.h"
	#include "ThePEG/Interface/Parameter.h"
	#include "ThePEG/Interface/Reference.h"
	#include "ThePEG/EventRecord/Particle.h"
	#include "ThePEG/Repository/UseRandom.h"
	#include "ThePEG/Repository/EventGenerator.h"
	#include "ThePEG/Utilities/DescribeClass.h"
	#include "Herwig/MatrixElement/Matchbox/Base/MatchboxMEBase.h"

	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"

	using namespace Herwig;

	MatchboxTopIndividualMTScale::MatchboxTopIndividualMTScale() :
	-theFactor(1.) {}
	+ theShowerScaleMode(1), theFactor(1.) {}

	MatchboxTopIndividualMTScale::~MatchboxTopIndividualMTScale() {}

	IBPtr MatchboxTopIndividualMTScale::clone() const {
	return new_ptr(*this);
	}

	IBPtr MatchboxTopIndividualMTScale::fullclone() const {
	return new_ptr(*this);
	}

	Energy2 MatchboxTopIndividualMTScale::renormalizationScale() const {

	size_t k = 2;
	int top = -1;
	int antitop = -1;

	while ( (top == -1 \|\| antitop == -1) && k < mePartonData().size() ){
	if ( mePartonData()[k]->id() == 6 ) {
	if ( top < 0 )
	top = k;
	else
	assert(false);
	} else if ( mePartonData()[k]->id() == -6 ) {
	if ( antitop < 0 )
	antitop = k;
	else
	assert(false);
	}
	k++;
	}

	if ( top < 2 \|\| antitop < 2 ){
	throw Exception() << "MatchboxTopIndividualMTScale: Could not find a top-antitop-pair in the final state!\n"
	<< Exception::runerror;
	}

	// Not using .mt() as this is signed and not what we want.
	Energy topMt = sqrt(meMomenta()[top].mt2());
	return sqr(topMt*theFactor);

	}

	Energy2 MatchboxTopIndividualMTScale::factorizationScale() const {
	return(renormalizationScale());
	}

	+Energy2 MatchboxTopIndividualMTScale::showerScale() const {
	+
	+ if ( theShowerScaleMode == 1 )
	+ return factorizationScale();
	+
	+ else {
	+ assert(theShowerScaleMode == 2);
	+ size_t k = 2;
	+ int top = -1;
	+ int antitop = -1;
	+ int emission = -1;
	+
	+ if ( mePartonData().size() > 5 )
	+ assert(false && "MatchboxTopIndividualMTScale.cc: mePartonData().size() > 5.\n");
	+
	+ while ( (top == -1 \|\| antitop == -1 \|\| emission == -1) && k < mePartonData().size() ){
	+ if ( mePartonData()[k]->id() == 6 ) {
	+ if ( top < 0 )
	+ top = k;
	+ else
	+ assert(false);
	+ } else if ( mePartonData()[k]->id() == -6 ) {
	+ if ( antitop < 0 )
	+ antitop = k;
	+ else
	+ assert(false);
	+ }
	+ else if ( abs(mePartonData()[k]->id()) < 6 \|\| mePartonData()[k]->id() == 21 ) {
	+ if ( emission < 0 )
	+ emission = k;
	+ else
	+ assert(false && "More than one emission in MatchboxTopIndividualMTScale.\n");
	+ }
	+
	+ k++;
	+ }
	+
	+ if ( top < 2 \|\| antitop < 2 ){
	+ throw Exception() << "MatchboxTopIndividualMTScale: Could not find a top-antitop-pair in the final state!\n"
	+ << Exception::runerror;
	+ }
	+
	+
	+ if ( mePartonData().size() > 4 && emission < 2 ){
	+ throw Exception() << "MatchboxTopIndividualMTScale: Could not find an emission in a state where mePartonData().size() > 4.\n"
	+ << Exception::runerror;
	+ }
	+
	+ if ( emission < 2 )
	+ return (1./2.)*(meMomenta()[top].mt2()+meMomenta()[antitop].mt2());
	+
	+ else
	+ return (1./3.)*(meMomenta()[top].mt2()+meMomenta()[antitop].mt2() + meMomenta()[emission].mt2());
	+ }
	+}
	+
	// If needed, insert default implementations of virtual function defined
	// in the InterfacedBase class here (using ThePEG-interfaced-impl in Emacs).


	void MatchboxTopIndividualMTScale::persistentOutput(PersistentOStream & os) const {
	-os << theFactor;
	+ os << theShowerScaleMode << theFactor;
	}

	void MatchboxTopIndividualMTScale::persistentInput(PersistentIStream & is, int) {
	-is >> theFactor;
	+ is >> theShowerScaleMode >> theFactor;
	}

	// * Attention * The following static variable is needed for the type
	// description system in ThePEG. Please check that the template arguments
	// are correct (the class and its base class), and that the constructor
	// arguments are correct (the class name and the name of the dynamically
	// loadable library where the class implementation can be found).
	DescribeClass<MatchboxTopIndividualMTScale,MatchboxScaleChoice>
	describeHerwigMatchboxTopIndividualMTScale("Herwig::MatchboxTopIndividualMTScale", "HwMatchboxScales.so");

	void MatchboxTopIndividualMTScale::Init() {

	static ClassDocumentation<MatchboxTopIndividualMTScale> documentation
	- ("MatchboxTopIndividualMTScale implements the linear sum of the transverse masses of the top and antitop quark as a scale choice.");
	+ ("MatchboxTopIndividualMTScale implements the transverse mass of the top quark as the scale choice.");
	+
	+ static Switch<MatchboxTopIndividualMTScale, unsigned int> interfaceShowerScaleMode
	+ ("ShowerScaleMode",
	+ "Choose the definition of the shower hard scale.",
	+ &MatchboxTopIndividualMTScale::theShowerScaleMode, 1, false, false);
	+ static SwitchOption FactorizationScale
	+ (interfaceShowerScaleMode,"FactorizationScale","Use the factorization scale.", 1);
	+ static SwitchOption MeanMT2
	+ (interfaceShowerScaleMode,"MeanMT2","Use the mean squared transverse mass of the outgoing particles.", 2);

	static Parameter<MatchboxTopIndividualMTScale,double> interfaceMultiplicationFactor
	("MultiplicationFactor",
	"Set a multiplicative factor to include in the scale choice definition.",
	&MatchboxTopIndividualMTScale::theFactor, 1., 0., 0,
	false, false, Interface::lowerlim);

	}

	diff --git a/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.h b/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.h
	--- a/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.h
	+++ b/MatrixElement/Matchbox/Scales/MatchboxTopIndividualMTScale.h
	@@ -1,125 +1,135 @@
	// -- C++ --
	//
	// MatchboxTopIndividualMTScale.h is a part of Herwig - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2014 The Herwig Collaboration
	//
	// Herwig is licenced under version 2 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	#ifndef Herwig_MatchboxTopIndividualMTScale_H
	#define Herwig_MatchboxTopIndividualMTScale_H
	//
	// This is the declaration of the MatchboxTopIndividualMTScale class.
	//

	#include "Herwig/MatrixElement/Matchbox/Utility/MatchboxScaleChoice.h"

	namespace Herwig {

	using namespace ThePEG;

	/**
	* \ingroup Matchbox
	- * \author Daniel Rauch
	+ * \author Stephen Webster
	*
	- * \brief MatchboxTopIndividualMTScale implements a scale choice related to the quadratic
	- * sum of transverse masses of a top and antitop quark.
	+ * \brief MatchboxTopIndividualMTScale implements the scale as
	+ * the transverse masses of the top quark in a top-pair process.
	*
	*/
	class MatchboxTopIndividualMTScale: public MatchboxScaleChoice {

	public:

	/** @name Standard constructors and destructors. */
	//@{
	/**
	* The default constructor.
	*/
	MatchboxTopIndividualMTScale();

	/**
	* The destructor.
	*/
	virtual ~MatchboxTopIndividualMTScale();
	//@}

	public:

	/**
	* Return the renormalization scale.
	*/
	virtual Energy2 renormalizationScale() const;

	/**
	* Return the factorization scale.
	*/
	virtual Energy2 factorizationScale() const;

	+ /**
	+ * Return the shower hard scale.
	+ */
	+ virtual Energy2 showerScale() const;
	+
	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);
	//@}

	/**
	* The standard Init function used to initialize the interfaces.
	* Called exactly once for each class by the class description system
	* before the main function starts or
	* when this class is dynamically loaded.
	*/
	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:


	// If needed, insert declarations of virtual function defined in the
	// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).


	private:

	/**
	+ * Switch to choose the definition of the shower hard scale.
	+ */
	+ unsigned int theShowerScaleMode;
	+
	+ /**
	* Multiplicative factor for the scale definition
	**/
	double theFactor;

	/**
	* The assignment operator is private and must never be called.
	* In fact, it should not even be implemented.
	*/
	MatchboxTopIndividualMTScale & operator=(const MatchboxTopIndividualMTScale &);

	};

	}

	#endif /* Herwig_MatchboxTopIndividualMTScale_H */
	diff --git a/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.cc b/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.cc
	--- a/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.cc
	+++ b/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.cc
	@@ -1,110 +1,178 @@
	// -- C++ --
	//
	// MatchboxTopLinearSumMTScale.cc is a part of Herwig - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2012 The Herwig Collaboration
	//
	// Herwig is licenced under version 2 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	//
	// This is the implementation of the non-inlined, non-templated member
	// functions of the MatchboxTopLinearSumMTScale class.
	//

	#include "MatchboxTopLinearSumMTScale.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	+#include "ThePEG/Interface/Switch.h"
	#include "ThePEG/Interface/Parameter.h"
	#include "ThePEG/Interface/Reference.h"
	#include "ThePEG/EventRecord/Particle.h"
	#include "ThePEG/Repository/UseRandom.h"
	#include "ThePEG/Repository/EventGenerator.h"
	#include "ThePEG/Utilities/DescribeClass.h"
	#include "Herwig/MatrixElement/Matchbox/Base/MatchboxMEBase.h"

	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"

	using namespace Herwig;

	MatchboxTopLinearSumMTScale::MatchboxTopLinearSumMTScale() :
	-theFactor(1.) {}
	+ theShowerScaleMode(1) , theFactor(1.) {}

	MatchboxTopLinearSumMTScale::~MatchboxTopLinearSumMTScale() {}

	IBPtr MatchboxTopLinearSumMTScale::clone() const {
	return new_ptr(*this);
	}

	IBPtr MatchboxTopLinearSumMTScale::fullclone() const {
	return new_ptr(*this);
	}

	Energy2 MatchboxTopLinearSumMTScale::renormalizationScale() const {

	size_t k = 2;
	int top = -1;
	int antitop = -1;

	while ( (top == -1 \|\| antitop == -1) && k < mePartonData().size() ){
	if ( mePartonData()[k]->id() == 6 ) {
	if ( top < 0 )
	top = k;
	else
	assert(false);
	} else if ( mePartonData()[k]->id() == -6 ) {
	if ( antitop < 0 )
	antitop = k;
	else
	assert(false);
	}
	k++;
	}

	if ( top < 2 \|\| antitop < 2 ){
	throw Exception() << "MatchboxTopLinearSumMTScale: Could not find a top-antitop-pair in the final state!\n"
	<< Exception::runerror;
	}

	// Not using .mt() as this is signed and not what we want.
	Energy sum = sqrt(meMomenta()[top].mt2()) + sqrt(meMomenta()[antitop].mt2());
	return sqr(sum*theFactor);

	}

	Energy2 MatchboxTopLinearSumMTScale::factorizationScale() const {
	return(renormalizationScale());
	}

	+Energy2 MatchboxTopLinearSumMTScale::showerScale() const {
	+
	+ if ( theShowerScaleMode == 1 )
	+ return factorizationScale();
	+
	+ else {
	+ assert(theShowerScaleMode == 2);
	+ size_t k = 2;
	+ int top = -1;
	+ int antitop = -1;
	+ int emission = -1;
	+
	+ if ( mePartonData().size() > 5 )
	+ assert(false && "MatchboxTopLinearSumMTScale.cc: mePartonData().size() > 5.\n");
	+
	+ while ( (top == -1 \|\| antitop == -1 \|\| emission == -1) && k < mePartonData().size() ){
	+ if ( mePartonData()[k]->id() == 6 ) {
	+ if ( top < 0 )
	+ top = k;
	+ else
	+ assert(false);
	+ } else if ( mePartonData()[k]->id() == -6 ) {
	+ if ( antitop < 0 )
	+ antitop = k;
	+ else
	+ assert(false);
	+ }
	+ else if ( abs(mePartonData()[k]->id()) < 6 \|\| mePartonData()[k]->id() == 21 ) {
	+ if ( emission < 0 )
	+ emission = k;
	+ else
	+ assert(false && "More than one emission in MatchboxTopLinearSumMTScale.\n");
	+ }
	+
	+ k++;
	+ }
	+
	+ if ( top < 2 \|\| antitop < 2 ){
	+ throw Exception() << "MatchboxTopLinearSumMTScale: Could not find a top-antitop-pair in the final state!\n"
	+ << Exception::runerror;
	+ }
	+
	+
	+ if ( mePartonData().size() > 4 && emission < 2 ){
	+ throw Exception() << "MatchboxTopLinearSumMTScale: Could not find an emission in a state where mePartonData().size() > 4.\n"
	+ << Exception::runerror;
	+ }
	+
	+ if ( emission < 2 )
	+ return (1./2.)*(meMomenta()[top].mt2()+meMomenta()[antitop].mt2());
	+
	+ else
	+ return (1./3.)*(meMomenta()[top].mt2()+meMomenta()[antitop].mt2() + meMomenta()[emission].mt2());
	+ }
	+}
	+
	// If needed, insert default implementations of virtual function defined
	// in the InterfacedBase class here (using ThePEG-interfaced-impl in Emacs).


	void MatchboxTopLinearSumMTScale::persistentOutput(PersistentOStream & os) const {
	-os << theFactor;
	+ os << theShowerScaleMode << theFactor;
	}

	void MatchboxTopLinearSumMTScale::persistentInput(PersistentIStream & is, int) {
	-is >> theFactor;
	+ is >> theShowerScaleMode >> theFactor;
	}

	// * Attention * The following static variable is needed for the type
	// description system in ThePEG. Please check that the template arguments
	// are correct (the class and its base class), and that the constructor
	// arguments are correct (the class name and the name of the dynamically
	// loadable library where the class implementation can be found).
	DescribeClass<MatchboxTopLinearSumMTScale,MatchboxScaleChoice>
	describeHerwigMatchboxTopLinearSumMTScale("Herwig::MatchboxTopLinearSumMTScale", "HwMatchboxScales.so");

	void MatchboxTopLinearSumMTScale::Init() {

	static ClassDocumentation<MatchboxTopLinearSumMTScale> documentation
	("MatchboxTopLinearSumMTScale implements the linear sum of the transverse masses of the top and antitop quark as a scale choice.");
	+

	+ static Switch<MatchboxTopLinearSumMTScale, unsigned int> interfaceShowerScaleMode
	+ ("ShowerScaleMode",
	+ "Choose the definition of the shower hard scale.",
	+ &MatchboxTopLinearSumMTScale::theShowerScaleMode, 1, false, false);
	+ static SwitchOption FactorizationScale
	+ (interfaceShowerScaleMode,"FactorizationScale","Use the factorization scale.", 1);
	+ static SwitchOption MeanMT2
	+ (interfaceShowerScaleMode,"MeanMT2","Use the mean squared transverse mass of the outgoing particles.", 2);
	+
	+
	static Parameter<MatchboxTopLinearSumMTScale,double> interfaceMultiplicationFactor
	("MultiplicationFactor",
	"Set a multiplicative factor to include in the scale choice definition.",
	&MatchboxTopLinearSumMTScale::theFactor, 1., 0., 0,
	false, false, Interface::lowerlim);

	}

	diff --git a/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.h b/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.h
	--- a/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.h
	+++ b/MatrixElement/Matchbox/Scales/MatchboxTopLinearSumMTScale.h
	@@ -1,125 +1,136 @@
	// -- C++ --
	//
	// MatchboxTopLinearSumMTScale.h is a part of Herwig - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2014 The Herwig Collaboration
	//
	// Herwig is licenced under version 2 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	#ifndef Herwig_MatchboxTopLinearSumMTScale_H
	#define Herwig_MatchboxTopLinearSumMTScale_H
	//
	// This is the declaration of the MatchboxTopLinearSumMTScale class.
	//

	#include "Herwig/MatrixElement/Matchbox/Utility/MatchboxScaleChoice.h"

	namespace Herwig {

	using namespace ThePEG;

	/**
	* \ingroup Matchbox
	- * \author Daniel Rauch
	+ * \author Stephen Webster
	*
	- * \brief MatchboxTopLinearSumMTScale implements a scale choice related to the linear
	- * sum of transverse masses of a top and antitop quark.
	+ * \brief MatchboxTopLinearSumMTScale implements the scale as the linear
	+ * sum of the transverse masses of the top and antitop quarks in
	+ * a top pair process.
	*
	*/
	class MatchboxTopLinearSumMTScale: public MatchboxScaleChoice {

	public:

	/** @name Standard constructors and destructors. */
	//@{
	/**
	* The default constructor.
	*/
	MatchboxTopLinearSumMTScale();

	/**
	* The destructor.
	*/
	virtual ~MatchboxTopLinearSumMTScale();
	//@}

	public:

	/**
	* Return the renormalization scale.
	*/
	virtual Energy2 renormalizationScale() const;

	/**
	* Return the factorization scale.
	*/
	virtual Energy2 factorizationScale() const;

	+ /**
	+ * Return the shower hard scale.
	+ */
	+ virtual Energy2 showerScale() const;
	+
	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);
	//@}

	/**
	* The standard Init function used to initialize the interfaces.
	* Called exactly once for each class by the class description system
	* before the main function starts or
	* when this class is dynamically loaded.
	*/
	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:


	// If needed, insert declarations of virtual function defined in the
	// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).


	private:

	/**
	+ * Switch to choose the definition of the shower hard scale.
	+ */
	+ unsigned int theShowerScaleMode;
	+
	+ /**
	* Multiplicative factor for the scale definition
	**/
	double theFactor;

	/**
	* The assignment operator is private and must never be called.
	* In fact, it should not even be implemented.
	*/
	MatchboxTopLinearSumMTScale & operator=(const MatchboxTopLinearSumMTScale &);

	};

	}

	#endif /* Herwig_MatchboxTopLinearSumMTScale_H */
	diff --git a/MatrixElement/Matchbox/Scales/MatchboxTopMassScale.cc b/MatrixElement/Matchbox/Scales/MatchboxTopMassScale.cc
	--- a/MatrixElement/Matchbox/Scales/MatchboxTopMassScale.cc
	+++ b/MatrixElement/Matchbox/Scales/MatchboxTopMassScale.cc
	@@ -1,175 +1,172 @@
	// -- C++ --
	//
	// MatchboxTopMassScale.cc is a part of Herwig - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2014 The Herwig Collaboration
	//
	// Herwig is licenced under version 2 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	//
	// This is the implementation of the non-inlined, non-templated member
	// functions of the MatchboxTopMassScale class.
	//

	#include "MatchboxTopMassScale.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	#include "ThePEG/Interface/Switch.h"
	#include "ThePEG/Interface/Parameter.h"
	#include "ThePEG/Interface/Reference.h"
	#include "ThePEG/EventRecord/Particle.h"
	#include "ThePEG/Repository/UseRandom.h"
	#include "ThePEG/Repository/EventGenerator.h"
	#include "ThePEG/Utilities/DescribeClass.h"
	#include "Herwig/MatrixElement/Matchbox/Base/MatchboxMEBase.h"

	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"

	using namespace Herwig;

	MatchboxTopMassScale::MatchboxTopMassScale() :
	theShowerScaleMode(1) {}

	MatchboxTopMassScale::~MatchboxTopMassScale() {}

	IBPtr MatchboxTopMassScale::clone() const {
	return new_ptr(*this);
	}

	IBPtr MatchboxTopMassScale::fullclone() const {
	return new_ptr(*this);
	}

	Energy2 MatchboxTopMassScale::renormalizationScale() const {

	size_t k = 2;
	int top = -1;
	int antitop = -1;

	while ( (top == -1 \|\| antitop == -1) && k < mePartonData().size() ){
	if ( mePartonData()[k]->id() == 6 ) {
	if ( top < 0 )
	top = k;
	else
	assert(false);
	} else if ( mePartonData()[k]->id() == -6 ) {
	if ( antitop < 0 )
	antitop = k;
	else
	assert(false);
	}
	k++;
	}

	if ( top < 2 \|\| antitop < 2 ){
	throw Exception() << "MatchboxTopMassScale: Could not find a top-antitop-pair in the final state!\n"
	<< Exception::runerror;
	}
	// cerr << " sqrt(TopMassScale) = "
	// << sqrt((meMomenta()[top]+meMomenta()[antitop]).m2())/GeV
	// << "\n" << flush;
	return((meMomenta()[top]+meMomenta()[antitop]).m2());

	}

	Energy2 MatchboxTopMassScale::factorizationScale() const {
	return(renormalizationScale());
	}

	-
	Energy2 MatchboxTopMassScale::showerScale() const {

	if ( theShowerScaleMode == 1 )
	return factorizationScale();

	else {
	assert(theShowerScaleMode == 2);
	size_t k = 2;
	int top = -1;
	int antitop = -1;
	int emission = -1;

	if ( mePartonData().size() > 5 )
	assert(false && "MatchboxTopMassScale.cc: mePartonData().size() > 5.\n");

	while ( (top == -1 \|\| antitop == -1 \|\| emission == -1) && k < mePartonData().size() ){
	if ( mePartonData()[k]->id() == 6 ) {
	if ( top < 0 )
	top = k;
	else
	assert(false);
	} else if ( mePartonData()[k]->id() == -6 ) {
	if ( antitop < 0 )
	antitop = k;
	else
	assert(false);
	}
	else if ( abs(mePartonData()[k]->id()) < 6 \|\| mePartonData()[k]->id() == 21 ) {
	if ( emission < 0 )
	emission = k;
	else
	assert(false && "More than one emission in MatchboxTopMassScale.\n");
	}

	k++;
	}

	if ( top < 2 \|\| antitop < 2 ){
	throw Exception() << "MatchboxTopMassScale: Could not find a top-antitop-pair in the final state!\n"
	<< Exception::runerror;
	}


	if ( mePartonData().size() > 4 && emission < 2 ){
	throw Exception() << "MatchboxTopMassScale: Could not find an emission in a state where mePartonData().size() > 4.\n"
	<< Exception::runerror;
	}

	if ( emission < 2 )
	return (1./2.)*(meMomenta()[top].mt2()+meMomenta()[antitop].mt2());

	else
	return (1./3.)*(meMomenta()[top].mt2()+meMomenta()[antitop].mt2() + meMomenta()[emission].mt2());
	}
	}

	-
	-
	// If needed, insert default implementations of virtual function defined
	// in the InterfacedBase class here (using ThePEG-interfaced-impl in Emacs).


	void MatchboxTopMassScale::persistentOutput(PersistentOStream & os) const {
	os << theShowerScaleMode;
	}

	void MatchboxTopMassScale::persistentInput(PersistentIStream & is, int) {
	is >> theShowerScaleMode;
	}


	// * Attention * The following static variable is needed for the type
	// description system in ThePEG. Please check that the template arguments
	// are correct (the class and its base class), and that the constructor
	// arguments are correct (the class name and the name of the dynamically
	// loadable library where the class implementation can be found).
	DescribeClass<MatchboxTopMassScale,MatchboxScaleChoice>
	describeHerwigMatchboxTopMassScale("Herwig::MatchboxTopMassScale", "HwMatchboxScales.so");

	void MatchboxTopMassScale::Init() {

	static ClassDocumentation<MatchboxTopMassScale> documentation
	("MatchboxTopMassScale implements invariant mass of top-antitop pair as scale choice.");


	static Switch<MatchboxTopMassScale, unsigned int> interfaceShowerScaleMode
	("ShowerScaleMode",
	"Choose the definition of the shower hard scale.",
	&MatchboxTopMassScale::theShowerScaleMode, 1, false, false);
	static SwitchOption FactorizationScale
	(interfaceShowerScaleMode,"FactorizationScale","Use the factorization scale.", 1);
	static SwitchOption MeanMT2
	(interfaceShowerScaleMode,"MeanMT2","Use the mean squared transverse mass of the outgoing particles.", 2);

	}

	diff --git a/src/defaults/MatchboxDefaults.in.in b/src/defaults/MatchboxDefaults.in.in
	--- a/src/defaults/MatchboxDefaults.in.in
	+++ b/src/defaults/MatchboxDefaults.in.in
	@@ -1,761 +1,764 @@
	# -- ThePEG-repository --

	################################################################################
	#
	# Default setup for Matchbox matrix element generation.
	# You do not need to make any change in here; processes of
	# interest can be chosen in the standard input files.
	#
	################################################################################

	################################################################################
	# Load libraries
	################################################################################

	library JetCuts.so
	library FastJetFinder.so
	library HwMatchboxScales.so
	library HwMatchboxCuts.so
	library HwSampling.so
	library HwColorFull.so
	library HwMatchboxBuiltin.so

	################################################################################
	# Integration/sampling
	################################################################################

	mkdir /Herwig/Samplers
	cd /Herwig/Samplers

	create Herwig::BinSampler FlatBinSampler
	set FlatBinSampler:InitialPoints 1000
	set FlatBinSampler:UseAllIterations No

	create Herwig::CellGridSampler CellGridSampler
	set CellGridSampler:InitialPoints 10000
	set CellGridSampler:ExplorationPoints 500
	set CellGridSampler:ExplorationSteps 4
	set CellGridSampler:Gain 0.3
	set CellGridSampler:Epsilon 1.0
	set CellGridSampler:MinimumSelection 0.000001
	set CellGridSampler:NIterations 1
	set CellGridSampler:EnhancementFactor 1
	set CellGridSampler:UseAllIterations No

	set CellGridSampler:RemapperPoints 50000
	set CellGridSampler:RemapperMinSelection 0.00001
	set CellGridSampler:RemapChannelDimension Yes
	set CellGridSampler:LuminosityMapperBins 20
	set CellGridSampler:GeneralMapperBins 0

	set CellGridSampler:HalfPoints No
	set CellGridSampler:MaxNewMax 30
	set CellGridSampler:NonZeroInPresampling Yes

	create Herwig::MonacoSampler MonacoSampler
	set MonacoSampler:InitialPoints 15000
	set MonacoSampler:NIterations 4
	set MonacoSampler:EnhancementFactor 1.2
	set MonacoSampler:UseAllIterations No
	set MonacoSampler:RemapChannelDimension No
	set MonacoSampler:LuminosityMapperBins 0
	set MonacoSampler:HalfPoints No
	set MonacoSampler:MaxNewMax 30
	set MonacoSampler:NonZeroInPresampling Yes

	create Herwig::GeneralSampler Sampler
	set Sampler:UpdateAfter 1000

	set Sampler:BinSampler CellGridSampler

	set Sampler:AddUpSamplers Off
	set Sampler:GlobalMaximumWeight Off
	set Sampler:FlatSubprocesses Off
	set Sampler:MinSelection 0.000001
	set Sampler:AlmostUnweighted Off
	set Sampler:RunCombinationData Off
	set Sampler:WriteGridsOnFinish No
	set Sampler:MaxEnhancement 1.1

	################################################################################
	# Setup the factory object
	################################################################################

	mkdir /Herwig/MatrixElements/Matchbox
	cd /Herwig/MatrixElements/Matchbox

	create Herwig::MatchboxFactory Factory

	do Factory:StartParticleGroup p
	insert Factory:ParticleGroup 0 /Herwig/Particles/b
	insert Factory:ParticleGroup 0 /Herwig/Particles/bbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/c
	insert Factory:ParticleGroup 0 /Herwig/Particles/cbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/s
	insert Factory:ParticleGroup 0 /Herwig/Particles/sbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/d
	insert Factory:ParticleGroup 0 /Herwig/Particles/dbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/u
	insert Factory:ParticleGroup 0 /Herwig/Particles/ubar
	insert Factory:ParticleGroup 0 /Herwig/Particles/g
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup pbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/b
	insert Factory:ParticleGroup 0 /Herwig/Particles/bbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/c
	insert Factory:ParticleGroup 0 /Herwig/Particles/cbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/s
	insert Factory:ParticleGroup 0 /Herwig/Particles/sbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/d
	insert Factory:ParticleGroup 0 /Herwig/Particles/dbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/u
	insert Factory:ParticleGroup 0 /Herwig/Particles/ubar
	insert Factory:ParticleGroup 0 /Herwig/Particles/g
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup j
	insert Factory:ParticleGroup 0 /Herwig/Particles/b
	insert Factory:ParticleGroup 0 /Herwig/Particles/bbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/c
	insert Factory:ParticleGroup 0 /Herwig/Particles/cbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/s
	insert Factory:ParticleGroup 0 /Herwig/Particles/sbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/d
	insert Factory:ParticleGroup 0 /Herwig/Particles/dbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/u
	insert Factory:ParticleGroup 0 /Herwig/Particles/ubar
	insert Factory:ParticleGroup 0 /Herwig/Particles/g
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup u
	insert Factory:ParticleGroup 0 /Herwig/Particles/u
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup ubar
	insert Factory:ParticleGroup 0 /Herwig/Particles/ubar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup d
	insert Factory:ParticleGroup 0 /Herwig/Particles/d
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup dbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/dbar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup s
	insert Factory:ParticleGroup 0 /Herwig/Particles/s
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup sbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/sbar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup c
	insert Factory:ParticleGroup 0 /Herwig/Particles/c
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup cbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/cbar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup b
	insert Factory:ParticleGroup 0 /Herwig/Particles/b
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup bbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/bbar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup t
	insert Factory:ParticleGroup 0 /Herwig/Particles/t
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup tbar
	insert Factory:ParticleGroup 0 /Herwig/Particles/tbar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup g
	insert Factory:ParticleGroup 0 /Herwig/Particles/g
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup gamma
	insert Factory:ParticleGroup 0 /Herwig/Particles/gamma
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup h0
	insert Factory:ParticleGroup 0 /Herwig/Particles/h0
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup W+
	insert Factory:ParticleGroup 0 /Herwig/Particles/W+
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup W-
	insert Factory:ParticleGroup 0 /Herwig/Particles/W-
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup Z0
	insert Factory:ParticleGroup 0 /Herwig/Particles/Z0
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup e+
	insert Factory:ParticleGroup 0 /Herwig/Particles/e+
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup e-
	insert Factory:ParticleGroup 0 /Herwig/Particles/e-
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup mu+
	insert Factory:ParticleGroup 0 /Herwig/Particles/mu+
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup mu-
	insert Factory:ParticleGroup 0 /Herwig/Particles/mu-
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup tau+
	insert Factory:ParticleGroup 0 /Herwig/Particles/tau+
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup tau-
	insert Factory:ParticleGroup 0 /Herwig/Particles/tau-
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu_e
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_e
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu_mu
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_mu
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu_tau
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_tau
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu_ebar
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_ebar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu_mubar
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_mubar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu_taubar
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_taubar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup l
	insert Factory:ParticleGroup 0 /Herwig/Particles/e+
	insert Factory:ParticleGroup 0 /Herwig/Particles/mu+
	insert Factory:ParticleGroup 0 /Herwig/Particles/e-
	insert Factory:ParticleGroup 0 /Herwig/Particles/mu-
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup nu
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_e
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_mu
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_ebar
	insert Factory:ParticleGroup 0 /Herwig/Particles/nu_mubar
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup l+
	insert Factory:ParticleGroup 0 /Herwig/Particles/e+
	insert Factory:ParticleGroup 0 /Herwig/Particles/mu+
	do Factory:EndParticleGroup

	do Factory:StartParticleGroup l-
	insert Factory:ParticleGroup 0 /Herwig/Particles/e-
	insert Factory:ParticleGroup 0 /Herwig/Particles/mu-
	do Factory:EndParticleGroup

	################################################################################
	# Default settings for hard process widths
	################################################################################

	set /Herwig/Particles/mu+:HardProcessWidth 0*GeV
	set /Herwig/Particles/mu-:HardProcessWidth 0*GeV
	set /Herwig/Particles/tau+:HardProcessWidth 0*GeV
	set /Herwig/Particles/tau-:HardProcessWidth 0*GeV

	################################################################################
	# Setup amplitudes
	################################################################################

	cd /Herwig/MatrixElements/Matchbox
	mkdir Amplitudes
	cd Amplitudes

	create ColorFull::TraceBasis TraceBasis

	create Herwig::MatchboxHybridAmplitude GenericProcesses

	@LOAD_MADGRAPH@ HwMatchboxMadGraph.so
	@CREATE_MADGRAPH@ Herwig::MadGraphAmplitude MadGraph
	@SET_MADGRAPH@ MadGraph:ColourBasis TraceBasis

	@LOAD_GOSAM@ HwMatchboxGoSam.so
	@CREATE_GOSAM@ Herwig::GoSamAmplitude GoSam

	@LOAD_NJET@ HwMatchboxNJet.so
	@CREATE_NJET@ Herwig::NJetsAmplitude NJet

	@DO_NJET@ NJet:Massless 5
	@DO_NJET@ NJet:Massless -5

	@LOAD_OPENLOOPS@ HwMatchboxOpenLoops.so
	@CREATE_OPENLOOPS@ Herwig::OpenLoopsAmplitude OpenLoops

	@LOAD_VBFNLO@ HwMatchboxVBFNLO.so
	@CREATE_VBFNLO@ Herwig::VBFNLOAmplitude VBFNLO

	mkdir Builtin
	cd Builtin

	create Herwig::SimpleColourBasis SimpleColourBasis
	create Herwig::SimpleColourBasis2 SimpleColourBasis2

	create Herwig::MatchboxAmplitudellbarqqbar Amplitudellbarqqbar
	set Amplitudellbarqqbar:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudellbarqqbarg Amplitudellbarqqbarg
	set Amplitudellbarqqbarg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudellbarqqbargg Amplitudellbarqqbargg
	set Amplitudellbarqqbargg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudellbarqqbarqqbar Amplitudellbarqqbarqqbar
	set Amplitudellbarqqbarqqbar:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudelnuqqbar Amplitudelnuqqbar
	set Amplitudelnuqqbar:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudelnuqqbarg Amplitudelnuqqbarg
	set Amplitudelnuqqbarg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudelnuqqbargg Amplitudelnuqqbargg
	set Amplitudelnuqqbargg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudelnuqqbarqqbar Amplitudelnuqqbarqqbar
	set Amplitudelnuqqbarqqbar:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudehgg Amplitudehgg
	set Amplitudehgg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudehggg Amplitudehggg
	set Amplitudehggg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudehqqbarg Amplitudehqqbarg
	set Amplitudehqqbarg:ColourBasis SimpleColourBasis

	create Herwig::MatchboxAmplitudeqqbarttbar Amplitudeqqbarttbar
	set Amplitudeqqbarttbar:ColourBasis SimpleColourBasis2

	create Herwig::MatchboxAmplitudeqqbarttbarg Amplitudeqqbarttbarg
	set Amplitudeqqbarttbarg:ColourBasis SimpleColourBasis2

	create Herwig::MatchboxAmplitudeggttbar Amplitudeggttbar
	set Amplitudeggttbar:ColourBasis SimpleColourBasis2

	create Herwig::MatchboxAmplitudeggttbarg Amplitudeggttbarg
	set Amplitudeggttbarg:ColourBasis SimpleColourBasis2


	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudellbarqqbar
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudellbarqqbarg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudellbarqqbargg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudellbarqqbarqqbar

	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudelnuqqbar
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudelnuqqbarg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudelnuqqbargg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudelnuqqbarqqbar

	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudehgg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudehggg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudehqqbarg

	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudeqqbarttbar
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudeqqbarttbarg
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudeggttbar
	insert /Herwig/MatrixElements/Matchbox/Factory:Amplitudes 0 Amplitudeggttbarg

	################################################################################
	# Setup phasespace generators
	################################################################################

	cd /Herwig/MatrixElements/Matchbox
	mkdir Phasespace
	cd Phasespace

	create Herwig::PhasespaceCouplings PhasespaceCouplings

	create Herwig::MatchboxRambo Rambo
	set Rambo:CouplingData PhasespaceCouplings

	create Herwig::FlatInvertiblePhasespace InvertiblePhasespace
	set InvertiblePhasespace:CouplingData PhasespaceCouplings

	create Herwig::FlatInvertibleLabframePhasespace InvertibleLabframePhasespace
	set InvertibleLabframePhasespace:CouplingData PhasespaceCouplings
	set InvertibleLabframePhasespace:LogSHat False

	create Herwig::TreePhasespaceChannels TreePhasespaceChannels
	create Herwig::TreePhasespace TreePhasespace
	set TreePhasespace:ChannelMap TreePhasespaceChannels
	set TreePhasespace:M0 0.0001*GeV
	set TreePhasespace:MC 0.00005*GeV
	set TreePhasespace:CouplingData PhasespaceCouplings

	do TreePhasespace:SetPhysicalCoupling 21 -1 1 0.059
	do TreePhasespace:SetPhysicalCoupling 21 -2 2 0.059
	do TreePhasespace:SetPhysicalCoupling 21 -3 3 0.059
	do TreePhasespace:SetPhysicalCoupling 21 -4 4 0.059
	do TreePhasespace:SetPhysicalCoupling 21 -5 5 0.059
	do TreePhasespace:SetPhysicalCoupling 21 -6 6 0.059
	do TreePhasespace:SetPhysicalCoupling 21 1 -1 0.059
	do TreePhasespace:SetPhysicalCoupling 21 2 -2 0.059
	do TreePhasespace:SetPhysicalCoupling 21 3 -3 0.059
	do TreePhasespace:SetPhysicalCoupling 21 4 -4 0.059
	do TreePhasespace:SetPhysicalCoupling 21 5 -5 0.059
	do TreePhasespace:SetPhysicalCoupling 21 6 -6 0.059
	do TreePhasespace:SetPhysicalCoupling 1 21 1 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 2 21 2 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 3 21 3 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 4 21 4 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 5 21 5 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 6 21 6 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -1 21 -1 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -2 21 -2 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -3 21 -3 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -4 21 -4 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -5 21 -5 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -6 21 -6 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 1 1 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 2 2 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 3 3 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 4 4 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 5 5 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling 6 6 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -1 -1 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -2 -2 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -3 -3 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -4 -4 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -5 -5 21 0.15733333333333333333
	do TreePhasespace:SetPhysicalCoupling -6 -6 21 0.15733333333333333333
	do TreePhasespace:SetCoupling 25 -1 1 0
	do TreePhasespace:SetCoupling 25 -2 2 0
	do TreePhasespace:SetCoupling 25 -3 3 0.00000001184279069851
	do TreePhasespace:SetCoupling 25 -4 4 0.00000205034465001885
	do TreePhasespace:SetCoupling 25 -5 5 0.00002314757096085280
	do TreePhasespace:SetCoupling 25 -6 6 0.03982017320025470767
	do TreePhasespace:SetCoupling 25 -11 11 0.00000000000034264835
	do TreePhasespace:SetCoupling 25 -12 12 0
	do TreePhasespace:SetCoupling 25 -13 13 0.00000001464912263400
	do TreePhasespace:SetCoupling 25 -14 14 0
	do TreePhasespace:SetCoupling 25 -15 15 0.00000414359033108195
	do TreePhasespace:SetCoupling 25 -16 16 0
	do TreePhasespace:SetCoupling 22 -1 1 0.00083932358497608365
	do TreePhasespace:SetCoupling 22 -2 2 0.00335729433990433461
	do TreePhasespace:SetCoupling 22 -3 3 0.00083932358497608365
	do TreePhasespace:SetCoupling 22 -4 4 0.00335729433990433461
	do TreePhasespace:SetCoupling 22 -5 5 0.00083932358497608365
	do TreePhasespace:SetCoupling 22 -6 6 0.00335729433990433461
	do TreePhasespace:SetCoupling 22 -11 11 0.00755391226478475287
	do TreePhasespace:SetCoupling 22 -13 13 0.00755391226478475287
	do TreePhasespace:SetCoupling 22 -15 15 0.00755391226478475287
	do TreePhasespace:SetCoupling 24 -2 1 0.01652748072644379386
	do TreePhasespace:SetCoupling 24 -4 1 0.00382028458188709739
	do TreePhasespace:SetCoupling 24 -6 1 0.00014707756360995175
	do TreePhasespace:SetCoupling 24 -2 3 0.00382265953677814621
	do TreePhasespace:SetCoupling 24 -4 3 0.01651340063673257587
	do TreePhasespace:SetCoupling 24 -6 3 0.00068534412570265868
	do TreePhasespace:SetCoupling 24 -2 5 0.00005954351191129535
	do TreePhasespace:SetCoupling 24 -4 5 0.00069891529650865192
	do TreePhasespace:SetCoupling 24 -6 5 0.01694947628265615369
	do TreePhasespace:SetCoupling 24 -12 11 0.01696396350749155147
	do TreePhasespace:SetCoupling 24 -14 13 0.01696396350749155147
	do TreePhasespace:SetCoupling 24 -16 15 0.01696396350749155147
	do TreePhasespace:SetCoupling -24 2 -1 0.01652748072644379386
	do TreePhasespace:SetCoupling -24 4 -1 0.00382028458188709739
	do TreePhasespace:SetCoupling -24 6 -1 0.00014707756360995175
	do TreePhasespace:SetCoupling -24 2 -3 0.00382265953677814621
	do TreePhasespace:SetCoupling -24 4 -3 0.01651340063673257587
	do TreePhasespace:SetCoupling -24 6 -3 0.00068534412570265868
	do TreePhasespace:SetCoupling -24 2 -5 0.00005954351191129535
	do TreePhasespace:SetCoupling -24 4 -5 0.00069891529650865192
	do TreePhasespace:SetCoupling -24 6 -5 0.01694947628265615369
	do TreePhasespace:SetCoupling -24 12 -11 0.01696396350749155147
	do TreePhasespace:SetCoupling -24 14 -13 0.01696396350749155147
	do TreePhasespace:SetCoupling -24 16 -15 0.01696396350749155147
	do TreePhasespace:SetCoupling 23 -1 1 0.00407649129960709158
	do TreePhasespace:SetCoupling 23 -2 2 0.00317809816318353030
	do TreePhasespace:SetCoupling 23 -3 3 0.00407649129960709158
	do TreePhasespace:SetCoupling 23 -4 4 0.00317809816318353030
	do TreePhasespace:SetCoupling 23 -5 5 0.00407649129960709158
	do TreePhasespace:SetCoupling 23 -6 6 0.00317809816318353030
	do TreePhasespace:SetCoupling 23 -11 11 0.00276049468148072129
	do TreePhasespace:SetCoupling 23 -12 12 0.00545567409075140513
	do TreePhasespace:SetCoupling 23 -13 13 0.00276049468148072129
	do TreePhasespace:SetCoupling 23 -14 14 0.00545567409075140513
	do TreePhasespace:SetCoupling 23 -15 15 0.00276049468148072129
	do TreePhasespace:SetCoupling 23 -16 16 0.00545567409075140513
	do TreePhasespace:SetCoupling 21 21 21 0.354
	do TreePhasespace:SetCoupling 25 21 21 0.00000000016160437564
	do TreePhasespace:SetCoupling 25 25 25 0.18719783125611995353
	do TreePhasespace:SetCoupling 25 22 22 0.00000000006295673620
	do TreePhasespace:SetCoupling 25 24 -24 219.30463760755686425818
	do TreePhasespace:SetCoupling 25 23 23 362.91922658249853887524
	do TreePhasespace:SetCoupling 22 24 -24 0.00755391226478475287
	do TreePhasespace:SetCoupling 23 24 -24 0.02637401475019835008

	@CREATE_VBFNLO@ Herwig::VBFNLOPhasespace VBFNLOPhasespace
	@SET_VBFNLO@ VBFNLOPhasespace:CouplingData PhasespaceCouplings

	set /Herwig/MatrixElements/Matchbox/Factory:Phasespace TreePhasespace


	################################################################################
	# Setup utilities for matching
	################################################################################

	cd /Herwig/MatrixElements/Matchbox

	create Herwig::HardScaleProfile HardScaleProfile

	create Herwig::MEMatching MEMatching
	set MEMatching:RestrictPhasespace On
	set MEMatching:HardScaleProfile /Herwig/MatrixElements/Matchbox/HardScaleProfile
	set MEMatching:BornScaleInSubtraction BornScale
	set MEMatching:RealEmissionScaleInSubtraction RealScale
	set MEMatching:EmissionScaleInSubtraction RealScale
	set MEMatching:BornScaleInSplitting ShowerScale
	set MEMatching:RealEmissionScaleInSplitting ShowerScale
	set MEMatching:EmissionScaleInSplitting ShowerScale
	set MEMatching:TruncatedShower Yes
	set MEMatching:MaxPtIsMuF Yes
	set MEMatching:FFPtCut 1.0*GeV
	set MEMatching:FIPtCut 1.0*GeV
	set MEMatching:IIPtCut 1.0*GeV
	set MEMatching:SafeCut 0.*GeV

	create Herwig::ShowerApproximationGenerator MECorrectionHandler
	set MECorrectionHandler:ShowerApproximation MEMatching
	set MECorrectionHandler:Phasespace /Herwig/MatrixElements/Matchbox/Phasespace/InvertiblePhasespace
	set MECorrectionHandler:PresamplingPoints 50000
	set MECorrectionHandler:FreezeGrid 100000

	create Herwig::DipoleMatching DipoleMatching HwDipoleMatching.so
	# set in DipoleShowerDefaults.in as not available at this point
	# set DipoleMatching:ShowerHandler /Herwig/DipoleShower/DipoleShowerHandler
	set DipoleMatching:BornScaleInSubtraction BornScale
	set DipoleMatching:RealEmissionScaleInSubtraction BornScale
	set DipoleMatching:EmissionScaleInSubtraction BornScale
	set DipoleMatching:FFPtCut 1.0*GeV
	set DipoleMatching:FIPtCut 1.0*GeV
	set DipoleMatching:IIPtCut 1.0*GeV
	set DipoleMatching:SafeCut 4.*GeV

	create Herwig::QTildeMatching QTildeMatching HwQTildeMatching.so
	set QTildeMatching:ShowerHandler /Herwig/Shower/ShowerHandler
	set QTildeMatching:BornScaleInSubtraction BornScale
	set QTildeMatching:RealEmissionScaleInSubtraction BornScale
	set QTildeMatching:EmissionScaleInSubtraction BornScale
	set QTildeMatching:QTildeFinder /Herwig/Shower/PartnerFinder
	set QTildeMatching:SafeCut 4.*GeV
	# just a dummy, since SudakovCommonn can't be used
	# it's only used to get the value of the kinCutoffScale
	set QTildeMatching:QTildeSudakov /Herwig/Shower/QtoQGSudakov

	################################################################################
	# Setup utilities for process generation
	################################################################################

	cd /Herwig/MatrixElements/Matchbox
	mkdir Utility
	cd Utility

	create Herwig::Tree2toNGenerator DiagramGenerator
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/FFGVertex
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/GGGVertex

	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/FFPVertex
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/FFZVertex

	cp /Herwig/Vertices/FFWVertex /Herwig/Vertices/FFWMatchboxVertex
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/FFWMatchboxVertex

	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/WWHVertex
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/WWWVertex
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/HGGVertex
	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/HHHVertex

	cp /Herwig/Vertices/FFHVertex /Herwig/Vertices/TTHVertex
	set /Herwig/Vertices/TTHVertex:Fermion 6

	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/TTHVertex

	cp /Herwig/Vertices/FFHVertex /Herwig/Vertices/BBHVertex
	set /Herwig/Vertices/BBHVertex:Fermion 5

	cp /Herwig/Vertices/FFHVertex /Herwig/Vertices/TauTauHVertex
	set /Herwig/Vertices/TauTauHVertex:Fermion 15

	insert DiagramGenerator:Vertices 0 /Herwig/Vertices/TauTauHVertex

	cp /Herwig/Vertices/FFHVertex /Herwig/Vertices/MuMuHVertex
	set /Herwig/Vertices/MuMuHVertex:Fermion 13

	create Herwig::ProcessData ProcessData

	set /Herwig/MatrixElements/Matchbox/Factory:DiagramGenerator DiagramGenerator
	set /Herwig/MatrixElements/Matchbox/Factory:ProcessData ProcessData

	################################################################################
	# Setup jet cuts
	################################################################################

	cd /Herwig/Cuts

	create Herwig::MatchboxFactoryMatcher MatchboxJetMatcher
	set MatchboxJetMatcher:Group j

	create ThePEG::FastJetFinder JetFinder
	set JetFinder:UnresolvedMatcher MatchboxJetMatcher
	set JetFinder:Variant AntiKt
	set JetFinder:RecombinationScheme E
	set JetFinder:Mode Inclusive
	set JetFinder:ConeRadius 0.7

	create ThePEG::JetRegion FirstJet
	set FirstJet:PtMin 20.*GeV
	do FirstJet:YRange -5.0 5.0
	set FirstJet:Fuzzy Yes
	set FirstJet:EnergyCutWidth 4.0*GeV
	set FirstJet:RapidityCutWidth 0.4
	insert FirstJet:Accepts[0] 1

	create ThePEG::JetRegion SecondJet
	set SecondJet:PtMin 20.*GeV
	do SecondJet:YRange -5.0 5.0
	set SecondJet:Fuzzy Yes
	set SecondJet:EnergyCutWidth 4.0*GeV
	set SecondJet:RapidityCutWidth 0.4
	insert SecondJet:Accepts[0] 2

	create ThePEG::JetRegion ThirdJet
	set ThirdJet:PtMin 20.*GeV
	do ThirdJet:YRange -5.0 5.0
	set ThirdJet:Fuzzy Yes
	set ThirdJet:EnergyCutWidth 4.0*GeV
	set ThirdJet:RapidityCutWidth 0.4
	insert ThirdJet:Accepts[0] 3

	create ThePEG::JetRegion FourthJet
	set FourthJet:PtMin 20.*GeV
	do FourthJet:YRange -5.0 5.0
	set FourthJet:Fuzzy Yes
	set FourthJet:EnergyCutWidth 4.0*GeV
	set FourthJet:RapidityCutWidth 0.4
	insert FourthJet:Accepts[0] 4

	create ThePEG::FuzzyTheta FuzzyTheta
	set FuzzyTheta:EnergyWidth 4.0*GeV
	set FuzzyTheta:RapidityWidth 0.4
	set FuzzyTheta:AngularWidth 0.4

	create ThePEG::NJetsCut NJetsCut
	set NJetsCut:UnresolvedMatcher MatchboxJetMatcher
	set NJetsCut:NJetsMin 2

	create ThePEG::JetCuts JetCuts
	set JetCuts:UnresolvedMatcher MatchboxJetMatcher
	set JetCuts:Ordering OrderPt

	create Herwig::IdentifiedParticleCut IdentifiedParticleCut

	cp IdentifiedParticleCut LeptonCut
	set LeptonCut:Matcher /Herwig/Matchers/Lepton

	cp IdentifiedParticleCut ChargedLeptonCut
	set ChargedLeptonCut:Matcher /Herwig/Matchers/ChargedLepton

	cp IdentifiedParticleCut BottomQuarkCut
	set BottomQuarkCut:Matcher /Herwig/Matchers/Bottom

	cp IdentifiedParticleCut TopQuarkCut
	set TopQuarkCut:Matcher /Herwig/Matchers/Top

	cp IdentifiedParticleCut WBosonCut
	set WBosonCut:Matcher /Herwig/Matchers/WBoson

	cp IdentifiedParticleCut ZBosonCut
	set ZBosonCut:Matcher /Herwig/Matchers/ZBoson

	cp IdentifiedParticleCut HiggsBosonCut
	set HiggsBosonCut:Matcher /Herwig/Matchers/HiggsBoson

	cp IdentifiedParticleCut PhotonCut
	set PhotonCut:Matcher /Herwig/Matchers/Photon

	create Herwig::FrixionePhotonSeparationCut PhotonIsolationCut
	set PhotonIsolationCut:UnresolvedMatcher MatchboxJetMatcher

	create Herwig::MatchboxDeltaRCut MatchboxDeltaRCut

	cp MatchboxDeltaRCut LeptonDeltaRCut
	set LeptonDeltaRCut:FirstMatcher /Herwig/Matchers/Lepton
	set LeptonDeltaRCut:SecondMatcher /Herwig/Matchers/Lepton

	cp MatchboxDeltaRCut ChargedLeptonDeltaRCut
	set ChargedLeptonDeltaRCut:FirstMatcher /Herwig/Matchers/ChargedLepton
	set ChargedLeptonDeltaRCut:SecondMatcher /Herwig/Matchers/ChargedLepton

	create Herwig::InvariantMassCut InvariantMassCut

	cp InvariantMassCut LeptonPairMassCut
	set LeptonPairMassCut:FirstMatcher /Herwig/Matchers/Lepton
	set LeptonPairMassCut:SecondMatcher /Herwig/Matchers/Lepton

	cp InvariantMassCut ChargedLeptonPairMassCut
	set ChargedLeptonPairMassCut:FirstMatcher /Herwig/Matchers/ChargedLepton
	set ChargedLeptonPairMassCut:SecondMatcher /Herwig/Matchers/ChargedLepton

	create Herwig::MissingPtCut MissingPtCut
	set MissingPtCut:Matcher /Herwig/Matchers/Neutrino



	################################################################################
	# Setup scale choices
	################################################################################

	cd /Herwig/MatrixElements/Matchbox
	mkdir Scales
	cd Scales

	create Herwig::MatchboxScaleChoice SHatScale
	cp SHatScale FixedScale
	set FixedScale:FixedScale 100.*GeV
	create Herwig::MatchboxPtScale MaxJetPtScale
	set MaxJetPtScale:JetFinder /Herwig/Cuts/JetFinder
	create Herwig::MatchboxLeptonMassScale LeptonPairMassScale
	create Herwig::MatchboxLeptonPtScale LeptonPairPtScale
	create Herwig::MatchboxHtScale HTScale
	create Herwig::MatchboxTopMassScale TopPairMassScale
	create Herwig::MatchboxTopMTScale TopPairMTScale
	+create Herwig::MatchboxTopLinearSumMTScale TopPairLinearMTScale
	+create Herwig::MatchboxTopIndividualMTScale TopPairIndividualMTScale
	+
	set HTScale:JetFinder /Herwig/Cuts/JetFinder
	set HTScale:IncludeMT No
	set HTScale:JetPtCut 15.*GeV
	cp HTScale HTPrimeScale
	set HTPrimeScale:IncludeMT Yes
	set HTPrimeScale:JetPtCut 15.*GeV
	cp LeptonPairMassScale LeptonQ2Scale

	set /Herwig/MatrixElements/Matchbox/Factory:ScaleChoice LeptonPairMassScale

	cd /

File Metadata

Mime Type: text/x-diff
Expires: Tue, Nov 19, 3:00 PM (1 d, 14 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 3798823
Default Alt Text: (55 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions