./code/analysis/LowPhotonAnalysis.f90

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/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/twoPi_output
2.12. LowPhotonAnalysis/sigmaMomEEEE.dat
2.13. LowPhotonAnalysis/sigmaRadiusEEEE.dat
2.14. LowPhotonAnalysis/sigmadm_mom_EEEE_ZZZ.dat
2.15. LowPhotonAnalysis/twoPi_output_hist
2.16. 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=.true.

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:

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

LowPhotonAnalysis/analyze_Photon [ Subroutines ]

[ Top ] [ LowPhotonAnalysis ] [ Subroutines ]

NAME

subroutine analyze_Photon(teilchen,finalFlag)

INPUTS

type(particle), dimension(:,:), intent(in) :: teilchen ! particle vector
logical, intent (in) :: 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

type(particle), target, dimension(:,:) :: p ! particle vector

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 3gamma 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

type(particle), target, dimension(:,:) :: p ! 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(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(:,:), intent(in) :: p ! particle vector
logical, intent (in) :: 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, intent(in) :: 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(:,:), intent(in) :: p ! particle vector
logical, intent (in) :: 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(:,:), intent(in) :: teilchen ! particle vector
logical, intent (in) :: finalFlag ! .true. if it is the last call for one specific energy, therefore final output must be made.
real, intent(in) ,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(:,:), intent(in) :: p ! particle vector
logical, intent (in) :: 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, INTENT(IN),dimension(1:3) :: p
real, INTENT(IN),dimension(1:3) :: q
integer, INTENT(IN) :: chargep
integer, INTENT(IN) :: chargeq

OUTPUT NONE

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

TABLE OF CONTENTS

/LowPhotonAnalysis [ Modules ]

lowPhotonAnalysis/outputEvents [ Global module-variables ]

lowPhotonAnalysis/outputEvents_onlyFree [ Global module-variables ]

lowPhotonAnalysis/FissumOutput [ Global module-variables ]

lowPhotonAnalysis/KruscheOutput [ Global module-variables ]

lowPhotonAnalysis/KruscheAnalyse_cut [ Global module-variables ]

lowPhotonAnalysis/TwoPiOutput [ Global module-variables ]

lowPhotonAnalysis/photonAnalyse [ Global module-variables ]

lowPhotonAnalysis/pi0gamma_analysis [ Global module-variables ]

lowPhotonAnalysis/pi0gamma_momcut [ Global module-variables ]

lowPhotonAnalysis/pi0gamma_masscut [ Global module-variables ]

lowPhotonAnalysis/pi0gamma_mombin [ Global module-variables ]

lowPhotonAnalysis/pi0gamma_massres_sigma [ Global module-variables ]

lowPhotonAnalysis/pi0gamma_massres_tau [ Global module-variables ]

lowPhotonAnalysis/Ekin_pi0_cut [ Global module-variables ]

lowPhotonAnalysis/readinput [ Subroutines ]

lowPhotonAnalysis/lowPhotonAnalysis [ Namelists ]

LowPhotonAnalysis/analyze_Photon [ Subroutines ]

LowPhotonAnalysis/photonFS_Analyse [ Subroutines ]

LowPhotonAnalysis/pi0gamma_analyze [ Subroutines ]

LowPhotonAnalysis/pi0gamma_dsigma_dm.dat [ Output files ]

LowPhotonAnalysis/pi0gamma_dsigma_dp.dat [ Output files ]

LowPhotonAnalysis/pi0gamma_dsigma_dE.dat [ Output files ]

LowPhotonAnalysis/pi0gamma_dsigma_dEkin.dat [ Output files ]

LowPhotonAnalysis/OmegaExcitFunc_pi0gamma.dat [ Output files ]

LowPhotonAnalysis/kruscheAnalyse [ Subroutines ]

kruscheAnalyse/KruscheAnalyse.dat [ Output files ]

kruscheAnalyse/KruscheAnalyse_prelim.dat [ Output files ]

kruscheAnalyse/get_deltaE [ Subroutines ]

kruscheAnalyse/get_EPion [ Subroutines ]

LowPhotonAnalysis/fissumAnalyse [ Subroutines ]

LowPhotonAnalysis/twoPi_output [ Subroutines ]

LowPhotonAnalysis/sigmaMomEEEE.dat [ Output files ]

LowPhotonAnalysis/sigmaRadiusEEEE.dat [ Output files ]

LowPhotonAnalysis/sigmadm_mom_EEEE_ZZZ.dat [ Output files ]

LowPhotonAnalysis/twoPi_output_hist [ Subroutines ]

LowPhotonAnalysis/massInv_TwoPi [ Subroutines ]