GiBUU documentation

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

TABLE OF CONTENTS

• 1. /yScalingAnalysis
• 1.1. yScalingAnalysis/debug
• 1.2. yScalingAnalysis/analyze
• 1.3. yScalingAnalysis/optionalOutput
• 1.4. yScalingAnalysis/variable
• 1.5. yScalingAnalysis/kFermi
• 1.6. yScalingAnalysis/E_shift
• 1.7. yScalingAnalysis/mA
• 1.8. yScalingAnalysis/mA_1
• 1.9. yScalingAnalysis/readinput
• 1.10. yScalingAnalysis/YScalingAnalysis
• 1.11. yScalingAnalysis/debug_prompt
• 1.12. yScalingAnalysis/calculateY
• 1.13. yScalingAnalysis/calculateYRes
• 1.14. yScalingAnalysis/mottXsection
• 1.15. yScalingAnalysis/cc1Xsection
• 1.16. yScalingAnalysis/scalingFunction
• 1.17. yScalingAnalysis/scaling_Analyze
• 1.18. yScalingAnalysis/scaling_info.dat
• 1.19. yScalingAnalysis/scaling_analysis.dat
• 1.20. yScalingAnalysis/scaling_analysis_fullRFG.dat
• 1.21. yScalingAnalysis/scaling_analysis_approxRFG.dat
• 1.22. yScalingAnalysis/scaling_analysis_fullPWIA.dat
• 1.23. yScalingAnalysis/single_nucleon.dat
• 1.24. yScalingAnalysis/scaling_delta.dat

---

/yScalingAnalysis [ Modules ]

[ Top ] [ Modules ]

NAME

module yScalingAnalysis

PURPOSE

• This module concludes the superscaling analysis of lepton induced processes
• For details on superscaling see Ivan Lappo-Danilevski diploma thesis at
• theorie.uni-giessen.de

---

yScalingAnalysis/debug [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: debug=.false.

PURPOSE

switches on and off debug infos: If .true. debug infos are produced, if .false. not

---

yScalingAnalysis/analyze [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: analyze=.false.

PURPOSE

Determines wether the y-scaling analysis is performed

---

yScalingAnalysis/optionalOutput [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: optionalOutput=.false.

PURPOSE

Determines wether in addition to the standard 'scaling_analysis.dat' other histograms will be generated. E.g. * 'single_nucleon.dat' - a table for comparing nucleon-knockout with fully inclusive

   x sections

* 'scaling_info.dat' - general parameters of the analysis, to be used for quick

   analyis

* 'scaling_delta.dat' - output to be used for analyis of scaling function in

   resonance excitation region

---

yScalingAnalysis/variable [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  integer, save :: variable=1

PURPOSE

determines which kind of scaling variable will be used (cf. Donnelly, Sick 1999) 1) RFG full variable Psi 2) RFG approximation Psi 3) PWIA full Upsilon (y/kf) 9) evaluation will be done for all variables, output written to seperate files

---

