n_particles:='3$ 
n_total_spin:='6$ 
http:='nil$ 
sub_model:='ok$ 
direct_model:='ok$ 
just_v_type:='nil$ 
nnvertex:='81$ 
pdl_add:='((pr1 0 p1m w1d 0 0 0 0 0 0 1 0 0 0) (pr2 2 0 0 0
0 0 0 0 0 1 0 0 0) (pr3 2 p3m w3d 0 0 0 0 0 0 1 0 0 0) (pr4
2 p4m w4d 0 0 0 0 0 0 1 0 0 0) (pr5 2 p5m w5d 0 0 0 0 0 0 1
0 0 0) (pr6 2 p6m w6d 0 0 0 0 0 0 1 0 0 0) (pr7 2 p7m w7d 0
0 0 0 0 0 1 0 0 0) (pr8 2 p8m w8d 0 0 0 0 0 0 1 0 0 0) (pr9
2 p9m w9d 0 0 0 0 0 0 1 0 0 0) (pr10 6 p10m w10d 0 0 0 0 0 0
1 0 0 0) (pr11 2 p11m w11d 0 0 0 0 0 0 1 0 0 0) (pr12 6 p12m
w12d 0 0 0 0 0 0 1 0 0 0) (pr13 0 p13m w13d 0 0 0 0 0 0 1 0
0 0))$ 
spflist_add:='(list)$ 
v_t_list:='nil$ 
chargelist:='((pr1 0) (pr2 0) (pr3 0) (pr4 0) (pr5 0) (pr6 0
) (pr7 0) (pr8 0) (pr9 0) (pr10 0) (pr11 0) (pr12 0) (pr13 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 p3m w3d p4m w4d p5m w5d p6m w6d p7m
w7d p8m w8d p9m w9d p10m w10d p11m w11d p12m w12d p13m w13d)
$ 
platexname_add:='( 
      (pr1  "\eta") 
      (pr2  "\gamma") 
      (pr3  "\omega[782]") 
      (pr4  "\phi[1020]") 
      (pr5  "h_1[1170]") 
      (pr6  "h_1[1380]") 
      (pr7  "\omega[1420]") 
      (pr8  "h_1[1595]") 
      (pr9  "\omega_[1650]") 
      (pr10  "\omega_3[1670]") 
      (pr11  "\phi[1680]") 
      (pr12  "\phi_3[1850]") 
      (pr13  "\eta_{c}[2980]") 
     )$ 
 
 
inputlist_add:='( 
      (p1m  0.5473) 
      (w1d  0.0001) 
      (p3m  0.78194) 
      (w3d  0.0084) 
      (p4m  1.01941) 
      (w4d  0.0043) 
      (p5m  1.17) 
      (w5d  0.36) 
      (p6m  1.386) 
      (w6d  0.17) 
      (p7m  1.419) 
      (w7d  0.174) 
      (p8m  1.594) 
      (w8d  0.384) 
      (p9m  1.652) 
      (w9d  0.05) 
      (p10m  1.667) 
      (w10d  0.168) 
      (p11m  1.68) 
      (w11d  0.15) 
      (p12m  1.854) 
      (w12d  0.087) 
      (p13m  2.9792) 
      (w13d  0.013) 
     )$ 
 
