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: 
3: 
4:  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.
***** End-of-file read in file ../model/model_state
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%---FDC---%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %  To build the model, the first step is to perform:                        %
 %     model_generator1();                                                   %
 %  If you have already run it or you are sure matter interaction terms in   %
 %  the file ../model/matterintinput is ok.
 %  Be carefull about the matter interaction terms, Then you can perform     %
 %     model_generator2();                                                   %
 %  if it already be run and you are sure that the vertices which you want   %
 %  to drop are listed in the file  %  ../model/vertices_drop_list
 %  Then you can go head to perform:                                         %
 %     model_generator3();                                                   %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5: 
Time: 0 ms

6: *** local variable n in procedure tpdm not used 
*** es already defined as operator 
*** local variable b in procedure sin_solve not used 
*** find_charge called with 6 instead of 3 arguments 
*** local variable n in procedure find_charge1 not used 
*** local variable a in procedure project_color_factor3 not used 
*** b declared operator 
*** a declared operator 
*** gs declared operator 
*** h0 declared operator 
*** nue declared operator 
*** ef declared operator 
*** qu declared operator 
*** qd declared operator 
*** w declared operator 
*** z declared operator 
*** p declared operator 

*** qu already defined as operator 

*** qd already defined as operator 

*** se1 already defined as operator 

*** s1 already defined as operator 

*** qd already defined as operator 

*** qu already defined as operator 

*** ef already defined as operator 

*** nue already defined as operator 

*** h0 already defined as operator 

*** gs already defined as operator 

*** w already defined as operator 

*** z already defined as operator 

*** p already defined as operator 

hfixd==pd(gs(ic10,va),va)

 

      - pd(gs(ic10,va),va)*pd(gs(ic10,v),v)
aa==----------------------------------------
                       2

 

hg==gsg(-1,)*gsg(1,ic11)*pd(gs(ic12,v),v)*cg(ic10,ic11,ic12,3)*g3

     + gsg(-1,)*pd(gsg(1,ic11),v)*cg(ic10,ic11,ic12,3)*gs(ic12,v)*g3

     - gsg(-1,)*pd(pd(gsg(1,ic10),v),v)

 

                 2      2                 2
