herwig is hosted by Hepforge, IPPP Durham
Herwig  7.2.1
Herwig::MEPP2Higgs Class Reference

The MEPP2Higgs class implements the matrix element for the process pp->Higgs with different Higgs shape prescriptions (see details in hep-ph/9505211) and the NLL corrected Higgs width (see details in the FORTRAN HERWIG manual). More...

#include <MEPP2Higgs.h>

Inheritance diagram for Herwig::MEPP2Higgs:

Public Member Functions

 MEPP2Higgs ()
 The default constructor.
 
virtual CrossSection dSigHatDR () const
 Return the matrix element for the kinematical configuation previously provided by the last call to setKinematics(). More...
 
virtual void setKinematics ()
 Set the typed and momenta of the incoming and outgoing partons to be used in subsequent calls to me() and colourGeometries() according to the associated XComb object.
 
Member functions for the generation of hard QCD radiation
virtual POWHEGType hasPOWHEGCorrection ()
 Has a POWHEG style correction.
 
virtual bool hasMECorrection ()
 Has an old fashioned ME correction.
 
virtual void initializeMECorrection (RealEmissionProcessPtr, double &, double &)
 Initialize the ME correction.
 
virtual RealEmissionProcessPtr applyHardMatrixElementCorrection (RealEmissionProcessPtr)
 Apply the hard matrix element correction to a given hard process or decay.
 
virtual bool softMatrixElementVeto (PPtr parent, PPtr progenitor, const bool &fs, const Energy &highestpT, const vector< tcPDPtr > &ids, const double &z, const Energy &scale, const Energy &pT)
 Apply the soft matrix element correction. More...
 
virtual RealEmissionProcessPtr generateHardest (RealEmissionProcessPtr, ShowerInteraction)
 Apply the POWHEG style correction.
 
Virtual functions required by the MEBase class.
virtual unsigned int orderInAlphaS () const
 Return the order in $\alpha_S$ in which this matrix element is given.
 
virtual unsigned int orderInAlphaEW () const
 Return the order in $\alpha_{EW}$ in which this matrix element is given.
 
virtual double me2 () const
 The matrix element for the kinematical configuration previously provided by the last call to setKinematics(), suitably scaled by sHat() to give a dimension-less number. More...
 
