TABLE OF CONTENTS
- 1. /hadronTensor_ResProd
- 1.1. hadronTensor_ResProd/hadronTensor_R
- 1.1.1. hadronTensor_R/writeError
- 1.2. hadronTensor_ResProd/hadronTensor_1_2
- 1.2.1. hadronTensor_1_2/j_pos
- 1.2.2. hadronTensor_1_2/j_neg
- 1.2.3. hadronTensor_1_2/j_dagger
- 1.3. hadronTensor_ResProd/hadronTensor_3_2
- 1.3.1. hadronTensor_3_2/j_pos
- 1.3.2. hadronTensor_3_2/j_neg
- 1.3.3. hadronTensor_3_2/j_dagger
- 1.4. hadronTensor_ResProd/parityError
/hadronTensor_ResProd [ Modules ]
NAME
module hadronTensor_ResProd
PURPOSE
This module defines the hadronic tensors for resonance production by EM- and weak currents
hadronTensor_ResProd/hadronTensor_R [ Functions ]
[ Top ] [ hadronTensor_ResProd ] [ Functions ]
NAME
function hadronTensor_R(pi,pf,resID,targetCharge,process,matrix) result(success)
PURPOSE
This function returns the hadronic tensor for resonance production of the resonance "resID" by a current type defined by "process" (e.g. Gamma* N -> R).
INPUTS
- real, dimension(0:3),intent(in) :: pi -- 4-momentum of incoming nucleon
- real, dimension(0:3),intent(in) :: pf -- 4-momentum of outgoing resonance
- integer,intent(in) :: resID -- ID of resonance
- integer,intent(in) :: process -- type of the process (EM, CC, NC)
- integer,intent(in) :: targetCharge -- Charge of target nucleon
OUTPUT
- complex, dimension(0:3,0:3) :: matrix -- H^mu nu
- logical, intent(out) :: success -- true if resonance has coupling strength to the process;
- false if there is no coupling, therefore H^mu nu=0
NOTES
- We use the Peskin notation where the spin sum(u(p) ubar(p))=slashed(p)+m !!
hadronTensor_R/writeError [ Subroutines ]
[ Top ] [ hadronTensor_R ] [ Subroutines ]
PURPOSE
Error Message
hadronTensor_ResProd/hadronTensor_1_2 [ Functions ]
[ Top ] [ hadronTensor_ResProd ] [ Functions ]
NAME
function hadronTensor_1_2(pi,pf,parity,G,GA,resID) result(matrix)
PURPOSE
This function returns the hadronic tensor for spin=1/2 resonances.
INPUTS
- real, dimension(1:2),intent(in) :: F -- Vector or EM-Form factors
- real, dimension(1:2),intent(in) :: FA -- Axial form factors
- real, dimension(0:3),intent(in) :: pi -- 4-momentum of incoming nucleon
- real, dimension(0:3),intent(in) :: pf -- 4-momentum of outgoing resonance
- integer,intent(in) :: resID -- ID of resonance
- integer,intent(in) :: parity -- Parity of resonance
- integer,intent(in) :: process -- type of the process (EM, CC, NC)
OUTPUT
hadronTensor_1_2/j_pos [ Functions ]
[ Top ] [ hadronTensor_1_2 ] [ Functions ]
NAME
function j_pos(mu) result(matrix)
PURPOSE
This function the hadronic flux "J^mu" for resonances of positive parity
INPUTS
integer, intent(in) :: mu
OUTPUT
- complex, dimension(0:3,0:3) :: matrix
hadronTensor_1_2/j_neg [ Functions ]
[ Top ] [ hadronTensor_1_2 ] [ Functions ]
NAME
function j_neg(mu) result(matrix)
PURPOSE
This function the hadronic flux "J^mu" for resonances of negative parity
INPUTS
integer, intent(in) :: mu
OUTPUT
- complex, dimension(0:3,0:3) :: matrix
hadronTensor_1_2/j_dagger [ Functions ]
[ Top ] [ hadronTensor_1_2 ] [ Functions ]
NAME
function j_neg_dagger(mu) result(matrix)
PURPOSE
This function returns "gamma_0 (J^mu)^dagger gamma_0" for resonances of negative parity
INPUTS
integer, intent(in) :: mu
OUTPUT
- complex, dimension(0:3,0:3) :: matrix
hadronTensor_ResProd/hadronTensor_3_2 [ Functions ]
[ Top ] [ hadronTensor_ResProd ] [ Functions ]
NAME
function hadronTensor_3_2(pi,pf,parity,formfactor,resID) result(matrix)
PURPOSE
This function returns the hadronic tensor for spin=3/2 resonances.
INPUTS
- real, dimension(0:3),intent(in) :: pi -- 4-momentum of incoming nucleon
- real, dimension(0:3),intent(in) :: pf -- 4-momentum of outgoing resonance
- integer,intent(in) :: parity -- Parity of resonance
- real, dimension(1:8),intent(in) :: FormFactor -- 1-4: Vector or EM, 5-8: Axial form factors
- integer,intent(in) :: resID -- ID of resonance
- integer,intent(in) :: process -- type of the process (EM, CC, NC)
OUTPUT
hadronTensor_3_2/j_pos [ Functions ]
[ Top ] [ hadronTensor_3_2 ] [ Functions ]
NAME
function j_pos(lambda,nu) result(matrix)
PURPOSE
- This function evaluates the hadronic flux "J^lambda nu" for resonances of positive parity
- Same as negative parity case times gamma(5).
INPUTS
- integer, intent(in) :: lambda,nu -- lambda is the lorentz index which is contracted with the 3/2 spinor
OUTPUT
- complex, dimension(0:3,0:3) :: matrix
hadronTensor_3_2/j_neg [ Functions ]
[ Top ] [ hadronTensor_3_2 ] [ Functions ]
NAME
function j_neg(alpha,nu) result(matrix)
PURPOSE
- This function evaluates the hadronic flux "J^lambda nu" for resonances of negative parity
- See Phys Rev D74, 014009 (2006), eq. 4.3
INPUTS
- integer, intent(in) :: alpha,nu -- lambda is the lorentz index which is contracted with the 3/2 spinor
OUTPUT
- complex, dimension(0:3,0:3) :: matrix
hadronTensor_3_2/j_dagger [ Functions ]
[ Top ] [ hadronTensor_3_2 ] [ Functions ]
NAME
function j_dagger(alpha,nu) result(matrix)
PURPOSE
- This function evaluates "gamma_o (J^lambda nu)^dagger gamm_o" for resonances of negative parity
- See Phys Rev D74, 014009 (2006), eq. 4.3
INPUTS
- integer, intent(in) :: alpha,nu -- lambda is the lorentz index which is contracted with the 3/2 spinor
OUTPUT
- complex, dimension(0:3,0:3) :: matrix
hadronTensor_ResProd/parityError [ Subroutines ]
[ Top ] [ hadronTensor_ResProd ] [ Subroutines ]
PURPOSE
Error Message