hh==(6*cos(theta) *h0(0) *w(1,v)*w(-1,v)*g

                    2      2       2  2
      + 3*cos(theta) *h0(0) *z(0,v) *g

                     2                       2
      + 12*cos(theta) *h0(0)*w(1,v)*w(-1,v)*g *v0

                    2             2  2
      + 6*cos(theta) *h0(0)*z(0,v) *g *v0

                     2       2         2  2
      + 12*cos(theta) *w(1,v) *w(-1,va) *g

                     2                                  2
      - 12*cos(theta) *w(1,v)*w(1,va)*w(-1,v)*w(-1,va)*g

                     2                       2  2
      - 24*cos(theta) *w(1,v)*w(-1,v)*z(0,va) *g

                    2                 2   2
      + 6*cos(theta) *w(1,v)*w(-1,v)*g *v0

                     2                                 2
      + 24*cos(theta) *w(1,v)*w(-1,va)*z(0,v)*z(0,va)*g

                    2       2  2   2
      + 3*cos(theta) *z(0,v) *g *v0

                          2                           2
      - 6*cos(theta)*h0(0) *p(0,v)*sin(theta)*z(0,v)*g

                                                      2
      - 12*cos(theta)*h0(0)*p(0,v)*sin(theta)*z(0,v)*g *v0

                                               2   2
      - 6*cos(theta)*p(0,v)*sin(theta)*z(0,v)*g *v0  - 24*ef(-1)*gg(v)

     *ef(1)*cos(theta)*p(0,v)*sin(theta)*te1(1,ef)*te1(-1,ef)*g - 12

                                        2
     *ef(-1)*gg(v)*l(1)*ef(1)*cos(theta) *se1(1,ef)*se1(-1,ef)*z(0,v)

     *g - 12*ef(-1)*gg(v)*l(1)*ef(1)*cos(theta)*p(0,v)*se1(1,ef)

     *se1(-1,ef)*sin(theta)*g + 24*ef(-1)*gg(v)*l(1)*ef(1)*cos(theta)

     *p(0,v)*sin(theta)*te1(1,ef)*te1(-1,ef)*g + 12*sqrt(2)*ef(-1)

                                  2
     *gg(v)*l(1)*nue(1)*cos(theta) *se1(1,nue)*se1(-1,ef)*w(-1,v)*g + 

                                                2
     24*ef(-1)*gg(v)*l(1)*pd(ef(1),v)*cos(theta) *se1(1,ef)*se1(-1,ef)

                                                     2
     *i - 24*ef(-1)*gg(v)*l(1)*pd(ef(1),v)*cos(theta) *te1(1,ef)

     *te1(-1,ef)*i

                                              2
      + 24*ef(-1)*gg(v)*pd(ef(1),v)*cos(theta) *te1(1,ef)*te1(-1,ef)*i

               2       2           2  2
      + 3*h0(0) *p(0,v) *sin(theta) *g

                      2           2  2
      + 6*h0(0)*p(0,v) *sin(theta) *g *v0 + 12*sqrt(2)*nue(-1)*gg(v)

                           2
     *l(1)*ef(1)*cos(theta) *se1(1,ef)*se1(-1,nue)*w(1,v)*g + 12

                                          2
     *nue(-1)*gg(v)*l(1)*nue(1)*cos(theta) *se1(1,nue)*se1(-1,nue)

     *z(0,v)*g - 12*nue(-1)*gg(v)*l(1)*nue(1)*cos(theta)*p(0,v)

     *se1(1,nue)*se1(-1,nue)*sin(theta)*g + 24*nue(-1)*gg(v)*l(1)

                             2
     *pd(nue(1),v)*cos(theta) *se1(1,nue)*se1(-1,nue)*i

                2           2  2   2
      + 3*p(0,v) *sin(theta) *g *v0

                             2           2
      - 12*pd(gs(0,v,ic1),va) *cos(theta)

                                                           2
      + 12*pd(gs(0,v,ic1),va)*pd(gs(0,va,ic1),v)*cos(theta)

                      2           2                   2           2
      + 12*pd(h0(0),v) *cos(theta)  - 12*pd(p(0,v),va) *cos(theta)

                                                 2
      + 12*pd(p(0,v),va)*pd(p(0,va),v)*cos(theta)

                                   2
      + 24*pd(w(1,v),va)*cos(theta) *w(-1,v)*z(0,va)*g*i

                                   2
      - 24*pd(w(1,v),va)*cos(theta) *w(-1,va)*z(0,v)*g*i

                                                  2
      - 12*pd(w(1,v),va)*pd(w(-1,v),va)*cos(theta)

                                                  2
      + 12*pd(w(1,v),va)*pd(w(-1,va),v)*cos(theta)

                                    2
      - 24*pd(w(-1,v),va)*cos(theta) *w(1,v)*z(0,va)*g*i

                                    2
      + 24*pd(w(-1,v),va)*cos(theta) *w(1,va)*z(0,v)*g*i

                                                  2
      - 12*pd(w(-1,v),va)*pd(w(1,v),va)*cos(theta)

                                                  2
      + 12*pd(w(-1,v),va)*pd(w(1,va),v)*cos(theta)

                        2           2
      - 12*pd(z(0,v),va) *cos(theta)

                                   2
      + 24*pd(z(0,v),va)*cos(theta) *w(1,v)*w(-1,va)*g*i

                                   2
      - 24*pd(z(0,v),va)*cos(theta) *w(1,va)*w(-1,v)*g*i

                                                 2
      + 12*pd(z(0,v),va)*pd(z(0,va),v)*cos(theta)  + 24*qd(-1,id1)

                                           2
     *gg(v)*l(1)*pd(qd(1,id1),v)*cos(theta) *s1(1,qd)*s1(-1,qd)*i - 24

                                                      2
     *qd(-1,id1)*gg(v)*l(1)*pd(qd(1,id1),v)*cos(theta) *t1(1,qd)

                                                                 2
     *t1(-1,qd)*i - 12*qd(-1,id1)*gg(v)*l(1)*qd(1,id1)*cos(theta)

     *s1(1,qd)*s1(-1,qd)*z(0,v)*g + 4*qd(-1,id1)*gg(v)*l(1)*qd(1,id1)

     *cos(theta)*p(0,v)*s1(1,qd)*s1(-1,qd)*sin(theta)*g + 8*qd(-1,id1)

     *gg(v)*l(1)*qd(1,id1)*cos(theta)*p(0,v)*sin(theta)*t1(1,qd)

     *t1(-1,qd)*g + 12*sqrt(2)*qd(-1,id1)*gg(v)*l(1)*qu(1,id1)

                2
     *cos(theta) *s1(1,qu)*s1(-1,qd)*w(-1,v)*g + 24*qd(-1,id1)*gg(v)

                                2
     *pd(qd(1,id1),v)*cos(theta) *t1(1,qd)*t1(-1,qd)*i - 8*qd(-1,id1)

     *gg(v)*qd(1,id1)*cos(theta)*p(0,v)*sin(theta)*t1(1,qd)*t1(-1,qd)

                                                             2
     *g + 24*qu(-1,id1)*gg(v)*l(1)*pd(qu(1,id1),v)*cos(theta)

     *s1(1,qu)*s1(-1,qu)*i - 24*qu(-1,id1)*gg(v)*l(1)*pd(qu(1,id1),v)

                2
     *cos(theta) *t2(1,qu)*t2(-1,qu)*i + 12*sqrt(2)*qu(-1,id1)*gg(v)

                               2
     *l(1)*qd(1,id1)*cos(theta) *s1(1,qd)*s1(-1,qu)*w(1,v)*g + 12

                                                2
     *qu(-1,id1)*gg(v)*l(1)*qu(1,id1)*cos(theta) *s1(1,qu)*s1(-1,qu)

     *z(0,v)*g + 4*qu(-1,id1)*gg(v)*l(1)*qu(1,id1)*cos(theta)*p(0,v)

     *s1(1,qu)*s1(-1,qu)*sin(theta)*g - 16*qu(-1,id1)*gg(v)*l(1)

     *qu(1,id1)*cos(theta)*p(0,v)*sin(theta)*t2(1,qu)*t2(-1,qu)*g + 24

                                                 2
     *qu(-1,id1)*gg(v)*pd(qu(1,id1),v)*cos(theta) *t2(1,qu)*t2(-1,qu)

     *i + 16*qu(-1,id1)*gg(v)*qu(1,id1)*cos(theta)*p(0,v)*sin(theta)

                                      2
     *t2(1,qu)*t2(-1,qu)*g)/cos(theta)

 

