Loading image file: /usr/lib/reduce/pslbuild/red/reduce.img 
Reduce (Free PSL version, revision 4715),  9-Aug-2018 ...

1: 
(drop model_state midmodel_add massolve equationlist feyrshow midmodel 
particle_table symbole_list_gen model_gen model_gen1 particle_table_gen 
feynman_rules model_input matterintinput symbole_list add_vertices 
physical_parameters physical_parameters_gen)

(sout alldiagramlist outadir fout amp_information option)

(showmore selfenergy justnumber ndiagfile ncolor nout nsy nsym state set_ext_sym
mkpsymboleok)

phys_cut := {}

phyinput_add := {}

ok

hpux

no

model_generator1

model_generator2

model_generator3

*** local variable n_coupling in procedure diagram_generator not used 

diagram_generator

write_diagram_information

write_next_amplitude

write_next_kinematic

generation_kinematics

process_chose23

gen_phyinput

get_x

amplitude_evaluation

model_chose

model_chose1

write_build_model_first

write_build_model

write_start_process

write_start_process1

write_model_information

write_process_information

process_chose

gmprocess

gen_wdlist

add_d0

add_d

d_alg

drop_zero_pa

mk_fpa_diag

gen_fit_parameter_list

find_n_it

gen_resonace_list

get_fasl_type

*** nonlocal use of undeclared variable new_parameter_list in procedure gmhtml 

gmhtml

ghtml

start_fdc
fdc_home:="/home/fdcpwa/fdc-pwa3.0/"$


op_system:='linux$


$

end$

model_state := model_state

"debug/massolve"

"debug/equationlist"

"debug/debug_input"

feyrshow := feynman_rule_show

drop := vertices_drop_list

add_vertices

"system/midmodel"

"system/midmodel_add"

"system/physical_parameters"

"debug/debug_input"

"diagram/"

"diagram/sout"

amp_information

option
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%---FDC---%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %  Now, you have just load the FDC package for Feynman Diagram Calculation  %
 %  The First step, you should choose a Physical model by perform:           %
 %  model_chose(aa);                                                         %
 %  Where aa refer to the file in which the model definition information is  %
 %  stored. The example of it could be found in file:                        %
 %  /home/fdcpwa/fdc-pwa3.0/example/model.def
 %  If the model is already built,  you can go directly to chose             %
 %  Physical process by perform:                                             %
 %  process_chose(aa);                                                       %
 %  Where aa refer to the file in which the process definition information   %
 %  is stored. The example of it could be found in file:                     %
 %  /home/fdcpwa/fdc-pwa3.0/example/process.def
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2:  You choose the The Phenomenology Lagrangian Model for $J/\Psi$ Physics
 in the directory ../model/, which is inputed by 
Jian-Xiong Wang and Rui-Chang Niu on April 15, 2022.
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%---FDC---%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %  The model is already built, so that you can go directly to chose         %
 %  Physical process by perform:                                             %
 %  process_chose(aa);                                                       %
 %  Where aa refer to the file in which the process definition information   %
 %  is stored. The example of it could be found in file:                     %
 %  /home/fdcpwa/fdc-pwa3.0/example/process.def
 %  If you want to rebuild the model with some other choice, You must delete %
 %  all lines in the file ../model/model_state
 %  Then run REDUCE and FDC again.                                           %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3: 
4: 
5: 
6: 6: 6: n_particles:='3$

 
n_total_spin:='6$

 
http:='nil$

 
sub_model:='ok$

 
direct_model:='ok$

 
just_v_type:='nil$

 
nnvertex:='144$

 
pdl_add:='((pr1 0 p1m w1d 0 0 0 0 0 0 1 0 0 0) (pr2 0 p2m 
w2d 0 1 0 0 0 0 1 -1 0 0) (pr3 0 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 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 4 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 2 p13m w13d 0 1 0 0 0
0 1 -1 0 0) (pr14 4 p14m w14d 0 1 0 0 0 0 1 -1 0 0) (pr15 2
p15m w15d 0 1 0 0 0 0 1 -1 0 0) (pr16 6 p16m w16d 0 1 0 0 0
0 1 -1 0 0) (pr17 8 p17m w17d 0 1 0 0 0 0 1 -1 0 0) (pr18 10
p18m w18d 0 1 0 0 0 0 1 -1 0 0) (pr19 0 p19m w19d 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
) (pr14 0) (pr15 0) (pr16 0) (pr17 0) (pr18 0) (pr19 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 p17m w17d p18m w18d
p19m w19d)$

 
platexname_add:='( 
      (pr1  "\pi^0") 
      (pr2  "K^0") 
      (pr3  "\eta^{'}[958]") 
      (pr4  "\omega^h[1900]") 
      (pr5  "f_2[1430]") 
      (pr6  "f^{'}_2[1525]") 
      (pr7  "f_2[1565]") 
      (pr8  "f_2[1640]") 
      (pr9  "f_2[1810]") 
      (pr10  "f_2[1910]") 
      (pr11  "f_2[1950]") 
      (pr12  "f_2[2010]") 
      (pr13  "K^{*}[892]^0") 
      (pr14  "K^{*}_2[1430]") 
      (pr15  "K^{*}[1680]") 
      (pr16  "K^{*}_3[1780]") 
      (pr17  "K^{*}_4[2045]") 
      (pr18  "K^{*}_5[2380]") 
      (pr19  "\eta_{c}[2980]") 
     )$

 
 
 
