TABLE OF CONTENTS
/AntibaryonWidth [ Modules ]
NAME
module AntibaryonWidth
PURPOSE
Includes the routines administrating the simultaneous annihilation and calculation of annihilation and total widths of an antibaryon.
AntibaryonWidth/decideOnAnnihilation [ Subroutines ]
[ Top ] [ AntibaryonWidth ] [ Subroutines ]
NAME
subroutine decideOnAnnihilation(parts,time,flag,widthAverageOut,widthTotalAverageOut, baryonDensityAverageOut,sigmaAnniAverageOut,vRelAverageOut, srtS_vacuumAverageOut)
PURPOSE
Decide whether to annihilate an antiproton at a given time step.
INPUTS
OUTPUT
- logical:: flag -- if .true. --- do annihilation
- real, optional :: widthAverageOut --- annihilation width of an antiproton averaged over all parallel ensembles
- real, optional :: widthTotalAverageOut --- total width of an antiproton averaged over all parallel ensembles
- real, optional :: baryonDensityAverageOut --- average density of baryons
- real, optional :: sigmaAnniAverageOut --- average annihilation cross section
- real, optional :: vRelAverageOut --- average relative velocity
- real, optional :: srtS_vacuumAverageOut --- average vacuum c.m. energy
NOTES
Presently suitable for the real particles and for the parallel ensemble mode.
decideOnAnnihilation/gamma [ Subroutines ]
[ Top ] [ decideOnAnnihilation ] [ Subroutines ]
NAME
subroutine gamma(antiBaryon,time,width,widthTotal,baryonDensityOut,sigmaAnniOut,vRelOut,srtS_vacuumOut)
PURPOSE
Determine the annihilation width of an antibaryon.
INPUTS
OUTPUT
- real :: width ! width w/r to annihilation (c/fm)
- real :: widthTotal ! total width (c/fm)
- real, optional :: baryonDensityOut ! density of baryons (fm**-3)
- real, optional :: sigmaAnniOut ! annihilation cross section (mb)
- real, optional :: vRelOut ! relative velocity (c)
- real, optional :: srtS_vacuumOut ! c.m. energy (GeV) used in calculation of annihilation cross section
decideOnAnnihilation/ResetPosition [ Subroutines ]
[ Top ] [ decideOnAnnihilation ] [ Subroutines ]
NAME
subroutine ResetPosition
PURPOSE
- calculate maxID, i.e. the number of finalState-particles
- set position of finalState-particles: It loops over the FS particles vector and sets for the first FS-particle, which is baryon/meson as the first IS-particle, the spatial position as the first IS particle. Looping continues and the same is done for the next FS particle, which is of the same type as the secon IS particle. All other particles get as position the position of the interaction.
decideOnAnnihilation/gamma_exact [ Subroutines ]
[ Top ] [ decideOnAnnihilation ] [ Subroutines ]
NAME
subroutine gamma_exact(antiBaryon,time,width)
PURPOSE
Determine the annihilation width of an antibaryon.
INPUTS
- type(particle), intent(in) :: antiBaryon ! antibaryon particle
- real, intent(in) :: time ! actual time step
OUTPUT
- real, intent(out) :: width ! width w/r to annihilation (c/fm)
AntibaryonWidth/DoAnnihilation [ Subroutines ]
[ Top ] [ AntibaryonWidth ] [ Subroutines ]
NAME
subroutine DoAnnihilation(parts,time)
PURPOSE
Looks for an antibaryon and annihilates it with the closest baryon.
INPUTS
- type(particle), dimension(:,:) :: parts --- real particle vector
- real :: time --- actual time step
- real :: annihilationTime --- time when annihilation is enforced
OUTPUT
NOTES
Presently suitable for the real particles and for the parallel ensemble mode.