ThePEG  1.8.0
Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
ThePEG::MEGroup Class Reference

The MEGroup class represents a 'head' matrix element in association with a group of dependent matrix elements. More...

#include <MEGroup.h>

Inheritance diagram for ThePEG::MEGroup:
Inheritance graph
[legend]

List of all members.

Public Member Functions

tMEPtr head () const
 Return the head matrix element.
const MEVectordependent () const
 Visit the dependent matrix elements.
void head (tMEPtr me)
 Set the head matrix element.
MEVectordependent ()
 Access the dependent matrix elements.
int dependentOffset (tMEPtr dep) const
 Return the random number offset to access the random numbers provided for the given matrix element to generate dependent kinematics.
StdDependentXCombPtr makeDependentXComb (tStdXCombPtr xcHead, const cPDVector &proc, tMEPtr depME, const PartonPairVec &allPBins) const
 Create an StdDependentXComb object to be used for the given process steered bythe head object and dependent matrix element.
virtual bool subProcessGroups () const
 Return true, if SubProcessGroups should be setup from this MEGroup.
virtual bool mcSumDependent () const
 Return true for MC summation of dependent matrix elements, if feasible.
tStdDependentXCombPtr lastDependentXComb () const
 Return a pointer to the dependent xcomb object selected if mcSumDependent returned true.
virtual bool uniformAdditional () const =0
 Return true, if the same additional random numbers should be presented to any of the dependent matrix elements.
virtual MEBase::DiagramVector dependentDiagrams (const cPDVector &proc, tMEPtr depME) const =0
 Given a process from the head matrix element, return a list of diagrams which should be considered for the given dependent matrix element.
Standard constructors and destructors.
 MEGroup ()
 Default constructor.
virtual ~MEGroup ()
 Destructor.
Virtual functions from MEBase.
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_{EM}$ in which this matrix element is given.
virtual double me2 () const
 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 Energy2 scale () const
 Return the scale associated with the phase space point provided by the last call to setKinematics().
virtual double alphaS () const
 Return the value of $\alpha_S$ associated with the phase space point provided by the last call to setKinematics().
virtual double alphaEM () const
 Return the value of $\alpha_EM$ associated with the phase space point provided by the last call to setKinematics().
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.
virtual void constructVertex (tSubProPtr sub)
 construct the spin information for the interaction
virtual int nDim () const
 The number of internal degreed of freedom used in the matrix element.
