gibuu is hosted by Hepforge, IPPP Durham
GiBUU

jobcards/T2K_0pi_water: eO_1080_32.job

File eO_1080_32.job, 15.9 KB (added by gallmei, 6 years ago)

Jobcard for electronO

Line 
1!-------------------------------------------------- -*- fortran -*- ----
2! sample jobcard for lepton induced
3!-----------------------------------------------------------------------
4! please visit the GiBUU homepage for further information:
5! https://gibuu.hepforge.org
6!-----------------------------------------------------------------------
7
8!-----------------------------------------------------------------------
9! This jobcard is for fixed Enu, fixed scattering angle of the outgoing
10! lepton and output produced for various energy of outgoing lepton
11! (dsigma/dcostheta/dEprime xsec nuXsectionMode=1)
12!
13! see comments to the namelist &neutrino_induced for further
14! differential xsec (nuXsectionMode=0-5,7)
15!-----------------------------------------------------------------------
16! This jobcard would be useful, e.g., for producing double differential
17! xsec for inclusive electron scattering
18!-----------------------------------------------------------------------
19! In this jobcard most of the analysis for final states is switched off,
20! cf. namelist &neutrinoAnalysis;
21! The propagation of particles out of the nucleus is also switched off
22! (cf. numTimeSteps=0),
23! thus we are doing inclusive analysis
24! The event output is switched OFF
25! for generation of events better use nuXsectionMode=6
26! (16 for experimental flux)  and another jobcard as a prototype
27!-----------------------------------------------------------------------
28
29************************************************************************
30! DON'T FORGET to change "path_to_input" in the namelist &input according
31! to your local GiBUU installation
32************************************************************************
33
34
35!################################################
36!init flags
37!################################################
38!****************
39! neutrino init:
40!****************
41
42! file: code/init/neutrino/initNeutrino.f90
43&neutrino_induced
44       process_ID      =  1 ! 1:EM, 2:CC, 3:NC, -1:antiEM, -2:antiCC, -3:antiNC
45       flavor_ID       =  1 ! 1:electron, 2:muon, 3:tau
46
47!       nuXsectionMode  =  0 !  0: integratedSigma
48      nuXsectionMode  =  1 !  1: dSigmadCosThetadElepton
49!       nuXsectionMode  =  2 !  2: dSigmadQsdElepton
50!       nuXsectionMode  =  3 !  3: dSigmadQs
51!       nuXsectionMode  =  4 !  4: dSigmadCosTheta
52!       nuXsectionMode  =  5 !  5: dSigmadElepton
53!       nuXsectionMode  =  6 !  6: dSigmaMC
54!       nuXsectionMode  =  7 !  6: dSigmadW
55!
56! calculation for flux-veraged cross sections (choose with flag nuExp):
57!        above (from 1 to 6) plus 10
58!
59!       nuXsectionMode  = 11 ! 11: EXP_dSigmadCosThetadElepton
60!       nuXsectionMode  = 12 ! 12: EXP_dSigmadQsdElepton
61!       nuXsectionMode  = 13 ! 13: EXP_dSigmadQs
62!       nuXsectionMode  = 14 ! 14: EXP_dSigmadCosTheta
63!       nuXsectionMode  = 15 ! 15: EXP_dSigmadElepton
64!       nuXsectionMode  = 16 ! 16: EXP_dSigmaMC
65!       nuXsectionMode  = 16 ! 16: EXP_dSigmadW
66
67!       
68! The values for nuExp are ! only relevant if nuXsectionmode.ge.10
69!         
70       nuExp           =  0 !  0: noExp
71!       nuExp           =  1 !  1: MiniBooNE-nu
72!       nuExp           =  2 !  2: ANL
73!       nuExp           =  3 !  3: K2K
74!       nuExp           =  4 !  4: BNL
75!       nuExp           =  5 !  5: MinibooNE-barnu
76!       nuExp           =  6 !  6: MINOS-numu-numode
77!       nuExp           =  7 !  7: MINOS-barnumu-numode
78!       nuExp           =  8 !  8: NOvA
79!       nuExp           =  9 !  9: T2K-OA2.5-ND280
80!       nuExp           = 10 ! 10: uniform-distribution
81!       nuExp           = 11 ! 11: MINOS-numu-barnumode
82!       nuExp           = 12 ! 12: MINOS-barnumu-barnumode
83
84       includeQE       = T
85       includeDELTA    = T
86       includeRES      = T
87       include1pi      = T
88       includeDIS      = T
89       include2p2hQE   = T
90       include2p2hDelta= F
91       include2pi      = T 
92       
93       printAbsorptionXS = T
94/
95
96
97! file: code/density/nucleus.f90
98&target     
99      target_Z=8
100      target_A=16
101      densitySwitch_Static=2    ! 0: density=0.0, 1: Woods-Saxon by Lenske, 2 : NPA 554, 3: Woods-Saxon by Lenske, different neutron and proton radii,
102                                ! 5: density distribution is a sphere with density according to the input value of "fermiMomentum_input".
103      fermiMomentum_input=0.225 ! Input value of the fermi momentum for densitySwitch_static=5.
104      fermiMotion=.true.
105      ReAdjustForConstBinding=.true.
106!   ConstBinding=-0.008             ! only valid for ReAdjustForConstBinding=true
107/
108
109!#############################################################
110! general input, controls among others the statistics of output
111!#############################################################
112
113
114! the number of generated events is proportional (but not equal, because some "events" can be in prohibited phase space)
115! to target_A * numEnsembles * num_runs_SameEnergy
116! file: code/inputOutput/input.f90
117&input
118      numEnsembles=2000      ! for C12 you can use 4000, for heavier nuclei should be lower; decrease it if there are problems with memory
119      eventtype=5               ! 5=neutrino or electron, 200=simple transport of a given particle
120      numTimeSteps=0            ! the distance numTimeSteps*delta_T (100*0.2=20 fm in this example) should significantly exceed
121                                ! the radius of the target nucleus
122      delta_T=0.2
123      fullensemble=.true.
124      localEnsemble=.true.
125      num_runs_SameEnergy=1     ! increase these if you want to increase statistics (= number of generated events)
126      num_Energies=30      ! for nuXsectionMode=6 or 16 keep this 1, for other nuXsectionMode this is the number of steps in the
127                           !running Q2/Elepton/cosTheta variables, resp.
128      printParticleVectors=.false.
129      freezeRealParticles=.true.
130      LRF_equals_CALC_frame=.true. !if .false.: no offshelltransport possible
131      path_to_input='/home/mosel/GiBUU/buuinput'  ! for local run cluster
132/
133
134! file: code/density/density.f90
135&initDensity
136      densitySwitch=2           ! 1=dynamic density according to testparticle density, 2=analytic density prescription
137/
138
139! file: code/density/pauliBlocking.f90
140&initPauli
141      pauliSwitch     = 2 ! 0: no Pauli blocking, 1: dynamic, 2: analytic
142/
143
144
145! file: code/init/neutrino/neutrinoXsection.f90
146&nl_neutrinoxsection
147      singlePiModel=1           ! 0 for HNV, 1 for MAID
148      invariantMasscut=5.     ! cut events with invariant mass above
149      REScutW1=2.0              ! these 4 parameters describe the transition region in Invariant mass W
150                                ! between the resonances and DIS, see "Lalakulich Gallmeister Mosel PRC86(2012)014607" for details
151      REScutW2=2.05
152      DIScutW1=1.6
153      DIScutW2=1.65
154!      mcutDIS = 0.6
155     
156!      DISmassless = .true.
157/
158
159
160! SPECFIC NEUTRINO FLAGS depending on nuXsectionMode: (choose one)
161
162
163! relevant for nuXsectionMode=0
164! file: code/init/neutrino/neutrinoXsection.f90
165!&nl_integratedSigma
166!   enu=
167!   delta_enu=
168!/
169
170! relevant for nuXsectionMode=1 and (except parameter enu) 11
171! file: code/init/neutrino/neutrinoXsection.f90
172&nl_dSigmadCosThetadElepton     
173      enu      = 1.080  ! neutrino energy is fixed to 1 GeV
174      costheta = 0.84805 ! cos(outgoing-lepton angle scattering angle with respect to neutrino direction) is fixed
175      elepton  = 0.35 ! start calculations with energy of the outgoing lepton
176      delta_elepton = 0.025      ! make steps of 0.020; 
177!     for number of steps set parameter num_Energies in the namelist &input
178/
179     
180
181! relevant for nuXsectionMode=2 and (except parameter enu) 12
182! file: code/init/neutrino/neutrinoXsection.f90
183!&nl_dSigmadQsdElepton
184!   enu=
185!   Qs=
186!   elepton=
187!   delta_elepton=
188!/
189
190
191! relevant for nuXsectionMode=3 and (except parameter enu) 13
192! file: code/init/neutrino/neutrinoXsection.f90
193!&nl_dSigmadQs
194!   enu=
195!   Qs=
196!   delta_Qs=
197!/
198
199! relevant for nuXsectionMode=4  and (except parameter enu) 14
200! file: code/init/neutrino/neutrinoXsection.f90
201!&nl_dSigmadCosTheta
202!   enu=
203!   costheta=
204!   delta_costheta=
205!/
206
207! relevant for nuXsectionMode=5 and (except parameter enu) 15
208! file: code/init/neutrino/neutrinoXsection.f90
209!&nl_dSigmadElepton
210!   enu=
211!   elepton=
212!   delta_elepton=
213!/
214
215! relevant for nuXsectionMode=6 and (except parameter enu) 16
216! file: code/init/neutrino/neutrinoXsection.f90
217&nl_SigmaMC
218      MC_xmax=2.0
219      nu=1.0
220!/
221
222! relevant for nuXsectionMode=7
223&nl_dSigmadW
224       enu = 4.629
225       W = 0.6
226       delta_W = 0.16
227
228
229!********************************************
230! Model ingredients
231!******************************************
232
233! file: code/init/lowElectron/lepton2p2h.f90
234&lepton2p2h
235      ME_Version=4            ! This is the Christy model
236/
237
238
239! file: code/init/neutrino/matrixelement.f90
240&neutrino_matrixelement
241      which_resonanceModel=0    !0=QE + matrixelements from MAID, 1=QE matrixelements + old Delta, 2=Rein-Sehgal
242/
243
244
245! file: code/init/lepton/formfactors_QE_nucleon/FF_QE_nucleonScattering.f90
246&ff_QE
247      parametrization=3         ! 1=BBA03, 2=BBBA05, 3=BBBA07
248      useNonStandardMA=.false.  ! if true, use value of MA_in for axial mass MA, if false, use best fit
249!   MA_in=1.0
250/
251
252!****************
253! form factors for resonance production
254!****************
255! file: code/init/lepton/formfactors_ResProd/formFactor_ResProd.f90
256&input_FF_ResProd
257      FF_ResProd=0              ! 0=MAID in CM-frame, 1=fit of Lalakulich ,2=MAID in LAB-frame
258  !MA=1.3                         ! axial mass in the Delta resonance form factors based on ANL=0.95  or BNL=1.3(default)
259/
260
261!*****************
262!background parameters for fit to ANL or BNL pion data
263!*****************
264! file: code/init/neutrino/singlePionProductionMAIDlike.f90
265&neutrino_MAIDlikeBG
266!  b_proton_pinull=6.0            ! parameters for 1-pion background  ANL=3  (def)      BNL=6
267!  b_neutron_piplus=3.0           !                                   ANL=1.5  (def)    BNL=3.0
268/
269
270
271
272!################################################
273!analysis flags
274!################################################
275
276
277!****************
278! Neutrino analysis
279!****************
280
281! file: code/analysis/neutrinoAnalysis.f90
282&neutrinoAnalysis
283      radialScale=0                   
284      detailed_diff_output=.false. ! differential xsec; see namelist detailed_diff for max values and bins of the histograms
285      kineticEnergyDetectionThreshold_nucleon=0.
286      kineticEnergyDetectionThreshold_chargedpion=0.
287      kineticEnergyDetectionThreshold_neutralpion=0.
288      kineticEnergyDetectionThreshold_lepton=0.0
289      calorimetric_analysis=.false. ! calorimetric reconstruction of energy; see namelist nl_calorimetric_analysis
290      ZeroPion_analysis=.false. ! extra cross sections  for events with 0 pions in the final state
291      reconstruct_neutrino_energy=.false. ! calculate neutrino energy reconstruction for varios final states (see nl_specificEvent)
292      specificEvent_analysis=.false. ! extra cross sections for specific final states; see namelist nl_specificEvent
293                                                                                 ! also provides Q2 and energy reconstruction
294      inclusiveAnalysis=.true.  !if .true. we don't care whether particles have made it out       
295      outputEvents=.false.      ! output list of events and all outgoing particles in each event to the file FinalEvents.dat
296/
297
298
299
300! file: code/analysis/neutrinoAnalysis.f90
301&nl_calorimetric_analysis
302      numax=30
303      nubin=0.02
304      Enumax=30
305      enubin=0.02
306/
307
308! file: code/analysis/neutrinoAnalysis.f90
309$nl_specificEvent          ! any .true. in this namelist must be combines with calorimetric_analysis=.true. in $neutrinoAnalysis
310      no_pi=.false.              ! specificEvent=1        set to .false. if you want produce e.g. one-diff xsec versus lepton varaibles (Ekin, Q2, costheta)
311      p_Xn_no_pi=.false.         ! specificEvent=2
312      piplus=.false.             ! specificEvent=3
313      pi0=.false.                ! specificEvent=4
314      pi0_MULTI=.false.          ! specificEvent=5
315      piplus_MULTI=.false.       ! specificEvent=6
316      pp_no_pi=.false.          ! specificEvent=7
317      pn_no_pi=.false.          ! specificEvent=8
318      nn_no_pi=.false.          ! specificEvent=9
319      pp_Xn_no_pi=.false.       ! specificEvent=10
320      nn_Xp_no_pi=.false.       ! specificEvent=11
321      ppp_Xn_no_pi=.false.      ! specificEvent=12
322      pppp_Xn_no_pi=.false.     ! specificEvent=13
323      p_no_pi=.false.           ! specificEvent=14
324      n_no_pi=.false.           ! specificEvent=15
325      Xn_no_pi=.false.          ! specificEvent=16
326/
327
328
329! file: code/analysis/neutrinoAnalysis.f90
330&detailed_diff
331      ekinMax=3.0               ! maximum value of Ekin in the output of distribution of hadrons versus kinetic energy
332      dEkin=0.02                ! binning of Ekin in the output versus kinetic energy
333      fornucleon=.false.         ! produce output versus kinetic energy, angle, ... for this outgoing hadron
334      forpion=.false.            ! to the files diff_....
335      foreta=.false.
336      forkaon=.false.
337      forkaonBar=.false.
338      forLambda=.false.
339      forSigmaResonance=.false.
340/
341
342
343
344
345
346
347!################################################
348! in-medium width
349!################################################
350
351! file: code/width/baryonWidthMedium.f90
352&width_Baryon
353      mediumSwitch=.false.       ! if .false. vacuum widths will be used for all resonances
354                                      ! and other switches are not meaningful/relevant
355      mediumSwitch_Delta=.false. ! Use Oset broadening of the Delta (only relevant for mediumSwitch_coll=.false.)
356      mediumSwitch_coll=.false. ! Use collisional broadening of all resonances
357/
358
359
360! file: code/spectralFunctions/spectralFunc.f90
361&spectralFunction
362      nuclwidth=0.001           ! if mediumSwitch_coll=.false.: this
363                                !parameter determines the "bare nucleon width", for numerical purpose
364      relativistic=.true.       ! if mediumSwitch_coll=.false.: relativistic or non-rel. spectral function
365/
366
367!################################################
368! propagation
369!################################################
370
371! file:  code/propagation/propagation.f90   
372&propagation
373      delta_P=0.01              ! Delta Momentum for derivatives
374      delta_E=0.01              ! Delta Energy for derivatives
375 
376      UseHadronic=.true.        ! Whether to use hadronic potentials in propagation
377      RungeKuttaOrder=2               ! 1=first order Runge-Kutta, 2=second order Runge-Kutta
378      Mode=2                          ! =0 Cascade, =1 Euler, =2 Predictor-Corrector   
379/
380
381
382! file: code/width/offShellPotential.f90
383&offShellPotential
384      useOffShellPotentialBaryons=.false. ! if inMediumSwitch=.false. also useOffShellpotential will be set to .false.
385      extrapolateBaryonWidth=.false. ! whether to extrapolate the baryon width below minimal mass
386/
387
388
389!################################################
390! potentials
391!################################################
392
393! file: code/potential/baryonPotential.f90
394&baryonPotential
395      EQS_Type=5                ! 1=soft mom-dep, 2=hard mom-dep, 3=soft non-mom-dep, 4=hard non-mom-dep, 5=medium mom-dep
396      DeltaPot=1                ! 1=2/3 of nucleon potential, 2=100 MeV*rho/rhoNull
397      noPerturbativePotential=.false.
398      symmetryPotFlag=0  ! Switch for the assymetry term in the nucleon potential
399/
400
401
402! file: code/potential/coulomb/coulomb.f90
403&Coulomb
404      CoulombFlag=.false.
405/
406
407
408!################################################
409! collision term
410!################################################
411
412! file: code/collisions/collisionTerm.f90
413&collisionTerm
414      oneBodyProcesses=.true.
415      twoBodyProcesses=.true.
416      threeBodyProcesses=.true.
417      DoJustAbsorptive=.false.
418/
419
420
421! file: code/collisions/twoBodyReactions/master_2Body.f90
422&master_2Body
423      baryonBaryonScattering=.true.
424      baryonMesonScattering=.true.
425      mesonMesonScattering=.false.
426/
427
428
429! file: code/collisions/twoBodyReactions/baryonMeson/resonanceCrossSections.f90
430&resonanceCrossSections
431      fullPropagator=.false.    ! Use self energies in resonance propagators (only when collisional broadening is switched on)
432/
433
434
435
436! file: code/collisions/insertion.f90
437&insertion
438      minimumEnergy=0.005       ! default is 0.005 GeV
439/
440
441
442! file: code/collisions/twoBodyReactions/HiEnergy/DoCollTools.f90
443&pythia
444      PARP(91)=0.44
445/