Page MenuHomeHEPForge
Files F11221491

MatchboxMEPP2llbar.h
No OneTemporary
Actions

MatchboxMEPP2llbar.h
View Options

// -*- C++ -*-
//
// MatchboxMEPP2llbar.h 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.
//
#ifndef HERWIG_MatchboxMEPP2llbar_H
#define HERWIG_MatchboxMEPP2llbar_H
//
// MatchboxMEPP2llbar
//
#include "Herwig++/MatrixElement/Matchbox/Base/MatchboxMEBase.h"
#include "Herwig++/MatrixElement/Matchbox/Builtin/Processes/MatchboxMEllbarqqbar.h"
namespace Herwig {
using namespace ThePEG;
/**
* \ingroup Matchbox
* \author Simon Platzer
*
* \brief MatchboxMEPP2llbar implements charged lepton
* pair production in hadron-hadron collisions.
*
* @see \ref MatchboxMEPP2llbarInterfaces "The interfaces"
* defined for MatchboxMEPP2llbar.
*/
class MatchboxMEPP2llbar: public MatchboxMEBase, public MatchboxMEllbarqqbar {
public:
/** @name Standard constructors and destructors. */
//@{
/**
* The default constructor.
*/
MatchboxMEPP2llbar();
/**
* The destructor.
*/
virtual ~MatchboxMEPP2llbar();
//@}
public:
/**
* Return the order in \f$\alpha_S\f$ in which this matrix element
* is given.
*/
virtual unsigned int orderInAlphaS() const { return 0; }
/**
* Return the order in \f$\alpha_{EM}\f$ in which this matrix
* element is given. Returns 0.
*/
virtual unsigned int orderInAlphaEW() const { return 2; }
/**
* Return the matrix element for the kinematical configuation
* previously provided by the last call to setKinematics(), suitably
* scaled by sHat() to give a dimension-less number.
*/
virtual double me2() const;
/**
* Evaluate the finite virtual correction for the
* variables supplied through the Born XComb object
* and possible additional random numbers.
*/
virtual double oneLoopInterference() const {
return
(lastAlphaS()/(2.*Constants::pi)) *
((SM().Nc()*SM().Nc()-1.0)/(2.*SM().Nc())) *
( sqr(Constants::pi) - 8. ) *
me2();
}
/**
* Return true, if this matrix element is capable of calculating
* one-loop (QCD) corrections.
*/
virtual bool haveOneLoop() const { return true; };
/**
* Return true, if one loop corrections have been calculated in
* dimensional reduction. Otherwise conventional dimensional
* regularization is assumed. Note that renormalization is always
* assumed to be MSbar.
*/
virtual bool isDR() const { return false; }
/**
* Return true, if one loop corrections are given in the conventions
* of the integrated dipoles.
*/
virtual bool isCS() const { return true; }
/**
* Return the renormalization scale for the last generated phasespace point.
*/
virtual Energy2 factorizationScale() const;
/**
* Return the renormalization scale for the last generated phasespace point.
*/
virtual Energy2 renormalizationScale() const;
/**
* Return the number of light flavours, this matrix
* element is calculated for.
*/
virtual unsigned int nLight() const { return theQuarkFlavours.size(); }
/**
* Return the colour correlated matrix element squared with
* respect to the given two partons as appearing in mePartonData(),
* suitably scaled by sHat() to give a dimension-less number.
*/
virtual double colourCorrelatedME2(pair<int,int>) const;
/**
* Return the colour and spin correlated matrix element squared for
* the gluon indexed by the first argument using the given
* correlation tensor.
*/
virtual double spinColourCorrelatedME2(pair<int,int> emitterSpectator,
const SpinCorrelationTensor& c) const;
/**
* Add all possible diagrams with the add() function.
*/
virtual void getDiagrams() const;
/**
* With the information previously supplied with the
* setKinematics(...) method, a derived class may optionally
* override this method to weight the given diagrams with their
* (although certainly not physical) relative probabilities.
*/
virtual Selector<DiagramIndex> diagrams(const DiagramVector &) const;
/**
* Return a Selector with possible colour geometries for the selected
* diagram weighted by their relative probabilities.
*/
virtual Selector<const ColourLines *>
colourGeometries(tcDiagPtr diag) 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;
//@}
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
protected:
/** @name Standard Interfaced functions. */
//@{
/**
* Initialize this object after the setup phase before saving an
* EventGenerator to disk.
* @throws InitException if object could not be initialized properly.
*/
virtual void doinit();
//@}
private:
/**
* The lepton flavours to be considered.
*/
PDVector theLeptonFlavours;
/**
* The quark flavours to be considered.
*/
PDVector theQuarkFlavours;
/**
* A user defined scale; if zero, the center of mass
* energy is used.
*/
Energy theUserScale;
private:
/**
* The static object used to initialize the description of this class.
* Indicates that this is a concrete class with persistent data.
*/
static ClassDescription<MatchboxMEPP2llbar> initMatchboxMEPP2llbar;
/**
* The assignment operator is private and must never be called.
* In fact, it should not even be implemented.
*/
MatchboxMEPP2llbar & operator=(const MatchboxMEPP2llbar &);
};
}
#include "ThePEG/Utilities/ClassTraits.h"
namespace ThePEG {
/** @cond TRAITSPECIALIZATIONS */
/** This template specialization informs ThePEG about the
* base classes of MatchboxMEPP2llbar. */
template <>
struct BaseClassTrait<Herwig::MatchboxMEPP2llbar,1> {
/** Typedef of the first base class of MatchboxMEPP2llbar. */
typedef Herwig::MatchboxMEBase NthBase;
};
/** This template specialization informs ThePEG about the name of
* the MatchboxMEPP2llbar class and the shared object where it is defined. */
template <>
struct ClassTraits<Herwig::MatchboxMEPP2llbar>
: public ClassTraitsBase<Herwig::MatchboxMEPP2llbar> {
/** Return a platform-independent class name */
static string className() { return "Herwig::MatchboxMEPP2llbar"; }
/**
* The name of a file containing the dynamic library where the class
* MatchboxMEPP2llbar is implemented. It may also include several, space-separated,
* libraries if the class MatchboxMEPP2llbar depends on other classes (base classes
* excepted). In this case the listed libraries will be dynamically
* linked in the order they are specified.
*/
static string library() { return "HwMatchbox.so"; }
};
/** @endcond */
}
#endif /* HERWIG_MatchboxMEPP2llbar_H */

File Metadata

Mime Type
text/x-c++
Expires
Wed, May 14, 10:30 AM (1 d, 14 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5111159
Default Alt Text
MatchboxMEPP2llbar.h (7 KB)

Event Timeline