gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/Coll_nuN [ Modules ]

[ Top ] [ 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

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