inputlist_add:='( 
      (p1m  0.13498) 
      (w1d  0.0001) 
      (p2m  0.49767) 
      (w2d  0.0001) 
      (p3m  0.95778) 
      (w3d  0.0002) 
      (p4m  1.9) 
      (w4d  0.15) 
      (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.815) 
      (w9d  0.197) 
      (p10m  1.915) 
      (w10d  0.163) 
      (p11m  1.934) 
      (w11d  0.453) 
      (p12m  2.011) 
      (w12d  0.202) 
      (p13m  0.89166) 
      (w13d  0.051) 
      (p14m  1.4256) 
      (w14d  0.099) 
      (p15m  1.717) 
      (w15d  0.332) 
      (p16m  1.776) 
      (w16d  0.159) 
      (p17m  2.045) 
      (w17d  0.198) 
      (p18m  2.382) 
      (w18d  0.178) 
      (p19m  2.9792) 
      (w19d  0.013) 
     )$

 
 
pa_decay_list:='((pr13 (24 (pr2 pr1) 2) (23 (pr2 pr1) 1)) (
pr14 (70 (pr13 pr1) 2) (69 (pr13 pr1) 1) (37 (pr2 pr1) 2) (
36 (pr2 pr1) 1)) (pr5 (86 (pr13 pr2) 2) (85 (pr13 pr2) 2) (
63 (pr1 pr1) 2) (62 (pr2 pr2) 2)) (pr6 (84 (pr13 pr2) 2) (83
(pr13 pr2) 2) (60 (pr1 pr1) 2) (59 (pr2 pr2) 2)) (pr7 (82 (
pr13 pr2) 2) (81 (pr13 pr2) 2) (57 (pr1 pr1) 2) (56 (pr2 pr2
) 2)) (pr8 (80 (pr13 pr2) 2) (79 (pr13 pr2) 2) (54 (pr1 pr1)
2) (53 (pr2 pr2) 2)) (pr15 (66 (pr14 pr1) 1) (65 (pr14 pr1)
2) (30 (pr13 pr1) 2) (29 (pr13 pr1) 1) (20 (pr2 pr1) 2) (19
(pr2 pr1) 1) (18 (pr3 pr2) 2) (17 (pr3 pr2) 1)) (pr16 (116 (
pr14 pr1) 2) (115 (pr14 pr1) 1) (108 (pr13 pr1) 2) (107 (
pr13 pr1) 1) (92 (pr2 pr1) 2) (91 (pr2 pr1) 1) (90 (pr3 pr2)
2) (89 (pr3 pr2) 1)) (pr9 (96 (pr13 pr13) 2) (78 (pr13 pr2)
2) (77 (pr13 pr2) 2) (51 (pr1 pr1) 2) (50 (pr2 pr2) 2)) (pr4
(64 (pr13 pr13) 2) (33 (pr13 pr2) 2) (32 (pr13 pr2) 2) (26 (
pr1 pr1) 2)) (pr10 (95 (pr13 pr13) 2) (76 (pr13 pr2) 2) (75
(pr13 pr2) 2) (48 (pr1 pr1) 2) (47 (pr2 pr2) 2)) (pr11 (101
(pr14 pr2) 2) (100 (pr14 pr2) 2) (94 (pr13 pr13) 2) (74 (
pr13 pr2) 2) (73 (pr13 pr2) 2) (45 (pr1 pr1) 2) (44 (pr2 pr2
) 2) (43 (pr3 pr3) 2)) (pr12 (99 (pr14 pr2) 2) (98 (pr14 pr2
) 2) (93 (pr13 pr13) 2) (72 (pr13 pr2) 2) (71 (pr13 pr2) 2)
(41 (pr1 pr1) 2) (40 (pr2 pr2) 2) (39 (pr3 pr3) 2)) (pr17 (
138 (pr5 pr2) 2) (137 (pr5 pr2) 1) (136 (pr6 pr2) 2) (135 (
pr6 pr2) 1) (134 (pr7 pr2) 2) (133 (pr7 pr2) 1) (132 (pr14 
pr1) 2) (131 (pr14 pr1) 1) (124 (pr13 pr1) 2) (123 (pr13 pr1
) 1) (122 (pr13 pr3) 2) (121 (pr13 pr3) 1) (118 (pr15 pr1) 2
) (117 (pr15 pr1) 1) (114 (pr2 pr1) 2) (113 (pr2 pr1) 1) (
112 (pr3 pr2) 2) (111 (pr3 pr2) 1)) (pr18 (144 (pr13 pr1) 2)
(143 (pr13 pr1) 1) (142 (pr13 pr3) 2) (141 (pr13 pr3) 1) (
140 (pr15 pr1) 2) (139 (pr15 pr1) 1) (130 (pr2 pr1) 2) (129
(pr2 pr1) 1) (128 (pr3 pr2) 2) (127 (pr3 pr2) 1)) (pr19 (126
(pr18 pr2) 2) (125 (pr18 pr2) 2) (120 (pr17 pr13) 2) (119 (
pr17 pr13) 2) (110 (pr17 pr2) 2) (109 (pr17 pr2) 2) (106 (
pr16 pr13) 2) (105 (pr16 pr13) 2) (104 (pr5 pr5) 2) (103 (
pr6 pr5) 2) (102 (pr7 pr5) 2) (97 (pr14 pr14) 2) (88 (pr16 
pr2) 2) (87 (pr16 pr2) 2) (68 (pr14 pr13) 2) (67 (pr14 pr13)
2) (61 (pr5 pr3) 2) (58 (pr6 pr3) 2) (55 (pr7 pr3) 2) (52 (
pr8 pr3) 2) (49 (pr9 pr3) 2) (46 (pr10 pr3) 2) (42 (pr11 pr3
) 2) (38 (pr12 pr3) 2) (35 (pr14 pr2) 2) (34 (pr14 pr2) 2) (
31 (pr13 pr13) 2) (28 (pr15 pr13) 2) (27 (pr15 pr13) 2) (25
(pr4 pr3) 2) (22 (pr13 pr2) 2) (21 (pr13 pr2) 2) (16 (pr15 
pr2) 2) (15 (pr15 pr2) 2)))$

 
htionvvlist_add:='((144 ((pr1 0 r3) (pr13 -1 r2 vf) (pr18 1
r1 v va vb vc vd)) (times (cons r2 va) (cons r2 vb) (cons r2
vc) (cons r2 vd) (eps r1 r2 v vf) f113 g i)) (143 ((pr18 -1
r1 v va vb vc vd) (pr13 1 r2 vf) (pr1 0 r3)) (minus (times (
cons r2 va) (cons r2 vb) (cons r2 vc) (cons r2 vd) (eps r1 
r2 v vf) f113 g i))) (142 ((pr3 0 r3) (pr13 -1 r2 vf) (pr18
1 r1 v va vb vc vd)) (minus (times (cons r2 va) (cons r2 vb)
(cons r2 vc) (cons r2 vd) (eps r1 r2 v vf) f112 g))) (141 ((
pr18 -1 r1 v va vb vc vd) (pr13 1 r2 vf) (pr3 0 r3)) (minus
(times (cons r2 va) (cons r2 vb) (cons r2 vc) (cons r2 vd) (
eps r1 r2 v vf) f112 g))) (140 ((pr1 0 r3) (pr15 -1 r2 vf) (
pr18 1 r1 v va vb vc vd)) (times (cons r2 va) (cons r2 vb) (
cons r2 vc) (cons r2 vd) (eps r1 r2 v vf) f111 g i)) (139 ((
pr18 -1 r1 v va vb vc vd) (pr15 1 r2 vf) (pr1 0 r3)) (minus
(times (cons r2 va) (cons r2 vb) (cons r2 vc) (cons r2 vd) (
eps r1 r2 v vf) f111 g i))) (138 ((pr2 -1 r3) (pr5 0 r2 vd 
vf) (pr17 1 r1 v va vb vc)) (plus (minus (times (cons r1 vf)
(cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f110
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f109 g)))) (137 ((pr17 -1 r1 v va vb vc) (
pr5 0 r2 vd vf) (pr2 1 r3)) (plus (minus (times (cons r1 vf)
(cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f110
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f109 g)))) (136 ((pr2 -1 r3) (pr6 0 r2 vd vf
) (pr17 1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (
cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f108
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f107 g)))) (135 ((pr17 -1 r1 v va vb vc) (
pr6 0 r2 vd vf) (pr2 1 r3)) (plus (minus (times (cons r1 vf)
(cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f108
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f107 g)))) (134 ((pr2 -1 r3) (pr7 0 r2 vd vf
) (pr17 1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (
cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f106
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f105 g)))) (133 ((pr17 -1 r1 v va vb vc) (
pr7 0 r2 vd vf) (pr2 1 r3)) (plus (minus (times (cons r1 vf)
(cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f106
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f105 g)))) (132 ((pr1 0 r3) (pr14 -1 r2 vd 
vf) (pr17 1 r1 v va vb vc)) (plus (times (cons r1 vf) (cons
r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f104 g i)
(times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1 r2 v 
vd) f103 g i))) (131 ((pr17 -1 r1 v va vb vc) (pr14 1 r2 vd
vf) (pr1 0 r3)) (plus (minus (times (cons r1 vf) (cons r2 va
) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f104 g i)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f103 g i)))) (130 ((pr1 0 r3) (pr2 -1 r2) (pr18 1 
r1 v va vb vc vd)) (minus (times (cons r2 v) (cons r2 va) (
cons r2 vb) (cons r2 vc) (cons r2 vd) f102 g))) (129 ((pr18
-1 r1 v va vb vc vd) (pr2 1 r2) (pr1 0 r3)) (times (cons r2
v) (cons r2 va) (cons r2 vb) (cons r2 vc) (cons r2 vd) f102
g)) (128 ((pr2 -1 r3) (pr3 0 r2) (pr18 1 r1 v va vb vc vd))
(times (cons r2 v) (cons r2 va) (cons r2 vb) (cons r2 vc) (
cons r2 vd) f101 g i)) (127 ((pr18 -1 r1 v va vb vc vd) (pr3
0 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va) (cons r2 
vb) (cons r2 vc) (cons r2 vd) f101 g i)) (126 ((pr2 -1 r3) (
pr19 0 r2) (pr18 1 r1 v va vb vc vd)) (times (cons r2 v) (
cons r2 va) (cons r2 vb) (cons r2 vc) (cons r2 vd) f100 g i)
) (125 ((pr18 -1 r1 v va vb vc vd) (pr19 0 r2) (pr2 1 r3)) (
times (cons r2 v) (cons r2 va) (cons r2 vb) (cons r2 vc) (
cons r2 vd) f100 g i)) (124 ((pr1 0 r3) (pr13 -1 r2 vd) (
pr17 1 r1 v va vb vc)) (times (cons r2 va) (cons r2 vb) (
cons r2 vc) (eps r1 r2 v vd) f99 g)) (123 ((pr17 -1 r1 v va
vb vc) (pr13 1 r2 vd) (pr1 0 r3)) (minus (times (cons r2 va)
(cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f99 g))) (122 ((
pr3 0 r3) (pr13 -1 r2 vd) (pr17 1 r1 v va vb vc)) (minus (
times (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd
) f98 g i))) (121 ((pr17 -1 r1 v va vb vc) (pr13 1 r2 vd) (
pr3 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (cons r2
vc) (eps r1 r2 v vd) f98 g i))) (120 ((pr19 0 r3) (pr13 -1 
r2 vd) (pr17 1 r1 v va vb vc)) (minus (times (cons r2 va) (
cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f97 g i))) (119 ((
pr17 -1 r1 v va vb vc) (pr13 1 r2 vd) (pr19 0 r3)) (minus (
times (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd
) f97 g i))) (118 ((pr1 0 r3) (pr15 -1 r2 vd) (pr17 1 r1 v 
va vb vc)) (times (cons r2 va) (cons r2 vb) (cons r2 vc) (
eps r1 r2 v vd) f96 g)) (117 ((pr17 -1 r1 v va vb vc) (pr15
1 r2 vd) (pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb)
(cons r2 vc) (eps r1 r2 v vd) f96 g))) (116 ((pr1 0 r3) (
pr14 -1 r2 vc vd) (pr16 1 r1 v va vb)) (plus (times (cons r1
vd) (cons r2 va) (cons r2 vb) (eps r1 r2 v vc) f95 g) (times
(cons r2 va) (cons vb vd) (eps r1 r2 v vc) f94 g))) (115 ((
pr16 -1 r1 v va vb) (pr14 1 r2 vc vd) (pr1 0 r3)) (plus (
minus (times (cons r1 vd) (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f95 g)) (minus (times (cons r2 va) (cons vb vd) (
eps r1 r2 v vc) f94 g)))) (114 ((pr1 0 r3) (pr2 -1 r2) (pr17
1 r1 v va vb vc)) (minus (times (cons r2 v) (cons r2 va) (
cons r2 vb) (cons r2 vc) f93 g i))) (113 ((pr17 -1 r1 v va 
vb vc) (pr2 1 r2) (pr1 0 r3)) (times (cons r2 v) (cons r2 va
) (cons r2 vb) (cons r2 vc) f93 g i)) (112 ((pr2 -1 r3) (pr3
0 r2) (pr17 1 r1 v va vb vc)) (times (cons r2 v) (cons r2 va
) (cons r2 vb) (cons r2 vc) f92 g)) (111 ((pr17 -1 r1 v va 
vb vc) (pr3 0 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va
) (cons r2 vb) (cons r2 vc) f92 g)) (110 ((pr2 -1 r3) (pr19
0 r2) (pr17 1 r1 v va vb vc)) (times (cons r2 v) (cons r2 va
) (cons r2 vb) (cons r2 vc) f91 g)) (109 ((pr17 -1 r1 v va 
vb vc) (pr19 0 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 
va) (cons r2 vb) (cons r2 vc) f91 g)) (108 ((pr1 0 r3) (pr13
-1 r2 vc) (pr16 1 r1 v va vb)) (times (cons r2 va) (cons r2
vb) (eps r1 r2 v vc) f90 g i)) (107 ((pr16 -1 r1 v va vb) (
pr13 1 r2 vc) (pr1 0 r3)) (minus (times (cons r2 va) (cons 
r2 vb) (eps r1 r2 v vc) f90 g i))) (106 ((pr19 0 r3) (pr13 
-1 r2 vc) (pr16 1 r1 v va vb)) (minus (times (cons r2 va) (
cons r2 vb) (eps r1 r2 v vc) f89 g))) (105 ((pr16 -1 r1 v va
vb) (pr13 1 r2 vc) (pr19 0 r3)) (minus (times (cons r2 va) (
cons r2 vb) (eps r1 r2 v vc) f89 g))) (104 ((pr5 0 r1 v va)
(pr5 0 r2 vb vc) (pr19 0 r3)) (plus (minus (times (cons r1 
vc) (cons r2 va) (eps r1 r2 v vb) f88 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f87 g)))) (103 ((pr6 0 r1 v va)
(pr5 0 r2 vb vc) (pr19 0 r3)) (plus (minus (times (cons r1 
vc) (cons r2 va) (eps r1 r2 v vb) f86 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f85 g)))) (102 ((pr7 0 r1 v va)
(pr5 0 r2 vb vc) (pr19 0 r3)) (plus (minus (times (cons r1 
vc) (cons r2 va) (eps r1 r2 v vb) f84 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f83 g)))) (101 ((pr2 -1 r3) (
pr11 0 r2 vb vc) (pr14 1 r1 v va)) (plus (minus (times (cons
r1 vc) (cons r2 va) (eps r1 r2 v vb) f82 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f81 g)))) (100 ((pr14 -1 r1 v 
va) (pr11 0 r2 vb vc) (pr2 1 r3)) (plus (minus (times (cons
r1 vc) (cons r2 va) (eps r1 r2 v vb) f82 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f81 g)))) (99 ((pr2 -1 r3) (
pr12 0 r2 vb vc) (pr14 1 r1 v va)) (plus (minus (times (cons
r1 vc) (cons r2 va) (eps r1 r2 v vb) f80 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f79 g)))) (98 ((pr14 -1 r1 v va
) (pr12 0 r2 vb vc) (pr2 1 r3)) (plus (minus (times (cons r1
vc) (cons r2 va) (eps r1 r2 v vb) f80 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f79 g)))) (97 ((pr14 -1 r1 v va
) (pr14 1 r2 vb vc) (pr19 0 r3)) (plus (minus (times (cons 
r1 vc) (cons r2 va) (eps r1 r2 v vb) f78 g)) (minus (times (
cons va vc) (eps r1 r2 v vb) f77 g)))) (96 ((pr9 0 r1 v va)
(pr13 -1 r2 vb) (pr13 1 r3 vc)) (plus (times (cons r1 vb) (
cons r1 vc) (cons r2 v) (cons r2 va) f76 g) (times (cons r1
vb) (cons r2 v) (cons va vc) f75 g) (minus (times (cons r1 
vc) (cons r2 v) (cons r2 va) (cons r3 vb) f76 g)) (times (
cons r1 vc) (cons r2 v) (cons va vb) f73 g) (minus (times (
cons r1 vc) (cons r2 v) (cons va vb) f75 g)) (times (cons r2
v) (cons r2 va) (cons vb vc) f74 g) (times (cons r2 v) (cons
r3 vb) (cons va vc) f73 g) (times (cons v vb) (cons va vc) 
f72 g) (times (cons v vc) (cons va vb) f72 g))) (95 ((pr10 0
r1 v va) (pr13 -1 r2 vb) (pr13 1 r3 vc)) (plus (times (cons
r1 vb) (cons r1 vc) (cons r2 v) (cons r2 va) f71 g) (times (
cons r1 vb) (cons r2 v) (cons va vc) f70 g) (minus (times (
cons r1 vc) (cons r2 v) (cons r2 va) (cons r3 vb) f71 g)) (
times (cons r1 vc) (cons r2 v) (cons va vb) f68 g) (minus (
times (cons r1 vc) (cons r2 v) (cons va vb) f70 g)) (times (
cons r2 v) (cons r2 va) (cons vb vc) f69 g) (times (cons r2
v) (cons r3 vb) (cons va vc) f68 g) (times (cons v vb) (cons
va vc) f67 g) (times (cons v vc) (cons va vb) f67 g))) (94 (
(pr11 0 r1 v va) (pr13 -1 r2 vb) (pr13 1 r3 vc)) (plus (
times (cons r1 vb) (cons r1 vc) (cons r2 v) (cons r2 va) f66
g) (times (cons r1 vb) (cons r2 v) (cons va vc) f65 g) (
minus (times (cons r1 vc) (cons r2 v) (cons r2 va) (cons r3
vb) f66 g)) (times (cons r1 vc) (cons r2 v) (cons va vb) f63
g) (minus (times (cons r1 vc) (cons r2 v) (cons va vb) f65 g
)) (times (cons r2 v) (cons r2 va) (cons vb vc) f64 g) (
times (cons r2 v) (cons r3 vb) (cons va vc) f63 g) (times (
cons v vb) (cons va vc) f62 g) (times (cons v vc) (cons va 
vb) f62 g))) (93 ((pr12 0 r1 v va) (pr13 -1 r2 vb) (pr13 1 
r3 vc)) (plus (times (cons r1 vb) (cons r1 vc) (cons r2 v) (
cons r2 va) f61 g) (times (cons r1 vb) (cons r2 v) (cons va
vc) f60 g) (minus (times (cons r1 vc) (cons r2 v) (cons r2 
va) (cons r3 vb) f61 g)) (times (cons r1 vc) (cons r2 v) (
cons va vb) f58 g) (minus (times (cons r1 vc) (cons r2 v) (
cons va vb) f60 g)) (times (cons r2 v) (cons r2 va) (cons vb
vc) f59 g) (times (cons r2 v) (cons r3 vb) (cons va vc) f58
g) (times (cons v vb) (cons va vc) f57 g) (times (cons v vc)
(cons va vb) f57 g))) (92 ((pr1 0 r3) (pr2 -1 r2) (pr16 1 r1
v va vb)) (minus (times (cons r2 v) (cons r2 va) (cons r2 vb
) f56 g))) (91 ((pr16 -1 r1 v va vb) (pr2 1 r2) (pr1 0 r3))
(times (cons r2 v) (cons r2 va) (cons r2 vb) f56 g)) (90 ((
pr2 -1 r3) (pr3 0 r2) (pr16 1 r1 v va vb)) (times (cons r2 v
) (cons r2 va) (cons r2 vb) f55 g i)) (89 ((pr16 -1 r1 v va
vb) (pr3 0 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va) (
cons r2 vb) f55 g i)) (88 ((pr2 -1 r3) (pr19 0 r2) (pr16 1 
r1 v va vb)) (times (cons r2 v) (cons r2 va) (cons r2 vb) 
f54 g i)) (87 ((pr16 -1 r1 v va vb) (pr19 0 r2) (pr2 1 r3))
(times (cons r2 v) (cons r2 va) (cons r2 vb) f54 g i)) (86 (
(pr2 -1 r3) (pr13 1 r2 vb) (pr5 0 r1 v va)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f53 g i))) (85 ((pr5 0 r1 v va)
(pr13 -1 r2 vb) (pr2 1 r3)) (minus (times (cons r2 va) (eps
r1 r2 v vb) f53 g i))) (84 ((pr2 -1 r3) (pr13 1 r2 vb) (pr6
0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v vb) f52
g i))) (83 ((pr6 0 r1 v va) (pr13 -1 r2 vb) (pr2 1 r3)) (
minus (times (cons r2 va) (eps r1 r2 v vb) f52 g i))) (82 ((
pr2 -1 r3) (pr13 1 r2 vb) (pr7 0 r1 v va)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f51 g i))) (81 ((pr7 0 r1 v va)
(pr13 -1 r2 vb) (pr2 1 r3)) (minus (times (cons r2 va) (eps
r1 r2 v vb) f51 g i))) (80 ((pr2 -1 r3) (pr13 1 r2 vb) (pr8
0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v vb) f50
g i))) (79 ((pr8 0 r1 v va) (pr13 -1 r2 vb) (pr2 1 r3)) (
minus (times (cons r2 va) (eps r1 r2 v vb) f50 g i))) (78 ((
pr2 -1 r3) (pr13 1 r2 vb) (pr9 0 r1 v va)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f49 g i))) (77 ((pr9 0 r1 v va)
(pr13 -1 r2 vb) (pr2 1 r3)) (minus (times (cons r2 va) (eps
r1 r2 v vb) f49 g i))) (76 ((pr2 -1 r3) (pr13 1 r2 vb) (pr10
0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v vb) f48
g i))) (75 ((pr10 0 r1 v va) (pr13 -1 r2 vb) (pr2 1 r3)) (
minus (times (cons r2 va) (eps r1 r2 v vb) f48 g i))) (74 ((
pr2 -1 r3) (pr13 1 r2 vb) (pr11 0 r1 v va)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f47 g i))) (73 ((pr11 0 r1 v va
) (pr13 -1 r2 vb) (pr2 1 r3)) (minus (times (cons r2 va) (
eps r1 r2 v vb) f47 g i))) (72 ((pr2 -1 r3) (pr13 1 r2 vb) (
pr12 0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v vb)
f46 g i))) (71 ((pr12 0 r1 v va) (pr13 -1 r2 vb) (pr2 1 r3))
(minus (times (cons r2 va) (eps r1 r2 v vb) f46 g i))) (70 (
(pr1 0 r3) (pr13 -1 r2 vb) (pr14 1 r1 v va)) (times (cons r2
va) (eps r1 r2 v vb) f45 g)) (69 ((pr14 -1 r1 v va) (pr13 1
r2 vb) (pr1 0 r3)) (minus (times (cons r2 va) (eps r1 r2 v 
vb) f45 g))) (68 ((pr19 0 r3) (pr13 -1 r2 vb) (pr14 1 r1 v 
va)) (minus (times (cons r2 va) (eps r1 r2 v vb) f44 g i)))
(67 ((pr14 -1 r1 v va) (pr13 1 r2 vb) (pr19 0 r3)) (minus (
times (cons r2 va) (eps r1 r2 v vb) f44 g i))) (66 ((pr1 0 
r3) (pr15 -1 r2 vb) (pr14 1 r1 v va)) (times (cons r2 va) (
eps r1 r2 v vb) f43 g)) (65 ((pr14 -1 r1 v va) (pr15 1 r2 vb
) (pr1 0 r3)) (minus (times (cons r2 va) (eps r1 r2 v vb) 
f43 g))) (64 ((pr13 -1 r1 v) (pr13 1 r2 va) (pr4 0 r3 vb)) (
plus (times (cons r1 va) (cons r1 vb) (cons r2 v) f42 g i) (
times (cons r1 va) (cons v vb) f41 g i) (times (cons r1 vb)
(cons r2 v) (cons r3 va) f42 g i) (times (cons r2 v) (cons 
va vb) f40 g i) (times (cons r2 v) (cons va vb) f41 g i) (
minus (times (cons r3 va) (cons v vb) f40 g i)))) (63 ((pr5
0 r1 v va) (pr1 0 r2) (pr1 0 r3)) (times (cons r2 v) (cons 
r2 va) f39 g)) (62 ((pr5 0 r1 v va) (pr2 -1 r2) (pr2 1 r3))
(times (cons r2 v) (cons r2 va) f38 g)) (61 ((pr5 0 r1 v va)
(pr19 0 r2) (pr3 0 r3)) (times (cons r2 v) (cons r2 va) f37
g)) (60 ((pr6 0 r1 v va) (pr1 0 r2) (pr1 0 r3)) (times (cons
r2 v) (cons r2 va) f36 g)) (59 ((pr6 0 r1 v va) (pr2 -1 r2)
(pr2 1 r3)) (times (cons r2 v) (cons r2 va) f35 g)) (58 ((
pr6 0 r1 v va) (pr19 0 r2) (pr3 0 r3)) (times (cons r2 v) (
cons r2 va) f34 g)) (57 ((pr7 0 r1 v va) (pr1 0 r2) (pr1 0 
r3)) (times (cons r2 v) (cons r2 va) f33 g)) (56 ((pr7 0 r1
v va) (pr2 -1 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va
) f32 g)) (55 ((pr7 0 r1 v va) (pr19 0 r2) (pr3 0 r3)) (
times (cons r2 v) (cons r2 va) f31 g)) (54 ((pr8 0 r1 v va)
(pr1 0 r2) (pr1 0 r3)) (times (cons r2 v) (cons r2 va) f30 g
)) (53 ((pr8 0 r1 v va) (pr2 -1 r2) (pr2 1 r3)) (times (cons
r2 v) (cons r2 va) f29 g)) (52 ((pr8 0 r1 v va) (pr19 0 r2)
(pr3 0 r3)) (times (cons r2 v) (cons r2 va) f28 g)) (51 ((
pr9 0 r1 v va) (pr1 0 r2) (pr1 0 r3)) (times (cons r2 v) (
cons r2 va) f27 g)) (50 ((pr9 0 r1 v va) (pr2 -1 r2) (pr2 1
r3)) (times (cons r2 v) (cons r2 va) f26 g)) (49 ((pr9 0 r1
v va) (pr19 0 r2) (pr3 0 r3)) (times (cons r2 v) (cons r2 va
) f25 g)) (48 ((pr10 0 r1 v va) (pr1 0 r2) (pr1 0 r3)) (
times (cons r2 v) (cons r2 va) f24 g)) (47 ((pr10 0 r1 v va)
(pr2 -1 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va) f23
g)) (46 ((pr10 0 r1 v va) (pr19 0 r2) (pr3 0 r3)) (times (
cons r2 v) (cons r2 va) f22 g)) (45 ((pr11 0 r1 v va) (pr1 0
r2) (pr1 0 r3)) (times (cons r2 v) (cons r2 va) f21 g)) (44
((pr11 0 r1 v va) (pr2 -1 r2) (pr2 1 r3)) (times (cons r2 v)
(cons r2 va) f20 g)) (43 ((pr11 0 r1 v va) (pr3 0 r2) (pr3 0
r3)) (times (cons r2 v) (cons r2 va) f19 g)) (42 ((pr11 0 r1
v va) (pr19 0 r2) (pr3 0 r3)) (times (cons r2 v) (cons r2 va
) f18 g)) (41 ((pr12 0 r1 v va) (pr1 0 r2) (pr1 0 r3)) (
times (cons r2 v) (cons r2 va) f17 g)) (40 ((pr12 0 r1 v va)
(pr2 -1 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va) f16
g)) (39 ((pr12 0 r1 v va) (pr3 0 r2) (pr3 0 r3)) (times (
cons r2 v) (cons r2 va) f15 g)) (38 ((pr12 0 r1 v va) (pr19
0 r2) (pr3 0 r3)) (times (cons r2 v) (cons r2 va) f14 g)) (
37 ((pr1 0 r3) (pr2 -1 r2) (pr14 1 r1 v va)) (minus (times (
cons r2 v) (cons r2 va) f13 g i))) (36 ((pr14 -1 r1 v va) (
pr2 1 r2) (pr1 0 r3)) (times (cons r2 v) (cons r2 va) f13 g
i)) (35 ((pr2 -1 r3) (pr19 0 r2) (pr14 1 r1 v va)) (times (
cons r2 v) (cons r2 va) f12 g)) (34 ((pr14 -1 r1 v va) (pr19
0 r2) (pr2 1 r3)) (times (cons r2 v) (cons r2 va) f12 g)) (
33 ((pr2 -1 r3) (pr4 0 r2 va) (pr13 1 r1 v)) (minus (times (
eps r1 r2 v va) f11 g))) (32 ((pr13 -1 r1 v) (pr4 0 r2 va) (
pr2 1 r3)) (minus (times (eps r1 r2 v va) f11 g))) (31 ((
pr13 -1 r1 v) (pr13 1 r2 va) (pr19 0 r3)) (minus (times (eps
r1 r2 v va) f10 g))) (30 ((pr1 0 r3) (pr13 -1 r2 va) (pr15 1
r1 v)) (times (eps r1 r2 v va) f9 g i)) (29 ((pr15 -1 r1 v)
(pr13 1 r2 va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f9
g i))) (28 ((pr19 0 r3) (pr13 -1 r2 va) (pr15 1 r1 v)) (
minus (times (eps r1 r2 v va) f8 g))) (27 ((pr15 -1 r1 v) (
pr13 1 r2 va) (pr19 0 r3)) (minus (times (eps r1 r2 v va) f8
g))) (26 ((pr4 0 r1 v) (pr1 0 r2) (pr1 0 r3)) (times (cons 
r2 v) f7 g i)) (25 ((pr4 0 r1 v) (pr19 0 r2) (pr3 0 r3)) (
times (cons r2 v) f6 g i)) (24 ((pr1 0 r3) (pr2 -1 r2) (pr13
1 r1 v)) (minus (times (cons r2 v) f5 g))) (23 ((pr13 -1 r1
v) (pr2 1 r2) (pr1 0 r3)) (times (cons r2 v) f5 g)) (22 ((
pr2 -1 r3) (pr19 0 r2) (pr13 1 r1 v)) (times (cons r2 v) f4
g i)) (21 ((pr13 -1 r1 v) (pr19 0 r2) (pr2 1 r3)) (times (
cons r2 v) f4 g i)) (20 ((pr1 0 r3) (pr2 -1 r2) (pr15 1 r1 v
)) (minus (times (cons r2 v) f3 g))) (19 ((pr15 -1 r1 v) (
pr2 1 r2) (pr1 0 r3)) (times (cons r2 v) f3 g)) (18 ((pr2 -1
r3) (pr3 0 r2) (pr15 1 r1 v)) (times (cons r2 v) f2 g i)) (
17 ((pr15 -1 r1 v) (pr3 0 r2) (pr2 1 r3)) (times (cons r2 v)
f2 g i)) (16 ((pr2 -1 r3) (pr19 0 r2) (pr15 1 r1 v)) (times
(cons r2 v) f1 g i)) (15 ((pr15 -1 r1 v) (pr19 0 r2) (pr2 1
r3)) (times (cons r2 v) f1 g i)))$

 
vertices_add_comment:='"
In this model, We have added some effective 
vertices, to do some calculation related to some resonance states
"$

 
end;

