./PaxHeaders.94862/pythia8052 0000644 0000000 0000000 00000000072 13466754007 012543 x ustar 00 30 atime=1619556553.820645587 28 ctime=1619556554.2886476 pythia8052/ 0000755 0002170 0000144 00000000000 13466754007 013421 5 ustar 00pythia users 0000000 0000000 pythia8052/PaxHeaders.94862/Makefile 0000644 0000000 0000000 00000000074 13466754007 014125 x ustar 00 30 atime=1619556553.820645587 30 ctime=1619556553.824645604 pythia8052/Makefile 0000644 0002170 0000144 00000007612 13466754007 015067 0 ustar 00pythia users 0000000 0000000 # # Libraries Makefile. Some ideas from Geant4 Makefiles # # M. Kirsanov 07.04.2006 SHELL = /bin/sh # Compilers and compiler options. # Default. FC = g77 CC = gcc CPP = g++ FFLAGS= -O CCFLAGS= -O CPPFLAGS= -g -ansi -pedantic -O -W -Wall # Linux platform with gcc3: as default. ifeq ($(ARCH), Linux) FC = g77 CC = gcc CPP = g++ FFLAGS= -O CCFLAGS= -O CPPFLAGS= -g -ansi -pedantic -O -W -Wall endif # Linux platform with gcc4: new Fortran90 compiler. ifeq ($(ARCH), Linux-gcc4) FC = gfortran CC = gcc CPP = g++ FFLAGS= -O CCFLAGS= -O CPPFLAGS= -g -ansi -pedantic -O -W -Wall endif # Location of directories. TMPDIR=tmp TOPDIR=$(shell \pwd) INCDIR=include SRCDIR=src LIBDIR=lib BINDIR=bin # Location of libraries to be built. targets=$(LIBDIR)/libpythia8.a ifneq (x$(HEPMCLOCATION),x) targets+=$(LIBDIR)/libhepmcinterface.a endif ifeq (x$(PYTHIA6LOCATION),x) targets+=$(LIBDIR)/libpythia6.a endif # Main part: build Pythia8 library. $(TMPDIR)/%.o : $(SRCDIR)/%.cc @mkdir -p $(TMPDIR) $(CPP) $(CPPFLAGS) -c -I$(INCDIR) $< -o $@ # Creating the dependency files *.d # The compiler with option -M is used to build the dependency strings. They # are further edited with sed (stream editor). The first sed command adds the # dependency for the *.d files themselves, the second one is needed because # object files are put in the directory different from src. The last line # removes empty *.d files produced in case of error. $(TMPDIR)/%.d : $(SRCDIR)/%.cc @echo Making dependency for file $<; \ mkdir -p $(TMPDIR); \ $(CC) -M -I$(INCDIR) $< | \ sed 's,\($*\)\.o[ :]*,\1.o $@ : ,g' | \ sed 's/$*.o/$(TMPDIR)\/$*.o/' > $@; \ [ -s $@ ] || rm -f $@ all: $(targets) objects := $(patsubst $(SRCDIR)/%.cc,$(TMPDIR)/%.o,$(wildcard $(SRCDIR)/*.cc)) $(LIBDIR)/libpythia8.a: $(objects) @mkdir -p $(LIBDIR) ar cru $(LIBDIR)/libpythia8.a $(objects) deps := $(patsubst $(SRCDIR)/%.cc,$(TMPDIR)/%.d,$(wildcard $(SRCDIR)/*.cc)) # The "if" below is needed in order to avoid producing the dependency files # when you want to just clean ifneq ($(MAKECMDGOALS),clean) -include $(deps) endif # Build Pythia6 library if a location with existing Pythia6 is not set ifeq (x$(PYTHIA6LOCATION),x) $(TMPDIR)/%.o : pythia6/%.f @mkdir -p $(TMPDIR) $(FC) $(FFLAGS) -c $< -o $@ objectsP6 := $(patsubst pythia6/%.f,$(TMPDIR)/%.o,$(wildcard pythia6/*.f)) $(LIBDIR)/libpythia6.a : $(objectsP6) @mkdir -p $(LIBDIR) ar cru $(LIBDIR)/libpythia6.a $(objectsP6) endif # Build HepMC interface part if HepMC and CLHEP locations are set. ifneq (x$(HEPMCLOCATION),x) ifneq (x$(CLHEPLOCATION),x) $(TMPDIR)/%.o : hepmcinterface/%.cc @mkdir -p $(TMPDIR) $(CPP) $(CPPFLAGS) -c -I$(INCDIR) -I$(HEPMCLOCATION)/include \ -I$(CLHEPLOCATION)/include $< -o $@ $(TMPDIR)/%.d : hepmcinterface/%.cc @echo Making dependency for file $<; \ mkdir -p $(TMPDIR); \ $(CC) -M -I$(INCDIR) -I$(HEPMCLOCATION)/include -I$(CLHEPLOCATION)/include $< | \ sed 's,\($*\)\.o[ :]*,\1.o $@ : ,g' | \ sed 's/$*.o/$(TMPDIR)\/$*.o/' > $@; \ [ -s $@ ] || rm -f $@ objectsI := $(patsubst hepmcinterface/%.cc,$(TMPDIR)/%.o,$(wildcard hepmcinterface/*.cc)) $(LIBDIR)/libhepmcinterface.a : $(objectsI) @mkdir -p $(LIBDIR) ar cru $(LIBDIR)/libhepmcinterface.a $(objectsI) depsI := $(patsubst hepmcinterface/%.cc,$(TMPDIR)/%.d,$(wildcard hepmcinterface/*.cc)) ifneq ($(MAKECMDGOALS),clean) -include $(depsI) endif else $(LIBDIR)/libhepmcinterface.a : hepmcinterface/I_Pythia8.cc @echo ERROR, CLHEPLOCATION should be defined with HEPMCLOCATION endif endif # Clean up: remove (almost?) everything that cannot be recreated. .PHONY: clean clean: rm -f *~; rm -f \#*; rm -rf $(TMPDIR) rm -rf $(LIBDIR) rm -rf $(BINDIR) cd $(SRCDIR); rm -f *~; rm -f \#*; cd - cd $(INCDIR); rm -f *~; rm -f \#*; cd - cd doc; rm -f *~; rm -f \#*; cd - cd pythia6; rm -f *~; rm -f \#*; cd - cd examples; rm -rf *.exe; rm -f *~; rm -f \#*; rm -f core*; cd - pythia8052/PaxHeaders.94862/README 0000644 0000000 0000000 00000000074 13466754007 013345 x ustar 00 30 atime=1619556553.824645604 30 ctime=1619556553.824645604 pythia8052/README 0000644 0002170 0000144 00000003547 13466754007 014312 0 ustar 00pythia users 0000000 0000000 1) First edit file "config.csh" if you use the csh or tcsh shells, else edit "config.sh". (If in doubt, use "echo $SHELL" to find out.) Put there the locations of available libraries. If some libraries are unavailable, comment out the respective lines. If a compiled Pythia6 is unavailable, the source code included in this distribution will be used to build the local library. If HepMC or CLHEP are unavailable, the HepMC interface library will not be built and some examples will be impossible to use. The program will work perfectly well standalone, however. 2) Type "source config.csh" for csh and tcsh shells, else type "source config.sh". This makes the above library paths take effect. 3) Type "gmake" or "make". This will create the libraries (up to three of them). On Linux systems gmake and make are usually equivalent. For others, in case of doubt, use gmake. This applies below as well. 4) The libraries should now be ready for use. To study some examples, go to the "examples" subdirectory. There you can use "gmake mainNN", for NN a two-digit number corresponding to one of the existing examples. This will build an executable "mainNN.exe", that you then can run with "./mainNN.exe > outputfile" (where the initial "./" may not be necessary). 5) Typing "gmake clean" will remove all temporary files, except any output files you may have created during your test runs. (If you do it in the main directory it will apply to all subdirectories, if in examples then only for that subdirectory.) There are no shared libraries for the moment. To learn more about the current program: A) A brief self-contained introduction is found in the pythia8051.pdf file. B) Details can be found by opening a web browser on the doc/Welcome.html file, and then navigating among the pages in the index there. pythia8052/PaxHeaders.94862/config.csh 0000644 0000000 0000000 00000000074 13466754007 014431 x ustar 00 30 atime=1619556553.824645604 30 ctime=1619556553.828645621 pythia8052/config.csh 0000644 0002170 0000144 00000002745 13466754007 015375 0 ustar 00pythia users 0000000 0000000 #!/bin/csh # # Find platform. # setenv ARCH "`uname`" set theGcc1=`g++ --version | awk '{print$3}'` set theGcc=`echo $theGcc1 | awk -F . '{print $1}'` if( $theGcc == 4 ) then if( -e /usr/bin/gfortran ) then setenv ARCH ${ARCH}-gcc4 else echo ERROR: compiler is gcc4 but gfortran is not installed endif endif echo Platform is $ARCH # # Environment variables managing the usage of Pythia6 library. If the # variables below are not set, the local Pythia6 library will be built # and used. # The default values here correspond to CERN AFS (but you may want to change # the version and/or platform). # #setenv PYTHIA6LOCATION /afs/cern.ch/sw/lcg/app/releases/GENSER/GENSER_1_3_0/slc3_ia32_gcc323/lib/archive #setenv PYTHIA6LIBNAME "-lpythia6_326 -ldummy_pythia6_326 -lpdfdummy_pythia6_326" # # Environment variables for building HepMC interface library. Comment out ALL # of them if this interface is not needed or if you don't have both CLHEP # and HepMC in your computer. Note that this library is used by the example # main11. # The default values here correspond to CERN AFS (but you may want to change # the version and/or platform). # #setenv HEPMCLOCATION /afs/cern.ch/sw/lcg/external/HepMC/1.26/slc3_ia32_gcc323 #setenv CLHEPLOCATION /afs/cern.ch/sw/lcg/external/clhep/1.9.2.2/slc3_ia32_gcc323 #if( $?LD_LIBRARY_PATH ) then # setenv LD_LIBRARY_PATH ${LD_LIBRARY_PATH}:${HEPMCLOCATION}/lib:${CLHEPLOCATION}/lib #else # setenv LD_LIBRARY_PATH ${HEPMCLOCATION}/lib:${CLHEPLOCATION}/lib #endif pythia8052/PaxHeaders.94862/config.sh 0000644 0000000 0000000 00000000074 13466754007 014266 x ustar 00 30 atime=1619556553.828645621 30 ctime=1619556553.828645621 pythia8052/config.sh 0000644 0002170 0000144 00000002407 13466754007 015225 0 ustar 00pythia users 0000000 0000000 #!/bin/sh # # Find platform. # export ARCH="`uname`" export theGcc1=`g++ --version | awk '{print$3}'` export theGcc=`echo ${theGcc1} | awk -F . '{print $1}'` if [ ${theGcc} = 4 ]; then export ARCH=${ARCH}-gcc4 fi echo Platform is $ARCH # # Environment variables managing the usage of Pythia6 library. If the # variables below are not set, the local Pythia6 library will be built # and used. # The default values here correspond to CERN AFS (but you may want to change # the version and/or platform). # #export PYTHIA6LOCATION=/afs/cern.ch/sw/lcg/app/releases/GENSER/GENSER_1_3_0/slc3_ia32_gcc323/lib/archive #export PYTHIA6LIBNAME="-lpythia6_326 -ldummypythia6_326 -lpdfdummypythia6_326" # # Environment variables for building HepMC interface library. Comment out ALL # of them if this interface is not needed or if you don't have both CLHEP # and HepMC in your computer. Note that this library is used by the example # main11. # The default values here correspond to CERN AFS (but you may want to change # the version and/or platform). # #export HEPMCLOCATION=/afs/cern.ch/sw/lcg/external/HepMC/1.26/slc3_ia32_gcc323 #export CLHEPLOCATION=/afs/cern.ch/sw/lcg/external/clhep/1.9.2.2/slc3_ia32_gcc323 #export LD_LIBRARY_PATH=${LD_LIBRARY_PATH}:${HEPMCLOCATION}/lib:${CLHEPLOCATION}/lib pythia8052/PaxHeaders.94862/doc 0000644 0000000 0000000 00000000074 13466754007 013155 x ustar 00 30 atime=1619556554.292647617 30 ctime=1619556554.012646414 pythia8052/doc/ 0000755 0002170 0000144 00000000000 13466754007 014166 5 ustar 00pythia users 0000000 0000000 pythia8052/doc/PaxHeaders.94862/Analysis.html 0000644 0000000 0000000 00000000074 13466754007 015703 x ustar 00 30 atime=1619556553.832645639 30 ctime=1619556553.832645639 pythia8052/doc/Analysis.html 0000644 0002170 0000144 00000017736 13466754007 016655 0 ustar 00pythia users 0000000 0000000
Event
format, and use a few basic facilities such
as four-vectors.
isInvisible()
particle method),
and
sph.analyze( event)where
event
is an object of the Event
class,
most likely the pythia.event
one. If the routine returns
false
the analysis failed, e.g. if too few particles are
present to analyze.
After the analysis has been performed, a few Sphericity
class methods are available to return the result of the analysis:
method name="sph()"
Vec4
with vanishing
time/energy component.
method name="list()"
CellJet
(a.k.a. PYCELL
) is a simple cone jet
finder in the UA1 spirit, see the PYTHIA 6 manual. It works in an
(eta, phi, eT) space, where eta is pseudorapidity,
phi azimuthal angle and eT transverse energy.
It will draw cones in R = sqrt(Delta-eta^2 + Delta-phi^2)
around seed cells. If the total eT inside the cone exceeds
the threshold, a jet is formed, and the cells are removed from further
analysis. There are no split or merge procedures, so later-found jet
may be missing some of the edge regions already used up by previous
ones.
A cell jet analysis object is declared by
class name="CellJet cellJet( eTjetMin, coneRadius, select, etaMax,
nEta, nPhi, eTseed, smear, resolution, upperCut, threshold)"
isInvisible()
particle
method),
and
cellJet.analyze( event)where event is an object of the
Event
class, most likely
the pythia.event
one. If the routine returns
false
the analysis failed, but currently this is not
foreseen ever to happen.
After the analysis has been performed, a few CellJet
class methods are available to return the result of the analysis:
method name="size()"
Vec4
corresponding to the four-momentum defined by
the sum of all the contributing cells to the i'th jet, where
each cell contributes a four-momentum as if all the eT is
deposited in the center of the cell,
method name="m(i)"
pMassive
above,
method name="list()"
pMassless
,
for reasons of space).
pythia8052/doc/PaxHeaders.94862/Analysis.xml 0000644 0000000 0000000 00000000074 13466754007 015537 x ustar 00 30 atime=1619556553.832645639
30 ctime=1619556553.832645639
pythia8052/doc/Analysis.xml 0000644 0002170 0000144 00000016615 13466754007 016504 0 ustar 00pythia users 0000000 0000000 Event
format, and use a few basic facilities such
as four-vectors.
isInvisible()
particle method),
andsph.analyze( event)where
event
is an object of the Event
class,
most likely the pythia.event
one. If the routine returns
false
the analysis failed, e.g. if too few particles are
present to analyze.
After the analysis has been performed, a few Sphericity
class methods are available to return the result of the analysis:
Vec4
with vanishing
time/energy component.
CellJet
(a.k.a. PYCELL
) is a simple cone jet
finder in the UA1 spirit, see the PYTHIA 6 manual. It works in an
isInvisible()
particle
method),
andcellJet.analyze( event)where event is an object of the
Event
class, most likely
the pythia.event
one. If the routine returns
false
the analysis failed, but currently this is not
foreseen ever to happen.
After the analysis has been performed, a few CellJet
class methods are available to return the result of the analysis:
Vec4
corresponding to the four-momentum defined by
the sum of all the contributing cells to the pMassive
above,
pMassless
,
for reasons of space).
BeamParticle
class contains information on all partons
extracted from a beam (so far). As each consecutive multiple interaction
defines its respective incoming parton to the hard scattering a
new slot is added to the list. This information is modified when
the backwards evolution of the spacelike shower defines a new
initiator parton. It is used, both for the multiple interactions
and the spacelike showers, to define rescaled parton densities based
on the x and flavours already extracted, and to distinguish
between valence, sea and companion quarks. Once the perturbative
evolution is finished, further beam remnants are added to obtain a
consistent set of flavours. The current physics framework is further
described in [Sjo04].
Much of the above information is stored in a vector of
ResolvedParton
objects, which each contains flavour and
momentum information, as well as valence/companion information and more.
The BeamParticle
method list()
shows the contents of
this vector, mainly for debug purposes.
The BeamRemnants
class takes over for the final step of adding
primordial kT to the initiators and remnants, assigning the
relative longitudinal momentum sharing among the remnants, and
constructing the overall kinematics. This step couples the two sides
of an event, and could therefore not be covered in the
BeamParticle
class, which only considers one beam at a time.
Neither of the methods of these classes are intended for general use,
and so are not described here.
Pythia
class, see the
Generic page. Then pointers to the pdf's are handed on to
BeamParticle
at initialization, for all subsequent usage.
valenceDiqEnhance
factor.
parameter name="Beams:valencePowerMeson" default="0.8" min="0."
off
option is intended for debug purposes only, as follows.
When more than one valence quark is kicked out of a baryon beam,
as part of the multiple interactions scenario, the subsequent
hadronization is described in terms of a junction string topology.
This description involves a number of technical complications that
may make the program more unstable. As an alternative, by switching
this option off, junction configurations are rejected, and the
multiple interactions and their showers are redone until a
junction-free topology is found.
diffLargeMassSuppress
parameter.
pythia8052/doc/PaxHeaders.94862/BeamRemnants.xml 0000644 0000000 0000000 00000000074 13466754007 016330 x ustar 00 30 atime=1619556553.836645656
30 ctime=1619556553.836645656
pythia8052/doc/BeamRemnants.xml 0000644 0002170 0000144 00000016427 13466754007 017276 0 ustar 00pythia users 0000000 0000000 BeamParticle
class contains information on all partons
extracted from a beam (so far). As each consecutive multiple interaction
defines its respective incoming parton to the hard scattering a
new slot is added to the list. This information is modified when
the backwards evolution of the spacelike shower defines a new
initiator parton. It is used, both for the multiple interactions
and the spacelike showers, to define rescaled parton densities based
on the ResolvedParton
objects, which each contains flavour and
momentum information, as well as valence/companion information and more.
The BeamParticle
method list()
shows the contents of
this vector, mainly for debug purposes.
The BeamRemnants
class takes over for the final step of adding
primordial BeamParticle
class, which only considers one beam at a time.
Neither of the methods of these classes are intended for general use,
and so are not described here.
Pythia
class, see the
Generic page. Then pointers to the pdf's are handed on to
BeamParticle
at initialization, for all subsequent usage.
valenceDiqEnhance
factor.
off
option is intended for debug purposes only, as follows.
When more than one valence quark is kicked out of a baryon beam,
as part of the multiple interactions scenario, the subsequent
hadronization is described in terms of a junction string topology.
This description involves a number of technical complications that
may make the program more unstable. As an alternative, by switching
this option off, junction configurations are rejected, and the
multiple interactions and their showers are redone until a
junction-free topology is found.
diffLargeMassSuppress
parameter.
PythiaStdlib
only exists as a header file, collecting all
the include
and using
statements that
are required by most other classes to access the C++
Stdlib
containers and methods.
In addition, the file contains inline functions pow2(x)
,
pow3(x)
, pow4(x)
and pow5(x)
,
for small integer powers, sqrtpos(x)
where a
max(0., x)
ensures one does not take the square root of
a negative number, and tolower(string)
that converts
a whole string to lowercase characters (extending on the
tolower
fuction for a single character).
pythia8052/doc/PaxHeaders.94862/Comments.xml 0000644 0000000 0000000 00000000074 13466754007 015541 x ustar 00 30 atime=1619556553.840645673
30 ctime=1619556553.840645673
pythia8052/doc/Comments.xml 0000644 0002170 0000144 00000001664 13466754007 016504 0 ustar 00pythia users 0000000 0000000 PythiaStdlib
only exists as a header file, collecting all
the include
and using
statements that
are required by most other classes to access the C++
Stdlib
containers and methods.
In addition, the file contains inline functions pow2(x)
,
pow3(x)
, pow4(x)
and pow5(x)
,
for small integer powers, sqrtpos(x)
where a
max(0., x)
ensures one does not take the square root of
a negative number, and tolower(string)
that converts
a whole string to lowercase characters (extending on the
tolower
fuction for a single character).
Event
class for event records basically is a vector of
Particle
s, so that it can expand to fit the event size.
The index operator is overloaded, so that event[i]
corresponds to the i
'th particle of an Event
object event
. Thus event[i].id()
returns the
identity of the i
'th particle. References to the first,
i
'th and last particle are obtained with
event.front()
, event.at(i)
and
event.back()
, respectively.
The event size can be found with size()
, i.e. valid
particles are stored in 0 <= i < event.size()
.
Line 0 is used to represent the event as a whole, with its total
four-momentum and invariant mass, but does not form part of the
event history. Lines 1 and 2 contains the two incoming beams, and
only from here on history tracing works as could be expected. That
way unassigned mother and daughter indices can be put 0 without
ambiguity. Depending on the task at hand, a loop may therefore start
at an index 1 without any loss. Specifically, for translation to other
event record formats such as HepMC [Dob01], where the first
index is 1, the Pythia entry 0 definitely ought to be skipped in order
to minimize the danger of errors.
New particles can be added to the end of the current event record
with append(Particle)
, or
append(id, status, mother1, mother2, daughter1, daughter2,
col, acol, p, m, scale)
where p
is the four-momentum vector, and everything except
id
defaults to 0. The append
method returns
the index of the new particle position.
The existing particle at index iCopy
can be copied to the end
with copy(iCopy, newStatus = 0)
. By default, i.e. with
newStatus = 0
, everything is copied precisely as it is,
which means that history information has to be modified further by hand
to make sense. With a positive newStatus
, the new copy is set
up to be the daughter of the old, with status code newStatus
,
and the status code of iCopy
is negated. With a negative
newStatus
, the new copy is instead set up to be the mother
of iCopy
.
A listing of the whole event is obtained with list()
. This
function takes an ostream
object as an optional argument.
The basic id, status, mother, daughter, colour, four-momentum
and mass data are always given, but the following switches can be
set to provide further information, or restrict the set of particles
listed:
flag name="Event:listFinalOnly" default="off"
motherList(i)
and daughterList(i)
methods
described below. It is purely intended for debug purposes,
e.g. when writing an interface to another event record format.
flag name="Event:extraBlankLine" default="off"
Event:listMothersAndDaughters
on, and a
longer listing.
flag name="Event:listJunctions" default="off"
Event::initStatic()
yourself for them to take effect.
The event record can be emptied for the next event by a
clear()
. The last n
entries can be removed by
popBack(n)
, where n = 1
by default.
The user would normally be concerned with the Event
object
that is a public member event
of the Pythia
class.
Thus, having declared a Pythia
object pythia
,
pythia.event[i].id()
would be used to return the identity
of the i
'th particle, and pythia.event.size()
to
give the size of the event record.
A Pythia
object contains a second event record for the
hard process alone, called process
, used as input for the
generation of the complete event. Thus one may e.g. call either
pythia.process.list()
or pythia.event.list()
.
To distinguish those two rapidly at visual inspection, the
"Pythia Event Listing"
header is printed out differently,
in one case adding "(hard process)"
and in the other
"(complete event)"
. This is set by a header(string) method.
One data member in an Event
object is used to keep track of the
largest col()
or acol()
tag set so far, so that
new ones do not clash. The lastcolTag()
method returns the
last tag assigned, i.e. largest value in the current event, and
nextColTag()
ups it by one before returing the value. The
latter method thus is used when a new colour tag is needed.
mode name="Event:startColTag" default="100" min="0" max="1000"
startColTag+1
, etc. The Les Houches accord [Boo01]
suggests this number to be 500, but 100 works equally well.
The scale()
methods can be used to set or get the scale
(in GeV) of the event as a whole. Further methods for event properties
may be added in the future.
There are also a few methods with an individual particle index
i
as input, but requiring some search operations in
the event record, and therefore not possible to define inside the
Particle
class:
method name="motherList(i)"
vector<int>
containing a list of all the
mothers of a particle. This list is empty for entry 0 , while
normally it contains one, two or many mothers. The latter case
applies e.g. to string fragmentation, where the whole fragmenting
system is counted as mothers. Mothers are listed in ascending order.
method name="daughterList(i)"
vector<int>
containing a list of all the
daughters of a particle. This list is empty for a particle that did
not decay (or, if the evolution is stopped early enough, a parton
that did not branch), while otherwise it can contain a list of
varying length, from one to many. Many partons may have the same
daughterList
, e.g. in the hard process and fragmentation steps.
For the two incoming beam particles, all shower initiators and beam
remnants are counted as daughters, with the one in slot 0 being
the one leading up to the hardest interaction.
method name="iTopCopy(i), iBotCopy(i)"
i
up
to its top mother or down to its bottom daughter. If there are no such
carbon copies, i
itself will be returned.
method name="iTopCopyId(i), iBotCopyId(i)"
id
code. The behaviour when trying to
trace a gluon through a shower, with its g -> g g branchings,
is then rather unpredictable. Similarly, a hard scattering such as
u u -> u u is better avoided. It should work well for "rare"
particles, not found in the beams, so that the program is not fooled
by ambiguities.
method name="sisterList(i)"
vector<int>
containing a list of all the
sisters of a particle, i.e. all the daughters of the first mother,
except the particle itself.
method name="sisterListTopBot(i)"
vector<int>
containing a list of all the
sisters of a particle, tracking up and back down through carbon copies
if required. That is, the particle is first traced up with
iTopCopy()
before its mother is found, and then all
the particles in the daughterList()
of this mother are
traced down with iBotCopy()
, omitting the original
particle itself.
method name="isAncestor(i, iAncestor)"
i
upwards through mother, grandmother,
and so on, until either iAncestor
is found or the top of
the record is reached. Normally one unique mother is required,
as is the case e.g. in decay chains or in parton showers, so that
e.g. the tracing through a hard scattering would not work. For
hadronization, currently first-rank hadrons are identified with the
respective string endpoint quark, which may be useful e.g. for b
physics. However, currently ministrings or junction topologies give
false
. (??)
kind = 1
: incoming colourless particle to three
outgoing colours (e.g. baryon beam remnant or
neutralino -> q q q
);kind = 2
: incoming colourless particle to three
outgoing anticolours;kind = 3
: one incoming anticolor (stored first)
and two outgoing colours (e.g. antisquark decaying to quark);kind = 4
: one incoming color (stored first) and two
outgoing anticolours;kind = 5
: incoming colour octet to three colours,
where the incoming colour passes through unchanged and so need not
be bokkept here, while the incoming anticolor (stored first) and the
two outgoing colours are (e.g. gluino decay to three quarks);kind = 6
: incoming colour octet to three anticolours,
where the incoming anticolour passes through unchanged and so need not
be bookkept here, while the incoming color (stored first) and the two
outgoing colours are.kind
codes corresponds to a +1 (-1) change in
baryon number across the junction.
kind = 1, 2
are
implemented.
The kind and colour information in the list of junctions can be set
or read with methods of the Event
class, but are not of
common interest and so not described here.
pythia8052/doc/PaxHeaders.94862/Event.xml 0000644 0000000 0000000 00000000072 13466754007 015033 x ustar 00 29 atime=1619556553.84464569
29 ctime=1619556553.84464569
pythia8052/doc/Event.xml 0000644 0002170 0000144 00000027720 13466754007 016001 0 ustar 00pythia users 0000000 0000000 Event
class for event records basically is a vector of
Particle
s, so that it can expand to fit the event size.
The index operator is overloaded, so that event[i]
corresponds to the i
'th particle of an Event
object event
. Thus event[i].id()
returns the
identity of the i
'th particle. References to the first,
i
'th and last particle are obtained with
event.front()
, event.at(i)
and
event.back()
, respectively.
The event size can be found with size()
, i.e. valid
particles are stored in 0 <= i < event.size()
.
Line 0 is used to represent the event as a whole, with its total
four-momentum and invariant mass, but does not form part of the
event history. Lines 1 and 2 contains the two incoming beams, and
only from here on history tracing works as could be expected. That
way unassigned mother and daughter indices can be put 0 without
ambiguity. Depending on the task at hand, a loop may therefore start
at an index 1 without any loss. Specifically, for translation to other
event record formats such as HepMC Dob01, where the first
index is 1, the Pythia entry 0 definitely ought to be skipped in order
to minimize the danger of errors.
New particles can be added to the end of the current event record
with append(Particle)
, or
append(id, status, mother1, mother2, daughter1, daughter2,
col, acol, p, m, scale)
where p
is the four-momentum vector, and everything except
id
defaults to 0. The append
method returns
the index of the new particle position.
The existing particle at index iCopy
can be copied to the end
with copy(iCopy, newStatus = 0)
. By default, i.e. with
newStatus = 0
, everything is copied precisely as it is,
which means that history information has to be modified further by hand
to make sense. With a positive newStatus
, the new copy is set
up to be the daughter of the old, with status code newStatus
,
and the status code of iCopy
is negated. With a negative
newStatus
, the new copy is instead set up to be the mother
of iCopy
.
A listing of the whole event is obtained with list()
. This
function takes an ostream
object as an optional argument.
The basic id, status, mother, daughter, colour, four-momentum
and mass data are always given, but the following switches can be
set to provide further information, or restrict the set of particles
listed:
motherList(i)
and daughterList(i)
methods
described below. It is purely intended for debug purposes,
e.g. when writing an interface to another event record format.
Event:listMothersAndDaughters
on, and a
longer listing.
Event::initStatic()
yourself for them to take effect.
The event record can be emptied for the next event by a
clear()
. The last n
entries can be removed by
popBack(n)
, where n = 1
by default.
The user would normally be concerned with the Event
object
that is a public member event
of the Pythia
class.
Thus, having declared a Pythia
object pythia
,
pythia.event[i].id()
would be used to return the identity
of the i
'th particle, and pythia.event.size()
to
give the size of the event record.
A Pythia
object contains a second event record for the
hard process alone, called process
, used as input for the
generation of the complete event. Thus one may e.g. call either
pythia.process.list()
or pythia.event.list()
.
To distinguish those two rapidly at visual inspection, the
"Pythia Event Listing"
header is printed out differently,
in one case adding "(hard process)"
and in the other
"(complete event)"
. This is set by a header(string) method.
One data member in an Event
object is used to keep track of the
largest col()
or acol()
tag set so far, so that
new ones do not clash. The lastcolTag()
method returns the
last tag assigned, i.e. largest value in the current event, and
nextColTag()
ups it by one before returing the value. The
latter method thus is used when a new colour tag is needed.
startColTag+1
, etc. The Les Houches accord Boo01
suggests this number to be 500, but 100 works equally well.
scale()
methods can be used to set or get the scale
(in GeV) of the event as a whole. Further methods for event properties
may be added in the future.
There are also a few methods with an individual particle index
i
as input, but requiring some search operations in
the event record, and therefore not possible to define inside the
Particle
class:
vector<int>
containing a list of all the
mothers of a particle. This list is empty for entry 0 , while
normally it contains one, two or many mothers. The latter case
applies e.g. to string fragmentation, where the whole fragmenting
system is counted as mothers. Mothers are listed in ascending order.
vector<int>
containing a list of all the
daughters of a particle. This list is empty for a particle that did
not decay (or, if the evolution is stopped early enough, a parton
that did not branch), while otherwise it can contain a list of
varying length, from one to many. Many partons may have the same
daughterList
, e.g. in the hard process and fragmentation steps.
For the two incoming beam particles, all shower initiators and beam
remnants are counted as daughters, with the one in slot 0 being
the one leading up to the hardest interaction.
i
up
to its top mother or down to its bottom daughter. If there are no such
carbon copies, i
itself will be returned.
id
code. The behaviour when trying to
trace a gluon through a shower, with its vector<int>
containing a list of all the
sisters of a particle, i.e. all the daughters of the first mother,
except the particle itself.
vector<int>
containing a list of all the
sisters of a particle, tracking up and back down through carbon copies
if required. That is, the particle is first traced up with
iTopCopy()
before its mother is found, and then all
the particles in the daughterList()
of this mother are
traced down with iBotCopy()
, omitting the original
particle itself.
i
upwards through mother, grandmother,
and so on, until either iAncestor
is found or the top of
the record is reached. Normally one unique mother is required,
as is the case e.g. in decay chains or in parton showers, so that
e.g. the tracing through a hard scattering would not work. For
hadronization, currently first-rank hadrons are identified with the
respective string endpoint quark, which may be useful e.g. for false
. (??)
kind = 1
: incoming colourless particle to three
outgoing colours (e.g. baryon beam remnant or
neutralino -> q q q
);kind = 2
: incoming colourless particle to three
outgoing anticolours;kind = 3
: one incoming anticolor (stored first)
and two outgoing colours (e.g. antisquark decaying to quark);kind = 4
: one incoming color (stored first) and two
outgoing anticolours;kind = 5
: incoming colour octet to three colours,
where the incoming colour passes through unchanged and so need not
be bokkept here, while the incoming anticolor (stored first) and the
two outgoing colours are (e.g. gluino decay to three quarks);kind = 6
: incoming colour octet to three anticolours,
where the incoming anticolour passes through unchanged and so need not
be bookkept here, while the incoming color (stored first) and the two
outgoing colours are.kind
codes corresponds to a +1 (-1) change in
baryon number across the junction.
kind = 1, 2
are
implemented.
The kind and colour information in the list of junctions can be set
or read with methods of the Event
class, but are not of
common interest and so not described here.
Info
class collects various one-of-a-kind information,
some relevant for all events and others for the current event.
An object info
is a public member of the Pythia
class, so if you e.g. have declared Pythia pythia
, the
Info
methods can be accessed by
pythia.info.method()
. Most of this is information that
could also be obtained e.g. from the event record, but is here more
directly available. It is intended for processes generated internally
in Pythia8, not for ones read in e.g. via the Les Houches Accord.
Here are the currently available methods:
method name="list()"
hasSub()
is true. For instance, for
a minimum-bias event the code
would always be 101,
while codeSub()
would vary depending on the actual
hardest interaction, e.g. 111 for g g -> g g. The methods
below would also provide information for this particular subcollision.
method name="id1(), id2()"
Info
class collects various one-of-a-kind information,
some relevant for all events and others for the current event.
An object info
is a public member of the Pythia
class, so if you e.g. have declared Pythia pythia
, the
Info
methods can be accessed by
pythia.info.method()
. Most of this is information that
could also be obtained e.g. from the event record, but is here more
directly available. It is intended for processes generated internally
in Pythia8, not for ones read in e.g. via the Les Houches Accord.
Here are the currently available methods:
hasSub()
is true. For instance, for
a minimum-bias event the code
would always be 101,
while codeSub()
would vary depending on the actual
hardest interaction, e.g. 111 for pythia.statistics();assuming
pythia
is an instance of the Pythia
class. The statistics
method in its turn calls on the
methods below.
ProcessLevel::statistics()
member will loop over the
list of existing processes, and for each write out name, code,
the number of tried and accepted events (in the sense that the
Monte Carlo integration involves an initial upper estimate of the
cross section, that then is compensated by the try-accept/reject
step), the cross section and the estimated error on the latter.
If PYTHIA 6 is used to generate the hard processes, instead the
Fortran PYSTAT
method is called.
ErrorMessages
class is rather small. It is handed any
warning or error messages during the event generation phase, and will
store each distinct message, with a counter for how many times it is
issued. Thus it is possible to limit the number of identical messages
issued. The summary table printed by pythia.statistics()
provides a table with all the different messages issued, in
alphabetical order, with the total number of times each was generated.
There is only one mode affecting its operation:
mode name="ErrorMessages:timesToPrint" default="1" min="0"
pythia.statistics();assuming
pythia
is an instance of the Pythia
class. The statistics
method in its turn calls on the
methods below.
ProcessLevel::statistics()
member will loop over the
list of existing processes, and for each write out name, code,
the number of tried and accepted events (in the sense that the
Monte Carlo integration involves an initial upper estimate of the
cross section, that then is compensated by the try-accept/reject
step), the cross section and the estimated error on the latter.
If PYTHIA 6 is used to generate the hard processes, instead the
Fortran PYSTAT
method is called.
ErrorMessages
class is rather small. It is handed any
warning or error messages during the event generation phase, and will
store each distinct message, with a counter for how many times it is
issued. Thus it is possible to limit the number of identical messages
issued. The summary table printed by pythia.statistics()
provides a table with all the different messages issued, in
alphabetical order, with the total number of times each was generated.
There is only one mode affecting its operation:
DecayHandler
is a base class for the external handling of
decays. The user-written derived class is called if a pointer to it has
been given with the pythia.decayPtr()
method, where it also
is specified which particles it will be called for.
There is only one pure virtual method in DecayHandler
,
to do the decay:
method name="virtual bool decay(vector<int>& idProd,
vector<double>& mProd, vector<Vec4>& pProd, int iDec,
const Event& event)"
idProd[0]
,
mProd[0]
and pProd[0]
contains information on the
particle that is to be decayed. At output, the vectors should have
increased by the addition of all the decay products. Even if initially
defined in the rest frame of the mother, the products should have been
boosted so that their four-momenta add up to the pProd[0]
of
the decaying particle.
Should it be of interest to know the prehistory of the decaying
particle, e.g. to set some helicity information affecting the
decay angular distribution, the full event record is available
read-only, with info in which slot iDec
the decaying particle
is stored.
The routine should return true
if it managed the decay and
false
otherwise, in which case Pythia
will try
to do the decay itself.
Note that the decay vertex is always set by Pythia
, and that
B-Bbar oscillations have already been taken into account,
if they were switched on. Thus idProd[0]
may be the opposite
of event[iDec].id()
, where the latter provides the code at
production.
A sample test program is available in main04.cc
, providing
a simple example of how to use thing facility.
pythia8052/doc/PaxHeaders.94862/ExternalDecays.xml 0000644 0000000 0000000 00000000074 13466754007 016667 x ustar 00 30 atime=1619556553.856645742
30 ctime=1619556553.856645742
pythia8052/doc/ExternalDecays.xml 0000644 0002170 0000144 00000004352 13466754007 017627 0 ustar 00pythia users 0000000 0000000 DecayHandler
is a base class for the external handling of
decays. The user-written derived class is called if a pointer to it has
been given with the pythia.decayPtr()
method, where it also
is specified which particles it will be called for.
There is only one pure virtual method in DecayHandler
,
to do the decay:
idProd[0]
,
mProd[0]
and pProd[0]
contains information on the
particle that is to be decayed. At output, the vectors should have
increased by the addition of all the decay products. Even if initially
defined in the rest frame of the mother, the products should have been
boosted so that their four-momenta add up to the pProd[0]
of
the decaying particle.
Should it be of interest to know the prehistory of the decaying
particle, e.g. to set some helicity information affecting the
decay angular distribution, the full event record is available
read-only, with info in which slot iDec
the decaying particle
is stored.
The routine should return true
if it managed the decay and
false
otherwise, in which case Pythia
will try
to do the decay itself.
Note that the decay vertex is always set by Pythia
, and that
idProd[0]
may be the opposite
of event[iDec].id()
, where the latter provides the code at
production.
A sample test program is available in main04.cc
, providing
a simple example of how to use thing facility.
.h
in the include
subdirectory, where the public interface of the class is declared,
and inline methods are defined..cc
in the src
subdirectory, where most of the methods are implemented..d
, and
compiled code, .o
are created in the tmp
subdirectory.
In part the .xml
documentation files in the doc
subdirectory have matching names, but the match is broken by the
desire to group topics more by user interaction than internal
operation. On these pages the function of the different code files
is summarized. Currently, each .xml
file is also roughly
translated into an .html
one, to allow easy viewing of
the contents in a web browser, but without the interactivity
eventually strived for.
file name="Pythia"
Pythia
,
so only a Pythia
object needs to be explicitly instantiated
and addressed by the user.
file name="Event"
Particle
class, used by
Event
. Pythia
uses two Event
objects, one for the process-level record (process
) and
one for the complete event (event
).
file name="Information"
Pythia6.cc
file but instead a Pythia6.f
one, located in the pythia6
subdirectory. The header file
contains the required wrappers to access the relevant parts of the
Fortran code from C++.
file name="LesHouches"
Pythia
. Should be linked to external programs or files.
file name="TimeShower"
xml
files, but these values can then
be changed by the user.
file name="ParticleData"
Rndm
, a four-vector class Vec4
, and a
histogram class Hist
.
file name="PythiaStdlib"
.h
file, containing all the Stdlib
headers used in Pythia 8, with using
directives. Also
inlines pow2(x)
- pow5(x)
for fast powers of
small integers, sqrtpos(x)
to avoid that roundoff gives the
square root of a negative number, and tolower(string)
for
producing a lowercase string.
pythia8052/doc/PaxHeaders.94862/Files.xml 0000644 0000000 0000000 00000000074 13466754007 015016 x ustar 00 30 atime=1619556553.856645742
30 ctime=1619556553.860645759
pythia8052/doc/Files.xml 0000644 0002170 0000144 00000016514 13466754007 015761 0 ustar 00pythia users 0000000 0000000 .h
in the include
subdirectory, where the public interface of the class is declared,
and inline methods are defined..cc
in the src
subdirectory, where most of the methods are implemented..d
, and
compiled code, .o
are created in the tmp
subdirectory.
In part the .xml
documentation files in the doc
subdirectory have matching names, but the match is broken by the
desire to group topics more by user interaction than internal
operation. On these pages the function of the different code files
is summarized. Currently, each .xml
file is also roughly
translated into an .html
one, to allow easy viewing of
the contents in a web browser, but without the interactivity
eventually strived for.
Pythia
,
so only a Pythia
object needs to be explicitly instantiated
and addressed by the user.
Particle
class, used by
Event
. Pythia
uses two Event
objects, one for the process-level record (process
) and
one for the complete event (event
).
Pythia6.cc
file but instead a Pythia6.f
one, located in the pythia6
subdirectory. The header file
contains the required wrappers to access the relevant parts of the
Fortran code from C++.
Pythia
. Should be linked to external programs or files.
xml
files, but these values can then
be changed by the user.
Rndm
, a four-vector class Vec4
, and a
histogram class Hist
.
.h
file, containing all the Stdlib
headers used in Pythia 8, with using
directives. Also
inlines pow2(x)
- pow5(x)
for fast powers of
small integers, sqrtpos(x)
to avoid that roundoff gives the
square root of a negative number, and tolower(string)
for
producing a lowercase string.
StringFlav
class handles the choice of a new flavour
in the fragmentation process, and the production of a new hadron
from a set of input flavours. It is mainly used by the string
fragmentation machinery (including ministrings), but also e.g.
in some particle decays and for some beam-remnant cases. The basic
concepts are in agreement with [And83].
Currently only simple diquark production of baryons is
included; the more sophisticated popcorn scenario remains to be
implemented. Therefore the number of parameters on this page will
be expanded later on.
StringFlav
class handles the choice of a new flavour
in the fragmentation process, and the production of a new hadron
from a set of input flavours. It is mainly used by the string
fragmentation machinery (including ministrings), but also e.g.
in some particle decays and for some beam-remnant cases. The basic
concepts are in agreement with And83.
Currently only simple diquark production of baryons is
included; the more sophisticated popcorn scenario remains to be
implemented. Therefore the number of parameters on this page will
be expanded later on.
Pythia
class should
be used in the user-supplied main program, further outlined in the
following.
#include "Pythia.h"To simplify typing, it also makes sense to declare
using namespace Pythia8;
Pythia pythia;It is this object that we will use from now on. Normally a run will only contain one
Pythia
object, but hypothetically
you could use several.Pythia
will be on the cout
stream. If this is not convenient, you can give a reference to another
stream as an optional argument
Pythia pythia(ostream);but this does not work so far! (??)
Pythia
constructor. A third static class, eventually to disappear, is the
interface to the old Fortran 77 PYTHIA 6, currently used for
generating hard processes.pythia.readString(string)method provides a covenient uniform interface to all three of them. The information in the string is case-insensitive, but upper- and lowercase can be combined for clarity. The rules are that
Pythia6:
, this part is peeled
off, and the rest is sent on to Fortran PYTHIA 6, using the
pygive
method in that package;pythia.particleData.readString(string)
;pythia.settings.readString(string)
.false
is used to
switch off warnings.pythia.readString("Pythia6:msel = 6"); pythia.readString("111:mayDecay = false"); pythia.readString("TimeShower:pTmin = 1.0");The methods in this paragraph are intended for small changes; for more extensive ones it is better to store all the changes in a file, see next.
pythia.readFile(fileName);Each line in this file with be processes by the
pythia.readString()
method introduced above. You can thus
freely mix comment lines and lines handed on to Settings
,
ParticleDataTable
and Pythia6
.
This would be the normal way to set up what a run is supposed
to do. Again, an optional second argument false
allows you to
switch off warning messages for unknown variables.readString
method above may be more appropriate.Pythia
, you can suppy your own by a call to the
PDFptr
method
pythia.PDFptr( pdfA, pdfB);where
pdfA
and pdfB
are pointers to two
Pythia
PDF objects
(further
instructions).
Note that pdfA
and pdfB
cannot point to
the same object; even if the PDF set is the same, two copies are
needed to keep track of two separate sets of x and density
values.decayPtr
method
pythia.decayPtr( decayHandler, particles)where the
decayHandler
derives from the
DecayHandler
base class and particles
is a
vector of particle codes to be handled
(further instructions).
rndmEnginePtr
method
pythia.rndmEnginePtr( rndmEngine)where
rndmEngine
derives from the RndmEngine
base class (further
instructions). The Pythia
default random number
generator is perfectly good, so this is only intended for consistency
in bigger frameworks.init
method allows a few different input formats,
so you can pick the one convenient for you:LHAinit
class object, and that LHA event information
will be provided by the LHAevnt
class object
(further instructions).pythia.settings.listChanged(); pythia.settings.listAll(); pythia.particleData.listChanged(); pythia.particleData.listAll();
next
method,
pythia.next();This method takes no arguments; everything has already been specified. It does return a bool value, however,
false
when the
generation failed. This can be a "programmed death" when the
supply of input parton-level configurations on file is exhausted,
but also caused by a failure of Pythia
to generate an event,
or that an event was generated but something strange was detected
in it.event
object, of type Event
, which is a public member of
pythia
. You therefore have access to all the tools described
on the pages under the "Study Output" header in the index. For instance,
an event can be listed with
pythia.event.list()
, the identity of the i'th
particle is given by pythia.event[i].id()
, and so on.process
, also of type
Event
.pythia.statistics();to get some run statistics, on cross sections and the number of errors and warnings encountered.
Pythia
class should
be used in the user-supplied main program, further outlined in the
following.
#include "Pythia.h"To simplify typing, it also makes sense to declare
using namespace Pythia8;
Pythia pythia;It is this object that we will use from now on. Normally a run will only contain one
Pythia
object, but hypothetically
you could use several.Pythia
will be on the cout
stream. If this is not convenient, you can give a reference to another
stream as an optional argument
Pythia pythia(ostream);but this does not work so far! (??)
Pythia
constructor. A third static class, eventually to disappear, is the
interface to the old Fortran 77 PYTHIA 6, currently used for
generating hard processes.pythia.readString(string)method provides a covenient uniform interface to all three of them. The information in the string is case-insensitive, but upper- and lowercase can be combined for clarity. The rules are that
Pythia6:
, this part is peeled
off, and the rest is sent on to Fortran PYTHIA 6, using the
pygive
method in that package;pythia.particleData.readString(string)
;pythia.settings.readString(string)
.false
is used to
switch off warnings.pythia.readString("Pythia6:msel = 6"); pythia.readString("111:mayDecay = false"); pythia.readString("TimeShower:pTmin = 1.0");The methods in this paragraph are intended for small changes; for more extensive ones it is better to store all the changes in a file, see next.
pythia.readFile(fileName);Each line in this file with be processes by the
pythia.readString()
method introduced above. You can thus
freely mix comment lines and lines handed on to Settings
,
ParticleDataTable
and Pythia6
.
This would be the normal way to set up what a run is supposed
to do. Again, an optional second argument false
allows you to
switch off warning messages for unknown variables.readString
method above may be more appropriate.Pythia
, you can suppy your own by a call to the
PDFptr
method
pythia.PDFptr( pdfA, pdfB);where
pdfA
and pdfB
are pointers to two
Pythia
PDF objects
(further
instructions).
Note that pdfA
and pdfB
cannot point to
the same object; even if the PDF set is the same, two copies are
needed to keep track of two separate sets of decayPtr
method
pythia.decayPtr( decayHandler, particles)where the
decayHandler
derives from the
DecayHandler
base class and particles
is a
vector of particle codes to be handled
(further instructions).
rndmEnginePtr
method
pythia.rndmEnginePtr( rndmEngine)where
rndmEngine
derives from the RndmEngine
base class (further
instructions). The Pythia
default random number
generator is perfectly good, so this is only intended for consistency
in bigger frameworks.init
method allows a few different input formats,
so you can pick the one convenient for you:LHAinit
class object, and that LHA event information
will be provided by the LHAevnt
class object
(further instructions).pythia.settings.listChanged(); pythia.settings.listAll(); pythia.particleData.listChanged(); pythia.particleData.listAll();
next
method,
pythia.next();This method takes no arguments; everything has already been specified. It does return a bool value, however,
false
when the
generation failed. This can be a "programmed death" when the
supply of input parton-level configurations on file is exhausted,
but also caused by a failure of Pythia
to generate an event,
or that an event was generated but something strange was detected
in it.event
object, of type Event
, which is a public member of
pythia
. You therefore have access to all the tools described
on the pages under the "Study Output" header in the index. For instance,
an event can be listed with
pythia.event.list()
, the identity of the pythia.event[i].id()
, and so on.process
, also of type
Event
.pythia.statistics();to get some run statistics, on cross sections and the number of errors and warnings encountered.
Vec4
class gives an implementation of four-vectors.
The member function names are based on the assumption that these
represent momentum vectors. Thus one can get or set
px()
, py()
, pz()
and
()e
, but not x, y, z or t. (When
production vertices are defined in the particle class, this is
partly circumvented by new methods that hide a Vec4
.)
Derived quantities like the pT()
, the pAbs()
,
and the theta()
and phi()
angles
may be read out. The names should be self-explanatory, so we refer
to the header class.
A set of overloaded operators are defined for four-vectors, so that
one may naturally add, subtract, multiply or divide four-vectors with
each other or with double numbers, for all the cases that are
meaningful.
The Particle
object contains a Vec4 p
that
stores the particle four-momentum, and another Vec4 vProd
for the production vertex. Therefore a user would not normally access the
Vec4
class directly, but by using the similarly-named methods
of the Particle
class. (The latter also stores the particle mass
separately, offering an element of redundancy, helpful in avoiding some
roundoff errors.) However, for simple analysis tasks it may be convenient
to use Vec4
, e.g., to define the four-vector sum of a set of
particles.
Simple rotations and boosts of the four-vectors are easily obtained
with member functions. For a longer sequence of rotations and boosts,
and where several Vec4
are involved for the same set of
operations, a more efficient approach is to define a
RotBstMatrix
, which forms a separate auxiliary class.
This matrix can be built up from the successive set of operations to be
performed and, once defined, it can be applied on as many
Vec4
as required.
pythia8052/doc/PaxHeaders.94862/FourVector.xml 0000644 0000000 0000000 00000000074 13466754007 016052 x ustar 00 30 atime=1619556553.868645794
30 ctime=1619556553.868645794
pythia8052/doc/FourVector.xml 0000644 0002170 0000144 00000004103 13466754007 017004 0 ustar 00pythia users 0000000 0000000 Vec4
class gives an implementation of four-vectors.
The member function names are based on the assumption that these
represent momentum vectors. Thus one can get or set
px()
, py()
, pz()
and
()e
, but not Vec4
.)
Derived quantities like the pT()
, the pAbs()
,
and the theta()
and phi()
angles
may be read out. The names should be self-explanatory, so we refer
to the header class.
A set of overloaded operators are defined for four-vectors, so that
one may naturally add, subtract, multiply or divide four-vectors with
each other or with double numbers, for all the cases that are
meaningful.
The Particle
object contains a Vec4 p
that
stores the particle four-momentum, and another Vec4 vProd
for the production vertex. Therefore a user would not normally access the
Vec4
class directly, but by using the similarly-named methods
of the Particle
class. (The latter also stores the particle mass
separately, offering an element of redundancy, helpful in avoiding some
roundoff errors.) However, for simple analysis tasks it may be convenient
to use Vec4
, e.g., to define the four-vector sum of a set of
particles.
Simple rotations and boosts of the four-vectors are easily obtained
with member functions. For a longer sequence of rotations and boosts,
and where several Vec4
are involved for the same set of
operations, a more efficient approach is to define a
RotBstMatrix
, which forms a separate auxiliary class.
This matrix can be built up from the successive set of operations to be
performed and, once defined, it can be applied on as many
Vec4
as required.
StringZ
class handles the choice of longitudinal
lightcone fraction z according to one of two possible
shape sets.
The Lund symmetric fragmentation function [And83] is the
only alternative for light quarks. It is of the form
StringPT
class handles the choice of fragmentation
pT. At each string breaking the quark and antiquark of the pair are
supposed to receive opposite and compensating pT kicks according
to a Gaussian distribution in p_x and p_y separately.
Call sigma_q the width of the p_x and p_y
distributions separately, i.e.
MiniStringFragmentation
machinery is only used when a
string system has so small invariant mass that normal string fragmentation
is difficult/impossible. Instead one or two particles are produced,
in the former case shuffling energy-momentum relative to another
colour singlet system in the event, while preserving the invariant
mass of that system. With one exception parameters are the same as
defined for normal string fragmentation, to the extent that they are
at all applicable in this case.
A discussion of the relevant physics is found in [Nor00].
The current implementation does not completely abide to the scheme
presented there, however, but has in part been simplified. (In part
for greater clarity, in part since the class is not quite finished yet.)
mode name="MiniStringFragmentation:nTry" default="2" min="1" max="10"
StringZ
class handles the choice of longitudinal
lightcone fraction StringPT
class handles the choice of fragmentation
MiniStringFragmentation
machinery is only used when a
string system has so small invariant mass that normal string fragmentation
is difficult/impossible. Instead one or two particles are produced,
in the former case shuffling energy-momentum relative to another
colour singlet system in the event, while preserving the invariant
mass of that system. With one exception parameters are the same as
defined for normal string fragmentation, to the extent that they are
at all applicable in this case.
A discussion of the relevant physics is found in Nor00.
The current implementation does not completely abide to the scheme
presented there, however, but has in part been simplified. (In part
for greater clarity, in part since the class is not quite finished yet.)
Pythia
class.
Pythia
class, namely which parton densities to use, a choice that then is
propagated through the program. The simplest option is to pick one
of the few distributions available internally:
mode name="Pythia:pPDFset" default="2" min="1" max="2"
Pythia::init
routine is called. If off, the random number
generator is initialized with its default seed at the beginning
of the run, and never again. If on, each new Pythia::init
call (should several be made in the same run) results in the random
number being re-initialized, thereby possibly starting over with the
same sequence, if you do not watch out.
mode name="Pythia:seed" default="-1" max="900000000"
setSeed
is on.pythia.statistics()
at the end of the run.
parameter name="Pythia:epTolerance" default="1e-5"
Pythia
class.
Pythia
class, namely which parton densities to use, a choice that then is
propagated through the program. The simplest option is to pick one
of the few distributions available internally:
Pythia::init
routine is called. If off, the random number
generator is initialized with its default seed at the beginning
of the run, and never again. If on, each new Pythia::init
call (should several be made in the same run) results in the random
number being re-initialized, thereby possibly starting over with the
same sequence, if you do not watch out.
setSeed
is on.pythia.statistics()
at the end of the run.
main11.c
main program.
pythia8052/doc/PaxHeaders.94862/HepMC.xml 0000644 0000000 0000000 00000000074 13466754007 014710 x ustar 00 30 atime=1619556553.884645863
30 ctime=1619556553.884645863
pythia8052/doc/HepMC.xml 0000644 0002170 0000144 00000000675 13466754007 015654 0 ustar 00pythia users 0000000 0000000 main11.c
main program.
Hist
class gives a simple implementation of
one-dimensional histograms, useful for quick-and-dirty testing,
without the need to link to more sophisticated packages.
A Histogram is declared by a
class name="Hist name( title, numberOfBins, xMin, xMax)"
Hist ZpT( "Z0 pT spectrum", 100, 0., 100.);Alternatively you can first declare it and later define it:
Hist ZpT; ZpT.book( "Z0 pT spectrum", 100, 0., 100.);Once declared, its contents can be added by repeated calls to fill method name="fill( xValue, weight)"
ZpT.fill( 22.7, 1.);Since the weight defaults to 1 the last argument could have been omitted in this case. A histogram can be printed by making use of the overloaded << operator, e.g.:
cout << ZpT;A set of overloaded operators have been defined, so that histograms can be added, divided by each other (bin by bin) and so on. Also overloaded operations with double real numbers are available, so that e.g. histograms easily can be rescaled. Thus one may write e.g.
allpT = ZpT + 2. * HpTassuming that
allpT
, ZpT
and HpT
have been booked with the same number of bins and x range. That
responsibility rests on the user; some checks are made for compatibility,
but not enough to catch all possible mistakes.
Some further possibilities are included, like writing out histogram
contents as a table, for plotting e.g. with Gnuplot.
pythia8052/doc/PaxHeaders.94862/Histogram.xml 0000644 0000000 0000000 00000000072 13466754007 015707 x ustar 00 29 atime=1619556553.88864588
29 ctime=1619556553.88864588
pythia8052/doc/Histogram.xml 0000644 0002170 0000144 00000004353 13466754007 016652 0 ustar 00pythia users 0000000 0000000 Hist
class gives a simple implementation of
one-dimensional histograms, useful for quick-and-dirty testing,
without the need to link to more sophisticated packages.
A Histogram is declared by a
Hist ZpT( "Z0 pT spectrum", 100, 0., 100.);Alternatively you can first declare it and later define it:
Hist ZpT; ZpT.book( "Z0 pT spectrum", 100, 0., 100.);Once declared, its contents can be added by repeated calls to fill
ZpT.fill( 22.7, 1.);Since the weight defaults to 1 the last argument could have been omitted in this case. A histogram can be printed by making use of the overloaded << operator, e.g.:
cout << ZpT;A set of overloaded operators have been defined, so that histograms can be added, divided by each other (bin by bin) and so on. Also overloaded operations with double real numbers are available, so that e.g. histograms easily can be rescaled. Thus one may write e.g.
allpT = ZpT + 2. * HpTassuming that
allpT
, ZpT
and HpT
have been booked with the same number of bins and LHAinit
and LHAevnt
classes are base classes,
containing reading and printout functions, plus a pure virtual function
each. Derived classes have to provide these two virtual functions to do
the actual work. Currently the only derived classes are for reading
information from the respective Fortran commonblock or for reading
it from PYTHIA 6.4-produced files.
Normally, pointers to objects of the derived classes should be handed
in with the pythia.init()
method. (If you use the LHA
interface to PYTHIA 6.4, this is taken care of
internally, so no pointers need to be handed in.)
LHAinit
class stores information equivalent to the
/HEPRUP/
commonblock, as required to initialize the event
generation chain. The main difference is that the vector container
now allows a flexible number of subprocesses to be defined. For the
rest, names have been modified, since the 6-character-limit does not
apply, and variables have been regrouped for clarity, but nothing
fundamental.
The pure virtual function set()
has to be implemented in the
derived class, to set relevant information when called. It should
return false
if it fails to set the info.
Inside set()
, such information can be set by the following
methods:
method name="beamA( identity, energy, pdfGroup, pdfSet)"
IDBMUP(1), EBMUP(1), PDFGUP(1), PDFSUP(1)
), and similarly
a beamB
method exists. The parton distribution information
defaults to zero, meaning that internal sets are used.
method name="strategy( choice)"
IDWTUP
).
method name="process( idProcess, crossSection, crossSectionError,
crossSectionMaximum)"
LPRUP, XSECUP, XERRUP,
XMAXUP
).
Each new call will append one more entry to the list of processes.
Information is handed back by the following methods:
method name="idBeamA(), eBeamA(), pdfGroupBeamA(), pdfSetBeamA()"
i
in the range 0 <= i <
size()
.
The information can also be printed using the overloaded
<< operator, e.g. cout << LHAinitObject
.
LHAevnt
class stores information equivalent to the
/HEPEUP/
commonblock, as required to hand in the next
parton-level configuration for complete event generation. The main
difference is that the vector container now allows a flexible number
of partons to be defined. For the rest, names have been modified,
since the 6-character-limit does not apply, and variables have been
regrouped for clarity, but nothing fundamental.
The Les Houches standard is based on Fortran arrays beginning with
index 1, and mother information is defined accordingly. In order to
be compatible with this convention, the zeroth line of the C++ particle
array is kept empty, so that index 1 also here corresponds to the first
particle. One small incompatibility is that the size()
method returns the full size of the particle array, including the
empty zeroth line, and thus is one larger than the true number of
particles (NUP
).
The pure virtual function set()
has to be implemented in
the derived class, to set relevant information when called. It should
return false
if it fails to set the info, e.g. if the
supply of events in a file is exhausted.
Inside set()
, cuch information can be set by the following
methods:
method name="process( idProcess, weight, scale, alphaQED, alphaQCD)"
IDPRUP, XWTGUP, SCALUP, AQEDUP, AQCDUP
). This method
also resets the size of the particle list, and adds the empty zeroth
line, so it has to be called before the particle method below.
method name="particle( id, status, mother1, mother2, colourTag1,
colourTag2, p_x, p_y, p_z, e, m, tau, spin)"
IDUP, ISTUP,
MOTHUP(1,..), MOTHUP(2,..), ICOLUP(1,..), ICOLUP(2,..), PUP(J,..),
VTIMUP, SPINUP
) .
Information is handed back by the following methods:
method name="idProc(), weight(), scale(), alphaQED(), alphaQCD()".
i
in the range
0 <= i < size()
. (But again note that
i = 0
is an empty line, so the true range begins at 1.)
The information can also be printed using the overloaded
<< operator, e.g. cout << LHAevntObject
.
LHAinitFortran
class derives from LHAinit
.
It reads initialization information from the Les Houches standard
Fortran commonblock, assuming this commonblock behaves like an
extern "C" struct
named heprup_
. (Note the final
underscore, to match how the gcc compiler internally names Fortran
files.)
Initialization is with
LHAinitFortran lhaInit();i.e. does not require any arguments. The
LHAevntFortran
class derives from LHAevnt
.
It reads information on the next event, stored in the Les Houches
standard Fortran commonblock, assuming this commonblock behaves like
an extern "C" struct
named hepeup_
.
Initialization is with
LHAevntFortran lhaEvnt();i.e. does not require any arguments. See further the page on "Accessing PYTHIA 6 Processes" for information how PYTHIA 6.4 can be linked to make use of this facility. Several of the example main programs illustrate how it can be used.
LHAinitPythia6
class derives from LHAinit
.
It reads initialization information, written from PYTHIA 6.4
with the PYUPIN
routine, from a given file. The file name
should be given as argument at the instantiation, e.g.
LHAinitPythia6 lhaInit("ttsample.init");Thw
LHAevntPythia6
class derives from LHAevnt
.
It reads information on the next event, written from PYTHIA 6.4 with
the PYUPEV
routine, from a given file. The file name should be
given as argument at the instantiation, e.g.
LHAevntPythia6 lhaEvnt("ttsample.evnt");The
PYUPEV
routine writes mother indices using the Fortran
(and Les Houches) convention of the first particle being number 1,
and so matches the C++ "fix" of an empty zeroth particle.
An example how to use events in such exernal files is found in
main02.cc
.
strategy = 10
(not present in the IDWTUP
specification) has been added. It takes a given partonic input,
no questions asked, and hadronizes it, i.e. does string fragmentation
and decay. Thereby the normal process-level and parton-level machineries
are bypassed, to the largest extent possible. (Some parts are used,
e.g. first to translate the Les Houches event to the process record
and later to the event record.) Such an option can therefore be used
to feed in ready-made parton-level configurations, without needing to
specify where these come from, i.e. there need be no beams or any such
explicit information, but of course the user must have taken care of it
beforehand.
An example how this can be used for toy studies is found in
main03.cc
.
pythia8052/doc/PaxHeaders.94862/LesHouches.xml 0000644 0000000 0000000 00000000074 13466754007 016016 x ustar 00 30 atime=1619556553.900645932
30 ctime=1619556553.900645932
pythia8052/doc/LesHouches.xml 0000644 0002170 0000144 00000021356 13466754007 016761 0 ustar 00pythia users 0000000 0000000 LHAinit
and LHAevnt
classes are base classes,
containing reading and printout functions, plus a pure virtual function
each. Derived classes have to provide these two virtual functions to do
the actual work. Currently the only derived classes are for reading
information from the respective Fortran commonblock or for reading
it from PYTHIA 6.4-produced files.
Normally, pointers to objects of the derived classes should be handed
in with the pythia.init()
method. (If you use the LHA
interface to PYTHIA 6.4, this is taken care of
internally, so no pointers need to be handed in.)
LHAinit
class stores information equivalent to the
/HEPRUP/
commonblock, as required to initialize the event
generation chain. The main difference is that the vector container
now allows a flexible number of subprocesses to be defined. For the
rest, names have been modified, since the 6-character-limit does not
apply, and variables have been regrouped for clarity, but nothing
fundamental.
The pure virtual function set()
has to be implemented in the
derived class, to set relevant information when called. It should
return false
if it fails to set the info.
Inside set()
, such information can be set by the following
methods:
IDBMUP(1), EBMUP(1), PDFGUP(1), PDFSUP(1)
), and similarly
a beamB
method exists. The parton distribution information
defaults to zero, meaning that internal sets are used.
IDWTUP
).
LPRUP, XSECUP, XERRUP,
XMAXUP
).
Each new call will append one more entry to the list of processes.
i
in the range 0 <= i <
size()
.
cout << LHAinitObject
.
LHAevnt
class stores information equivalent to the
/HEPEUP/
commonblock, as required to hand in the next
parton-level configuration for complete event generation. The main
difference is that the vector container now allows a flexible number
of partons to be defined. For the rest, names have been modified,
since the 6-character-limit does not apply, and variables have been
regrouped for clarity, but nothing fundamental.
The Les Houches standard is based on Fortran arrays beginning with
index 1, and mother information is defined accordingly. In order to
be compatible with this convention, the zeroth line of the C++ particle
array is kept empty, so that index 1 also here corresponds to the first
particle. One small incompatibility is that the size()
method returns the full size of the particle array, including the
empty zeroth line, and thus is one larger than the true number of
particles (NUP
).
The pure virtual function set()
has to be implemented in
the derived class, to set relevant information when called. It should
return false
if it fails to set the info, e.g. if the
supply of events in a file is exhausted.
Inside set()
, cuch information can be set by the following
methods:
IDPRUP, XWTGUP, SCALUP, AQEDUP, AQCDUP
). This method
also resets the size of the particle list, and adds the empty zeroth
line, so it has to be called before the particle method below.
IDUP, ISTUP,
MOTHUP(1,..), MOTHUP(2,..), ICOLUP(1,..), ICOLUP(2,..), PUP(J,..),
VTIMUP, SPINUP
) .
i
in the range
0 <= i < size()
. (But again note that
i = 0
is an empty line, so the true range begins at 1.)
cout << LHAevntObject
.
LHAinitFortran
class derives from LHAinit
.
It reads initialization information from the Les Houches standard
Fortran commonblock, assuming this commonblock behaves like an
extern "C" struct
named heprup_
. (Note the final
underscore, to match how the gcc compiler internally names Fortran
files.)
Initialization is with
LHAinitFortran lhaInit();i.e. does not require any arguments. The
LHAevntFortran
class derives from LHAevnt
.
It reads information on the next event, stored in the Les Houches
standard Fortran commonblock, assuming this commonblock behaves like
an extern "C" struct
named hepeup_
.
Initialization is with
LHAevntFortran lhaEvnt();i.e. does not require any arguments. See further the page on "Accessing PYTHIA 6 Processes" for information how PYTHIA 6.4 can be linked to make use of this facility. Several of the example main programs illustrate how it can be used.
LHAinitPythia6
class derives from LHAinit
.
It reads initialization information, written from PYTHIA 6.4
with the PYUPIN
routine, from a given file. The file name
should be given as argument at the instantiation, e.g.
LHAinitPythia6 lhaInit("ttsample.init");Thw
LHAevntPythia6
class derives from LHAevnt
.
It reads information on the next event, written from PYTHIA 6.4 with
the PYUPEV
routine, from a given file. The file name should be
given as argument at the instantiation, e.g.
LHAevntPythia6 lhaEvnt("ttsample.evnt");The
PYUPEV
routine writes mother indices using the Fortran
(and Les Houches) convention of the first particle being number 1,
and so matches the C++ "fix" of an empty zeroth particle.
An example how to use events in such exernal files is found in
main02.cc
.
strategy = 10
(not present in the IDWTUP
specification) has been added. It takes a given partonic input,
no questions asked, and hadronizes it, i.e. does string fragmentation
and decay. Thereby the normal process-level and parton-level machineries
are bypassed, to the largest extent possible. (Some parts are used,
e.g. first to translate the Les Houches event to the process record
and later to the event record.) Such an option can therefore be used
to feed in ready-made parton-level configurations, without needing to
specify where these come from, i.e. there need be no beams or any such
explicit information, but of course the user must have taken care of it
beforehand.
An example how this can be used for toy studies is found in
main03.cc
.
Settings
machinery. They can thus be put among
the other cards, and then read back in by the main program.
numberOfEvents/numberToShow
events.
mode name="Main:timesAllowErrors" default="10" min = "0"
pythia.next()
returns false,
i.e. that an event is flawed, before aborting the run.
flag name="Main:showChangedSettings" default="true"
id
code for the first incoming particle.
mode name="Main:idBeamB" default="2212"
id
code for the second incoming particle.
flag name="Main:inCMframe" default="true"
Main:inCMframe
is true.
parameter name="Main:eBeamA" default="7000." min="0."
main01.cmnd
"cards file" illustrating many of the
parameters listed on these pages.ttsample.init
and
ttsample.evnt
. These currently only contain 100 events,
so this program is mainly a demonstration of principles.
The Fortran program main02for.f
has been used to generate
the input files. This program can be modified and linked to
Pythia6.o
to generate other files, always in pairs.main07.cmnd
"cards file" illustrating many of the
parameters listed below. (Brief example given in talks.)ttsample.init
and ttsample.evnt
. These again were generated with main02.f.
(Brief example given in talks.)CellJet
cone jet finder.Settings
machinery. They can thus be put among
the other cards, and then read back in by the main program.
numberOfEvents/numberToShow
events.
pythia.next()
returns false,
i.e. that an event is flawed, before aborting the run.
id
code for the first incoming particle.
id
code for the second incoming particle.
Main:inCMframe
is true.
main01.cmnd
"cards file" illustrating many of the
parameters listed on these pages.ttsample.init
and
ttsample.evnt
. These currently only contain 100 events,
so this program is mainly a demonstration of principles.
The Fortran program main02for.f
has been used to generate
the input files. This program can be modified and linked to
Pythia6.o
to generate other files, always in pairs.main07.cmnd
"cards file" illustrating many of the
parameters listed below. (Brief example given in talks.)ttsample.init
and ttsample.evnt
. These again were generated with main02.f.
(Brief example given in talks.)CellJet
cone jet finder.MultipleInteractions
class.
In addition there is the possibility of a global rescaling of
cross sections (which could not easily be accommodated by a
changed alpha_strong, since alpha_strong runs)
parameter name="MultipleInteractions:Kfactor" default="1.0" min="0.5"
max="4.0"
MultipleInteractions
class.
In addition there is the possibility of a global rescaling of
cross sections (which could not easily be accommodated by a
changed ProcessLevel
class administrates the initial step of
the event generation, wherein the basic process is selected. Currently
this is done either using some of the internal processes, or with
Les Houches Accord input, or from the Fortran PYTHIA 6.4 code,
the latter also using the Les Houches Accord standard.
Since there cannot be any event at all without an initial process,
there is no possibility to switch off this part of the story. It is
possible, however, to stop the generation immediately after the
basic process has been selected:
flag name="Pythia:afterProcessLevel" default="on"
process
record is filled, but the event
one is not.
PartonLevel
class administrates the middle step of the
event generation, i.e. the evolution from an input (hard) process from
ProcessLevel
, containing a few partons only, to a complete
parton-level configuration to be handed on to HadronLevel
.
This step involves the application of initial- and final-state radiation,
multiple interactions and the structure of beam remnants.
Some parts of the event generation on this level may be switched off
individually:
flag name="PartonLevel:MI" default="on"
HadronLevel
class administrates the final step of the
event generation, wherein the partonic configuration from
PartonLevel
is hadronized, including string fragmentation
and secondary decays.
Most of the code in this class deals with subdividing the partonic
content of the event into separate colour singlets, that can be
treated individually by the string fragmentation machinery. When a
junction and an antijunction are directly connected, it also breaks
the string between the two, so that the topology can be reduced back
to two separate one-junction systems, while still preserving the
expected particle flow in the junction-junction string region(s).
Some parts of the event generation on this level may be switched off
individually:
flag name="HadronLevel:Hadronize" default="on"
ProcessLevel
class administrates the initial step of
the event generation, wherein the basic process is selected. Currently
this is done either using some of the internal processes, or with
Les Houches Accord input, or from the Fortran PYTHIA 6.4 code,
the latter also using the Les Houches Accord standard.
Since there cannot be any event at all without an initial process,
there is no possibility to switch off this part of the story. It is
possible, however, to stop the generation immediately after the
basic process has been selected:
process
record is filled, but the event
one is not.
PartonLevel
class administrates the middle step of the
event generation, i.e. the evolution from an input (hard) process from
ProcessLevel
, containing a few partons only, to a complete
parton-level configuration to be handed on to HadronLevel
.
This step involves the application of initial- and final-state radiation,
multiple interactions and the structure of beam remnants.
Some parts of the event generation on this level may be switched off
individually:
HadronLevel
class administrates the final step of the
event generation, wherein the partonic configuration from
PartonLevel
is hadronized, including string fragmentation
and secondary decays.
Most of the code in this class deals with subdividing the partonic
content of the event into separate colour singlets, that can be
treated individually by the string fragmentation machinery. When a
junction and an antijunction are directly connected, it also breaks
the string between the two, so that the topology can be reduced back
to two separate one-junction systems, while still preserving the
expected particle flow in the junction-junction string region(s).
Some parts of the event generation on this level may be switched off
individually:
modeME()
above, almost agrees with the Pythia6 ones, but several obsolete have
been removed, and a few moved for better consistency. Here is the list
of currently allowed codes:
code
- 10 number
of particles in the decay channelcode
- 20particle: 1 d dbar -1 1 0.33000 0.00000 0.00000 0.00000E+00 0 particle: 2 u ubar 2 1 0.33000 0.00000 0.00000 0.00000E+00 0 particle: 3 s sbar -1 1 0.50000 0.00000 0.00000 0.00000E+00 0 particle: 4 c cbar 2 1 1.50000 0.00000 0.00000 0.00000E+00 0 particle: 5 b bbar -1 1 4.80000 0.00000 0.00000 0.00000E+00 0 particle: 6 t tbar 2 1 175.00000 1.39816 13.98156 0.00000E+00 1 particle: 7 b' b'bar -1 1 400.00000 0.00000 0.00000 0.00000E+00 1 particle: 8 t' t'bar 2 1 400.00000 0.00000 0.00000 0.00000E+00 1 particle: 11 e- e+ -3 0 0.00051 0.00000 0.00000 0.00000E+00 0 particle: 12 nu_e nu_ebar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 13 mu- mu+ -3 0 0.10566 0.00000 0.00000 6.58654E+05 0 channel: 1.0000000 42 -12 11 14 particle: 14 nu_mu nu_mubar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 15 tau- tau+ -3 0 1.77700 0.00000 0.00000 8.72000E-02 1 channel: 0.1783000 42 -12 11 16 channel: 0.1735000 42 -14 13 16 channel: 0.1131000 0 16 -211 channel: 0.2494000 0 16 -213 channel: 0.0030000 41 16 -211 111 channel: 0.0900000 41 16 -213 111 channel: 0.0027000 41 16 -211 111 111 channel: 0.0100000 41 16 -213 111 111 channel: 0.0014000 41 16 -211 111 111 111 channel: 0.0012000 41 16 -213 111 111 111 channel: 0.0002500 41 16 -211 310 channel: 0.0002500 41 16 -211 130 channel: 0.0071000 0 16 -321 channel: 0.0120000 0 16 -323 channel: 0.0004000 41 16 -321 111 channel: 0.0007500 41 16 -323 111 channel: 0.0000600 41 16 -323 111 111 channel: 0.0007800 41 16 -321 310 channel: 0.0007800 41 16 -321 130 channel: 0.0034000 41 16 -321 321 -211 channel: 0.0800000 41 16 -211 113 channel: 0.0110000 41 16 -211 211 -211 channel: 0.0191000 41 16 -211 223 channel: 0.0000600 41 16 -211 221 channel: 0.0050000 41 16 -213 113 channel: 0.0133000 41 16 -211 113 111 channel: 0.0067000 41 16 -213 211 -211 channel: 0.0005000 41 16 -211 211 -211 111 channel: 0.0035000 41 16 -213 223 channel: 0.0006000 41 16 -211 223 111 channel: 0.0015000 41 16 -213 221 channel: 0.0002100 41 16 -211 221 111 channel: 0.0002000 41 16 -213 113 111 channel: 0.0007500 41 16 -211 113 113 channel: 0.0001000 41 16 -211 221 221 channel: 0.0002000 41 16 -211 113 111 111 channel: 0.0011000 41 16 -213 113 111 111 channel: 0.0002000 41 16 -211 213 -213 channel: 0.0002000 41 16 -211 213 -211 111 channel: 0.0002000 41 16 -211 -213 211 111 channel: 0.0002200 41 16 -211 113 113 111 channel: 0.0004000 41 16 -323 111 111 channel: 0.0001000 41 16 -321 111 111 111 channel: 0.0020500 41 16 -211 310 111 channel: 0.0020500 41 16 -211 130 111 channel: 0.0006900 41 16 -321 310 111 channel: 0.0006900 41 16 -321 130 111 channel: 0.0002500 41 16 -211 310 310 channel: 0.0005100 41 16 -211 310 130 channel: 0.0002500 41 16 -211 130 130 particle: 16 nu_tau nu_taubar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 17 tau'- tau'+ -3 0 400.00000 0.00000 0.00000 0.00000E+00 1 particle: 18 nu'_tau nu'_taubar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 21 g void 0 2 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 22 gamma void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 23 Z0 void 0 0 91.18800 2.47813 24.78129 0.00000E+00 1 particle: 24 W+ W- 3 0 80.45000 2.07115 20.71149 0.00000E+00 1 particle: 25 h0 void 0 0 115.00000 0.00367 0.03669 0.00000E+00 1 particle: 32 Z'0 void 0 0 500.00000 14.54029 145.40294 0.00000E+00 1 particle: 33 Z"0 void 0 0 900.00000 0.00000 0.00000 0.00000E+00 0 particle: 34 W'+ W'- 3 0 500.00000 16.66099 166.60993 0.00000E+00 1 particle: 35 H0 void 0 0 300.00000 8.38842 83.88423 0.00000E+00 1 particle: 36 A0 void 0 0 300.00000 3.37520 33.75195 0.00000E+00 1 particle: 37 H+ H- 3 0 300.00000 4.17669 41.76694 0.00000E+00 1 particle: 39 Graviton void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 41 R0 Rbar0 0 0 5000.00000 417.29147 4172.91467 0.00000E+00 1 particle: 42 LQ_ue LQ_uebar -1 1 200.00000 0.39162 3.91621 0.00000E+00 1 particle: 81 specflav void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 82 rndmflav rndmflavbar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 88 junction void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 90 system void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 91 cluster void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 92 string void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 93 indep. void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 94 CMshower void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 95 SPHEaxis void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 96 THRUaxis void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 97 CLUSjet void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 98 CELLjet void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 99 table void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 110 reggeon void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 111 pi0 void 0 0 0.13498 0.00000 0.00000 3.00000E-05 1 channel: 0.9880000 0 22 22 channel: 0.0120000 2 22 11 -11 particle: 113 rho0 void 0 0 0.76850 0.15100 0.40000 0.00000E+00 1 channel: 0.9987390 3 211 -211 channel: 0.0007900 0 111 22 channel: 0.0003800 0 221 22 channel: 0.0000460 0 13 -13 channel: 0.0000450 0 11 -11 particle: 115 a_20 void 0 0 1.31800 0.10700 0.25000 0.00000E+00 1 channel: 0.3472500 0 213 -211 channel: 0.3472500 0 -213 211 channel: 0.1440000 0 221 111 channel: 0.1040000 0 223 211 -211 channel: 0.0245000 0 321 -321 channel: 0.0122500 0 130 130 channel: 0.0122500 0 310 310 channel: 0.0028000 0 111 22 channel: 0.0057000 0 331 111 particle: 130 K_L0 void 0 0 0.49767 0.00000 0.00000 1.55000E+04 0 channel: 0.2112000 0 111 111 111 channel: 0.1256000 0 211 -211 111 channel: 0.1939000 42 -12 11 211 channel: 0.1939000 42 12 -11 -211 channel: 0.1359000 42 -14 13 211 channel: 0.1359000 42 14 -13 -211 channel: 0.0020000 0 211 -211 channel: 0.0010000 0 111 111 channel: 0.0006000 0 22 22 particle: 211 pi+ pi- 3 0 0.13957 0.00000 0.00000 7.80450E+03 0 channel: 0.9998770 0 -13 14 channel: 0.0001230 0 -11 12 particle: 213 rho+ rho- 3 0 0.76690 0.14900 0.40000 0.00000E+00 1 channel: 0.9995500 3 211 111 channel: 0.0004500 0 211 22 particle: 215 a_2+ a_2- 3 0 1.31800 0.10700 0.25000 0.00000E+00 1 channel: 0.3472500 0 213 111 channel: 0.3472500 0 113 211 channel: 0.1440000 0 221 211 channel: 0.1040000 0 223 211 111 channel: 0.0490000 0 321 -311 channel: 0.0028000 0 211 22 channel: 0.0057000 0 331 211 particle: 221 eta void 0 0 0.54745 0.00000 0.00000 0.00000E+00 1 channel: 0.3923000 0 22 22 channel: 0.3210000 0 111 111 111 channel: 0.2317000 0 211 -211 111 channel: 0.0478000 0 22 211 -211 channel: 0.0049000 2 22 11 -11 channel: 0.0013000 0 211 -211 11 -11 channel: 0.0003000 0 22 13 -13 channel: 0.0007000 0 111 22 22 particle: 223 omega void 0 0 0.78194 0.00843 0.10000 0.00000E+00 1 channel: 0.8900000 1 211 -211 111 channel: 0.0869300 0 22 111 channel: 0.0221000 3 211 -211 channel: 0.0008300 0 221 22 channel: 0.0000700 0 111 111 22 channel: 0.0000700 0 11 -11 particle: 225 f_2 void 0 0 1.27500 0.18500 0.17000 0.00000E+00 1 channel: 0.5640000 0 211 -211 channel: 0.2820000 0 111 111 channel: 0.0720000 0 211 -211 111 111 channel: 0.0280000 0 211 -211 211 -211 channel: 0.0230000 0 321 -321 channel: 0.0115000 0 130 130 channel: 0.0115000 0 310 310 channel: 0.0050000 0 221 221 channel: 0.0030000 0 111 111 111 111 particle: 310 K_S0 void 0 0 0.49767 0.00000 0.00000 2.67620E+01 1 channel: 0.6861000 0 211 -211 channel: 0.3139000 0 111 111 particle: 311 K0 Kbar0 0 0 0.49767 0.00000 0.00000 0.00000E+00 1 channel: 0.5000000 0 130 channel: 0.5000000 0 310 particle: 313 K*0 K*bar0 0 0 0.89610 0.05050 0.20000 0.00000E+00 1 channel: 0.6650000 3 321 -211 channel: 0.3330000 3 311 111 channel: 0.0020000 0 311 22 particle: 315 K*_20 K*_2bar0 0 0 1.43200 0.10900 0.12000 0.00000E+00 1 channel: 0.3330000 0 321 -211 channel: 0.1660000 0 311 111 channel: 0.1680000 0 323 -211 channel: 0.0840000 0 313 111 channel: 0.0870000 0 323 -211 111 channel: 0.0430000 0 313 211 -211 channel: 0.0590000 0 321 -213 channel: 0.0290000 0 311 113 channel: 0.0290000 0 311 223 channel: 0.0020000 0 311 221 particle: 321 K+ K- 3 0 0.49360 0.00000 0.00000 3.70900E+03 0 channel: 0.6352000 0 -13 14 channel: 0.2116000 0 211 111 channel: 0.0559000 0 211 211 -211 channel: 0.0173000 0 211 111 111 channel: 0.0482000 42 12 -11 111 channel: 0.0318000 42 14 -13 111 particle: 323 K*+ K*- 3 0 0.89160 0.04980 0.20000 0.00000E+00 1 channel: 0.6660000 3 311 211 channel: 0.3330000 3 321 111 channel: 0.0010000 0 321 22 particle: 325 K*_2+ K*_2- 3 0 1.42500 0.09800 0.12000 0.00000E+00 1 channel: 0.3320000 0 311 211 channel: 0.1660000 0 321 111 channel: 0.1680000 0 313 211 channel: 0.0840000 0 323 111 channel: 0.0860000 0 313 211 111 channel: 0.0430000 0 323 211 -211 channel: 0.0590000 0 311 213 channel: 0.0290000 0 321 113 channel: 0.0290000 0 321 223 channel: 0.0020000 0 321 221 channel: 0.0020000 0 321 22 particle: 331 eta' void 0 0 0.95777 0.00020 0.00200 0.00000E+00 1 channel: 0.4370000 0 211 -211 221 channel: 0.2080000 0 111 111 221 channel: 0.3020000 0 22 113 channel: 0.0302000 0 22 223 channel: 0.0212000 0 22 22 channel: 0.0016000 0 111 111 111 particle: 333 phi void 0 0 1.01940 0.00443 0.01500 0.00000E+00 1 channel: 0.4894700 3 321 -321 channel: 0.3400000 3 130 310 channel: 0.0430000 0 -213 211 channel: 0.0430000 0 113 111 channel: 0.0430000 0 213 -211 channel: 0.0270000 1 211 -211 111 channel: 0.0126000 0 22 221 channel: 0.0013000 0 111 22 channel: 0.0003000 0 11 -11 channel: 0.0002500 0 13 -13 channel: 0.0000800 0 211 -211 particle: 335 f'_2 void 0 0 1.52500 0.07600 0.20000 0.00000E+00 1 channel: 0.4440000 0 321 -321 channel: 0.2220000 0 130 130 channel: 0.2220000 0 310 310 channel: 0.1040000 0 221 221 channel: 0.0040000 0 211 -211 channel: 0.0040000 0 111 111 particle: 411 D+ D- 3 0 1.86930 0.00000 0.00000 3.17000E-01 1 channel: 0.0700000 42 -11 12 -311 channel: 0.0650000 42 -11 12 -313 channel: 0.0050000 42 -11 12 -311 111 channel: 0.0050000 42 -11 12 -321 211 channel: 0.0110000 42 -11 12 -313 111 channel: 0.0110000 42 -11 12 -323 211 channel: 0.0010000 42 -11 12 111 channel: 0.0010000 42 -11 12 221 channel: 0.0010000 42 -11 12 331 channel: 0.0010000 42 -11 12 113 channel: 0.0010000 42 -11 12 223 channel: 0.0700000 42 -13 14 -311 channel: 0.0650000 42 -13 14 -313 channel: 0.0050000 42 -13 14 -311 111 channel: 0.0050000 42 -13 14 -321 211 channel: 0.0110000 42 -13 14 -313 111 channel: 0.0110000 42 -13 14 -323 211 channel: 0.0010000 42 -13 14 111 channel: 0.0010000 42 -13 14 221 channel: 0.0010000 42 -13 14 331 channel: 0.0010000 42 -13 14 113 channel: 0.0010000 42 -13 14 223 channel: 0.0260000 0 -311 211 channel: 0.0190000 0 -313 211 channel: 0.0660000 0 -311 213 channel: 0.0410000 0 -313 213 channel: 0.0450000 0 -20313 211 channel: 0.0760000 0 -311 20213 channel: 0.0073000 0 -311 321 channel: 0.0047000 0 -313 321 channel: 0.0047000 0 -311 323 channel: 0.0260000 0 -313 323 channel: 0.0010000 0 111 211 channel: 0.0006000 0 111 213 channel: 0.0066000 0 221 211 channel: 0.0050000 0 221 213 channel: 0.0030000 0 331 211 channel: 0.0030000 0 331 213 channel: 0.0006000 0 113 211 channel: 0.0006000 0 113 213 channel: 0.0010000 0 223 211 channel: 0.0010000 0 223 213 channel: 0.0060000 0 333 211 channel: 0.0050000 0 333 213 channel: 0.0120000 0 -311 211 111 channel: 0.0057000 0 -313 211 113 channel: 0.0670000 0 -321 211 211 channel: 0.0080000 0 -321 213 211 channel: 0.0022000 0 211 211 -211 channel: 0.0270000 0 -311 321 -311 channel: 0.0040000 0 -321 321 211 channel: 0.0190000 0 333 211 111 channel: 0.0120000 0 -311 211 211 -211 channel: 0.0020000 0 -313 211 211 -211 channel: 0.0090000 0 -321 211 211 111 channel: 0.0218000 0 211 211 -211 111 channel: 0.0010000 0 -321 211 211 211 -211 channel: 0.0220000 0 -321 211 211 111 111 channel: 0.0870000 0 -311 211 211 -211 111 channel: 0.0010000 0 -311 113 211 211 -211 channel: 0.0019000 0 -321 113 211 211 111 channel: 0.0015000 0 211 211 211 -211 -211 channel: 0.0028000 0 113 211 211 -211 111 particle: 413 D*+ D*- 3 0 2.01000 0.00000 0.00000 0.00000E+00 1 channel: 0.6830000 3 421 211 channel: 0.3060000 3 411 111 channel: 0.0110000 0 411 22 particle: 415 D*_2+ D*_2- 3 0 2.46000 0.02300 0.12000 0.00000E+00 1 channel: 0.3000000 0 421 211 channel: 0.1500000 0 411 111 channel: 0.1600000 0 423 211 channel: 0.0800000 0 413 111 channel: 0.1300000 0 423 211 111 channel: 0.0600000 0 413 211 -211 channel: 0.0800000 0 421 211 111 channel: 0.0400000 0 411 211 -211 particle: 421 D0 Dbar0 0 0 1.86450 0.00000 0.00000 1.24400E-01 1 channel: 0.0340000 42 -11 12 -321 channel: 0.0270000 42 -11 12 -323 channel: 0.0020000 42 -11 12 -311 -211 channel: 0.0020000 42 -11 12 -321 111 channel: 0.0040000 42 -11 12 -313 -211 channel: 0.0040000 42 -11 12 -323 111 channel: 0.0020000 42 -11 12 -211 channel: 0.0020000 42 -11 12 -213 channel: 0.0340000 42 -13 14 -321 channel: 0.0270000 42 -13 14 -323 channel: 0.0020000 42 -13 14 -311 -211 channel: 0.0020000 42 -13 14 -321 111 channel: 0.0040000 42 -13 14 -313 -211 channel: 0.0040000 42 -13 14 -323 111 channel: 0.0020000 42 -13 14 -211 channel: 0.0020000 42 -13 14 -213 channel: 0.0365000 0 -321 211 channel: 0.0450000 0 -323 211 channel: 0.0730000 0 -321 213 channel: 0.0620000 0 -323 213 channel: 0.0210000 0 -311 111 channel: 0.0210000 0 -313 111 channel: 0.0210000 0 -313 221 channel: 0.0061000 0 -311 113 channel: 0.0150000 0 -313 113 channel: 0.0250000 0 -311 223 channel: 0.0088000 0 -311 333 channel: 0.0740000 0 -321 20213 channel: 0.0109000 0 -10323 211 channel: 0.0041000 0 -321 321 channel: 0.0020000 0 -323 321 channel: 0.0035000 0 -321 323 channel: 0.0011000 0 -311 311 channel: 0.0010000 0 -313 311 channel: 0.0027000 0 -313 313 channel: 0.0016000 0 211 -211 channel: 0.0016000 0 111 111 channel: 0.0018000 0 333 113 channel: 0.0110000 0 -321 211 111 channel: 0.0063000 0 -321 211 113 channel: 0.0052000 0 -321 321 -311 channel: 0.0180000 0 -311 211 -211 channel: 0.0160000 0 -313 211 -211 channel: 0.0034000 0 -321 311 211 channel: 0.0036000 0 -313 321 -211 channel: 0.0009000 0 310 310 310 channel: 0.0006000 0 333 211 -211 channel: 0.0150000 0 211 -211 111 channel: 0.0923000 0 -321 211 111 111 channel: 0.0180000 0 -321 211 211 -211 channel: 0.0220000 0 -311 211 -211 111 channel: 0.0077000 0 -313 211 -211 111 channel: 0.0090000 0 -311 321 -321 111 channel: 0.0075000 0 211 211 -211 -211 channel: 0.0240000 0 -321 211 211 -211 111 channel: 0.0085000 0 -311 211 211 -211 -211 channel: 0.0670000 0 -311 211 -211 111 111 channel: 0.0511000 0 -311 113 111 111 111 channel: 0.0170000 0 211 211 -211 -211 111 channel: 0.0004000 0 113 211 211 -211 -211 channel: 0.0028000 0 321 -321 211 -211 111 particle: 423 D*0 D*bar0 0 0 2.00670 0.00000 0.00000 0.00000E+00 1 channel: 0.6190000 3 421 111 channel: 0.3810000 0 421 22 particle: 425 D*_20 D*_2bar0 0 0 2.46000 0.02300 0.12000 0.00000E+00 1 channel: 0.3000000 0 411 -211 channel: 0.1500000 0 421 111 channel: 0.1600000 0 413 -211 channel: 0.0800000 0 423 111 channel: 0.1300000 0 413 -211 111 channel: 0.0600000 0 423 211 -211 channel: 0.0800000 0 411 -211 111 channel: 0.0400000 0 421 211 -211 particle: 431 D_s+ D_s- 3 0 1.96850 0.00000 0.00000 1.40000E-01 1 channel: 0.0100000 0 -15 16 channel: 0.0200000 42 -11 12 221 channel: 0.0200000 42 -11 12 331 channel: 0.0300000 42 -11 12 333 channel: 0.0050000 42 -11 12 321 -321 channel: 0.0050000 42 -11 12 311 -311 channel: 0.0200000 42 -13 14 221 channel: 0.0200000 42 -13 14 331 channel: 0.0300000 42 -13 14 333 channel: 0.0050000 42 -13 14 321 -321 channel: 0.0050000 42 -13 14 311 -311 channel: 0.0150000 0 221 211 channel: 0.0370000 0 331 211 channel: 0.0280000 0 333 211 channel: 0.0790000 0 221 213 channel: 0.0950000 0 331 213 channel: 0.0520000 0 333 213 channel: 0.0078000 0 10221 211 channel: 0.0010000 0 211 111 channel: 0.0010000 0 213 111 channel: 0.0010000 0 211 113 channel: 0.0010000 0 213 113 channel: 0.0280000 0 321 -311 channel: 0.0330000 0 323 -311 channel: 0.0260000 0 321 -313 channel: 0.0500000 0 323 -313 channel: 0.0100000 0 2212 -2112 channel: 0.0050000 0 221 321 channel: 0.0050000 0 331 321 channel: 0.0050000 0 333 321 channel: 0.0050000 0 221 323 channel: 0.2500000 13 2 -1 3 -3 channel: 0.0952000 13 2 -1 particle: 433 D*_s+ D*_s- 3 0 2.11240 0.00000 0.00000 0.00000E+00 1 channel: 0.9400000 0 431 22 channel: 0.0600000 0 431 111 particle: 435 D*_2s+ D*_2s- 3 0 2.57350 0.01500 0.05000 0.00000E+00 1 channel: 0.4000000 0 421 321 channel: 0.4000000 0 411 311 channel: 0.1000000 0 423 321 channel: 0.1000000 0 413 311 particle: 441 eta_c void 0 0 2.97980 0.00130 0.00500 0.00000E+00 1 channel: 1.0000000 12 82 -82 particle: 443 J/psi void 0 0 3.09688 0.00000 0.00000 0.00000E+00 1 channel: 0.0602000 0 11 -11 channel: 0.0601000 0 13 -13 channel: 0.8797000 12 82 -82 particle: 445 chi_2c void 0 0 3.55620 0.00200 0.01000 0.00000E+00 1 channel: 0.1350000 0 443 22 channel: 0.8650000 12 82 -82 particle: 511 B0 Bbar0 0 0 5.27920 0.00000 0.00000 4.68000E-01 1 channel: 0.0200000 42 12 -11 -411 channel: 0.0550000 42 12 -11 -413 channel: 0.0050000 42 12 -11 -10413 channel: 0.0050000 42 12 -11 -10411 channel: 0.0080000 42 12 -11 -20413 channel: 0.0120000 42 12 -11 -415 channel: 0.0200000 42 14 -13 -411 channel: 0.0550000 42 14 -13 -413 channel: 0.0050000 42 14 -13 -10413 channel: 0.0050000 42 14 -13 -10411 channel: 0.0080000 42 14 -13 -20413 channel: 0.0120000 42 14 -13 -415 channel: 0.0100000 42 16 -15 -411 channel: 0.0300000 42 16 -15 -413 channel: 0.0035000 0 -411 211 channel: 0.0110000 0 -411 213 channel: 0.0055000 0 -411 20213 channel: 0.0042000 0 -413 211 channel: 0.0090000 0 -413 213 channel: 0.0180000 0 -413 20213 channel: 0.0150000 0 -411 431 channel: 0.0185000 0 -411 433 channel: 0.0135000 0 -413 431 channel: 0.0250000 0 -413 433 channel: 0.0004000 0 441 311 channel: 0.0007000 0 441 313 channel: 0.0008000 0 443 311 channel: 0.0014000 0 443 313 channel: 0.0019000 0 20443 311 channel: 0.0025000 0 20443 313 channel: 0.4291000 43 2 -1 -4 1 channel: 0.0800000 13 2 -4 -1 1 channel: 0.0700000 13 4 -3 -4 1 channel: 0.0200000 13 4 -4 -3 1 channel: 0.0150000 42 2 -1 -2 1 channel: 0.0050000 42 4 -3 -2 1 particle: 513 B*0 B*bar0 0 0 5.32480 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 511 22 particle: 515 B*_20 B*_2bar0 0 0 5.83000 0.02000 0.05000 0.00000E+00 1 channel: 0.3000000 0 521 -211 channel: 0.1500000 0 511 111 channel: 0.1600000 0 523 -211 channel: 0.0800000 0 513 111 channel: 0.1300000 0 523 -211 111 channel: 0.0600000 0 513 211 -211 channel: 0.0800000 0 521 -211 111 channel: 0.0400000 0 511 211 -211 particle: 521 B+ B- 3 0 5.27890 0.00000 0.00000 4.62000E-01 1 channel: 0.0200000 42 12 -11 -421 channel: 0.0550000 42 12 -11 -423 channel: 0.0050000 42 12 -11 -10423 channel: 0.0050000 42 12 -11 -10421 channel: 0.0080000 42 12 -11 -20423 channel: 0.0120000 42 12 -11 -425 channel: 0.0200000 42 14 -13 -421 channel: 0.0550000 42 14 -13 -423 channel: 0.0050000 42 14 -13 -10423 channel: 0.0050000 42 14 -13 -10421 channel: 0.0080000 42 14 -13 -20423 channel: 0.0120000 42 14 -13 -425 channel: 0.0100000 42 16 -15 -421 channel: 0.0300000 42 16 -15 -423 channel: 0.0035000 0 -421 211 channel: 0.0110000 0 -421 213 channel: 0.0055000 0 -421 20213 channel: 0.0042000 0 -423 211 channel: 0.0090000 0 -423 213 channel: 0.0180000 0 -423 20213 channel: 0.0150000 0 -421 431 channel: 0.0185000 0 -421 433 channel: 0.0135000 0 -423 431 channel: 0.0250000 0 -423 433 channel: 0.0004000 0 441 321 channel: 0.0007000 0 441 323 channel: 0.0008000 0 443 321 channel: 0.0014000 0 443 323 channel: 0.0019000 0 20443 321 channel: 0.0025000 0 20443 323 channel: 0.4291000 43 2 -1 -4 2 channel: 0.0800000 13 2 -4 -1 2 channel: 0.0700000 13 4 -3 -4 2 channel: 0.0200000 13 4 -4 -3 2 channel: 0.0150000 42 2 -1 -2 2 channel: 0.0050000 42 4 -3 -2 2 particle: 523 B*+ B*- 3 0 5.32480 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 521 22 particle: 525 B*_2+ B*_2- 3 0 5.83000 0.02000 0.05000 0.00000E+00 1 channel: 0.3000000 0 511 211 channel: 0.1500000 0 521 111 channel: 0.1600000 0 513 211 channel: 0.0800000 0 523 111 channel: 0.1300000 0 513 211 111 channel: 0.0600000 0 523 211 -211 channel: 0.0800000 0 511 211 111 channel: 0.0400000 0 521 211 -211 particle: 531 B_s0 B_sbar0 0 0 5.36930 0.00000 0.00000 4.83000E-01 1 channel: 0.0200000 42 12 -11 -431 channel: 0.0550000 42 12 -11 -433 channel: 0.0050000 42 12 -11 -10433 channel: 0.0050000 42 12 -11 -10431 channel: 0.0080000 42 12 -11 -20433 channel: 0.0120000 42 12 -11 -435 channel: 0.0200000 42 14 -13 -431 channel: 0.0550000 42 14 -13 -433 channel: 0.0050000 42 14 -13 -10433 channel: 0.0050000 42 14 -13 -10431 channel: 0.0080000 42 14 -13 -20433 channel: 0.0120000 42 14 -13 -435 channel: 0.0100000 42 16 -15 -431 channel: 0.0300000 42 16 -15 -433 channel: 0.0035000 0 -431 211 channel: 0.0110000 0 -431 213 channel: 0.0055000 0 -431 20213 channel: 0.0042000 0 -433 211 channel: 0.0090000 0 -433 213 channel: 0.0180000 0 -433 20213 channel: 0.0150000 0 -431 431 channel: 0.0185000 0 -431 433 channel: 0.0135000 0 -433 431 channel: 0.0250000 0 -433 433 channel: 0.0002000 0 441 221 channel: 0.0002000 0 441 331 channel: 0.0007000 0 441 333 channel: 0.0004000 0 443 221 channel: 0.0004000 0 443 331 channel: 0.0014000 0 443 333 channel: 0.0010000 0 20443 221 channel: 0.0009000 0 20443 331 channel: 0.0025000 0 20443 333 channel: 0.4291000 43 2 -1 -4 3 channel: 0.0800000 13 2 -4 -1 3 channel: 0.0700000 13 4 -3 -4 3 channel: 0.0200000 13 4 -4 -3 3 channel: 0.0150000 42 2 -1 -2 3 channel: 0.0050000 42 4 -3 -2 3 particle: 533 B*_s0 B*_sbar0 0 0 5.41630 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 531 22 particle: 535 B*_2s0 B*_2sbar0 0 0 6.07000 0.02000 0.05000 0.00000E+00 1 channel: 0.3000000 0 521 -321 channel: 0.3000000 0 511 -311 channel: 0.2000000 0 523 -321 channel: 0.2000000 0 513 -311 particle: 541 B_c+ B_c- 3 0 6.59400 0.00000 0.00000 1.50000E-01 1 channel: 0.0470000 0 16 -15 channel: 0.1220000 11 4 -3 channel: 0.0060000 11 4 -1 channel: 0.0120000 42 12 -11 441 channel: 0.0350000 42 12 -11 443 channel: 0.0120000 42 14 -13 441 channel: 0.0350000 42 14 -13 443 channel: 0.0030000 42 16 -15 441 channel: 0.0070000 42 16 -15 443 channel: 0.1500000 42 2 -1 -4 4 channel: 0.0370000 42 2 -4 -1 4 channel: 0.0080000 42 2 -3 -4 4 channel: 0.0020000 42 2 -4 -3 4 channel: 0.0500000 42 4 -3 -4 4 channel: 0.0150000 42 4 -4 -3 4 channel: 0.0030000 42 4 -1 -4 4 channel: 0.0010000 42 4 -4 -1 4 channel: 0.0140000 42 -11 12 531 channel: 0.0420000 42 -11 12 533 channel: 0.0140000 42 -13 14 531 channel: 0.0420000 42 -13 14 533 channel: 0.2400000 42 -1 2 3 -5 channel: 0.0650000 42 -1 3 2 -5 channel: 0.0120000 42 -3 2 3 -5 channel: 0.0030000 42 -3 3 2 -5 channel: 0.0010000 42 -11 12 511 channel: 0.0020000 42 -11 12 513 channel: 0.0010000 42 -13 14 511 channel: 0.0020000 42 -13 14 513 channel: 0.0140000 42 -1 2 1 -5 channel: 0.0030000 42 -1 1 2 -5 particle: 543 B*_c+ B*_c- 3 0 6.60200 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 541 22 particle: 545 B*_2c+ B*_2c- 3 0 7.35000 0.02000 0.05000 0.00000E+00 1 channel: 0.3000000 0 511 411 channel: 0.3000000 0 521 421 channel: 0.2000000 0 513 411 channel: 0.2000000 0 523 421 particle: 551 eta_b void 0 0 9.40000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 32 21 21 particle: 553 Upsilon void 0 0 9.46030 0.00000 0.00000 0.00000E+00 1 channel: 0.0252000 0 11 -11 channel: 0.0248000 0 13 -13 channel: 0.0267000 0 15 -15 channel: 0.0150000 32 1 -1 channel: 0.0450000 32 2 -2 channel: 0.0150000 32 3 -3 channel: 0.0450000 32 4 -4 channel: 0.7743000 33 21 21 21 channel: 0.0290000 33 22 21 21 particle: 555 chi_2b void 0 0 9.91320 0.00000 0.00000 0.00000E+00 1 channel: 0.2200000 0 553 22 channel: 0.7800000 32 21 21 particle: 990 pomeron void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 1103 dd_1 dd_1bar -2 -1 0.77133 0.00000 0.00000 0.00000E+00 0 particle: 1114 Delta- Deltabar+ -3 0 1.23400 0.12000 0.14000 0.00000E+00 1 channel: 1.0000000 0 2112 -211 particle: 2101 ud_0 ud_0bar 1 -1 0.57933 0.00000 0.00000 0.00000E+00 0 particle: 2103 ud_1 ud_1bar 1 -1 0.77133 0.00000 0.00000 0.00000E+00 0 particle: 2112 n0 nbar0 0 0 0.93957 0.00000 0.00000 0.00000E+00 0 particle: 2114 Delta0 Deltabar0 0 0 1.23300 0.12000 0.14000 0.00000E+00 1 channel: 0.3310000 0 2212 -211 channel: 0.6630000 0 2112 111 channel: 0.0060000 0 2112 22 particle: 2203 uu_1 uu_1bar 4 -1 0.77133 0.00000 0.00000 0.00000E+00 0 particle: 2212 p+ pbar- 3 0 0.93827 0.00000 0.00000 0.00000E+00 0 particle: 2214 Delta+ Deltabar- 3 0 1.23200 0.12000 0.14000 0.00000E+00 1 channel: 0.6630000 0 2212 111 channel: 0.3310000 0 2112 211 channel: 0.0060000 0 2212 22 particle: 2224 Delta++ Deltabar-- 6 0 1.23100 0.12000 0.14000 0.00000E+00 1 channel: 1.0000000 0 2212 211 particle: 3101 sd_0 sd_0bar -2 -1 0.80473 0.00000 0.00000 0.00000E+00 0 particle: 3103 sd_1 sd_1bar -2 -1 0.92953 0.00000 0.00000 0.00000E+00 0 particle: 3112 Sigma- Sigmabar+ -3 0 1.19744 0.00000 0.00000 4.43400E+01 1 channel: 0.9990000 0 2112 -211 channel: 0.0010000 0 -12 11 2112 particle: 3114 Sigma*- Sigma*bar+ -3 0 1.38720 0.03940 0.04000 0.00000E+00 1 channel: 0.8800000 0 3122 -211 channel: 0.0600000 0 3212 -211 channel: 0.0600000 0 3112 111 particle: 3122 Lambda0 Lambdabar0 0 0 1.11568 0.00000 0.00000 7.88800E+01 1 channel: 0.6390000 0 2212 -211 channel: 0.3580000 0 2112 111 channel: 0.0020000 0 2112 22 channel: 0.0010000 0 -12 11 2212 particle: 3201 su_0 su_0bar 1 -1 0.80473 0.00000 0.00000 0.00000E+00 0 particle: 3203 su_1 su_1bar 1 -1 0.92953 0.00000 0.00000 0.00000E+00 0 particle: 3212 Sigma0 Sigmabar0 0 0 1.19255 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 3122 22 particle: 3214 Sigma*0 Sigma*bar0 0 0 1.38370 0.03600 0.03500 0.00000E+00 1 channel: 0.8800000 0 3122 111 channel: 0.0600000 0 3222 -211 channel: 0.0600000 0 3112 211 particle: 3222 Sigma+ Sigmabar- 3 0 1.18937 0.00000 0.00000 2.39600E+01 1 channel: 0.5160000 0 2212 111 channel: 0.4830000 0 2112 211 channel: 0.0010000 0 2212 22 particle: 3224 Sigma*+ Sigma*bar- 3 0 1.38280 0.03580 0.03500 0.00000E+00 1 channel: 0.8800000 0 3122 211 channel: 0.0600000 0 3222 111 channel: 0.0600000 0 3212 211 particle: 3303 ss_1 ss_1bar -2 -1 1.09361 0.00000 0.00000 0.00000E+00 0 particle: 3312 Xi- Xibar+ -3 0 1.32130 0.00000 0.00000 4.91000E+01 1 channel: 0.9988000 0 3122 -211 channel: 0.0001000 0 3112 22 channel: 0.0006000 0 -12 11 3122 channel: 0.0004000 0 -14 13 3122 channel: 0.0001000 0 -12 11 3212 particle: 3314 Xi*- Xi*bar+ -3 0 1.53500 0.00990 0.05000 0.00000E+00 1 channel: 0.6670000 0 3322 -211 channel: 0.3330000 0 3312 111 particle: 3322 Xi0 Xibar0 0 0 1.31490 0.00000 0.00000 8.71000E+01 1 channel: 0.9954000 0 3122 111 channel: 0.0011000 0 3122 22 channel: 0.0035000 0 3212 22 particle: 3324 Xi*0 Xi*bar0 0 0 1.53180 0.00910 0.05000 0.00000E+00 1 channel: 0.3330000 0 3322 111 channel: 0.6670000 0 3312 211 particle: 3334 Omega- Omegabar+ -3 0 1.67245 0.00000 0.00000 2.46000E+01 1 channel: 0.6760000 0 3122 -321 channel: 0.2340000 0 3322 -211 channel: 0.0850000 0 3312 111 channel: 0.0050000 0 -12 11 3322 particle: 4101 cd_0 cd_0bar 1 -1 1.96908 0.00000 0.00000 0.00000E+00 0 particle: 4103 cd_1 cd_1bar 1 -1 2.00808 0.00000 0.00000 0.00000E+00 0 particle: 4112 Sigma_c0 Sigma_cbar0 0 0 2.45210 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4122 -211 particle: 4114 Sigma*_c0 Sigma*_cbar0 0 0 2.50000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4122 -211 particle: 4122 Lambda_c+ Lambda_cbar- 3 0 2.28490 0.00000 0.00000 6.18000E-02 1 channel: 0.0180000 42 -11 12 3122 channel: 0.0050000 42 -11 12 3212 channel: 0.0050000 42 -11 12 3214 channel: 0.0030000 42 -11 12 2112 channel: 0.0020000 42 -11 12 2114 channel: 0.0060000 42 -11 12 2212 -211 channel: 0.0060000 42 -11 12 2112 111 channel: 0.0180000 42 -13 14 3122 channel: 0.0050000 42 -13 14 3212 channel: 0.0050000 42 -13 14 3214 channel: 0.0030000 42 -13 14 2112 channel: 0.0020000 42 -13 14 2114 channel: 0.0060000 42 -13 14 2212 -211 channel: 0.0060000 42 -13 14 2112 111 channel: 0.0066000 0 2224 -321 channel: 0.0250000 0 2224 -323 channel: 0.0160000 0 2212 -311 channel: 0.0088000 0 2212 -313 channel: 0.0050000 0 2214 -311 channel: 0.0050000 0 2214 -313 channel: 0.0058000 0 3122 211 channel: 0.0050000 0 3122 213 channel: 0.0055000 0 3212 211 channel: 0.0040000 0 3212 213 channel: 0.0040000 0 3214 211 channel: 0.0040000 0 3214 213 channel: 0.0040000 0 3222 111 channel: 0.0020000 0 3222 221 channel: 0.0020000 0 3222 331 channel: 0.0040000 0 3222 113 channel: 0.0040000 0 3222 223 channel: 0.0030000 0 3224 111 channel: 0.0020000 0 3224 221 channel: 0.0030000 0 3224 113 channel: 0.0030000 0 3224 223 channel: 0.0020000 0 3322 321 channel: 0.0020000 0 3322 323 channel: 0.0020000 0 3324 321 channel: 0.0010000 0 2224 -211 channel: 0.0010000 0 2224 -213 channel: 0.0020000 0 2212 111 channel: 0.0010000 0 2212 221 channel: 0.0010000 0 2212 331 channel: 0.0020000 0 2212 113 channel: 0.0020000 0 2212 223 channel: 0.0013000 0 2212 333 channel: 0.0018000 0 2212 10221 channel: 0.0010000 0 2214 111 channel: 0.0010000 0 2214 221 channel: 0.0010000 0 2214 331 channel: 0.0010000 0 2214 113 channel: 0.0010000 0 2214 223 channel: 0.0030000 0 2112 211 channel: 0.0030000 0 2112 213 channel: 0.0030000 0 2114 211 channel: 0.0030000 0 2114 213 channel: 0.0050000 0 3122 321 channel: 0.0050000 0 3122 323 channel: 0.0020000 0 3212 321 channel: 0.0020000 0 3212 323 channel: 0.0010000 0 3214 321 channel: 0.0010000 0 3214 323 channel: 0.0020000 0 3222 311 channel: 0.0020000 0 3222 313 channel: 0.0010000 0 3224 311 channel: 0.0010000 0 3224 313 channel: 0.2432000 13 2 -1 3 2101 channel: 0.0570000 13 2 -1 3 2103 channel: 0.0350000 13 2 -3 3 2101 channel: 0.0350000 13 2 -1 1 2101 channel: 0.1500000 13 3 2203 channel: 0.0750000 13 2 3201 channel: 0.0750000 13 2 3203 channel: 0.0300000 13 1 2203 channel: 0.0150000 13 2 2101 channel: 0.0150000 13 2 2103 particle: 4132 Xi_c0 Xi_cbar0 0 0 2.47030 0.00000 0.00000 2.90000E-02 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4201 cu_0 cu_0bar 4 -1 1.96908 0.00000 0.00000 0.00000E+00 0 particle: 4203 cu_1 cu_1bar 4 -1 2.00808 0.00000 0.00000 0.00000E+00 0 particle: 4212 Sigma_c+ Sigma_cbar- 3 0 2.45350 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4122 111 particle: 4214 Sigma*_c+ Sigma*_cbar- 3 0 2.50000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4122 111 particle: 4222 Sigma_c++ Sigma_cbar-- 6 0 2.45290 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4122 211 particle: 4224 Sigma*_c++ Sigma*_cbar-- 6 0 2.50000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4122 211 particle: 4232 Xi_c+ Xi_cbar- 3 0 2.46560 0.00000 0.00000 1.06000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4301 cs_0 cs_0bar 1 -1 2.15432 0.00000 0.00000 0.00000E+00 0 particle: 4303 cs_1 cs_1bar 1 -1 2.17967 0.00000 0.00000 0.00000E+00 0 particle: 4312 Xi'_c0 Xi'_cbar0 0 0 2.55000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4132 22 particle: 4314 Xi*_c0 Xi*_cbar0 0 0 2.63000 0.00000 0.00000 0.00000E+00 1 channel: 0.5000000 0 4132 111 channel: 0.5000000 0 4132 22 particle: 4322 Xi'_c+ Xi'_cbar- 3 0 2.55000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4232 22 particle: 4324 Xi*_c+ Xi*_cbar- 3 0 2.63000 0.00000 0.00000 0.00000E+00 1 channel: 0.5000000 0 4232 111 channel: 0.5000000 0 4232 22 particle: 4332 Omega_c0 Omega_cbar0 0 0 2.70400 0.00000 0.00000 1.90000E-02 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4334 Omega*_c0 Omega*_cbar0 0 0 2.80000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 4332 22 particle: 4403 cc_1 cc_1bar 4 -1 3.27531 0.00000 0.00000 0.00000E+00 0 particle: 4412 Xi_cc+ Xi_ccbar- 3 0 3.59798 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4414 Xi*_cc+ Xi*_ccbar- 3 0 3.65648 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4422 Xi_cc++ Xi_ccbar-- 6 0 3.59798 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4424 Xi*_cc++ Xi*_ccbar-- 6 0 3.65648 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4432 Omega_cc+ Omega_ccbar- 3 0 3.78663 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4434 Omega*_cc+ Omega*_ccbar- 3 0 3.82466 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 4444 Omega*_ccc++ Omega*_cccbar- 6 0 4.91594 0.00000 0.00000 1.00000E-01 1 channel: 0.0800000 42 -11 12 3 81 channel: 0.0800000 42 -13 14 3 81 channel: 0.7600000 11 2 -1 3 81 channel: 0.0800000 11 2 -3 3 81 particle: 5101 bd_0 bd_0bar -2 -1 5.38897 0.00000 0.00000 0.00000E+00 0 particle: 5103 bd_1 bd_1bar -2 -1 5.40145 0.00000 0.00000 0.00000E+00 0 particle: 5112 Sigma_b- Sigma_bbar+ -3 0 5.80000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5122 -211 particle: 5114 Sigma*_b- Sigma*_bbar+ -3 0 5.81000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5122 -211 particle: 5122 Lambda_b0 Lambda_bbar0 0 0 5.64100 0.00000 0.00000 3.42000E-01 1 channel: 0.1050000 42 -12 11 4122 channel: 0.1050000 42 -14 13 4122 channel: 0.0400000 42 -16 15 4122 channel: 0.0077000 0 4122 -211 channel: 0.0200000 0 4122 -213 channel: 0.0235000 0 4122 -20213 channel: 0.0285000 0 4122 -431 channel: 0.0435000 0 4122 -433 channel: 0.0011000 0 441 3122 channel: 0.0022000 0 443 3122 channel: 0.0044000 0 20443 3122 channel: 0.4291000 43 -2 1 4 2101 channel: 0.0800000 13 -2 4 1 2101 channel: 0.0700000 13 -4 3 4 2101 channel: 0.0200000 13 -4 4 3 2101 channel: 0.0150000 42 -2 1 2 2101 channel: 0.0050000 42 -4 3 2 2101 particle: 5132 Xi_b- Xi_bbar+ -3 0 5.84000 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5142 Xi_bc0 Xi_bcbar0 0 0 7.00575 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5201 bu_0 bu_0bar 1 -1 5.38897 0.00000 0.00000 0.00000E+00 0 particle: 5203 bu_1 bu_1bar 1 -1 5.40145 0.00000 0.00000 0.00000E+00 0 particle: 5212 Sigma_b0 Sigma_bbar0 0 0 5.80000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5122 111 particle: 5214 Sigma*_b0 Sigma*_bbar0 0 0 5.81000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5122 111 particle: 5222 Sigma_b+ Sigma_bbar- 3 0 5.80000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5122 211 particle: 5224 Sigma*_b+ Sigma*_bbar- 3 0 5.81000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5122 211 particle: 5232 Xi_b0 Xi_bbar0 0 0 5.84000 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5242 Xi_bc+ Xi_bcbar- 3 0 7.00575 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5301 bs_0 bs_0bar -2 -1 5.56725 0.00000 0.00000 0.00000E+00 0 particle: 5303 bs_1 bs_1bar -2 -1 5.57536 0.00000 0.00000 0.00000E+00 0 particle: 5312 Xi'_b- Xi'_bbar+ -3 0 5.96000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5132 22 particle: 5314 Xi*_b- Xi*_bbar+ -3 0 5.97000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5132 22 particle: 5322 Xi'_b0 Xi'_bbar0 0 0 5.96000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5232 22 particle: 5324 Xi*_b0 Xi*_bbar0 0 0 5.97000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5232 22 particle: 5332 Omega_b- Omega_bbar+ -3 0 6.12000 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5334 Omega*_b- Omega*_bbar+ -3 0 6.13000 0.00000 0.00000 0.00000E+00 1 channel: 1.0000000 0 5332 22 particle: 5342 Omega_bc0 Omega_bcbar0 0 0 7.19099 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5401 bc_0 bc_0bar 1 -1 6.67143 0.00000 0.00000 0.00000E+00 0 particle: 5403 bc_1 bc_1bar 1 -1 6.67397 0.00000 0.00000 0.00000E+00 0 particle: 5412 Xi'_bc0 Xi'_bcbar0 0 0 7.03724 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5414 Xi*_bc0 Xi*_bcbar0 0 0 7.04850 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5422 Xi'_bc+ Xi'_bcbar- 3 0 7.03724 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5424 Xi*_bc+ Xi*_bcbar- 3 0 7.04850 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5432 Omega'_bc0 Omega'_bcba 0 0 7.21101 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5434 Omega*_bc0 Omega*_bcbar0 0 0 7.21900 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5442 Omega_bcc+ Omega_bccbar- 3 0 8.30945 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5444 Omega*_bcc+ Omega*_bccbar- 3 0 8.31325 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5503 bb_1 bb_1bar -2 -1 10.07354 0.00000 0.00000 0.00000E+00 0 particle: 5512 Xi_bb- Xi_bbbar+ -3 0 10.42272 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5514 Xi*_bb- Xi*_bbbar+ -3 0 10.44144 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5522 Xi_bb0 Xi_bbbar0 0 0 10.42272 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5524 Xi*_bb0 Xi*_bbbar0 0 0 10.44144 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5532 Omega_bb- Omega_bbbar+ -3 0 10.60209 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5534 Omega*_bb- Omega*_bbbar+ -3 0 10.61426 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5542 Omega_bbc0 Omega_bbcbar0 0 0 11.70767 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5544 Omega*_bbc0 Omega*_bbcbar0 0 0 11.71147 0.00000 0.00000 3.87000E-01 1 channel: 0.1050000 42 -12 11 4 81 channel: 0.1050000 42 -14 13 4 81 channel: 0.0400000 42 -16 15 4 81 channel: 0.5000000 42 -2 1 4 81 channel: 0.0800000 42 -2 4 1 81 channel: 0.1400000 42 -4 3 4 81 channel: 0.0100000 42 -4 4 3 81 channel: 0.0150000 42 -2 1 2 81 channel: 0.0050000 42 -4 3 2 81 particle: 5554 Omega*_bbb- Omega*_bbbbar+ -3 0 15.11061 0.00000 0.00000 3.87000E-01 1 particle: 10111 a_00 void 0 0 0.98350 0.06000 0.05000 0.00000E+00 1 channel: 1.0000000 0 221 111 particle: 10113 b_10 void 0 0 1.23100 0.14200 0.25000 0.00000E+00 1 channel: 1.0000000 0 223 111 particle: 10211 a_0+ a_0- 3 0 0.98350 0.06000 0.05000 0.00000E+00 1 channel: 1.0000000 0 221 211 particle: 10213 b_1+ b_1- 3 0 1.23100 0.14200 0.25000 0.00000E+00 1 channel: 1.0000000 0 223 211 particle: 10221 f_0 void 0 0 1.00000 0.00000 0.00000 0.00000E+00 1 channel: 0.5200000 0 211 -211 channel: 0.2600000 0 111 111 channel: 0.1100000 0 321 -321 channel: 0.0550000 0 130 130 channel: 0.0550000 0 310 310 particle: 10223 h_1 void 0 0 1.17000 0.36000 0.20000 0.00000E+00 1 channel: 0.3330000 0 213 -211 channel: 0.3340000 0 113 111 channel: 0.3330000 0 -213 211 particle: 10311 K*_00 K*_0bar0 0 0 1.42900 0.28700 0.40000 0.00000E+00 1 channel: 0.6670000 0 321 -211 channel: 0.3330000 0 311 111 particle: 10313 K_10 K_1bar0 0 0 1.29000 0.09000 0.00500 0.00000E+00 1 channel: 0.2800000 0 321 -213 channel: 0.1400000 0 311 113 channel: 0.3130000 0 323 -211 channel: 0.1570000 0 313 111 channel: 0.1100000 0 311 223 particle: 10321 K*_0+ K*_0- 3 0 1.42900 0.28700 0.40000 0.00000E+00 1 channel: 0.6670000 0 311 211 channel: 0.3330000 0 321 111 particle: 10323 K_1+ K_1- 3 0 1.29000 0.09000 0.01000 0.00000E+00 1 channel: 0.2800000 0 311 213 channel: 0.1400000 0 321 113 channel: 0.3130000 0 313 211 channel: 0.1570000 0 323 111 channel: 0.1100000 0 321 223 particle: 10331 f'_0 void 0 0 1.40000 0.25000 0.35000 0.00000E+00 1 channel: 0.3600000 0 211 -211 channel: 0.1800000 0 111 111 channel: 0.0300000 0 321 -321 channel: 0.0150000 0 130 130 channel: 0.0150000 0 310 310 channel: 0.2000000 0 211 -211 211 -211 channel: 0.2000000 0 211 -211 111 111 particle: 10333 h'_1 void 0 0 1.40000 0.08000 0.00100 0.00000E+00 1 channel: 0.2500000 0 313 -311 channel: 0.2500000 0 -313 311 channel: 0.2500000 0 323 -321 channel: 0.2500000 0 -323 321 particle: 10411 D*_0+ D*_0- 3 0 2.27200 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 421 211 channel: 0.3330000 0 411 111 particle: 10413 D_1+ D_1- 3 0 2.42400 0.02000 0.08000 0.00000E+00 1 channel: 0.6670000 0 423 211 channel: 0.3330000 0 413 111 particle: 10421 D*_00 D*_0bar0 0 0 2.27200 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 411 -211 channel: 0.3330000 0 421 111 particle: 10423 D_10 D_1bar0 0 0 2.42400 0.02000 0.08000 0.00000E+00 1 channel: 0.6670000 0 413 -211 channel: 0.3330000 0 423 111 particle: 10431 D*_0s+ D*_0s- 3 0 2.50000 0.05000 0.10000 0.00000E+00 1 channel: 0.5000000 0 411 311 channel: 0.5000000 0 421 321 particle: 10433 D_1s+ D_1s- 3 0 2.53600 0.00000 0.00000 0.00000E+00 1 channel: 0.5000000 0 423 321 channel: 0.5000000 0 413 311 particle: 10441 chi_0c void 0 0 3.41510 0.01400 0.05000 0.00000E+00 1 channel: 0.0070000 0 443 22 channel: 0.9930000 12 82 -82 particle: 10443 h_1c void 0 0 3.46000 0.01000 0.02000 0.00000E+00 1 channel: 1.0000000 12 82 -82 particle: 10511 B*_00 B*_0bar0 0 0 5.68000 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 521 -211 channel: 0.3330000 0 511 111 particle: 10513 B_10 B_1bar0 0 0 5.73000 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 523 -211 channel: 0.3330000 0 513 111 particle: 10521 B*_0+ B*_0- 3 0 5.68000 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 511 211 channel: 0.3330000 0 521 111 particle: 10523 B_1+ B_1- 3 0 5.73000 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 513 211 channel: 0.3330000 0 523 111 particle: 10531 B*_0s0 B*_0sbar0 0 0 5.92000 0.05000 0.10000 0.00000E+00 1 channel: 0.5000000 0 521 -321 channel: 0.5000000 0 511 -311 particle: 10533 B_1s0 B_1sbar0 0 0 5.97000 0.05000 0.10000 0.00000E+00 1 channel: 0.5000000 0 523 -321 channel: 0.5000000 0 513 -311 particle: 10541 B*_0c+ B*_0c- 3 0 7.25000 0.05000 0.05000 0.00000E+00 1 channel: 0.5000000 0 511 411 channel: 0.5000000 0 521 421 particle: 10543 B_1c+ B_1c- 3 0 7.30000 0.05000 0.10000 0.00000E+00 1 channel: 0.5000000 0 513 411 channel: 0.5000000 0 523 421 particle: 10551 chi_0b void 0 0 9.85980 0.00000 0.00000 0.00000E+00 1 channel: 0.0200000 0 553 22 channel: 0.9800000 32 21 21 particle: 10553 h_1b void 0 0 9.87500 0.01000 0.02000 0.00000E+00 1 channel: 1.0000000 32 21 21 particle: 20113 a_10 void 0 0 1.23000 0.40000 0.30000 0.00000E+00 1 channel: 0.5000000 0 213 -211 channel: 0.5000000 0 -213 211 particle: 20213 a_1+ a_1- 3 0 1.23000 0.40000 0.30000 0.00000E+00 1 channel: 0.5000000 0 113 211 channel: 0.5000000 0 213 111 particle: 20223 f_1 void 0 0 1.28200 0.02500 0.05000 0.00000E+00 1 channel: 0.1460000 0 10211 -211 channel: 0.1460000 0 10111 111 channel: 0.1460000 0 -10211 211 channel: 0.0500000 0 221 211 -211 channel: 0.0500000 0 221 111 111 channel: 0.0500000 0 213 -211 111 channel: 0.1500000 0 113 211 -211 channel: 0.0500000 0 113 111 111 channel: 0.0500000 0 -213 211 111 channel: 0.0240000 0 321 -321 111 channel: 0.0240000 0 321 -311 -211 channel: 0.0240000 0 311 -311 111 channel: 0.0240000 0 311 -321 211 channel: 0.0660000 0 113 22 particle: 20313 K*_10 K*_1bar0 0 0 1.40200 0.17400 0.30000 0.00000E+00 1 channel: 0.6670000 0 323 -211 channel: 0.3330000 0 313 111 particle: 20323 K*_1+ K*_1- 3 0 1.40200 0.17400 0.30000 0.00000E+00 1 channel: 0.6670000 0 313 211 channel: 0.3330000 0 323 111 particle: 20333 f'_1 void 0 0 1.42700 0.05300 0.02000 0.00000E+00 1 channel: 0.2500000 0 313 -311 channel: 0.2500000 0 -313 311 channel: 0.2500000 0 323 -321 channel: 0.2500000 0 -323 321 particle: 20413 D*_1+ D*_1- 3 0 2.37200 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 423 211 channel: 0.3330000 0 413 111 particle: 20423 D*_10 D*_1bar0 0 0 2.37200 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 413 -211 channel: 0.3330000 0 423 111 particle: 20433 D*_1s+ D*_1s- 3 0 2.56000 0.05000 0.03000 0.00000E+00 1 channel: 0.5000000 0 423 321 channel: 0.5000000 0 413 311 particle: 20443 chi_1c void 0 0 3.51060 0.00090 0.00100 0.00000E+00 1 channel: 0.2730000 0 443 22 channel: 0.7270000 12 82 -82 particle: 20513 B*_10 B*_1bar0 0 0 5.78000 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 523 -211 channel: 0.3330000 0 513 111 particle: 20523 B*_1+ B*_1- 3 0 5.78000 0.05000 0.10000 0.00000E+00 1 channel: 0.6670000 0 513 211 channel: 0.3330000 0 523 111 particle: 20533 B*_1s0 B*_1sbar0 0 0 6.02000 0.05000 0.10000 0.00000E+00 1 channel: 0.5000000 0 523 -321 channel: 0.5000000 0 513 -311 particle: 20543 B*_1c+ B*_1c- 3 0 7.30000 0.05000 0.10000 0.00000E+00 1 channel: 0.5000000 0 513 411 channel: 0.5000000 0 523 421 particle: 20553 chi_1b void 0 0 9.89190 0.00000 0.00000 0.00000E+00 1 channel: 0.3500000 0 553 22 channel: 0.6500000 32 21 21 particle: 100443 psi' void 0 0 3.68600 0.00000 0.00000 0.00000E+00 1 channel: 0.0083000 0 11 -11 channel: 0.0083000 0 13 -13 channel: 0.1866000 12 82 -82 channel: 0.3240000 0 443 211 -211 channel: 0.1840000 0 443 111 111 channel: 0.0270000 0 443 221 channel: 0.0010000 0 443 111 channel: 0.0930000 0 10441 22 channel: 0.0870000 0 20443 22 channel: 0.0780000 0 445 22 channel: 0.0028000 0 441 22 particle: 100553 Upsilon' void 0 0 10.02330 0.00000 0.00000 0.00000E+00 1 channel: 0.0140000 0 11 -11 channel: 0.0140000 0 13 -13 channel: 0.0140000 0 15 -15 channel: 0.0080000 32 1 -1 channel: 0.0240000 32 2 -2 channel: 0.0080000 32 3 -3 channel: 0.0240000 32 4 -4 channel: 0.4250000 33 21 21 21 channel: 0.0200000 33 22 21 21 channel: 0.1850000 0 553 211 -211 channel: 0.0880000 0 553 111 111 channel: 0.0430000 0 10551 22 channel: 0.0670000 0 20553 22 channel: 0.0660000 0 555 22 particle: 1000001 ~d_L ~d_Lbar -1 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000002 ~u_L ~u_Lbar 2 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000003 ~s_L ~s_Lbar -1 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000004 ~c_L ~c_Lbar 2 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000005 ~b_1 ~b_1bar -1 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000006 ~t_1 ~t_1bar 2 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000011 ~e_L- ~e_L+ -3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000012 ~nu_eL ~nu_eLbar 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000013 ~mu_L- ~mu_L+ -3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000014 ~nu_muL ~nu_muLbar 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000015 ~tau_1- ~tau_1+ -3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000016 ~nu_tauL ~nu_tauLbar 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000021 ~g void 0 2 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000022 ~chi_10 void 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000023 ~chi_20 void 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000024 ~chi_1+ ~chi_1- 3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000025 ~chi_30 void 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000035 ~chi_40 void 0 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000037 ~chi_2+ ~chi_2- 3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 1000039 ~Gravitino void 0 0 500.00000 0.00000 0.00001 0.00000E+00 0 particle: 2000001 ~d_R ~d_Rbar -1 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000002 ~u_R ~u_Rbar 2 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000003 ~s_R ~s_Rbar -1 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000004 ~c_R ~c_Rbar 2 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000005 ~b_2 ~b_2bar -1 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000006 ~t_2 ~t_2bar 2 1 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000011 ~e_R- ~e_R+ -3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000012 ~nu_eR ~nu_eRbar 0 0 500.00000 0.00000 0.00001 0.00000E+00 0 particle: 2000013 ~mu_R- ~mu_R+ -3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000014 ~nu_muR ~nu_muRbar 0 0 500.00000 0.00000 0.00001 0.00000E+00 0 particle: 2000015 ~tau_2- ~tau_2+ -3 0 500.00000 1.00000 10.00000 0.00000E+00 1 particle: 2000016 ~nu_tauR ~nu_tauRbar 0 0 500.00000 0.00000 0.00001 0.00000E+00 0 particle: 3000111 pi_tc0 void 0 0 110.00000 0.02911 0.29108 0.00000E+00 1 particle: 3000211 pi_tc+ pi_tc- 3 0 110.00000 0.01741 0.17412 0.00000E+00 1 particle: 3000221 pi'_tc0 void 0 0 110.00000 0.04536 0.45362 0.00000E+00 1 particle: 3000331 eta_tc0 void 0 2 350.00000 0.09511 0.95114 0.00000E+00 1 particle: 3000113 rho_tc0 void 0 0 210.00000 0.86860 8.68604 0.00000E+00 1 particle: 3000213 rho_tc+ rho_tc- 3 0 210.00000 0.62395 6.23946 0.00000E+00 1 particle: 3000223 omega_tc void 0 0 210.00000 0.19192 1.91923 0.00000E+00 1 particle: 3100021 V8_tc void 0 2 500.00000 123.27638 450.00000 0.00000E+00 1 particle: 3100111 pi_22_1_tc void 0 0 125.00000 0.02296 0.22959 0.00000E+00 1 particle: 3200111 pi_22_8_tc void 0 2 250.00000 0.18886 1.88863 0.00000E+00 1 particle: 3100113 rho_11_tc void 0 2 400.00000 23.26819 232.68185 0.00000E+00 1 particle: 3200113 rho_12_tc void 0 2 350.00000 2.86306 28.63059 0.00000E+00 1 particle: 3300113 rho_21_tc void 0 2 350.00000 0.00000 0.00000 0.00000E+00 1 particle: 3400113 rho_22_tc void 0 2 300.00000 3.45903 34.59032 0.00000E+00 1 particle: 4000001 d* d*bar -1 1 400.00000 2.59359 25.93594 0.00000E+00 1 particle: 4000002 u* u*bar 2 1 400.00000 2.59687 25.96873 0.00000E+00 1 particle: 4000011 e*- e*bar+ -3 0 400.00000 0.42896 4.28961 0.00000E+00 1 particle: 4000012 nu*_e0 nu*_ebar0 0 0 400.00000 0.41912 4.19124 0.00000E+00 1 particle: 5000039 Graviton* void 0 0 1000.00000 0.14153 1.41528 0.00000E+00 1 particle: 9900012 nu_Re void 0 0 500.00000 0.00098 0.00977 0.00000E+00 1 particle: 9900014 nu_Rmu void 0 0 500.00000 0.00098 0.00976 0.00000E+00 1 particle: 9900016 nu_Rtau void 0 0 500.00000 0.00097 0.00973 0.00000E+00 1 particle: 9900023 Z_R0 void 0 0 1200.00000 26.72450 267.24501 0.00000E+00 1 particle: 9900024 W_R+ W_R- 3 0 750.00000 21.74916 217.49162 0.00000E+00 1 particle: 9900041 H_L++ H_L-- 6 0 200.00000 0.88159 8.81592 0.00000E+00 1 particle: 9900042 H_R++ H_R-- 6 0 200.00000 0.88001 8.80013 0.00000E+00 1 particle: 9900110 rho_diff0 void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 9900210 pi_diffr+ pi_diffr- 3 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 9900220 omega_di void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 9900330 phi_diff void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 9900440 J/psi_di void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 9902110 n_diffr0 n_diffrbar0 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 particle: 9902210 p_diffr+ p_diffrbar- 3 0 0.00000 0.00000 0.00000 0.00000E+00 0pythia8052/doc/PaxHeaders.94862/ParticleData.xml 0000644 0000000 0000000 00000000074 13466754007 016311 x ustar 00 30 atime=1619556553.936646086 30 ctime=1619556553.940646103 pythia8052/doc/ParticleData.xml 0000644 0002170 0000144 00000273216 13466754007 017260 0 ustar 00pythia users 0000000 0000000
modeME()
above, almost agrees with the Pythia6 ones, but several obsolete have
been removed, and a few moved for better consistency. Here is the list
of currently allowed codes:
code
- 10 number
of particles in the decay channelcode
- 201 d dbar -1 1 0.33000 0.00000 0.00000 0.00000E+00 0 2 u ubar 2 1 0.33000 0.00000 0.00000 0.00000E+00 0 3 s sbar -1 1 0.50000 0.00000 0.00000 0.00000E+00 0 4 c cbar 2 1 1.50000 0.00000 0.00000 0.00000E+00 0 5 b bbar -1 1 4.80000 0.00000 0.00000 0.00000E+00 0 6 t tbar 2 1 175.00000 1.39816 13.98156 0.00000E+00 1 7 b' b'bar -1 1 400.00000 0.00000 0.00000 0.00000E+00 1 8 t' t'bar 2 1 400.00000 0.00000 0.00000 0.00000E+00 1 11 e- e+ -3 0 0.00051 0.00000 0.00000 0.00000E+00 0 12 nu_e nu_ebar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 13 mu- mu+ -3 0 0.10566 0.00000 0.00000 6.58654E+05 0 1.0000000 42 -12 11 14 14 nu_mu nu_mubar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 15 tau- tau+ -3 0 1.77700 0.00000 0.00000 8.72000E-02 1 0.1783000 42 -12 11 16 0.1735000 42 -14 13 16 0.1131000 0 16 -211 0.2494000 0 16 -213 0.0030000 41 16 -211 111 0.0900000 41 16 -213 111 0.0027000 41 16 -211 111 111 0.0100000 41 16 -213 111 111 0.0014000 41 16 -211 111 111 111 0.0012000 41 16 -213 111 111 111 0.0002500 41 16 -211 310 0.0002500 41 16 -211 130 0.0071000 0 16 -321 0.0120000 0 16 -323 0.0004000 41 16 -321 111 0.0007500 41 16 -323 111 0.0000600 41 16 -323 111 111 0.0007800 41 16 -321 310 0.0007800 41 16 -321 130 0.0034000 41 16 -321 321 -211 0.0800000 41 16 -211 113 0.0110000 41 16 -211 211 -211 0.0191000 41 16 -211 223 0.0000600 41 16 -211 221 0.0050000 41 16 -213 113 0.0133000 41 16 -211 113 111 0.0067000 41 16 -213 211 -211 0.0005000 41 16 -211 211 -211 111 0.0035000 41 16 -213 223 0.0006000 41 16 -211 223 111 0.0015000 41 16 -213 221 0.0002100 41 16 -211 221 111 0.0002000 41 16 -213 113 111 0.0007500 41 16 -211 113 113 0.0001000 41 16 -211 221 221 0.0002000 41 16 -211 113 111 111 0.0011000 41 16 -213 113 111 111 0.0002000 41 16 -211 213 -213 0.0002000 41 16 -211 213 -211 111 0.0002000 41 16 -211 -213 211 111 0.0002200 41 16 -211 113 113 111 0.0004000 41 16 -323 111 111 0.0001000 41 16 -321 111 111 111 0.0020500 41 16 -211 310 111 0.0020500 41 16 -211 130 111 0.0006900 41 16 -321 310 111 0.0006900 41 16 -321 130 111 0.0002500 41 16 -211 310 310 0.0005100 41 16 -211 310 130 0.0002500 41 16 -211 130 130 16 nu_tau nu_taubar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 17 tau'- tau'+ -3 0 400.00000 0.00000 0.00000 0.00000E+00 1 18 nu'_tau nu'_taubar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 21 g void 0 2 0.00000 0.00000 0.00000 0.00000E+00 0 22 gamma void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 23 Z0 void 0 0 91.18800 2.47813 24.78129 0.00000E+00 1 24 W+ W- 3 0 80.45000 2.07115 20.71149 0.00000E+00 1 25 h0 void 0 0 115.00000 0.00367 0.03669 0.00000E+00 1 32 Z'0 void 0 0 500.00000 14.54029 145.40294 0.00000E+00 1 33 Z"0 void 0 0 900.00000 0.00000 0.00000 0.00000E+00 0 34 W'+ W'- 3 0 500.00000 16.66099 166.60993 0.00000E+00 1 35 H0 void 0 0 300.00000 8.38842 83.88423 0.00000E+00 1 36 A0 void 0 0 300.00000 3.37520 33.75195 0.00000E+00 1 37 H+ H- 3 0 300.00000 4.17669 41.76694 0.00000E+00 1 39 Graviton void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 41 R0 Rbar0 0 0 5000.00000 417.29147 4172.91467 0.00000E+00 1 42 LQ_ue LQ_uebar -1 1 200.00000 0.39162 3.91621 0.00000E+00 1 81 specflav void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 82 rndmflav rndmflavbar 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 88 junction void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 90 system void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 91 cluster void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 92 string void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 93 indep. void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 94 CMshower void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 95 SPHEaxis void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 96 THRUaxis void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 97 CLUSjet void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 98 CELLjet void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 99 table void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 110 reggeon void 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 111 pi0 void 0 0 0.13498 0.00000 0.00000 3.00000E-05 1 0.9880000 0 22 22 0.0120000 2 22 11 -11 113 rho0 void 0 0 0.76850 0.15100 0.40000 0.00000E+00 1 0.9987390 3 211 -211 0.0007900 0 111 22 0.0003800 0 221 22 0.0000460 0 13 -13 0.0000450 0 11 -11 115 a_20 void 0 0 1.31800 0.10700 0.25000 0.00000E+00 1 0.3472500 0 213 -211 0.3472500 0 -213 211 0.1440000 0 221 111 0.1040000 0 223 211 -211 0.0245000 0 321 -321 0.0122500 0 130 130 0.0122500 0 310 310 0.0028000 0 111 22 0.0057000 0 331 111 130 K_L0 void 0 0 0.49767 0.00000 0.00000 1.55000E+04 0 0.2112000 0 111 111 111 0.1256000 0 211 -211 111 0.1939000 42 -12 11 211 0.1939000 42 12 -11 -211 0.1359000 42 -14 13 211 0.1359000 42 14 -13 -211 0.0020000 0 211 -211 0.0010000 0 111 111 0.0006000 0 22 22 211 pi+ pi- 3 0 0.13957 0.00000 0.00000 7.80450E+03 0 0.9998770 0 -13 14 0.0001230 0 -11 12 213 rho+ rho- 3 0 0.76690 0.14900 0.40000 0.00000E+00 1 0.9995500 3 211 111 0.0004500 0 211 22 215 a_2+ a_2- 3 0 1.31800 0.10700 0.25000 0.00000E+00 1 0.3472500 0 213 111 0.3472500 0 113 211 0.1440000 0 221 211 0.1040000 0 223 211 111 0.0490000 0 321 -311 0.0028000 0 211 22 0.0057000 0 331 211 221 eta void 0 0 0.54745 0.00000 0.00000 0.00000E+00 1 0.3923000 0 22 22 0.3210000 0 111 111 111 0.2317000 0 211 -211 111 0.0478000 0 22 211 -211 0.0049000 2 22 11 -11 0.0013000 0 211 -211 11 -11 0.0003000 0 22 13 -13 0.0007000 0 111 22 22 223 omega void 0 0 0.78194 0.00843 0.10000 0.00000E+00 1 0.8900000 1 211 -211 111 0.0869300 0 22 111 0.0221000 3 211 -211 0.0008300 0 221 22 0.0000700 0 111 111 22 0.0000700 0 11 -11 225 f_2 void 0 0 1.27500 0.18500 0.17000 0.00000E+00 1 0.5640000 0 211 -211 0.2820000 0 111 111 0.0720000 0 211 -211 111 111 0.0280000 0 211 -211 211 -211 0.0230000 0 321 -321 0.0115000 0 130 130 0.0115000 0 310 310 0.0050000 0 221 221 0.0030000 0 111 111 111 111 310 K_S0 void 0 0 0.49767 0.00000 0.00000 2.67620E+01 1 0.6861000 0 211 -211 0.3139000 0 111 111 311 K0 Kbar0 0 0 0.49767 0.00000 0.00000 0.00000E+00 1 0.5000000 0 130 0.5000000 0 310 313 K*0 K*bar0 0 0 0.89610 0.05050 0.20000 0.00000E+00 1 0.6650000 3 321 -211 0.3330000 3 311 111 0.0020000 0 311 22 315 K*_20 K*_2bar0 0 0 1.43200 0.10900 0.12000 0.00000E+00 1 0.3330000 0 321 -211 0.1660000 0 311 111 0.1680000 0 323 -211 0.0840000 0 313 111 0.0870000 0 323 -211 111 0.0430000 0 313 211 -211 0.0590000 0 321 -213 0.0290000 0 311 113 0.0290000 0 311 223 0.0020000 0 311 221 321 K+ K- 3 0 0.49360 0.00000 0.00000 3.70900E+03 0 0.6352000 0 -13 14 0.2116000 0 211 111 0.0559000 0 211 211 -211 0.0173000 0 211 111 111 0.0482000 42 12 -11 111 0.0318000 42 14 -13 111 323 K*+ K*- 3 0 0.89160 0.04980 0.20000 0.00000E+00 1 0.6660000 3 311 211 0.3330000 3 321 111 0.0010000 0 321 22 325 K*_2+ K*_2- 3 0 1.42500 0.09800 0.12000 0.00000E+00 1 0.3320000 0 311 211 0.1660000 0 321 111 0.1680000 0 313 211 0.0840000 0 323 111 0.0860000 0 313 211 111 0.0430000 0 323 211 -211 0.0590000 0 311 213 0.0290000 0 321 113 0.0290000 0 321 223 0.0020000 0 321 221 0.0020000 0 321 22 331 eta' void 0 0 0.95777 0.00020 0.00200 0.00000E+00 1 0.4370000 0 211 -211 221 0.2080000 0 111 111 221 0.3020000 0 22 113 0.0302000 0 22 223 0.0212000 0 22 22 0.0016000 0 111 111 111 333 phi void 0 0 1.01940 0.00443 0.01500 0.00000E+00 1 0.4894700 3 321 -321 0.3400000 3 130 310 0.0430000 0 -213 211 0.0430000 0 113 111 0.0430000 0 213 -211 0.0270000 1 211 -211 111 0.0126000 0 22 221 0.0013000 0 111 22 0.0003000 0 11 -11 0.0002500 0 13 -13 0.0000800 0 211 -211 335 f'_2 void 0 0 1.52500 0.07600 0.20000 0.00000E+00 1 0.4440000 0 321 -321 0.2220000 0 130 130 0.2220000 0 310 310 0.1040000 0 221 221 0.0040000 0 211 -211 0.0040000 0 111 111 411 D+ D- 3 0 1.86930 0.00000 0.00000 3.17000E-01 1 0.0700000 42 -11 12 -311 0.0650000 42 -11 12 -313 0.0050000 42 -11 12 -311 111 0.0050000 42 -11 12 -321 211 0.0110000 42 -11 12 -313 111 0.0110000 42 -11 12 -323 211 0.0010000 42 -11 12 111 0.0010000 42 -11 12 221 0.0010000 42 -11 12 331 0.0010000 42 -11 12 113 0.0010000 42 -11 12 223 0.0700000 42 -13 14 -311 0.0650000 42 -13 14 -313 0.0050000 42 -13 14 -311 111 0.0050000 42 -13 14 -321 211 0.0110000 42 -13 14 -313 111 0.0110000 42 -13 14 -323 211