PYTHIA
8.312
|
Final-final electroweak antenna. More...
#include <VinciaEW.h>
Public Member Functions | |
virtual bool | init (Event &event, int iMotIn, int iRecIn, int iSysIn, vector< EWBranching > &branchings, Settings *settingsPtr) override |
Overridden virtual functions. More... | |
virtual double | generateTrial (double q2Start, double q2End, double alphaIn) override |
Generate a trial. More... | |
virtual bool | acceptTrial (Event &event) override |
Accept trial. More... | |
virtual void | updateEvent (Event &event) override |
Update an event. More... | |
Public Member Functions inherited from EWAntenna | |
EWAntenna () | |
Constructor and destructor. | |
void | print () |
Print, must be implemented by base classes. | |
void | initPtr (Info *infoPtrIn, VinciaCommon *vinComPtrIn, AlphaEM *alphaPtrIn, AmpCalculator *ampCalcPtrIn) |
Initialize pointers. | |
void | setVerbose (int verboseIn) |
Set verbosity level. | |
virtual void | updatePartonSystems (Event &event) |
Base class for an electroweak antenna. More... | |
int | getIndexMot () |
Return index. | |
int | getIndexRec () |
bool | isSplitToFermions () |
Check if splitting to fermions, inital, or resoance. | |
virtual bool | isInitial () |
virtual bool | isResonanceDecay () |
bool | selectChannel (int idx, const double &cSum, const map< double, int > &cSumSoFar, int &idi, int &idj, double &mi2, double &mj2) |
Select a channel. | |
Additional Inherited Members | |
Protected Attributes inherited from EWAntenna | |
int | iMot |
Indices, PID, and polarization of I, K in Pythia event record. | |
int | iRec |
int | idMot |
int | idRec |
int | polMot |
Vec4 | pMot |
Momenta of antenna constituents. | |
Vec4 | pRec |
double | sAnt |
Masses and invariants of antenna constituents. | |
double | mMot |
double | mMot2 |
double | mRec |
double | mRec2 |
double | alpha |
Overestimate of QED coupling. | |
int | iSys |
Parton system this antenna is part of. | |
vector< EWBranching > | brVec |
EW branchings. | |
bool | hasTrial |
Trial variables. | |
double | q2Trial |
double | sijTrial |
double | sjkTrial |
int | poliTrial |
int | poljTrial |
vector< Vec4 > | pNew |
Outgoing momenta after branching. | |
double | c0Sum |
Info on coefficents. | |
double | c1Sum |
double | c2Sum |
double | c3Sum |
map< double, int > | c0SumSoFar |
map< double, int > | c1SumSoFar |
map< double, int > | c2SumSoFar |
map< double, int > | c3SumSoFar |
double | q2Match |
Matching scale. | |
int | jNew |
Information for partonSystems. | |
unordered_map< int, int > | iReplace |
double | shat |
EWBranching * | brTrial {} |
Pointers. | |
Info * | infoPtr {} |
Rndm * | rndmPtr {} |
Logger * | loggerPtr {} |
PartonSystems * | partonSystemsPtr {} |
VinciaCommon * | vinComPtr {} |
AlphaEM * | alphaPtr {} |
AmpCalculator * | ampCalcPtr {} |
bool | doBosonInterference |
Settings. | |
int | verbose |
Verbosity. | |
Final-final electroweak antenna.
|
overridevirtual |
Accept trial.
Mark trial as used.
Get some variables.
Check on-shell phase space.
Make sure any new qqbar pair has at least the invariant mass of the lightest meson. sijMin is 0 if these are not quarks.
Construct veto probability, alpha ratio.
Trial kernel - uses kinematic mass of pMot.
Physical kernels per final state spin config - width = 0.
Sums and cumulants for spin states.
Antenna ratio.
Check if veto probability is ok.
Accept/Reject step.
Passed veto - now select a spin state.
Sample Breit-Wigner masses for the kinematics. Don't change resonance mass in an emission.
Recompute invariants to conserve Q2.
Check off-shell phase space.
Check if we veto resonance production with off-shellness larger than q2.
Proceed to do the kinematics. Pre- and post-branching momenta, phi, invariants, masses.
Do kinematics and return.
Apply bosonic interference factor.
Branching accepted!
Implements EWAntenna.
Reimplemented in EWAntennaFFres.
|
overridevirtual |
Generate a trial.
Auxiliary phase space variable ranges. All branching are generated in one go. This implementation overestimates z ranges in massless approx.
Antenna mass is too small for radiation.
Expand square root to prevent numerical issues.
Zeta integrals and trial generation weights.
Find the highest new trial scale and generate scale from c0.
Compute the trial invariants. Use the kinematic mass of pMot here.
Generate scale from c1.
Compute the trial invariants. Use the kinematic mass of pMot here.
Generate scale from c2.
Local veto is required for this trial.
Compute the trial invariants. Use the kinematic mass of pMot here.
Local veto probability to fix the trial distribution. Note that pAccept may be > 1 if sijNow < 0. This is not an issue because these trials get vetoed in acceptTrial.
Save the trial.
Generate scale from c3.
Compute the trial invariants. Use the kinematic mass of pMot here.
Return trial.
Implements EWAntenna.
Reimplemented in EWAntennaFFres.
|
overridevirtual |
Overridden virtual functions.
Final-final electroweak antenna.
Initialize.
Initialize variables.
The shower always uses the on-shell mass of mI. mRec only acts as recoiler, so we just use the kinematic mass.
This function is part of phase space, It is computed with the kinematic masses.
Store the system and branchings.
Find coefficients for overestimates.
Implements EWAntenna.
Reimplemented in EWAntennaFFres.
|
overridevirtual |
Update an event.
Clear information for replacing later in partonSystems.
Fetch ids and kinematic masses.
If new particles are qqbar pair, we need new color indices.
Otherwise carry I color indices to i.
Adjust old particles.
Save information for parton systems.
Implements EWAntenna.
Reimplemented in EWAntennaFFres.