C++ Interface to Tauola
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
Generated by doxygen 1.9.1