GiBUU documentation

Table of Contents Sourcefiles Index Functions Global module-variables Modules Namelists Output files Program Subroutines Types

TABLE OF CONTENTS

• 1. lowPhotonAnalysis/lowPhotonAnalysis
• 1.1. lowPhotonAnalysis/outputEvents
• 1.2. lowPhotonAnalysis/outputEvents_onlyFree
• 1.3. lowPhotonAnalysis/FissumOutput
• 1.4. lowPhotonAnalysis/KruscheOutput
• 1.5. lowPhotonAnalysis/KruscheAnalyse_cut
• 1.6. lowPhotonAnalysis/TwoPiOutput
• 1.7. lowPhotonAnalysis/photonAnalyse
• 1.8. lowPhotonAnalysis/pi0gamma_analysis
• 1.9. lowPhotonAnalysis/pi0gamma_momcut
• 1.10. lowPhotonAnalysis/pi0gamma_masscut
• 1.11. lowPhotonAnalysis/pi0gamma_mombin
• 1.12. lowPhotonAnalysis/pi0gamma_massres_sigma
• 1.13. lowPhotonAnalysis/pi0gamma_massres_tau
• 1.14. lowPhotonAnalysis/Ekin_pi0_cut
• 1.15. lowPhotonAnalysis/eta_analysis
• 1.16. lowPhotonAnalysis/MissingMass_analysis
• 1.17. lowPhotonAnalysis/DoOutChannels
• 1.18. lowPhotonAnalysis/readinput
• 2. /LowPhotonAnalysis
• 2.1. LowPhotonAnalysis/analyze_Photon
• 2.2. LowPhotonAnalysis/photonFS_Analyse
• 2.3. LowPhotonAnalysis/pi0gamma_analyze
• 2.4. LowPhotonAnalysis/pi0gamma_dsigma_dm.dat
• 2.5. LowPhotonAnalysis/pi0gamma_dsigma_dp.dat
• 2.6. LowPhotonAnalysis/pi0gamma_dsigma_dE.dat
• 2.7. LowPhotonAnalysis/pi0gamma_dsigma_dEkin.dat
• 2.8. LowPhotonAnalysis/OmegaExcitFunc_pi0gamma.dat
• 2.9. LowPhotonAnalysis/kruscheAnalyse
• 2.9.1. kruscheAnalyse/KruscheAnalyse.dat
• 2.9.2. kruscheAnalyse/KruscheAnalyse_prelim.dat
• 2.9.3. kruscheAnalyse/get_deltaE
• 2.9.4. kruscheAnalyse/get_EPion
• 2.10. LowPhotonAnalysis/fissumAnalyse
• 2.11. LowPhotonAnalysis/etaAnalyse
• 2.12. LowPhotonAnalysis/twoPi_output
• 2.13. LowPhotonAnalysis/sigmaMomEEEE.dat
• 2.14. LowPhotonAnalysis/sigmaRadiusEEEE.dat
• 2.15. LowPhotonAnalysis/sigmadm_mom_EEEE_ZZZ.dat
• 2.16. LowPhotonAnalysis/twoPi_output_hist
• 2.17. LowPhotonAnalysis/massInv_TwoPi

---

/LowPhotonAnalysis [ Modules ]

[ Top ] [ Modules ]

NAME

module LowPhotonAnalysis

PURPOSE

Includes analysis routines for photon induced runs

---

