gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/pionXsection [ Modules ]

[ Top ] [ Modules ]

NAME

pionXsection

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.

pionXsection/debug [ Global module-variables ]

[ Top ] [ pionXsection ] [ Global module-variables ]

SOURCE

  logical, save :: debug=.false.

PURPOSE

If .true. debug infos are produced.

pionXsection/CMFrame [ Global module-variables ]

[ Top ] [ pionXsection ] [ 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.

pionXsection/dsigma_dOmegadE_switch [ Global module-variables ]

[ Top ] [ pionXsection ] [ Global module-variables ]

SOURCE

  logical, save :: dsigma_dOmegadE_switch=.false.

PURPOSE

If .true. then dsigma/dOmega and dSigma/dOmega/dE are evaluated.

pionXsection/twoPi_switch [ Global module-variables ]

[ Top ] [ pionXsection ] [ Global module-variables ]

SOURCE

  logical, save :: twoPi_switch=.false.

PURPOSE

If .true. then 2Pi output is evaluated.

pionXsection/readinput [ Subroutines ]

[ Top ] [ pionXsection ] [ 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.

pionXsection/pionAnalysis [ Namelists ]

[ Top ] [ pionXsection ] [ Namelists ]

NAME

NAMELIST /pionAnalysis/

PURPOSE

Includes the input switches:

pionXsection/pionXsectionAnalysis [ Subroutines ]

[ Top ] [ pionXsection ] [ Subroutines ]

NAME

subroutine pionXsectionAnalysis(pertParticles,finalFlag)

INPUTS

type(particle), intent(in),dimension(:,:) :: pertParticles ! Particles which shall be analyzed logical, intent(in) :: finalFlag ! if .true. than 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

AUTHOR

oliver buss

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 und Sigma_abs for Pions in file "pionInduced_Xsections.dat"

USES

To generate 2pi Xsection output we use the routine: 

  * use LowPhotonAnalysis, only : twoPi_Output

pionXsection/dSigmadOmegadE [ Subroutines ]

[ Top ] [ pionXsection ] [ Subroutines ]

NAME

subroutine dSigmadOmegadE(pertParticles,flag)

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*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.

pionXsection/dSigmadOmegadENN [ Subroutines ]

[ Top ] [ pionXsection ] [ 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 impact Parameter integration

Gives dSigma/dOmega/dE(pion) for Pions in file. Only those pions are considered which are produced in NN collisions.

USES

use statistics, only : getInfo Makes only sense if useStatistics=.true. in collisionTerm

pionXsection/dSigmadOmegadE_Resonance [ Subroutines ]

[ Top ] [ pionXsection ] [ Subroutines ]

NAME

subroutine dSigmadOmegadE_Resonance(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_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.

USES

use statistics, only : getInfo Makes only sense if useStatistics=.true. in collisionTerm

pionXsection/dSigmadOmegadE_pionNuk [ Subroutines ]

[ Top ] [ pionXsection ] [ Subroutines ]

NAME

subroutine dSigmadOmegadE_pionNuk(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_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.

USES

use statistics, only : getInfo Makes only sense if useStatistics=.true. in collisionTerm