TEvtGen/Tauola/tauola-fortran/new-currents/RChL-currents/rcht_common/value_parameter.f

   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_RPT*FPI_RPT- FV_RPT*GV_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.*pi*pi*sqrt(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**2*FV_RPT*MRO/SQRT(2.)/FPI_RPT**2
 338       C1235 =  0.
 339       C4    = -0.07
 340       D123  =  0.05
 341       D3    = -MRO**2/(64.*PI*PI*FPI_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_kpi**2)/2.d0/mpic**4
 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.77*mmpi_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.62*mmpi_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