PYTHIA
8.316
|
The ThermalStringFlav class is used to select quark and hadron flavours. More...
#include <ThermalFragmentation.h>
Public Member Functions | |
ThermalStringFlav () | |
Constructor. | |
~ThermalStringFlav () | |
Destructor. | |
void | init () override |
Initialize data members. More... | |
FlavContainer | pick (FlavContainer &flavOld, double pT, double kappaModifier, bool allowPop) override |
Pick a new flavour (including diquarks) given an incoming one. More... | |
virtual int | getHadronIDwin () |
Return chosen hadron in case of thermal model. | |
double | getHadronMassWin (int idHad) override |
Return hadron mass. Used one if present, pick otherwise. | |
int | combineLastThermal (FlavContainer &flav1, FlavContainer &flav2, double pT, double kappaModifier) |
int | getHadronID (FlavContainer &flav1, FlavContainer &flav2, double pT=-1.0, double kappaModifier=-1.0, bool finalTwo=false) override |
Return already set hadron id or combination of the two flavours. | |
void | addQuarkDiquark (vector< pair< int, int > > &quarkCombis, int qID, int diqID, int hadronID) |
Check if quark-diquark combination should be added. If so add. | |
![]() | |
StringFlav () | |
Constructor. | |
virtual | ~StringFlav () |
Destructor. | |
virtual void | init (double kappaModifier, double strangeJunc, double probQQmod) |
Initialise parameters when using close packing. More... | |
int | pickLightQ () |
Pick a light d, u or s quark according to fixed ratios. | |
virtual int | combine (FlavContainer &flav1, FlavContainer &flav2) |
Combine two flavours (including diquarks) to produce a hadron. More... | |
virtual int | combineId (int id1, int id2, bool keepTrying=true) |
Ditto, simplified input argument for simple configurations. | |
virtual pair< int, int > | combineDiquarkJunction (int id1, int id2, int id3) |
Combine three (di-) quark flavours into two hadrons. More... | |
virtual int | combineToLightest (int id1, int id2) |
Combine two flavours to produce a hadron with lowest possible mass. More... | |
virtual int | idLightestNeutralMeson () |
Lightest flavour-neutral meson. | |
void | assignPopQ (FlavContainer &flav) |
Assign popcorn quark inside an original (= rank 0) diquark. More... | |
int | makeDiquark (int id1, int id2, int idHad=0) |
Combine two quarks to produce a diquark. More... | |
void | addQuarkDiquark (vector< pair< int, int > > &quarkCombis, int qID, int diqID, int hadronID) |
Check if quark-diquark combination should be added. If so add. | |
int | getMesonSpinCounter (int hadronID) |
Get spin counter for mesons. | |
double | getFlavourSpinRatios (int i, int j) |
![]() | |
void | initInfoPtr (Info &infoPtrIn) |
This function is called from above for physics objects used in a run. More... | |
virtual | ~PhysicsBase () |
Empty virtual destructor. | |
bool | flag (string key) const |
Shorthand to read settings values. | |
int | mode (string key) const |
double | parm (string key) const |
string | word (string key) const |
vector< bool > | fvec (string key) const |
vector< int > | mvec (string key) const |
vector< double > | pvec (string key) const |
vector< string > | wvec (string key) const |
Protected Attributes | |
bool | mesonNonetL1 |
Settings for thermal model. | |
double | temperature |
double | tempPreFactor |
int | nNewQuark |
double | mesMixRate1 [2][6] |
double | mesMixRate2 [2][6] |
double | mesMixRate3 [2][6] |
double | baryonOctWeight [6][6][6][2] |
double | baryonDecWeight [6][6][6][2] |
map< int, vector< pair< int, int > > > | hadronConstIDs |
Key = hadron id, value = list of constituent ids. | |
map< int, vector< pair< int, int > > > | possibleHadrons |
map< int, vector< double > > | possibleRatePrefacs |
Key = initial (di)quark id, value = prefactor to multiply rate. | |
map< pair< int, int >, vector< pair< int, int > > > | possibleHadronsLast |
Similar, but for combining the last two (di)quarks. Key = (di)quark pair. | |
map< pair< int, int >, vector< double > > | possibleRatePrefacsLast |
int | hadronIDwin |
Selection in thermal model. | |
int | quarkIDwin |
double | hadronMassWin |
![]() | |
bool | suppressLeadingB |
Settings for default Gaussian model. | |
double | probQQtoQ |
double | probStoUD |
double | probSQtoQQ |
double | probQQ1toQQ0 |
double | probQandQQ |
double | probQandS |
double | probQandSinQQ |
double | probQQ1corr |
double | probQQ1corrInv |
double | probQQ1norm |
double | probQQ1join [4] |
double | mesonRate [4][6] |
double | mesonRateSum [4] |
double | mesonMix1 [2][6] |
double | mesonMix2 [2][6] |
double | etaSup |
double | etaPrimeSup |
double | decupletSup |
double | baryonCGSum [6] |
double | baryonCGMax [6] |
double | popcornRate |
double | popcornSpair |
double | popcornSmeson |
double | barCGMax [8] |
double | scbBM [3] |
double | popFrac |
double | popS [3] |
double | dWT [3][7] |
double | lightLeadingBSup |
double | heavyLeadingBSup |
bool | qqKappa |
double | probStoUDSav |
double | probQQtoQSav |
double | probSQtoQQSav |
double | probQQ1toQQ0Sav |
double | alphaQQSav |
bool | closePacking |
Settings for closepacking. | |
bool | doEnhanceDiquark |
double | enhanceStrange |
double | enhancePT |
double | enhanceDiquark |
double | exponentMPI |
double | exponentNSP |
WeightsFragmentation * | wgtsPtr {} |
Fragmentation weights container. | |
![]() | |
Info * | infoPtr = {} |
Settings * | settingsPtr = {} |
Pointer to the settings database. | |
ParticleData * | particleDataPtr = {} |
Pointer to the particle data table. | |
Logger * | loggerPtr = {} |
Pointer to logger. | |
HadronWidths * | hadronWidthsPtr = {} |
Pointer to the hadron widths data table. | |
Rndm * | rndmPtr = {} |
Pointer to the random number generator. | |
CoupSM * | coupSMPtr = {} |
Pointers to SM and SUSY couplings. | |
CoupSUSY * | coupSUSYPtr = {} |
BeamSetup * | beamSetupPtr = {} |
BeamParticle * | beamAPtr = {} |
BeamParticle * | beamBPtr = {} |
BeamParticle * | beamPomAPtr = {} |
BeamParticle * | beamPomBPtr = {} |
BeamParticle * | beamGamAPtr = {} |
BeamParticle * | beamGamBPtr = {} |
BeamParticle * | beamVMDAPtr = {} |
BeamParticle * | beamVMDBPtr = {} |
PartonSystems * | partonSystemsPtr = {} |
Pointer to information on subcollision parton locations. | |
SigmaTotal * | sigmaTotPtr = {} |
Pointers to the total/elastic/diffractive cross sections. | |
SigmaCombined * | sigmaCmbPtr = {} |
set< PhysicsBase * > | subObjects |
UserHooksPtr | userHooksPtr |
mutex * | mutexPtr |
Mutex that should be locked for thread-unsafe code. | |
Additional Inherited Members | |
![]() | |
enum | Status { INCOMPLETE = -1, COMPLETE = 0, CONSTRUCTOR_FAILED, INIT_FAILED, LHEF_END, LOWENERGY_FAILED, PROCESSLEVEL_FAILED, PROCESSLEVEL_USERVETO, MERGING_FAILED, PARTONLEVEL_FAILED, PARTONLEVEL_USERVETO, HADRONLEVEL_FAILED, CHECK_FAILED, OTHER_UNPHYSICAL, HEAVYION_FAILED, HADRONLEVEL_USERVETO } |
Enumerate the different status codes the event generation can have. | |
![]() | |
virtual void | initDerived () |
Initialise derived parameters. More... | |
![]() | |
PhysicsBase () | |
Default constructor. | |
virtual void | onInitInfoPtr () |
virtual void | onBeginEvent () |
This function is called in the very beginning of each Pythia::next call. | |
virtual void | onEndEvent (Status) |
virtual void | onStat () |
This function is called from the Pythia::stat() call. | |
virtual void | onStat (vector< PhysicsBase * >, Pythia *) |
void | registerSubObject (PhysicsBase &pb) |
Register a sub object that should have its information in sync with this. | |
![]() | |
static const int | mesonMultipletCode [6] = { 1, 3, 10003, 10001, 20003, 5} |
Constants: could only be changed in the code itself. More... | |
static const double | baryonCGOct [6] = { 0.75, 0.5, 0., 0.1667, 0.0833, 0.1667} |
static const double | baryonCGDec [6] = { 0., 0., 1., 0.3333, 0.6667, 0.3333} |
The ThermalStringFlav class is used to select quark and hadron flavours.
int combineLastThermal | ( | FlavContainer & | flav1, |
FlavContainer & | flav2, | ||
double | pT, | ||
double | kappaModifier | ||
) |
Combine two flavours into hadron for last two remaining flavours for thermal model.
Combine two flavours (including diquarks) to produce a hadron. Function called in case of combining the two remaining flavours into last hadron.
Determine close-packing scaling.
Decide randomly on whether to treat flav1 or flav2 as incoming.
Temperature increase to work against asymmetry. Apply for s/c/b and diquarks.
Enhanced-rate prefactor for MPIs and/or nearby string pieces.
Get the list of allowed hadrons and constituents for that combination of (di)quarks. First parameter of pair is hadron ID, second is nr of hadron constituents in the list.
Vector with hadron masses. Is -1.0 if m0 is use for calculating the suppression rate and mSel if mSel is used.
Calculate rates/suppression factors for given pT.
Pick mass and calculate suppression factor.
Multiply rate with prefactor.
Save rate and add to sum
Normalize rates
Get accumulated rates
Random number to decide which hadron to pick
Save hadron.
Done.
|
overridevirtual |
Initialize data members.
The ThermalStringFlav class.
Initialize data members of the flavour generation.
Temperature parameters for thermal model.
Hadron multiplets in thermal model.
Parameters for uubar - ddbar - ssbar meson mixing.
Fill in (flavour, spin)-dependent probability of producing the lightest or the lightest two mesons of the nonet.
Fill in rates for multiplication.
Fill list of possible hadrons that are allowed to be produced. Also include a list of "emergency" hadrons that are needed to get rid of all possible endpoint (di)quarks.
Baryon octet and decuplet.
Check how many heavy baryons to include.
Only include lightest combinations.
Only include lightest combinations.
Antibaryons.
Mesons nonets. Take pseudoscalar PDG codes as basis.
Check how many heavy mesons to include. If not included in ordinary production, fill minimal list with "emergency" hadrons
Include all possible combinations, only pseudoscalar as they are the lightest ones.
Include all possible combinations, only pseudoscalar as they are the lightest ones.
Pseudoscalar nonet J=0, S=0, L=0.
Vector nonet J=1, S=1, L=0.
Include L=1 nonets?
Pseudovector nonet J=1, S=0, L=1.
Scalar nonet J=0, S=1, L=1.
Pseudovector nonet J=1, S=1, L=1.
Tensor nonet J=2, S=1, L=1.
Fill list of all hadrons ids (ordinary and "emergency").
Fill map with IDs of hadron constituents for all hadrons.
Baryon can be split into q + qq in several different ways.
Baryon octet J=1/2.
Add (q2+q3)_0/1 + q1. if (q2 < q3) (q2+q3)_0 and if (q2 > q3) (q2+q3)_1.
Add other combinations. Can be both, J=0 or J=1.
(q1+q3)j + q2
(q1+q2)j + q3
Quarks with the same flavour form J=1, all other combinations can be both, J=0 or J=1.
(q1+q2)1 + q3
(q1+q3)1 + q2
(q2+q3)1 + q1
Baryon decuplet J=3/2.
All quark pairs form diquarks with J=1. (q1+q2)1 + q3
(q1+q3)1 + q2
(q2+q3)1 + q1
Mesons usually have a trivial subdivision into quark + antiquark. Mixing of diagonal mesons is taken into account later.
id > 0: downtype+uptype: up = quark, down = antiquark (default) id > 0: same type -> larger id decides
Copy into smaller versions (one for ordinary production, two for "emergency")
List with all possible initial (di)quarks we could get.
If we include heavy quark hadrons we include the following diquarks in addition.
Loop over list with all possible initial (di)quarks. Fill map possibleHadrons with key = initial (di)quark id, value = list of possible hadron ids
For heavy quarks add "emergency" list, if needed.
Fill list: first parameter of pair is hadron ID, second is nr of hadron constituents in the list.
Loop through list with hadrons and their (di)quark content, check if possible to produce given the choice of initial (di)quark.
Loop over constituent IDs.
To include uubar-ddbar-ssbar mixing include all diagonal mesons.
Calculate baryon octet and decuplet weighting factors based on Clebsch-Gordan coefficients and spin counting. Parameters: qDi1 qDi2 q3 spin. Zero for flavour=0 and same flavour diquarks with J=0.
qq0 + r
qq0 + r
Clebsch-Gordon for the rest.
qq1 + q
qq1 + r
qr0 + q
rq0 + q
qr1 + q
rq1 + q
qr0 + s
qr1 + s
Spin 1 diquarks get extra factor of 3. And all factors get relative baryon-to-meson ratio.
Go through the list of possible hadrons and calculate the prefactor that will multiply the rate.
Get hadron and constituents.
Extra suppression factor for s/c/b quarks.
Extra factor according to last digit for spin counting.
Include correct uubar-ddbar-ssbar mixing factor;
Get spin used as counter for the different multiplets
Check if baryon is octet or decuplet.
Make sure ID2 is diquark.
Extract quark flavours and spin from diquark.
Single quark.
Find Clebsch-Gordan: q1 in DQ | q2 in DQ | q3 | S of DQ
Special cases for Lamda (312) and Sigma (321) or the like.
Extract the two lightest quarks from hadron.
Extract the two quarks from the diquark.
Don't do anything if (12) or (21) is diquark.
Sigma (321)
Lamda (312)
Save prefactor.
Now the same again for joining the last two (di)quarks into hadron.
Loop over possible partners, start with next quark.
Skip all combinations with two diquarks.
Skip all combinations with two quarks or two antiquarks.
Skip all combinations with quark-antidiquark and antiquark-diquark. (1 = diquark, 2 = quark not possible).
If we are not including heavy quarks skip combinations of heavy quark - diquark with heavy quark.
Now decide which list of possible hadrons to use. As we might have to use the special list for heavy quarks we use the maximum of the absolute ids in case of two quarks and check the maximum flavour in case of quark - diquark pair.
quark - quark
quark - diquark
Check if diquark contains a heavier flavour than the quark.
New list to fill.
Now loop over possible hadrons and check if other (di)quark in constituents matches idIn2.
Get constituents.
Can take this combination.
Save.
Enhanced-rate prefactor for MPIs and/or nearby string pieces.
Initialize winning parameters.
Reimplemented from StringFlav.
|
overridevirtual |
Pick a new flavour (including diquarks) given an incoming one.
Pick a hadron, based on generated pT value and initial (di)quark. Check all possible hadrons and calculate their relative suppression based on exp(-mThadron/T), possibly multiplied by spin counting, meson mixing or baryon weighting factors. First return value is hadron ID, second new (di)quark ID.
Determine close-packing scaling.
Initial values for new flavour.
Temperature increase to work against asymmetry. Apply for s/c/b and diquarks.
Enhanced-rate prefactor for MPIs and/or nearby string pieces.
Get the list of allowed hadrons and constituents for that initial (di)quark. First parameter of pair is hadron ID, second is nr of hadron constituents in the list.
Vector with hadron masses. Is -1.0 if m0 is use for calculating the suppression rate and mSel if mSel is used.
Calculate rates/suppression factors for given pT.
Pick mass and calculate suppression factor.
Multiply rate with prefactor.
Save rate and add to sum
Normalize rates
Get accumulated rates
Random number to decide which hadron to pick
Get flavour of (di)quark to use next time.
Mesons
Special case for diagonal meson, flavour remains
Baryons
Save new flavour and hadron.
id used to build hadron
id used in next step
Done.
Reimplemented from StringFlav.
|
protected |
Key = initial (di)quark id, value = list of possible hadron ids