lowPhotonAnalysis/outputEvents [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical,save  :: outputEvents=.false.

PURPOSE

If .true. then all events are printed to file.

---

lowPhotonAnalysis/outputEvents_onlyFree [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical,save  :: outputEvents_onlyFree=.false.

PURPOSE

If outputEvents=.true. then only particles which may leave the nucleus, i.e. may become "free", are printed to file.

---

lowPhotonAnalysis/FissumOutput [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical,save  ::  FissumOutput=.false.

PURPOSE

If .true. then we perform an analysis as in PRC 53,#3 pages 1278 ff. (1996) Produces dsigma/dOmega/dT_pi for pi^+

---

lowPhotonAnalysis/KruscheOutput [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical,save  ::  KruscheOutput=.false.

PURPOSE

If .true. then we perform an analysis as in EPJA22 347-351 (2004)

---

lowPhotonAnalysis/KruscheAnalyse_cut [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: KruscheAnalyse_cut=0. ! GeV

PURPOSE

Value of the cut for the deltaE cut in EPJA22 347-351 (2004).

---

lowPhotonAnalysis/TwoPiOutput [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical,save :: TwoPiOutput = .false.

PURPOSE

If .true. then we perform an analysis for 2pi production, including statistics for the mass of the pi-pi pair.

---

lowPhotonAnalysis/photonAnalyse [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical,save :: photonAnalyse=.false.

PURPOSE

Special analysis for final state photons

---

lowPhotonAnalysis/pi0gamma_analysis [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: pi0gamma_analysis = .false.

PURPOSE

Do analysis of pi0 gamma pairs (dsigma/dm), to reconstruct invariant mass spectrum of omega mesons.

---

lowPhotonAnalysis/pi0gamma_momcut [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: pi0gamma_momcut = 0.5

PURPOSE

Cut on the absolute omega three momentum in GeV, being applied to the pi0 gamma spectrum.

---

lowPhotonAnalysis/pi0gamma_masscut [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, dimension(1:2), save :: pi0gamma_masscut = (/0.,2./)

PURPOSE

Cuts on the pi0-gamma invariant mass in GeV, being applied to all pi0-gamma spectra (except the mass spectrum). First component is lower limit, second component is upper limit.

---

lowPhotonAnalysis/pi0gamma_mombin [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: pi0gamma_mombin = 0.050

PURPOSE

Bin size for the pi0 gamma momenentum spectrum in GeV.

---

lowPhotonAnalysis/pi0gamma_massres_sigma [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: pi0gamma_massres_sigma = 0.025

PURPOSE

Sigma parameter for the exp. resolution smearing (width of the Gauss or Novosibirsk function in GeV). See also pi0gamma_massres_tau.

---

lowPhotonAnalysis/pi0gamma_massres_tau [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: pi0gamma_massres_tau = -0.090

PURPOSE

Skewness parameter tau of the Novosibirsk function (for exp. resolution smearing). See also pi0gamma_massres_sigma.

---

lowPhotonAnalysis/Ekin_pi0_cut [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: Ekin_pi0_cut = 0.

PURPOSE

Cut on the kinetic energy of neutral pions in the pi0gamma_analysis. Only pions with kinetic energies larger than this cutoff are used for the analysis.

---

lowPhotonAnalysis/eta_analysis [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: eta_analysis = .false.

PURPOSE

Do analysis of eta mesons

---

lowPhotonAnalysis/MissingMass_analysis [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: MissingMass_analysis = .false.

PURPOSE

Do analysis of missing mass distributions for 1 nucleon and 2 nucleon final states

---

lowPhotonAnalysis/DoOutChannels [ Global module-variables ]

[ Top ] [ lowPhotonAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: DoOutChannels = .false.

PURPOSE

switch on/off: reporting of all final state channels

---

lowPhotonAnalysis/readinput [ Subroutines ]

[ Top ] [ lowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine readinput

INPUTS

NONE

OUTPUT

NONE

PURPOSE

This subroutine reads out the jobcard "photonAnalysis". Only called once to initialize the module.

---

lowPhotonAnalysis/lowPhotonAnalysis [ Namelists ]

[ Top ] [ Namelists ]

NAME

NAMELIST /lowPhotonAnalysis/

PURPOSE

Includes the switches:

• outputEvents
• outputEvents_onlyFree
• KruscheOutput
• KruscheAnalyse_cut
• FissumOutput
• photonAnalyse
• TwoPiOutput
• pi0gamma_analysis
• pi0gamma_momcut
• pi0gamma_masscut
• pi0gamma_mombin
• pi0gamma_massres_sigma
• pi0gamma_massres_tau
• Ekin_pi0_cut
• eta_analysis
• MissingMass_analysis
• DoOutChannels

---

LowPhotonAnalysis/analyze_Photon [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine analyze_Photon(Parts,finalFlag)

INPUTS

• type(particle), dimension(:,:) :: Parts -- particle vector
• logical, :: finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

NONE

PURPOSE

Main routine for the analysis. Calls all subroutines.

---

LowPhotonAnalysis/photonFS_Analyse [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine photonFS_Analyse(Parts)

INPUTS

• type(particle), target, dimension(:,:) :: Parts --- particle vector
• logical, finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

if "outputEvents" is true:

  * Prints all events to a file called 'FinalEvents.dat'
  * Events stemming from omega, rho or phi production are printed to
    seperate files 'FinalEvents.rho.dat', ...

if "photonAnalyse" is true:

  * Produces files called Photons_m2_vs_m3*.dat and Photons_dsigma_dm*.dat.
  * The first one contains the distribution of the invariant mass of all
    3 gamma combinations vs the inv. mass of all 2gamma combinations.
  * The second one contains the cross section dSigma/dm vs the mass of
    3 gammas.

PURPOSE

Groups all final state particles into events and writes them to file

---

LowPhotonAnalysis/pi0gamma_analyze [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine pi0gamma_analyze(Parts)

INPUTS

• type(particle), target, dimension(:,:) :: Parts -- particle vector

OUTPUT

• Produces an invariant mass spectrum of all pi0 gamma pairs in the final state and writes it to a file called 'pi0gamma_dsigma_dm.dat'

---

LowPhotonAnalysis/pi0gamma_dsigma_dm.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file pi0gamma_dsigma_dm.dat

PURPOSE

Pi0-gamma invariant-mass spectrum in microbarn/GeV/A.

---

LowPhotonAnalysis/pi0gamma_dsigma_dp.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file pi0gamma_dsigma_dp.dat

PURPOSE

Pi0-gamma momentum spectrum in microbarn/GeV/A.

---

LowPhotonAnalysis/pi0gamma_dsigma_dE.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file pi0gamma_dsigma_dE.dat

PURPOSE

Pi0-gamma energy spectrum in microbarn/GeV/A.

---

LowPhotonAnalysis/pi0gamma_dsigma_dEkin.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file pi0gamma_dsigma_dEkin.dat

PURPOSE

Pi0-gamma kinetic energy spectrum in microbarn/GeV/A.

---

LowPhotonAnalysis/OmegaExcitFunc_pi0gamma.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file OmegaExcitFunc_pi0gamma.dat

PURPOSE

Contains the omega excitation function, i.e. the energy-dependent inclusive omega photo-production cross section, as observed through the decay into pi0-gamma. In contrast to OmegaExcitFunc.dat, this does include FSI effects.

---

LowPhotonAnalysis/kruscheAnalyse [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine kruscheAnalyse(Parts,finalFlag)

PURPOSE

Evaluates the cross section for pi^0 production in the same way as Krusche in EPJA 22, 347 (2004)

INPUTS

• type(particle), dimension(:,:) :: Parts -- particle vector
• logical, finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

• See 'KruscheAnalyse.dat' for Xsections and delta(E) histograms are given in 'deltaE_E_gamma_'//realToChar4(energy_gamma)//'_GeV.dat'

---

kruscheAnalyse/KruscheAnalyse.dat [ Output files ]

[ Top ] [ kruscheAnalyse ] [ Output files ]

NAME

file KruscheAnalyse.dat

PURPOSE

Evaluate the cross section for pi^0 production in the same way as Krusche in EPJA 22, 347 (2004).

PURPOSE

Columns:

• 1 : E_gamma [Gev]
• 2 : sigma pi0 [mb]
• 3 : sigma pi0 after cut [mb]
• 4 : error sigma pi0 [mb]
• 5 : error sigma pi0 after cut [mb]

---

kruscheAnalyse/KruscheAnalyse_prelim.dat [ Output files ]

[ Top ] [ kruscheAnalyse ] [ Output files ]

NAME

file KruscheAnalyse_prelim.dat

PURPOSE

Same as KruscheAnalyse.dat, but output is written afer each run and not only the averaged result after each energy.

---

kruscheAnalyse/get_deltaE [ Subroutines ]

[ Top ] [ kruscheAnalyse ] [ Subroutines ]

NAME

real function get_deltaE(part)

PURPOSE

Evaluates for a given particle "delta E" according to eq. 4 in arXiv:nucl-ex/0406002v1 (June 2004).

INPUTS

• type(particle), intent(in) :: part

---

kruscheAnalyse/get_EPion [ Subroutines ]

[ Top ] [ kruscheAnalyse ] [ Subroutines ]

NAME

real function get_EPion(E_gamma)

PURPOSE

Evaluates the energy of a pion in the CM-frame of a gamma Nucleon collision.

INPUTS

• real:: E_gamma -- Energy of gamma in GeV in frame where nucleon rests

---

LowPhotonAnalysis/fissumAnalyse [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine fissumAnalyse(Parts,finalFlag)

PURPOSE

Evaluates the cross section for pi^+ production in the same binnings as Fissum PRC 53,3 pages 1278ff. (19996)

INPUTS

• type(particle), dimension(:,:) :: Parts -- particle vector
• logical :: finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

---

LowPhotonAnalysis/etaAnalyse [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine etaAnalyse(Parts,finalFlag)

PURPOSE

...

INPUTS

• type(particle), dimension(:,:) :: Parts -- particle vector
• logical :: finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

NOTES

the inclusive production cross sections are averaged over the different runs.

---

LowPhotonAnalysis/twoPi_output [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine twoPi_output(Parts, finalFlag,elab)

INPUTS

• type(particle), dimension(:,:) :: Parts -- particle vector
• logical :: finalFlag -- .true. if it is the last call for one specific & energy, therefore final output must be made.
• real, optional :: elab

OUTPUT

NONE

PURPOSE

Makes output for two-pion production processes. Writes to the following files:

• Total cross sections for pi pi production : "twoPi_Total.dat"
• dsigma/dm : "sigmadm*.dat"
• dsigma/dp "sigmaMom*.dat"
• dsigma/dm/dp "sigmadm_mom*.dat"

where * is the given energy elab or the photon energy if elab is not present.

Only for photon induced events also dsigm/dr is available in file 'sigmaRadius*.dat'

Here r is the radius of the production point of the measured pi pi pair.

NOTES

If you want to use this routine for non photon-induced events, then you must provide "elab"!!!

---

LowPhotonAnalysis/sigmaMomEEEE.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file sigmaMomEEEE.dat, EEEE= energy of the incoming photon in MeV

PURPOSE

The files are produced in the runs with eventtype=3=LowPhoton They show xsec (microbarn) for the __two-pion events__ versus pion momentum (GeV)

Columns:

• #1: momentum (absolute value of the pion 3-momentum)
• #2: PiMinus (xsec fo pi- )
• #3: PiNull (xsec fo pi0 )
• #4: PiPlus (xsec fo pi+ )

---

LowPhotonAnalysis/sigmaRadiusEEEE.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file sigmaRadiusEEEE.dat, EEEE= energy of the incoming photon in MeV

PURPOSE

The files are produced in the runs with eventtype=3=LowPhoton They show xsec (microbarn) for the __two-pion events__ versus radius of the nucleus for various final states

Columns:

• #1: radius[fm]
• #2: baryon density[fm^-3]
• #3: dsigma/dRadius (1:6) [10^-6 b/fm/A]
• #4: xsec for PiNull PiNull [microbarn]
• #5 xsec for PiNull PiPlus
• #6 xsec for piNull piMinus
• #7 xsec for piPlus piMinus
• #8 xsec for PiPlus PiPlus
• #9 xsec for piMinus piMinus

---

LowPhotonAnalysis/sigmadm_mom_EEEE_ZZZ.dat [ Output files ]

[ Top ] [ LowPhotonAnalysis ] [ Output files ]

NAME

file sigmadm_mom_EEEE_ZZZ.dat, EEEE= energy of the incoming photon in MeV ZZZ = channel (predefined two-photon final state)

PURPOSE

The files are produced in the runs with eventtype=3=LowPhoton They show xsec (microbarn) for the __two-pion events__ versus pion momentum for various final states

Columns:

• #1: mass [GeV] = invariant mass of the two outgoing pions
• #2: momentum [GeV]
• #3: dsigma/dm/dp [10^-6 b/GeV^2/A]

---

LowPhotonAnalysis/twoPi_output_hist [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine twoPi_output_hist(Parts,finalFlag)

PURPOSE

Makes dsigma/dm output for two-pion production processes.

INPUTS

• type(particle), dimension(:,:) :: Parts -- particle vector
• logical :: finalFlag -- .true. if it is the last call for one specific

energy, therefore final output must be made.

OUTPUT

Writes to the following files:

• twopi_dsigma_dm_"energy_gamma"_MeV.dat'

"energy_gamma" is the photon energy

---

LowPhotonAnalysis/massInv_TwoPi [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

function massInv_TwoPi(p,q,chargep,chargeq) Result(m)

INPUTS

• real, dimension(1:3) :: p
• real, dimension(1:3) :: q
• integer :: chargep
• integer :: chargeq

PURPOSE

Evaluates invariant mass of two pions with momentum p,q and charge chargep,chargeq

OUTPUT

real :: m -- m=SQRT((p+q)*(p+q)) where p and q are the 4-vectors

NOTES

real, parameter :: massPiZero=0.138 real, parameter :: massPiCharged=0.138

Generated from ./code/analysis/LowPhotonAnalysis.f90 with ROBODoc V4.99.42 on Sun Dec 22 2024 05:33:09