gibuu is hosted by Hepforge, IPPP Durham
GiBUU

Changes between Initial Version and Version 1 of jobCards/omegaTaps


Ignore:
Timestamp:
Feb 18, 2009, 3:34:33 PM (15 years ago)
Author:
oliver
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • jobCards/omegaTaps

    v1 v1  
     1= Omega production at TAPS =
     2
     3A sample job card:
     4
     5{{{
     6!****************************************************
     7!       Input specific for the reactions
     8!****************************************************
     9
     10$target      !Only relevant for eventtypes 0,1 and 3 (??????)
     11target_Z=20,
     12target_A=40,
     13fermimotion=.true.
     14densitySwitch_Static=1     !1=wood-saxon, 2=according NPA 554,554 (Oset)
     15$end
     16
     17$low_photo_induced
     18energy_gamma=1.1
     19debugFlag=.false.
     20pascalTwoPi=.true.
     21delta_energy=0.01
     22equalDistribution_twoPi=.true.
     23vecmes=.true.
     24resonances=.true.
     25singlePi=.true.
     26pi0eta=.true.
     27twopi=.true.
     28$end
     29
     30$lowPhotonAnalysis
     31outputEvents=.true.
     32KruscheOutput=.false.
     33fissumOutput=.false.
     34photonAnalyse=.true.
     35twoPiOutput=.false.
     36$end
     37
     38
     39!***************************************************
     40! General Input :
     41!***************************************************
     42
     43$initDatabase
     44propagationSwitch=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
     50stabilityFlag(101) = 4
     51$end
     52
     53$pythia
     54MDCY(102,1)=1 ! KC code of pi0, not KF!
     55$end
     56
     57
     58$input
     59numEnsembles= 5      ! number of ensembles
     60eventtype   = 3       ! 2=PionA, 3=photon A
     61numTimeSteps= 0       ! number of time steps
     62delta_T     = 0.4     ! time step size
     63fullensemble=.false.
     64num_runs_SameEnergy=1
     65num_Energies=1
     66FinalCoulombCorrection=.false.
     67set_length_perturbative=.true.
     68!length_perturbative=5000
     69length_perturbative=800
     70PrintParticleVectors=.false.
     71DoPRLevel(1)=.false.
     72DoPRLevel(2)=.false.
     73DoPRLevel(3)=.false.
     74DoPRLevel(4)=.false.
     75$end
     76
     77$initRandom
     78SEED=45678         !Seed for the random number
     79$end
     80
     81$initDensity
     82densitySwitch=2            !1=dynamic density according to testparticle density, 2=analytic density prescription
     83splineExtraPolation=.true. !Switch for linear spline extrapolation for dynamically calculated density: Extrapolates density between
     84$end
     85
     86$initPauli
     87pauliSwitch=2           !1=dynamic, 2=analytic
     88$end
     89
     90$propagation
     91delta_P=0.01                ! Delta Momentum for derivatives
     92coulomb=.true.             ! Whether to use coulomb in propagation
     93hadronic=.true.             ! Whether to use hadronic potentials in propagation
     94DerivativeType=1            ! 1=first order Range-Kutta, 2=second order Range-Kutta
     95predictorCorrector=.false.  ! Whether to use a predictor/corrector algorithm to do the propagation
     96$end
     97
     98!************************************************************
     99!       Input for potentials
     100!************************************************************
     101$Coulomb
     102CoulombFlag=.false.
     103$
     104
     105$mesonPotential
     106pionPot_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
     115EQS_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
     122DeltaPot=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
     133tauForma=0.8             ! formation proper time in restframe of hadron
     134$end
     135
     136$collisionTerm
     137energyCheck=0.01          ! accuracy of energy check in GeV
     138oneBodyProcesses=.true.
     139oneBodyAdditional=.true.
     140twoBodyProcesses=.true.
     141threeBodyProcesses=.true.
     142$end
     143
     144$insertion
     145minimumEnergy=0.005
     146propagateNoPhoton=.false.
     147$end
     148
     149$master_1Body
     150$end
     151
     152$master_2Body
     153baryonBaryonScattering=.true.
     154baryonMesonScattering=.true.
     155mesonMesonScattering=.false.
     156$end
     157
     158$master_3Body
     159$end
     160!****************************************************************
     161!   The widths of the particles
     162!***************************************************************
     163
     164$width_Baryon
     165mediumSwitch=.false.       ! Switch on/off in-medium width of all baryons at once -> The vacuum width are used.
     166mediumSwitch_coll=.false.  ! Use consistent collisional broadening
     167mediumSwitch_Delta=.true. ! Switch on/off in-medium width of the delta. .false.=vacuum width
     168$end
     169
     170$width_Meson
     171mediumSwitch=.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
     180temperatureSwitch=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}}}