TABLE OF CONTENTS
/pi0eta_photoproduction [ Modules ]
NAME
module pi0eta_photoproduction
PURPOSE
Includes a parametrization for gamma p -> p pi0 eta cross sections.
pi0eta_photoproduction/event_gamma_p_to_p_pi0_eta [ Subroutines ]
[ Top ] [ pi0eta_photoproduction ] [ Subroutines ]
NAME
subroutine event_gamma_p_to_p_pi0_eta(nuc,gamma,sigTot,srts,srtFree,betaCM,mediumAtPosition,fState,flag,perturbative)
PURPOSE
Generates a pi^0 eta nucleon event.
INPUTS
- type(particle),intent(in):: nuc,gamma ! incoming nucleon, photon
- real,intent(in):: sigTot ! total cross section
- real,intent(in):: srts ! sqrt(s)
- real,intent(in):: srtFree ! sqrt(s) in vacuum
- real,intent(in):: betaCM(3) ! beta for boost from CM to Lab frame
- logical , intent(in ) :: perturbative ! .true.= final state consists of perturbative particles
- type(medium),intent(in) :: mediumAtPosition
OUTPUT
- type(particle),intent(out):: fState(3) ! final state particles
- logical,intent(out):: flag ! true => success
pi0eta_photoproduction/sigma_gamma_p_to_p_pi0_eta [ Functions ]
[ Top ] [ pi0eta_photoproduction ] [ Functions ]
NAME
function sigma_gamma_p_to_p_pi0_eta(W) result(sigma)
PURPOSE
Returns a splined curve which goes through the data given in the file trim(path_to_input)//"/photo_pi0eta/gamma_p_TO_p_pi0_eta.dat"
INPUTS
real, intent(in) :: W -- Center of mass energy in vacuum
OUTPUT
real :: sigma -- Xsection in mu b
NOTES
If the point W lies outside the input data, then we set sigma=0 for W<smallest energy and assume linear interpolation for W>W_max.
pi0eta_photoproduction/readData [ Subroutines ]
[ Top ] [ pi0eta_photoproduction ] [ Subroutines ]
NAME
subroutine readData()
PURPOSE
Reads the data in the file trim(path_to_input)//"/photo_pi0eta/gamma_p_TO_p_pi0_eta.dat" to the field sigma_field. Max_index is set to the number of data points.