7: 
matterinteraction:=el(-1)*er(1)*hig(1)*c(1,4)
+hig(-1)*er(-1)*el(1)*c(-1,4)
+hig(-1)*hig(1)*c(0,1)
+hig(-1)*hig(1)*hig(-1)*hig(1)*c(0,5)
+hig(-1)*q1dr(-1)*q1l(1)*c(-1,3)
+higtid(-1)*q1ur(-1)*q1l(1)*c(-1,2)
+q1l(-1)*q1dr(1)*hig(1)*c(1,3)
+q1l(-1)*q1ur(1)*higtid(1)*c(1,2)$


*** el declared operator 

*** er declared operator 

*** hig declared operator 

*** q1dr declared operator 

*** q1l declared operator 

*** higtid declared operator 

*** q1ur declared operator 


 $


gl_s_list:='(    )$


 $



end$

8: ;

end;

9: platexname:='(
        (p "\gamma") 
        (pa1 "\gamma~Z^0") 
        (z "Z") 
        (w w) 
        (gs g) 
        (nue "{\nu}_e") 
        (numu "{\nu}_{\mu}") 
        (nut "{\nu}_{\tau}") 
        (ef e)
        (mu "\mu")
        (tau "\tau")
        (qu u) 
        (qc c) 
        (qt t) 
        (qd d) 
        (qs s) 
        (qb b) 
        (h0 "h") 
        (gsg "{g_g}")
      );


