Omega production at TAPS

This job card is suited to analyze photon production in photon-proton or photon-nucleus reactions. It is divided into several parts

Target definition

$target      
! Proton:
target_Z=1,   
target_A=1,
! Nb_93
!target_Z=41,   
!target_A=93,
! Calcium_40
!target_Z=20,   
!target_A=40,
fermimotion=.true.
densitySwitch_Static=1 	   !1=wood-saxon, 2=according NPA 554,554 (Oset)
$end

Initialization and final Analysis

$low_photo_induced
energy_gamma=1.3
delta_energy=0.01
! Switch for specific initial channels
vecmes    =.true.
resonances=.true.
singlePi  =.true.
pi0eta    =.true.
twopi     =.true.
$end

$lowPhotonAnalysis
! Analysis flags
outputEvents   =.true.     ! Print events to file
photonAnalyse  =.true.     ! Generate analysis for final state photons
! Switch off unnessary analysis
KruscheOutput  =.false.
fissumOutput   =.false.
twoPiOutput    =.false.
$end

General Input

$initDatabase
propagationSwitch=0        ! 0=all resonances with more than 1 star, 1=only delta&nucleon 2=only nucleons
                           ! used to initialize the resonances' database
$end


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

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


$input
numEnsembles= 1           ! number of ensembles
eventtype   = 3           ! 3=photon A
numTimeSteps= 0           ! number of time steps
delta_T     = 0.4         ! time step size
num_runs_SameEnergy = 1   ! Number of runs with the same energy
num_Energies        = 1   ! Number of different energies
set_length_perturbative=.true.
! 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.
PrintParticleVectors=.false.
$end

$initRandom
SEED=45678         !Seed for the random number
$end

Numerical details

$initDensity
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 
$end

$initPauli
pauliSwitch=2           !1=dynamic, 2=analytic
$end

$propagation
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
$end

Input for potentials

$Coulomb
CoulombFlag=.false.
$

$mesonPotential
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
$end

$baryonPotential
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
$end

The collision term

$hadronFormation
tauForma=0.8             ! formation proper time in restframe of hadron
$end

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

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

$master_1Body
$end

$master_2Body
baryonBaryonScattering = .true.
baryonMesonScattering  = .true.
mesonMesonScattering   = .false.
$end

$master_3Body
$end

The widths of the particles

$width_Baryon
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
$end

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

Temperature and thermodynamics

Don't touch this! Otherwise computation time blows up

$initThermoDynamics
temperatureSwitch=1 
! 1=groundstate calculations (T=0,mu=E_F)
! 2=the full procedure according to testparticle density (real particles only!)
$end