gibuu is hosted by Hepforge, IPPP Durham
GiBUU

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


Last modified 9 years ago Last modified on Jun 9, 2010, 2:49:32 PM