Changes between Initial Version and Version 1 of jobcards/jobPhotoProd

Timestamp:

Jun 9, 2010, 2:49:32 PM (14 years ago)

Author:

gallmei

Comment:

--

jobcards/jobPhotoProd

@@ +1 @@
+ = 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
+
+
+}}}