TEvtGen/Tauola/tauola-fortran/new-currents/RChL-currents/rcht_3pi/funct_3pi.f

   1 C **********************************************************
   2 C library of functions used in calculation of currents
   3 C references:
   4 C [1]   arXiv:0911.4436 (hep-ph)  D. Gomez Dumm et al. (tau -> 3pi nu)
   5 C [2]   arXiv:0911.2640 (hep-ph) D. Gomezz Dumm et al. (tau -> KKpi nu)
   6 C [3]   arXiv:0807.4883 (hep-ph) D. R. Boito et al. (tau -> Kpi nu)
   7 C [4]   P. Roig, talk at (tau -> 2 pi nu)
   8 C [5]   arXiv:0803.2039  (hep-ph) E. Arganda et al., Appendix B (tau -> KK nu)
   9 C **********************************************************
  10
  11       FUNCTION ALP1_RPT(Q,S,T,M1SQ,M2SQ,M3SQ,
  12      $                           MMM,MMM2,GGG)
  13       IMPLICIT NONE
  14       REAL                      Q,S,T
  15       DOUBLE PRECISION                M1SQ,M2SQ,M3SQ,
  16      &                           MMM,MMM2,GGG
  17 C **********************************************************
  18 C     COMPLEX FUNCTION ALP1_RPT;
  19 C     one-resonance contribution to
  20 C     F3pi_rcht(iform=1,2) at f3pi_rcht.f;
  21 C     corrected  formula (8) of REF [1] including rho1,
  22 C     the factor sqrt(2)*FV_RPT*GV_RPT/(3.*FPI_RPT**3) is included in
  23 C     F3pi_rcht(iform1=1,2) at f3pi_rcht.f
  24 C **********************************************************
  25       REAL U
  26       include '../parameter.inc'
  27       include '../funct_declar.inc'
  28
  29        U= Q-S-T+M1SQ+M2SQ+M3SQ
  30        ALP1_RPT = - 3.*S/(1.+BETA_RHO)*
  31      $     (1./(S-MMM**2+i*MMM*GRHO_RCHT(S,MMM))+
  32      $      BETA_RHO/(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG)))
  33      $      +(2.*GV_RPT/FV_RPT-1.)*
  34      $      ((2.*Q-2.*S-U)/(1.+BETA_RHO)*
  35      $      (1./(S-MMM**2+i*MMM*GRHO_RCHT(S,MMM))+
  36      $      BETA_RHO/(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG)))+
  37      $      (U-S)/(1.+BETA_RHO)*
  38      $      (1./(T-MMM**2+i*MMM*GRHO_RCHT(T,MMM))+
  39      $       BETA_RHO/(T-MMM2**2+i*MMM2*GRHO1_RCHT(T,MMM2,GGG))))
  40
  41       RETURN
  42
  43       END
  44
  45
  46
  47       FUNCTION BETA_RPT(Q,S,T,M1SQ,M2SQ,M3SQ,MMM,MMM2,GGG)
  48       IMPLICIT NONE
  49       REAL                      Q,S,T
  50       DOUBLE PRECISION                M1SQ,M2SQ,M3SQ,MMM,MMM2,GGG
  51 C **********************************************************
  52 C     COMPLEX FUNCTION BETA_RPT ; two-resonance contribution to
  53 C     F3pi_rcht(iform=1,2) at f3pi_rcht.f;
  54 C     corrected  formula (8) of REF [1] including rho1,
  55 C     the factor 4*FA_RPT*GV_RPT/(3.*FPI_RPT**3)*QQ/D_a1(QQ)
  56 C     is included in F3pi_rcht(iform1=1,2) at f3pi_rcht.f
  57 C **********************************************************
  58       include '../parameter.inc'
  59       include '../funct_declar.inc'
  60
  61       REAL U,LAM0_RPT,LAM1_RPT,LAM2_RPT,FF1_RPT,FF2_RPT,FF_REL
  62
  63       FF_REL = FPI_RPT*FPI_RPT/(FV_RPT*FV_RPT)
  64
  65        LAM1_RPT = FPI_RPT*FPI_RPT/(2.D0*SQRT(2.D0)*FA_RPT*GV_RPT)
  66        LAM2_RPT = -(1.-2.*FF_REL)*LAM1_RPT
  67        LAM0_RPT = (LAM1_RPT + LAM2_RPT)/4.
  68
  69         U= Q-S-T+M1SQ+M2SQ+M3SQ
  70
  71       FF1_RPT = -LAM0_RPT*M1SQ/Q +LAM1_RPT*S/Q+LAM2_RPT
  72       FF2_RPT = -LAM0_RPT*M1SQ/Q +LAM1_RPT*T/Q+LAM2_RPT
  73
  74       BETA_RPT = -3.*(LAM1_RPT+LAM2_RPT)*S/(1.+BETA_RHO)*
  75      $     (1./(S-MMM**2+i*MMM*GRHO_RCHT(S,MMM))+
  76      $      BETA_RHO/(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG)))+
  77      $         FF1_RPT*(2.*Q+S-U)/(1.+BETA_RHO)*
  78      $     (1./(S-MMM**2+i*MMM*GRHO_RCHT(S,MMM))+
  79      $      BETA_RHO/(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG)))+
  80      $         FF2_RPT*(U-S)/(1.+BETA_RHO)*
  81      $     (1./(T-MMM**2+i*MMM*GRHO_RCHT(T,MMM))+
  82      $      BETA_RHO/(T-MMM2**2+i*MMM2*GRHO1_RCHT(T,MMM2,GGG)))
  83
  84       RETURN
  85
  86       END
  87
  88       FUNCTION ALP1_RPT_RHO1(Q,S,T,M1SQ,M2SQ,M3SQ,
  89      $                           MMM,MMM2,GGG)
  90       IMPLICIT NONE
  91       REAL                      Q,S,T
  92       DOUBLE PRECISION                M1SQ,M2SQ,M3SQ,
  93      &                           MMM,MMM2,GGG
  94 C **********************************************************
  95 C     COMPLEX FUNCTION ALP1_RPT_RHO1;
  96 C     one-resonance contribution to
  97 C     F3pi_rcht(iform=1,2) at f3pi_rcht.f;
  98 C     corrected  formula (8) of REF [1] including rho1,
  99 C     the factor sqrt(2)*FV_RPT*GV_RPT/(3.*FPI_RPT**3) is included in
 100 C     F3pi_rcht(iform1=1,2) at f3pi_rcht.f
 101 C **********************************************************
 102       REAL U
 103       include '../parameter.inc'
 104       include '../funct_declar.inc'
 105
 106        U= Q-S-T+M1SQ+M2SQ+M3SQ
 107        ALP1_RPT_RHO1 = - 3.*S
 108      $      /(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG))
 109      $      +(2.*GV1_RPT/FV1_RPT-1.)*
 110      $      (  (2.*Q-2.*S-U)
 111      $      /(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG))+
 112      $      (U-S)/(T-MMM2**2+i*MMM2*GRHO1_RCHT(T,MMM2,GGG))  )
 113
 114       RETURN
 115
 116       END
 117
 118
 119
 120       FUNCTION BETA_RPT_RHO1(Q,S,T,M1SQ,M2SQ,M3SQ,MMM,MMM2,GGG)
 121       IMPLICIT NONE
 122       REAL                      Q,S,T
 123       DOUBLE PRECISION                M1SQ,M2SQ,M3SQ,MMM,MMM2,GGG
 124 C **********************************************************
 125 C     COMPLEX FUNCTION BETA_RPT_RHO1 ; two-resonance contribution to
 126 C     F3pi_rcht(iform=1,2) at f3pi_rcht.f;
 127 C     corrected  formula (8) of REF [1] including rho1,
 128 C     the factor 4*FA_RPT*GV_RPT/(3.*FPI_RPT**3)*QQ/D_a1(QQ)
 129 C     is included in F3pi_rcht(iform1=1,2) at f3pi_rcht.f
 130 C **********************************************************
 131       include '../parameter.inc'
 132       include '../funct_declar.inc'
 133
 134       REAL U,LAM0_RPT,LAM1_RPT,LAM2_RPT,FF1_RPT,FF2_RPT,FF_REL
 135
 136       FF_REL = FPI_RPT*FPI_RPT/(FV1_RPT*FV1_RPT)
 137
 138        LAM1_RPT = FPI_RPT*FPI_RPT/(2.D0*SQRT(2.D0)*FA_RPT*GV1_RPT)
 139        LAM2_RPT = -(1.-2.*FF_REL)*LAM1_RPT
 140        LAM0_RPT = (LAM1_RPT + LAM2_RPT)/4.
 141
 142         U= Q-S-T+M1SQ+M2SQ+M3SQ
 143
 144       FF1_RPT = -LAM0_RPT*M1SQ/Q +LAM1_RPT*S/Q+LAM2_RPT
 145       FF2_RPT = -LAM0_RPT*M1SQ/Q +LAM1_RPT*T/Q+LAM2_RPT
 146
 147       BETA_RPT_RHO1 = -3.*(LAM1_RPT+LAM2_RPT)*S
 148      $     /(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG)) +
 149      $         FF1_RPT*(2.*Q+S-U)
 150      $     /(S-MMM2**2+i*MMM2*GRHO1_RCHT(S,MMM2,GGG))+
 151      $         FF2_RPT*(U-S)
 152      $     /(T-MMM2**2+i*MMM2*GRHO1_RCHT(T,MMM2,GGG))
 153
 154       RETURN
 155
 156       END
 157
 158
 159
 160