hh== - 2*gsg(-1,)*pd(pd(gsg(1,ic1),v),v)

     - pd(gs(0,v,ic1),v)*pd(gs(0,va,ic1),va)

 

hh==c(1,4)*se1(-1,ef)*te1(1,ef)*(ef(-1)*ef(1)*h0(0) + ef(-1)*ef(1)*v0

        - ef(-1)*l(1)*ef(1)*h0(0) - ef(-1)*l(1)*ef(1)*v0)

 

hh==c(-1,4)*se1(1,ef)*te1(-1,ef)

    *(ef(-1)*l(1)*ef(1)*h0(0) + ef(-1)*l(1)*ef(1)*v0)

 

                 2                  2
hh==c(0,1)*(h0(0)  + 2*h0(0)*v0 + v0 )

 

                 4          3             2   2             3     4
hh==c(0,5)*(h0(0)  + 4*h0(0) *v0 + 6*h0(0) *v0  + 4*h0(0)*v0  + v0 )

 

hh==c(-1,3)*s1(1,qd)*t1(-1,qd)

    *(qd(-1,id1)*l(1)*qd(1,id1)*h0(0) + qd(-1,id1)*l(1)*qd(1,id1)*v0)

 

hh==c(-1,2)*s1(1,qu)*t2(-1,qu)

    *(qu(-1,id1)*l(1)*qu(1,id1)*h0(0) + qu(-1,id1)*l(1)*qu(1,id1)*v0)

 

hh==c(1,3)*s1(-1,qd)*t1(1,qd)*( - qd(-1,id1)*l(1)*qd(1,id1)*h0(0)

        - qd(-1,id1)*l(1)*qd(1,id1)*v0 + qd(-1,id1)*qd(1,id1)*h0(0)

        + qd(-1,id1)*qd(1,id1)*v0)

 