yScalingAnalysis/kFermi [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: kFermi=0.2251

PURPOSE

Nucleon Fermi momentum in nucleus. If none specified 0.2251 will be used, except if densitySwitch_static is set to 5, then fermiMomentum_input is used. The 0.225_1_ aims at preventing confusion whith delibaretely set differences between kFermi and fermiMomentum_input

---

yScalingAnalysis/E_shift [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: E_shift=0.020

PURPOSE

Energy correction to account for binding effects, otherwise neglected in RFG model

---

yScalingAnalysis/mA [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: mA=-1

PURPOSE

Tabulated mass of target nucleus

---

yScalingAnalysis/mA_1 [ Global module-variables ]

[ Top ] [ yScalingAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: mA_1=-1

PURPOSE

Tabulated mass of rest nucleus (target nucleus with knocked out nucleon)

---

yScalingAnalysis/readinput [ Subroutines ]

[ Top ] [ yScalingAnalysis ] [ Subroutines ]

NAME

subroutine readinput

INPUTS

NONE

OUTPUT

NONE

PURPOSE

This subroutine reads out the namelist "YScalingAnalysis". Only called once to initialize the module

---

yScalingAnalysis/YScalingAnalysis [ Namelists ]

[ Top ] [ yScalingAnalysis ] [ Namelists ]

NAME

NAMELIST /YScalingAnalysis/

PURPOSE

This namelist includes:

• debug
• analyze
• optionalOutput
• variable
• kFermi
• E_shift

---

yScalingAnalysis/debug_prompt [ Subroutines ]

[ Top ] [ yScalingAnalysis ] [ Subroutines ]

NAME

subroutine debug_prompt

INPUTS

NONE

OUTPUT

NONE

PURPOSE

This subroutine prints out basic debug informations for the y-scaling analysis (Not in use at the moment)

---

yScalingAnalysis/calculateY [ Functions ]

[ Top ] [ yScalingAnalysis ] [ Functions ]

NAME

function calculateY

INPUTS

• integer, intent(in) :: var -- type of scaling variable (cf. variable)
• real, intent(in) :: Li -- initial lepton energy
• real, intent(in) :: Lf -- final lepton energy
• real, intent(in) :: Qsquared -- four-momentum transfer

OUTPUT

• real :: calculateY -- resulting scaling variable

PURPOSE

This function calculates the scaling variable

---

yScalingAnalysis/calculateYRes [ Functions ]

[ Top ] [ yScalingAnalysis ] [ Functions ]

NAME

function calculateYRes

INPUTS

• real, intent(in) :: Li -- initial lepton energy
• real, intent(in) :: Lf -- final lepton energy
• real, intent(in) :: Qsquared -- four-momentum transfer
• real, intent(in) :: mass -- resonance mass

OUTPUT

• real :: calculateYRes -- resulting scaling variable

PURPOSE

This function calculates the scaling variable for the resonance region

---

yScalingAnalysis/mottXsection [ Functions ]

[ Top ] [ yScalingAnalysis ] [ Functions ]

NAME

function mottXsection

INPUTS

• real, intent(in) :: Ei -- incoming lepton energy
• real, intent(in) :: theta -- scattering angle
• integer, intent(in) :: Z -- target charge

OUTPUT

• real :: mottXsection

PURPOSE

This function calculates the Mott cross section for an ultrarelativistic spin 1/2 lepton. Result dimension is 1/Gev^2

---

yScalingAnalysis/cc1Xsection [ Subroutines ]

[ Top ] [ yScalingAnalysis ] [ Subroutines ]

NAME

function cc1Xsection

INPUTS

• real, intent(in) :: q -- 3-momentum transfer [GeV]
• real, intent(in) :: p -- momentum of struck nucleon [GeV]
• real, intent(in) :: Eps -- excitation energy [GeV]
• real, intent(in) :: mAi -- rest mass of target nucleus [GeV]
• real, intent(in) :: mA_1i -- rest mass of rest nucleus [GeV]

OUTPUT

• real, intent(out) :: wl -- longitudinal response
• real, intent(out) :: wt -- transverse response

PURPOSE

This subroutine calculates the off-shell single-nucleon responses according to de Forests cc1 perscription (1983, Nucl. Phys. A 392, 232-248). But following the conventions of Donnelly and Sick (in Phys. Rev. C 60, 065502).

---

yScalingAnalysis/scalingFunction [ Functions ]

[ Top ] [ yScalingAnalysis ] [ Functions ]

NAME

function scalingFunction

INPUTS

• integer :: var -- type of scaling variable (cf. variable)
• real :: Li -- initial lepton energy
• real :: Lf -- final lepton energy
• real :: Qsquared -- four-momentum transfer
• integer :: process_ID -- 1,2,3 for EM,CC,NC see initNeutrino.f90
• integer :: flavor_ID -- 1,2,3 for e,mu,tau

OUTPUT

• real :: scalingFunction

PURPOSE

This function calculates the scaling function f for different models

---

yScalingAnalysis/scaling_Analyze [ Subroutines ]

[ Top ] [ yScalingAnalysis ] [ Subroutines ]

NAME

subroutine scaling_Analyze(Particles,finalFlag,num_runs_sameEnergy,eventType)

INPUTS

• particles
• finalFlag
• num_runs_sameEnergy
• eventType

This subroutine produces output for lepton-nucleus y-scaling analysis

---

yScalingAnalysis/scaling_info.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file scaling_info.dat

PURPOSE

• Collection of information concerning the scaling analysis
• order of entries: kFermi, E_shift, N ,Z

NOTE

---

yScalingAnalysis/scaling_analysis.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file scaling_analysis.dat

PURPOSE

Evaluate scaling variable and scaling function (as described by Donnelly and Sick in Phys. Rev. C 60, 065502).

PURPOSE

Columns:

• 1 : scaling variable
• 2 : scaling function
• 3 : initial lepton energy [Gev]
• 4 : final lepton energy [Gev]
• 5 : momentum transfer Q^2 [Gev^2]
• 6 : angle of outgoing lepton [°]
• 7 : inclusive cross section [mb/GeV]
• 8 : mott cross section [mb/GeV]
• 9 : Longitudinal single-nucleon response GL
• 10: Transversal single-nucleon response GT

---

yScalingAnalysis/scaling_analysis_fullRFG.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file scaling_analysis_fullRFG.dat

PURPOSE

Specific analysis with full RFG scaling variable computation See : yScalingAnalysis/scaling_analysis.dat

---

yScalingAnalysis/scaling_analysis_approxRFG.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file scaling_analysis_approxRFG.dat

PURPOSE

Specific analysis with approximate RFG scaling variable computation See : yScalingAnalysis/scaling_analysis.dat

---

yScalingAnalysis/scaling_analysis_fullPWIA.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file scaling_analysis_fullPWIA.dat

PURPOSE

Specific analysis with full PWIA scaling variable computation See : yScalingAnalysis/scaling_analysis.dat

---

yScalingAnalysis/single_nucleon.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file single_nucleon.dat

PURPOSE

• A table for comparing nucleon-knockout with fully inclusive x sections
• order of entries: incoming lepton energy, outgoing lepton energy, sigma single-neutron, sigma single-proton, sigma single-nucleon sigma inclusive, ratio sig_sn/sig_inc
• Units are [mb/Gev]

NOTE should only differ significantly for timesteps > 0

---

yScalingAnalysis/scaling_delta.dat [ Output files ]

[ Top ] [ yScalingAnalysis ] [ Output files ]

NAME

file scaling_delta.dat

PURPOSE

Evaluate scaling variable for the resonance region and particulary the delta (as described by Maieron, Donnelly and Sick in Phys. Rev. C 65, 025502).

PURPOSE

Columns:

• 1 : initial lepton energy [Gev]
• 2 : final lepton energy [Gev]
• 3 : momentum transfer Q^2 [Gev^2]
• 4 : angle of outgoing lepton [°]
• 5 : inclusive cross section [mb/GeV]
• 6 : mott cross section [mb/GeV]
• 7 : scaling variable
• [8] : scaling function (to be filled in by 'delta_scaling.py'
• [9] : cross section as estimated by simple RFG model (to be filled in by 'delta_scaling.py')

Generated from ./code/analysis/yScalingAnalysis.f90 with ROBODoc V4.99.42 on Wed Feb 10 2016 21:00:10