gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/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/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:

LowPhotonAnalysis/analyze_Photon [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine analyze_Photon(teilchen,finalFlag)

INPUTS

  • type(particle), dimension(:,:) :: teilchen -- 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(p)

INPUTS

OUTPUT for "outputEvents" * 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', ... OUTPUT for "photonAnalyse" * 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(p)

INPUTS

OUTPUT

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(p,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(:,:) :: p -- particle vector
  • logical, finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

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

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(p,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(:,:) :: p -- particle vector
  • logical :: finalFlag -- .true. if it is the last call for one specific energy, therefore final output must be made.

OUTPUT

LowPhotonAnalysis/twoPi_output [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine twoPi_output(teilchen, finalFlag,elab)

INPUTS

  • type(particle), dimension(:,:) :: teilchen -- 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(p,finalFlag)

PURPOSE

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

INPUTS

  • type(particle), dimension(:,:) :: p -- 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:

"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