TABLE OF CONTENTS
- 1. /LoPionAnalysis
- 1.1. LoPionAnalysis/debug
- 1.2. LoPionAnalysis/CMFrame
- 1.3. LoPionAnalysis/dsigma_dOmegadE_switch
- 1.4. LoPionAnalysis/twoPi_switch
- 1.5. LoPionAnalysis/readinput
- 1.6. LoPionAnalysis/pionAnalysis
- 1.7. LoPionAnalysis/DoLoPionAnalysis
- 1.8. LoPionAnalysis/dSigmadOmegadE
- 1.9. LoPionAnalysis/dSigmadOmegadENN
- 1.10. LoPionAnalysis/dSigmadOmegadE_Resonance
- 1.11. LoPionAnalysis/dSigmadOmegadE_pionNuk
/LoPionAnalysis [ Modules ]
NAME
LoPionAnalysis
PURPOSE
Module with all routines necessary to make output for pion induced processes. Valid up to 180 MeV incoming pion energy (above that 2pi production becomes important).Therefore there could be a pipi finalstate, which could lead to double counting.
LoPionAnalysis/debug [ Global module-variables ]
[ Top ] [ LoPionAnalysis ] [ Global module-variables ]
SOURCE
logical, parameter :: debug=.false.
PURPOSE
If .true. debug infos are produced.
LoPionAnalysis/CMFrame [ Global module-variables ]
[ Top ] [ LoPionAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: CMFrame = .false.
PURPOSE
If .true. Xsection is evaluated in CM-Frame of the incoming pion and a resting nucleon, else in calculation frame.
LoPionAnalysis/dsigma_dOmegadE_switch [ Global module-variables ]
[ Top ] [ LoPionAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: dsigma_dOmegadE_switch=.false.
PURPOSE
If .true. then dsigma/dOmega and dSigma/dOmega/dE are evaluated.
LoPionAnalysis/twoPi_switch [ Global module-variables ]
[ Top ] [ LoPionAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: twoPi_switch=.false.
PURPOSE
If .true. then 2Pi output is evaluated.
LoPionAnalysis/readinput [ Subroutines ]
[ Top ] [ LoPionAnalysis ] [ Subroutines ]
NAME
subroutine readinput
INPUTS
NONE
OUTPUT
NONE
PURPOSE
This subroutine reads the namelist "pionAnalysis" out of the jobcard. Only called once to initialize the module.
LoPionAnalysis/pionAnalysis [ Namelists ]
[ Top ] [ LoPionAnalysis ] [ Namelists ]
NAME
NAMELIST /pionAnalysis/
PURPOSE
Includes the input switches:
LoPionAnalysis/DoLoPionAnalysis [ Subroutines ]
[ Top ] [ LoPionAnalysis ] [ Subroutines ]
NAME
subroutine DoLoPionAnalysis(pertParticles,finalFlag)
INPUTS
type(particle), intent(in),dimension(:,:) :: pertParticles --- Particles which shall be analyzed logical, intent(in) :: finalFlag --- if .true. then the final output for a series of calls will be done
OUTPUT
- File './pionInduced_xSections.dat' with absorption, reaction and quasielastic cross sections
- File 'pionInduced_dTheta'//RealToChar(getEkin()*1000.) //'MeV.dat' with quasielastic cross section : dsigma/dtheta in mB/radian
NOTES
This routine is based on "pionout" in inpion.f by M.Effenberger
PURPOSE
- This subroutine produces output for pion-nucleus scattering, including impact parameter integration.
- Gives dSigma/dOmega for Pions in file "pionInduced_dTheta... .dat"
- Gives Sigma_reac and Sigma_abs for Pions in "pionInduced_Xsections.dat"
USES
To generate 2pi Xsection output we use the routine:
- use LowPhotonAnalysis, only : twoPi_Output
LoPionAnalysis/dSigmadOmegadE [ Subroutines ]
[ Top ] [ LoPionAnalysis ] [ Subroutines ]
NAME
subroutine dSigmadOmegadE(pertParticles,flag)
INPUTS
type(particle), dimension(:,:) :: pertParticles ! The particle vector logical :: flag ! Whether it's the final run
OUTPUT
The files:
- pionInduced_dOmegadE_piMinus*MeV.dat
- pionInduced_dOmegadE_piNull*MeV.dat
- pionInduced_dOmegadE_piPlus*MeV.dat
where * is the kinetic energy of the incoming pions.
PURPOSE
This subroutine produces output for pion-nucleus scattering, including impact parameter integration. Gives dSigma/dOmega/dE(pion) for Pions in file.
LoPionAnalysis/dSigmadOmegadENN [ Subroutines ]
[ Top ] [ LoPionAnalysis ] [ Subroutines ]
NAME
subroutine dSigmadOmegadE_NN(pertParticles,finalFlag)
INPUTS
type(particle), intent(in),dimension(:,:) :: pertParticles ! The particle vector logical,intent(in) :: flag ! Whether it's the final run
OUTPUT
The files:
- pionInduced_dOmegadE_piMinus_NN*MeV.dat
- pionInduced_dOmegadE_piNull_NN*MeV.dat
- pionInduced_dOmegadE_piPlus_NN*MeV.dat
where * is the kinetic energy of the incoming pions.
PURPOSE
This subroutine produces output for pion-nucleus scattering, including impactp arameter integration
Gives dSigma/dOmega/dE(pion) for Pions in file. Only those pions are considered which are produced in NN collisions.
NOTES
Makes only sense if useStatistics=.true. in collisionTerm
LoPionAnalysis/dSigmadOmegadE_Resonance [ Subroutines ]
[ Top ] [ LoPionAnalysis ] [ Subroutines ]
NAME
subroutine dSigmadOmegadE_Resonance(pertParticles,finalFlag)
INPUTS
- type(particle), intent(in),dimension(:,:) :: pertParticles --- The particle vector
- logical,intent(in) :: finalFlag __ Whether it's the final run
OUTPUT
The files:
- pionInduced_dOmegadE_piMinus_R*MeV.dat
- pionInduced_dOmegadE_piNull_R*MeV.dat
- pionInduced_dOmegadE_piPlus_R*MeV.dat
where * is the kinetic energy of the incoming pions.
PURPOSE
This subroutine produces output for pion-nucleus scattering, including impact parameter integration
Gives dSigma/dOmega/dE(pion) for Pions in file. Only those pions are considered which are produced in Resonance decays.
NOTES
Makes only sense if useStatistics=.true. in collisionTerm
LoPionAnalysis/dSigmadOmegadE_pionNuk [ Subroutines ]
[ Top ] [ LoPionAnalysis ] [ Subroutines ]
NAME
subroutine dSigmadOmegadE_pionNuk(pertParticles,finalFlag)
INPUTS
- type(particle), intent(in),dimension(:,:) :: pertParticles --- The particle vector
- logical,intent(in) :: finalFlag --- Whether it's the final run
OUTPUT
The files:
- pionInduced_dOmegadE_piMinus_piN*MeV.dat
- pionInduced_dOmegadE_piNull_piN*MeV.dat
- pionInduced_dOmegadE_piPlus_piN*MeV.dat
where * is the kinetic energy of the incoming pions.
PURPOSE
This subroutine produces output for pion-nucleus scattering, including impact parameter integration
Gives dSigma/dOmega/dE(pion) for Pions in file. Only those pions are considered which are produced in pi N collisions.
NOTES
Makes only sense if useStatistics=.true. in collisionTerm