fRE Function

  1. korc_experimental_pdf.f90
  2. korc_experimental_pdf
  3. fRE

private function fRE(eta, p)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in) :: eta
real(kind=rp), intent(in) :: p

Return Value real(kind=rp)


Calls

proc~~fre~4~~CallsGraph proc~fre~4 fRE proc~fgamma~2 fGamma proc~fre~4->proc~fgamma~2 proc~deg2rad~5 deg2rad proc~fre~4->proc~deg2rad~5

Called by

proc~~fre~4~~CalledByGraph proc~fre~4 fRE proc~frexpr fRExPR proc~frexpr->proc~fre~4 proc~sample_distribution~2 sample_distribution proc~sample_distribution~2->proc~fre~4 proc~get_experimentalg_distribution get_experimentalG_distribution proc~get_experimentalg_distribution->proc~sample_distribution~2

Contents

Source Code

fRE

Source Code

  FUNCTION fRE(eta,p)
    REAL(rp), INTENT(IN) 	:: eta ! pitch angle in degrees
    REAL(rp), INTENT(IN) 	:: p ! momentum in units of mc
    REAL(rp) 				:: fRE
    REAL(rp) 				:: fE
    REAL(rp) 				:: feta
    REAL(rp) 				:: A
    REAL(rp) 				:: Eo

    Eo = SQRT(p**2.0_rp + 1.0_rp)

    A = (2.0_rp*pdf_params%E/(pdf_params%Zeff + 1.0_rp))*(p**2/SQRT(p**2.0_rp + 1.0_rp))
    A = A*pdf_params%A_fact
    feta = 0.5_rp*A*EXP(A*COS(deg2rad(eta)))/SINH(A)

    fE = fGamma(Eo,pdf_params%k,pdf_params%t/xo)/pdf_params%fGo

    fRE = fE*feta
  END FUNCTION fRE