Changes between Version 1 and Version 2 of jobCards/omegaTaps

Timestamp:

Feb 25, 2009, 2:19:50 PM (16 years ago)

Author:

oliver

Comment:

--

jobCards/omegaTaps

@@ -1 +1 @@
 = Omega production at TAPS =
 
-A sample job card:
+This job card is suited to analyze photon production in photon-proton or photon-nucleus reactions. It is divided into several parts
 
+== Target definition ==
 {{{
-!****************************************************
-!       Input specific for the reactions
-!****************************************************
-
-$target      !Only relevant for eventtypes 0,1 and 3 (??????)
-target_Z=20,
-target_A=40,
+$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)
+densitySwitch_Static=1     !1=wood-saxon, 2=according NPA 554,554 (Oset)
 $end
-
+}}}
+== Initialization and final Analysis ==
+{{{
 $low_photo_induced
-energy_gamma=1.1
-debugFlag=.false.
-pascalTwoPi=.true.
+energy_gamma=1.3
 delta_energy=0.01
-equalDistribution_twoPi=.true.
-vecmes=.true.
+! Switch for specific initial channels
+vecmes    =.true.
 resonances=.true.
-singlePi=.true.
-pi0eta=.true.
-twopi=.true.
+singlePi  =.true.
+pi0eta    =.true.
+twopi     =.true.
 $end
 
 $lowPhotonAnalysis
-outputEvents=.true.
-KruscheOutput=.false.
-fissumOutput=.false.
-photonAnalyse=.true.
-twoPiOutput=.false.
+! 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 :
-!***************************************************
-
+}}}
+== General Input ==
+{{{
 $initDatabase
 propagationSwitch=0        ! 0=all resonances with more than 1 star, 1=only delta&nucleon 2=only nucleons
@@ -48 +51 @@
 
 $initStability
-stabilityFlag(101) = 4
+stabilityFlag(101) = 4 ! Let Pi^0 Decay
 $end
 
 $pythia
-MDCY(102,1)=1 ! KC code of pi0, not KF!
-$end
+MDCY(102,1)=1 ! KC code of pi0, not KF!  !Pi^0 unstable in Pythia
+$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.
+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_perturbative=5000
-length_perturbative=800
+! 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.
-DoPRLevel(1)=.false.
-DoPRLevel(2)=.false.
-DoPRLevel(3)=.false.
-DoPRLevel(4)=.false.
 $end
 
@@ -79 +84 @@
 $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
+splineExtraPolation=.true. !Switch for linear spline extrapolation for dynamically calculated density: Extrapolates density between 
 $end
 
@@ -89 +100 @@
 
 $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
+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
-!************************************************************
+
+}}}
+== Input for potentials ==
+{{{
+
 $Coulomb
 CoulombFlag=.false.
@@ -125 +138 @@
 $end
 
-
-!*********************************************************
-!          The collision term
-!**********************************************************
+}}}
+== The collision term ==
+{{{
 
 $hadronFormation
@@ -135 +147 @@
 
 $collisionTerm
-energyCheck=0.01          ! accuracy of energy check in GeV
-oneBodyProcesses=.true.
-oneBodyAdditional=.true.
-twoBodyProcesses=.true.
-threeBodyProcesses=.true.
+energyCheck        =0.01          ! accuracy of energy check in GeV
+oneBodyProcesses   =.true.
+oneBodyAdditional  =.true.
+twoBodyProcesses   =.true.
+threeBodyProcesses =.true.
 $end
 
 $insertion
-minimumEnergy=0.005
-propagateNoPhoton=.false.
+minimumEnergy      =0.005       ! Minimal kinetic energy for a proton
+propagateNoPhoton  =.false.     ! Photons are propagated 
 $end
 
@@ -151 +163 @@
 
 $master_2Body
-baryonBaryonScattering=.true.
-baryonMesonScattering=.true.
-mesonMesonScattering=.false.
+baryonBaryonScattering = .true.
+baryonMesonScattering  = .true.
+mesonMesonScattering   = .false.
 $end
 
 $master_3Body
 $end
-!****************************************************************
-!   The widths of the particles
-!***************************************************************
+
+}}}
+== 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
+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.
+mediumSwitch=.false.         ! Switch on/off in-medium width of all mesons at once -> The vacuum width are used.
 $end
-
-
-!******************************************************************
-! Temperature and thermodynamics
-!****************************************************
-
+}}}
+== Temperature and thermodynamics ==
+Don't touch this! Otherwise computation time blows up
+{{{
 $initThermoDynamics
-temperatureSwitch=1
+temperatureSwitch=1 
 ! 1=groundstate calculations (T=0,mu=E_F)
 ! 2=the full procedure according to testparticle density (real particles only!)