[u/mrichter/AliRoot.git] / TEvtGen / Tauola / tauola-fortran / new-currents / RChL-currents / rcht_common / value_parameter.f

      SUBROUTINE CH3PISET(JJ)
C information on 3 pion sub-channel under construction obtained
C J3PI=1 means 3 prong 
C J3PI=2 means 1 prong
C to be initialized in routine DPHSAA of tauola.f
      COMMON /CHANOPT/ J3PI
      INTEGER          J3PI
      J3PI=JJ
      end
     
      SUBROUTINE CH3PIGET(JJ)
C information on 3 pion sub-channel under construction obtained
C J3PI=1 means 3 prong 
C J3PI=2 means 1 prong
C to be initialized in routine DPHSAA of tauola.f
      COMMON /CHANOPT/ J3PI
      INTEGER          J3PI
      IF (J3PI.EQ.1.OR.J3PI.EQ.2) THEN
        JJ=J3PI
      ELSE
       write(*,*) 'FROM value_parameter.f CH3PIGET, wrong J3PI=',J3PI
       stop
      ENDIF
      end

      SUBROUTINE GETFF2PIRHO(JJ)
      IMPLICIT NONE
      include '../parameter.inc'
      INTEGER JJ
      JJ = FF2PIRHO
      END

      SUBROUTINE OLACHNL(SIGN)
C provides sign of tau, to be used in CP dependent parts of current.
      COMMON / JAKI   /  JAK1,JAK2,JAKP,JAKM,KTOM
      INTEGER            JAK1,JAK2,JAKP,JAKM,KTOM
      COMMON / IDFC  / IDFF
      INTEGER KTO
      REAL    SIGN
      IF     (KTOM.EQ.1.OR.KTOM.EQ.-1) THEN
        SIGN= IDFF/ABS(IDFF)
      ELSEIF (KTOM.EQ.2) THEN
        SIGN=-IDFF/ABS(IDFF)
      ELSE
        PRINT *, 'STOP IN OLACHNL: KTOM=',KTOM
        STOP
      ENDIF
      END

      FUNCTION COEFrr(I,J)
C clebsh gordan (or so ...)  coefs for 3 scalar final states
      implicit none
C TAUOLA RChL COEF(I,J) =  COEFr(I,J)
      REAL COEFr(1:5,0:7)
      REAL COEFrr
      DATA PI /3.141592653589793238462643/
      REAL PI
      DATA ICONT /0/
      INTEGER ICONT
      INTEGER I,J
      REAL FPIr

C initialization of FPI matrix defined in ...
C FPIc is to be used with cleo initialization
C FPIr is to be used with RChL initialization
C actual choice is made in ???


      DATA  FPIr /92.4E-3/


C initialization of COEF matrix defined in ...
C COEFc is to be used with cleo initialization
C COEFr is to be used with RChL initialization
      IF (ICONT.EQ.0) THEN
       ICONT=1
C
C********* COEFr(I,J) *******

       COEFr(1,0)= 1.
       COEFr(2,0)= -1.
       COEFr(3,0)= 0.
       COEFr(4,0)= 1.
       COEFr(5,0)= 0.

       COEFr(1,1)= 1.
       COEFr(2,1)= -1.
       COEFr(3,1)= 0.
       COEFr(4,1)= 1.
       COEFr(5,1)= 1.
C
       COEFr(1,2)=1.
       COEFr(2,2)= -1.
       COEFr(3,2)= 0.0
       COEFr(4,2)= 1.
       COEFr(5,2)=1.
C
       COEFr(1,3)= 0.
       COEFr(2,3)= 1.
       COEFr(3,3)= -1.
       COEFr(4,3)= 1.
       COEFr(5,3)= - 1.
C
       COEFr(1,4)= 1.0/SQRT(2.)/3.0
       COEFr(2,4)=-1.0/SQRT(2.)/3.0
       COEFr(3,4)= 0.0
       COEFr(4,4)= 0.0
       COEFr(5,4)= 0.0
C
       COEFr(1,5)=-SQRT(2.)/3.0
       COEFr(2,5)= SQRT(2.)/3.0
       COEFr(3,5)= 0.0
       COEFr(4,5)= 0.0
       COEFr(5,5)=-SQRT(2.)
C
       COEFr(1,6)= 1./3.
       COEFr(2,6)=-2./3.
       COEFr(3,6)= 2./3.
       COEFr(4,6)= 0.0
       COEFr(5,6)=-2.0
C
       COEFr(1,7)= 0.0
       COEFr(2,7)= 0.0
       COEFr(3,7)= 0.0
       COEFr(4,7)= 0.0
       COEFr(5,7)=-SQRT(2.0/3.0)
      ENDIF

      COEFrr=COEFr(I,J)
      END

      subroutine rchl_parameters(KAK)
      implicit none
