PYTHIA  8.303
Public Member Functions | Public Attributes | List of all members
Dire_fsr_qcd_Q2QGG Class Reference
Inheritance diagram for Dire_fsr_qcd_Q2QGG:
DireSplittingQCD DireSplitting

Public Member Functions

 Dire_fsr_qcd_Q2QGG (string idIn, int softRS, Settings *settings, ParticleData *particleData, Rndm *rndm, BeamParticle *beamA, BeamParticle *beamB, CoupSM *coupSM, Info *info, DireInfo *direInfo)
bool canRadiate (const Event &, pair< int, int >, unordered_map< string, bool >=unordered_map< string, bool >(), Settings *=NULL, PartonSystems *=NULL, BeamParticle *=NULL)
 Return true if this kernel should partake in the evolution.
bool canRadiate (const Event &, int iRadBef, int iRecBef, Settings *=NULL, PartonSystems *=NULL, BeamParticle *=NULL)
vector< int > radAndEmt (int idDaughter, int)
int nEmissions ()
int kinMap ()
bool canUseForBranching ()
vector< pair< int, int > > radAndEmtCols (int iRad, int colType, Event state)
pair< int, int > radBefCols (int colRadAfter, int acolRadAfter, int colEmtAfter, int acolEmtAfter)
 Return colours of recombined radiator (before splitting!)
int radBefID (int idRadAfter, int idEmtAfter)
 Return id of recombined radiator (before splitting!)
double gaugeFactor (int=0, int=0)
 Return color factor for splitting.
double symmetryFactor (int=0, int=0)
 Return symmetry factor for splitting.
double zSplit (double zMinAbs, double zMaxAbs, double m2dip)
 Pick z for new splitting. More...
double overestimateInt (double zMinAbs, double zMaxAbs, double pT2Old, double m2dip, int order=-1)
 New overestimates, z-integrated versions. More...
double overestimateDiff (double z, double m2dip, int order=-1)
 Return kernel for new splitting. More...
bool calc (const Event &state=Event(), int order=-1)
 Functions to calculate the kernel from SplitInfo information. More...
double counterTerm (double si1, double si2, double sj1, double sj2, double sij, double s12)
bool allow_sai_endpoint_for_kinematics ()
 Treatment of additional virtual corrections.
bool allow_xa_endpoint_for_kinematics ()
void try_sai_endpoint ()
 Functions to set if kernel should contribute to a kinematical endpoint.
void try_xa_endpoint ()
bool is_sai_endpoint ()
 Return endpoint information.
bool is_xa_endpoint ()
- Public Member Functions inherited from DireSplittingQCD
 DireSplittingQCD (string idIn, int softRS, Settings *settings, ParticleData *particleData, Rndm *rndm, BeamParticle *beamA, BeamParticle *beamB, CoupSM *coupSM, Info *info, DireInfo *direInfo)
 Constructor and destructor.
void init ()
double getNF (double pT2)
double GammaQCD2 (double NF=5.)
double GammaQCD3 (double NF=5.)
double betaQCD0 (double NF=5.)
double betaQCD1 (double NF=5.)
double betaQCD2 (double NF=5.)
double as2Pi (double pT2, int orderNow=-1, double renormMultFacNow=-1.)
double softRescaleInt (int order)
double softRescaleDiff (int order, double pT2, double renormMultFacNow=-1.)
double beta0Endpoint (int order, double m2dip, double pT2, double z, double renormMultFacNow=-1.)
double polevl (double x, double *coef, int N)
 Helper function to calculate dilogarithm.
double DiLog (double x)
 Function to calculate dilogarithm.
vector< int > sharedColor (const Event &event, int iRad, int iRec)
bool hasSharedColor (const Event &event, int iRad, int iRec)
int findCol (int col, vector< int > iExc, const Event &, int type)
bool useFastFunctions ()
virtual bool isPartial ()
virtual int couplingType (int, int)
virtual double coupling (double z, double pT2, double m2dip, double renormMultFacNow=-1., pair< int, bool > radBef=pair< int, bool >(), pair< int, bool > recBef=pair< int, bool >())
virtual double couplingScale2 (double z, double pT2, double m2dip, pair< int, bool > radBef, pair< int, bool > recBef)
virtual double getJacobian (const Event &=Event(), PartonSystems *partonSystems=0)
virtual unordered_map< string, double > getPhasespaceVars (const Event &=Event(), PartonSystems *=0)
double sCoef (int powz)
double sExp (int powz)
double kCoef (int powz)
double kExp (int powz)
double cCoef (int powz)
double cExp (int powz)
double fCoef ()
bool hasMECBef (const Event &state, double pT2)
bool hasMECAft (const Event &state, double pT2)
- Public Member Functions inherited from DireSplitting
 DireSplitting ()
 Constructor and destructor.
 DireSplitting (string idIn, int softRS, Settings *settings, ParticleData *particleData, Rndm *rndm, BeamParticle *beamA, BeamParticle *beamB, CoupSM *coupSMPtrIn, Info *infoPtrIn, DireInfo *direInfo)
void init ()
 The Splitting class.
bool is (ulong pattern)
string name ()
virtual bool aboveCutoff (double, const Particle &, const Particle &, int, PartonSystems *=NULL)
 Discard below the cut-off for the splitting.
virtual int motherID (int)
 Return id of mother after splitting.
virtual int sisterID (int)
 Return id of emission.
