TABLE OF CONTENTS
- 1. /gamma2Pi_Xsections
- 1.1. gamma2Pi_Xsections/gamma_2Pi_Xsections
- 1.2. gamma2Pi_Xsections/experimentalXsections
- 1.3. gamma2Pi_Xsections/init
- 1.4. gamma2Pi_Xsections/gamma2pi
- 1.5. gamma2Pi_Xsections/gP_ProtonPiPlusPiMinus
- 1.6. gamma2Pi_Xsections/gP_NeutronPiPlusPiNull
- 1.7. gamma2Pi_Xsections/gP_ProtonPiNullPiNul
- 1.8. gamma2Pi_Xsections/gN_NeutronPiPlusPiMinus
- 1.9. gamma2Pi_Xsections/gN_ProtonPiMinusPiNull
- 1.10. gamma2Pi_Xsections/gN_NeutronPiNullPiNull
- 1.11. gamma2Pi_Xsections/thres2pi
/gamma2Pi_Xsections [ Modules ]
NAME
module gamma2Pi_Xsections
PURPOSE
Includes the cross sections for gamma Nucleon -> Nucleon Pion Pion processes.
gamma2Pi_Xsections/gamma_2Pi_Xsections [ Namelists ]
[ Top ] [ gamma2Pi_Xsections ] [ Namelists ]
NAME
NAMELIST gamma_2Pi_Xsections
PURPOSE
Includes:
gamma2Pi_Xsections/experimentalXsections [ Global module-variables ]
[ Top ] [ gamma2Pi_Xsections ] [ Global module-variables ]
SOURCE
logical , save :: experimentalXsections=.true.
PURPOSE
- If .true. then the Xsections are taken from the experiment
- If .false. then the theoretical values are given
gamma2Pi_Xsections/init [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NAME
subroutine init
PURPOSE
Read out job card to initialize module parameters.
gamma2Pi_Xsections/gamma2pi [ Subroutines ]
[ Top ] [ gamma2Pi_Xsections ] [ Subroutines ]
NAME
subroutine gamma2pi(qnuk,srts,sig2pi,betaToLRF,mediumAtPosition,position)
PURPOSE
Calculation of 2 pion cross sections
INPUTS
OUTPUT
- sig2pi(0:3) Cross sections for gamma+nucleon->nucleon+2Pi production
- sig2pi(1) -> nucleon piMinus piPlus
- sig2pi(2) -> nucleon piPlus piNull or nucleon piMinus piNull
- sig2pi(3) -> nucleon piNull piNull
- sig2pi(0) -> Total Xsection into nucleon+2 Pions
NOTES
- Threshold improvement by Pascal Muehlich
- Units of cross sections : [10^-6 barn]
gamma2Pi_Xsections/gP_ProtonPiPlusPiMinus [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NOTES
Reads out data file and uses a spline with it.
OUTPUT
- real, optional :: minPlab,minSig ! plab and sigma of lowest lying data point in the data file
gamma2Pi_Xsections/gP_NeutronPiPlusPiNull [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NOTES
Reads out data file and uses a spline with it.
OUTPUT
- real, optional :: minPlab,minSig ! plab and sigma of lowest lying data point in the data file
gamma2Pi_Xsections/gP_ProtonPiNullPiNul [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NOTES
Reads out data file and uses a spline with it.
OUTPUT
- real, optional :: minPlab,minSig ! plab and sigma of lowest lying data point in the data file
gamma2Pi_Xsections/gN_NeutronPiPlusPiMinus [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NOTES
Reads out data file and uses a spline with it.
OUTPUT
- real, optional :: minPlab,minSig ! plab and sigma of lowest lying data point in the data file
gamma2Pi_Xsections/gN_ProtonPiMinusPiNull [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NOTES
Reads out data file and uses a spline with it.
OUTPUT
- real, optional :: minPlab,minSig ! plab and sigma of lowest lying data point in the data file
gamma2Pi_Xsections/gN_NeutronPiNullPiNull [ Functions ]
[ Top ] [ gamma2Pi_Xsections ] [ Functions ]
NOTES
Reads out data file and uses a spline with it.
OUTPUT
- real, optional :: minPlab,minSig ! plab and sigma of lowest lying data point in the data file
gamma2Pi_Xsections/thres2pi [ Subroutines ]
[ Top ] [ gamma2Pi_Xsections ] [ Subroutines ]
NAME
subroutine thres2pi(srts,sigma,qnuk,inThresholdRegion)
NOTES
- This subroutine provides a reasonable threshold behaviour for the two pion photoproduction cross sections. gamma nucleon -> nucleon+2pi
* Retrieves the lowest data point of all the parametrizations and fits a phase space curve to it.
AUTHOR
- Oliver Buss
INPUTS
OUTPUT
- real :: sigma(1:3) cross sections (1=double charged, 2= single charged, 3= pi^0 pi^0)
* logical, dimension(1:3), intent (out) :: inThresholdRegion
.true. if sqrt(s) is smaller than lowest data point in the channel (1=double charged, 2= single charged, 3= pi^0 pi^0)