C==============================================================================
C  Initialization, of '../parameter.inc' common block group 
C 
C  KAK may be equal to JAK of TAUOLA namespace, but it is not always the case 
C  Hard-coded  fit parameters:
C  rho, rhoprime, f2(1275), f0(1186), sigma(made up!)
C  The value of both the mass and width of resonances are taken 
C  from  fit to ALEPH data  (ref [1], Set 1) 
C  References: [1] arXiv: 0911.4436 [hep-ph] D.  Gomez Dumm et al
C              [2] arXiv: 0911.2640 [hep-ph] D.  Gomez Dumm et al. 
C              [3] P Roig, talk PhiPsi2011, Novosibirsk
C              [4] arXiv:0807.4883 [hep-ph] Diogo R. Boito et al.
C              [5] arXiv:0803.1786 [hep-ph] M. Jamin et al.
C  WARNING:    some of parameters require RERUN of da1wid_tot_rho1_gauss.f
C              pretabulating Q dependent a1 width,
C              directory RChL-currents/tabler/a1
C==============================================================================
      include '../parameter.inc'
      INTEGER KAK
      DATA IWARM/0/
      INTEGER IWARM
      INTEGER          J3PI
      COMMON /CHANOPT/ J3PI

         IF(KAK.EQ.4) THEN  
C  /MASS_RES/; resonances parameters initialization: 
C                           ! at present only for two pion mode non-default 
c                           ! values are used:
      mro = 0.77554d0
      mrho1 = 1.453d0
      grho1 = 0.50155D0
c   /PAR_RHOPRIME/; parameters of rho' and rho'' 
C                   used for 2 pion form factor, reference [3]
      COEF_GA =  0.14199D0    
      COEF_DE = -0.12623D0    
      phi_1   = -0.17377D0 
      phi_2   =  0.27632D0  
      grho2   =  0.41786D0  
      mrho2   =  1.8105d0
         ELSE IF(KAK.EQ.5) THEN 
      MRO   = 0.771849d0      !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MRHO1 = 1.35d0     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      GRHO1 = 0.448379d0 !0.473287d0       !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f

      ELSE
      MRO   = 0.775     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MRHO1 = 1.465     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      GRHO1 = 0.4       !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f

c   /PAR_RHOPRIME/; parameters of rho' and rho'' 
C                   used for 2 kaon form factor, reference [3]
c                   FOR THE MOMENT THEIR NUMERICAL VALUES COINCIDE WITH
c                   ONES FOR THE TWO PION MODE !!!!
      COEF_GA =  0.14199D0    
      COEF_DE = -0.12623D0    
      phi_1   = -0.17377D0 
      phi_2   =  0.27632D0  
      grho2   =  0.41786D0  
      mrho2   =  1.8105d0
         ENDIF


        IF(KAK.EQ.70)  THEN    ! non default values to be used
                               ! for KPI MODE NO FSR INTERACTION
c   /PAR_KPI/; parameters for Kpi mode, reference [4], table 4, row2
      MKST        = 0.943d0
      MKSTPR      = 1.374D0
      GAMMA_KST   = 0.06672d0  
      GAMMA_KSTPR = 0.240d0
      GAMMA_RCHT  =-0.039d0
        ELSE IF(KAK.EQ.71)  THEN     ! non default values to be used 
                                     ! for KPI MODE WITH FSR INTERACTION
c   parameters for Kpi mode, reference [5]
      MKST = 0.8953d0
      GAMMA_KST = 0.0475d0
      MKSTPR = 1.307d0
      GAMMA_KSTPR = 0.206d0
      GAMMA_RCHT = -0.043d0
        ELSE
C   /MASS_SCAL/; stable particles - final scalars 
      Mksp  = 0.89166d0   !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      Mks0  = 0.89610d0   !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MKST = (Mksp +Mks0)/2.
      MKSTPR = 1.374d0    !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      GAMMA_KST   = 0.06672  
      GAMMA_KSTPR = 0.240
      GAMMA_RCHT  = -0.043 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
        ENDIF

C  /RCHT_3PI/; model parameters; their value are from fit,
c             reference [1], set 1
C              CHANGE OF THEIR VALUES REQUIRES 
C              RERUN /tabler/a1/da1wid_tot_rho1_gauss.f
 
      IF(KAK.EQ.5) THEN
      FPI_RPT  = 0.091337d0
      FV_RPT   = 0.168652d0   
      FA_RPT   = 0.131425d0                 
      BETA_RHO = -0.318551d0 
      ELSE
      FPI_RPT  = 0.0924
      FV_RPT   = 0.18  
      FA_RPT   = 0.149                
      BETA_RHO =  -0.25  
      ENDIF
 
      FK_RPT = FPI_RPT*1.198d0  
      GV_RPT   = FPI_RPT*FPI_RPT/FV_RPT

c$$$c   It has to be used for a new parametrization of rho1 for 3pions, 
C$$$c   that is not checked yet
c$$$c        IF(KAK.EQ.5)  THEN  ! high energy behaviour imposes these relations   
c$$$c      GV_RPT   = 0.066  
c$$$c      FV1_RPT = 0.18D0
c$$$c      GV1_RPT = (FPI_RPT*FPI_RPT- FV_RPT*GV_RPT)/FV1_RPT
c$$$c        ELSE
c$$$c      GV_RPT   = FPI_RPT*FPI_RPT/FV_RPT
c$$$c      ENDIF

