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 14 years ago
Last modified on Jun 9, 2010, 2:49:32 PM