virtual bool swapRadEmt ()
virtual bool isSymmetric (const Particle *=NULL, const Particle *=NULL)
virtual vector< int > recPositions (const Event &, int, int)
virtual double getKernel (string="")
 Functions to store and retrieve all the variants of the kernel.
virtual unordered_map< string, double > getKernelVals ()
virtual void clearKernels ()
void setTimesPtr (shared_ptr< DireTimes > fsrIn)
void setTimesDecPtr (shared_ptr< DireTimes > fsrIn)
void setSpacePtr (shared_ptr< DireSpace > isrIn)
virtual bool allow_z_endpoint_for_kinematics ()
 Treatment of additional virtual corrections.
virtual bool allow_pT2_endpoint_for_kinematics ()
virtual void try_z_endpoint ()
 Functions to set if kernel should contribute to a kinematical endpoint.
virtual void try_pT2_endpoint ()
virtual bool is_z_endpoint ()
 Return endpoint information.
virtual bool is_pT2_endpoint ()
virtual double tdire_ff (double, double t, double)
 Functions to calculate Dire variables from the evolution variables.
virtual double zdire_ff (double z, double, double)
virtual double tdire_fi (double, double t, double)
virtual double zdire_fi (double z, double, double)
virtual double tdire_if (double, double t, double)
virtual double zdire_if (double z, double, double)
virtual double tdire_ii (double, double t, double)
virtual double zdire_ii (double z, double, double)
virtual void storeOverhead (double pT2, double x, int radid, int nf, double val)
virtual double overhead (double pT2, int idd, int nf)

Public Attributes

bool is_sai_endpoint_save
- Public Attributes inherited from DireSplittingQCD
double CA
double TR
double CF
double pTmin
double pT2minVariations
int NF_qcd_fsr
int orderSave
bool usePDFalphas
bool doVariations
bool doCorrelations
bool doMECs
double alphaSorder
double alphaS2pi
AlphaStrong alphaS
int asSchemeISR
int asSchemeFSR
bool useBackboneGluons
bool doGeneralizedKernel
- Public Attributes inherited from DireSplitting
double renormMultFac
string id
int correctionOrder
bool is_qcd
 Some short-cuts and string hashes to help avoid string comparisons.
bool is_qed
bool is_ewk
bool is_fsr
bool is_isr
bool is_dire
ulong nameHash
unordered_map< string, double > kernelVals
DireSplitInfo splitInfo
shared_ptr< DireSpaceisr
shared_ptr< DireTimesfsr
shared_ptr< DireTimesfsrDec
multimap< double, OverheadInfooverhead_map

Additional Inherited Members

- Static Public Attributes inherited from DireSplittingQCD
static const double SMALL_TEVOL = 2.0
 The SplittingQCD class.

Member Function Documentation

bool calc ( const Event state = Event(),
int  order = -1 

Functions to calculate the kernel from SplitInfo information.

Return kernel for new splitting.

Dummy statement to avoid compiler warnings.

Read all splitting variables.

Calculate argument of alphaS.

Do nothing without other NLO kernels!

Generate state after branching to extract momenta.

Get invariants.

Use only massless for now!

Get invariants.


Gluon contribution to cusp terms.

Additional contribution from si1 = 0. (might want separate endpoint eventually).

Multiply with LO kernel.


Do nothing below PS cut-off on kT of intermediate gluon, since never produced by PS, hence no underlying Born to correct.

Multiplicative weight factor to force first eikonal to current outgoing momenta.


Unsubtracted kernel.

Subtractions for all possible histories.

Avoid numerical issues if all subtractions are active (-> integrand should be exactly zero)

Overall color factor.

From xa integration volume.

Multiply by 2 since we randomly chose endpoint or fully differential.

Weight of sai-selection.

Remember that this might be an endpoint with vanishing sai.

Insert value of kernel into kernel list.

Create muR-variations.

Multiply with z to project out part where emitted gluon pair is soft, and quark is identified.

Store higher order correction separately.

Store kernel values.

Reimplemented from DireSplitting.

double counterTerm ( double  si1,
double  si2,
double  sj1,
double  sj2,
double  sij,
double  s12 

Counter-term is zero in unordered phase space.

Multiplicative weight factor to force first eikonal to current outgoing momenta.

Quark kernel.

Gluon kernel.

If this is an ordered region, only the subleading color dipoles remain.

Full counter-term


int kinMap ( )

Function to return an identifier for the phase space mapping that is used for setting up this splitting. return values: 1 –> Default Dire mapping. 2 –> Dire 1->3 mapping.

Reimplemented from DireSplitting.

double overestimateDiff ( double  z,
double  m2dip,
int  order = -1 

Return kernel for new splitting.

Return overestimate for new splitting.

Overestimate by soft + 1/(z+kappa2)

Reimplemented from DireSplitting.

double overestimateInt ( double  zMinAbs,
double  zMaxAbs,
double  pT2Old,
double  m2dip,
int  order = -1 

New overestimates, z-integrated versions.

Overestimate by soft + 1/(z+kappa2)

Reimplemented from DireSplitting.

vector<int> radAndEmt ( int  ,

Return a pair of ids for the radiator and emission after the splitting.

Reimplemented from DireSplittingQCD.

vector< pair< int, int > > radAndEmtCols ( int  iRad,
int  colType,
Event  state 

Also remember colors for "intermediate" particles in 1->3 splitting.

Reimplemented from DireSplitting.

double zSplit ( double  zMinAbs,
double  zMaxAbs,
double  m2dip 

Pick z for new splitting.

Pick according to soft + 1/(z+kappa2)

Reimplemented from DireSplitting.

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