TABLE OF CONTENTS
- 1. /Coll_nuN
/Coll_nuN [ Modules ]
NAME
module Coll_nuN
PURPOSE
neutrino Nucleon scatterings with PYTHIA
Coll_nuN/DoColl_nuN_Py [ Subroutines ]
[ Top ] [ Coll_nuN ] [ Subroutines ]
NAME
subroutine DoColl_nuN_Py(eNev,outPart,flagOK, doMassless,cross)
PURPOSE
generate a high energy neutrino event with PYTHIA
returned cross section is dsigma/dE'dcost in mb/GeV
Here the outgoing kinematics is fixed via PYTHIA-Option MSTP(23)=1
INPUTS
- type(electronNucleon_event) :: eNev -- kinematics
- logical :: doMassless -- neglect masses of quarks and diquarks
OUTPUT
- type(particle),dimension(:) :: outPart -- outgoing particles
- logical :: flagOK -- .TRUE., if everything okay
- real :: cross -- cross section
NOTES
With the PYTHIA-Option MSTP(23)=1, the routine PYREMN tries to keep the chosen x and Q2 values of the partonic process also for the outer hadronic process.
Coll_nuN/DoColl_nuN_Py_Free [ Subroutines ]
[ Top ] [ Coll_nuN ] [ Subroutines ]
NAME
subroutine DoColl_nuN_Py_Free(eNev,outPart,flagOK, doMassless,cross)
PURPOSE
generate a high energy neutrino event with PYTHIA
Here the outgoing kinematics is free at input and is set in this routine
INPUTS
- type(electronNucleon_event) :: eNev -- kinematics
- logical :: doMassless -- neglect masses of quarks and diquarks
OUTPUT
- type(electronNucleon_event) :: eNev -- kinematics
- type(particle),dimension(:) :: outPart -- outgoing particles
- logical :: flagOK -- .TRUE., if everything okay
- real :: cross -- cross section
NOTES
We give as input the neutrino and the nucleon in medium and free nucleon kinematics. Since we can not ask for any cuts on the outgoing lepton/neutrino, we set the momentum of this particle as output
Coll_nuN/AnaEstimate [ Functions ]
[ Top ] [ Coll_nuN ] [ Functions ]
NAME
real function CalcAnaEstimate(eNev)
PURPOSE
Calculate the cross sections according eqs. (C.37)-(C.39) of Phys.Rep.
OUTPUT
- function value = dsigma/dE'dcost in mb/GeV
Coll_nuN/AnaEstimatePythia [ Functions ]
[ Top ] [ Coll_nuN ] [ Functions ]
NAME
real function CalcAnaEstimatePythia(eNev)
PURPOSE
Calculate the cross sections according Pythia for the em case
OUTPUT
- function value = dsigma/dE'dcost in mb/GeV