Changes between Version 2 and Version 3 of jobCards/omegaTaps
- Timestamp:
- Feb 25, 2009, 2:41:20 PM (15 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
jobCards/omegaTaps
v2 v3 1 1 = Omega production at TAPS = 2 2 3 This job card is suited to analyze photon production in photon-proton or photon-nucleus reactions. It is divided into several parts 3 This job card is suited to analyze photon production in photon-proton or photon-nucleus reactions. 4 5 Each fortran name list given below represents input for one specific module in the code. Note that everything after an exclamation mark (!) is a comment in fortran name lists. 4 6 5 7 == Target definition == 8 The namelist ''target'' controls the target. Do not modify the ''fermimotion'' and ''densitySwitch_Static'' switches besides for testing. 6 9 {{{ 7 10 $target … … 20 23 }}} 21 24 == Initialization and final Analysis == 25 The next two namelists control the initialization step and the final analysis step: 22 26 {{{ 23 27 $low_photo_induced … … 43 47 }}} 44 48 == General Input == 45 {{{ 46 $initDatabase 47 propagationSwitch=0 ! 0=all resonances with more than 1 star, 1=only delta&nucleon 2=only nucleons 48 ! used to initialize the resonances' database 49 $end 50 51 52 $initStability 53 stabilityFlag(101) = 4 ! Let Pi^0 Decay 54 $end 55 56 $pythia 57 MDCY(102,1)=1 ! KC code of pi0, not KF! !Pi^0 unstable in Pythia 58 $end 59 60 49 These job cards control the general input for each run: 50 * 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. 51 * The variables ''num_runs_SameEnergy'' and ''num_Energies'' define how many subsequent runs are performed. 52 * Number of time steps: 53 * For proton targets we use ''numTimeSteps=0'' since there is no transport step necessary 54 * For nuclear targets, ''numTimeSteps'' should be chosen such that ''delta_T*numTimeSteps>40''. 55 * Please adjust ''path_to_input'' to your directory structure. 56 * The length of the perturbative vector should be adjusted to your needs to save memory. 57 {{{ 61 58 $input 62 59 numEnsembles= 1 ! number of ensembles 63 60 eventtype = 3 ! 3=photon A 64 61 numTimeSteps= 0 ! number of time steps 65 delta_T = 0. 4! time step size62 delta_T = 0.2 ! time step size 66 63 num_runs_SameEnergy = 1 ! Number of runs with the same energy 67 64 num_Energies = 1 ! Number of different energies … … 84 81 $end 85 82 86 83 $initDatabase 84 propagationSwitch=0 ! 0=all resonances with more than 1 star, 1=only delta&nucleon 2=only nucleons 85 ! used to initialize the resonances' database 86 $end 87 87 }}} 88 88 == Numerical details == 89 {{{ 90 91 89 Don't touch! 90 {{{ 92 91 $initDensity 93 92 densitySwitch=2 !1=dynamic density according to testparticle density, 2=analytic density prescription … … 110 109 }}} 111 110 == Input for potentials == 111 Influences the hadronic potentials. In the scenario below we use no potential for the mesons and our standard Skyrme-type potential for the baryons. 112 112 {{{ 113 113 … … 140 140 }}} 141 141 == The collision term == 142 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. 142 143 {{{ 143 144 … … 171 172 $end 172 173 174 175 176 $initStability 177 stabilityFlag(101) = 4 ! Let Pi^0 Decay 178 $end 179 180 $pythia 181 MDCY(102,1)=1 ! KC code of pi0, not KF! !Pi^0 unstable in Pythia 182 $end 183 184 173 185 }}} 174 186 == The widths of the particles == 175 {{{ 176 187 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 188 ''mediumSwitch_coll=.true.'' in the namelist width_Baryon and/or ''mediumSwitch=.true.'' in the namelist width_Meson. 189 {{{ 177 190 $width_Baryon 178 191 mediumSwitch =.true. ! Switch on/off in-medium width of all baryons at once -> The vacuum width are used.