c   /SCAL_3PI/; parameters of sigma meson for 3 pion mode
C* Parameteres for the sigma contribution, using BW for sigma
      IF(KAK.EQ.5)  THEN
          IF (J3PI.EQ.1) THEN
      alpsig = -8.795938d0 
      betasig = 9.763701d0 
      gamsig =  1.264263d0 
      delsig =  0.656762d0 
      rsigma =  1.866913d0 
           ELSE IF (J3PI.EQ.2) THEN
      alpsig = 1.139486d0*0.63d0
      betasig = 1.139486d0*0.63d0
      gamsig =  0.889769d0*0.63d0
      delsig =  0.889769d0*0.63d0 
      rsigma =  0.000013d0
           ENDIF
       ENDIF

C  /MASS_RES/ 
      IF(KAK.EQ.5)  THEN
      MMA1  = 1.091865d0      !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
          IF (J3PI.EQ.1) THEN
      MSIG  = 0.487512d0 
      GSIG  = 0.70d0     
          ELSE IF(J3PI.EQ.2) THEN
      MSIG  = 0.55d0 
      GSIG  = 0.7d0 
          ENDIF
      ELSE
      MMA1 = 1.12
      MSIG = 0.475
      GSIG = 0.550
      ENDIF
        call rchl_REparam(0,IWARM,KAK)     
         IF (IWARM.EQ.1) RETURN  ! parameters below do not need 
         IWARM=1                 ! re-initialization

C  /MASS_RES/ 
      GRO   = 0.149d0     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MF2   = 1.275d0
      GF2   = 0.185d0
      MF0   = 1.186d0
      GF0   = 0.350d0
      MSG   = 0.860d0
      GSG   = 0.880d0
      MPHI  = 1.019d0     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      GPHI  = 0.0042d0    !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MOM   = 0.781940d0  !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      GOM   = 0.00843d0   !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f

C    /RES_MIXING_RCHT/; a parameter defines w-phi angle mixing
      THETA = 35.*PI/180. 

C    /FF0SCKPI/ a parameter normalized FFSC_KPI 
      F00 = 0.972


C   /MASS_SCAL/; stable particles - final scalars 
C              CHANGE OF THEIR VALUES (useful for some tests) REQUIRES, 
C              RERUN /tabler/a1/da1wid_tot_rho1_gauss.f

      MPIZ    = 0.1349766d0           !PKORB(1,7)    ! NEUTRAL PION MASS
      MPIC    = 0.13957018d0          !PKORB(1,8)    ! CHARGED PION MASS
      MMPI_AV = (MPIZ+2.*MPIC)/3.d0
      MKZ     = 0.497648d0            !PKORB(1,12)    ! NEUTRAL KAON MASS
      MKC     = 0.493677d0            !PKORB(1,11)    ! CHARGED KAON MASS
      MMK     = (MKC+MKZ)/2.d0
      MTAU    = 1.777
      MNUTA   = 0.001
      META    = 0.547d0
 

c   /PAR_KKPI/; parameters to describe KKpi modes, reference [2]
C              CHANGE OF THEIR VALUES REQUIRES 
C              RERUN /tabler/a1/da1wid_tot_rho1_gauss.f

      G2    =  mro/(192.*pi*pi*sqrt(2.)*FV_RPT)*3.
      G13   = -2.*g2 
      G4    = -0.72
      G5    = -0.6-2.*g4 
      C125  =  0.
      C1256 = -3/96./pi**2*FV_RPT*MRO/SQRT(2.)/FPI_RPT**2 
      C1235 =  0.
      C4    = -0.07
      D123  =  0.05 
      D3    = -MRO**2/(64.*PI*PI*FPI_RPT**2) 


c   /PAR_KPI/; parameters to describe Kpi mode, reference [4]
      Ht0 = -1.2400398216503017D-2
