Example jobcard for High energy photon induced
!************************************************************ -*- fortran -*-
!* A JobCard for BUU:
!* ==================
!*
!* This JobCard tries to be a template for all
!* "HiEnergy photon+Nucleus"
!* collisions.
!* Here: on NUCLEON/DEUTERIUM
!* -- fermi motion : YES
!* -- shadowing : no
!* -- propagation : no
!* -- decays : YES
!*
!* This is just a special case of jobHiLepton_N.
!* You just give a small, but nonvanishing value for Q2.
!* Then nu is the photon energy and the actual value of
!* epsilon influences the total cross section in the percent
!* region (but may have some influence, whether the code
!* may 'reconstruct' the kinematics of a scattered electron.
!************************************************************
!************************************************************
! General Input :
!************************************************************
$input
eventtype = 14 ! HiLepton
numEnsembles= 200 ! number of ensembles
numTimeSteps= 0 ! number of time steps
set_length_perturbative = .TRUE.
length_perturbative = 50
! num_runs_SameEnergy=10000 ! number of runs per energy
num_runs_SameEnergy= 1
path_To_Input = '~/buuinput'
DoPrLevel(1) = .FALSE.
DoPrLevel(2) = .FALSE.
$end
$initRandom
SEED=45678 ! Seed for the random number
$end
$initDensity
densitySwitch=2 ! 1=dynamic, 2=analytic
$end
$initPauli
pauliSwitch=2 ! 1=dynamic, 2=analytic
$end
$propagation
coulomb=.false. ! Whether to use coulomb in propagation
predictorCorrector=.false. ! Whether to use a predictor/corrector algorithm
$end
!************************************************************
! Input for potentials
!************************************************************
$Coulomb
CoulombFlag=.false.
$end
$mesonPotential
noPerturbativePotential=.true. ! perturbative mesons feel no potential
$end
$baryonPotential
noPerturbativePotential=.true. ! perturbative baryons feel no potential
$end
$Yukawa
yukawaFlag=.false. ! whether Yukawa is switched on/off
$end
!************************************************************
! Input specific for the reactions
!************************************************************
$target
! fermiMotion = .FALSE.
target_Z= 1, target_A= 1 ! deuterium (proton,neutron)
$end
!************************************************************
$HiLeptonNucleus ! EVENTTYPE = 14
shadow=.false.
iExperiment=0 ! no experiment/fixed kinematics
$end
$HiGammaNucleus
$end
$HiPhotonKinematics
eps = 0.99, Q2 = 0.01, nu = 5.75
! eps = 0.1, Q2 = 0.01, nu = 5.75
$end
!************************************************************
! The collision term
!************************************************************
$hadronFormation
useJetSetVec = .FALSE.
$end
$collisionTerm
energyCheck=0.100 ! accuracy of energy check in GeV
twoBodyProcessesRealReal = .FALSE.
! oneBodyProcesses=.false.
! twoBodyProcesses=.false.
! threeBodyProcesses=.false.
$end
$insertion
minimumEnergy=0.100 ! minimal kinetic energy of produced nucleons (GeV)
$end
$master_2Body
correctEnergy_message=.false.
$end
!************************************************************
! Temperature and thermodynamics
!************************************************************
$initThermoDynamics
temperatureSwitch=1 ! 1=groundstate calculations (T=0,mu=E_F)
$end
!************************************************************
! Collision Statistics
!************************************************************
$collReporter
UseCollReporter=.FALSE.
$end
!************************************************************
! Additional Pythia Parameters
!************************************************************
$pythia
PARP( 91)=0.44 ! width intrinsic kT
$end