Index: trunk/omega/tests/keystones_omegalib_bispinors_generate.ml =================================================================== --- trunk/omega/tests/keystones_omegalib_bispinors_generate.ml (revision 8488) +++ trunk/omega/tests/keystones_omegalib_bispinors_generate.ml (revision 8489) @@ -1,53 +1,53 @@ (* keystones_omegalib_bispinors_generate.ml -- Copyright (C) 2019-2020 by Wolfgang Kilian Thorsten Ohl Juergen Reuter WHIZARD is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2, or (at your option) any later version. WHIZARD is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. *) open Coupling open Keystones let vector_spinor_current tag = { tag = Printf.sprintf "vector_spinor_current__%s_ff" tag; keystones = [ { bra = (Majorana, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Vector, 1); F (Majorana, 2)] }; { bra = (Vector, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (Majorana, 0); F (Majorana, 2)] } ] } let scalar_spinor_current tag = { tag = Printf.sprintf "scalar_spinor_current__%s_ff" tag; keystones = [ { bra = (Majorana, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Scalar, 1); F (Majorana, 2)] }; { bra = (Scalar, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (Majorana, 0); F (Majorana, 2)] } ] } let vertices = List.concat [ List.map vector_spinor_current ["v"; "a"; "vl"; "vr"]; List.map scalar_spinor_current ["s"; "p"; "sl"; "sr"] ] let _ = Keystones.generate - ~reps:10000 ~threshold:0.70 ~omega_module:"omega95_bispinors" + ~reps:1000 ~threshold:0.70 ~omega_module:"omega95_bispinors" vertices; exit 0 Index: trunk/omega/tests/comparisons_recola.list =================================================================== --- trunk/omega/tests/comparisons_recola.list (revision 8488) +++ trunk/omega/tests/comparisons_recola.list (revision 8489) @@ -1,38 +1,38 @@ ### comparisons_recola.list -- #! ### -------------------------------------------------------------------- -### thr abs n roots model mode process ... +### thr abs n roots model mode process ... ### -------------------------------------------------------------------- ### QED -eeee 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> e+ e- -eemm 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> mu+ mu- -eeaa 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> A A -eeaaa 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> A A A -eeaaaa 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> A A A A +eeee 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> e+ e- +eemm 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> mu+ mu- +eeaa 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> A A +eeaaa 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> A A A +eeaaaa 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> A A A A ### -------------------------------------------------------------------- ### NC -eezz 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> Z Z +eezz 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> Z Z ### -------------------------------------------------------------------- ### CC -eeww 0.25 1E-30 1000 1000 SM_Higgs scatter e+ e- -> W+ W- +eeww 0.25 1E-30 100 1000 SM_Higgs scatter e+ e- -> W+ W- ### VBS -#wwww 0.25 1E-30 1000 1000 SM_Higgs scatter W+ W- -> W+ W- -#zzww 0.25 1E-30 1000 1000 SM_Higgs scatter Z Z -> W+ W- +#wwww 0.25 1E-30 100 1000 SM_Higgs scatter W+ W- -> W+ W- +#zzww 0.25 1E-30 100 1000 SM_Higgs scatter Z Z -> W+ W- ### -------------------------------------------------------------------- ### QCD -#gggg 0.25 1E-30 1000 1000 SM_Higgs scatter g g -> g g +#gggg 0.25 1E-30 100 1000 SM_Higgs scatter g g -> g g ### -------------------------------------------------------------------- ### Misc -wwhh 0.10 1E-30 100 10000 SM_Higgs scatter mu+ e- -> numubar nue H H | mu+ e- -> nu_mu~ nu_e H H -wwhhh 0.05 1E-30 100 10000 SM_Higgs scatter mu+ e- -> numubar nue H H H | mu+ e- -> nu_mu~ nu_e H H H +wwhh 0.10 1E-30 50 10000 SM_Higgs scatter mu+ e- -> numubar nue H H | mu+ e- -> nu_mu~ nu_e H H +wwhhh 0.05 1E-30 50 10000 SM_Higgs scatter mu+ e- -> numubar nue H H H | mu+ e- -> nu_mu~ nu_e H H H ### -------------------------------------------------------------------- ### recola particles ### scalars: 'H', 'p0', 'p+', 'p-' ### vector bosons: 'g', 'A', 'Z', 'W+', 'W-' ### leptons: 'nu_e', 'nu_e~', 'e-', 'e+', ### 'nu_mu', 'nu_mu~', 'mu-', 'mu+', ### 'nu_tau', 'nu_tau~', 'tau-', 'tau+' ### quarks: 'u', 'u~', 'd', 'd~', ### 'c', 'c~', 's', 's~', ### 't', 't~', 'b', 'b~' ### ---------------------------------------------------------------------- Index: trunk/omega/tests/keystones_omegalib_generate.ml =================================================================== --- trunk/omega/tests/keystones_omegalib_generate.ml (revision 8488) +++ trunk/omega/tests/keystones_omegalib_generate.ml (revision 8489) @@ -1,79 +1,79 @@ (* keystones_omegalib_generate.ml -- Copyright (C) 2019-2020 by Wolfgang Kilian Thorsten Ohl Juergen Reuter WHIZARD is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2, or (at your option) any later version. WHIZARD is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. *) open Coupling open Keystones let vector_spinor_current tag = { tag = Printf.sprintf "vector_spinor_current__%s_ff" tag; keystones = [ { bra = (ConjSpinor, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Vector, 1); F (Spinor, 2)] }; { bra = (Vector, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (ConjSpinor, 0); F (Spinor, 2)] }; { bra = (Spinor, 2); name = Printf.sprintf "f_f%s" tag; args = [G (0); F (ConjSpinor, 0); F (Vector, 1)] } ] } let scalar_spinor_current tag = { tag = Printf.sprintf "scalar_spinor_current__%s_ff" tag; keystones = [ { bra = (ConjSpinor, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Scalar, 1); F (Spinor, 2)] }; { bra = (Scalar, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (ConjSpinor, 0); F (Spinor, 2)] }; { bra = (Spinor, 2); name = Printf.sprintf "f_f%s" tag; args = [G (0); F (ConjSpinor, 0); F (Scalar, 1)] } ] } (* NB: the vertex is anti-symmetric in the scalars and we need to use a cyclic permutation. *) let vector_scalar_current = { tag = "vector_scalar_current__v_ss"; keystones = [ { bra = (Vector, 0); name = "v_ss"; args = [G (0); F (Scalar, 1); P (1); F (Scalar, 2); P (2)] }; { bra = (Scalar, 2); name = "s_vs"; args = [G (0); F (Vector, 0); P (0); F (Scalar, 1); P (1)] } ] } let scalar_vector_current tag = { tag = Printf.sprintf "transversal_vector_current__s_vv_%s" tag; keystones = [ { bra = (Scalar, 0); name = Printf.sprintf "s_vv_%s" tag; args = [G (0); F (Vector, 1); P (1); F (Vector, 2); P (2)] }; { bra = (Vector, 1); name = Printf.sprintf "v_sv_%s" tag; args = [G (0); F (Scalar, 0); P (0); F (Vector, 2); P (2)] } ] } let vertices = List.concat [ List.map vector_spinor_current ["v"; "a"; "vl"; "vr"]; List.map scalar_spinor_current ["s"; "p"; "sl"; "sr"]; [ vector_scalar_current ]; List.map scalar_vector_current ["t"; "6D"; "6DP"] ] let _ = - Keystones.generate ~reps:10000 ~threshold:0.70 vertices; + Keystones.generate ~reps:1000 ~threshold:0.70 vertices; exit 0 Index: trunk/omega/tests/ward_identities.list =================================================================== --- trunk/omega/tests/ward_identities.list (revision 8488) +++ trunk/omega/tests/ward_identities.list (revision 8489) @@ -1,92 +1,92 @@ # ward_identities.list -- # ---------------------------------------------------------------------- # thr n roots model i process ... # ---------------------------------------------------------------------- -eeaa 0.75 1000 1000 QED 3 scatter e+ e- -> A A -eeaaa 0.75 1000 1000 QED 3 scatter e+ e- -> A A A -eeaee 0.35 1000 1000 QED 3 scatter e+ e- -> A e+ e- -aeae 0.35 1000 1000 QED 3 scatter A e+ -> A e+ -aeaea 0.35 1000 1000 QED 3 scatter A e+ -> A e+ A -# -uugg 0.70 1000 1000 QCD 3 scatter u ubar -> gl gl -gggg 0.70 1000 1000 QCD 3 scatter gl gl -> gl gl -uuguu 0.35 500 1000 QCD 3 scatter u ubar -> gl u ubar -gugu 0.30 1000 1000 QCD 3 scatter gl u -> gl u +eeaa 0.75 100 1000 QED 3 scatter e+ e- -> A A +eeaaa 0.75 100 1000 QED 3 scatter e+ e- -> A A A +eeaee 0.35 100 1000 QED 3 scatter e+ e- -> A e+ e- +aeae 0.35 100 1000 QED 3 scatter A e+ -> A e+ +aeaea 0.35 100 1000 QED 3 scatter A e+ -> A e+ A +# +uugg 0.70 100 1000 QCD 3 scatter u ubar -> gl gl +gggg 0.70 100 1000 QCD 3 scatter gl gl -> gl gl +uuguu 0.35 50 1000 QCD 3 scatter u ubar -> gl u ubar +gugu 0.30 100 1000 QCD 3 scatter gl u -> gl u # numerical noise improves w/quad: -uugguu 0.25 200 1000 QCD 3 scatter u ubar -> gl gl u ubar +uugguu 0.25 50 1000 QCD 3 scatter u ubar -> gl gl u ubar # -qqgg1 0.70 1000 1000 SYM 3 scatter q1 Q1 -> g1 g1 -ssgg1 0.70 1000 1000 SYM 3 scatter sq1 sQ1 -> g1 g1 -ppgg1 0.70 1000 1000 SYM 3 scatter phi phi -> g1 g1 -gggg1 0.70 1000 1000 SYM 3 scatter g1 g1 -> g1 g1 -# -ggtt 0.35 1000 1000 SM 1 scatter gl gl -> t tbar -gtgt 0.35 1000 1000 SM 3 scatter gl t -> gl t -ggtt_a 0.35 1000 1000 SM_top_anom 1 scatter gl gl -> t tbar -gtgt_a 0.35 1000 1000 SM_top_anom 3 scatter gl t -> gl t -# -aatt 0.30 1000 1000 SM 3 scatter t tbar -> A A -ttaa 0.30 1000 1000 SM 1 scatter A A -> t tbar -tata 0.30 1000 1000 SM 3 scatter A t -> A t -atat 0.30 1000 1000 SM 3 scatter A tbar -> A tbar -aatt_a 0.30 1000 1000 SM_top_anom 3 scatter t tbar -> A A -ttaa_a 0.30 1000 1000 SM_top_anom 1 scatter A A -> t tbar -tata_a 0.30 1000 1000 SM_top_anom 3 scatter A t -> A t -atat_a 0.30 1000 1000 SM_top_anom 3 scatter A tbar -> A tbar -# -#ttaz 0.30 1000 1000 SM 3 scatter t tbar -> A Z -#aztt 0.30 1000 1000 SM 1 scatter A Z -> t tbar -#tzta 0.30 1000 1000 SM 3 scatter Z t -> A t -#ztat 0.30 1000 1000 SM 3 scatter Z tbar -> A tbar -#ttaz_a 0.30 1000 1000 SM_top_anom 3 scatter t tbar -> A Z -#aztt_a 0.30 1000 1000 SM_top_anom 1 scatter A Z -> t tbar -#tzta_a 0.30 1000 1000 SM_top_anom 3 scatter Z t -> A t -#ztat_a 0.30 1000 1000 SM_top_anom 3 scatter Z tbar -> A tbar -# -#ttza 0.30 1000 1000 SM 3 scatter t tbar -> Z A -#zatt 0.30 1000 1000 SM 1 scatter Z A -> t tbar -#tatz 0.30 1000 1000 SM 3 scatter A t -> Z t -#atzt 0.30 1000 1000 SM 3 scatter A tbar -> Z tbar -#ttza_a 0.30 1000 1000 SM_top_anom 3 scatter t tbar -> Z A -#zatt_a 0.30 1000 1000 SM_top_anom 1 scatter Z A -> t tbar -#tatz_a 0.30 1000 1000 SM_top_anom 3 scatter A t -> Z t -#atzt_a 0.30 1000 1000 SM_top_anom 3 scatter A tbar -> Z tbar -# -#btwz 0.35 1000 1000 SM 4 scatter t bbar -> W+ Z -#wbzt 0.35 1000 1000 SM 3 scatter W+ b -> Z t -#btwz_a 0.35 1000 1000 SM_top_anom 4 scatter t bbar -> W+ Z -#wbzt_a 0.35 1000 1000 SM_top_anom 3 scatter W+ b -> Z t -# -#tbwz 0.35 1000 1000 SM 4 scatter b tbar -> W- Z -#wtzb 0.35 1000 1000 SM 3 scatter W- t -> Z b -#tbwz_a 0.35 1000 1000 SM_top_anom 4 scatter b tbar -> W- Z -#wtzb_a 0.35 1000 1000 SM_top_anom 3 scatter W- t -> Z b -# -#btwa 0.35 1000 1000 SM 4 scatter t bbar -> W+ A -#wbat 0.35 1000 1000 SM 3 scatter W+ b -> A t -#btwa_a 0.35 1000 1000 SM_top_anom 4 scatter t bbar -> W+ A -#wbat_a 0.35 1000 1000 SM_top_anom 3 scatter W+ b -> A t -# -#tbaw 0.35 1000 1000 SM 3 scatter b tbar -> A W- -#wtab 0.35 1000 1000 SM 3 scatter W- t -> A b -#tbaw_a 0.35 1000 1000 SM_top_anom 3 scatter b tbar -> A W- -#wtab_a 0.35 1000 1000 SM_top_anom 3 scatter W- t -> A b -# -#ttzz 0.30 1000 1000 SM 3 scatter t tbar -> Z Z -#zztt 0.30 1000 1000 SM 1 scatter Z Z -> t tbar -#tztz 0.30 1000 1000 SM 3 scatter Z t -> Z t -#ztzt 0.30 1000 1000 SM 3 scatter Z tbar -> Z tbar -#ttzz_a 0.30 1000 1000 SM_top_anom 3 scatter t tbar -> Z Z -#zztt_a 0.30 1000 1000 SM_top_anom 1 scatter Z Z -> t tbar -#tztz_a 0.30 1000 1000 SM_top_anom 3 scatter Z t -> Z t -#ztzt_a 0.30 1000 1000 SM_top_anom 3 scatter Z tbar -> Z tbar -# -#wwtt 0.35 1000 1000 SM 3 scatter t tbar -> W+ W- -#wtwt 0.35 1000 1000 SM 3 scatter W+ t -> W+ t -#wwtt_a 0.35 1000 1000 SM_top_anom 3 scatter t tbar -> W+ W- -#wtwt_a 0.35 1000 1000 SM_top_anom 3 scatter W+ t -> W+ t -# -#wwbb 0.35 1000 1000 SM 3 scatter b bbar -> W+ W- -#wbwb 0.35 1000 1000 SM 3 scatter W+ b -> W+ b -#wwbb_a 0.35 1000 1000 SM_top_anom 3 scatter b bbar -> W+ W- -#wbwb_a 0.35 1000 1000 SM_top_anom 3 scatter W+ b -> W+ b +qqgg1 0.70 100 1000 SYM 3 scatter q1 Q1 -> g1 g1 +ssgg1 0.70 100 1000 SYM 3 scatter sq1 sQ1 -> g1 g1 +ppgg1 0.70 100 1000 SYM 3 scatter phi phi -> g1 g1 +gggg1 0.70 100 1000 SYM 3 scatter g1 g1 -> g1 g1 +# +ggtt 0.35 100 1000 SM 1 scatter gl gl -> t tbar +gtgt 0.35 100 1000 SM 3 scatter gl t -> gl t +ggtt_a 0.35 100 1000 SM_top_anom 1 scatter gl gl -> t tbar +gtgt_a 0.35 100 1000 SM_top_anom 3 scatter gl t -> gl t +# +aatt 0.30 100 1000 SM 3 scatter t tbar -> A A +ttaa 0.30 100 1000 SM 1 scatter A A -> t tbar +tata 0.30 100 1000 SM 3 scatter A t -> A t +atat 0.30 100 1000 SM 3 scatter A tbar -> A tbar +aatt_a 0.30 100 1000 SM_top_anom 3 scatter t tbar -> A A +ttaa_a 0.30 100 1000 SM_top_anom 1 scatter A A -> t tbar +tata_a 0.30 100 1000 SM_top_anom 3 scatter A t -> A t +atat_a 0.30 100 1000 SM_top_anom 3 scatter A tbar -> A tbar +# +#ttaz 0.30 100 1000 SM 3 scatter t tbar -> A Z +#aztt 0.30 100 1000 SM 1 scatter A Z -> t tbar +#tzta 0.30 100 1000 SM 3 scatter Z t -> A t +#ztat 0.30 100 1000 SM 3 scatter Z tbar -> A tbar +#ttaz_a 0.30 100 1000 SM_top_anom 3 scatter t tbar -> A Z +#aztt_a 0.30 100 1000 SM_top_anom 1 scatter A Z -> t tbar +#tzta_a 0.30 100 1000 SM_top_anom 3 scatter Z t -> A t +#ztat_a 0.30 100 1000 SM_top_anom 3 scatter Z tbar -> A tbar +# +#ttza 0.30 100 1000 SM 3 scatter t tbar -> Z A +#zatt 0.30 100 1000 SM 1 scatter Z A -> t tbar +#tatz 0.30 100 1000 SM 3 scatter A t -> Z t +#atzt 0.30 100 1000 SM 3 scatter A tbar -> Z tbar +#ttza_a 0.30 100 1000 SM_top_anom 3 scatter t tbar -> Z A +#zatt_a 0.30 100 1000 SM_top_anom 1 scatter Z A -> t tbar +#tatz_a 0.30 100 1000 SM_top_anom 3 scatter A t -> Z t +#atzt_a 0.30 100 1000 SM_top_anom 3 scatter A tbar -> Z tbar +# +#btwz 0.35 100 1000 SM 4 scatter t bbar -> W+ Z +#wbzt 0.35 100 1000 SM 3 scatter W+ b -> Z t +#btwz_a 0.35 100 1000 SM_top_anom 4 scatter t bbar -> W+ Z +#wbzt_a 0.35 100 1000 SM_top_anom 3 scatter W+ b -> Z t +# +#tbwz 0.35 100 1000 SM 4 scatter b tbar -> W- Z +#wtzb 0.35 100 1000 SM 3 scatter W- t -> Z b +#tbwz_a 0.35 100 1000 SM_top_anom 4 scatter b tbar -> W- Z +#wtzb_a 0.35 100 1000 SM_top_anom 3 scatter W- t -> Z b +# +#btwa 0.35 100 1000 SM 4 scatter t bbar -> W+ A +#wbat 0.35 100 1000 SM 3 scatter W+ b -> A t +#btwa_a 0.35 100 1000 SM_top_anom 4 scatter t bbar -> W+ A +#wbat_a 0.35 100 1000 SM_top_anom 3 scatter W+ b -> A t +# +#tbaw 0.35 100 1000 SM 3 scatter b tbar -> A W- +#wtab 0.35 100 1000 SM 3 scatter W- t -> A b +#tbaw_a 0.35 100 1000 SM_top_anom 3 scatter b tbar -> A W- +#wtab_a 0.35 100 1000 SM_top_anom 3 scatter W- t -> A b +# +#ttzz 0.30 100 1000 SM 3 scatter t tbar -> Z Z +#zztt 0.30 100 1000 SM 1 scatter Z Z -> t tbar +#tztz 0.30 100 1000 SM 3 scatter Z t -> Z t +#ztzt 0.30 100 1000 SM 3 scatter Z tbar -> Z tbar +#ttzz_a 0.30 100 1000 SM_top_anom 3 scatter t tbar -> Z Z +#zztt_a 0.30 100 1000 SM_top_anom 1 scatter Z Z -> t tbar +#tztz_a 0.30 100 1000 SM_top_anom 3 scatter Z t -> Z t +#ztzt_a 0.30 100 1000 SM_top_anom 3 scatter Z tbar -> Z tbar +# +#wwtt 0.35 100 1000 SM 3 scatter t tbar -> W+ W- +#wtwt 0.35 100 1000 SM 3 scatter W+ t -> W+ t +#wwtt_a 0.35 100 1000 SM_top_anom 3 scatter t tbar -> W+ W- +#wtwt_a 0.35 100 1000 SM_top_anom 3 scatter W+ t -> W+ t +# +#wwbb 0.35 100 1000 SM 3 scatter b bbar -> W+ W- +#wbwb 0.35 100 1000 SM 3 scatter W+ b -> W+ b +#wwbb_a 0.35 100 1000 SM_top_anom 3 scatter b bbar -> W+ W- +#wbwb_a 0.35 100 1000 SM_top_anom 3 scatter W+ b -> W+ b Index: trunk/omega/tests/benchmark_processes.list =================================================================== --- trunk/omega/tests/benchmark_processes.list (revision 8488) +++ trunk/omega/tests/benchmark_processes.list (revision 8489) @@ -1,38 +1,38 @@ # processes for benchmarking that can be used in different scripts # # Use `n` points for a fixed `roots` for the scattering `process` in their # respective `model` -# ---------------------------------------------------------------------- -# n roots model process -# ---------------------------------------------------------------------- -#4g 10000 1000 QCD scatter gl gl -> gl gl -#5g 5000 1000 QCD scatter gl gl -> gl gl gl -6g 2310 1000 QCD scatter gl gl -> gl gl gl gl -#7g 2310 1000 QCD scatter gl gl -> gl gl gl gl gl +# --------------------------------------------------------------------- +# n roots model process +# --------------------------------------------------------------------- +#4g 1000 1000 QCD scatter gl gl -> gl gl +#5g 500 1000 QCD scatter gl gl -> gl gl gl +6g 231 1000 QCD scatter gl gl -> gl gl gl gl +#7g 231 1000 QCD scatter gl gl -> gl gl gl gl gl -#DY0j 10000 1000 SM scatter u ubar -> e+ e- -#DY1j 10000 1000 SM scatter u ubar -> e+ e- u:ubar:gl -#DY2j 10000 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl -#DY3j 10000 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl -#DY4j 2310 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl -#DY5j 2310 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl -#DY6j 210 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl +#DY0j 1000 1000 SM scatter u ubar -> e+ e- +#DY1j 1000 1000 SM scatter u ubar -> e+ e- u:ubar:gl +#DY2j 1000 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl +#DY3j 1000 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl +#DY4j 231 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl +#DY5j 231 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl +#DY6j 21 1000 SM scatter u ubar -> e+ e- u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl -#4jets 10000 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl -#5jets 10000 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl -#6jets 2310 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl -#7jets 210 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl +#4jets 1000 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl +#5jets 1000 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl +#6jets 231 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl +#7jets 21 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl -#ee2a 10000 1000 QED scatter e+ e- -> A A -#ee3a 10000 1000 QED scatter e+ e- -> A A A -#ee4a 10000 1000 QED scatter e+ e- -> A A A A -#ee5a 10000 1000 QED scatter e+ e- -> A A A A A -#ee6a 10000 1000 QED scatter e+ e- -> A A A A A A -#ee7a 5000 1000 QED scatter e+ e- -> A A A A A A A -#ee8a 2310 1000 QED scatter e+ e- -> A A A A A A A A -#ee9a 2310 1000 QED scatter e+ e- -> A A A A A A A A A -#ee10a 2310 1000 QED scatter e+ e- -> A A A A A A A A A A +#ee2a 1000 1000 QED scatter e+ e- -> A A +#ee3a 1000 1000 QED scatter e+ e- -> A A A +#ee4a 1000 1000 QED scatter e+ e- -> A A A A +#ee5a 1000 1000 QED scatter e+ e- -> A A A A A +#ee6a 1000 1000 QED scatter e+ e- -> A A A A A A +#ee7a 500 1000 QED scatter e+ e- -> A A A A A A A +#ee8a 231 1000 QED scatter e+ e- -> A A A A A A A A +#ee9a 231 1000 QED scatter e+ e- -> A A A A A A A A A +#ee10a 231 1000 QED scatter e+ e- -> A A A A A A A A A A # Doesn't compile/link with gfortran or ifort, only possible with OVM -#8g 210 1000 QCD scatter gl gl -> gl gl gl gl gl gl -#8jets 200 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl +#8g 21 1000 QCD scatter gl gl -> gl gl gl gl gl gl +#8jets 20 1000 QCD scatter u:ubar:gl u:ubar:gl -> u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl u:ubar:gl Index: trunk/omega/tests/ward_identities_UFO.list =================================================================== --- trunk/omega/tests/ward_identities_UFO.list (revision 8488) +++ trunk/omega/tests/ward_identities_UFO.list (revision 8489) @@ -1,16 +1,16 @@ # ward_identities_UFO.list -- #! # ---------------------------------------------------------------------- -# thr n roots model i process ... +# thr n roots model i process ... # ---------------------------------------------------------------------- -eeaa 0.30 1000 1000 SM 3 scatter e+ e- -> a a -eeaaa 0.30 1000 1000 SM 3 scatter e+ e- -> a a a -uuaa 0.25 1000 1000 SM 3 scatter u u~ -> a a -uuga 0.25 1000 1000 SM 3 scatter u u~ -> g a -uuag 0.25 1000 1000 SM 3 scatter u u~ -> a g -uugg 0.25 1000 1000 SM 3 scatter u u~ -> g g -uugga 0.25 1000 1000 SM 3 scatter u u~ -> g g a -uuagg 0.20 1000 1000 SM 3 scatter u u~ -> a g g -uuggg 0.20 1000 1000 SM 3 scatter u u~ -> g g g -gggg 0.50 1000 1000 SM 3 scatter g g -> g g -ggggg 0.50 1000 1000 SM 3 scatter g g -> g g g +eeaa 0.30 100 1000 SM 3 scatter e+ e- -> a a +eeaaa 0.30 100 1000 SM 3 scatter e+ e- -> a a a +uuaa 0.25 100 1000 SM 3 scatter u u~ -> a a +uuga 0.25 100 1000 SM 3 scatter u u~ -> g a +uuag 0.25 100 1000 SM 3 scatter u u~ -> a g +uugg 0.25 100 1000 SM 3 scatter u u~ -> g g +uugga 0.25 100 1000 SM 3 scatter u u~ -> g g a +uuagg 0.20 100 1000 SM 3 scatter u u~ -> a g g +uuggg 0.20 100 1000 SM 3 scatter u u~ -> g g g +gggg 0.50 100 1000 SM 3 scatter g g -> g g +ggggg 0.50 100 1000 SM 3 scatter g g -> g g g Index: trunk/omega/tests/comparisons.list =================================================================== --- trunk/omega/tests/comparisons.list (revision 8488) +++ trunk/omega/tests/comparisons.list (revision 8489) @@ -1,6 +1,6 @@ # comparisons.list -- # ---------------------------------------------------------------------- -# thr abs_thr n roots model process ... +# thr abs_thr n roots model process ... # ---------------------------------------------------------------------- -gggg 0.60 1E-11 1000 1000 QCD scatter gl gl -> gl gl -ggggg 0.60 1E-11 200 1000 QCD scatter gl gl -> gl gl gl +gggg 0.60 1E-11 100 1000 QCD scatter gl gl -> gl gl +ggggg 0.60 1E-11 50 1000 QCD scatter gl gl -> gl gl gl Index: trunk/omega/tests/ward_identities_long.list =================================================================== --- trunk/omega/tests/ward_identities_long.list (revision 8488) +++ trunk/omega/tests/ward_identities_long.list (revision 8489) @@ -1,10 +1,10 @@ # ward_identities_long.list -- # ---------------------------------------------------------------------- -# thr n roots model i process ... +# thr n roots model i process ... # ---------------------------------------------------------------------- -uuggg 0.7 500 1000 QCD 3 scatter u ubar -> gl gl gl -uggggg 0.7 100 1000 QCD 3 scatter u ubar -> gl gl gl gl -ggggg 0.7 500 1000 QCD 3 scatter gl gl -> gl gl gl -gggggg 0.7 100 1000 QCD 3 scatter gl gl -> gl gl gl gl +uuggg 0.7 50 1000 QCD 3 scatter u ubar -> gl gl gl +uggggg 0.7 10 1000 QCD 3 scatter u ubar -> gl gl gl gl +ggggg 0.7 50 1000 QCD 3 scatter gl gl -> gl gl gl +gggggg 0.7 10 1000 QCD 3 scatter gl gl -> gl gl gl gl Index: trunk/omega/tests/comparisons_UFO.list =================================================================== --- trunk/omega/tests/comparisons_UFO.list (revision 8488) +++ trunk/omega/tests/comparisons_UFO.list (revision 8489) @@ -1,105 +1,105 @@ ### comparisons_UFO.list -- #! ### -------------------------------------------------------------------- -### thr abs n roots model mode process ... +### thr abs n roots model mode process ... ### -------------------------------------------------------------------- ### QCD -uudd 0.75 1E-11 1000 1000 SM scatter u ubar -> d dbar | u u~ -> d d~ -uucc 0.75 1E-11 1000 1000 SM scatter u ubar -> c cbar | u u~ -> c c~ -uuga 0.30 1E-11 1000 1000 SM scatter u ubar -> g A | u u~ -> g a -ddga 0.30 1E-11 1000 1000 SM scatter d dbar -> g A | d d~ -> g a -uugg 0.25 1E-11 1000 1000 SM scatter u ubar -> g g | u u~ -> g g -ddgg 0.25 1E-11 1000 1000 SM scatter d dbar -> g g | d d~ -> g g -uugga 0.30 1E-11 1000 1000 SM scatter u ubar -> g g A | u u~ -> g g a -uuggg 0.30 1E-11 1000 1000 SM scatter u ubar -> g g g | u u~ -> g g g +uudd 0.75 1E-11 100 1000 SM scatter u ubar -> d dbar | u u~ -> d d~ +uucc 0.75 1E-11 100 1000 SM scatter u ubar -> c cbar | u u~ -> c c~ +uuga 0.30 1E-11 100 1000 SM scatter u ubar -> g A | u u~ -> g a +ddga 0.30 1E-11 100 1000 SM scatter d dbar -> g A | d d~ -> g a +uugg 0.25 1E-11 100 1000 SM scatter u ubar -> g g | u u~ -> g g +ddgg 0.25 1E-11 100 1000 SM scatter d dbar -> g g | d d~ -> g g +uugga 0.30 1E-11 100 1000 SM scatter u ubar -> g g A | u u~ -> g g a +uuggg 0.30 1E-11 100 1000 SM scatter u ubar -> g g g | u u~ -> g g g ### -------------------------------------------------------------------- ### VB Higgs -wwhh 0.65 1E-11 1000 1000 SM scatter W+ W- -> H H -whwh 0.65 1E-11 1000 1000 SM scatter W+ H -> W+ H -hhww 0.70 1E-11 1000 1000 SM scatter H H -> W+ W- -zzhh 0.70 1E-11 1000 1000 SM scatter Z Z -> H H -zhzh 0.70 1E-11 1000 1000 SM scatter Z H -> Z H -hhzz 0.70 1E-11 1000 1000 SM scatter H H -> Z Z +wwhh 0.65 1E-11 100 1000 SM scatter W+ W- -> H H +whwh 0.65 1E-11 100 1000 SM scatter W+ H -> W+ H +hhww 0.70 1E-11 100 1000 SM scatter H H -> W+ W- +zzhh 0.70 1E-11 100 1000 SM scatter Z Z -> H H +zhzh 0.70 1E-11 100 1000 SM scatter Z H -> Z H +hhzz 0.70 1E-11 100 1000 SM scatter H H -> Z Z ### -------------------------------------------------------------------- ### Higgs -hhhh 0.75 1E-11 1000 1000 SM scatter H H -> H H -hhhhh 0.75 1E-11 500 1000 SM scatter H H -> H H H -hhhhhh 0.75 1E-11 200 1000 SM scatter H H -> H H H H +hhhh 0.75 1E-11 100 1000 SM scatter H H -> H H +hhhhh 0.75 1E-11 50 1000 SM scatter H H -> H H H +hhhhhh 0.75 1E-11 50 1000 SM scatter H H -> H H H H ### -------------------------------------------------------------------- ### VBS -wwww 0.50 1E-11 1000 1000 SM scatter W+ W- -> W+ W- -wwww_1 0.50 1E-11 1000 1000 SM scatter W+ W+ -> W+ W+ -wwww_2 0.50 1E-11 1000 1000 SM scatter W+ W- -> W- W+ -wwww_3 0.50 1E-11 1000 1000 SM scatter W- W- -> W- W- -aaww 0.55 1E-11 1000 1000 SM scatter A A -> W+ W- | a a -> W+ W- -awaw 0.55 1E-11 1000 1000 SM scatter A W- -> A W- | a W- -> a W- -awwa 0.55 1E-11 1000 1000 SM scatter A W- -> W- A | a W- -> W- a -wwaa 0.55 1E-11 1000 1000 SM scatter W+ W- -> A A | W+ W- -> a a -wwzz 0.50 1E-11 1000 1000 SM scatter W+ W- -> Z Z -wzwz 0.55 1E-11 1000 1000 SM scatter W+ Z -> W+ Z -wzwz_1 0.55 1E-11 1000 1000 SM scatter W- Z -> W- Z -wzzw 0.55 1E-11 1000 1000 SM scatter W+ Z -> Z W+ -wzzw_1 0.55 1E-11 1000 1000 SM scatter W+ Z -> W+ Z -wwaz 0.50 1E-11 1000 1000 SM scatter W+ W- -> A Z | W+ W- -> a Z -wwza 0.50 1E-11 1000 1000 SM scatter W+ W- -> Z A | W+ W- -> Z a -zaww 0.55 1E-11 1000 1000 SM scatter A Z -> W+ W- | a Z -> W+ W- -azww 0.55 1E-11 1000 1000 SM scatter Z A -> W+ W- | Z a -> W+ W- -wawz 0.55 1E-11 1000 1000 SM scatter W+ A -> W+ Z | W+ a -> W+ Z -wazw 0.55 1E-11 1000 1000 SM scatter W+ A -> Z W+ | W+ a -> Z W+ -wzwa 0.55 1E-11 1000 1000 SM scatter W+ Z -> W+ A | W+ Z -> W+ a -wzaw 0.55 1E-11 1000 1000 SM scatter W+ Z -> A W+ | W+ Z -> a W+ +wwww 0.50 1E-11 100 1000 SM scatter W+ W- -> W+ W- +wwww_1 0.50 1E-11 100 1000 SM scatter W+ W+ -> W+ W+ +wwww_2 0.50 1E-11 100 1000 SM scatter W+ W- -> W- W+ +wwww_3 0.50 1E-11 100 1000 SM scatter W- W- -> W- W- +aaww 0.55 1E-11 100 1000 SM scatter A A -> W+ W- | a a -> W+ W- +awaw 0.55 1E-11 100 1000 SM scatter A W- -> A W- | a W- -> a W- +awwa 0.55 1E-11 100 1000 SM scatter A W- -> W- A | a W- -> W- a +wwaa 0.55 1E-11 100 1000 SM scatter W+ W- -> A A | W+ W- -> a a +wwzz 0.50 1E-11 100 1000 SM scatter W+ W- -> Z Z +wzwz 0.55 1E-11 100 1000 SM scatter W+ Z -> W+ Z +wzwz_1 0.55 1E-11 100 1000 SM scatter W- Z -> W- Z +wzzw 0.55 1E-11 100 1000 SM scatter W+ Z -> Z W+ +wzzw_1 0.55 1E-11 100 1000 SM scatter W+ Z -> W+ Z +wwaz 0.50 1E-11 100 1000 SM scatter W+ W- -> A Z | W+ W- -> a Z +wwza 0.50 1E-11 100 1000 SM scatter W+ W- -> Z A | W+ W- -> Z a +zaww 0.55 1E-11 100 1000 SM scatter A Z -> W+ W- | a Z -> W+ W- +azww 0.55 1E-11 100 1000 SM scatter Z A -> W+ W- | Z a -> W+ W- +wawz 0.55 1E-11 100 1000 SM scatter W+ A -> W+ Z | W+ a -> W+ Z +wazw 0.55 1E-11 100 1000 SM scatter W+ A -> Z W+ | W+ a -> Z W+ +wzwa 0.55 1E-11 100 1000 SM scatter W+ Z -> W+ A | W+ Z -> W+ a +wzaw 0.55 1E-11 100 1000 SM scatter W+ Z -> A W+ | W+ Z -> a W+ ### -------------------------------------------------------------------- ### NC -eemm 0.60 1E-11 1000 1000 SM scatter e+ e- -> mu+ mu- -emem 0.60 1E-11 1000 1000 SM scatter e- mu- -> e- mu- -eemmmm 0.60 1E-11 200 1000 SM scatter e+ e- -> mu+ mu- mu+ mu- -eeee 0.60 1E-11 1000 1000 SM scatter e+ e- -> e+ e- -eeee_2 0.60 1E-11 1000 1000 SM scatter e- e- -> e- e- -eeee_3 0.60 1E-11 1000 1000 SM scatter e+ e+ -> e+ e+ -eeeeee 0.60 1E-11 200 1000 SM scatter e+ e- -> e+ e- e+ e- -eaea 0.60 1E-11 1000 1000 SM scatter e- A -> e- A | e- a -> e- a -eeaa 0.60 1E-11 1000 1000 SM scatter e+ e- -> A A | e+ e- -> a a -aaee 0.60 1E-11 1000 1000 SM scatter A A -> e+ e- | a a -> e+ e- -eeaaa 0.60 1E-11 500 1000 SM scatter e+ e- -> A A A | e+ e- -> a a a -eeaaaa 0.60 1E-11 200 1000 SM scatter e+ e- -> A A A A | e+ e- -> a a a a -ezez 0.25 1E-11 1000 1000 SM scatter e- Z -> e- Z -eezz 0.30 1E-11 1000 1000 SM scatter e+ e- -> Z Z -zzee 0.25 1E-11 1000 1000 SM scatter Z Z -> e+ e- -eaez 0.25 1E-11 1000 1000 SM scatter e- A -> e- Z | e- a -> e- Z -eeaz 0.25 1E-11 1000 1000 SM scatter e+ e- -> A Z | e+ e- -> a Z -azee 0.25 1E-11 1000 1000 SM scatter A Z -> e+ e- | a Z -> e+ e- -enuenu 0.60 1E-11 200 1000 SM scatter e- numu -> e- numu | e- vm -> e- vm -eenunu 0.60 1E-11 200 1000 SM scatter e+ e- -> nue nuebar | e+ e- -> ve ve~ -nu4 0.60 1E-11 1000 1000 SM scatter nue nue -> nue nue | ve ve -> ve ve -eedd 0.60 1E-11 1000 1000 SM scatter e+ e- -> d dbar | e+ e- -> d d~ -eded 0.60 1E-11 1000 1000 SM scatter e+ d -> e+ d | e+ d -> e+ d -eeuu 0.60 1E-11 1000 1000 SM scatter e+ e- -> u ubar | e+ e- -> u u~ -eueu 0.60 1E-11 1000 1000 SM scatter e+ u -> e+ u | e+ u -> e+ u -uuee 0.60 1E-11 1000 1000 SM scatter u ubar -> e+ e-| u u~ -> e+ e- -uuaa 0.45 1E-11 1000 1000 SM scatter u ubar -> A A | u u~ -> a a -uuaz 0.25 1E-11 1000 1000 SM scatter u ubar -> A Z | u u~ -> a Z +eemm 0.60 1E-11 100 1000 SM scatter e+ e- -> mu+ mu- +emem 0.60 1E-11 100 1000 SM scatter e- mu- -> e- mu- +eemmmm 0.60 1E-11 50 1000 SM scatter e+ e- -> mu+ mu- mu+ mu- +eeee 0.60 1E-11 100 1000 SM scatter e+ e- -> e+ e- +eeee_2 0.60 1E-11 100 1000 SM scatter e- e- -> e- e- +eeee_3 0.60 1E-11 100 1000 SM scatter e+ e+ -> e+ e+ +eeeeee 0.60 1E-11 50 1000 SM scatter e+ e- -> e+ e- e+ e- +eaea 0.60 1E-11 100 1000 SM scatter e- A -> e- A | e- a -> e- a +eeaa 0.60 1E-11 100 1000 SM scatter e+ e- -> A A | e+ e- -> a a +aaee 0.60 1E-11 100 1000 SM scatter A A -> e+ e- | a a -> e+ e- +eeaaa 0.60 1E-11 50 1000 SM scatter e+ e- -> A A A | e+ e- -> a a a +eeaaaa 0.60 1E-11 50 1000 SM scatter e+ e- -> A A A A | e+ e- -> a a a a +ezez 0.25 1E-11 100 1000 SM scatter e- Z -> e- Z +eezz 0.30 1E-11 100 1000 SM scatter e+ e- -> Z Z +zzee 0.25 1E-11 100 1000 SM scatter Z Z -> e+ e- +eaez 0.25 1E-11 100 1000 SM scatter e- A -> e- Z | e- a -> e- Z +eeaz 0.25 1E-11 100 1000 SM scatter e+ e- -> A Z | e+ e- -> a Z +azee 0.25 1E-11 100 1000 SM scatter A Z -> e+ e- | a Z -> e+ e- +enuenu 0.60 1E-11 50 1000 SM scatter e- numu -> e- numu | e- vm -> e- vm +eenunu 0.60 1E-11 50 1000 SM scatter e+ e- -> nue nuebar | e+ e- -> ve ve~ +nu4 0.60 1E-11 100 1000 SM scatter nue nue -> nue nue | ve ve -> ve ve +eedd 0.60 1E-11 100 1000 SM scatter e+ e- -> d dbar | e+ e- -> d d~ +eded 0.60 1E-11 100 1000 SM scatter e+ d -> e+ d | e+ d -> e+ d +eeuu 0.60 1E-11 100 1000 SM scatter e+ e- -> u ubar | e+ e- -> u u~ +eueu 0.60 1E-11 100 1000 SM scatter e+ u -> e+ u | e+ u -> e+ u +uuee 0.60 1E-11 100 1000 SM scatter u ubar -> e+ e-| u u~ -> e+ e- +uuaa 0.45 1E-11 100 1000 SM scatter u ubar -> A A | u u~ -> a a +uuaz 0.25 1E-11 100 1000 SM scatter u ubar -> A Z | u u~ -> a Z ### -------------------------------------------------------------------- ### CC -enumunu 0.80 1E-11 200 1000 SM scatter e- nuebar -> mu- numubar | e- ve~ -> mu- vm~ -emununu 0.75 1E-11 200 1000 SM scatter e- mu+ -> nue numubar | e- mu+ -> ve vm~ -uuww 0.35 1E-11 200 1000 SM scatter u ubar -> W+ W- | u u~ -> W+ W- +enumunu 0.80 1E-11 50 1000 SM scatter e- nuebar -> mu- numubar | e- ve~ -> mu- vm~ +emununu 0.75 1E-11 50 1000 SM scatter e- mu+ -> nue numubar | e- mu+ -> ve vm~ +uuww 0.35 1E-11 50 1000 SM scatter u ubar -> W+ W- | u u~ -> W+ W- ### -------------------------------------------------------------------- ### CC/NC interferences -enuwa 0.60 1E-11 500 1000 SM scatter e- nuebar -> W- A | e- ve~ -> W- a -ewnua 0.75 1E-11 500 1000 SM scatter e- W+ -> nue A | e- W+ -> ve a -wenua 0.75 1E-11 500 1000 SM scatter e+ W- -> nuebar A | e+ W- -> ve~ a -enuwz 0.60 1E-11 500 1000 SM scatter e- nuebar -> W- Z | e- ve~ -> W- Z -ewnuz 0.75 1E-11 500 1000 SM scatter e- W+ -> nue Z | e- W+ -> ve Z +enuwa 0.60 1E-11 50 1000 SM scatter e- nuebar -> W- A | e- ve~ -> W- a +ewnua 0.75 1E-11 50 1000 SM scatter e- W+ -> nue A | e- W+ -> ve a +wenua 0.75 1E-11 50 1000 SM scatter e+ W- -> nuebar A | e+ W- -> ve~ a +enuwz 0.60 1E-11 50 1000 SM scatter e- nuebar -> W- Z | e- ve~ -> W- Z +ewnuz 0.75 1E-11 50 1000 SM scatter e- W+ -> nue Z | e- W+ -> ve Z ### -------------------------------------------------------------------- ### Top Yukawa -tthh 0.80 1E-11 1000 1000 SM scatter t tbar -> H H | t t~ -> H H -ttgh 0.80 1E-11 1000 1000 SM scatter t tbar -> g H | t t~ -> g H -tbwh 0.80 1E-11 1000 1000 SM scatter t bbar -> W+ H | t b~ -> W+ H +tthh 0.80 1E-11 100 1000 SM scatter t tbar -> H H | t t~ -> H H +ttgh 0.80 1E-11 100 1000 SM scatter t tbar -> g H | t t~ -> g H +tbwh 0.80 1E-11 100 1000 SM scatter t bbar -> W+ H | t b~ -> W+ H ### -------------------------------------------------------------------- ### Pure QCD -ggg 0.70 1E-11 1000 1000 SM decay g -> g g -gggg 0.60 1E-11 500 1000 SM scatter g g -> g g -ggggg 0.30 1E-11 200 1000 SM scatter g g -> g g g +ggg 0.70 1E-11 100 1000 SM decay g -> g g +gggg 0.60 1E-11 50 1000 SM scatter g g -> g g +ggggg 0.30 1E-11 50 1000 SM scatter g g -> g g g ### Other -eenu2H 0.75 1E-11 500 1000 SM scatter e+ e- -> nue nuebar H | e+ e- -> ve ve~ H +eenu2H 0.75 1E-11 50 1000 SM scatter e+ e- -> nue nuebar H | e+ e- -> ve ve~ H ### ---------------------------------------------------------------------- Index: trunk/omega/tests/comparisons_majorana_legacy.list =================================================================== --- trunk/omega/tests/comparisons_majorana_legacy.list (revision 8488) +++ trunk/omega/tests/comparisons_majorana_legacy.list (revision 8489) @@ -1,23 +1,23 @@ # comparisons_majorana_legacy.list -- ######################################################################## # NB: the *overall* sign of an apmplitude may change depending on wether # we compute with Fusion.Mixed23 or Fusion_vintage.Mixed23_Majorana ######################################################################## # ---------------------------------------------------------------------- -# thr abs_thr n roots model sign process ... +# thr abs_thr n roots model sign process ... # ---------------------------------------------------------------------- -eeee 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- -eeeea 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- A -eeeeaa 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- A A -eeeeep 0.80 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- e+ e- -aaep 0.90 1E-11 1000 1000 SM .false. scatter A A -> e+ e- -aaepep 0.75 1E-11 1000 1000 SM .false. scatter A A -> e+ e- e+ e- -eaea 0.90 1E-11 1000 1000 SM .false. scatter e- A -> e- A -epaa 0.90 1E-11 1000 1000 SM .true. scatter e+ e- -> A A -epep 0.90 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- -epepa 0.75 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- A -epepaa 0.80 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- A A -epepep 0.80 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- -epepepa 0.70 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A -epepepaa 0.70 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A A -epepepep 0.70 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- e+ e- +eeee 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- +eeeea 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- A +eeeeaa 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- A A +eeeeep 0.80 1E-11 100 1000 SM .true. scatter e- e- -> e- e- e+ e- +aaep 0.90 1E-11 100 1000 SM .false. scatter A A -> e+ e- +aaepep 0.75 1E-11 100 1000 SM .false. scatter A A -> e+ e- e+ e- +eaea 0.90 1E-11 100 1000 SM .false. scatter e- A -> e- A +epaa 0.90 1E-11 100 1000 SM .true. scatter e+ e- -> A A +epep 0.90 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- +epepa 0.75 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- A +epepaa 0.80 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- A A +epepep 0.80 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- +epepepa 0.70 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A +epepepaa 0.70 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A A +epepepep 0.70 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- e+ e- Index: trunk/omega/tests/keystones_UFO_bispinors_generate.ml =================================================================== --- trunk/omega/tests/keystones_UFO_bispinors_generate.ml (revision 8488) +++ trunk/omega/tests/keystones_UFO_bispinors_generate.ml (revision 8489) @@ -1,155 +1,155 @@ (* keystones_UFO_generate.ml -- Copyright (C) 2019-2020 by Wolfgang Kilian Thorsten Ohl Juergen Reuter WHIZARD is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2, or (at your option) any later version. WHIZARD is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. *) open Coupling open Keystones let qed = equivalent_tensors [| Majorana; Vector; Majorana |] [ ("qed", "Gamma(2,1,3)") ] let axial = equivalent_tensors [| Majorana; Vector; Majorana |] [ ("axial1", "Gamma5(1,-1)*Gamma(2,-1,3)"); ("axial2", "-Gamma(2,1,-3)*Gamma5(-3,3)") ] let left = equivalent_tensors [| Majorana; Vector; Majorana |] [ ("left1", "(Identity(1,-1)+Gamma5(1,-1))*Gamma(2,-1,3)"); ("left2", "2*ProjP(1,-1)*Gamma(2,-1,3)"); ("left3", "Gamma(2,1,-3)*(Identity(-3,3)-Gamma5(-3,3))"); ("left4", "2*Gamma(2,1,-3)*ProjM(-3,3)") ] let right = equivalent_tensors [| Majorana; Vector; Majorana |] [ ("right1", "(Identity(1,-1)-Gamma5(1,-1))*Gamma(2,-1,3)"); ("right2", "2*ProjM(1,-1)*Gamma(2,-1,3)"); ("right3", "Gamma(2,1,-3)*(Identity(-3,3)+Gamma5(-3,3))"); ("right4", "2*Gamma(2,1,-3)*ProjP(-3,3)") ] let vector_spinor_current tag = { tag = Printf.sprintf "vector_spinor_current__%s_ff" tag; keystones = [ { bra = (Majorana, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Vector, 1); F (Majorana, 2)] }; { bra = (Vector, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (Majorana, 0); F (Majorana, 2)] } ] } let scalar = equivalent_tensors [| Majorana; Scalar; Majorana |] [ ("scalar_current", "Identity(1,3)") ] let pseudo = equivalent_tensors [| Majorana; Scalar; Majorana |] [ ("pseudo_current", "Gamma5(1,3)") ] let left_scalar = equivalent_tensors [| Majorana; Scalar; Majorana |] [ ("left_scalar1", "Identity(1,3)-Gamma5(1,3)"); ("left_scalar2", "2*ProjM(1,3)") ] let right_scalar = equivalent_tensors [| Majorana; Scalar; Majorana |] [ ("right_scalar1", "Identity(1,3)+Gamma5(1,3)"); ("right_scalar2", "2*ProjP(1,3)") ] let scalar_spinor_current tag = { tag = Printf.sprintf "scalar_spinor_current__%s_ff" tag; keystones = [ { bra = (Majorana, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Scalar, 1); F (Majorana, 2)] }; { bra = (Scalar, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (Majorana, 0); F (Majorana, 2)] } ] } let empty = { tag = "empty"; keystones = [ ] } let vertices = [ (qed, vector_spinor_current "v"); (axial, vector_spinor_current "a"); (left, vector_spinor_current "vl"); (right, vector_spinor_current "vr"); (scalar, scalar_spinor_current "s"); (pseudo, scalar_spinor_current "p"); (left_scalar, scalar_spinor_current "sl"); (right_scalar, scalar_spinor_current "sr"); ] let parse_propagator (p_tag, p_omega, p_spins, numerator, denominator) = let p = UFO.Propagator.of_propagator_UFO ~majorana:true { UFO.Propagator_UFO.name = p_tag; UFO.Propagator_UFO.numerator = UFOx.Lorentz.of_string numerator; UFO.Propagator_UFO.denominator = UFOx.Lorentz.of_string denominator } in { p_tag; p_omega; p_spins; p_propagator = p } let default_denominator = "P('mu', id) * P('mu', id) - Mass(id) * Mass(id) \ + complex(0,1) * Mass(id) * Width(id)" let majorana_propagator = ( "majorana", "pr_psi", (Majorana, Majorana), "Gamma('mu', 1, 2) * P('mu', id) + Mass(id) * Identity(1, 2)", default_denominator ) let gravitino_propagator = ( "vectorspinor", "pr_grav", (Vectorspinor, Vectorspinor), "(Gamma(-1,1,2)*P(-1,id) - Mass(id)*Identity(1,2)) \ * (Metric(1,2) - P(1,id)*P(2,id)/Mass(id)**2) \ + 1/3 * (Gamma(1,1,-1) - P(1,id)/Mass(id)*Identity(1,-1)) \ * (Gamma(-3,-1,-2)*P(-3,id) + Mass(id)*Identity(-1,-2)) \ * (Gamma(2,-2,2) - P(2,id)/Mass(id)*Identity(-2,2)) ", default_denominator ) let gravitino_propagator = ( "vectorspinor", "pr_grav", (Vectorspinor, Vectorspinor), "(Gamma(-1,2001,2002)*P(-1,id) - Mass(id)*Identity(2001,2002)) \ * (Metric(1001,1002) - P(1001,id)*P(1002,id)/Mass(id)**2) \ + 1/3 * (Gamma(1001,2001,-1) - P(1001,id)/Mass(id)*Identity(2001,-1)) \ * (Gamma(-3,-1,-2)*P(-3,id) + Mass(id)*Identity(-1,-2)) \ * (Gamma(1002,-2,2002) - P(1002,id)/Mass(id)*Identity(-2,2002)) ", default_denominator ) let propagators = List.map parse_propagator [ majorana_propagator; (* [gravitino_propagator] *) ] let _ = generate_ufo - ~reps:10000 ~threshold:0.70 ~omega_module:"omega95_bispinors" + ~reps:1000 ~threshold:0.70 ~omega_module:"omega95_bispinors" "fusions_UFO_bispinors" vertices propagators; exit 0 Index: trunk/omega/tests/comparisons_VM.list =================================================================== --- trunk/omega/tests/comparisons_VM.list (revision 8488) +++ trunk/omega/tests/comparisons_VM.list (revision 8489) @@ -1,93 +1,93 @@ # processes to cross-check the OVM with compiled code # # For these comparisons, use the relative threshold `thr` except for amplitudes # smaller than `abs_thr` with `n` points for a fixed `roots` and decay or # scattering `process` in their respective `model` # # ---------------------------------------------------------------------- # thr abs_thr n roots model process # ---------------------------------------------------------------------- # Check QCD and color flows -gggg 0.30 1E-12 200 1000 QCD scatter gl gl -> gl gl -ggggg 0.30 1E-12 100 1000 QCD scatter gl gl -> gl gl gl -uubddb 0.30 1E-15 200 1000 QCD scatter u ubar -> d dbar -udud 0.60 1E-15 200 1000 QCD scatter u d -> u d -gguub 0.60 1E-15 200 1000 QCD scatter gl gl -> u ubar -uubgg 0.60 1E-15 200 1000 QCD scatter u ubar -> gl gl +gggg 0.30 1E-12 100 1000 QCD scatter gl gl -> gl gl +ggggg 0.30 1E-12 50 1000 QCD scatter gl gl -> gl gl gl +uubddb 0.30 1E-15 100 1000 QCD scatter u ubar -> d dbar +udud 0.60 1E-15 100 1000 QCD scatter u d -> u d +gguub 0.60 1E-15 100 1000 QCD scatter gl gl -> u ubar +uubgg 0.60 1E-15 100 1000 QCD scatter u ubar -> gl gl # Check multiple level parallelization with simple QED -eeaa 0.60 1E-15 200 1000 QED scatter e+ e- -> A A -ee4a 0.40 1E-15 200 50 QED scatter e+ e- -> A A A A -eeaee 0.60 1E-15 200 500 QED scatter e+ e- -> A e+ e- -aeae 0.60 1E-15 200 1000 QED scatter A e+ -> A e+ -aeaea 0.60 1E-15 200 50 QED scatter A e+ -> A e+ A +eeaa 0.60 1E-15 100 1000 QED scatter e+ e- -> A A +ee4a 0.40 1E-15 100 50 QED scatter e+ e- -> A A A A +eeaee 0.60 1E-15 100 500 QED scatter e+ e- -> A e+ e- +aeae 0.60 1E-15 100 1000 QED scatter A e+ -> A e+ +aeaea 0.60 1E-15 100 50 QED scatter A e+ -> A e+ A # Check some SM vertices -aatt 0.60 1E-15 200 1000 SM scatter t tbar -> A A -ttaa 0.60 1E-15 200 1000 SM scatter A A -> t tbar -tata 0.60 1E-15 200 1000 SM scatter A t -> A t -atat 0.60 1E-15 200 1000 SM scatter A tbar -> A tbar +aatt 0.60 1E-15 100 1000 SM scatter t tbar -> A A +ttaa 0.60 1E-15 100 1000 SM scatter A A -> t tbar +tata 0.60 1E-15 100 1000 SM scatter A t -> A t +atat 0.60 1E-15 100 1000 SM scatter A tbar -> A tbar # Check multi-flavor algorithm 2uubgg 0.60 1E-12 200 1000 QCD scatter u:ubar u:ubar -> gl gl 4qj 0.60 1E-12 100 1000 QCD scatter u:ubar:d:dbar u:ubar:d:dbar -> u:ubar:d:dbar u:ubar:d:dbar eeall 0.60 1E-12 200 50 QED scatter e+ e- -> e+:e-:m+:m-:A e+:e-:m+:m-:A e+:e-:m+:m-:A # Check V-A couplings DY 0.60 1E-15 200 1000 SM scatter u ubar -> e+ e- DYj 0.60 1E-15 200 1000 SM scatter u ubar -> e+ e- u:ubar:gl # Check scalar couplings and 1 -> 2, 2 -> 1 -SSS 0.60 1E-15 200 1000 SM decay H -> H H -3Sx 0.60 1E-15 1 1000 SM scatter H H -> H -4H 0.60 1E-15 200 1000 SM scatter H H -> H H -5H 0.60 1E-15 200 1000 SM scatter H H -> H H H -6H 0.60 1E-15 200 50 SM scatter H H -> H H H H +SSS 0.60 1E-15 100 1000 SM decay H -> H H +3Sx 0.60 1E-15 100 1000 SM scatter H H -> H +4H 0.60 1E-15 100 1000 SM scatter H H -> H H +5H 0.60 1E-15 100 1000 SM scatter H H -> H H H +6H 0.60 1E-15 100 50 SM scatter H H -> H H H H # Check scalar-vector couplings -SVV 0.60 1E-15 200 1000 SM decay H -> W+ W- -VSV 0.60 1E-15 200 1000 SM decay Z -> Z H -VSSV 0.60 1E-15 200 1000 SM scatter W+ W- -> H H -SSVV 0.60 1E-15 200 1000 SM scatter H H -> W+ W- -WpWm 0.60 1E-15 200 1000 SM scatter W+ W- -> W+ W- -ZZZZ 0.60 1E-15 200 1000 SM scatter Z Z -> Z Z -VBF 0.60 1E-12 100 1000 SM scatter u:ubar:d:dbar u:ubar:d:dbar -> u:ubar:d:dbar u:ubar:d:dbar H +SVV 0.60 1E-15 100 1000 SM decay H -> W+ W- +VSV 0.60 1E-15 100 1000 SM decay Z -> Z H +VSSV 0.60 1E-15 100 1000 SM scatter W+ W- -> H H +SSVV 0.60 1E-15 100 1000 SM scatter H H -> W+ W- +WpWm 0.60 1E-15 100 1000 SM scatter W+ W- -> W+ W- +ZZZZ 0.60 1E-15 100 1000 SM scatter Z Z -> Z Z +VBF 0.60 1E-12 50 1000 SM scatter u:ubar:d:dbar u:ubar:d:dbar -> u:ubar:d:dbar u:ubar:d:dbar H # Check zero amplitude -#uudd 0.60 1E-12 200 1000 QCD scatter u u -> d d +#uudd 0.60 1E-12 100 1000 QCD scatter u u -> d d # Check SM_CKM -ttbar 0.60 1E-12 100 50 SM_CKM scatter u ubar -> e+ nue b bbar mu- numubar -pidec 0.60 1E-12 200 1000 SM_CKM scatter u dbar -> mu+ numu -betdec 0.60 1E-12 200 1000 SM_CKM decay d -> u e- nuebar -Bdec 0.60 1E-12 200 1000 SM_CKM decay b -> c e- nuebar -scat 0.60 1E-12 200 1000 SM_CKM scatter u g -> d mu+ numu +ttbar 0.60 1E-12 50 50 SM_CKM scatter u ubar -> e+ nue b bbar mu- numubar +pidec 0.60 1E-12 100 1000 SM_CKM scatter u dbar -> mu+ numu +betdec 0.60 1E-12 100 1000 SM_CKM decay d -> u e- nuebar +Bdec 0.60 1E-12 100 1000 SM_CKM decay b -> c e- nuebar +scat 0.60 1E-12 100 1000 SM_CKM scatter u g -> d mu+ numu # Check SM_Higgs -ggH 0.60 1E-12 1 1000 SM_Higgs scatter g g -> H -Hgg 0.60 1E-12 200 1000 SM_Higgs decay H -> g g -ggAA 0.60 1E-12 200 1000 SM_Higgs scatter g g -> A A -ZHgg 0.60 1E-12 200 1000 SM_Higgs scatter Z A -> g g -ggAAH 0.60 1E-12 200 1000 SM_Higgs scatter g g -> A A H -ggAAHH 0.60 1E-12 200 1000 SM_Higgs scatter g g -> A A H H +ggH 0.60 1E-12 100 1000 SM_Higgs scatter g g -> H +Hgg 0.60 1E-12 100 1000 SM_Higgs decay H -> g g +ggAA 0.60 1E-12 100 1000 SM_Higgs scatter g g -> A A +ZHgg 0.60 1E-12 100 1000 SM_Higgs scatter Z A -> g g +ggAAH 0.60 1E-12 100 1000 SM_Higgs scatter g g -> A A H +ggAAHH 0.60 1E-12 100 1000 SM_Higgs scatter g g -> A A H H # Check Zprime -eeuu 0.60 1E-15 200 10000 Zprime scatter e+ e- -> u ubar -Zhuub 0.60 1E-15 200 10000 Zprime decay Zh -> u ubar -ZhddbA 0.60 1E-15 200 10000 Zprime decay Zh -> d dbar A +eeuu 0.60 1E-15 100 10000 Zprime scatter e+ e- -> u ubar +Zhuub 0.60 1E-15 100 10000 Zprime decay Zh -> u ubar +ZhddbA 0.60 1E-15 100 10000 Zprime decay Zh -> d dbar A # Check THDM -eeWW 0.60 1E-15 200 10000 THDM scatter e+ e- -> W+ W- -eeWWx 0.60 1E-15 200 10000 THDM scatter e+ W- -> e+ W- -eeZZ 0.60 1E-15 200 10000 THDM scatter e+ e- -> Z Z -eehH 0.60 1E-15 200 10000 THDM scatter e+ e- -> h0 H0 -eehA 0.60 1E-15 200 10000 THDM scatter e+ e- -> h0 A0 -WWZZZ 0.30 1E-12 200 10000 THDM scatter W+ W- -> Z Z Z +eeWW 0.60 1E-15 100 10000 THDM scatter e+ e- -> W+ W- +eeWWx 0.60 1E-15 100 10000 THDM scatter e+ W- -> e+ W- +eeZZ 0.60 1E-15 100 10000 THDM scatter e+ e- -> Z Z +eehH 0.60 1E-15 100 10000 THDM scatter e+ e- -> h0 H0 +eehA 0.60 1E-15 100 10000 THDM scatter e+ e- -> h0 A0 +WWZZZ 0.30 1E-12 100 10000 THDM scatter W+ W- -> Z Z Z # Check THDM_CKM -utbdbb 0.60 1E-15 1 1000 THDM_CKM scatter u tbar -> d bbar +utbdbb 0.60 1E-15 50 1000 THDM_CKM scatter u tbar -> d bbar # Check HSExt -ttbS 0.60 1E-15 1 1000 HSExt scatter t tbar -> S -ttbHS 0.60 1E-15 200 1000 HSExt scatter t tbar -> S H +ttbS 0.60 1E-15 50 1000 HSExt scatter t tbar -> S +ttbHS 0.60 1E-15 200 1000 HSExt scatter t tbar -> S H Index: trunk/omega/tests/fermi.list =================================================================== --- trunk/omega/tests/fermi.list (revision 8488) +++ trunk/omega/tests/fermi.list (revision 8489) @@ -1,29 +1,29 @@ # fermi.list -- # ---------------------------------------------------------------------- -# thr abs_thr n roots model i j eps process ... +# thr abs_thr n roots model i j eps process ... # ---------------------------------------------------------------------- -eeee 0.75 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- -eeee 0.80 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- -eeeea 0.65 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- A -eeeea 0.80 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- A -eeeeaa 0.72 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- A A -eeeeaa 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- A A -eeeeaa 0.75 1E-11 1000 1000 SM 5 6 1 scatter e- e- -> e- e- A A -eeeeee 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e+ e- -> e+ e+ e- e- -eeeeee 0.75 1E-11 1000 1000 SM 5 6 -1 scatter e+ e- -> e+ e+ e- e- -Seeee 0.70 1E-11 1000 1000 MSSM 1 2 -1 scatter e- e- -> e- e- -Seeee 0.75 1E-11 1000 1000 MSSM 3 4 -1 scatter e- e- -> e- e- -Seeeea 0.75 1E-11 1000 1000 MSSM 3 4 -1 scatter e- e- -> e- e- A -Seenn 0.75 1E-11 1000 1000 MSSM 3 4 -1 scatter e+ e- -> neu1 neu1 -Seen4 0.75 1E-11 1000 1000 MSSM 3 4 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 -Seen4 0.75 1E-11 1000 1000 MSSM 3 5 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 -Seen4 0.75 1E-11 1000 1000 MSSM 3 6 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 -Seen4 0.75 1E-11 1000 1000 MSSM 4 5 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 -Seen4 0.75 1E-11 1000 1000 MSSM 4 6 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 -Seen4 0.75 1E-11 1000 1000 MSSM 5 6 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 -Seess 0.75 1E-11 1000 1000 MSSM 1 2 -1 scatter e- e- -> se1- se1- -Seess 0.75 1E-11 1000 1000 MSSM 3 4 1 scatter e- e- -> se1- se1- -Sees4 0.75 1E-11 1000 1000 MSSM 3 4 1 scatter e+ e- -> se1+ se1+ se1- se1- -Sees4 0.75 1E-11 1000 1000 MSSM 5 6 1 scatter e+ e- -> se1+ se1+ se1- se1- +eeee 0.75 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- +eeee 0.80 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- +eeeea 0.65 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- A +eeeea 0.80 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- A +eeeeaa 0.72 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- A A +eeeeaa 0.75 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- A A +eeeeaa 0.75 1E-11 100 1000 SM 5 6 1 scatter e- e- -> e- e- A A +eeeeee 0.75 1E-11 100 1000 SM 3 4 -1 scatter e+ e- -> e+ e+ e- e- +eeeeee 0.75 1E-11 100 1000 SM 5 6 -1 scatter e+ e- -> e+ e+ e- e- +Seeee 0.70 1E-11 100 1000 MSSM 1 2 -1 scatter e- e- -> e- e- +Seeee 0.75 1E-11 100 1000 MSSM 3 4 -1 scatter e- e- -> e- e- +Seeeea 0.75 1E-11 100 1000 MSSM 3 4 -1 scatter e- e- -> e- e- A +Seenn 0.75 1E-11 100 1000 MSSM 3 4 -1 scatter e+ e- -> neu1 neu1 +Seen4 0.75 1E-11 100 1000 MSSM 3 4 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 +Seen4 0.75 1E-11 100 1000 MSSM 3 5 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 +Seen4 0.75 1E-11 100 1000 MSSM 3 6 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 +Seen4 0.75 1E-11 100 1000 MSSM 4 5 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 +Seen4 0.75 1E-11 100 1000 MSSM 4 6 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 +Seen4 0.75 1E-11 100 1000 MSSM 5 6 -1 scatter e+ e- -> neu1 neu1 neu1 neu1 +Seess 0.75 1E-11 100 1000 MSSM 1 2 -1 scatter e- e- -> se1- se1- +Seess 0.75 1E-11 100 1000 MSSM 3 4 1 scatter e- e- -> se1- se1- +Sees4 0.75 1E-11 100 1000 MSSM 3 4 1 scatter e+ e- -> se1+ se1+ se1- se1- +Sees4 0.75 1E-11 100 1000 MSSM 5 6 1 scatter e+ e- -> se1+ se1+ se1- se1- ! -uuuu 0.75 1E-11 1000 1000 SM 3 4 -1 scatter u u -> u u +uuuu 0.75 1E-11 100 1000 SM 3 4 -1 scatter u u -> u u Index: trunk/omega/tests/keystones_UFO_generate.ml =================================================================== --- trunk/omega/tests/keystones_UFO_generate.ml (revision 8488) +++ trunk/omega/tests/keystones_UFO_generate.ml (revision 8489) @@ -1,379 +1,379 @@ (* keystones_UFO_generate.ml -- Copyright (C) 2019-2020 by Wolfgang Kilian Thorsten Ohl Juergen Reuter WHIZARD is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2, or (at your option) any later version. WHIZARD is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. *) open Coupling open Keystones let qed = equivalent_tensors [| ConjSpinor; Vector; Spinor |] [ ("qed", "Gamma(2,1,3)") ] let axial = equivalent_tensors [| ConjSpinor; Vector; Spinor |] [ ("axial1", "Gamma5(1,-1)*Gamma(2,-1,3)"); ("axial2", "-Gamma(2,1,-3)*Gamma5(-3,3)") ] let left = equivalent_tensors [| ConjSpinor; Vector; Spinor |] [ ("left1", "(Identity(1,-1)+Gamma5(1,-1))*Gamma(2,-1,3)"); ("left2", "2*ProjP(1,-1)*Gamma(2,-1,3)"); ("left3", "Gamma(2,1,-3)*(Identity(-3,3)-Gamma5(-3,3))"); ("left4", "2*Gamma(2,1,-3)*ProjM(-3,3)") ] let right = equivalent_tensors [| ConjSpinor; Vector; Spinor |] [ ("right1", "(Identity(1,-1)-Gamma5(1,-1))*Gamma(2,-1,3)"); ("right2", "2*ProjM(1,-1)*Gamma(2,-1,3)"); ("right3", "Gamma(2,1,-3)*(Identity(-3,3)+Gamma5(-3,3))"); ("right4", "2*Gamma(2,1,-3)*ProjP(-3,3)") ] let vector_spinor_current tag = { tag = Printf.sprintf "vector_spinor_current__%s_ff" tag; keystones = [ { bra = (ConjSpinor, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Vector, 1); F (Spinor, 2)] }; { bra = (Vector, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (ConjSpinor, 0); F (Spinor, 2)] }; { bra = (Spinor, 2); name = Printf.sprintf "f_f%s" tag; args = [G (0); F (ConjSpinor, 0); F (Vector, 1)] } ] } let scalar = equivalent_tensors [| ConjSpinor; Scalar; Spinor |] [ ("scalar_current", "Identity(1,3)") ] let pseudo = equivalent_tensors [| ConjSpinor; Scalar; Spinor |] [ ("pseudo_current", "Gamma5(1,3)") ] let left_scalar = equivalent_tensors [| ConjSpinor; Scalar; Spinor |] [ ("left_scalar1", "Identity(1,3)-Gamma5(1,3)"); ("left_scalar2", "2*ProjM(1,3)") ] let right_scalar = equivalent_tensors [| ConjSpinor; Scalar; Spinor |] [ ("right_scalar1", "Identity(1,3)+Gamma5(1,3)"); ("right_scalar2", "2*ProjP(1,3)") ] let scalar_spinor_current tag = { tag = Printf.sprintf "scalar_spinor_current__%s_ff" tag; keystones = [ { bra = (ConjSpinor, 0); name = Printf.sprintf "f_%sf" tag; args = [G (0); F (Scalar, 1); F (Spinor, 2)] }; { bra = (Scalar, 1); name = Printf.sprintf "%s_ff" tag; args = [G (0); F (ConjSpinor, 0); F (Spinor, 2)] }; { bra = (Spinor, 2); name = Printf.sprintf "f_f%s" tag; args = [G (0); F (ConjSpinor, 0); F (Scalar, 1)] } ] } let fermi_ss = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_ss", "Identity(1,2)*Identity(3,4)"); ("fermi_ss_f", " (1/4) * Identity(1,4)*Identity(3,2)" ^ " + (1/4) * Gamma(-1,1,4)*Gamma(-1,3,2)" ^ " + (1/8) * Sigma(-1,-2,1,4)*Sigma(-1,-2,3,2)" ^ " - (1/4) * Gamma(-1,1,-4)*Gamma5(-4,4)*Gamma(-1,3,-2)*Gamma5(-2,2)" ^ " + (1/4) * Gamma5(1,4)*Gamma5(3,2)") ] let fermi_vv = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_vv", "Gamma(-1,1,2)*Gamma(-1,3,4)"); ("fermi_vv_f", " Identity(1,4)*Identity(3,2)" ^ " - (1/2) * Gamma(-1,1,4)*Gamma(-1,3,2)" ^ " - (1/2) * Gamma(-1,1,-4)*Gamma5(-4,4)*Gamma(-1,3,-2)*Gamma5(-2,2)" ^ " - Gamma5(1,4)*Gamma5(3,2)") ] let fermi_tt = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_tt1", " Sigma(-1,-2,1,2)*Sigma(-1,-2,3,4)"); ("fermi_tt2", " - Sigma(-1,-2,1,2)*Sigma(-2,-1,3,4)"); ("fermi_tt3", " - Sigma(-2,-1,1,2)*Sigma(-1,-2,3,4)"); ("fermi_tt_f", " 3 * Identity(1,4)*Identity(3,2)" ^ " - (1/2) * Sigma(-1,-2,1,4)*Sigma(-1,-2,3,2)" ^ " + 3 * Gamma5(1,4)*Gamma5(3,2)") ] let fermi_aa = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_aa", "Gamma5(1,-2)*Gamma(-1,-2,2)*Gamma5(3,-3)*Gamma(-1,-3,4)"); ("fermi_aa_f", " - Identity(1,4)*Identity(3,2)" ^ " - (1/2) * Gamma(-1,1,4)*Gamma(-1,3,2)" ^ " - (1/2) * Gamma(-1,1,-4)*Gamma5(-4,4)*Gamma(-1,3,-2)*Gamma5(-2,2)" ^ " + Gamma5(1,4)*Gamma5(3,2)") ] let fermi_pp = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_pp", "Gamma5(1,2)*Gamma5(3,4)"); ("fermi_pp_f", " (1/4) * Identity(1,4)*Identity(3,2)" ^ " - (1/4) * Gamma(-1,1,4)*Gamma(-1,3,2)" ^ " + (1/8) * Sigma(-1,-2,1,4)*Sigma(-1,-2,3,2)" ^ " + (1/4) * Gamma(-1,1,-4)*Gamma5(-4,4)*Gamma(-1,3,-2)*Gamma5(-2,2)" ^ " + (1/4) * Gamma5(1,4)*Gamma5(3,2)") ] let fermi_ll = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_ll", " Gamma(-1,1,-2)*ProjM(-2,2)*Gamma(-1,3,-4)*ProjM(-4,4)"); ("fermi_ll_f", " - Gamma(-1,1,-2)*ProjM(-2,4)*Gamma(-1,3,-4)*ProjM(-4,2)") ] let fermi_va = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_va", "Gamma(-1,1,2)*Gamma5(3,-3)*Gamma(-1,-3,4)") ] let fermi_av = equivalent_tensors [| ConjSpinor; Spinor; ConjSpinor; Spinor |] [ ("fermi_av", "Gamma5(1,-2)*Gamma(-1,-2,2)*Gamma(-1,3,4)") ] let sqed = equivalent_tensors [| Scalar; Vector; Scalar |] [ ("sqed1", "P(2,3)-P(2,1)"); ("sqed2", "2*P(2,3)+P(2,2)"); ("sqed3", "-P(2,2)-2*P(2,1)") ] let vector_scalar_current = { tag = "vector_scalar_current__v_ss"; keystones = [ { bra = (Vector, 1); name = "v_ss"; args = [G (0); F (Scalar, 2); P (2); F (Scalar, 0); P (0)] }; { bra = (Scalar, 0); name = "s_vs"; args = [G (0); F (Vector, 1); P (1); F (Scalar, 2); P (2)] } ] } let svv_t = equivalent_tensors [| Scalar; Vector; Vector |] [ ("svv_t", "P(-1,2)*P(-1,3)*Metric(2,3)-P(2,3)*P(3,2)") ] let scalar_vector_current tag = { tag = Printf.sprintf "transversal_vector_current__s_vv_%s" tag; keystones = [ { bra = (Scalar, 0); name = Printf.sprintf "s_vv_%s" tag; args = [G (0); F (Vector, 1); P (1); F (Vector, 2); P (2)] }; { bra = (Vector, 1); name = Printf.sprintf "v_sv_%s" tag; args = [G (0); F (Scalar, 0); P (0); F (Vector, 2); P (2)] } ] } let gauge = equivalent_tensors [| Vector; Vector; Vector |] [ ("gauge", " Metric(1,2)*P(3,1) - Metric(1,2)*P(3,2) \ + Metric(3,1)*P(2,3) - Metric(3,1)*P(2,1) \ + Metric(2,3)*P(1,2) - Metric(2,3)*P(1,3)") ] let gauge_omega = { tag = "g_gg"; keystones = [ { bra = (Vector, 0); name = "(0,1)*g_gg"; args = [G (0); F (Vector, 1); P (1); F (Vector, 2); P (2)] } ] } (* Note that $C^{-1}=-C$ for the charge conjugation matrix.*) let charge_conjugate_s = equivalent_tensors [| Scalar; ConjSpinor; Spinor |] [ ("gamma1", "Identity(2,3)"); ("gamma1_cc", "C(3,-3)*Identity(-3,-2)*(-C(-2,2))"); ("gamma1_cx", "C(3,-1)*(-C(-1,2))") ] (* $C \gamma_5 C^{-1} = \gamma_5^T$ *) let charge_conjugate_p = equivalent_tensors [| Scalar; ConjSpinor; Spinor |] [ ("gamma5", "Gamma5(2,3)"); ("gamma5_cc", "C(3,-3)*Gamma5(-3,-2)*(-C(-2,2))") ] (* $C \gamma_\mu C^{-1} = - \gamma_\mu^T$ *) let charge_conjugate_v = equivalent_tensors [| Vector; ConjSpinor; Spinor |] [ ("gamma_mu", "Gamma(1,2,3)"); ("gamma_mu_cc", "-C(3,-3)*Gamma(1,-3,-2)*(-C(-2,2))") ] (* $C \gamma_5\gamma_\mu C^{-1} = (\gamma_5\gamma_\mu)^T$ *) let charge_conjugate_a = equivalent_tensors [| Vector; ConjSpinor; Spinor |] [ ("gamma_5mu", "Gamma5(2,-2)*Gamma(1,-2,3)"); ("gamma_5mu_cc", "C(3,-3)*Gamma5(-3,-1)*Gamma(1,-1,-2)*(-C(-2,2))") ] (* $C \sigma_{\mu\nu} C^{-1} = - \sigma_{\mu\nu}^T$ *) let charge_conjugate_t = equivalent_tensors [| Vector; Vector; ConjSpinor; Spinor |] [ ("sigma_munu", "Sigma(1,2,3,4)"); ("sigma_munu_cc", "-C(4,-4)*Sigma(1,2,-4,-3)*(-C(-3,3))") ] (* $C \gamma_\mu \gamma_\nu C^{-1} = \gamma_\nu^T \gamma_\mu^T$ *) let charge_conjugate_vv = equivalent_tensors [| Vector; Vector; ConjSpinor; Spinor |] [ ("gamma_mu_nu", "Gamma(1,3,-1)*Gamma(2,-1,4)"); ("gamma_mu_nu_cc", "C(4,-4)*Gamma(2,-4,-1)*Gamma(1,-1,-3)*(-C(-3,3))") ] let empty = { tag = "empty"; keystones = [ ] } let vertices = [ (qed, vector_spinor_current "v"); (axial, vector_spinor_current "a"); (left, vector_spinor_current "vl"); (right, vector_spinor_current "vr"); (scalar, scalar_spinor_current "s"); (pseudo, scalar_spinor_current "p"); (left_scalar, scalar_spinor_current "sl"); (right_scalar, scalar_spinor_current "sr"); (sqed, vector_scalar_current); (fermi_ss, empty); (fermi_vv, empty); (fermi_tt, empty); (fermi_aa, empty); (fermi_pp, empty); (fermi_ll, empty); (fermi_va, empty); (fermi_av, empty); (svv_t, scalar_vector_current "t"); (gauge, gauge_omega); (charge_conjugate_s, empty); (charge_conjugate_p, empty); (charge_conjugate_v, empty); (charge_conjugate_a, empty); (charge_conjugate_t, empty); (charge_conjugate_vv, empty) ] let parse_propagator (p_tag, p_omega, p_spins, numerator, denominator) = let p = UFO.Propagator.of_propagator_UFO { UFO.Propagator_UFO.name = p_tag; UFO.Propagator_UFO.numerator = UFOx.Lorentz.of_string numerator; UFO.Propagator_UFO.denominator = UFOx.Lorentz.of_string denominator } in { p_tag; p_omega; p_spins; p_propagator = p } let default_denominator = "P('mu', id) * P('mu', id) - Mass(id) * Mass(id) \ + complex(0,1) * Mass(id) * Width(id)" let scalar_propagator = ( "scalar", "pr_phi", (Scalar, Scalar), "1", default_denominator ) let spinor_propagator = ( "spinor", "pr_psi", (ConjSpinor, Spinor), "Gamma('mu', 1, 2) * P('mu', id) + Mass(id) * Identity(1, 2)", default_denominator ) let conjspinor_propagator = ( "conjspinor", "pr_psibar", (ConjSpinor, Spinor), "Gamma('mu', 1, 2) * P('mu', id) + Mass(id) * Identity(1, 2)", default_denominator ) let feynman_propagator = ( "feynman", "pr_feynman", (Vector, Vector), " - Metric(1, 2)", "P('mu', id) * P('mu', id)" ) let gauge_propagator = ( "gauge_propagator", "pr_gauge", (Vector, Vector), " - Metric(1, 2) + (1 - 42) * P(1,id) * P(2,id) / " ^ "( P('mu', id) * P('mu', id) )", "P('mu', id) * P('mu', id)" ) let rxi_propagator = ( "rxi_propagator", "pr_rxi", (Vector, Vector), " - Metric(1, 2) + (1 - 42) * P(1,id) * P(2,id) / " ^ "( P('mu', id) * P('mu', id) - 42 * Mass(id)**2 )", default_denominator ) let unitarity_propagator = ( "unitarity", "pr_unitarity", (Massive_Vector, Massive_Vector), "- Metric(1, 2) + Metric(1,'mu')*P('mu', id)*P(2, id)/Mass(id)**2", default_denominator ) let tensor_propagator = ( "tensor", "pr_tensor", (Tensor_2, Tensor_2), " 1/2 * (Metric(1001,1002) - P(1001,id)*P(1002,id)/Mass(id)**2) \ * (Metric(2001,2002) - P(2001,id)*P(2002,id)/Mass(id)**2) \ + 1/2 * (Metric(1001,2002) - P(1001,id)*P(2002,id)/Mass(id)**2) \ * (Metric(2001,1002) - P(2001,id)*P(1002,id)/Mass(id)**2) \ - 1/3 * (Metric(1001,2001) - P(1001,id)*P(2001,id)/Mass(id)**2) \ * (Metric(1002,2002) - P(1002,id)*P(2002,id)/Mass(id)**2) ", default_denominator ) let tensor_propagator_51_52 = ( "tensor_51_52", "pr_tensor", (Tensor_2, Tensor_2), " 1/2 * (Metric( 1, 2) - P( 1,id)*P( 2,id)/Mass(id)**2) \ * (Metric(51,52) - P(51,id)*P(52,id)/Mass(id)**2) \ + 1/2 * (Metric( 1,52) - P( 1,id)*P(52,id)/Mass(id)**2) \ * (Metric(51, 2) - P(51,id)*P( 2,id)/Mass(id)**2) \ - 1/3 * (Metric( 1,51) - P( 1,id)*P(51,id)/Mass(id)**2) \ * (Metric( 2,52) - P( 2,id)*P(52,id)/Mass(id)**2) ", default_denominator ) let propagators = List.map parse_propagator [ scalar_propagator; spinor_propagator; feynman_propagator; gauge_propagator; rxi_propagator; unitarity_propagator; tensor_propagator; tensor_propagator_51_52 ] let conjugate_propagators = List.map (fun p -> transpose (parse_propagator p)) [ conjspinor_propagator ] let all_propagators = propagators @ conjugate_propagators let _ = generate_ufo - ~reps:10000 ~threshold:0.70 "fusions_UFO" vertices all_propagators; + ~reps:1000 ~threshold:0.70 "fusions_UFO" vertices all_propagators; exit 0 Index: trunk/omega/tests/comparisons_majorana.list =================================================================== --- trunk/omega/tests/comparisons_majorana.list (revision 8488) +++ trunk/omega/tests/comparisons_majorana.list (revision 8489) @@ -1,23 +1,23 @@ # comparisons_majorana.list -- ######################################################################## # NB: the *overall* sign of an apmplitude may change depending on whether # we compute with Fusion.Mixed23 or Fusion.Mixed23_Majorana ######################################################################## # ---------------------------------------------------------------------- -# thr abs_thr n roots model sign process ... +# thr abs_thr n roots model sign process ... # ---------------------------------------------------------------------- -eeee 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- -eeeea 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- A -eeeeaa 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- A A -eeeeep 0.80 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- e+ e- -aaep 0.90 1E-11 1000 1000 SM .false. scatter A A -> e+ e- -aaepep 0.75 1E-11 1000 1000 SM .false. scatter A A -> e+ e- e+ e- -eaea 0.90 1E-11 1000 1000 SM .false. scatter e- A -> e- A -epaa 0.90 1E-11 1000 1000 SM .true. scatter e+ e- -> A A -epep 0.90 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- -epepa 0.75 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- A -epepaa 0.80 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- A A -epepep 0.80 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- -epepepa 0.70 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A -epepepaa 0.70 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A A -epepepep 0.70 1E-11 1000 1000 SM .true. scatter e+ e- -> e+ e- e+ e- e+ e- +eeee 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- +eeeea 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- A +eeeeaa 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- A A +eeeeep 0.80 1E-11 100 1000 SM .true. scatter e- e- -> e- e- e+ e- +aaep 0.90 1E-11 100 1000 SM .false. scatter A A -> e+ e- +aaepep 0.75 1E-11 100 1000 SM .false. scatter A A -> e+ e- e+ e- +eaea 0.90 1E-11 100 1000 SM .false. scatter e- A -> e- A +epaa 0.90 1E-11 100 1000 SM .true. scatter e+ e- -> A A +epep 0.90 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- +epepa 0.75 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- A +epepaa 0.80 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- A A +epepep 0.80 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- +epepepa 0.70 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A +epepepaa 0.70 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- A A +epepepep 0.70 1E-11 100 1000 SM .true. scatter e+ e- -> e+ e- e+ e- e+ e- Index: trunk/omega/tests/fermi_UFO.list =================================================================== --- trunk/omega/tests/fermi_UFO.list (revision 8488) +++ trunk/omega/tests/fermi_UFO.list (revision 8489) @@ -1,30 +1,30 @@ # fermi_UFO.list -- # ---------------------------------------------------------------------- -# thr abs_thr n roots model i j eps process ... +# thr abs_thr n roots model i j eps process ... # ---------------------------------------------------------------------- -eeee 0.70 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- -eeee 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- -eeeea 0.60 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- a -eeeea 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- a -eeeeaa 0.70 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- a a -eeeeaa 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- a a -eeeeaa 0.75 1E-11 1000 1000 SM 5 6 1 scatter e- e- -> e- e- a a -eeeepe 0.70 1E-11 1000 1000 SM 1 2 -1 scatter e- e- -> e- e- e+ e- -eeeepe 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e- e- -> e- e- e+ e- -eeeepe 0.75 1E-11 1000 1000 SM 3 6 -1 scatter e- e- -> e- e- e+ e- -eeeepe 0.75 1E-11 1000 1000 SM 4 6 -1 scatter e- e- -> e- e- e+ e- -peppee 0.75 1E-11 1000 1000 SM 3 4 -1 scatter e+ e- -> e+ e+ e- e- -peppee 0.75 1E-11 1000 1000 SM 5 6 -1 scatter e+ e- -> e+ e+ e- e- -meeee 0.70 1E-11 1000 1000 SM/M 1 2 -1 scatter e- e- -> e- e- -meeee 0.75 1E-11 1000 1000 SM/M 3 4 -1 scatter e- e- -> e- e- -#eeeea 0.40 1E-11 1000 1000 SM/M 1 2 -1 scatter e- e- -> e- e- a -meeeea 0.70 1E-11 1000 1000 SM/M 3 4 -1 scatter e- e- -> e- e- a -#eeeeaa 0.50 1E-11 1000 1000 SM/M 1 2 -1 scatter e- e- -> e- e- a a -meeeeaa 0.75 1E-11 1000 1000 SM/M 3 4 -1 scatter e- e- -> e- e- a a -meeeeaa 0.75 1E-11 1000 1000 SM/M 5 6 1 scatter e- e- -> e- e- a a -meeeepe 0.70 1E-11 1000 1000 SM/M 1 2 -1 scatter e- e- -> e- e- e+ e- -meeeepe 0.75 1E-11 1000 1000 SM/M 3 4 -1 scatter e- e- -> e- e- e+ e- -#eeeepe 0.70 1E-11 1000 1000 SM/M 3 6 -1 scatter e- e- -> e- e- e+ e- -#eeeepe 0.75 1E-11 1000 1000 SM/M 4 6 -1 scatter e- e- -> e- e- e+ e- -mpeppee 0.75 1E-11 1000 1000 SM/M 3 4 -1 scatter e+ e- -> e+ e+ e- e- -mpeppee 0.75 1E-11 1000 1000 SM/M 5 6 -1 scatter e+ e- -> e+ e+ e- e- +eeee 0.70 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- +eeee 0.75 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- +eeeea 0.60 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- a +eeeea 0.75 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- a +eeeeaa 0.70 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- a a +eeeeaa 0.75 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- a a +eeeeaa 0.75 1E-11 100 1000 SM 5 6 1 scatter e- e- -> e- e- a a +eeeepe 0.70 1E-11 100 1000 SM 1 2 -1 scatter e- e- -> e- e- e+ e- +eeeepe 0.75 1E-11 100 1000 SM 3 4 -1 scatter e- e- -> e- e- e+ e- +eeeepe 0.75 1E-11 100 1000 SM 3 6 -1 scatter e- e- -> e- e- e+ e- +eeeepe 0.75 1E-11 100 1000 SM 4 6 -1 scatter e- e- -> e- e- e+ e- +peppee 0.75 1E-11 100 1000 SM 3 4 -1 scatter e+ e- -> e+ e+ e- e- +peppee 0.75 1E-11 100 1000 SM 5 6 -1 scatter e+ e- -> e+ e+ e- e- +meeee 0.70 1E-11 100 1000 SM/M 1 2 -1 scatter e- e- -> e- e- +meeee 0.75 1E-11 100 1000 SM/M 3 4 -1 scatter e- e- -> e- e- +#eeeea 0.40 1E-11 100 1000 SM/M 1 2 -1 scatter e- e- -> e- e- a +meeeea 0.70 1E-11 100 1000 SM/M 3 4 -1 scatter e- e- -> e- e- a +#eeeeaa 0.50 1E-11 100 1000 SM/M 1 2 -1 scatter e- e- -> e- e- a a +meeeeaa 0.75 1E-11 100 1000 SM/M 3 4 -1 scatter e- e- -> e- e- a a +meeeeaa 0.75 1E-11 100 1000 SM/M 5 6 1 scatter e- e- -> e- e- a a +meeeepe 0.70 1E-11 100 1000 SM/M 1 2 -1 scatter e- e- -> e- e- e+ e- +meeeepe 0.75 1E-11 100 1000 SM/M 3 4 -1 scatter e- e- -> e- e- e+ e- +#eeeepe 0.70 1E-11 100 1000 SM/M 3 6 -1 scatter e- e- -> e- e- e+ e- +#eeeepe 0.75 1E-11 100 1000 SM/M 4 6 -1 scatter e- e- -> e- e- e+ e- +mpeppee 0.75 1E-11 100 1000 SM/M 3 4 -1 scatter e+ e- -> e+ e+ e- e- +mpeppee 0.75 1E-11 100 1000 SM/M 5 6 -1 scatter e+ e- -> e+ e+ e- e- Index: trunk/omega/tests/comparisons_majorana_UFO.list =================================================================== --- trunk/omega/tests/comparisons_majorana_UFO.list (revision 8488) +++ trunk/omega/tests/comparisons_majorana_UFO.list (revision 8489) @@ -1,26 +1,26 @@ # comparisons_majorana_UFO.list -- ######################################################################## # NB: the *overall* sign of an apmplitude may change depending on whether # we compute with Fusion.Mixed23 or Fusion.Mixed23_Majorana ######################################################################## # -# ---------------------------------------------------------------------- -# thr abs_thr n roots model sign process ... -# ---------------------------------------------------------------------- -eeee 0.90 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- -#eeea 0.40 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- a -#eeeaa 0.90 1E-11 1000 1000 SM .false. scatter e- e- -> e- e- a a -#eeeep 0.80 1E-11 1000 1000 SM .true. scatter e- e- -> e- e- e+ e- -#eeeeep 0.80 1E-11 1000 1000 SM .false. scatter e- e- -> e- e- e+ e- -aaep 0.90 1E-11 1000 1000 SM .false. scatter a a -> e+ e- -#aepep 0.90 1E-11 1000 1000 SM .true. scatter a a -> e+ e- e+ e- -eaea 0.90 1E-11 1000 1000 SM .false. scatter e- a -> e- a -epaa 0.90 1E-11 1000 1000 SM .true. scatter e+ e- -> a a -epep 0.65 1E-11 1000 1000 SM .false. scatter e+ e- -> e+ e- -#pepa 0.90 1E-11 1000 1000 SM .false. scatter e+ e- -> e+ e- a -#pepaa 0.90 1E-11 1000 1000 SM .false. scatter e+ e- -> e+ e- a a -#pepep 0.80 1E-11 1000 1000 SM .false. scatter e+ e- -> e+ e- e+ e- -epvv 0.90 1E-11 1000 1000 SM .false. scatter e+ e- -> ve~ ve -epvva 0.90 1E-11 1000 1000 SM .false. scatter e+ e- -> ve~ ve a -vvep 0.90 1E-11 1000 1000 SM .false. scatter ve~ ve -> e+ e- -#vepa 0.40 1E-11 1000 1000 SM .false. scatter ve~ ve -> e+ e- a +# --------------------------------------------------------------------- +# thr abs_thr n roots model sign process ... +# --------------------------------------------------------------------- +eeee 0.90 1E-11 100 1000 SM .true. scatter e- e- -> e- e- +#eeea 0.40 1E-11 100 1000 SM .true. scatter e- e- -> e- e- a +#eeeaa 0.90 1E-11 100 1000 SM .false. scatter e- e- -> e- e- a a +#eeeep 0.80 1E-11 100 1000 SM .true. scatter e- e- -> e- e- e+ e- +#eeeeep 0.80 1E-11 100 1000 SM .false. scatter e- e- -> e- e- e+ e- +aaep 0.90 1E-11 100 1000 SM .false. scatter a a -> e+ e- +#aepep 0.90 1E-11 100 1000 SM .true. scatter a a -> e+ e- e+ e- +eaea 0.90 1E-11 100 1000 SM .false. scatter e- a -> e- a +epaa 0.90 1E-11 100 1000 SM .true. scatter e+ e- -> a a +epep 0.65 1E-11 100 1000 SM .false. scatter e+ e- -> e+ e- +#pepa 0.90 1E-11 100 1000 SM .false. scatter e+ e- -> e+ e- a +#pepaa 0.90 1E-11 100 1000 SM .false. scatter e+ e- -> e+ e- a a +#pepep 0.80 1E-11 100 1000 SM .false. scatter e+ e- -> e+ e- e+ e- +epvv 0.90 1E-11 100 1000 SM .false. scatter e+ e- -> ve~ ve +epvva 0.90 1E-11 100 1000 SM .false. scatter e+ e- -> ve~ ve a +vvep 0.90 1E-11 100 1000 SM .false. scatter ve~ ve -> e+ e- +#vepa 0.40 1E-11 100 1000 SM .false. scatter ve~ ve -> e+ e- a