C     Ht0 = !!!!! TO ADD A FORMULAE FOR Ht0 (Jamin's email) !!!!!
      lap_KPI = 24.66e-3
      lapp_KPI = 11.99e-4
      c1_KPI = lap_KPI/mpic**2
      c2_KPI = (lapp_KPI - lap_kpi**2)/2.d0/mpic**4

c /KPISC_EM/; parameters for Kpi scalar FF from 
c  http://arxiv.org/pdf/1103.4855.pdf
      lnC = 0.20193d0
      lambda0 = 0.013139d0

c   /SCAL_3PI/; parameters of sigma meson for 3 pion mode
      a00_3piscal = 0.220
      b00_3piscal = 0.268/mmpi_av**2
      c00_3piscal = -0.0139/mmpi_av**4
      d00_3piscal = -0.00139/mmpi_av**6
      x00_3piscal = 36.77*mmpi_av**2
      a02_3piscal = -0.0444
      b02_3piscal = -0.0857/mmpi_av**2
      c02_3piscal = -0.00221/mmpi_av**4
      d02_3piscal = -0.000129/mmpi_av**6
      x02_3piscal = -21.62*mmpi_av**2
      MMF0 = 0.441

c /SCAL_3PI/; parameters for the scalar part 3 pion modes
c   Pablo private
      ALPHA0_3PI = 1.
      ALPHA1_3PI = 1.
      GAMMA0_3PI = 1.
      GAMMA1_3PI = 1.


C FFVEC: dipswitch for Final State interaction in two scalar modes
C     with FSI (default FFVEC =1) and
C     without FSI (FFVEC =0)

      FFVEC = 1

C  FFKPIVEC : parameter to choose the parametrization for 
C             vector Kpi form factor with FSI effects
C    FFKPIVEC = 0 parametrization Eqs.(17),(18) of [4]
C    FFKPIVEC = 1 parametrization Eq.(5) of [5]
C    FFKPIVEC = 2 parmetrization [4], total result
      FFKPIVEC = 2
C  FFKPISCAL : parameter to choose the parametrization for 
C             scalar Kpi form factor with FSI effects
C  FFKPISCAL = 0 no scalar contribution
C  FFKPISCAL = 1 parametrization of Mathias Jamin,adopted his private code
C  FFKPISCAL = 2 parametrization of Emilie Passerman,
C               adopted her private code []
      FFKPISCAL = 1 

C FFKKVEC: dipswitch for K0K- mode
C     with rho' and rho'' (FFKKVEC =1) and
C     without rho' and rho''  (default FFKKVEC =0)
      FFKKVEC = 0

C FF3PISCAL: dipswitch for the scalar contribution for 3 pion modes
C     with the scalar contribution ( default FF3PISCAL = 2) 
c     FF3PISCAL = 2 BW parametrization for sigma meson
c     FF3PISCAL = 1 simplified RCHT results
C     FF3PISCAL =0  no sigma contribution
      FF3PISCAL = 2

C  Implemetation of another parametrization rho1, not checked yet by tests 
C FF3PIRHOPR: dipswitch for the parametrization for rho' contribution 
C     For 3 pion modes
C     general parametrization  ( default FF3RHOPR =1) and
C     simplified  (FF3PIRHOPR =0)
      FF3PIRHOPR = 0


C FF2PIRHO: dipswitch for the two pion form factor (default FF2PIRHO = 1) 
C FF2PIRHO =1   RChL parametrization 
C FF2PIRHO = 2  Belle parametrization, 
C               all parameters par (1...11) of fit are free 
C FF2PIRHO = 3  Belle parametrization, 
C               parameters of fit are free 
C               except for fixed par(1)=F_pi(0)=1
      FF2PIRHO =2


C FCOUL: dipswitch for the Coulomb interaction 
C FCOUL = 1  with 
C FCOUL = 0  without
      FCOUL = 0

      call rchl_REparam(1,IWARM,KAK)
      return
      end

      subroutine rchl_REparam(IMODE,IWARM,KAK)
      include '../parameter.inc'
      common / PARAMS / P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12,P13,P14,P15,P16,IUSE
      INTEGER                                                                  IUSE
      DOUBLE PRECISION  P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12,P13,P14,P15,P16
      DATA IUSE /0/

      IF(IUSE.EQ.0) RETURN
      
      IF (IMODE.EQ.-1) THEN
        IWARM=IWARM
      ELSE

C FF3PISCAL: dipswitch for the scalar contribution for 3 pion modes
C     with the scalar contribution ( default FF3PISCAL = 2) 
c     FF3PISCAL = 2 BW parametrization for sigma meson
c     FF3PISCAL = 1 simplified RCHT results
C     FF3PISCAL =0  no sigma contribution
c      FF3PISCAL = 2

C     CANDIDATES FOR PARAMETERS TO FIT with default values
C* Parameteres for the sigma contribution, using BW for sigma
      alpsig  = P1
      betasig = P2
      gamsig  = P3
      delsig  = P4
      rsigma  = P5


      MRO   = P6     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MRHO1 = P7     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      GRHO1 = P8       !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
C  /MASS_RES/ 
      GRO   = P9     !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MMA1  = P10    !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
      MSIG  = P11
      GSIG  = P12

C  /RCHT_3PI/; model parameters; their value are from fit,
c             reference [1], set 1
C              CHANGE OF THEIR VALUES REQUIRES 
C              RERUN /tabler/a1/da1wid_tot_rho1_gauss.f
      FPI_RPT  = P13
      FV_RPT   = P14
      FA_RPT   = P15
      BETA_RHO = P16
      FK_RPT = FPI_RPT*1.198d0  
      GV_RPT   = FPI_RPT*FPI_RPT/FV_RPT

      ENDIF
      return
      end
Commit	Line	Data
0ca57c2f	1	SUBROUTINE CH3PISET(JJ)
	2	C information on 3 pion sub-channel under construction obtained
	3	C J3PI=1 means 3 prong
	4	C J3PI=2 means 1 prong
	5	C to be initialized in routine DPHSAA of tauola.f
	6	COMMON /CHANOPT/ J3PI
	7	INTEGER J3PI
	8	J3PI=JJ
	9	end
	10
	11	SUBROUTINE CH3PIGET(JJ)
	12	C information on 3 pion sub-channel under construction obtained
	13	C J3PI=1 means 3 prong
	14	C J3PI=2 means 1 prong
	15	C to be initialized in routine DPHSAA of tauola.f
	16	COMMON /CHANOPT/ J3PI
	17	INTEGER J3PI
	18	IF (J3PI.EQ.1.OR.J3PI.EQ.2) THEN
	19	JJ=J3PI
	20	ELSE
	21	write(,) 'FROM value_parameter.f CH3PIGET, wrong J3PI=',J3PI
	22	stop
	23	ENDIF
	24	end
	25
	26	SUBROUTINE GETFF2PIRHO(JJ)
	27	IMPLICIT NONE
	28	include '../parameter.inc'
	29	INTEGER JJ
	30	JJ = FF2PIRHO
	31	END
	32
	33	SUBROUTINE OLACHNL(SIGN)
	34	C provides sign of tau, to be used in CP dependent parts of current.
	35	COMMON / JAKI / JAK1,JAK2,JAKP,JAKM,KTOM
	36	INTEGER JAK1,JAK2,JAKP,JAKM,KTOM
	37	COMMON / IDFC / IDFF
	38	INTEGER KTO
	39	REAL SIGN
	40	IF (KTOM.EQ.1.OR.KTOM.EQ.-1) THEN
	41	SIGN= IDFF/ABS(IDFF)
	42	ELSEIF (KTOM.EQ.2) THEN
	43	SIGN=-IDFF/ABS(IDFF)
	44	ELSE
	45	PRINT *, 'STOP IN OLACHNL: KTOM=',KTOM
	46	STOP
	47	ENDIF
	48	END
	49
	50	FUNCTION COEFrr(I,J)
	51	C clebsh gordan (or so ...) coefs for 3 scalar final states
	52	implicit none
	53	C TAUOLA RChL COEF(I,J) = COEFr(I,J)
	54	REAL COEFr(1:5,0:7)
	55	REAL COEFrr
	56	DATA PI /3.141592653589793238462643/
	57	REAL PI
	58	DATA ICONT /0/
	59	INTEGER ICONT
	60	INTEGER I,J
	61	REAL FPIr
	62
	63	C initialization of FPI matrix defined in ...
	64	C FPIc is to be used with cleo initialization
65	C FPIr is to be used with RChL initialization
66	C actual choice is made in ???
67
68
69	DATA FPIr /92.4E-3/
70
71
72	C initialization of COEF matrix defined in ...
73	C COEFc is to be used with cleo initialization
74	C COEFr is to be used with RChL initialization
75	IF (ICONT.EQ.0) THEN
76	ICONT=1
77	C
78	C******* COEFr(I,J) *****
79
80	COEFr(1,0)= 1.
81	COEFr(2,0)= -1.
82	COEFr(3,0)= 0.
83	COEFr(4,0)= 1.
84	COEFr(5,0)= 0.
85
86	COEFr(1,1)= 1.
87	COEFr(2,1)= -1.
88	COEFr(3,1)= 0.
89	COEFr(4,1)= 1.
90	COEFr(5,1)= 1.
91	C
92	COEFr(1,2)=1.
93	COEFr(2,2)= -1.
94	COEFr(3,2)= 0.0
95	COEFr(4,2)= 1.
96	COEFr(5,2)=1.
97	C
98	COEFr(1,3)= 0.
99	COEFr(2,3)= 1.
100	COEFr(3,3)= -1.
101	COEFr(4,3)= 1.
102	COEFr(5,3)= - 1.
103	C
104	COEFr(1,4)= 1.0/SQRT(2.)/3.0
105	COEFr(2,4)=-1.0/SQRT(2.)/3.0
106	COEFr(3,4)= 0.0
107	COEFr(4,4)= 0.0
108	COEFr(5,4)= 0.0
109	C
110	COEFr(1,5)=-SQRT(2.)/3.0
111	COEFr(2,5)= SQRT(2.)/3.0
112	COEFr(3,5)= 0.0
113	COEFr(4,5)= 0.0
114	COEFr(5,5)=-SQRT(2.)
115	C
116	COEFr(1,6)= 1./3.
117	COEFr(2,6)=-2./3.
118	COEFr(3,6)= 2./3.
119	COEFr(4,6)= 0.0
120	COEFr(5,6)=-2.0
121	C
122	COEFr(1,7)= 0.0
123	COEFr(2,7)= 0.0
124	COEFr(3,7)= 0.0
125	COEFr(4,7)= 0.0
126	COEFr(5,7)=-SQRT(2.0/3.0)
127	ENDIF
128
129	COEFrr=COEFr(I,J)
130	END
131
132	subroutine rchl_parameters(KAK)
133	implicit none
134	C==============================================================================
135	C Initialization, of '../parameter.inc' common block group
136	C
137	C KAK may be equal to JAK of TAUOLA namespace, but it is not always the case
138	C Hard-coded fit parameters:
139	C rho, rhoprime, f2(1275), f0(1186), sigma(made up!)
140	C The value of both the mass and width of resonances are taken
141	C from fit to ALEPH data (ref [1], Set 1)
142	C References: [1] arXiv: 0911.4436 [hep-ph] D. Gomez Dumm et al
143	C [2] arXiv: 0911.2640 [hep-ph] D. Gomez Dumm et al.
144	C [3] P Roig, talk PhiPsi2011, Novosibirsk
145	C [4] arXiv:0807.4883 [hep-ph] Diogo R. Boito et al.
146	C [5] arXiv:0803.1786 [hep-ph] M. Jamin et al.
147	C WARNING: some of parameters require RERUN of da1wid_tot_rho1_gauss.f
148	C pretabulating Q dependent a1 width,
149	C directory RChL-currents/tabler/a1
150	C==============================================================================
151	include '../parameter.inc'
152	INTEGER KAK
153	DATA IWARM/0/
154	INTEGER IWARM
155	INTEGER J3PI
156	COMMON /CHANOPT/ J3PI
157
158	IF(KAK.EQ.4) THEN
159	C /MASS_RES/; resonances parameters initialization:
160	C ! at present only for two pion mode non-default
161	c ! values are used:
162	mro = 0.77554d0
163	mrho1 = 1.453d0
164	grho1 = 0.50155D0
165	c /PAR_RHOPRIME/; parameters of rho' and rho''
166	C used for 2 pion form factor, reference [3]
167	COEF_GA = 0.14199D0
168	COEF_DE = -0.12623D0
169	phi_1 = -0.17377D0
170	phi_2 = 0.27632D0
171	grho2 = 0.41786D0
172	mrho2 = 1.8105d0
173	ELSE IF(KAK.EQ.5) THEN
174	MRO = 0.771849d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
175	MRHO1 = 1.35d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
176	GRHO1 = 0.448379d0 !0.473287d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
177
178	ELSE
179	MRO = 0.775 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
180	MRHO1 = 1.465 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
181	GRHO1 = 0.4 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
182
183	c /PAR_RHOPRIME/; parameters of rho' and rho''
184	C used for 2 kaon form factor, reference [3]
185	c FOR THE MOMENT THEIR NUMERICAL VALUES COINCIDE WITH
186	c ONES FOR THE TWO PION MODE !!!!
187	COEF_GA = 0.14199D0
188	COEF_DE = -0.12623D0
189	phi_1 = -0.17377D0
190	phi_2 = 0.27632D0
191	grho2 = 0.41786D0
192	mrho2 = 1.8105d0
193	ENDIF
194
195
196	IF(KAK.EQ.70) THEN ! non default values to be used
197	! for KPI MODE NO FSR INTERACTION
198	c /PAR_KPI/; parameters for Kpi mode, reference [4], table 4, row2
199	MKST = 0.943d0
200	MKSTPR = 1.374D0
201	GAMMA_KST = 0.06672d0
202	GAMMA_KSTPR = 0.240d0
203	GAMMA_RCHT =-0.039d0
204	ELSE IF(KAK.EQ.71) THEN ! non default values to be used
205	! for KPI MODE WITH FSR INTERACTION
206	c parameters for Kpi mode, reference [5]
207	MKST = 0.8953d0
208	GAMMA_KST = 0.0475d0
209	MKSTPR = 1.307d0
210	GAMMA_KSTPR = 0.206d0
211	GAMMA_RCHT = -0.043d0
212	ELSE
213	C /MASS_SCAL/; stable particles - final scalars
214	Mksp = 0.89166d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
215	Mks0 = 0.89610d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
216	MKST = (Mksp +Mks0)/2.
217	MKSTPR = 1.374d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
218	GAMMA_KST = 0.06672
219	GAMMA_KSTPR = 0.240
220	GAMMA_RCHT = -0.043 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
221	ENDIF
222
223	C /RCHT_3PI/; model parameters; their value are from fit,
224	c reference [1], set 1
225	C CHANGE OF THEIR VALUES REQUIRES
226	C RERUN /tabler/a1/da1wid_tot_rho1_gauss.f
227
228	IF(KAK.EQ.5) THEN
229	FPI_RPT = 0.091337d0
230	FV_RPT = 0.168652d0
231	FA_RPT = 0.131425d0
232	BETA_RHO = -0.318551d0
233	ELSE
234	FPI_RPT = 0.0924
235	FV_RPT = 0.18
236	FA_RPT = 0.149
237	BETA_RHO = -0.25
238	ENDIF
239
240	FK_RPT = FPI_RPT*1.198d0
241	GV_RPT = FPI_RPT*FPI_RPT/FV_RPT
242
243	c$$$c It has to be used for a new parametrization of rho1 for 3pions,
244	C$$$c that is not checked yet
245	c$$$c IF(KAK.EQ.5) THEN ! high energy behaviour imposes these relations
246	c$$$c GV_RPT = 0.066
247	c$$$c FV1_RPT = 0.18D0
248	c$$$c GV1_RPT = (FPI_RPTFPI_RPT- FV_RPTGV_RPT)/FV1_RPT
249	c$$$c ELSE
250	c$$$c GV_RPT = FPI_RPT*FPI_RPT/FV_RPT
251	c$$$c ENDIF
252
253	c /SCAL_3PI/; parameters of sigma meson for 3 pion mode
254	C* Parameteres for the sigma contribution, using BW for sigma
255	IF(KAK.EQ.5) THEN
256	IF (J3PI.EQ.1) THEN
257	alpsig = -8.795938d0
258	betasig = 9.763701d0
259	gamsig = 1.264263d0
260	delsig = 0.656762d0
261	rsigma = 1.866913d0
262	ELSE IF (J3PI.EQ.2) THEN
263	alpsig = 1.139486d0*0.63d0
264	betasig = 1.139486d0*0.63d0
265	gamsig = 0.889769d0*0.63d0
266	delsig = 0.889769d0*0.63d0
267	rsigma = 0.000013d0
268	ENDIF
269	ENDIF
270
271	C /MASS_RES/
272	IF(KAK.EQ.5) THEN
273	MMA1 = 1.091865d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
274	IF (J3PI.EQ.1) THEN
275	MSIG = 0.487512d0
276	GSIG = 0.70d0
277	ELSE IF(J3PI.EQ.2) THEN
278	MSIG = 0.55d0
279	GSIG = 0.7d0
280	ENDIF
281	ELSE
282	MMA1 = 1.12
283	MSIG = 0.475
284	GSIG = 0.550
285	ENDIF
286	call rchl_REparam(0,IWARM,KAK)
287	IF (IWARM.EQ.1) RETURN ! parameters below do not need
288	IWARM=1 ! re-initialization
289
290	C /MASS_RES/
291	GRO = 0.149d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
292	MF2 = 1.275d0
293	GF2 = 0.185d0
294	MF0 = 1.186d0
295	GF0 = 0.350d0
296	MSG = 0.860d0
297	GSG = 0.880d0
298	MPHI = 1.019d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
299	GPHI = 0.0042d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
300	MOM = 0.781940d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
301	GOM = 0.00843d0 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
302
303	C /RES_MIXING_RCHT/; a parameter defines w-phi angle mixing
304	THETA = 35.*PI/180.
305
306	C /FF0SCKPI/ a parameter normalized FFSC_KPI
307	F00 = 0.972
308
309
310
311
312	C /MASS_SCAL/; stable particles - final scalars
313	C CHANGE OF THEIR VALUES (useful for some tests) REQUIRES,
314	C RERUN /tabler/a1/da1wid_tot_rho1_gauss.f
315
316	MPIZ = 0.1349766d0 !PKORB(1,7) ! NEUTRAL PION MASS
317	MPIC = 0.13957018d0 !PKORB(1,8) ! CHARGED PION MASS
318	MMPI_AV = (MPIZ+2.*MPIC)/3.d0
319	MKZ = 0.497648d0 !PKORB(1,12) ! NEUTRAL KAON MASS
320	MKC = 0.493677d0 !PKORB(1,11) ! CHARGED KAON MASS
321	MMK = (MKC+MKZ)/2.d0
322	MTAU = 1.777
323	MNUTA = 0.001
324	META = 0.547d0
325
326
327
328	c /PAR_KKPI/; parameters to describe KKpi modes, reference [2]
329	C CHANGE OF THEIR VALUES REQUIRES
330	C RERUN /tabler/a1/da1wid_tot_rho1_gauss.f
331
332	G2 = mro/(192.pipisqrt(2.)FV_RPT)*3.
333	G13 = -2.*g2
334	G4 = -0.72
335	G5 = -0.6-2.*g4
336	C125 = 0.
337	C1256 = -3/96./pi*2FV_RPTMRO/SQRT(2.)/FPI_RPT*2
338	C1235 = 0.
339	C4 = -0.07
340	D123 = 0.05
341	D3 = -MRO*2/(64.PIPIFPI_RPT**2)
342
343
344	c /PAR_KPI/; parameters to describe Kpi mode, reference [4]
345	Ht0 = -1.2400398216503017D-2
346	C Ht0 = !!!!! TO ADD A FORMULAE FOR Ht0 (Jamin's email) !!!!!
347	lap_KPI = 24.66e-3
348	lapp_KPI = 11.99e-4
349	c1_KPI = lap_KPI/mpic**2
350	c2_KPI = (lapp_KPI - lap_kpi2)/2.d0/mpic4
351
352	c /KPISC_EM/; parameters for Kpi scalar FF from
353	c http://arxiv.org/pdf/1103.4855.pdf
354	lnC = 0.20193d0
355	lambda0 = 0.013139d0
356
357	c /SCAL_3PI/; parameters of sigma meson for 3 pion mode
358	a00_3piscal = 0.220
359	b00_3piscal = 0.268/mmpi_av**2
360	c00_3piscal = -0.0139/mmpi_av**4
361	d00_3piscal = -0.00139/mmpi_av**6
362	x00_3piscal = 36.77mmpi_av*2
363	a02_3piscal = -0.0444
364	b02_3piscal = -0.0857/mmpi_av**2
365	c02_3piscal = -0.00221/mmpi_av**4
366	d02_3piscal = -0.000129/mmpi_av**6
367	x02_3piscal = -21.62mmpi_av*2
368	MMF0 = 0.441
369
370	c /SCAL_3PI/; parameters for the scalar part 3 pion modes
371	c Pablo private
372	ALPHA0_3PI = 1.
373	ALPHA1_3PI = 1.
374	GAMMA0_3PI = 1.
375	GAMMA1_3PI = 1.
376
377
378	C FFVEC: dipswitch for Final State interaction in two scalar modes
379	C with FSI (default FFVEC =1) and
380	C without FSI (FFVEC =0)
381
382	FFVEC = 1
383
384	C FFKPIVEC : parameter to choose the parametrization for
385	C vector Kpi form factor with FSI effects
386	C FFKPIVEC = 0 parametrization Eqs.(17),(18) of [4]
387	C FFKPIVEC = 1 parametrization Eq.(5) of [5]
388	C FFKPIVEC = 2 parmetrization [4], total result
389	FFKPIVEC = 2
390	C FFKPISCAL : parameter to choose the parametrization for
391	C scalar Kpi form factor with FSI effects
392	C FFKPISCAL = 0 no scalar contribution
393	C FFKPISCAL = 1 parametrization of Mathias Jamin,adopted his private code
394	C FFKPISCAL = 2 parametrization of Emilie Passerman,
395	C adopted her private code []
396	FFKPISCAL = 1
397
398	C FFKKVEC: dipswitch for K0K- mode
399	C with rho' and rho'' (FFKKVEC =1) and
400	C without rho' and rho'' (default FFKKVEC =0)
401	FFKKVEC = 0
402
403	C FF3PISCAL: dipswitch for the scalar contribution for 3 pion modes
404	C with the scalar contribution ( default FF3PISCAL = 2)
405	c FF3PISCAL = 2 BW parametrization for sigma meson
406	c FF3PISCAL = 1 simplified RCHT results
407	C FF3PISCAL =0 no sigma contribution
408	FF3PISCAL = 2
409
410	C Implemetation of another parametrization rho1, not checked yet by tests
411	C FF3PIRHOPR: dipswitch for the parametrization for rho' contribution
412	C For 3 pion modes
413	C general parametrization ( default FF3RHOPR =1) and
414	C simplified (FF3PIRHOPR =0)
415	FF3PIRHOPR = 0
416
417
418	C FF2PIRHO: dipswitch for the two pion form factor (default FF2PIRHO = 1)
419	C FF2PIRHO =1 RChL parametrization
420	C FF2PIRHO = 2 Belle parametrization,
421	C all parameters par (1...11) of fit are free
422	C FF2PIRHO = 3 Belle parametrization,
423	C parameters of fit are free
424	C except for fixed par(1)=F_pi(0)=1
425	FF2PIRHO =2
426
427
428	C FCOUL: dipswitch for the Coulomb interaction
429	C FCOUL = 1 with
430	C FCOUL = 0 without
431	FCOUL = 0
432
433	call rchl_REparam(1,IWARM,KAK)
434	return
435	end
436
437	subroutine rchl_REparam(IMODE,IWARM,KAK)
438	include '../parameter.inc'
439	common / PARAMS / P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12,P13,P14,P15,P16,IUSE
440	INTEGER IUSE
441	DOUBLE PRECISION P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12,P13,P14,P15,P16
442	DATA IUSE /0/
443
444	IF(IUSE.EQ.0) RETURN
445
446	IF (IMODE.EQ.-1) THEN
447	IWARM=IWARM
448	ELSE
449
450	C FF3PISCAL: dipswitch for the scalar contribution for 3 pion modes
451	C with the scalar contribution ( default FF3PISCAL = 2)
452	c FF3PISCAL = 2 BW parametrization for sigma meson
453	c FF3PISCAL = 1 simplified RCHT results
454	C FF3PISCAL =0 no sigma contribution
455	c FF3PISCAL = 2
456
457	C CANDIDATES FOR PARAMETERS TO FIT with default values
458	C* Parameteres for the sigma contribution, using BW for sigma
459	alpsig = P1
460	betasig = P2
461	gamsig = P3
462	delsig = P4
463	rsigma = P5
464
465
466	MRO = P6 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
467	MRHO1 = P7 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
468	GRHO1 = P8 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
469	C /MASS_RES/
470	GRO = P9 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
471	MMA1 = P10 !CHANGE REQUIRES RERUN da1wid_tot_rho1_gauss.f
472	MSIG = P11
473	GSIG = P12
474
475	C /RCHT_3PI/; model parameters; their value are from fit,
476	c reference [1], set 1
477	C CHANGE OF THEIR VALUES REQUIRES
478	C RERUN /tabler/a1/da1wid_tot_rho1_gauss.f
479	FPI_RPT = P13
480	FV_RPT = P14
481	FA_RPT = P15
482	BETA_RHO = P16
483	FK_RPT = FPI_RPT*1.198d0
484	GV_RPT = FPI_RPT*FPI_RPT/FV_RPT
485
486	ENDIF
487	return
488	end