| 1 | = Omega production at TAPS = |
| 2 | |
| 3 | A sample job card: |
| 4 | |
| 5 | {{{ |
| 6 | !**************************************************** |
| 7 | ! Input specific for the reactions |
| 8 | !**************************************************** |
| 9 | |
| 10 | $target !Only relevant for eventtypes 0,1 and 3 (??????) |
| 11 | target_Z=20, |
| 12 | target_A=40, |
| 13 | fermimotion=.true. |
| 14 | densitySwitch_Static=1 !1=wood-saxon, 2=according NPA 554,554 (Oset) |
| 15 | $end |
| 16 | |
| 17 | $low_photo_induced |
| 18 | energy_gamma=1.1 |
| 19 | debugFlag=.false. |
| 20 | pascalTwoPi=.true. |
| 21 | delta_energy=0.01 |
| 22 | equalDistribution_twoPi=.true. |
| 23 | vecmes=.true. |
| 24 | resonances=.true. |
| 25 | singlePi=.true. |
| 26 | pi0eta=.true. |
| 27 | twopi=.true. |
| 28 | $end |
| 29 | |
| 30 | $lowPhotonAnalysis |
| 31 | outputEvents=.true. |
| 32 | KruscheOutput=.false. |
| 33 | fissumOutput=.false. |
| 34 | photonAnalyse=.true. |
| 35 | twoPiOutput=.false. |
| 36 | $end |
| 37 | |
| 38 | |
| 39 | !*************************************************** |
| 40 | ! General Input : |
| 41 | !*************************************************** |
| 42 | |
| 43 | $initDatabase |
| 44 | propagationSwitch=0 ! 0=all resonances with more than 1 star, 1=only delta&nucleon 2=only nucleons |
| 45 | ! used to initialize the resonances' database |
| 46 | $end |
| 47 | |
| 48 | |
| 49 | $initStability |
| 50 | stabilityFlag(101) = 4 |
| 51 | $end |
| 52 | |
| 53 | $pythia |
| 54 | MDCY(102,1)=1 ! KC code of pi0, not KF! |
| 55 | $end |
| 56 | |
| 57 | |
| 58 | $input |
| 59 | numEnsembles= 5 ! number of ensembles |
| 60 | eventtype = 3 ! 2=PionA, 3=photon A |
| 61 | numTimeSteps= 0 ! number of time steps |
| 62 | delta_T = 0.4 ! time step size |
| 63 | fullensemble=.false. |
| 64 | num_runs_SameEnergy=1 |
| 65 | num_Energies=1 |
| 66 | FinalCoulombCorrection=.false. |
| 67 | set_length_perturbative=.true. |
| 68 | !length_perturbative=5000 |
| 69 | length_perturbative=800 |
| 70 | PrintParticleVectors=.false. |
| 71 | DoPRLevel(1)=.false. |
| 72 | DoPRLevel(2)=.false. |
| 73 | DoPRLevel(3)=.false. |
| 74 | DoPRLevel(4)=.false. |
| 75 | $end |
| 76 | |
| 77 | $initRandom |
| 78 | SEED=45678 !Seed for the random number |
| 79 | $end |
| 80 | |
| 81 | $initDensity |
| 82 | densitySwitch=2 !1=dynamic density according to testparticle density, 2=analytic density prescription |
| 83 | splineExtraPolation=.true. !Switch for linear spline extrapolation for dynamically calculated density: Extrapolates density between |
| 84 | $end |
| 85 | |
| 86 | $initPauli |
| 87 | pauliSwitch=2 !1=dynamic, 2=analytic |
| 88 | $end |
| 89 | |
| 90 | $propagation |
| 91 | delta_P=0.01 ! Delta Momentum for derivatives |
| 92 | coulomb=.true. ! Whether to use coulomb in propagation |
| 93 | hadronic=.true. ! Whether to use hadronic potentials in propagation |
| 94 | DerivativeType=1 ! 1=first order Range-Kutta, 2=second order Range-Kutta |
| 95 | predictorCorrector=.false. ! Whether to use a predictor/corrector algorithm to do the propagation |
| 96 | $end |
| 97 | |
| 98 | !************************************************************ |
| 99 | ! Input for potentials |
| 100 | !************************************************************ |
| 101 | $Coulomb |
| 102 | CoulombFlag=.false. |
| 103 | $ |
| 104 | |
| 105 | $mesonPotential |
| 106 | pionPot_Switch=0 ! Switch for pionPotential |
| 107 | ! 1=Oset potential (NPA 554), which is valid up to 50Mev kinetic energy |
| 108 | ! 2=Kapusta suggestion for pion potential (rather unusual) |
| 109 | ! 3=Delta Hole potential, which is valid up to 130 MeV kinetic energy |
| 110 | ! 4=Smooth spline transition between switch 1 and 3. |
| 111 | ! else=no pion potential |
| 112 | $end |
| 113 | |
| 114 | $baryonPotential |
| 115 | EQS_Type=5, ! Switch for equation of state for nucleon resonances spin=1/2 |
| 116 | ! Parameters for nucleon potentials: |
| 117 | ! 1=soft mom-dep lambda = 2.130 |
| 118 | ! 2=hard mom-dep lambda = 2.126 |
| 119 | ! 3=soft non-mom-dep |
| 120 | ! 4=hard non-mom-dep |
| 121 | ! 5=medium mom-dep |
| 122 | DeltaPot=1, ! Switch for potential of spin=3/2 resonances |
| 123 | ! 1=nucleon (spin=1/2) potential times 3/5 [according to ericson/Weise book] |
| 124 | ! 2= 100 MeV *rho/rhoNull |
| 125 | $end |
| 126 | |
| 127 | |
| 128 | !********************************************************* |
| 129 | ! The collision term |
| 130 | !********************************************************** |
| 131 | |
| 132 | $hadronFormation |
| 133 | tauForma=0.8 ! formation proper time in restframe of hadron |
| 134 | $end |
| 135 | |
| 136 | $collisionTerm |
| 137 | energyCheck=0.01 ! accuracy of energy check in GeV |
| 138 | oneBodyProcesses=.true. |
| 139 | oneBodyAdditional=.true. |
| 140 | twoBodyProcesses=.true. |
| 141 | threeBodyProcesses=.true. |
| 142 | $end |
| 143 | |
| 144 | $insertion |
| 145 | minimumEnergy=0.005 |
| 146 | propagateNoPhoton=.false. |
| 147 | $end |
| 148 | |
| 149 | $master_1Body |
| 150 | $end |
| 151 | |
| 152 | $master_2Body |
| 153 | baryonBaryonScattering=.true. |
| 154 | baryonMesonScattering=.true. |
| 155 | mesonMesonScattering=.false. |
| 156 | $end |
| 157 | |
| 158 | $master_3Body |
| 159 | $end |
| 160 | !**************************************************************** |
| 161 | ! The widths of the particles |
| 162 | !*************************************************************** |
| 163 | |
| 164 | $width_Baryon |
| 165 | mediumSwitch=.false. ! Switch on/off in-medium width of all baryons at once -> The vacuum width are used. |
| 166 | mediumSwitch_coll=.false. ! Use consistent collisional broadening |
| 167 | mediumSwitch_Delta=.true. ! Switch on/off in-medium width of the delta. .false.=vacuum width |
| 168 | $end |
| 169 | |
| 170 | $width_Meson |
| 171 | mediumSwitch=.false. ! Switch on/off in-medium width of all mesons at once -> The vacuum width are used. |
| 172 | $end |
| 173 | |
| 174 | |
| 175 | !****************************************************************** |
| 176 | ! Temperature and thermodynamics |
| 177 | !**************************************************** |
| 178 | |
| 179 | $initThermoDynamics |
| 180 | temperatureSwitch=1 |
| 181 | ! 1=groundstate calculations (T=0,mu=E_F) |
| 182 | ! 2=the full procedure according to testparticle density (real particles only!) |
| 183 | $end |
| 184 | }}} |