gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


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


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

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

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')