virtual bool generateKinematics (const double *r)
 Generate internal degrees of freedom given nDim() uniform random numbers in the interval $ ]0,1[ $. More...
 
virtual Energy2 scale () const
 Return the scale associated with the last set phase space point.
 
virtual int nDim () const
 The number of internal degrees of freedom used in the matrix element.
 
virtual void getDiagrams () const
 Add all possible diagrams with the add() function.
 
virtual Selector< DiagramIndexdiagrams (const DiagramVector &dv) const
 Get diagram selector. More...
 
virtual Selector< const ColourLines * > colourGeometries (tcDiagPtr diag) const
 Return a Selector with possible colour geometries for the selected diagram weighted by their relative probabilities. More...
 
virtual void constructVertex (tSubProPtr)
 Construct the vertex of spin correlations.
 
Functions used by the persistent I/O system.
void persistentOutput (PersistentOStream &os) const
 Function used to write out object persistently. More...
 
void persistentInput (PersistentIStream &is, int version)
 Function used to read in object persistently. More...
 
- Public Member Functions inherited from Herwig::HwMEBase
 HwMEBase ()
 Default constructor.
 
void persistentOutput (PersistentOStream &os) const
 Function used to write out object persistently. More...
 
void persistentInput (PersistentIStream &is, int version)
 Function used to read in object persistently. More...
 
- Public Member Functions inherited from ThePEG::MEBase
virtual double alphaS () const
 
virtual double alphaEM () const
 
void setKinematics (tPPair in, const PVector &out)
 
virtual void constructVertex (tSubProPtr sub, const ColourLines *cl)
 
virtual bool wantCMS () const
 
virtual bool headCuts () const
 
virtual bool ignoreCuts () const
 
virtual bool apply () const
 
virtual map< string, double > generateOptionalWeights ()
 
virtual bool haveX1X2 () const
 
virtual bool havePDFWeight1 () const
 
virtual bool havePDFWeight2 () const
 
virtual bool keepRandomNumbers () const
 
virtual void generateSubCollision (SubProcess &)
 
virtual void clearKinematics ()
 
virtual bool noMirror () const
 
const DiagramVectordiagrams () const
 
virtual const ColourLinesselectColourGeometry (tcDiagPtr diag) const
 
virtual DiagramIndex diagram (const DiagramVector &) const
 
bool reweighted () const
 
double reWeight () const
 
double preWeight () const
 
void addReweighter (tReweightPtr rw)
 
void addPreweighter (tReweightPtr rw)
 
Ptr< Amplitude >::pointer amplitude () const
 
void amplitude (Ptr< Amplitude >::pointer amp)
 
Energy2 sHat () const
 
double preweight () const
 
virtual void flushCaches ()
 
virtual StdXCombPtr makeXComb (Energy newMaxEnergy, const cPDPair &inc, tEHPtr newEventHandler, tSubHdlPtr newSubProcessHandler, tPExtrPtr newExtractor, tCascHdlPtr newCKKW, const PBPair &newPartonBins, tCutsPtr newCuts, const DiagramVector &newDiagrams, bool mir, const PartonPairVec &allPBins, tStdXCombPtr newHead=tStdXCombPtr(), tMEPtr newME=tMEPtr())
 
virtual StdXCombPtr makeXComb (tStdXCombPtr newHead, const PBPair &newPartonBins, const DiagramVector &newDiagrams, tMEPtr newME=tMEPtr())
 
virtual void fillProjectors ()
 
virtual void setXComb (tStdXCombPtr)
 
const DVectormeInfo () const
 
void meInfo (const DVector &info) const
 
virtual int maxMultCKKW () const
 
virtual int minMultCKKW () const
 
virtual void maxMultCKKW (int mult)
 
virtual void minMultCKKW (int mult)
 
virtual void setVetoScales (tSubProPtr) const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
virtual double alphaS () const
 
virtual double alphaEM () const
 
void setKinematics (tPPair in, const PVector &out)
 
virtual void constructVertex (tSubProPtr sub, const ColourLines *cl)
 
virtual bool wantCMS () const
 
virtual bool headCuts () const
 
virtual bool ignoreCuts () const
 
virtual bool apply () const
 
virtual map< string, double > generateOptionalWeights ()
 
virtual bool haveX1X2 () const
 
virtual bool havePDFWeight1 () const
 
virtual bool havePDFWeight2 () const
 
virtual bool keepRandomNumbers () const
 
virtual void generateSubCollision (SubProcess &)
 
virtual void clearKinematics ()
 
virtual bool noMirror () const
 
const DiagramVectordiagrams () const
 
virtual const ColourLinesselectColourGeometry (tcDiagPtr diag) const
 
virtual DiagramIndex diagram (const DiagramVector &) const
 
bool reweighted () const
 
double reWeight () const
 
double preWeight () const
 
void addReweighter (tReweightPtr rw)
 
void addPreweighter (tReweightPtr rw)
 
Ptr< Amplitude >::pointer amplitude () const
 
void amplitude (Ptr< Amplitude >::pointer amp)
 
Energy2 sHat () const
 
double preweight () const
 
virtual void flushCaches ()
 
virtual StdXCombPtr makeXComb (Energy newMaxEnergy, const cPDPair &inc, tEHPtr newEventHandler, tSubHdlPtr newSubProcessHandler, tPExtrPtr newExtractor, tCascHdlPtr newCKKW, const PBPair &newPartonBins, tCutsPtr newCuts, const DiagramVector &newDiagrams, bool mir, const PartonPairVec &allPBins, tStdXCombPtr newHead=tStdXCombPtr(), tMEPtr newME=tMEPtr())
 
virtual StdXCombPtr makeXComb (tStdXCombPtr newHead, const PBPair &newPartonBins, const DiagramVector &newDiagrams, tMEPtr newME=tMEPtr())
 
virtual void fillProjectors ()
 
virtual void setXComb (tStdXCombPtr)
 
const DVectormeInfo () const
 
void meInfo (const DVector &info) const
 
virtual int maxMultCKKW () const
 
virtual int minMultCKKW () const
 
virtual void maxMultCKKW (int mult)
 
virtual void minMultCKKW (int mult)
 
virtual void setVetoScales (tSubProPtr) const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
- Public Member Functions inherited from HandlerBaseT< UseRandom >
double rnd () const
 
double rnd (double xu) const
 
double rnd (double xl, double xu) const
 
bool rndbool () const
 
bool rndbool (double p) const
 
bool rndbool (double p1, double p2) const
 
int rndsign (double p1, double p2, double p3) const
 
int rnd2 (double p0, double p1) const
 
int rnd3 (double p0, double p1, double p2) const
 
int rnd4 (double p0, double p1, double p2, double p3) const
 
long irnd (long xu=2) const
 
long irnd (long xl, long xu) const
 
const StandardModelBase & SM () const
 
tSMPtr standardModel () const
 
- Public Member Functions inherited from ThePEG::Interfaced
virtual bool defaultInit ()
 
PPtr getParticle (PID) const
 
PDPtr getParticleData (PID) const
 
bool used () const
 
void useMe () const
 
tEGPtr generator () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
PPtr getParticle (PID) const
 
PDPtr getParticleData (PID) const
 
bool used () const
 
void useMe () const
 
tEGPtr generator () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
- Public Member Functions inherited from ThePEG::InterfacedBase
string fullName () const
 
string name () const
 
string path () const
 
string comment () const
 
void setup (istream &is)
 
virtual void debugme () const
 
void update ()
 
void init ()
 
virtual bool preInitialize () const
 
void initrun ()
 
void finish ()
 
void touch ()
 
void reset ()
 
void clear ()
 
InitState state () const
 
bool locked () const
 
bool touched () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
void update ()
 
void init ()
 
virtual bool preInitialize () const
 
void initrun ()
 
void finish ()
 
void touch ()
 
void reset ()
 
void clear ()
 
InitState state () const
 
bool locked () const
 
bool touched () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
- Public Member Functions inherited from ThePEG::Base
void debug () const
 
- Public Member Functions inherited from ThePEG::Pointer::ReferenceCounted
CounterType referenceCount () const
 
- Public Member Functions inherited from ThePEG::Named
 Named (const string &newName=string())
 
 Named (const Named &)=default
 
const string & name () const
 
bool operator== (const Named &other) const
 
bool operator< (const Named &other) const
 
- Public Member Functions inherited from ThePEG::LastXCombInfo< StandardXComb >
const StandardXComblastXComb () const
 
const StandardXComblastXComb () const
 
tXCPtr lastXCombPtr () const
 
tXCPtr lastXCombPtr () const
 
const StandardXComblastHeadXComb () const
 
const StandardXComblastHeadXComb () const
 
tXCPtr lastHeadXCombPtr () const
 
tXCPtr lastHeadXCombPtr () const
 
const EventHandlerlastEventHandler () const
 
const EventHandlerlastEventHandler () const
 
tPExtrPtr lastExtractor () const
 
tPExtrPtr lastExtractor () const
 
PDFT pdf (tcPPtr parton) const
 
PDFT pdf (tcPPtr parton) const
 
const CutslastCuts () const
 
const CutslastCuts () const
 
tCutsPtr lastCutsPtr () const
 
tCutsPtr lastCutsPtr () const
 
const PPairlastParticles () const
 
const PPairlastParticles () const
 
Energy2 lastS () const
 
Energy2 lastS () const
 
const PPairlastPartons () const
 
const PPairlastPartons () const
 
Energy2 lastSHat () const
 
Energy2 lastSHat () const
 
double lastTau () const
 
double lastTau () const
 
double lastY () const
 
double lastY () const
 
double lastP1 () const
 
double lastP1 () const
 
double lastP2 () const
 
double lastP2 () const
 
double lastL1 () const
 
double lastL1 () const
 
double lastL2 () const
 
double lastL2 () const
 
double lastX1 () const
 
double lastX1 () const
 
double lastX2 () const
 
double lastX2 () const
 
double lastE1 () const
 
double lastE1 () const
 
double lastE2 () const
 
double lastE2 () const
 
double lastFL1L2 () const
 
double lastFL1L2 () const
 
Energy2 lastScale () const
 
Energy2 lastScale () const
 
double lastAlphaS () const
 
double lastAlphaS () const
 
double lastAlphaEM () const
 
double lastAlphaEM () const
 
const vector< Lorentz5Momentum > & meMomenta () const
 
const vector< Lorentz5Momentum > & meMomenta () const
 
double lastME2 () const
 
double lastME2 () const
 
double lastPreweight () const
 
double lastPreweight () const
 
double jacobian () const
 
double jacobian () const
 
CrossSection lastMECrossSection () const
 
CrossSection lastMECrossSection () const
 
double lastMEPDFWeight () const
 
double lastMEPDFWeight () const
 
double lastMECouplings () const
 
double lastMECouplings () const
 
tSubProPtr subProcess () const
 
tSubProPtr subProcess () const
 
const cPDVectormePartonData () const
 
const cPDVectormePartonData () const
 

Static Public Member Functions

static void Init ()
 The standard Init function used to initialize the interfaces. More...
 
- Static Public Member Functions inherited from Herwig::HwMEBase
static void Init ()
 The standard Init function used to initialize the interfaces. More...
 
- Static Public Member Functions inherited from ThePEG::MEBase
static void Init ()
 
- Static Public Member Functions inherited from ThePEG::HandlerBase
static void Init ()
 
- Static Public Member Functions inherited from ThePEG::Interfaced
static void Init ()
 
- Static Public Member Functions inherited from ThePEG::InterfacedBase
static void Init ()
 
- Static Public Member Functions inherited from ThePEG::Base
static void Init ()
 

Protected Member Functions

tPDPtr quarkFlavour (tcPDFPtr pdf, Energy2 scale, double x, tcBeamPtr beam, double &pdfweight, bool anti)
 Method to extract the PDF weight for quark/antiquark initiated processes and select the quark flavour.
 
bool applyHard (ParticleVector gluons, vector< tcBeamPtr > beams, PPtr higgs, unsigned int &iemit, unsigned int &itype, vector< Lorentz5Momentum > &pnew, pair< double, double > &xnew, tPDPtr &out)
 Return the momenta and type of hard matrix element correction. More...
 
bool getEvent (vector< Lorentz5Momentum > &pnew, int &emissiontype)
 generates the hardest emission (yj,p) More...
 
double getResult (int emis_type, Energy pt, double yj, tcPDPtr &outParton)
 Returns the matrix element for a given type of process, rapidity of the jet $y_j$ and transverse momentum $p_T$. More...
 
Clone Methods.
virtual IBPtr clone () const
 Make a simple clone of this object. More...
 
virtual IBPtr fullclone () const
 Make a clone of this object, possibly modifying the cloned object to make it sane. More...
 
Standard Interfaced functions.
virtual void doinit ()
 Initialize this object after the setup phase before saving an EventGenerator to disk. More...
 
virtual void dofinish ()
 Finalize this object. More...
 
Energy4 loME () const
 Members to calculate the real emission matrix elements. More...
 
Energy2 ggME (Energy2 s, Energy2 t, Energy2 u)
 The matrix element for $gg\to H g$.
 
Energy2 qgME (Energy2 s, Energy2 t, Energy2 u)
 The matrix element for $qg\to H q$.
 
Energy2 qbargME (Energy2 s, Energy2 t, Energy2 u)
 The matrix element for $qbarg\to H qbar$.
 
Complex B (Energy2 s, Energy2 mf2) const
 Members to calculate the functions for the loop diagrams. More...
 
complex< InvEnergy2 > C (Energy2 s, Energy2 mf2) const
 The $C(s)$ function of NBP339 (1990) 38-66. More...
 
complex< InvEnergy4 > D (Energy2 s, Energy2 t, Energy2 u, Energy2 mf2) const
 The $C(s)$ function of NBP339 (1990) 38-66. More...
 
Complex dIntegral (Energy2 a, Energy2 b, double y0) const
 The integral $\int\frac{dy}{y-y_0}\log(a-i\epsilon-b y(1-y))$ from NBP339 (1990) 38-66. More...
 
complex< Energy > me1 (Energy2 s, Energy2 t, Energy2 u, Energy2 mf2, unsigned int i, unsigned int j, unsigned int k, unsigned int i1, unsigned int j1, unsigned int k1) const
 The $M_{+++}$ matrix element of NBP339 (1990) 38-66. More...
 
complex< Energy > me2 (Energy2 s, Energy2 t, Energy2 u, Energy2 mf2) const
 The $M_{++-}$ matrix element of NBP339 (1990) 38-66. More...
 
Complex F (double x) const
 The $F(x)$ function for the leading-order result.
 
- Protected Member Functions inherited from Herwig::HwMEBase
void massOption (vector< unsigned int > iopt)
 Set the treatment of the outgoing masses. More...
 
bool generateMasses (vector< Energy > &masses, double &mjac, const double *r)
 Generate the masses of the particles.
 
virtual double getCosTheta (double cthmin, double cthmax, const double r)
 Used internally by generateKinematics, after calculating the limits on cos(theta).
 
Energy2 tHat () const
 Return the $\hat{t}$ of the last set phase space point.
 
Energy2 uHat () const
 Return the $\hat{u}$ of the last set phase space point.
 
double phi () const
 Return the azimuth angle of the last set phase space point.
 
void tHat (Energy2 e2)
 Set the $\hat{t}$ of the last set phase space point.
 
void uHat (Energy2 e2)
 Set the $\hat{u}$ of the last set phase space point.
 
void phi (double phi)
 Set the azimuth angle of the last set phase space point.
 
void rescalingOption (unsigned int iopt)
 Rescaled momenta for the helicity ME. More...
 
bool rescaleMomenta (const vector< Lorentz5Momentum > &, const cPDVector &)
 rescale the momenta for the computation of the helicity matrix element
 
const vector< Lorentz5Momentum > & rescaledMomenta () const
 Access to the rescaled momenta.
 
- Protected Member Functions inherited from ThePEG::MEBase
void add (DiagPtr dp) const
 
vector< Lorentz5Momentum > & meMomenta ()
 
void lastME2 (double v) const
 
void lastPreweight (double w) const
 
void lastMECrossSection (CrossSection v) const
 
void lastMEPDFWeight (double v) const
 
void lastMECouplings (double v) const
 
void jacobian (double j)
 
void use (tcMEPtr other)
 
void useDiagrams (tcMEPtr other) const
 
virtual void doinitrun ()
 
virtual void doinitrun ()
 
- Protected Member Functions inherited from ThePEG::Interfaced
void reporeg (IBPtr object, string name) const
 
bool setDefaultReference (PtrT &ptr, string classname, string objectname)
 
 Interfaced (const string &newName)
 
 Interfaced (const Interfaced &i)
 
void setGenerator (tEGPtr generator)
 
- Protected Member Functions inherited from ThePEG::InterfacedBase
 InterfacedBase (string newName)
 
 InterfacedBase (const InterfacedBase &i)
 
virtual void readSetup (istream &is)
 
virtual void doupdate ()
 
virtual IVector getReferences ()
 
virtual void rebind (const TranslationMap &)
 
virtual void readSetup (istream &is)
 
virtual void doupdate ()
 
virtual IVector getReferences ()
 
virtual void rebind (const TranslationMap &)
 
- Protected Member Functions inherited from ThePEG::Pointer::ReferenceCounted
 ReferenceCounted (const ReferenceCounted &)
 
ReferenceCountedoperator= (const ReferenceCounted &)
 
 ReferenceCounted (const ReferenceCounted &)
 
ReferenceCountedoperator= (const ReferenceCounted &)
 
- Protected Member Functions inherited from ThePEG::Named
const Namedoperator= (const Named &other)
 
const string & name (const string &newName)
 

Private Member Functions

MEPP2Higgsoperator= (const MEPP2Higgs &)=delete
 The assignment operator is private and must never be called. More...
 
double ggME (vector< VectorWaveFunction > g1, vector< VectorWaveFunction > g2, ScalarWaveFunction &, bool calc) const
 Members to return the matrix elements for the different subprocesses. More...
 
double qqME (vector< SpinorWaveFunction > &fin, vector< SpinorBarWaveFunction > &ain, ScalarWaveFunction &, bool calc) const
 Calculates the matrix element for the process q,qbar->h. More...
 

Private Attributes

unsigned int scaleopt_
 Selects a dynamic (sHat) or fixed factorization scale.
 
Energy mu_F_
 The value associated to the fixed factorization scale option.
 
unsigned int shapeOption_
 Defines the Higgs resonance shape.
 
unsigned int processOption_
 The processes to be included (GG->H and/or qq->H)
 
int minFlavour_
 Minimum flavour of incoming quarks.
 
int maxFlavour_
 Maximum flavour of incoming quarks.
 
ProductionMatrixElement me_
 Matrix element for spin correlations.
 
AbstractVVSVertexPtr HGGVertex_
 Pointer to the H-> 2 gluon vertex (used in gg->H)
 
AbstractFFSVertexPtr HFFVertex_
 Pointer to the fermion-fermion Higgs vertex (used in qq->H)
 
GenericMassGeneratorPtr hmass_
 The mass generator for the Higgs.
 
Energy mh_
 On-shell mass for the higgs.
 
Energy wh_
 On-shell width for the higgs.
 
ShowerAlphaPtr alpha_
 Pointer to the object calculating the strong coupling.
 
Energy2 mh2_
 Mass squared of Higgs.
 
double channelwgtA_
 Relative weight of the $qg$ to the $gg$ channel.
 
double channelwgtB_
 Relative weight for the $\bar{q}g$ to the $gg$ channel.
 
vector< double > channelWeights_
 Weights for the channels as a vector.
 
double ggPow_
 Power for the $\frac{{\rm d}\hat{s}}{\hat{s}^n}$ importance sampling of the $gg$ component.
 
double qgPow_
 Power for the $\frac{{\rm d}\hat{s}}{\hat{s}^n}$ importance sampling of the $qg$ and $\bar{q}g$ components.
 
double enhance_
 The enhancement factor for initial-state radiation.
 
unsigned int nover_
 Number of weights greater than 1.
 
unsigned int ntry_
 Number of attempts.
 
unsigned int ngen_
 Number which suceed.
 
double maxwgt_
 Maximum weight.
 
Energy minpT_
 The transverse momentum of the jet.
 
bool spinCorrelations_
 Whether of not to construct the vertex for spin correlations.
 
unsigned int minLoop_
 Stuff for the ME correction. More...
 
unsigned int maxLoop_
 Maximum flavour of quarks to include in the loops.
 
unsigned int massOption_
 Option for treatment of the fermion loops.
 
unsigned int mu_R_opt_
 Option for dynamic scale choice in alpha_S (0=mT,>0=pT)
 
unsigned int mu_F_opt_
 Option for dynamic scale choice in PDFs (0=mT,>0=pT)
 
Complex bi_ [5]
 Storage of the loop functions. More...
 
complex< InvEnergy2 > ci_ [8]
 C functions.
 
complex< InvEnergy4 > di_ [4]
 D functions.
 
double power_
 Constants for the sampling. More...
 
double pregg_
 The prefactor, $c$ for the $gg$ channel.
 
double preqg_
 The prefactor, $c$ for the $qg$ channel.
 
double pregqbar_
 The prefactor, $c$ for the $g\bar{q}$ channel.
 
vector< double > prefactor_
 The prefactors as a vector for easy use.
 
vector< tcBeamPtrbeams_
 Properties of the incoming particles. More...
 
vector< tcPDPtrpartons_
 Pointers to the ParticleDataObjects for the partons.
 
double yh_
 Properties of the boson and jets. More...
 
Energy mass_
 The mass of the Higgs boson.
 
double yj_
 the rapidity of the jet
 
Energy pt_
 The transverse momentum of the jet.
 
tcPDPtr out_
 The outgoing parton.
 

Static Private Attributes

static const complex< Energy2 > epsi_
 Small complex number to regularize some integrals.
 

Additional Inherited Members

- Public Types inherited from Herwig::HwMEBase
enum  POWHEGType { No, ISR, FSR, Both }
 Virtual members to be overridden by inheriting classes which implement hard corrections. More...
 
- Public Types inherited from ThePEG::MEBase
typedef vector< DiagPtr > DiagramVector
 
typedef DiagramVector::size_type DiagramIndex
 
typedef vector< ReweightPtrReweightVector
 
- Public Types inherited from ThePEG::InterfacedBase
enum  InitState
 
- Public Types inherited from ThePEG::Pointer::ReferenceCounted
typedef unsigned int CounterType
 
- Public Types inherited from ThePEG::LastXCombInfo< StandardXComb >
typedef ThePEG::Ptr< StandardXComb >::pointer XCPtr
 
typedef ThePEG::Ptr< StandardXComb >::const_pointer cXCPtr
 
typedef ThePEG::Ptr< StandardXComb >::transient_pointer tXCPtr
 
typedef ThePEG::Ptr< StandardXComb >::transient_const_pointer tcXCPtr
 
- Public Attributes inherited from ThePEG::InterfacedBase
 initializing
 
 uninitialized
 
 initialized
 
 runready
 
- Public Attributes inherited from ThePEG::Pointer::ReferenceCounted
const unsigned long uniqueId
 
- Static Protected Member Functions inherited from ThePEG::Interfaced
static void registerRepository (IBPtr)
 
static void registerRepository (IBPtr, string newName)
 
- Protected Attributes inherited from ThePEG::LastXCombInfo< StandardXComb >
XCPtr theLastXComb
 

Detailed Description

The MEPP2Higgs class implements the matrix element for the process pp->Higgs with different Higgs shape prescriptions (see details in hep-ph/9505211) and the NLL corrected Higgs width (see details in the FORTRAN HERWIG manual).

See also
The interfaces defined for MEPP2Higgs.

Definition at line 38 of file MEPP2Higgs.h.

Member Function Documentation

◆ applyHard()

bool Herwig::MEPP2Higgs::applyHard ( ParticleVector  gluons,
vector< tcBeamPtr beams,
PPtr  higgs,
unsigned int &  iemit,
unsigned int &  itype,
vector< Lorentz5Momentum > &  pnew,
pair< double, double > &  xnew,
tPDPtr out 
)
protected

Return the momenta and type of hard matrix element correction.

Parameters
gluonsThe original incoming particles.
beamsThe BeamParticleData objects
higgsThe original outgoing higgs
iemitWhether the first (0) or second (1) particle emitted the radiation
itypeThe type of radiated particle (0 is gluon, 1 is quark and 2 is antiquark)
pnewThe momenta of the new particles
xnewThe new values of the momentuym fractions
outThe ParticleData object for the outgoing parton
Returns
Whether or not the matrix element correction needs to be applied

◆ B()

Complex Herwig::MEPP2Higgs::B ( Energy2  s,
Energy2  mf2 
) const
protected

Members to calculate the functions for the loop diagrams.

The $B(s)$ function of NBP339 (1990) 38-66

Parameters
sThe scale
mf2The fermion mass squared.

◆ C()

complex<InvEnergy2> Herwig::MEPP2Higgs::C ( Energy2  s,
Energy2  mf2 
) const
protected

The $C(s)$ function of NBP339 (1990) 38-66.

Parameters
sThe scale
mf2The fermion mass squared.

◆ clone()

virtual IBPtr Herwig::MEPP2Higgs::clone ( ) const
inlineprotectedvirtual

Make a simple clone of this object.

Returns
a pointer to the new object.

Implements ThePEG::InterfacedBase.

Reimplemented in Herwig::MEPP2HiggsPowheg.

Definition at line 227 of file MEPP2Higgs.h.

◆ colourGeometries()

virtual Selector<const ColourLines *> Herwig::MEPP2Higgs::colourGeometries ( tcDiagPtr  diag) const
virtual

Return a Selector with possible colour geometries for the selected diagram weighted by their relative probabilities.

Parameters
diagthe diagram chosen.
Returns
the possible colour geometries weighted by their relative probabilities.

Implements ThePEG::MEBase.

◆ D()

complex<InvEnergy4> Herwig::MEPP2Higgs::D ( Energy2  s,
Energy2  t,
Energy2  u,
Energy2  mf2 
) const
protected

The $C(s)$ function of NBP339 (1990) 38-66.

Parameters
sThe $s$ invariant
tThe $t$ invariant
uThe $u$ invariant
mf2The fermion mass squared

◆ diagrams()

virtual Selector<DiagramIndex> Herwig::MEPP2Higgs::diagrams ( const DiagramVector dv) const
virtual

Get diagram selector.

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.

Parameters
dvthe diagrams to be weighted.
Returns
a Selector relating the given diagrams to their weights.

Reimplemented from ThePEG::MEBase.

◆ dIntegral()

Complex Herwig::MEPP2Higgs::dIntegral ( Energy2  a,
Energy2  b,
double  y0 
) const
protected

The integral $\int\frac{dy}{y-y_0}\log(a-i\epsilon-b y(1-y))$ from NBP339 (1990) 38-66.

Parameters
aThe parameter $a$.
bThe parameter $b$.
y0The parameter $y_0$.

◆ dofinish()

virtual void Herwig::MEPP2Higgs::dofinish ( )
protectedvirtual

Finalize this object.

Called in the run phase just after a run has ended. Used eg. to write out statistics.

Reimplemented from ThePEG::InterfacedBase.

◆ doinit()

virtual void Herwig::MEPP2Higgs::doinit ( )
protectedvirtual

Initialize this object after the setup phase before saving an EventGenerator to disk.

Exceptions
InitExceptionif object could not be initialized properly.

Reimplemented from ThePEG::MEBase.

◆ dSigHatDR()

virtual CrossSection Herwig::MEPP2Higgs::dSigHatDR ( ) const
virtual

Return the matrix element for the kinematical configuation previously provided by the last call to setKinematics().

Uses me().

Reimplemented from Herwig::HwMEBase.

◆ fullclone()

virtual IBPtr Herwig::MEPP2Higgs::fullclone ( ) const
inlineprotectedvirtual

Make a clone of this object, possibly modifying the cloned object to make it sane.

Returns
a pointer to the new object.

Reimplemented from ThePEG::InterfacedBase.

Reimplemented in Herwig::MEPP2HiggsPowheg.

Definition at line 233 of file MEPP2Higgs.h.

◆ generateKinematics()

virtual bool Herwig::MEPP2Higgs::generateKinematics ( const double *  r)
virtual

Generate internal degrees of freedom given nDim() uniform random numbers in the interval $ ]0,1[ $.

To help the phase space generator, the dSigHatDR should be a smooth function of these numbers, although this is not strictly necessary.

Parameters
ra pointer to the first of nDim() consecutive random numbers.
Returns
true if the generation succeeded, otherwise false.

Reimplemented from Herwig::HwMEBase.

Reimplemented in Herwig::MEPP2HiggsPowheg.

◆ getEvent()

bool Herwig::MEPP2Higgs::getEvent ( vector< Lorentz5Momentum > &  pnew,
int &  emissiontype 
)
protected

generates the hardest emission (yj,p)

Parameters
pnewThe momenta of the new particles
emissiontypeThe type of emission, as for getResult
Returns
Whether not an emission was generated

◆ getResult()

double Herwig::MEPP2Higgs::getResult ( int  emis_type,
Energy  pt,
double  yj,
tcPDPtr outParton 
)
protected

Returns the matrix element for a given type of process, rapidity of the jet $y_j$ and transverse momentum $p_T$.

Parameters
emis_typethe type of emission, (0 is $gg\to h^0g$, 1 is $qg\to h^0q$ and 2 is $g\bar{q}\to h^0\bar{q}$)
ptThe transverse momentum of the jet
yjThe rapidity of the jet
outPartonthe outgoing parton

◆ ggME()

double Herwig::MEPP2Higgs::ggME ( vector< VectorWaveFunction g1,
vector< VectorWaveFunction g2,
ScalarWaveFunction ,
bool  calc 
) const
private

Members to return the matrix elements for the different subprocesses.

Calculates the matrix element for the process g,g->h (via quark loops)

Parameters
g1a vector of wave functions of the first incoming gluon
g2a vector of wave functions of the second incoming gluon
calcWhether or not to calculate the matrix element for spin correlations
Returns
the amlitude value.

◆ Init()

static void Herwig::MEPP2Higgs::Init ( )
static

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.

◆ loME()

Energy4 Herwig::MEPP2Higgs::loME ( ) const
protected

Members to calculate the real emission matrix elements.

The leading-order matrix element for $gg\to H$

◆ me1()

complex<Energy> Herwig::MEPP2Higgs::me1 ( Energy2  s,
Energy2  t,
Energy2  u,
Energy2  mf2,
unsigned int  i,
unsigned int  j,
unsigned int  k,
unsigned int  i1,
unsigned int  j1,
unsigned int  k1 
) const
protected

The $M_{+++}$ matrix element of NBP339 (1990) 38-66.

Parameters
sThe $s$ invariant
tThe $t$ invariant
uThe $u$ invariant
mf2The fermion mass squared.
iWhich of the stored values to use for $D(u,t)$.
jWhich of the stored values to use for $D(u,s)$.
kWhich of the stored values to use for $D(s,t)$.
i1Which of the stored values to use for $C_1(s)$.
j1Which of the stored values to use for $C_1(t)$.
k1Which of the stored values to use for $C_1(u)$.

◆ me2() [1/2]

virtual double Herwig::MEPP2Higgs::me2 ( ) const
virtual

The matrix element for the kinematical configuration previously provided by the last call to setKinematics(), suitably scaled by sHat() to give a dimension-less number.

Returns
the matrix element scaled with sHat() to give a dimensionless number.

Implements ThePEG::MEBase.

Reimplemented in Herwig::MEPP2HiggsPowheg.

◆ me2() [2/2]

complex<Energy> Herwig::MEPP2Higgs::me2 ( Energy2  s,
Energy2  t,
Energy2  u,
Energy2  mf2 
) const
protected

The $M_{++-}$ matrix element of NBP339 (1990) 38-66.

Parameters
sThe $s$ invariant
tThe $t$ invariant
uThe $u$ invariant
mf2The fermion mass squared.

◆ operator=()

MEPP2Higgs& Herwig::MEPP2Higgs::operator= ( const MEPP2Higgs )
privatedelete

The assignment operator is private and must never be called.

In fact, it should not even be implemented.

◆ persistentInput()

void Herwig::MEPP2Higgs::persistentInput ( PersistentIStream is,
int  version 
)

Function used to read in object persistently.

Parameters
isthe persistent input stream read from.
versionthe version number of the object when written.

◆ persistentOutput()

void Herwig::MEPP2Higgs::persistentOutput ( PersistentOStream os) const

Function used to write out object persistently.

Parameters
osthe persistent output stream written to.

◆ qqME()

double Herwig::MEPP2Higgs::qqME ( vector< SpinorWaveFunction > &  fin,
vector< SpinorBarWaveFunction > &  ain,
ScalarWaveFunction ,
bool  calc 
) const
private

Calculates the matrix element for the process q,qbar->h.

Parameters
fina vector of quark spinors
aina vector of anti-quark spinors
calcWhether or not to calculate the matrix element for spin correlations
Returns
the amlitude value.

◆ softMatrixElementVeto()

virtual bool Herwig::MEPP2Higgs::softMatrixElementVeto ( PPtr  parent,
PPtr  progenitor,
const bool &  fs,
const Energy &  highestpT,
const vector< tcPDPtr > &  ids,
const double &  z,
const Energy &  scale,
const Energy &  pT 
)
virtual

Apply the soft matrix element correction.

Parameters
parentThe initial particle in the current branching
progenitorThe progenitor particle of the jet
fsWhether the emission is initial or final-state
highestpTThe highest pT so far in the shower
idsids of the particles produced in the branching
zThe momentum fraction of the branching
scalethe evolution scale of the branching
pTThe transverse momentum of the branching
Returns
If true the emission should be vetoed

Reimplemented from Herwig::HwMEBase.

Member Data Documentation

◆ beams_

vector<tcBeamPtr> Herwig::MEPP2Higgs::beams_
private

Properties of the incoming particles.

Pointers to the BeamParticleData objects

Definition at line 664 of file MEPP2Higgs.h.

◆ bi_

Complex Herwig::MEPP2Higgs::bi_[5]
mutableprivate

Storage of the loop functions.

B functions

Definition at line 545 of file MEPP2Higgs.h.

◆ minLoop_

unsigned int Herwig::MEPP2Higgs::minLoop_
private

Stuff for the ME correction.

Parameters for the evaluation of the loops for the matrix elements Minimum flavour of quarks to include in the loops

Definition at line 508 of file MEPP2Higgs.h.

◆ power_

double Herwig::MEPP2Higgs::power_
private

Constants for the sampling.

The distribution is assumed to have the form $\frac{c}{{\rm GeV}}\times\left(\frac{{\rm GeV}}{p_T}\right)^n$ The power, $n$, for the sampling

Definition at line 629 of file MEPP2Higgs.h.

◆ yh_

double Herwig::MEPP2Higgs::yh_
private

Properties of the boson and jets.

The rapidity of the Higgs boson

Definition at line 679 of file MEPP2Higgs.h.


The documentation for this class was generated from the following file: