n_particles:='3$ 
n_total_spin:='6$ 
http:='nil$ 
sub_model:='ok$ 
direct_model:='ok$ 
just_v_type:='nil$ 
nnvertex:='76$ 
pdl_add:='((pr1 0 p1m w1d 0 1 0 0 0 0 1 -1 0 0) (pr2 1 p2m 
w2d 3 1 0 0 0 0 1 0 1 0) (pr3 2 p3m w3d 0 0 0 0 0 0 1 0 0 0)
(pr4 4 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 0 p7m 
w7d 0 0 0 0 0 0 1 0 0 0) (pr8 8 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 7 p10m w10d 3 1 0
0 0 0 1 0 1 0) (pr11 5 p11m w11d 3 1 0 0 0 0 1 0 1 0) (pr12
3 p12m w12d 3 1 0 0 0 0 1 0 1 0) (pr13 1 p13m w13d 3 1 0 0 0
0 1 0 1 0) (pr14 1 p14m w14d 3 1 0 0 0 0 1 0 1 0) (pr15 5 
p15m w15d 3 1 0 0 0 0 1 0 1 0) (pr16 1 p16m w16d 3 1 0 0 0 0
1 0 1 0) (pr17 3 p17m w17d 3 1 0 0 0 0 1 -1 1 0) (pr18 1 
p18m w18d 3 1 0 0 0 0 1 -1 1 0))$ 
spflist_add:='(list)$ 
v_t_list:='nil$ 
chargelist:='((pr1 0) (pr2 1) (pr3 0) (pr4 0) (pr5 0) (pr6 0
) (pr7 0) (pr8 0) (pr9 0) (pr10 1) (pr11 1) (pr12 1) (pr13 1
) (pr14 1) (pr15 1) (pr16 1) (pr17 1) (pr18 1) (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 p17m w17d p18m w18d)
$ 
platexname_add:='( 
      (pr1  "K^0") 
      (pr2  "P") 
      (pr3  "\phi[1020]") 
      (pr4  "f_2[1910]") 
      (pr5  "f_2[1950]") 
      (pr6  "f_2[2010]") 
      (pr7  "f_0[2020]") 
      (pr8  "f_4[2050]") 
      (pr9  "J/\psi") 
      (pr10  "N[1990]") 
      (pr11  "N[2000]") 
      (pr12  "N[2080]") 
      (pr13  "N[2090]") 
      (pr14  "N[2100]") 
      (pr15  "\Delta[2000][+1]") 
      (pr16  "\Delta[2150][+1]") 
      (pr17  "\Sigma[1580][+1]") 
      (pr18  "\Sigma[1620][+1]") 
     )$ 
 
 
inputlist_add:='( 
      (p1m  0.49767) 
      (w1d  0.0001) 
      (p2m  0.93827) 
      (w2d  0.0001) 
      (p3m  1.01941) 
      (w3d  0.0043) 
      (p4m  1.915) 
      (w4d  0.163) 
      (p5m  1.934) 
      (w5d  0.453) 
      (p6m  2.011) 
      (w6d  0.202) 
      (p7m  2.02) 
      (w7d  0.41) 
      (p8m  2.025) 
      (w8d  0.194) 
      (p9m  3.09688) 
      (w9d  0.0001) 
      (p10m  1.99) 
      (w10d  0.38) 
      (p11m  2) 
      (w11d  0.3) 
      (p12m  2.08) 
      (w12d  0.27) 
      (p13m  2.09) 
      (w13d  0.3) 
      (p14m  2.1) 
      (w14d  0.26) 
      (p15m  2) 
      (w15d  0.325) 
      (p16m  2.15) 
      (w16d  0.16) 
      (p17m  1.58) 
      (w17d  0.015) 
      (p18m  1.62) 
      (w18d  0.04) 
     )$ 
 
