TABLE OF CONTENTS
- 1. /Coll_Pythia
- 1.1. Coll_Pythia/DoColl_Pythia
- 1.2. Coll_Pythia/Polym
- 1.3. Coll_Pythia/CalcMomentum
/Coll_Pythia [ Modules ]
NAME
module Coll_Pythia
FUNCTION
Implement a + b -> X processes/reactions done by PYTHIA
NOTES
There is a PreProcessor statement, which steers the weighting factor.
Coll_Pythia/DoColl_Pythia [ Subroutines ]
[ Top ] [ Coll_Pythia ] [ Subroutines ]
NAME
subroutine DoColl_Pythia (inPart, outPart, flagOK, sqrtS, pcm, beta, DoReweight, weight)
PURPOSE
Perform a collision of particles given in "inPart" with energy "sqrtS" and return outgoing particles in "outPart". "pcm" and "beta" are vectors used for Boost and Rotation of the event. if "flagOK" is false, no event happened, the output in "outPart" should be neglected! You can call Pythia in the "Reweighting mode" by setting "DoReweight" to true. Then "weight" has to be respected.
INPUTS
- type(particle),dimension(:) :: inPart -- incoming particles
- real :: sqrtS -- energy of ollision
- real, dimension(0:3) :: pcm -- boost-vector
- real, dimension(1:3) :: beta -- boost vector
- logical :: DoReweight -- Flag
OUTPUT
- type(particle),dimension(:) :: outPart ! outgoing particles
- real :: weight ! weight of event (default:1.0)
- logical :: flagOK ! event okay ?
NOTES
in order to understand the meaning of "pcm" and "beta": The (Pythia-)event is done in the restframe of the two given particles. Then a call to PYROBO according
phi = atan2(pcm(2),pcm(1))
theta = atan2(sqrt(pcm(1)**2+pcm(2)**2),pcm(3))
call PYROBO(1,N, theta,phi, beta(1),beta(2),beta(3))
is performed in order to transform the system into the desired (Lab-) system.
Coll_Pythia/Polym [ Functions ]
[ Top ] [ Coll_Pythia ] [ Functions ]
NAME
function Polym(facs,x)
PURPOSE
Calculate the value of a polynomial.
INPUTS
- real, dimension(:) :: facs -- coefficients
- real :: x -- value
OUTPUT
function value -- f(x)
NOTES
This is very elementary. Check code in order not to duplicate routines!
Coll_Pythia/CalcMomentum [ Subroutines ]
[ Top ] [ Coll_Pythia ] [ Subroutines ]
NAME
subroutine CalcMomentum(sqrts,mA,mB, EA,EB,pA)
PURPOSE
calculate the momentum in a two particle system (fundamental kinematics)
INPUTS
- sqrts,mA,mB, EA,EB -- cm energy, masses, energies
OUTPUT
NOTES
This is very elementary. Check code in order not to duplicate routines!