n_particles:='3$ 
n_total_spin:='6$ 
http:='nil$ 
sub_model:='ok$ 
direct_model:='ok$ 
just_v_type:='nil$ 
nnvertex:='57$ 
pdl_add:='((pr1 0 p1m w1d 3 1 0 0 0 0 1 -1 0 0) (pr2 0 p2m 
w2d 0 0 0 0 0 0 1 0 0 0) (pr3 4 p3m w3d 0 0 0 0 0 0 1 0 0 0)
(pr4 0 p4m w4d 0 0 0 0 0 0 1 0 0 0) (pr5 4 p5m w5d 0 0 0 0 0
0 1 0 0 0) (pr6 4 p6m w6d 0 0 0 0 0 0 1 0 0 0) (pr7 4 p7m 
w7d 0 0 0 0 0 0 1 0 0 0) (pr8 4 p8m w8d 0 0 0 0 0 0 1 0 0 0)
(pr9 0 p9m w9d 0 0 0 0 0 0 1 0 0 0) (pr10 4 p10m w10d 0 0 0
0 0 0 1 0 0 0) (pr11 4 p11m w11d 0 0 0 0 0 0 1 0 0 0) (pr12
4 p12m w12d 0 0 0 0 0 0 1 0 0 0) (pr13 4 p13m w13d 0 0 0 0 0
0 1 0 0 0) (pr14 0 p14m w14d 0 0 0 0 0 0 1 0 0 0) (pr15 4 
p15m w15d 3 1 0 0 0 0 1 -1 0 0) (pr16 0 p16m w16d 0 0 0 0 0
0 1 0 0 0))$ 
spflist_add:='(list)$ 
v_t_list:='nil$ 
chargelist:='((pr1 1) (pr2 0) (pr3 0) (pr4 0) (pr5 0) (pr6 0
) (pr7 0) (pr8 0) (pr9 0) (pr10 0) (pr11 0) (pr12 0) (pr13 0
) (pr14 0) (pr15 1) (pr16 0) (gsg 0) (gs 0) (p 0) (z 0) (w 1
) (w -1) (z 0) (p 0) (qd (quotient (minus 1) 3)) (qu (
quotient 2 3)) (qu (quotient 2 3)) (qd (quotient (minus 1) 3
)) (ef -1) (nue 0) (ef -1) (h0 0) (numu 0) (mu -1) (nut 0) (
tau -1) (mu -1) (tau -1) (qc (quotient 2 3)) (qs (quotient (
minus 1) 3)) (qt (quotient 2 3)) (qb (quotient (minus 1) 3))
(qc (quotient 2 3)) (qt (quotient 2 3)) (qs (quotient (minus
1) 3)) (qb (quotient (minus 1) 3)))$ 
phyinput_add:='(p1m w1d p2m w2d p3m w3d p4m w4d p5m w5d p6m
w6d p7m w7d p8m w8d p9m w9d p10m w10d p11m w11d p12m w12d 
p13m w13d p14m w14d p15m w15d p16m w16d)$ 
platexname_add:='( 
      (pr1  "K") 
      (pr2  "\eta^{'}[958]") 
      (pr3  "f_2[1270]") 
      (pr4  "f_0[1370]") 
      (pr5  "f_2[1430]") 
      (pr6  "f^{'}_2[1525]") 
      (pr7  "f_2[1565]") 
      (pr8  "f_2[1640]") 
      (pr9  "f_0[1710]") 
      (pr10  "f_2[1810]") 
      (pr11  "f_2[1910]") 
      (pr12  "f_2[1950]") 
      (pr13  "f_2[2010]") 
      (pr14  "f_0[2020]") 
      (pr15  "K^{*}_2[1980]") 
      (pr16  "\eta_{c}[2980]") 
     )$ 
 
 
inputlist_add:='( 
      (p1m  0.49368) 
      (w1d  0.0001) 
      (p2m  0.95778) 
      (w2d  0.0002) 
      (p3m  1.275) 
      (w3d  0.1855) 
      (p4m  1.37) 
      (w4d  0.35) 
      (p5m  1.43) 
      (w5d  0.081) 
      (p6m  1.525) 
      (w6d  0.076) 
      (p7m  1.542) 
      (w7d  1.131) 
      (p8m  1.638) 
      (w8d  0.099) 
      (p9m  1.712) 
      (w9d  0.133) 
      (p10m  1.815) 
      (w10d  0.197) 
      (p11m  1.915) 
      (w11d  0.163) 
      (p12m  1.934) 
      (w12d  0.453) 
      (p13m  2.011) 
      (w13d  0.202) 
      (p14m  2.02) 
      (w14d  0.41) 
      (p15m  1.973) 
      (w15d  0.373) 
      (p16m  2.9792) 
      (w16d  0.013) 
     )$ 
 