pa_decay_list:='((pr3 (67 (pr2 pr1) 2) (15 (ef ef) 2)) (pr4
(64 (pr2 pr1) 2) (16 (ef ef) 2)) (pr5 (60 (pr2 pr1) 2)) (pr6
(55 (pr2 pr1) 2) (53 (pr3 pr1) 2)) (pr7 (48 (pr2 pr1) 2) (46
(pr3 pr1) 2) (17 (ef ef) 2)) (pr8 (40 (pr2 pr1) 2) (38 (pr3
pr1) 2) (36 (pr4 pr1) 2)) (pr9 (33 (pr2 pr1) 2) (31 (pr3 pr1
) 2) (29 (pr4 pr1) 2) (18 (ef ef) 2)) (pr10 (81 (pr2 pr1) 2)
(79 (pr3 pr1) 2) (77 (pr4 pr1) 2)) (pr11 (26 (pr2 pr1) 2) (
24 (pr3 pr1) 2) (22 (pr4 pr1) 2) (19 (ef ef) 2)) (pr12 (74 (
pr2 pr1) 2) (72 (pr3 pr1) 2) (70 (pr4 pr1) 2) (68 (pr5 pr1)
2)) (pr13 (80 (pr10 pr2) 2) (78 (pr10 pr3) 2) (76 (pr10 pr4)
2) (75 (pr10 pr5) 2) (73 (pr12 pr2) 2) (71 (pr12 pr3) 2) (69
(pr12 pr4) 2) (66 (pr3 pr2) 2) (65 (pr3 pr3) 2) (63 (pr4 pr2
) 2) (62 (pr4 pr3) 2) (61 (pr4 pr4) 2) (59 (pr5 pr2) 2) (58
(pr5 pr3) 2) (57 (pr5 pr4) 2) (56 (pr5 pr5) 2) (54 (pr6 pr2)
2) (52 (pr6 pr3) 2) (51 (pr6 pr4) 2) (50 (pr6 pr5) 2) (49 (
pr6 pr6) 2) (47 (pr7 pr2) 2) (45 (pr7 pr3) 2) (44 (pr7 pr4)
2) (43 (pr7 pr5) 2) (42 (pr7 pr6) 2) (41 (pr7 pr7) 2) (39 (
pr8 pr2) 2) (37 (pr8 pr3) 2) (35 (pr8 pr4) 2) (34 (pr8 pr5)
2) (32 (pr9 pr2) 2) (30 (pr9 pr3) 2) (28 (pr9 pr4) 2) (27 (
pr9 pr5) 2) (25 (pr11 pr2) 2) (23 (pr11 pr3) 2) (21 (pr11 
pr4) 2) (20 (pr11 pr5) 2)))$ 
htionvvlist_add:='((81 ((pr10 0 r1 v va vb) (pr2 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f82 g))) (80 ((pr10 0 r1 v va vb) (pr2 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f81 g))) (79 ((pr10 0 r1 v va vb) (pr3 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f80 g))) (78 ((pr10 0 r1 v va vb) (pr3 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f79 g))) (77 ((pr10 0 r1 v va vb) (pr4 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f78 g))) (76 ((pr10 0 r1 v va vb) (pr4 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f77 g))) (75 ((pr10 0 r1 v va vb) (pr5 0 r2 vc) (
pr13 0 r3)) (plus (times (cons r1 vc) (cons r2 v) (cons r2 
va) (cons r2 vb) f76 g) (times (cons r2 v) (cons r2 va) (
cons vb vc) f75 g))) (74 ((pr12 0 r1 v va vb) (pr2 0 r2 vc)
(pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f74 g))) (73 ((pr12 0 r1 v va vb) (pr2 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f73 g))) (72 ((pr12 0 r1 v va vb) (pr3 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f72 g))) (71 ((pr12 0 r1 v va vb) (pr3 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f71 g))) (70 ((pr12 0 r1 v va vb) (pr4 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f70 g))) (69 ((pr12 0 r1 v va vb) (pr4 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f69 g))) (68 ((pr12 0 r1 v va vb) (pr5 0 r2 vc) (
pr1 0 r3)) (plus (times (cons r1 vc) (cons r2 v) (cons r2 va
) (cons r2 vb) f68 g) (times (cons r2 v) (cons r2 va) (cons
vb vc) f67 g))) (67 ((pr3 0 r1 v) (pr2 0 r2 va) (pr1 0 r3))
(minus (times (eps r1 r2 v va) f66 g))) (66 ((pr3 0 r1 v) (
pr2 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f65
g))) (65 ((pr3 0 r1 v) (pr3 0 r2 va) (pr13 0 r3)) (minus (
times (eps r1 r2 v va) f64 g))) (64 ((pr4 0 r1 v) (pr2 0 r2
va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f63 g))) (63
((pr4 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)) (minus (times (eps
r1 r2 v va) f62 g))) (62 ((pr4 0 r1 v) (pr3 0 r2 va) (pr13 0
r3)) (minus (times (eps r1 r2 v va) f61 g))) (61 ((pr4 0 r1
v) (pr4 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va)
f60 g))) (60 ((pr5 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)) (
quotient (plus (times (cons r1 r1) (cons v va) f59 g) (times
2 (cons r1 va) (cons r2 v) f59 g) (minus (times (cons r3 r3)
(cons v va) f59 g))) 2)) (59 ((pr5 0 r1 v) (pr2 0 r2 va) (
pr13 0 r3)) (quotient (plus (times (cons r1 r1) (cons v va)
f58 g) (times 2 (cons r1 va) (cons r2 v) f58 g) (minus (
times (cons r3 r3) (cons v va) f58 g))) 2)) (58 ((pr5 0 r1 v
) (pr3 0 r2 va) (pr13 0 r3)) (plus (times (cons r1 va) (cons
r2 v) f57 g) (times (cons v va) f56 g))) (57 ((pr5 0 r1 v) (
pr4 0 r2 va) (pr13 0 r3)) (plus (times (cons r1 va) (cons r2
v) f55 g) (times (cons v va) f54 g))) (56 ((pr5 0 r1 v) (pr5
0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f53 g))
) (55 ((pr6 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)) (quotient (
plus (times (cons r1 r1) (cons v va) f52 g) (times 2 (cons 
r1 va) (cons r2 v) f52 g) (minus (times (cons r3 r3) (cons v
va) f52 g))) 2)) (54 ((pr6 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)
) (quotient (plus (times (cons r1 r1) (cons v va) f51 g) (
times 2 (cons r1 va) (cons r2 v) f51 g) (minus (times (cons
r3 r3) (cons v va) f51 g))) 2)) (53 ((pr6 0 r1 v) (pr3 0 r2
va) (pr1 0 r3)) (plus (times (cons r1 va) (cons r2 v) f50 g)
(times (cons v va) f49 g))) (52 ((pr6 0 r1 v) (pr3 0 r2 va)
(pr13 0 r3)) (plus (times (cons r1 va) (cons r2 v) f48 g) (
times (cons v va) f47 g))) (51 ((pr6 0 r1 v) (pr4 0 r2 va) (
pr13 0 r3)) (plus (times (cons r1 va) (cons r2 v) f46 g) (
times (cons v va) f45 g))) (50 ((pr6 0 r1 v) (pr5 0 r2 va) (
pr13 0 r3)) (minus (times (eps r1 r2 v va) f44 g))) (49 ((
pr6 0 r1 v) (pr6 0 r2 va) (pr13 0 r3)) (minus (times (eps r1
r2 v va) f43 g))) (48 ((pr7 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)
) (minus (times (eps r1 r2 v va) f42 g))) (47 ((pr7 0 r1 v)
(pr2 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) 
f41 g))) (46 ((pr7 0 r1 v) (pr3 0 r2 va) (pr1 0 r3)) (minus
(times (eps r1 r2 v va) f40 g))) (45 ((pr7 0 r1 v) (pr3 0 r2
va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f39 g))) (44
((pr7 0 r1 v) (pr4 0 r2 va) (pr13 0 r3)) (minus (times (eps
r1 r2 v va) f38 g))) (43 ((pr7 0 r1 v) (pr5 0 r2 va) (pr13 0
r3)) (plus (times (cons r1 va) (cons r2 v) f37 g) (times (
cons v va) f36 g))) (42 ((pr7 0 r1 v) (pr6 0 r2 va) (pr13 0
r3)) (plus (times (cons r1 va) (cons r2 v) f35 g) (times (
cons v va) f34 g))) (41 ((pr7 0 r1 v) (pr7 0 r2 va) (pr13 0
r3)) (minus (times (eps r1 r2 v va) f33 g))) (40 ((pr8 0 r1
v) (pr2 0 r2 va) (pr1 0 r3)) (quotient (plus (times (cons r1
r1) (cons v va) f32 g) (times 2 (cons r1 va) (cons r2 v) f32
g) (minus (times (cons r3 r3) (cons v va) f32 g))) 2)) (39 (
(pr8 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)) (quotient (plus (
times (cons r1 r1) (cons v va) f31 g) (times 2 (cons r1 va)
(cons r2 v) f31 g) (minus (times (cons r3 r3) (cons v va) 
f31 g))) 2)) (38 ((pr8 0 r1 v) (pr3 0 r2 va) (pr1 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f30 g) (times (cons v 
va) f29 g))) (37 ((pr8 0 r1 v) (pr3 0 r2 va) (pr13 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f28 g) (times (cons v 
va) f27 g))) (36 ((pr8 0 r1 v) (pr4 0 r2 va) (pr1 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f26 g) (times (cons v 
va) f25 g))) (35 ((pr8 0 r1 v) (pr4 0 r2 va) (pr13 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f24 g) (times (cons v 
va) f23 g))) (34 ((pr8 0 r1 v) (pr5 0 r2 va) (pr13 0 r3)) (
minus (times (eps r1 r2 v va) f22 g))) (33 ((pr9 0 r1 v) (
pr2 0 r2 va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f21
g))) (32 ((pr9 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)) (minus (
times (eps r1 r2 v va) f20 g))) (31 ((pr9 0 r1 v) (pr3 0 r2
va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f19 g))) (30
((pr9 0 r1 v) (pr3 0 r2 va) (pr13 0 r3)) (minus (times (eps
r1 r2 v va) f18 g))) (29 ((pr9 0 r1 v) (pr4 0 r2 va) (pr1 0
r3)) (minus (times (eps r1 r2 v va) f17 g))) (28 ((pr9 0 r1
v) (pr4 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va)
f16 g))) (27 ((pr9 0 r1 v) (pr5 0 r2 va) (pr13 0 r3)) (plus
(times (cons r1 va) (cons r2 v) f15 g) (times (cons v va) 
f14 g))) (26 ((pr11 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)) (minus
(times (eps r1 r2 v va) f13 g))) (25 ((pr11 0 r1 v) (pr2 0 
r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f12 g)))
(24 ((pr11 0 r1 v) (pr3 0 r2 va) (pr1 0 r3)) (minus (times (
eps r1 r2 v va) f11 g))) (23 ((pr11 0 r1 v) (pr3 0 r2 va) (
pr13 0 r3)) (minus (times (eps r1 r2 v va) f10 g))) (22 ((
pr11 0 r1 v) (pr4 0 r2 va) (pr1 0 r3)) (minus (times (eps r1
r2 v va) f9 g))) (21 ((pr11 0 r1 v) (pr4 0 r2 va) (pr13 0 r3
)) (minus (times (eps r1 r2 v va) f8 g))) (20 ((pr11 0 r1 v)
(pr5 0 r2 va) (pr13 0 r3)) (plus (times (cons r1 va) (cons 
r2 v) f7 g) (times (cons v va) f6 g))) (19 ((pr11 0 r1 v) (
ef -1 r2) (ef 1 r3)) (times (g l v) f5 g)) (18 ((pr9 0 r1 v)
(ef -1 r2) (ef 1 r3)) (times (g l v) f4 g)) (17 ((pr7 0 r1 v
) (ef -1 r2) (ef 1 r3)) (times (g l v) f3 g)) (16 ((pr4 0 r1
v) (ef -1 r2) (ef 1 r3)) (times (g l v) f2 g)) (15 ((pr3 0 
r1 v) (ef -1 r2) (ef 1 r3)) (times (g l v) f1 g)))$ 
vertices_add_comment:='"
In this model, We have added some effective 
vertices, to do some calculation related to some resonance states
"$ 
end; 