((p "\gamma") (pa1 "\gamma~Z^0") (z "Z") (w w) (gs g) (nue "{\nu}_e") (numu 
"{\nu}_{\mu}") (nut "{\nu}_{\tau}") (ef e) (mu "\mu") (tau "\tau") (qu u) (qc c)
(qt t) (qd d) (qs s) (qb b) (h0 "h") (gsg "{g_g}"))

 
ptypelist:='( 
        (psi 1) 
        (psi3 1) 
        (psi5 1) 
        (b1 (gs 43) 42) 
        (b0 (gsg 121) (h0 121) 41)
      );


((psi 1) (psi3 1) (psi5 1) (b1 (gs 43) 42) (b0 (gsg 121) (h0 121) 41))
 

latex_s_list:='(
        (theta "{\theta_w}")  
        (xi "{\xi}")  
        (g3 "{g_s}") 
        (s0 "S_0")
        (s1 "V_1")
        (s2 "V_2")
        (s3 "V_3")
        (s4 "V_4")
        (s5 "V_5")
        (s6 "V_6")
        (s7 "V_7")
        (s8 "V_8")
        (sf "F")
        (sf1 "F_{1/2}")
        (sf3 "F_{3/2}")
        (sf5 "F_{5/2}")
        (sf7 "F_{7/2}")
        (sf9 "F_{9/2}")
        (sf11 "F_{11/2}")
        (sf13 "F_{13/2}")
        (sf15 "F_{15/2}")
     );


((theta "{\theta_w}") (xi "{\xi}") (g3 "{g_s}") (s0 "S_0") (s1 "V_1") (s2 "V_2")
(s3 "V_3") (s4 "V_4") (s5 "V_5") (s6 "V_6") (s7 "V_7") (s8 "V_8") (sf "F") (sf1
"F_{1/2}") (sf3 "F_{3/2}") (sf5 "F_{5/2}") (sf7 "F_{7/2}") (sf9 "F_{9/2}") (sf11
"F_{11/2}") (sf13 "F_{13/2}") (sf15 "F_{15/2}"))
 
end;

10: 
     *  ......The latex file model.tex was generated 
     *        Please latex it 
11: 
Time: 10 ms

12: 
Quitting