pa_decay_list:='((pr3 (45 (pr1 pr1) 2)) (pr4 (21 (pr1 pr1) 2
)) (pr5 (43 (pr1 pr1) 2)) (pr6 (41 (pr1 pr1) 2)) (pr7 (39 (
pr1 pr1) 2)) (pr8 (37 (pr1 pr1) 2)) (pr9 (20 (pr1 pr1) 2)) (
pr10 (35 (pr1 pr1) 2)) (pr11 (33 (pr1 pr1) 2)) (pr12 (31 (
pr1 pr1) 2) (30 (pr2 pr2) 2)) (pr15 (49 (pr3 pr1) 2) (48 (
pr3 pr1) 1) (47 (pr5 pr1) 2) (46 (pr5 pr1) 1) (25 (pr2 pr1)
2) (24 (pr2 pr1) 1)) (pr13 (28 (pr1 pr1) 2) (27 (pr2 pr2) 2)
) (pr14 (19 (pr1 pr1) 2) (18 (pr2 pr2) 2)) (pr16 (57 (pr3 
pr3) 2) (56 (pr5 pr3) 2) (55 (pr5 pr5) 2) (54 (pr6 pr3) 2) (
53 (pr6 pr5) 2) (52 (pr7 pr3) 2) (51 (pr7 pr5) 2) (50 (pr8 
pr3) 2) (44 (pr3 pr2) 2) (42 (pr5 pr2) 2) (40 (pr6 pr2) 2) (
38 (pr7 pr2) 2) (36 (pr8 pr2) 2) (34 (pr10 pr2) 2) (32 (pr11
pr2) 2) (29 (pr12 pr2) 2) (26 (pr13 pr2) 2) (23 (pr15 pr1) 2
) (22 (pr15 pr1) 2) (17 (pr4 pr2) 2) (16 (pr9 pr2) 2) (15 (
pr14 pr2) 2)))$ 
htionvvlist_add:='((57 ((pr3 0 r1 v va) (pr3 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f49 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f48 g)))) (56 ((pr5 0 r1 v va) (pr3 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f47 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f46 g)))) (55 ((pr5 0 r1 v va) (pr5 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f45 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f44 g)))) (54 ((pr6 0 r1 v va) (pr3 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f43 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f42 g)))) (53 ((pr6 0 r1 v va) (pr5 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f41 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f40 g)))) (52 ((pr7 0 r1 v va) (pr3 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f39 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f38 g)))) (51 ((pr7 0 r1 v va) (pr5 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f37 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f36 g)))) (50 ((pr8 0 r1 v va) (pr3 0 r2 vb vc) (
pr16 0 r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (
eps r1 r2 v vb) f35 g)) (minus (times (cons va vc) (eps r1 
r2 v vb) f34 g)))) (49 ((pr1 -1 r3) (pr3 0 r2 vb vc) (pr15 1
r1 v va)) (plus (minus (times (cons r1 vc) (cons r2 va) (eps
r1 r2 v vb) f33 g)) (minus (times (cons va vc) (eps r1 r2 v
vb) f32 g)))) (48 ((pr15 -1 r1 v va) (pr3 0 r2 vb vc) (pr1 1
r3)) (plus (minus (times (cons r1 vc) (cons r2 va) (eps r1 
r2 v vb) f33 g)) (minus (times (cons va vc) (eps r1 r2 v vb)
f32 g)))) (47 ((pr1 -1 r3) (pr5 0 r2 vb vc) (pr15 1 r1 v va)
) (plus (minus (times (cons r1 vc) (cons r2 va) (eps r1 r2 v
vb) f31 g)) (minus (times (cons va vc) (eps r1 r2 v vb) f30
g)))) (46 ((pr15 -1 r1 v va) (pr5 0 r2 vb vc) (pr1 1 r3)) (
plus (minus (times (cons r1 vc) (cons r2 va) (eps r1 r2 v vb
) f31 g)) (minus (times (cons va vc) (eps r1 r2 v vb) f30 g)
))) (45 ((pr3 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (times (
cons r2 v) (cons r2 va) f29 g)) (44 ((pr3 0 r1 v va) (pr16 0
r2) (pr2 0 r3)) (times (cons r2 v) (cons r2 va) f28 g)) (43
((pr5 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v)
(cons r2 va) f27 g)) (42 ((pr5 0 r1 v va) (pr16 0 r2) (pr2 0
r3)) (times (cons r2 v) (cons r2 va) f26 g)) (41 ((pr6 0 r1
v va) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va
) f25 g)) (40 ((pr6 0 r1 v va) (pr16 0 r2) (pr2 0 r3)) (
times (cons r2 v) (cons r2 va) f24 g)) (39 ((pr7 0 r1 v va)
(pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va) f23
g)) (38 ((pr7 0 r1 v va) (pr16 0 r2) (pr2 0 r3)) (times (
cons r2 v) (cons r2 va) f22 g)) (37 ((pr8 0 r1 v va) (pr1 -1
r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va) f21 g)) (36
((pr8 0 r1 v va) (pr16 0 r2) (pr2 0 r3)) (times (cons r2 v)
(cons r2 va) f20 g)) (35 ((pr10 0 r1 v va) (pr1 -1 r2) (pr1
1 r3)) (times (cons r2 v) (cons r2 va) f19 g)) (34 ((pr10 0
r1 v va) (pr16 0 r2) (pr2 0 r3)) (times (cons r2 v) (cons r2
va) f18 g)) (33 ((pr11 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (
times (cons r2 v) (cons r2 va) f17 g)) (32 ((pr11 0 r1 v va)
(pr16 0 r2) (pr2 0 r3)) (times (cons r2 v) (cons r2 va) f16
g)) (31 ((pr12 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (times (
cons r2 v) (cons r2 va) f15 g)) (30 ((pr12 0 r1 v va) (pr2 0
r2) (pr2 0 r3)) (times (cons r2 v) (cons r2 va) f14 g)) (29
((pr12 0 r1 v va) (pr16 0 r2) (pr2 0 r3)) (times (cons r2 v)
(cons r2 va) f13 g)) (28 ((pr13 0 r1 v va) (pr1 -1 r2) (pr1
1 r3)) (times (cons r2 v) (cons r2 va) f12 g)) (27 ((pr13 0
r1 v va) (pr2 0 r2) (pr2 0 r3)) (times (cons r2 v) (cons r2
va) f11 g)) (26 ((pr13 0 r1 v va) (pr16 0 r2) (pr2 0 r3)) (
times (cons r2 v) (cons r2 va) f10 g)) (25 ((pr1 -1 r3) (pr2
0 r2) (pr15 1 r1 v va)) (times (cons r2 v) (cons r2 va) f9 g
)) (24 ((pr15 -1 r1 v va) (pr2 0 r2) (pr1 1 r3)) (times (
cons r2 v) (cons r2 va) f9 g)) (23 ((pr1 -1 r3) (pr16 0 r2)
(pr15 1 r1 v va)) (times (cons r2 v) (cons r2 va) f8 g)) (22
((pr15 -1 r1 v va) (pr16 0 r2) (pr1 1 r3)) (times (cons r2 v
) (cons r2 va) f8 g)) (21 ((pr4 0 r1) (pr1 -1 r2) (pr1 1 r3)
) (times f7 g)) (20 ((pr9 0 r1) (pr1 -1 r2) (pr1 1 r3)) (
times f6 g)) (19 ((pr14 0 r1) (pr1 -1 r2) (pr1 1 r3)) (times
f5 g)) (18 ((pr14 0 r1) (pr2 0 r2) (pr2 0 r3)) (times f4 g))
(17 ((pr16 0 r1) (pr4 0 r2) (pr2 0 r3)) (times f3 g)) (16 ((
pr16 0 r1) (pr9 0 r2) (pr2 0 r3)) (times f2 g)) (15 ((pr16 0
r1) (pr14 0 r2) (pr2 0 r3)) (times f1 g)))$ 
vertices_add_comment:='"
In this model, We have added some effective 
vertices, to do some calculation related to some resonance states
"$ 
end; 
