ThePEG
1.8.0
|
The MEGroup class represents a 'head' matrix element in association with a group of dependent matrix elements. More...
#include <MEGroup.h>
Public Member Functions | |
tMEPtr | head () const |
Return the head matrix element. | |
const MEVector & | dependent () const |
Visit the dependent matrix elements. | |
void | head (tMEPtr me) |
Set the head matrix element. | |
MEVector & | dependent () |
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 ![]() | |
virtual unsigned int | orderInAlphaEW () const |
Return the order in ![]() | |
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 ![]() | |
virtual double | alphaEM () const |
Return the value of ![]() | |
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 ColourLines & | selectColourGeometry (tcDiagPtr diag) const |
Select a ColpurLines geometry. | |
virtual Selector< DiagramIndex > | diagrams (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. | |
![]() | |
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 DiagramVector & | diagrams () 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 DVector & | meInfo () 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. | |
![]() | |
const StandardXComb & | lastXComb () const |
Return a reference to the last selected XComb. | |
tXCPtr | lastXCombPtr () const |
Return a pointer to the last selected XComb. | |
const StandardXComb & | lastHeadXComb () 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 EventHandler & | lastEventHandler () 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 Cuts & | lastCuts () const |
A reference to the currently used kinematical cuts. | |
tCutsPtr | lastCutsPtr () const |
A pointer to the currently used kinematical cuts. | |
const PPair & | lastParticles () const |
Return the pair of incoming parton instances. | |
Energy2 | lastS () const |
The last generated total energy squared of the incoming particles. | |
const PPair & | lastPartons () 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 ![]() | |
double | lastE2 () const |
Return 1-lastX2() to highest possible precision for x ![]() | |
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 ![]() | |
double | lastAlphaEM () const |
Get the ![]() | |
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 cPDVector & | mePartonData () 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. | |
![]() | |
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 | |
MEGroup & | operator= (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 | |
![]() | |
typedef vector< DiagPtr > | DiagramVector |
A vector of pointers to DiagramBase objects. | |
typedef DiagramVector::size_type | DiagramIndex |
The size_type used in the DiagramVector. | |
typedef vector< ReweightPtr > | ReweightVector |
A vector of pointers to ReweightBase objects. | |
![]() | |
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 . |
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.
|
inlinevirtual |
Return the value of 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().
|
inlinevirtual |
Return the value of 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().
|
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().
|
protectedvirtual |
Initialize this object after the setup phase before saving an EventGenerator to disk.
InitException | if object could not be initialized properly. |
Reimplemented from ThePEG::InterfacedBase.
|
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.
|
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().
|
protectedvirtual |
Return a vector of all pointers to Interfaced objects used in this object.
Reimplemented from ThePEG::InterfacedBase.
|
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().
|
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.
|
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.
|
inlinevirtual |
Return the order in 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.
is | the persistent input stream read from. |
version | the 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.
os | the persistent output stream written to. |
Reimplemented from ThePEG::MEBase.
|
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.
trans | a TranslationMap relating the original objects to their respective clones. |
RebindException | if no cloned object was found for a given pointer. |
Reimplemented from ThePEG::InterfacedBase.
|
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 |
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.
|
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.