TABLE OF CONTENTS
- 1. /neutrinoAnalysis
- 1.1. neutrinoAnalysis/detailed_diff_output
- 1.2. neutrinoAnalysis/forPion
- 1.3. neutrinoAnalysis/forEta
- 1.4. neutrinoAnalysis/forKaon
- 1.5. neutrinoAnalysis/forKaonBar
- 1.6. neutrinoAnalysis/forDmeson
- 1.7. neutrinoAnalysis/forDbar
- 1.8. neutrinoAnalysis/forDs_plus
- 1.9. neutrinoAnalysis/forDs_minus
- 1.10. neutrinoAnalysis/forNucleon
- 1.11. neutrinoAnalysis/forLambda
- 1.12. neutrinoAnalysis/forSigmaResonance
- 1.13. neutrinoAnalysis/forXi
- 1.14. neutrinoAnalysis/forOmegaResonance
- 1.15. neutrinoAnalysis/include_W_dist
- 1.16. neutrinoAnalysis/dW_Npi
- 1.17. neutrinoAnalysis/Wmin_Npi
- 1.18. neutrinoAnalysis/Wmax_Npi
- 1.19. neutrinoAnalysis/dW_mupi
- 1.20. neutrinoAnalysis/Wmin_mupi
- 1.21. neutrinoAnalysis/Wmax_mupi
- 1.22. neutrinoAnalysis/dW_muN
- 1.23. neutrinoAnalysis/Wmin_muN
- 1.24. neutrinoAnalysis/Wmax_muN
- 1.25. neutrinoAnalysis/numin
- 1.26. neutrinoAnalysis/numax
- 1.27. neutrinoAnalysis/nubin
- 1.28. neutrinoAnalysis/Enumin
- 1.29. neutrinoAnalysis/Enumax
- 1.30. neutrinoAnalysis/Enubin
- 1.31. neutrinoAnalysis/EkinMin
- 1.32. neutrinoAnalysis/EkinMax
- 1.33. neutrinoAnalysis/dEkin
- 1.34. neutrinoAnalysis/EkinMin_lepton
- 1.35. neutrinoAnalysis/EkinMax_lepton
- 1.36. neutrinoAnalysis/dEkin_lepton
- 1.37. neutrinoAnalysis/cost_min
- 1.38. neutrinoAnalysis/cost_max
- 1.39. neutrinoAnalysis/delta_cost
- 1.40. neutrinoAnalysis/kineticEnergyDetectionThreshold_nucleon
- 1.41. neutrinoAnalysis/kineticEnergyDetectionThreshold_chargedpion
- 1.42. neutrinoAnalysis/kineticEnergyDetectionThreshold_neutralpion
- 1.43. neutrinoAnalysis/kineticEnergyDetectionThreshold_lepton
- 1.44. neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_nucleon
- 1.45. neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_chargedpion
- 1.46. neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_neutralpion
- 1.47. neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_lepton
- 1.48. neutrinoAnalysis/radialScale
- 1.49. neutrinoAnalysis/inclusiveAnalysis
- 1.50. neutrinoAnalysis/Fissum_analysis
- 1.51. neutrinoAnalysis/ZeroPion_analysis
- 1.52. neutrinoAnalysis/calorimetric_analysis
- 1.53. neutrinoAnalysis/reconstruct_neutrino_energy
- 1.54. neutrinoAnalysis/specificEvent_analysis
- 1.55. neutrinoAnalysis/no_pi
- 1.56. neutrinoAnalysis/p_Xn_no_pi
- 1.57. neutrinoAnalysis/piplus
- 1.58. neutrinoAnalysis/piplus_MULTI
- 1.59. neutrinoAnalysis/pi0
- 1.60. neutrinoAnalysis/pi0_MULTI
- 1.61. neutrinoAnalysis/piminus
- 1.62. neutrinoAnalysis/piminus_MULTI
- 1.63. neutrinoAnalysis/pp_no_pi
- 1.64. neutrinoAnalysis/pn_no_pi
- 1.65. neutrinoAnalysis/nn_no_pi
- 1.66. neutrinoAnalysis/pp_Xn_no_pi
- 1.67. neutrinoAnalysis/nn_Xp_no_pi
- 1.68. neutrinoAnalysis/ppp_Xn_no_pi
- 1.69. neutrinoAnalysis/pppp_Xn_no_pi
- 1.70. neutrinoAnalysis/p_no_pi
- 1.71. neutrinoAnalysis/n_no_pi
- 1.72. neutrinoAnalysis/Xn_no_pi
- 1.73. neutrinoAnalysis/excl_hadron
- 1.74. neutrinoAnalysis/QEp
- 1.75. neutrinoAnalysis/pcut
- 1.76. neutrinoAnalysis/binsizeQ2
- 1.77. neutrinoAnalysis/binsizeEnu
- 1.78. neutrinoAnalysis/maxQ2
- 1.79. neutrinoAnalysis/maxEnu
- 1.80. neutrinoAnalysis/outputEvents
- 1.81. neutrinoAnalysis/Xsection_analysis
- 1.82. neutrinoAnalysis/readinput
- 1.83. neutrinoAnalysis/NeutrinoAnalysis
- 1.84. neutrinoAnalysis/W_distributions
- 1.85. neutrinoAnalysis/nl_calorimetric_analysis
- 1.86. neutrinoAnalysis/nl_specificEvent
- 1.87. neutrinoAnalysis/detailed_diff
- 1.88. neutrinoAnalysis/neutrino_Analyze
- 1.89. neutrinoAnalysis/neutrino_total_Xsection_multiplicities.dat
- 1.90. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_QE.dat
- 1.91. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_Delta.dat
- 1.92. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_highRES.dat
- 1.93. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_1piBG.dat
- 1.94. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_DIS.dat
- 1.95. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_2p2h.dat
- 1.96. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen0.dat
- 1.97. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen1.dat
- 1.98. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen2.dat
- 1.99. neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen3ormore.dat
- 1.100. neutrinoAnalysis/neutrino_kinetic_energy_Xsection.dat
- 1.101. neutrinoAnalysis/neutrino_Xsection_numbers.dat
- 1.102. neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat
- 1.103. neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_1X.dat
- 1.104. neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_2X.dat
- 1.105. neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_MULTI.dat
- 1.106. neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE.dat
- 1.107. neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dTheta_HADRON_charge_CHARGE.dat
- 1.108. neutrinoAnalysis/diff_XXX_dSigma_dEkin_lepton_PPP.dat
- 1.109. neutrinoAnalysis/diff_XXX_dSigma_dQ2_lepton_PPP.dat
- 1.110. neutrinoAnalysis/diff_XXX_dSigma_dQ2p_lepton_PPP.dat
- 1.111. neutrinoAnalysis/diff_XXX_dSigma_dcos_lepton_PPP.dat
- 1.112. neutrinoAnalysis/diff_XXX_d2Sigma_dEdcost_lepton_no_pi.dat
- 1.113. neutrinoAnalysis/diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat
- 1.114. neutrinoAnalysis/diff_XXX_dSigma_dW_muon_nucleon_charge_CHARGE.dat
- 1.115. neutrinoAnalysis/diff_XXX_dSigma_dW_muon_pion_charge_CHARGE.dat
- 1.116. neutrinoAnalysis/FinalEvents.dat
- 1.117. neutrinoAnalysis/reconstruction_Enureal_PPP.ZZZ.dat
- 1.118. neutrinoAnalysis/reconstruction_Enurec_PPP.ZZZ.dat
- 1.119. neutrinoAnalysis/reconstruction_Enu_rec_versus_real_PPP.ZZZ.dat
- 1.120. neutrinoAnalysis/reconstruction_Q2real_PPP.ZZZ.dat
- 1.121. neutrinoAnalysis/reconstruction_Q2rec_PPP.ZZZ.dat
- 1.122. neutrinoAnalysis/reconstruction_Q2_rec_versus_real_PPP.ZZZ.dat
- 1.123. neutrinoAnalysis/oscillations_UUU_real_PPP.ZZZ.dat
- 1.124. neutrinoAnalysis/oscillations_UUU_rec_PPP.ZZZ.dat
- 1.125. neutrinoAnalysis/IsBound
- 1.126. neutrinoAnalysis/IsBelowThreshold
- 1.127. neutrinoAnalysis/oscillationProbability
- 1.128. neutrinoAnalysis/PrintVals10
- 1.129. neutrinoAnalysis/PrintHeader10
/neutrinoAnalysis [ Modules ]
NAME
module neutrinoAnalysis
PURPOSE
neutrinoAnalysis/detailed_diff_output [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: detailed_diff_output = .false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forPion [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forPion =.true.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forEta [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forEta =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forKaon [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forKaon =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forKaonBar [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forKaonBar =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forDmeson [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forDmeson =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forDbar [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forDbar =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forDs_plus [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forDs_plus =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forDs_minus [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forDs_minus =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forNucleon [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forNucleon =.true.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forLambda [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forLambda =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forSigmaResonance [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forSigmaResonance =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forXi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forXi =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/forOmegaResonance [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: forOmegaResonance =.false.
PURPOSE
If .true. then also the detailed output of differential cross sections is produced
neutrinoAnalysis/include_W_dist [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: include_W_dist=.false.
PURPOSE
If .true. then the invariant mass distributions for events with 1 pion and 1 nucleon in the final state are produced
neutrinoAnalysis/dW_Npi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: dW_Npi=0.02
PURPOSE
for dsigma/d(InvariantMass); only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/Wmin_Npi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Wmin_Npi=1.08
PURPOSE
for dsigma/d(InvariantMass); only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/Wmax_Npi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Wmax_Npi=1.6
PURPOSE
for dsigma/d(InvariantMass); only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/dW_mupi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: dW_mupi=0.04
PURPOSE
only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/Wmin_mupi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Wmin_mupi=0.24
PURPOSE
only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/Wmax_mupi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Wmax_mupi=1.2
PURPOSE
only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/dW_muN [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: dW_muN=0.04
PURPOSE
only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/Wmin_muN [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Wmin_muN=1.04
PURPOSE
only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/Wmax_muN [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Wmax_muN=2.12
PURPOSE
only work if include_W_dist is .true. set the min, max and steps for various W-distributions
neutrinoAnalysis/numin [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: numin=0.
PURPOSE
for calorimetric analysis: values for transferred energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions
neutrinoAnalysis/numax [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: numax=10.0
PURPOSE
for calorimetric analysis: values for transferred energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions
neutrinoAnalysis/nubin [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: nubin=0.1
PURPOSE
for calorimetric analysis: values for transferred energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions
neutrinoAnalysis/Enumin [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Enumin=0.
PURPOSE
for calorimetric analysis: values for neutrino energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions
neutrinoAnalysis/Enumax [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Enumax=10.0
PURPOSE
for calorimetric analysis: values for neutrino energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions
neutrinoAnalysis/Enubin [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: Enubin=0.1
PURPOSE
for calorimetric analysis: values for neutrino energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions
neutrinoAnalysis/EkinMin [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: EkinMin=0.
PURPOSE
if detailed_diff_output is TRUE: Maximal kinetic energy for dsigma/dEkin for hadrons
neutrinoAnalysis/EkinMax [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: EkinMax=2.
PURPOSE
if detailed_diff_output is TRUE: Maximal kinetic energy for dsigma/dEkin for hadrons
neutrinoAnalysis/dEkin [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: dEkin=0.01
PURPOSE
if detailed_diff_output is TRUE: Delta(eKin) for dsigma/dEKin for hadrons
neutrinoAnalysis/EkinMin_lepton [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: EkinMin_lepton=0.
PURPOSE
if detailed_diff_output is TRUE: Maximal kinetic energy for dsigma/dEkin for leptons
neutrinoAnalysis/EkinMax_lepton [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: EkinMax_lepton=2.
PURPOSE
Maximal kinetic energy for dsigma/dEkin for leptons
neutrinoAnalysis/dEkin_lepton [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: dEkin_lepton=0.01
PURPOSE
if detailed_diff_output is TRUE: Delta(eKin) for dsigma/dEKin for leptons
neutrinoAnalysis/cost_min [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: cost_min= -1.
PURPOSE
Minimal cos(theta) of outgoing leptons
neutrinoAnalysis/cost_max [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: cost_max= +1
PURPOSE
Maximal cos(theta) of outgoing leptons
neutrinoAnalysis/delta_cost [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: delta_cost = 0.02
PURPOSE
stepsize of cos(theta) of outgoing leptons
neutrinoAnalysis/kineticEnergyDetectionThreshold_nucleon [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: kineticEnergyDetectionThreshold_nucleon=0.0
PURPOSE
kineticEnergyDetectionThreshold lower detection threshold for nucleon kinetic energies
neutrinoAnalysis/kineticEnergyDetectionThreshold_chargedpion [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: kineticEnergyDetectionThreshold_chargedpion=0.0
PURPOSE
kineticEnergyDetectionThreshold
neutrinoAnalysis/kineticEnergyDetectionThreshold_neutralpion [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: kineticEnergyDetectionThreshold_neutralpion=0.0
PURPOSE
kineticEnergyDetectionThreshold
neutrinoAnalysis/kineticEnergyDetectionThreshold_lepton [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: kineticEnergyDetectionThreshold_lepton=0.0
PURPOSE
kineticEnergyDetectionThreshold only lepton kinetic energies above this threshold can be detected
neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_nucleon [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: AngleUpperDetectionThresholdDegrees_nucleon=180.0
PURPOSE
nucleon angles up to this value can be detected
neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_chargedpion [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: AngleUpperDetectionThresholdDegrees_chargedpion=180.0
PURPOSE
charged pion angles up to this value can be detected
neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_neutralpion [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: AngleUpperDetectionThresholdDegrees_neutralpion=180.0
PURPOSE
neutral pion angles angles up to this value can be detected
neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_lepton [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: AngleUpperDetectionThresholdDegrees_lepton=180.0
PURPOSE
lepton angles up to this value can be detected
neutrinoAnalysis/radialScale [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: radialScale=0.0
PURPOSE
If radial position of nucleon < radialScale, then the nucleon is assumed to be bound
neutrinoAnalysis/inclusiveAnalysis [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: inclusiveAnalysis=.false.
PURPOSE
If .true. then we don't care whether particle has left the nucleus or not
neutrinoAnalysis/Fissum_analysis [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: Fissum_analysis=.false.
PURPOSE
do analysis with cuts as needed for Fig 25 in Fissum et al, PRC 70, 034606 (2004)
neutrinoAnalysis/ZeroPion_analysis [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: ZeroPion_analysis=.false.
PURPOSE
produce output of xsec for various final states with 0 pions and 2 pions see file see sigma_0pions.dat for the list of the final states
see files neutrino_0pions.dat, neutrino_0pions_QE.dat, neutrino_0pions_Delta.dat, ... for output
neutrinoAnalysis/calorimetric_analysis [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: calorimetric_analysis=.false.
PURPOSE
do calorimetric energy-transfer and neutrino-energy reconstruction (for each QE, Delta, ...) as in the MINOS experiment
neutrinoAnalysis/reconstruct_neutrino_energy [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: reconstruct_neutrino_energy=.false.
PURPOSE
reconstruct neutrino energy for final state in "specificEvent_analysis"
NOTES
.true. must be combined with specificEvent_analysis=.true. and at least one specific event .true.
neutrinoAnalysis/specificEvent_analysis [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: specificEvent_analysis=.false.
PURPOSE
do analysis for specific final states
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=1, no_pi (for example, for QE-like MiniBooNE)
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/p_Xn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: p_Xn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=2
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/piplus [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: piplus=.false.
PURPOSE
do analysis for specific final states: specificEvent=3, 1 pi+ X nucleons mesons of other flavor
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/piplus_MULTI [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: piplus_MULTI=.false.
PURPOSE
do analysis for specific final states: specificEvent=4 >=1 pi+ X other pions (incl pi+) X nucleons
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pi0 [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: pi0=.false.
PURPOSE
do analysis for specific final states: specificEvent=5, 1 pi0 X nucleons, plus mesons of other flavor
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pi0_MULTI [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: pi0_MULTI=.false.
PURPOSE
do analysis for specific final states: specificEvent=6, >=1 pi0 X other pions X nucleons, (pi0 K2K)
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/piminus [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: piminus=.false.
PURPOSE
do analysis for specific final states: specificEvent=7 1 pi- X other pions X nucleons
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/piminus_MULTI [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: piminus_MULTI=.false.
PURPOSE
do analysis for specific final states: specificEvent=8 >=1 pi- X other pions X nucleons
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pp_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: pp_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=9 2 protons, X neutrons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: pn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=10 1 neutron, 1 proton, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/nn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: nn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=11 2 neutrons, X protons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pp_Xn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: pp_Xn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=12 2 protons, X neutrons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/nn_Xp_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: nn_Xp_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=13 2 neutrons, X protons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/ppp_Xn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: ppp_Xn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=14 3 protons, X neutrons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pppp_Xn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: pppp_Xn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=15 4 protons, X neutrons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/p_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: p_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=16 1 proton, 0 neutron, 0 pion
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/n_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: n_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=17 1 neutron, 0 proton, 0 pion
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/Xn_no_pi [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: Xn_no_pi=.false.
PURPOSE
do analysis for specific final states: specificEvent=18, 0 proton, X neutrons, 0 pions
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/excl_hadron [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: excl_hadron=.false.
PURPOSE
do analysis for specific final states: specificEvent=19,20,21 exclusive 1 pion, no other pions or other mesons of different flavor There could be still other mesons which are heavier than the D, Such events (very rare at DUNE energies) could be counted as exclusive single-meson cross section. This could be cured by extending the list of stable mesons
value can be changed in the namelist nl_specificEvent
neutrinoAnalysis/QEp [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical :: QEp=.false.
PURPOSE
if .true, do analysis for specific analysis for QE-like event with 1 mu, 0 pi, X p
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/pcut [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: pcut = 0.0
PURPOSE
determines lower acceptance cut for outgoing protons
values can be changed in the namelist neutrinoAnalysis via the variable kineticEnergyDetectionThreshold_nucleon
neutrinoAnalysis/binsizeQ2 [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: binsizeQ2=0.01
PURPOSE
do analysis for specific final states: binning for reconstruction of Q2 and Enu
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/binsizeEnu [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: binsizeEnu=0.02
PURPOSE
do analysis for specific final states: binning for reconstruction of Q2 and Enu
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/maxQ2 [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: maxQ2=5.0
PURPOSE
do analysis for specific final states: max values for reconstruction of Q2 and Enu
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/maxEnu [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
real, save :: maxEnu=5.0
PURPOSE
do analysis for specific final states: max values for reconstruction of Q2 and Enu
values can be changed in the namelist nl_specificEvent
neutrinoAnalysis/outputEvents [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: outputEvents = .false.
PURPOSE
If .true. then all events are printed to the file 'FinalEvents.dat'.
neutrinoAnalysis/Xsection_analysis [ Global module-variables ]
[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]
SOURCE
logical, save :: Xsection_analysis = .false.
PURPOSE
If .true. then files "..._total_Xsection_..." are printed.
neutrinoAnalysis/readinput [ Subroutines ]
[ Top ] [ neutrinoAnalysis ] [ Subroutines ]
NAME
subroutine readinput
INPUTS
NONE
OUTPUT
NONE
PURPOSE
This subroutine reads the namelist "neutrinoAnalysis". Only called once to initialize the module.
neutrinoAnalysis/NeutrinoAnalysis [ Namelists ]
[ Top ] [ neutrinoAnalysis ] [ Namelists ]
NAME
NAMELIST /NeutrinoAnalysis/
PURPOSE
This namelist includes:
- detailed_diff_output
- include_W_dist
- kineticEnergyDetectionThreshold_lepton
- AngleUpperDetectionThresholdDegrees_lepton
- kineticEnergyDetectionThreshold_nucleon
- AngleUpperDetectionThresholdDegrees_nucleon
- kineticEnergyDetectionThreshold_chargedpion
- AngleUpperDetectionThresholdDegrees_chargedpion
- kineticEnergyDetectionThreshold_neutralpion
- AngleUpperDetectionThresholdDegrees_neutralpion
- inclusiveAnalysis
- Fissum_analysis
- ZeroPion_analysis
- calorimetric_analysis
- radialScale
- reconstruct_neutrino_energy
- outputEvents
- specificEvent_analysis
- Xsection_analysis
neutrinoAnalysis/W_distributions [ Namelists ]
[ Top ] [ neutrinoAnalysis ] [ Namelists ]
NAME
NAMELIST /W_distributions/
PURPOSE
This Namelist includes:
neutrinoAnalysis/nl_calorimetric_analysis [ Namelists ]
[ Top ] [ neutrinoAnalysis ] [ Namelists ]
NAME
NAMELIST /nl_calorimetric_analysis/
PURPOSE
This Namelist includes:
neutrinoAnalysis/nl_specificEvent [ Namelists ]
[ Top ] [ neutrinoAnalysis ] [ Namelists ]
NAME
NAMELIST /nl_specificEvent/
PURPOSE
This namelist includes:
- no_pi
- p_Xn_no_pi
- piplus
- piplus_MULTI
- pi0
- pi0_MULTI
- piminus
- piminus_MULTI
- pp_no_pi
- pn_no_pi
- nn_no_pi
- pp_Xn_no_pi
- nn_Xp_no_pi
- ppp_Xn_no_pi
- pppp_Xn_no_pi
- p_no_pi
- n_no_pi
- Xn_no_pi
- binsizeQ2
- binsizeEnu
- maxQ2
- maxEnu
- excl_hadron
- QEp
neutrinoAnalysis/detailed_diff [ Namelists ]
[ Top ] [ neutrinoAnalysis ] [ Namelists ]
NAME
NAMELIST /detailed_diff/
PURPOSE
This namelist includes:
- EkinMin
- EkinMax
- dEkin
- EkinMin_lepton
- EkinMax_lepton
- dEkin_lepton
- forPion
- forEta
- forKaon
- forKaonBar
- forDmeson
- forDbar
- forDs_plus
- forDs_minus
- forNucleon
- forLambda
- forSigmaResonance
- forXi
- forOmegaResonance
neutrinoAnalysis/neutrino_Analyze [ Subroutines ]
[ Top ] [ neutrinoAnalysis ] [ Subroutines ]
NAME
subroutine neutrino_Analyze(Particles,finalFlag,num_runs_sameEnergy)
INPUTS
- type(particle), intent(in),dimension(:,:) :: Particles -- Particles which shall be analyzed
- logical, intent(in) :: finalFlag -- if .true. than the final output for a series of calls will be done
- integer :: num_runs_sameEnergy
NOTES
neutrinoAnalysis/neutrino_total_Xsection_multiplicities.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities.dat
PURPOSE
The file is produced in the runs with eventtype=5=neutrino .
The file shows the cross sections for preselected final states
Units:
- For process_ID=CC and NC the units 10^{-38} cm^2 for integrated xsec (10^{-38)cm^2/GeV for dsigma/dElepton, 10^{-38)cm^2/GeV^2 for dsigma/dQ^2, and so on)
- For process_ID=EM the units are nanobarns=10^{-33}cm^2
Columns:
- #1: variable which was raised (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
- #2-#119: see description in AnaEvent.f90, subroutine event_sigma OR in the output file sigma.dat The description of some columns is given below In each channel the outgoing lepton is presupposed (unless explicitely stated otherwise)
Some columns in detail:
- #2: 1 pi-, no other pions of any charge and anything else
- #3: 1 pi0, no other pions of any charge and anything else
- #4: 1 pi+, no other pions of any charge and anything else
- #5: 1 eta and anything else
- #14: 1 neutron, no other nucleons and anything else
- #15: 1 proton, no other nucleons and anything else
- #68: 1 nucleon and 1 pion, no other pions or nucleons, anything else
- #68: 1 proton (no other nucleons) and 0 pions and anything else (QE-like in Cherenkov detector)
- #69: 0 pions and anything else
- #70: at least 1 pi-, any number of pi0 and/or pi+ , anything else, each event is counted once
- #71: at least 1 pi0, any number of pi- and/or pi+ , anything else, each event is counted once
- #72: at least 1 pi+, any number of pi- and/or pi0 , anything else, each event is counted once
- #91: at least 1 pi0, any number of pi- and/or pi+ , anything else, each pi0 is counted
- #96: no nucleons, anything else
- #97: 5 or more nucleons, anything else
- #120-#239: errors to Columns #2-#119
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_QE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_QE.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat but for QE events (=the first interaction was quasielastic or elastic scattering)
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_Delta.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_Delta.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat but for Delta events (=the first interaction was production of the Delta resonance)
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_highRES.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_highRES.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat but for highRES events (=the first interaction was production any resonance beyond Delta)
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_1piBG.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_1piBG.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat but for 1piBG events (=the first interaction was background production of 1-pion final state)
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_DIS.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_DIS.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat but for DIS events (=the first interaction was DIS)
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_2p2h.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_2p2h.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat but for 2particle-2hole events
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen0.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_gen0.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat for particles of the 0th generation The definition of generation is given in history.f90, description of the module "history"
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen1.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_gen1.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat for particles of the 1st generation The definition of generation is given in history.f90, description of the module "history"
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen2.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_gen2.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat for particles of the 2nd generation The definition of generation is given in history.f90, description of the module "history"
neutrinoAnalysis/neutrino_total_Xsection_multiplicities_gen3ormore.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_total_Xsection_multiplicities_gen3ormore.dat
PURPOSE
The same as neutrino_total_Xsection_multiplicities.dat for particles of the 3rd-or-more generation The definition of generation is given in history.f90, description of the module "history"
neutrinoAnalysis/neutrino_kinetic_energy_Xsection.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_kinetic_energy_Xsection.dat
PURPOSE
The file is produced in the runs with eventtype=5=neutrino .
The file shows the kinetic energy differential cross sections after final state interactions for pions, protons and neutrons
Units:
- For process_ID=CC and NC: 10^{-38} cm^2/GeV for integrated xsec and so on ..
- For process_ID=EM: nanobarns=10^{-33}cm^2/ GeV and so on.
Columns:
- #1: variable which was raised (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
- #2: kinetic energy in GeV
- #3: xsec for events with at least one pi- in the final state (coincides with column 2 in diff_XXX_dSigma_dEkin_pi_charge_-1_MULTI.dat)
- #4: xsec for events with at least one pi0 (coincides with column 2 in diff_XXX_dSigma_dEkin_pi_charge_+0_MULTI.dat)
- #5: xsec for events with at least one pi+ (coincides with column 2 in diff_XXX_dSigma_dEkin_pi_charge_+1_MULTI.dat)
- #6: xsec for events with at least one neutron (coincides with column 2 in diff_XXX_dSigma_dEkin_N_charge_+0_MULTI.dat)
- #7: xsec for events with at least one proton (coincides with column 2 in diff_XXX_dSigma_dEkin_N_charge_+1_MULTI.dat)
HERE XXX=000,001, ... is the count of the raise variable
neutrinoAnalysis/neutrino_Xsection_numbers.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file neutrino_Xsection_numbers.dat
PURPOSE
The file is produced in the runs with eventtype=5=neutrino .
The file shows number of nucleons produced in one run before and after final state interactions (should be devided by numEnsembles to obtain the average multiplicity per target nucleus and divided further by target_A to obtain the average multiplicity per target nucleon)
Columns:
- #1: variable which was raised (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
- #2: all the protons produced before final state interactions (FSI)
- #3: protons that made it out of the nucleus after FSI
- #4: protons produced in generation 0 and made it out of the nucleus after FSI
- #5: protons produced in generation 1 and made it out of the nucleus after FSI
- #6: protons produced in generation 2 and made it out of the nucleus after FSI
- #7: protons produced in generation 3-or-more and made it out of the nucleus after FSI
- #8: all the neutrons produced before FSI
- #9: neutrons that made it out of the nucleus after FSI
- #10: neutrons produced in generation 0 and made it out of the nucleus after FSI
- #11: neutrons produced in generation 1 and made it out of the nucleus after FSI
- #12: neutrons produced in generation 2 and made it out of the nucleus after FSI
- #13: neutrons produced in generation 3-or-more and made it out of the nucleus after FSI
- #14: all the antiprotons produced before FSI
- #15: antiprotons that made it out of the nucleus after FSI
- #16: antiprotons produced in generation 0 and made it out of the nucleus after FSI
- #17: antiprotons produced in generation 1 and made it out of the nucleus after FSI
- #18: antiprotons produced in generation 2 and made it out of the nucleus after FSI
- #19: antiprotons produced in generation 3-or-more and made it out of the nucleus after FSI
neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat
with:
- ZZZ denotes the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneSX.ZZZ.dat).
If ZZZ is missing the total cross section (sum over all primary events) is given, otherwise ZZZ=Delta, DIS, ....
* XXX=000, 001, 002 --- the first, second, third and so on values of
the "raised variable" (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
* HADRON = pion, N, K, K~, Lambda, SigmaResonance, eta
only those are shown, for which the switches in the namelist "detailed_diff" are set to .true. possible hadrons are contained in the field particleIDs defined in AnaEvent.f90
* CHARGE = charge of the outgoing hadron
PURPOSE
The file is produced in runs with:
- eventtype=5=neutrino if switch "detailed_diff_output" in namelist "neutrinoAnalysis" the is set to .true.
- eventtype=3=LowPhoto if switch "dE_switch" in namelist "LowElePhoto_Analysis" is set to .true.
The file shows the cross sections for ___1-HADRON___ final state versus kinetic energy of the outgoing hadron (one HADRON of a given CHARGE and no HADRONs with same flavor, but different charges;
there could be additional hadrons with different flavor)
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2/GeV
- For eventtype=3: microbarns=10^{-30}cm^2/GeV
Columns:
- #1: kinetic energy of the outgoing HADRON of given CHARGE [GeV]
- #2: dsi/dEkin xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_1X.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE_1X.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but for ___1-HADRON-X___ final state (one HADRON of a given CHARGE and any number of HADRONs of different charges)
neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_2X.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE_2X.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but for ___2-HADRON-X___ final state (two HADRONs of a given CHARGE and any number of HADRONs of different charges)
neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_MULTI.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE_MULTI.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but for ___MULTI-HADRON___ final state (at least one HADRON of a given CHARGE and any number of HADRONs of different charges)
neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE.dat
file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE_1X.dat
file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE_2X.dat
file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE_MULTI.dat)
PURPOSE
The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but dsigma/d(E(1-cosTheta))
The file is produced in runs with:
- eventtype=5=neutrino if switch "detailed_diff_output" in namelist "neutrinoAnalysis" the is set to .true.
The file shows the cross sections dsigma/d(E(1-cosTheta)) versus E*(1-costheta), where E is the energy (full energy, not kinetic) of the outgoing hadron, costheta its polar scattering angle (recall here that in neutrino runs neutrinos are moving along z-direction)
Columns:
- #1: E*(1-cosTheta) [GeV]
- #2: dsi/d(E(1-costheta)) xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dTheta_HADRON_charge_CHARGE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE.dat
file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE_1X.dat
file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE_2X.dat
file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE_MULTI.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but dsigma/dTheta
The file is produced in runs with:
- eventtype=5=neutrino if switch "detailed_diff_output" in namelist "neutrinoAnalysis" the is set to .true.
- eventtype=3=LowPhoto if switch "dTheta_switch" in namelist "LowElePhoto_Analysis" is set to .true.
The file shows the cross sections dsigma/dTheta versus Theta, where Theta is the polar scattering angle (in radians) of the outgoing hadron (recall here that in neutrino runs neutrinos are moving along z-direction)
Columns:
- #1: Theta [radians]
- #2: dsi/dTheta xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_XXX_dSigma_dEkin_lepton_PPP.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dEkin_lepton_PPP.dat
with:
- XXX=000, 001, 002 --- the first, second, third and so on values of the "raised variable" (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
- PPP = no_pi, p_Xn_no_pi, piplus, pi0, .... (see namelist nl_specificEvent for the full list) standing for the specific final state under consideration
PURPOSE
The file is produced in runs with:
- eventtype=5=neutrino if switch "specificEvent_Analysis" in namelist "neutrinoAnalysis" the is set to .true. and specific final states in the namelist "nl_specificEvent" are set to .true.
The file shows the cross sections for a specific final state versus kinetic energy of the outgoing lepton
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2/GeV
Columns:
- #1: kinetic energy of the outgoing lepton [GeV]
- #2: dsi/dEkin xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_XXX_dSigma_dQ2_lepton_PPP.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dQ2_lepton_PPP.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_lepton_PPP.dat but for dsigma/dQ2
The file shows the cross sections for a specific final state versus Q2
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV^2
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2/GeV^2
Columns:
- #1: Q2 [GeV^2] squared transfer momentum (Q2 = -q_mu \cdot q^\mu)
- #2: dsi/dQ2 xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_XXX_dSigma_dQ2p_lepton_PPP.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dQ2p_lepton_PPP.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_lepton_PPP.dat but for dsigma/dQ2
The file shows the cross sections for a final state with 1 mu, 0 pi , and (at least) 1 p versus Q2 here Q2 is calculated from the kinematics of the outgoing leading proton: Q2prot = (mProt - epsB)**2 - mProt**2 &
& + 2*(mNeut - epsB)*(Tp + mProt - mNeut + epsB)
This distribution can be converted into a kinetic energy distribution of 0pion events
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV^2
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2/GeV^2
Columns:
- #1: Q2 [GeV^2] squared transfer momentum, from proton kinematics
- #2: dsi/dQ2 xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_XXX_dSigma_dcos_lepton_PPP.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dcos_lepton_PPP.dat
PURPOSE
The same as diff_XXX_dSigma_dEkin_lepton_PPP.dat but for dsigma/dcos(theta_l)
The file shows the cross sections for a specific final state versus cos of the scattering angle of the outgoin lepton (with respect to neutrino direction)
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2
Columns:
- #1: cos(theta_l) cos of the scattering angle of the outgoing lepton
- #2: dsi/dcos(theta_l) xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/diff_XXX_d2Sigma_dEdcost_lepton_no_pi.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dEdcost_lepton_no_pi.dat
PURPOSE
Double differential cross section d2sigma/(dEkin dcostheta) for outgoing lepton
The file shows the cross section for the outgoing lepton versus cos of the scattering angle of the outgoin lepton (with respect to neutrino direction) and its kinetic energy
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2
Columns:
- #1: Ekin of the outgoing lepton (in GeV)
- #2: cos(theta_l) cos of the scattering angle of the outgoing lepton
- #3: dsi/dcos(theta_l) xsec
- #4: number of events contributed (only for internal use, you can safely neglect it)
- #5: xsec-statistical-error
neutrinoAnalysis/diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat
PURPOSE
Similar to diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat (see notations XXX and CHARGE there) but gives pion-nucleon invariant mass (W) distribution (true invariant mass with the sum of the 4-momenta of the outgoing particles, W^2=(p_1out+p_2out) as opposed to widely used W2=mN2+2*mN*nu-Q2 from lepton kinematics ) for events with ___1pion___ AND ___1nucleon___ in the final state
neutrinoAnalysis/diff_XXX_dSigma_dW_muon_nucleon_charge_CHARGE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dW_muon_nucleon_charge_CHARGE.dat
PURPOSE
The same as diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat but for the muon-nucleon invariant mass distribution
NOTES
Note that CHARGE is still a pion charge
neutrinoAnalysis/diff_XXX_dSigma_dW_muon_pion_charge_CHARGE.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file diff_XXX_dSigma_dW_muon_pion_charge_CHARGE.dat
PURPOSE
The same as diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat but for the muon-pion invariant mass distribution
NOTES
Note that CHARGE is still a pion charge
neutrinoAnalysis/FinalEvents.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file FinalEvents.dat
PURPOSE
The file contains the positions and four-vectors of all outgoing particles from an event For a cross-section construction all of the events have to be weighted with the 'perweight' values in col. 5 The output of this file is switched on by the switch 'outputEvents'.
Columns:
- #1: run number (from 1 to num_runs_SameEnergy)
- #2: event number (from 1 to ... less then target_A*numEnsembles)
- #3: ID of the outgoing particle
- #4: charge of the outgoing particle
- #5: perweight of the event
- #6-#8: position of the outgoing particle (set to 0 for outgoing lepton)
- #9-#12: 4-momentum of the outgoing particle in GeV
- #13: history of the outgoing particle (see history.f90 for the definition)
- #14: production_ID (type of the first event: 1=QE, 2=Delta, 34=DIS)
- #15: incoming neutrino energy
NOTES
There is always in each event a particle with weight 0. This is the nucleon on which the initial interaction happened.
In the case of an initial 2p2h process the second initial-state nucleon is not written out. It is chosen to be at the same place as the first initial nucleon (the one with weight=0), with a randomly chosen momentum in the Fermi-sea.
For large target_A, numEnsembles and num_runs_SameEnergy, this file may become very large, e.g.:
- target_A=12, 1000ens x 20runs, QE and Delta results in 41 Mb
- target_A=56, 2000ens x 5runs, QE,Delta,highRES,1piBG,DIS results in 500 Mb
neutrinoAnalysis/reconstruction_Enureal_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file reconstruction_Enureal_PPP.ZZZ.dat
PURPOSE
Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for true neutrino energy
The file shows the event distribution (flux times xsec) versus true neutrino energy for a specific final state The file is produced only for the runs with a neutrino flux (in the namelist "neutrino_induced" nuXsectionMode > 10, nuExp>0)
Columns:
- #1: true neutrino energy [GeV]
- #2: event distribution: normalized flux times xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: statistical-error of #2
neutrinoAnalysis/reconstruction_Enurec_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file reconstruction_Enurec_PPP.ZZZ.dat
with:
- PPP = no_pi, p_Xn_no_pi, piplus, pi0, .... (see namelist nl_specificEvent for the full list) standing for the specific final state under consideration
- ZZZ=000 - 008 is the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneXS.ZZZ.dat)
PURPOSE
The file shows the event distribution (normalized flux times xsec) versus reconstructed neutrino energy (see arXiv:1208.3678 [nucl-th]) for a specific final state
Reconstruction method depends on a specific process; * generally for no_pion events it is based on QE-like kinematics * and for pion events it is based on on-shell-Delta-creation assumption
The file is produced in runs with:
- eventtype=5=neutrino if switch "reconstruct_neutrino_energy" is set to .true. and switch "specificEvent_Analysis" in namelist "neutrinoAnalysis" is set to .true. and specific final states in the namelist "nl_specificEvent" are set to .true.
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV
- For eventtype=5 and process_ID=EM: nanobarns=10^{-33}cm^2/GeV
- All x-sec per particle (1/A)
Columns:
- #1: reconstructed neutrino energy [GeV]
- #2: event distribution: flux times xsec)
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: statistical-error of #2
neutrinoAnalysis/reconstruction_Enu_rec_versus_real_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file reconstruction_Enu_rec_versus_real_PPP.ZZZ.dat
PURPOSE
Similar to reconstruction_Enurec_PPP.ZZZ.dat, but double-differential
The file shows the 2-D density of the flux times xsec (see arXiv:1208.3678 [nucl-th]) versus true and reconstructed neutrino energies for a specific final state The file is produced only for the runs with a neutrino flux (in the namelist "neutrino_induced" nuXsectionMode > 10, nuExp>0)
Units:
- For event_type=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV^2
- All xsec per particle (1/A)
Columns:
- #1: true neutrino energy [GeV]
- #2: reconstructed neutrino energy [GeV]
- #3: flux-folded xsec
- #4: number of events contributed (only for internal use, you can safely neglect it)
- #5: xsec-statistical-error
NOTES
up to normalization this is "migration matrix" between true and reconstructed energies
neutrinoAnalysis/reconstruction_Q2real_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file reconstruction_Q2real_PPP.ZZZ.dat
PURPOSE
Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for the true Q2
The file contains the flux averaged cross section dsigma/dQ2, i.e. integral \int \Phi(E) dsigma/dQ2 (E) dE, versus true Q2 for a specific final state
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV^2
- For eventtype=5 and process_ID=EM (one can run it, but makes no physical sense): nanobarns=10^{-33}cm^2/GeV^2
- All xsec per particle (1/A)
Columns:
- #1: true Q2 [GeV^2] (squared momentum transfer Q2 = - q_mu \cdot q^mu)
- #2: flux-averaged dsigma/dQ2true
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/reconstruction_Q2rec_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file reconstruction_Q2rec_PPP.ZZZ.dat
PURPOSE
Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for the reconstructed Q2 (see arXiv:1208.3678 [nucl-th])
The file contains the flux averaged cross section dsigma/dQ2, i.e. integral \int \Phi(E) dsigma/dQ2 (E) dE, versus reconstructed Q2 for a specific final state
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV^2
- For eventtype=5 and process_ID=EM (one can run it, but makes no physical sense): nanobarns=10^{-33}cm^2/GeV^2
- All xsec per particle (1/A)
Columns:
- #1: reconstructed Q2 [GeV^2] (squared momentum transfer Q2 = - q_mu \cdot q^mu)
- #2: flux-averaged dsigma/dQ2rec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/reconstruction_Q2_rec_versus_real_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file reconstruction_Q2_rec_versus_real_PPP.ZZZ.dat
PURPOSE
Similar to reconstruction_Enurec_PPP.ZZZ.dat, but double-differential in Q2
The file shows the dsigma/dQ2true dQ2rec distribution versus true and reconstructed Q^2 for a specific final state
Units:
- For eventtype=5 and process_ID=CC and NC: 10^{-38} cm^2/GeV^4
- For eventtype=5 and process_ID=EM (one can run it, but makes no physical sense): nanobarns=10^{-33}cm^2/GeV^4
Columns:
- #1: true Q2 [GeV^2]
- #2: reconstructed Q2 [GeV^2]
- #3: dsigma/dQ2true dQ2rec
- #4: number of events contributed (only for internal use, you can safely neglect it)
- #5: xsec-statistical-error
neutrinoAnalysis/oscillations_UUU_real_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file oscillations_UUU_real_PPP.ZZZ.dat
with:
- UUU = "mumu", "mue", "mue_max", "mue_antimax" for oscillation of a given experimental flux for muon_neutrino survival ("mumu"), mu_e appearance with delta_CP=0 ("mue"), delta_CP=pi ("mue_max"), delta_CP=-pi ("mue_antimax")
- PPP = no_pi, p_Xn_no_pi, piplus, pi0, .... (see namelist nl_specificEvent for the full list) standing for the specific final state under consideration
- ZZZ=000 - 008 is the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneXS.ZZZ.dat)
PURPOSE
Similar to reconstruction_Enureal_PPP.ZZZ.dat, but for oscillated neutrino flux
Columns:
- #1: true neutrino energy [GeV]
- #2: flux-folded xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/oscillations_UUU_rec_PPP.ZZZ.dat [ Output files ]
[ Top ] [ neutrinoAnalysis ] [ Output files ]
NAME
file oscillations_UUU_rec_PPP.ZZZ.dat
with:
- UUU = "mumu", "mue", "mue_max", "mue_antimax" for oscillation of a given experimental flux for muon_neutrino survival ("mumu"), mu_e appearance with delta_CP=0 ("mue"), delta_CP=pi ("mue_max"), delta_CP=-pi ("mue_antimax")
- PPP = no_pi, p_Xn_no_pi, piplus, pi0, .... (see namelist nl_specificEvent for the full list) standing for the specific final state under consideration
- ZZZ=000 - 008 is the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneXS.ZZZ.dat)
PURPOSE
Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for oscillated neutrino flux
Columns:
- #1: reconstructed neutrino energy [GeV]
- #2: event rate: flux times xsec
- #3: number of events contributed (only for internal use, you can safely neglect it)
- #4: xsec-statistical-error
neutrinoAnalysis/IsBound [ Functions ]
[ Top ] [ neutrinoAnalysis ] [ Functions ]
NAME
logical function IsBound(part)
PURPOSE
return whether particle is bound or not
neutrinoAnalysis/IsBelowThreshold [ Functions ]
[ Top ] [ neutrinoAnalysis ] [ Functions ]
NAME
logical function IsBelowThreshold(part)
PURPOSE
Returns .true. when a particle is below a given detection threshold This routine can be used to remove all events with kinetic energies and angles below specified detection thresholds, if an event is below threshold it is not added to the full event in subroutine neutrino_Analyze
neutrinoAnalysis/oscillationProbability [ Subroutines ]
[ Top ] [ neutrinoAnalysis ] [ Subroutines ]
NAME
subroutine oscillationProbability(Enu,L,deltaCP,Posc_mumu,Posc_mue,Posc_mue_max, &
& Posc_mue_antimax)
PURPOSE
Calculate Oscillation Probability
neutrinoAnalysis/PrintVals10 [ Subroutines ]
[ Top ] [ neutrinoAnalysis ] [ Subroutines ]
NAME
PURPOSE
abbreviation used for many outputfiles: printing the vals stored in sigma, including its error
neutrinoAnalysis/PrintHeader10 [ Subroutines ]
[ Top ] [ neutrinoAnalysis ] [ Subroutines ]
NAME
PURPOSE
abbreviation used for many outputfiles: printing a header line