virtual bool generateKinematics (const double *r)
 Generate internal degrees of freedom given nDim() uniform random numbers in the interval ]0,1[.
virtual bool wantCMS () const
 Return true, if this matrix element expects the incoming partons in their center-of-mass system.
virtual CrossSection dSigHatDR () const
 Return the matrix element squared differential in the variables given by the last call to generateKinematics().
virtual bool haveX1X2 () const
 Return true, if this matrix element will generate momenta for the incoming partons itself.
virtual bool havePDFWeight1 () const
 Return true, if this matrix element provides the PDF weight for the first incoming parton itself.
virtual bool havePDFWeight2 () const
 Return true, if this matrix element provides the PDF weight for the second incoming parton itself.
virtual bool keepRandomNumbers () const
 Return true, if the XComb steering this matrix element should keep track of the random numbers used to generate the last phase space point.
virtual void generateSubCollision (SubProcess &sub)
 Comlete a SubProcess object using the internal degrees of freedom generated in the last generateKinematics() (and possible other degrees of freedom which was intergated over in dSigHatDR().
virtual void clearKinematics ()
 Clear the information previously provided by a call to setKinematics(...).
virtual void getDiagrams () const
 Add all possible diagrams with the add() function.
virtual bool noMirror () const
 Return true, if this matrix element does not want to make use of mirroring processes; in this case all possible partonic subprocesses with a fixed assignment of incoming particles need to be provided through the diagrams added with the add(...) method.
virtual Selector< const
ColourLines * > 
colourGeometries (tcDiagPtr diag) const
 Return a Selector with possible colour geometries for the selected diagram weighted by their relative probabilities.
virtual const ColourLinesselectColourGeometry (tcDiagPtr diag) const
 Select a ColpurLines geometry.
virtual Selector< DiagramIndexdiagrams (const DiagramVector &dv) 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 DiagramIndex diagram (const DiagramVector &dv) const
 Select a diagram.
virtual void setXComb (tStdXCombPtr)
 Set the XComb object to be used in the next call to generateKinematics() and dSigHatDR().
virtual void flushCaches ()
 Inform this matrix element that a new phase space point is about to be generated, so all caches should be flushed.
virtual void lastEventStatistics ()
 Collect information on the last evaluated phasespace point for verification or debugging purposes.
Functions used by the persistent I/O system.
void persistentOutput (PersistentOStream &os) const
 Function used to write out object persistently.
void persistentInput (PersistentIStream &is, int version)
 Function used to read in object persistently.
- Public Member Functions inherited from ThePEG::MEBase
 MEBase ()
 Default constructor.
virtual ~MEBase ()
 Destructor.
void setKinematics (tPPair in, const PVector &out)
 Set the typed and momenta of the incoming and outgoing partons to be used in subsequent calls to me() and colourGeometries().
virtual bool headCuts () const
 If this is a dependent matrix element in a ME group, return true, if cuts should be inherited from the head matrix element, i.e.
virtual bool ignoreCuts () const
 If this is a dependent matrix element in a ME group, return true, if cuts should be ignored.
virtual bool apply () const
 If this is a dependent matrix element in a ME group, return true, if it applies to the process set in lastXComb()
const DiagramVectordiagrams () const
 Return all possible diagrams.
bool reweighted () const
 Return true if this matrix element has associated (p)reWeight objects assigned.
double reWeight () const
 With the information previously supplied with the setKinematics(...) methods, return the combined effects of the reweighters.
double preWeight () const
 With the information previously supplied with the setKinematics(...) methods, return the comined effects of the peweighters.
void addReweighter (tReweightPtr rw)
 Add objects to the list of reweighters.
void addPreweighter (tReweightPtr rw)
 Add objects to the list of preweighters.
Ptr< Amplitude >::pointer amplitude () const
 Return the amplitude associated with this matrix element.
void amplitude (Ptr< Amplitude >::pointer amp)
 Set the amplitude associated with this matrix element.
Energy2 sHat () const
 Return the last set invariant mass squared.
double preweight () const
 Return the factor with which this matrix element was last pre-weighted.
const DVectormeInfo () const
 Retrieve information obtained in the calculation of the cross section to be used later when selecting diagrams and colour flow.
void meInfo (const DVector &info) const
 Save information obtained in the calculation of the cross section to be used later when selecting diagrams and colour flow.
int maxMultCKKW () const
 If this matrix element is to be used together with others for CKKW reweighting and veto, this should give the multiplicity of outgoing particles in the highest multiplicity matrix element in the group.
int minMultCKKW () const
 If this matrix element is to be used together with others for CKKW reweighting and veto, this should give the multiplicity of outgoing particles in the lowest multiplicity matrix element in the group.
virtual void setVetoScales (tSubProPtr) const
 Set veto scales on the particles at the given SubProcess which has been generated using this matrix element.
- Public Member Functions inherited from ThePEG::LastXCombInfo< StandardXComb >
const StandardXComblastXComb () const
 Return a reference to the last selected XComb.
tXCPtr lastXCombPtr () const
 Return a pointer to the last selected XComb.
const StandardXComblastHeadXComb () const
 If the last selected XComb object belongs to a group of XComb's return a reference to the head XComb object for this group.
tXCPtr lastHeadXCombPtr () const
 If the last selected XComb object belongs to a group of XComb's return a pointer to the head XComb object for this group.
const EventHandlerlastEventHandler () const
 Return a reference to the currently used EventHandler.
tPExtrPtr lastExtractor () const
 A pointer to the currently used parton extractor.
PDFT pdf (tcPPtr parton) const
 Return the parton density used to extract the given parton.
const CutslastCuts () const
 A reference to the currently used kinematical cuts.
tCutsPtr lastCutsPtr () const
 A pointer to the currently used kinematical cuts.
const PPairlastParticles () const
 Return the pair of incoming parton instances.
Energy2 lastS () const
 The last generated total energy squared of the incoming particles.
const PPairlastPartons () const
 Return the pair of incoming parton instances.
Energy2 lastSHat () const
 The last used interval in total parton-parton energy squared.
double lastTau () const
 Return lastSHat()/lastS().
double lastY () const
 The generated rapidity of the hard scattering sub-system.
double lastP1 () const
 Log of one over the momentum fraction of the first incoming particle w.r.t.
double lastP2 () const
 Log of one over the momentum fraction of the second incoming particle w.r.t.
double lastL1 () const
 Log of one over the first incoming parton momentum fraction w.r.t.
double lastL2 () const
 Log of one over the second incoming parton momentum fraction w.r.t.
double lastX1 () const
 The first incoming parton momentum fraction w.r.t.
double lastX2 () const
 The second incoming parton momentum fraction w.r.t.
double lastE1 () const
 Return 1-lastX1() to highest possible precision for x $\rightarrow$ 1.
double lastE2 () const
 Return 1-lastX2() to highest possible precision for x $\rightarrow$ 1.
double lastFL1L2 () const
 The product of the parton density functions at the last generated phase-space point.
Energy2 lastScale () const
 The chosen scale of the hard scattering.
double lastAlphaS () const
 Get the $\alpha_S$ used in the hard scattering.
double lastAlphaEM () const
 Get the $\alpha_{EM}$ used in the hard scattering.
const vector< Lorentz5Momentum > & meMomenta () const
 Return the momenta of the incoming and outgoing partons to be used by the matrix element object, in the order specified by the TreeDiagram objects given by the matrix element.
double lastME2 () const
 Return the matrix element squared as calculated for the last phase space point.
double jacobian () const
 Get the last jacobian obtained when generating the kinematics for the call to dSigHatDR.
CrossSection lastMECrossSection () const
 Return the partonic cross section as calculated for the last phase space point.
double lastMEPDFWeight () const
 Return the PDF weight as calculated for the last phase space point, if the matrix element does supply PDF weights.
tSubProPtr subProcess () const
 Return the SubProcess object corresponding to the last generated sub-process.
const cPDVectormePartonData () const
 Return the incoming and outgoing parton types to be used by the matrix element object, in the order specified by the TreeDiagram objects given by the matrix element.

Static Public Member Functions

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

Protected Member Functions

void lastDependentXComb (tStdDependentXCombPtr xc)
 Set a pointer to the dependent xcomb object selected if mcSumDependent returned true.
Standard Interfaced functions.
virtual void doinit ()
 Initialize this object after the setup phase before saving an EventGenerator to disk.
virtual void rebind (const TranslationMap &trans)
 Rebind pointer to other Interfaced objects.
virtual IVector getReferences ()
 Return a vector of all pointers to Interfaced objects used in this object.
- Protected Member Functions inherited from ThePEG::MEBase
void add (DiagPtr dp) const
 To be used by sub classes in the getDiagrams() method to add included diagrams.
vector< Lorentz5Momentum > & meMomenta ()
 Access the momenta set by the last call to generateKinematics().
void lastME2 (double v) const
 Set the matrix element squared as calculated for the last phase space point.
void lastMECrossSection (CrossSection v) const
 Set the partonic cross section as calculated for the last phase space point.
void lastMEPDFWeight (double v) const
 Set the PDF weight as calculated for the last phase space point, if the matrix element does supply PDF weights.
void jacobian (double j)
 Set the last jacobian obtained when generating the kinematics for the call to dSigHatDR.
void use (tcMEPtr other)
 Initialize all member variables from another MEBase object.
void useDiagrams (tcMEPtr other) const
 Initialize the diagrams from another MEBase object.

Private Member Functions

MEGroupoperator= (const MEGroup &)
 Private and non-existent assignment operator.

Private Attributes

MEPtr theHead
 The head matrix element.
MEVector theDependent
 The dependent matrix elements.
map< tMEPtr, int > theNDimMap
 Offsets to access additional random numbers required by the dependent matrix elements.
int theNDim
 The total number of random numbers required.
tStdDependentXCombPtr theLastDependentXComb
 A pointer to the dependent xcomb object selected if mcSumDependent returned true.

Static Private Attributes

static
AbstractClassDescription
< MEGroup
initMEGroup
 Describe a class with persistent data.

Additional Inherited Members

- Public Types inherited from ThePEG::MEBase
typedef vector< DiagPtrDiagramVector
 A vector of pointers to DiagramBase objects.
typedef DiagramVector::size_type DiagramIndex
 The size_type used in the DiagramVector.
typedef vector< ReweightPtrReweightVector
 A vector of pointers to ReweightBase objects.
- Public Types inherited from ThePEG::LastXCombInfo< StandardXComb >
typedef ThePEG::Ptr
< StandardXComb >::pointer 
XCPtr
 Alias for a reference counted pointer to XC .
typedef ThePEG::Ptr
< StandardXComb >
::const_pointer 
cXCPtr
 Alias for a reference counted pointer to a const XC .
typedef ThePEG::Ptr
< StandardXComb >
::transient_pointer 
tXCPtr
 Alias for a transient pointer to XC .
typedef ThePEG::Ptr
< StandardXComb >
::transient_const_pointer 
tcXCPtr
 Alias for a transient pointer to a const XC .

Detailed Description

The MEGroup class represents a 'head' matrix element in association with a group of dependent matrix elements.

It basically acts as a wrapper around its head matrix element however supplying additional information to the corresponding StdXCombGroup object.

See also:
StdXCombGroup
StdDependentXComb

Definition at line 32 of file MEGroup.h.


Member Function Documentation

virtual double ThePEG::MEGroup::alphaEM ( ) const
inlinevirtual

Return the value of $\alpha_EM$ associated with the phase space point provided by the last call to setKinematics().

This versions returns SM().alphaEM(scale()).

Reimplemented from ThePEG::MEBase.

Definition at line 90 of file MEGroup.h.

References head().

virtual double ThePEG::MEGroup::alphaS ( ) const
inlinevirtual

Return the value of $\alpha_S$ associated with the phase space point provided by the last call to setKinematics().

This versions returns SM().alphaS(scale()).

Reimplemented from ThePEG::MEBase.

Definition at line 83 of file MEGroup.h.

References head().

virtual DiagramIndex ThePEG::MEGroup::diagram ( const DiagramVector dv) const
inlinevirtual

Select a diagram.

Default version uses diagrams(const DiagramVector &) to select a diagram according to the weights. This is the only method used that should be outside of MEBase.

Reimplemented from ThePEG::MEBase.

Definition at line 226 of file MEGroup.h.

References head().

virtual void ThePEG::MEGroup::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::InterfacedBase.

virtual bool ThePEG::MEGroup::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. The return value should be true of the generation succeeded. If so the generated momenta should be stored in the meMomenta() vector.

Implements ThePEG::MEBase.

virtual void ThePEG::MEGroup::generateSubCollision ( SubProcess sub)
inlinevirtual

Comlete a SubProcess object using the internal degrees of freedom generated in the last generateKinematics() (and possible other degrees of freedom which was intergated over in dSigHatDR().

This default version does nothing. Will be made purely virtual in the future.

Reimplemented from ThePEG::MEBase.

Definition at line 170 of file MEGroup.h.

References head().

virtual IVector ThePEG::MEGroup::getReferences ( )
protectedvirtual

Return a vector of all pointers to Interfaced objects used in this object.

Returns:
a vector of pointers.

Reimplemented from ThePEG::InterfacedBase.

virtual bool ThePEG::MEGroup::haveX1X2 ( ) const
inlinevirtual

Return true, if this matrix element will generate momenta for the incoming partons itself.

The matrix element is required to store the incoming parton momenta in meMomenta()[0,1]. No mapping in tau and y is performed by the PartonExtractor object, if a derived class returns true here. The phase space jacobian is to include a factor 1/(x1 x2).

Reimplemented from ThePEG::MEBase.

Definition at line 142 of file MEGroup.h.

References head().

virtual void ThePEG::MEGroup::lastEventStatistics ( )
inlinevirtual

Collect information on the last evaluated phasespace point for verification or debugging purposes.

This only called, if the StdXCombGroup did accumulate a non-zero cross section from this ME group.

Definition at line 250 of file MEGroup.h.

virtual int ThePEG::MEGroup::nDim ( ) const
inlinevirtual

The number of internal degreed of freedom used in the matrix element.

This default version returns 0;

Reimplemented from ThePEG::MEBase.

Definition at line 110 of file MEGroup.h.

References theNDim.

virtual unsigned int ThePEG::MEGroup::orderInAlphaEW ( ) const
inlinevirtual

Return the order in $\alpha_{EM}$ in which this matrix element is given.

Returns 0.

Implements ThePEG::MEBase.

Definition at line 63 of file MEGroup.h.

References head().

void ThePEG::MEGroup::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.

Reimplemented from ThePEG::MEBase.

void ThePEG::MEGroup::persistentOutput ( PersistentOStream os) const

Function used to write out object persistently.

Parameters:
osthe persistent output stream written to.

Reimplemented from ThePEG::MEBase.

virtual void ThePEG::MEGroup::rebind ( const TranslationMap trans)
protectedvirtual

Rebind pointer to other Interfaced objects.

Called in the setup phase after all objects used in an EventGenerator has been cloned so that the pointers will refer to the cloned objects afterwards.

Parameters:
transa TranslationMap relating the original objects to their respective clones.
Exceptions:
RebindExceptionif no cloned object was found for a given pointer.

Reimplemented from ThePEG::InterfacedBase.

virtual const ColourLines& ThePEG::MEGroup::selectColourGeometry ( tcDiagPtr  diag) const
inlinevirtual

Select a ColpurLines geometry.

The default version returns a colour geometry selected among the ones returned from colourGeometries(tcDiagPtr).

Reimplemented from ThePEG::MEBase.

Definition at line 208 of file MEGroup.h.

References head().

virtual void ThePEG::MEGroup::setKinematics ( )
virtual

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.

If the function is overridden in a sub class the new function must call the base class one first.

Reimplemented from ThePEG::MEBase.

virtual bool ThePEG::MEGroup::subProcessGroups ( ) const
inlinevirtual

Return true, if SubProcessGroups should be setup from this MEGroup.

If not, a single SubProcess is constructed from the data provided by the head matrix element.

Definition at line 298 of file MEGroup.h.


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