hh==c(1,2)*s1(-1,qu)*t2(1,qu)*( - qu(-1,id1)*l(1)*qu(1,id1)*h0(0)

        - qu(-1,id1)*l(1)*qu(1,id1)*v0 + qu(-1,id1)*qu(1,id1)*h0(0)

        + qu(-1,id1)*qu(1,id1)*v0)

 

                        2           2
                      v0 *(c(0,5)*v0  + 2*c(0,1))
The vancum energy is=-----------------------------
                                   4

 

The particle is =h0(0,r1)

 

                              2
The equation is =v0*(c(0,5)*v0  + c(0,1))

 
This is mix in particles =(p z)

mkequation equationl bb=0

 

mkequation equationl bb=0

 

mkequation equationl bb=

                      2   2               2             2
 gm(v,va)*sin(theta)*g *v0 *( - cos(theta)  - sin(theta)  + 1)
---------------------------------------------------------------
                                     2
                         8*cos(theta)

 

mkequation equationl bb=

           2   2            2             2
 gm(v,va)*g *v0 *(cos(theta)  + sin(theta)  - 1)
-------------------------------------------------
                  8*cos(theta)

 
The particle list is =(p z)
something wrong for a's progaptor a==nue
something wrong for a's progaptor a==numu
something wrong for a's progaptor a==nut
something wrong for a's progaptor a==qd
something wrong for a's progaptor a==qs
something wrong for a's progaptor a==qb
something wrong for a's progaptor a==ef
something wrong for a's progaptor a==mu
something wrong for a's progaptor a==tau
something wrong for a's progaptor a==qu
something wrong for a's progaptor a==qc
something wrong for a's progaptor a==qt

     - sqrt(2)*fmt*g + 2*wangm9*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fmt*g + 2*wangm9*wm

 

     - sqrt(2)*fmc*g + 2*wangm8*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fmc*g + 2*wangm8*wm

 

     - sqrt(2)*fmu*g + 2*wangm7*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fmu*g + 2*wangm7*wm

 

     - sqrt(2)*fmtau*g + 2*wangm6*wm
0==----------------------------------
                   g

 

0== - sqrt(2)*fmtau*g + 2*wangm6*wm

 

     - sqrt(2)*fmmu*g + 2*wangm5*wm
0==---------------------------------
                   g

 

0== - sqrt(2)*fmmu*g + 2*wangm5*wm

 

     - sqrt(2)*fme*g + 2*wangm4*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fme*g + 2*wangm4*wm

 

     - sqrt(2)*fmb*g + 2*wangm3*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fmb*g + 2*wangm3*wm

 

     - sqrt(2)*fms*g + 2*wangm2*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fms*g + 2*wangm2*wm

 

     - sqrt(2)*fmd*g + 2*wangm1*wm
0==--------------------------------
                  g

 

0== - sqrt(2)*fmd*g + 2*wangm1*wm

 

                2           2    2   2
    12*c(0,5)*wm  + c(0,1)*g  + g *hm
0==------------------------------------
                     2
                    g

 

               2           2    2   2
0==12*c(0,5)*wm  + c(0,1)*g  + g *hm

 

                     2           2
    2*wm*(4*c(0,5)*wm  + c(0,1)*g )
0==---------------------------------
                   3
                  g

 

                    2           2
0==2*wm*(4*c(0,5)*wm  + c(0,1)*g )

 
couldn't found the higgs be eated by a=w
couldn't found the higgs be eated by a=z
couldn't found the higgs be eated by a=p
htionlist=((g (2 47 8 7 6 5 2 1) (1 54 53 52 51 50 49 48 46
45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26
25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 4 3)))

Time: 0 ms
 You are building the 
The Phenomenology Lagrangian Model for $J/\Psi$ Physics
 in the directory ../model/.
 The Feynman rules are generated and placed in the file 
../model/feynman_rule_show.
 Please go to look at it and find out which vertices you want to drop.
 Then you should created a file ../model/vertices_drop_list
 in which you write a list like:
 $
   droplist:='(2 3 7 ...... 19); 
 $
   ;end; 
 $
 This means that vertices 2,3,7,......,19 will be droped
 You can 'bye' REDUCE and do it and
 restat the procedure without performing model_generator1() and model_generator2(),
 and directly perform model_generator3(). You can also continue to perform 
 model_generator3() without 'bye' REDUCE if there is another way for you to
 revise the file
     *  ...... Successed in second step
7: 
Quitting
