Omega production at TAPS

A sample job card:

!****************************************************
!       Input specific for the reactions
!****************************************************

$target      !Only relevant for eventtypes 0,1 and 3 (??????)
target_Z=20,
target_A=40,
fermimotion=.true.
densitySwitch_Static=1     !1=wood-saxon, 2=according NPA 554,554 (Oset)
$end

$low_photo_induced
energy_gamma=1.1
debugFlag=.false.
pascalTwoPi=.true.
delta_energy=0.01
equalDistribution_twoPi=.true.
vecmes=.true.
resonances=.true.
singlePi=.true.
pi0eta=.true.
twopi=.true.
$end

$lowPhotonAnalysis
outputEvents=.true.
KruscheOutput=.false.
fissumOutput=.false.
photonAnalyse=.true.
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
$end

$pythia
MDCY(102,1)=1 ! KC code of pi0, not KF!
$end


$input
numEnsembles= 5      ! number of ensembles
eventtype   = 3       ! 2=PionA, 3=photon A
numTimeSteps= 0       ! number of time steps
delta_T     = 0.4     ! time step size
fullensemble=.false.
num_runs_SameEnergy=1
num_Energies=1
FinalCoulombCorrection=.false.
set_length_perturbative=.true.
!length_perturbative=5000
length_perturbative=800
PrintParticleVectors=.false.
DoPRLevel(1)=.false.
DoPRLevel(2)=.false.
DoPRLevel(3)=.false.
DoPRLevel(4)=.false.
$end

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

$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=.false.  ! 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
propagateNoPhoton=.false.
$end

$master_1Body
$end

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

$master_3Body
$end
!****************************************************************
!   The widths of the particles
!***************************************************************

$width_Baryon
mediumSwitch=.false.       ! 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
!****************************************************

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