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 2pi production becomes important. Therefore

there could be a pipi finalstate, which could lead to double counting.

pionXsection/pionAnalysis [ Namelists ]

[ Top ] [ pionXsection ] [ Namelists ]

NAME

NAMELIST /pionAnalysis/

PURPOSE

Includes the input switches:

pionXsection/debug [ Global module-variables ]

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

NAME

debug

SOURCE

  logical, save :: debug=.false.

PURPOSE

If .true. debug infos are produced, if .false. not..

pionXsection/CMFrame [ Global module-variables ]

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

NAME

CMFrame

SOURCE

  logical, save :: CMFrame=.false. ! Evaluate Xsection in CM frame of the incoming pion and a resting nucleon

PURPOSE

If .true. Xsection is evaluated in CM-Frame, if .false. in calculation frame.

pionXsection/dsigma_dOmegadE_switch [ Global module-variables ]

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

NAME

CMFrame

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 ]

NAME

CMFrame

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 out the jobcard "pionAnalysis" which may contain 

  * debug
  * CMFrame.

Only called once to initialize the module.

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