pa_decay_list:='((pr3 (19 (pr1 pr1) 2) (15 (ef ef) 2)) (pr17
(41 (pr2 pr1) 2) (40 (pr2 pr1) 1)) (pr18 (32 (pr2 pr1) 2) (
31 (pr2 pr1) 1)) (pr4 (63 (pr2 pr2) 2) (23 (pr1 pr1) 2)) (
pr5 (62 (pr2 pr2) 2) (22 (pr1 pr1) 2)) (pr10 (75 (pr2 pr3) 2
) (74 (pr2 pr3) 1)) (pr11 (71 (pr2 pr3) 2) (70 (pr2 pr3) 1))
(pr15 (67 (pr2 pr3) 2) (66 (pr2 pr3) 1)) (pr6 (61 (pr2 pr2)
2) (21 (pr1 pr1) 2)) (pr7 (33 (pr2 pr2) 2) (17 (pr1 pr1) 2))
(pr8 (76 (pr2 pr2) 2) (27 (pr1 pr1) 2)) (pr12 (60 (pr2 pr3)
2) (59 (pr2 pr3) 1) (56 (pr17 pr1) 1) (55 (pr17 pr1) 2)) (
pr13 (53 (pr2 pr3) 2) (52 (pr2 pr3) 1) (39 (pr17 pr1) 1) (38
(pr17 pr1) 2)) (pr14 (49 (pr2 pr3) 2) (48 (pr2 pr3) 1) (37 (
pr17 pr1) 1) (36 (pr17 pr1) 2)) (pr16 (45 (pr2 pr3) 2) (44 (
pr2 pr3) 1) (35 (pr17 pr1) 1) (34 (pr17 pr1) 2) (30 (pr18 
pr1) 1) (29 (pr18 pr1) 2)) (pr9 (73 (pr10 pr2) 2) (72 (pr10
pr2) 2) (69 (pr11 pr2) 2) (68 (pr11 pr2) 2) (65 (pr15 pr2) 2
) (64 (pr15 pr2) 2) (58 (pr12 pr2) 2) (57 (pr12 pr2) 2) (54
(pr2 pr2) 2) (51 (pr13 pr2) 2) (50 (pr13 pr2) 2) (47 (pr14 
pr2) 2) (46 (pr14 pr2) 2) (43 (pr16 pr2) 2) (42 (pr16 pr2) 2
) (28 (pr8 pr3) 2) (26 (pr4 pr3) 2) (25 (pr5 pr3) 2) (24 (
pr6 pr3) 2) (20 (pr3 pr7) 2) (18 (pr1 pr1) 2) (16 (ef ef) 2)
))$ 
htionvvlist_add:='((76 ((pr2 -1 r1) (pr2 1 r2) (pr8 0 r3 v 
va vb vc)) (plus (times (cons r2 v) (cons r2 va) (cons r2 vb
) (cons r2 vc) f85 g) (times (g l v) (cons r2 va) (cons r2 
vb) (cons r2 vc) f86 g))) (75 ((pr3 0 r3 vc) (pr2 -1 r2) (
pr10 1 r1 v va vb)) (plus (minus (times (g l a vc) (cons r2
v) (cons r2 va) (cons r2 vb) f84 g)) (minus (times (g l a) (
cons r2 v) (cons r2 va) (cons r2 vb) (cons r2 vc) f82 g)) (
minus (times (g l a) (cons r2 v) (cons r2 va) (cons vb vc) 
f83 g)))) (74 ((pr10 -1 r1 v va vb) (pr2 1 r2) (pr3 0 r3 vc)
) (plus (times (g l a vc) (cons r2 v) (cons r2 va) (cons r2
vb) f84 g) (times (g l a) (cons r2 v) (cons r2 va) (cons r2
vb) (cons r2 vc) f82 g) (times (g l a) (cons r2 v) (cons r2
va) (cons vb vc) f83 g))) (73 ((pr9 0 r3 vc) (pr2 -1 r2) (
pr10 1 r1 v va vb)) (plus (minus (times (g l a vc) (cons r2
v) (cons r2 va) (cons r2 vb) f81 g)) (minus (times (g l a) (
cons r2 v) (cons r2 va) (cons r2 vb) (cons r2 vc) f79 g)) (
minus (times (g l a) (cons r2 v) (cons r2 va) (cons vb vc) 
f80 g)))) (72 ((pr10 -1 r1 v va vb) (pr2 1 r2) (pr9 0 r3 vc)
) (plus (times (g l a vc) (cons r2 v) (cons r2 va) (cons r2
vb) f81 g) (times (g l a) (cons r2 v) (cons r2 va) (cons r2
vb) (cons r2 vc) f79 g) (times (g l a) (cons r2 v) (cons r2
va) (cons vb vc) f80 g))) (71 ((pr3 0 r3 vb) (pr2 -1 r2) (
pr11 1 r1 v va)) (plus (minus (times (cons r2 v) (cons r2 va
) (cons r2 vb) f76 g)) (minus (times (cons r2 v) (cons va vb
) f77 g)) (times (g l vb) (cons r2 v) (cons r2 va) f78 g)))
(70 ((pr11 -1 r1 v va) (pr2 1 r2) (pr3 0 r3 vb)) (plus (
times (cons r2 v) (cons r2 va) (cons r2 vb) f76 g) (times (
cons r2 v) (cons va vb) f77 g) (times (g l vb) (cons r2 v) (
cons r2 va) f78 g))) (69 ((pr9 0 r3 vb) (pr2 -1 r2) (pr11 1
r1 v va)) (plus (minus (times (cons r2 v) (cons r2 va) (cons
r2 vb) f73 g)) (minus (times (cons r2 v) (cons va vb) f74 g)
) (times (g l vb) (cons r2 v) (cons r2 va) f75 g))) (68 ((
pr11 -1 r1 v va) (pr2 1 r2) (pr9 0 r3 vb)) (plus (times (
cons r2 v) (cons r2 va) (cons r2 vb) f73 g) (times (cons r2
v) (cons va vb) f74 g) (times (g l vb) (cons r2 v) (cons r2
va) f75 g))) (67 ((pr3 0 r3 vb) (pr2 -1 r2) (pr15 1 r1 v va)
) (plus (minus (times (cons r2 v) (cons r2 va) (cons r2 vb)
f70 g)) (minus (times (cons r2 v) (cons va vb) f71 g)) (
times (g l vb) (cons r2 v) (cons r2 va) f72 g))) (66 ((pr15
-1 r1 v va) (pr2 1 r2) (pr3 0 r3 vb)) (plus (times (cons r2
v) (cons r2 va) (cons r2 vb) f70 g) (times (cons r2 v) (cons
va vb) f71 g) (times (g l vb) (cons r2 v) (cons r2 va) f72 g
))) (65 ((pr9 0 r3 vb) (pr2 -1 r2) (pr15 1 r1 v va)) (plus (
minus (times (cons r2 v) (cons r2 va) (cons r2 vb) f67 g)) (
minus (times (cons r2 v) (cons va vb) f68 g)) (times (g l vb
) (cons r2 v) (cons r2 va) f69 g))) (64 ((pr15 -1 r1 v va) (
pr2 1 r2) (pr9 0 r3 vb)) (plus (times (cons r2 v) (cons r2 
va) (cons r2 vb) f67 g) (times (cons r2 v) (cons va vb) f68
g) (times (g l vb) (cons r2 v) (cons r2 va) f69 g))) (63 ((
pr2 -1 r1) (pr2 1 r2) (pr4 0 r3 v va)) (plus (times (cons r2
v) (cons r2 va) f65 g) (times (g l v) (cons r2 va) f66 g)))
(62 ((pr2 -1 r1) (pr2 1 r2) (pr5 0 r3 v va)) (plus (times (
cons r2 v) (cons r2 va) f63 g) (times (g l v) (cons r2 va) 
f64 g))) (61 ((pr2 -1 r1) (pr2 1 r2) (pr6 0 r3 v va)) (plus
(times (cons r2 v) (cons r2 va) f61 g) (times (g l v) (cons
r2 va) f62 g))) (60 ((pr3 0 r3 va) (pr2 -1 r2) (pr12 1 r1 v)
) (plus (minus (times (cons r2 v) (cons r2 va) f58 g i)) (
minus (times (cons v va) f59 g i)) (times (g l va) (cons r2
v) f60 g i))) (59 ((pr12 -1 r1 v) (pr2 1 r2) (pr3 0 r3 va))
(plus (times (cons r2 v) (cons r2 va) f58 g i) (times (cons
v va) f59 g i) (times (g l va) (cons r2 v) f60 g i))) (58 ((
pr9 0 r3 va) (pr2 -1 r2) (pr12 1 r1 v)) (plus (minus (times
(cons r2 v) (cons r2 va) f55 g i)) (minus (times (cons v va)
f56 g i)) (times (g l va) (cons r2 v) f57 g i))) (57 ((pr12
-1 r1 v) (pr2 1 r2) (pr9 0 r3 va)) (plus (times (cons r2 v)
(cons r2 va) f55 g i) (times (cons v va) f56 g i) (times (g
l va) (cons r2 v) f57 g i))) (56 ((pr1 -1 r3) (pr12 -1 r2 va
) (pr17 1 r1 v)) (plus (times (g l a) (cons r2 v) (cons r3 
va) f53 g i) (times (g l a) (cons v va) f54 g i))) (55 ((
pr17 -1 r1 v) (pr12 1 r2 va) (pr1 1 r3)) (plus (times (g l a
) (cons r2 v) (cons r3 va) f53 g i) (times (g l a) (cons v 
va) f54 g i))) (54 ((pr2 -1 r1) (pr2 1 r2) (pr9 0 r3 v)) (
plus (minus (times (g l r3 v) f52 g)) (times (g l v r3) f52
g) (times (g l v) f51 g))) (53 ((pr3 0 r3 v) (pr2 -1 r2) (
pr13 1 r1)) (plus (times (g l a r3 v) f50 g i) (minus (times
(g l a v r3) f50 g i)) (minus (times (g l a v) f49 g i)))) (
52 ((pr13 -1 r1) (pr2 1 r2) (pr3 0 r3 v)) (plus (times (g l
a r3 v) f50 g i) (minus (times (g l a v r3) f50 g i)) (times
(g l a v) f49 g i))) (51 ((pr9 0 r3 v) (pr2 -1 r2) (pr13 1 
r1)) (plus (times (g l a r3 v) f48 g i) (minus (times (g l a
v r3) f48 g i)) (minus (times (g l a v) f47 g i)))) (50 ((
pr13 -1 r1) (pr2 1 r2) (pr9 0 r3 v)) (plus (times (g l a r3
v) f48 g i) (minus (times (g l a v r3) f48 g i)) (times (g l
a v) f47 g i))) (49 ((pr3 0 r3 v) (pr2 -1 r2) (pr14 1 r1)) (
plus (minus (times (g l r3 v) f46 g)) (times (g l v r3) f46
g) (times (g l v) f45 g))) (48 ((pr14 -1 r1) (pr2 1 r2) (pr3
0 r3 v)) (plus (minus (times (g l r3 v) f46 g)) (times (g l
v r3) f46 g) (times (g l v) f45 g))) (47 ((pr9 0 r3 v) (pr2
-1 r2) (pr14 1 r1)) (plus (minus (times (g l r3 v) f44 g)) (
times (g l v r3) f44 g) (times (g l v) f43 g))) (46 ((pr14 
-1 r1) (pr2 1 r2) (pr9 0 r3 v)) (plus (minus (times (g l r3
v) f44 g)) (times (g l v r3) f44 g) (times (g l v) f43 g)))
(45 ((pr3 0 r3 v) (pr2 -1 r2) (pr16 1 r1)) (plus (times (g l
a r3 v) f42 g i) (minus (times (g l a v r3) f42 g i)) (minus
(times (g l a v) f41 g i)))) (44 ((pr16 -1 r1) (pr2 1 r2) (
pr3 0 r3 v)) (plus (times (g l a r3 v) f42 g i) (minus (
times (g l a v r3) f42 g i)) (times (g l a v) f41 g i))) (43
((pr9 0 r3 v) (pr2 -1 r2) (pr16 1 r1)) (plus (times (g l a 
r3 v) f40 g i) (minus (times (g l a v r3) f40 g i)) (minus (
times (g l a v) f39 g i)))) (42 ((pr16 -1 r1) (pr2 1 r2) (
pr9 0 r3 v)) (plus (times (g l a r3 v) f40 g i) (minus (
times (g l a v r3) f40 g i)) (times (g l a v) f39 g i))) (41
((pr1 -1 r3) (pr2 -1 r2) (pr17 1 r1 v)) (times (g l a) (cons
r2 v) f38 g)) (40 ((pr17 -1 r1 v) (pr2 1 r2) (pr1 1 r3)) (
times (g l a) (cons r2 v) f38 g)) (39 ((pr1 -1 r3) (pr13 -1
r2) (pr17 1 r1 v)) (times (cons r2 v) f37 g i)) (38 ((pr17 
-1 r1 v) (pr13 1 r2) (pr1 1 r3)) (times (cons r2 v) f37 g i)
) (37 ((pr1 -1 r3) (pr14 -1 r2) (pr17 1 r1 v)) (times (g l a
) (cons r2 v) f36 g)) (36 ((pr17 -1 r1 v) (pr14 1 r2) (pr1 1
r3)) (times (g l a) (cons r2 v) f36 g)) (35 ((pr1 -1 r3) (
pr16 -1 r2) (pr17 1 r1 v)) (times (cons r2 v) f35 g i)) (34
((pr17 -1 r1 v) (pr16 1 r2) (pr1 1 r3)) (times (cons r2 v) 
f35 g i)) (33 ((pr2 -1 r1) (pr2 1 r2) (pr7 0 r3)) (times f34
g)) (32 ((pr1 -1 r3) (pr2 -1 r2) (pr18 1 r1)) (times f33 g))
(31 ((pr18 -1 r1) (pr2 1 r2) (pr1 1 r3)) (times f33 g)) (30
((pr1 -1 r3) (pr16 -1 r2) (pr18 1 r1)) (times (g l a) f32 g
i)) (29 ((pr18 -1 r1) (pr16 1 r2) (pr1 1 r3)) (times (g l a)
f32 g i)) (28 ((pr8 0 r1 v va vb vc) (pr9 0 r2 vd) (pr3 0 r3
vf)) (plus (times (cons r1 vd) (cons r1 vf) (cons r2 v) (
cons r2 va) (cons r2 vb) (cons r2 vc) f31 g) (times (cons r1
vd) (cons r2 v) (cons r2 va) (cons r2 vb) (cons vc vf) f30 g
) (times (cons r1 vf) (cons r2 v) (cons r2 va) (cons r2 vb)
(cons vc vd) f28 g) (times (cons r2 v) (cons r2 va) (cons r2
vb) (cons r2 vc) (cons vd vf) f29 g) (times (cons r2 v) (
cons r2 va) (cons vb vd) (cons vc vf) f27 g))) (27 ((pr8 0 
r1 v va vb vc) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (
cons r2 va) (cons r2 vb) (cons r2 vc) f26 g)) (26 ((pr4 0 r1
v va) (pr9 0 r2 vb) (pr3 0 r3 vc)) (plus (times (cons r1 vb)
(cons r1 vc) (cons r2 v) (cons r2 va) f25 g) (times (cons r1
vb) (cons r2 v) (cons va vc) f24 g) (times (cons r1 vc) (
cons r2 v) (cons va vb) f22 g) (times (cons r2 v) (cons r2 
va) (cons vb vc) f23 g) (times (cons v vb) (cons va vc) f21
g))) (25 ((pr5 0 r1 v va) (pr9 0 r2 vb) (pr3 0 r3 vc)) (plus
(times (cons r1 vb) (cons r1 vc) (cons r2 v) (cons r2 va) 
f20 g) (times (cons r1 vb) (cons r2 v) (cons va vc) f19 g) (
times (cons r1 vc) (cons r2 v) (cons va vb) f17 g) (times (
cons r2 v) (cons r2 va) (cons vb vc) f18 g) (times (cons v 
vb) (cons va vc) f16 g))) (24 ((pr6 0 r1 v va) (pr9 0 r2 vb)
(pr3 0 r3 vc)) (plus (times (cons r1 vb) (cons r1 vc) (cons
r2 v) (cons r2 va) f15 g) (times (cons r1 vb) (cons r2 v) (
cons va vc) f14 g) (times (cons r1 vc) (cons r2 v) (cons va
vb) f12 g) (times (cons r2 v) (cons r2 va) (cons vb vc) f13
g) (times (cons v vb) (cons va vc) f11 g))) (23 ((pr4 0 r1 v
va) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va)
f10 g)) (22 ((pr5 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (times
(cons r2 v) (cons r2 va) f9 g)) (21 ((pr6 0 r1 v va) (pr1 -1
r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va) f8 g)) (20 (
(pr9 0 r1 v) (pr3 0 r2 va) (pr7 0 r3)) (plus (times (cons r1
va) (cons r2 v) f7 g) (times (cons v va) f6 g))) (19 ((pr3 0
r1 v) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) f5 g)) (18
((pr9 0 r1 v) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) f4
g)) (17 ((pr7 0 r1) (pr1 -1 r2) (pr1 1 r3)) (times f3 g)) (
16 ((pr9 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; 
