gibuu is hosted by Hepforge, IPPP Durham
close Warning:
  • Error with navigation contributor "BrowserModule"
  • Failed to sync with repository "(default)": instance.__dict__ not accessible in restricted mode; repository information may be out of date. Look in the Trac log for more information including mitigation strategies.

Omega production at TAPS

The job card explained below is suited to analyze photon production in photon-proton or photon-nucleus reactions. Each fortran namelist given below represents input for one specific module in the code. Note: everything after an exclamation mark (!) is a comment in fortran namelists.

This jobcard can be found in the GiBUU repository as 003_metagExperiment.job

Target definition

The namelist target controls the target. Do not modify the fermimotion and densitySwitch_Static switches besides for testing.

! Proton:
! Nb_93
! Calcium_40
densitySwitch_Static=1 	   !1=wood-saxon, 2=according NPA 554,554 (Oset)

Initialization and final Analysis

The next two namelists control the initialization step and the final analysis step:

! Switch for specific initial channels
vecmes    =.true.
singlePi  =.true.
pi0eta    =.true.
twopi     =.true.

! Analysis flags
outputEvents          = .true.     ! Print events to file
outputEvents_onlyFree = .true.     ! Prints only "free" nucleons to file (.false.=print all nucleons)
photonAnalyse         = .true.     ! Generate analysis for final state photons
! Switch off unnessary analysis
KruscheOutput = .false.
fissumOutput  = .false.
twoPiOutput   = .false.

General Input

These namelists control the general input for each run:

  • The variable numEnsembles defines how many ensembles of test-particles are used to model a physical one, so it defines the granularity of our numerical realization. Use numEnsembles=1000...10000 for a real calculation.
  • The variables num_runs_SameEnergy and num_Energies define how many subsequent runs are performed.
  • Number of time steps:
    • For proton targets we use numTimeSteps=0 since there is no transport step necessary
    • For nuclear targets, numTimeSteps should be chosen such that delta_T*numTimeSteps>40.
  • Please adjust path_to_input to your directory structure.
  • The length of the perturbative vector should be adjusted to your needs to save memory.
numEnsembles= 1           ! number of ensembles
eventtype   = 3           ! 3=photon A
numTimeSteps= 0           ! number of time steps
delta_T     = 0.2         ! time step size
num_runs_SameEnergy = 1   ! Number of runs with the same energy
num_Energies        = 1   ! Number of different energies
! Length of particle vector. Must be adjusted to final state particle yield
! Proton
length_perturbative    =50   
! Calcium
!length_perturbative   =1000
! Niob
!length_perturbative    =3000
path_to_input='/home/hadron/oliver/buuinput_metag'     ! Path to input directory
fullensemble           =.false.
FinalCoulombCorrection =.false.

Numerical details

Don't touch!

densitySwitch=2            !1=dynamic density according to testparticle density, 2=analytic density prescription
splineExtraPolation=.true. !Switch for linear spline extrapolation for dynamically calculated density: Extrapolates density between 

pauliSwitch=2           !1=dynamic, 2=analytic

delta_P           =0.01    ! Delta Momentum for derivatives
coulomb           =.true.  ! Whether to use coulomb in propagation
hadronic          =.true.  ! Whether to use hadronic potentials in propagation
DerivativeType    =1       ! 1=first order Range-Kutta, 2=second order Range-Kutta
predictorCorrector=.true.  ! Whether to use a predictor/corrector algorithm to do the propagation

Input for potentials

Influences the hadronic potentials. In the scenario below we use no potential for the mesons and our standard Skyrme-type potential for the baryons.


pionPot_Switch=0  ! Switch for pionPotential
                  ! 1=Oset potential (NPA 554),  which is valid up to 50Mev kinetic energy
                  ! 2=Kapusta suggestion for pion potential (rather unusual)
                  ! 3=Delta Hole potential, which is valid up to 130 MeV kinetic energy
                  ! 4=Smooth spline transition between switch 1 and 3.
                  ! else=no pion potential

EQS_Type=5,  ! Switch for equation of state for nucleon resonances spin=1/2
             !     Parameters for nucleon potentials:
             !  1=soft mom-dep  lambda = 2.130
             !  2=hard mom-dep  lambda = 2.126
             !  3=soft  non-mom-dep
             !  4=hard  non-mom-dep
             !  5=medium  mom-dep
DeltaPot=1,  ! Switch for potential of spin=3/2 resonances
             ! 1=nucleon (spin=1/2) potential times  3/5   [according to ericson/Weise book]
             ! 2= 100 MeV *rho/rhoNull

The collision term

Here one can modify the collision term, e.g. by switching off three-body interactions. The scenario below is standard, so don't modify besides for testing. Note that the parameter minimumEnergy =0.005 in the namelist insertion removes all final state nucleons which have kinetic energies less than 5 MeV.

tauForma=0.8             ! formation proper time in restframe of hadron

energyCheck        =0.01          ! accuracy of energy check in GeV
oneBodyProcesses   =.true.
oneBodyAdditional  =.true.
twoBodyProcesses   =.true.
threeBodyProcesses =.true.

minimumEnergy      =0.005       ! Minimal kinetic energy for a proton
propagateNoPhoton  =.false.     ! Photons are propagated 

baryonBaryonScattering = .true.
baryonMesonScattering  = .true.
mesonMesonScattering   = .false.

stabilityFlag(101) = 4 ! Let Pi^0 Decay

MDCY(102,1)=1 ! KC code of pi0, not KF!  !Pi^0 unstable in Pythia

The widths of the particles

The scenario below corresponds to a broadening of the Delta, but to no broadening of any other particle. Such a broadening can be included by setting mediumSwitch_coll=.true. in the namelist width_Baryon and/or mediumSwitch=.true. in the namelist width_Meson.

mediumSwitch       =.true.   ! Switch on/off in-medium width of all baryons at once -> The vacuum width are used.
mediumSwitch_coll  =.false.  ! Use consistent collisional broadening
mediumSwitch_Delta =.true.   ! Switch on/off in-medium width of the delta. .false.=vacuum width

mediumSwitch=.false.         ! Switch on/off in-medium width of all mesons at once -> The vacuum width are used.

Temperature and thermodynamics

Don't touch this! Otherwise computation time blows up

! 1=groundstate calculations (T=0,mu=E_F)
! 2=the full procedure according to testparticle density (real particles only!)
Last modified 17 months ago Last modified on Jun 24, 2018, 11:14:00 AM