Index: trunk/share/tests/functional_tests/ref-output/model_test.ref =================================================================== --- trunk/share/tests/functional_tests/ref-output/model_test.ref (revision 8479) +++ trunk/share/tests/functional_tests/ref-output/model_test.ref (revision 8480) @@ -1,1810 +1,1810 @@ ?openmp_logging = false ?vis_history = false ?integration_timer = false seed = 0 | Switching to model 'THDM' | Process library 'model_test_1_lib': recorded process 'models_ccww_thdm' | Switching to model 'THDM_CKM' | Process library 'model_test_1_lib': recorded process 'models_ccww_thdm_ckm' | Switching to model 'GravTest' GravTest.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_grav_test' | Switching to model 'HSExt' HSExt.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_hsext' | Switching to model 'Littlest' Littlest.st => 1.000000000000E-01 Littlest.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_littlest' | Switching to model 'Littlest_Eta' Littlest_Eta.st => 5.000000000000E-01 Littlest_Eta.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_littlest_eta' | Switching to model 'Littlest_Tpar' Littlest_Tpar.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_littlest_tpar' | Switching to model 'MSSM' | Reading SLHA input file 'sps1ap_decays.slha' | SLHA: Initializing model 'MSSM' | SLHA: SUSY spectrum program info: | SLHA: SPheno | SLHA: 2.2.3 | Process library 'model_test_1_lib': recorded process 'models_ccww_mssm' | Switching to model 'MSSM_CKM' | Reading SLHA input file 'sps1ap_decays.slha' | SLHA: Initializing model 'MSSM_CKM' | SLHA: SUSY spectrum program info: | SLHA: SPheno | SLHA: 2.2.3 | Process library 'model_test_1_lib': recorded process 'models_ccww_mssm_ckm' | Switching to model 'MSSM_Grav' | Reading SLHA input file 'sps1ap_decays.slha' | SLHA: Initializing model 'MSSM_Grav' | SLHA: SUSY spectrum program info: | SLHA: SPheno | SLHA: 2.2.3 | Process library 'model_test_1_lib': recorded process 'models_ccww_mssm_grav' | Switching to model 'MSSM_Hgg' | Reading SLHA input file 'sps1ap_decays.slha' | SLHA: Initializing model 'MSSM_Hgg' | SLHA: SUSY spectrum program info: | SLHA: SPheno | SLHA: 2.2.3 | Process library 'model_test_1_lib': recorded process 'models_ccww_mssm_hgg' | Switching to model 'NMSSM' | Reading SLHA input file 'nmssm.slha' | SLHA: Initializing model 'NMSSM' | SLHA: SUSY spectrum program info: | SLHA: NMSSMTools | SLHA: 2 | Process library 'model_test_1_lib': recorded process 'models_ccww_nmssm' | Switching to model 'NMSSM_CKM' | Reading SLHA input file 'nmssm.slha' | SLHA: Initializing model 'NMSSM_CKM' | SLHA: SUSY spectrum program info: | SLHA: NMSSMTools | SLHA: 2 | Process library 'model_test_1_lib': recorded process 'models_ccww_nmssm_ckm' | Switching to model 'NMSSM_Hgg' | Reading SLHA input file 'nmssm.slha' | SLHA: Initializing model 'NMSSM_Hgg' | SLHA: SUSY spectrum program info: | SLHA: NMSSMTools | SLHA: 2 | Process library 'model_test_1_lib': recorded process 'models_ccww_nmssm_hgg' | Switching to model 'PSSSM' | Process library 'model_test_1_lib': recorded process 'models_ccww_psssm' | Switching to model 'QCD', scheme 'default' | Process library 'model_test_1_lib': recorded process 'models_uudd_qcd' | Switching to model 'QED' | Process library 'model_test_1_lib': recorded process 'models_eemm_qed' | Switching to model 'SM', scheme 'default' SM.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm' | Switching to model 'SM_CKM' SM_CKM.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_ckm' | Switching to model 'SM_ac' SM_ac.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_ac' | Switching to model 'SM_ac_CKM' SM_ac_CKM.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_ac_ckm' | Switching to model 'SM_dim6' SM_dim6.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_dim6' | Switching to model 'SM_rx' SM_rx.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_wwww_sm_rx' | Switching to model 'SM_ul' SM_ul.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_wwww_sm_ul' | Switching to model 'NoH_rx' | Process library 'model_test_1_lib': recorded process 'models_wwww_noh_rx' | Switching to model 'SM_top' SM_top.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_top' | Switching to model 'SM_top_anom' SM_top_anom.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_top_anom' | Switching to model 'SM_Higgs' SM_Higgs.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_higgs' | Switching to model 'SM_Higgs_CKM' SM_Higgs_CKM.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_higgs_ckm' | Switching to model 'SM_tt_threshold' SM_tt_threshold.mH => 1.250000000000E+02 SM_tt_threshold.test => 4.200000000000E+01 | Process library 'model_test_1_lib': recorded process 'models_ccww_sm_tt_thresh' | Switching to model 'Simplest' Simplest.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_simplest' | Switching to model 'Simplest_univ' Simplest_univ.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_simplest_univ' | Switching to model 'Template' Template.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_template' | Switching to model 'Threeshl' | Process library 'model_test_1_lib': recorded process 'models_ccww_threeshl' | Switching to model 'Threeshl_nohf' | Process library 'model_test_1_lib': recorded process 'models_ccww_threeshl_nohf' | Switching to model 'UED' UED.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_ued' | Switching to model 'SSC' SSC.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_wwww_ssc' | Switching to model 'SSC_2' SSC_2.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_wwww_ssc_2' | Switching to model 'SSC_AltT' SSC_AltT.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_wwww_ssc_altt' | Switching to model 'AltH' | Process library 'model_test_1_lib': recorded process 'models_wwww_alth' | Switching to model 'Xdim' Xdim.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_xdim' | Switching to model 'WZW' WZW.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_wzw' | Switching to model 'Zprime' Zprime.mH => 1.250000000000E+02 | Process library 'model_test_1_lib': recorded process 'models_ccww_zprime' | Process library 'model_test_lib': compiling ... | Process library 'model_test_lib': ... success. | Process library 'model_test_1_lib': compiling ... | Process library 'model_test_1_lib': writing makefile | Process library 'model_test_1_lib': removing old files | Process library 'model_test_1_lib': writing driver | Process library 'model_test_1_lib': creating source code | Process library 'model_test_1_lib': compiling sources | Process library 'model_test_1_lib': linking | Process library 'model_test_1_lib': loading | Process library 'model_test_1_lib': ... success. sqrts = 2.000000000000E+02 openmp_num_threads = 1 seed = 0 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 0 | Initializing integration for process models_ccww_thdm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_thdm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_thdm' | Library name = 'model_test_1_lib' | Process index = 1 | Process components: | 1: 'models_ccww_thdm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_thdm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 7.5144851E+03 6.76E+02 8.99 0.89* 40.19 |-----------------------------------------------------------------------------| 1 99 7.5144851E+03 6.76E+02 8.99 0.89 40.19 |=============================================================================| seed = 1 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 1 | Initializing integration for process models_ccww_thdm_ckm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_thdm_ckm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_thdm_ckm' | Library name = 'model_test_1_lib' | Process index = 2 | Process components: | 1: 'models_ccww_thdm_ckm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_thdm_ckm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 7.3770184E+03 6.73E+02 9.12 0.91* 39.73 |-----------------------------------------------------------------------------| 1 99 7.3770184E+03 6.73E+02 9.12 0.91 39.73 |=============================================================================| seed = 2 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 2 | Initializing integration for process models_ccww_grav_test: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_grav_test.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_grav_test' | Library name = 'model_test_1_lib' | Process index = 3 | Process components: | 1: 'models_ccww_grav_test_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_grav_test' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.5673184E+03 5.70E+02 10.24 1.02* 29.75 |-----------------------------------------------------------------------------| 1 99 5.5673184E+03 5.70E+02 10.24 1.02 29.75 |=============================================================================| seed = 3 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 3 | Initializing integration for process models_ccww_hsext: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_hsext.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_hsext' | Library name = 'model_test_1_lib' | Process index = 4 | Process components: | 1: 'models_ccww_hsext_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_hsext' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 8.3080339E+05 5.36E+04 6.46 0.64* 47.18 |-----------------------------------------------------------------------------| 1 99 8.3080339E+05 5.36E+04 6.46 0.64 47.18 |=============================================================================| seed = 4 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 4 | Initializing integration for process models_ccww_littlest: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_littlest.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_littlest' | Library name = 'model_test_1_lib' | Process index = 5 | Process components: | 1: 'models_ccww_littlest_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_littlest' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.7081379E+03 5.93E+02 8.84 0.88* 35.88 |-----------------------------------------------------------------------------| 1 99 6.7081379E+03 5.93E+02 8.84 0.88 35.88 |=============================================================================| seed = 5 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 5 | Initializing integration for process models_ccww_littlest_eta: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_littlest_eta.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_littlest_eta' | Library name = 'model_test_1_lib' | Process index = 6 | Process components: | 1: 'models_ccww_littlest_eta_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_littlest_eta' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.9060694E+03 5.73E+02 9.70 0.96* 31.86 |-----------------------------------------------------------------------------| 1 99 5.9060694E+03 5.73E+02 9.70 0.96 31.86 |=============================================================================| seed = 6 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 6 | Initializing integration for process models_ccww_littlest_tpar: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_littlest_tpar.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_littlest_tpar' | Library name = 'model_test_1_lib' | Process index = 7 | Process components: | 1: 'models_ccww_littlest_tpar_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_littlest_tpar' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.2761051E+03 5.73E+02 9.13 0.91* 35.11 |-----------------------------------------------------------------------------| 1 99 6.2761051E+03 5.73E+02 9.13 0.91 35.11 |=============================================================================| seed = 7 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 7 | Initializing integration for process models_ccww_mssm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_mssm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_mssm' | Library name = 'model_test_1_lib' | Process index = 8 | Process components: | 1: 'models_ccww_mssm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_mssm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.9494213E+03 6.28E+02 10.55 1.05* 31.92 |-----------------------------------------------------------------------------| 1 99 5.9494213E+03 6.28E+02 10.55 1.05 31.92 |=============================================================================| seed = 8 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 8 | Initializing integration for process models_ccww_mssm_ckm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_mssm_ckm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_mssm_ckm' | Library name = 'model_test_1_lib' | Process index = 9 | Process components: | 1: 'models_ccww_mssm_ckm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_mssm_ckm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| - 1 99 6.9673584E+03 6.29E+02 9.03 0.90* 37.22 + 1 99 6.9673757E+03 6.29E+02 9.03 0.90* 37.22 |-----------------------------------------------------------------------------| - 1 99 6.9673584E+03 6.29E+02 9.03 0.90 37.22 + 1 99 6.9673757E+03 6.29E+02 9.03 0.90 37.22 |=============================================================================| seed = 9 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 9 | Initializing integration for process models_ccww_mssm_grav: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_mssm_grav.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_mssm_grav' | Library name = 'model_test_1_lib' | Process index = 10 | Process components: | 1: 'models_ccww_mssm_grav_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_mssm_grav' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 7.0900237E+03 6.59E+02 9.29 0.92* 38.09 |-----------------------------------------------------------------------------| 1 99 7.0900237E+03 6.59E+02 9.29 0.92 38.09 |=============================================================================| seed = 10 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 10 | Initializing integration for process models_ccww_nmssm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_nmssm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_nmssm' | Library name = 'model_test_1_lib' | Process index = 12 | Process components: | 1: 'models_ccww_nmssm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_nmssm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.3772942E+03 5.41E+02 10.06 1.00* 29.82 |-----------------------------------------------------------------------------| 1 99 5.3772942E+03 5.41E+02 10.06 1.00 29.82 |=============================================================================| seed = 11 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 11 | Initializing integration for process models_ccww_nmssm_ckm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_nmssm_ckm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_nmssm_ckm' | Library name = 'model_test_1_lib' | Process index = 13 | Process components: | 1: 'models_ccww_nmssm_ckm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_nmssm_ckm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| - 1 99 5.6207771E+03 5.42E+02 9.64 0.96* 31.18 + 1 99 5.6207898E+03 5.42E+02 9.64 0.96* 31.18 |-----------------------------------------------------------------------------| - 1 99 5.6207771E+03 5.42E+02 9.64 0.96 31.18 + 1 99 5.6207898E+03 5.42E+02 9.64 0.96 31.18 |=============================================================================| seed = 12 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 12 | Initializing integration for process models_ccww_nmssm_hgg: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_nmssm_hgg.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_nmssm_hgg' | Library name = 'model_test_1_lib' | Process index = 14 | Process components: | 1: 'models_ccww_nmssm_hgg_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_nmssm_hgg' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.8020707E+03 5.14E+02 8.86 0.88* 32.08 |-----------------------------------------------------------------------------| 1 99 5.8020707E+03 5.14E+02 8.86 0.88 32.08 |=============================================================================| seed = 13 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 13 | Initializing integration for process models_ccww_psssm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_psssm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_psssm' | Library name = 'model_test_1_lib' | Process index = 15 | Process components: | 1: 'models_ccww_psssm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_psssm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.4182194E+08 5.96E+07 9.28 0.92* 34.08 |-----------------------------------------------------------------------------| 1 99 6.4182194E+08 5.96E+07 9.28 0.92 34.08 |=============================================================================| seed = 14 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 14 | Initializing integration for process models_uudd_qcd: | Beam structure: [any particles] | Beam data (collision): | u (mass = 0.0000000E+00 GeV) | ubar (mass = 0.0000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_uudd_qcd.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_uudd_qcd' | Library name = 'model_test_1_lib' | Process index = 16 | Process components: | 1: 'models_uudd_qcd_i1': u, ubar => d, dbar [omega] | ------------------------------------------------------------------------ | Phase space: 1 channels, 2 dimensions | Phase space: found 1 channel, collected in 1 grove. | Phase space: Using 1 equivalence between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_uudd_qcd' | Integrate: iterations = 1:100 | Integrator: 1 chains, 1 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 100 1.2587507E+05 2.94E+03 2.34 0.23* 67.96 |-----------------------------------------------------------------------------| 1 100 1.2587507E+05 2.94E+03 2.34 0.23 67.96 |=============================================================================| seed = 15 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 15 | Initializing integration for process models_eemm_qed: | Beam structure: [any particles] | Beam data (collision): | e- (mass = 5.1100000E-04 GeV) | e+ (mass = 5.1100000E-04 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_eemm_qed.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_eemm_qed' | Library name = 'model_test_1_lib' | Process index = 17 | Process components: | 1: 'models_eemm_qed_i1': e-, e+ => m-, m+ [omega] | ------------------------------------------------------------------------ | Phase space: 1 channels, 2 dimensions | Phase space: found 1 channel, collected in 1 grove. | Phase space: Using 1 equivalence between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_eemm_qed' | Integrate: iterations = 1:100 | Integrator: 1 chains, 1 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 100 2.1772937E+03 4.57E+01 2.10 0.21* 66.89 |-----------------------------------------------------------------------------| 1 100 2.1772937E+03 4.57E+01 2.10 0.21 66.89 |=============================================================================| seed = 16 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 16 | Initializing integration for process models_ccww_sm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm' | Library name = 'model_test_1_lib' | Process index = 18 | Process components: | 1: 'models_ccww_sm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.4136006E+03 5.84E+02 10.79 1.07* 29.81 |-----------------------------------------------------------------------------| 1 99 5.4136006E+03 5.84E+02 10.79 1.07 29.81 |=============================================================================| seed = 17 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 17 | Initializing integration for process models_ccww_sm_ckm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_ckm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_ckm' | Library name = 'model_test_1_lib' | Process index = 19 | Process components: | 1: 'models_ccww_sm_ckm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_ckm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| - 1 99 5.8798557E+03 5.74E+02 9.76 0.97* 33.05 + 1 99 5.8797654E+03 5.74E+02 9.76 0.97* 33.05 |-----------------------------------------------------------------------------| - 1 99 5.8798557E+03 5.74E+02 9.76 0.97 33.05 + 1 99 5.8797654E+03 5.74E+02 9.76 0.97 33.05 |=============================================================================| seed = 18 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 18 | Initializing integration for process models_ccww_sm_ac: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_ac.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_ac' | Library name = 'model_test_1_lib' | Process index = 20 | Process components: | 1: 'models_ccww_sm_ac_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_ac' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.7751152E+03 5.94E+02 10.28 1.02* 31.38 |-----------------------------------------------------------------------------| 1 99 5.7751152E+03 5.94E+02 10.28 1.02 31.38 |=============================================================================| seed = 19 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 19 | Initializing integration for process models_ccww_sm_ac_ckm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_ac_ckm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_ac_ckm' | Library name = 'model_test_1_lib' | Process index = 21 | Process components: | 1: 'models_ccww_sm_ac_ckm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_ac_ckm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| - 1 99 5.6718111E+03 6.23E+02 10.98 1.09* 31.58 + 1 99 5.6717247E+03 6.23E+02 10.98 1.09* 31.58 |-----------------------------------------------------------------------------| - 1 99 5.6718111E+03 6.23E+02 10.98 1.09 31.58 + 1 99 5.6717247E+03 6.23E+02 10.98 1.09 31.58 |=============================================================================| seed = 20 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 20 | Initializing integration for process models_ccww_sm_dim6: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_dim6.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_dim6' | Library name = 'model_test_1_lib' | Process index = 22 | Process components: | 1: 'models_ccww_sm_dim6_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_dim6' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.3233409E+03 5.61E+02 8.87 0.88* 34.07 |-----------------------------------------------------------------------------| 1 99 6.3233409E+03 5.61E+02 8.87 0.88 34.07 |=============================================================================| seed = 21 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 21 | Initializing integration for process models_wwww_sm_rx: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_sm_rx.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_sm_rx' | Library name = 'model_test_1_lib' | Process index = 23 | Process components: | 1: 'models_wwww_sm_rx_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 4 channels, 2 dimensions | Phase space: found 4 channels, collected in 2 groves. | Phase space: Using 4 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_sm_rx' | Integrate: iterations = 1:100 | Integrator: 2 chains, 4 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 100 3.0227479E+07 9.65E+06 31.94 3.19* 10.86 |-----------------------------------------------------------------------------| 1 100 3.0227479E+07 9.65E+06 31.94 3.19 10.86 |=============================================================================| seed = 22 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 22 | Initializing integration for process models_wwww_sm_ul: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_sm_ul.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_sm_ul' | Library name = 'model_test_1_lib' | Process index = 24 | Process components: | 1: 'models_wwww_sm_ul_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 4 channels, 2 dimensions | Phase space: found 4 channels, collected in 2 groves. | Phase space: Using 4 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_sm_ul' | Integrate: iterations = 1:100 | Integrator: 2 chains, 4 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 100 9.7080025E+08 7.49E+08 77.19 7.72* 4.95 |-----------------------------------------------------------------------------| 1 100 9.7080025E+08 7.49E+08 77.19 7.72 4.95 |=============================================================================| seed = 23 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 23 | Initializing integration for process models_wwww_noh_rx: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_noh_rx.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_noh_rx' | Library name = 'model_test_1_lib' | Process index = 25 | Process components: | 1: 'models_wwww_noh_rx_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 2 channels, 2 dimensions | Phase space: found 2 channels, collected in 1 grove. | Phase space: Using 2 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_noh_rx' | Integrate: iterations = 1:100 | Integrator: 1 chains, 2 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 100 2.1032248E+10 1.99E+10 94.76 9.48* 2.06 |-----------------------------------------------------------------------------| 1 100 2.1032248E+10 1.99E+10 94.76 9.48 2.06 |=============================================================================| seed = 24 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 24 | Initializing integration for process models_ccww_sm_top: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_top.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_top' | Library name = 'model_test_1_lib' | Process index = 26 | Process components: | 1: 'models_ccww_sm_top_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_top' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.7123251E+03 6.25E+02 9.30 0.93* 37.02 |-----------------------------------------------------------------------------| 1 99 6.7123251E+03 6.25E+02 9.30 0.93 37.02 |=============================================================================| seed = 25 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 25 | Initializing integration for process models_ccww_sm_top_anom: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_top_anom.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_top_anom' | Library name = 'model_test_1_lib' | Process index = 27 | Process components: | 1: 'models_ccww_sm_top_anom_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_top_anom' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.0274506E+03 5.41E+02 10.76 1.07* 28.07 |-----------------------------------------------------------------------------| 1 99 5.0274506E+03 5.41E+02 10.76 1.07 28.07 |=============================================================================| seed = 26 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 26 | Initializing integration for process models_ccww_sm_higgs: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_higgs.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_higgs' | Library name = 'model_test_1_lib' | Process index = 28 | Process components: | 1: 'models_ccww_sm_higgs_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_higgs' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.1901689E+03 5.99E+02 9.68 0.96* 33.39 |-----------------------------------------------------------------------------| 1 99 6.1901689E+03 5.99E+02 9.68 0.96 33.39 |=============================================================================| seed = 27 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 27 | Initializing integration for process models_ccww_sm_higgs_ckm: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_higgs_ckm.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_higgs_ckm' | Library name = 'model_test_1_lib' | Process index = 29 | Process components: | 1: 'models_ccww_sm_higgs_ckm_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_higgs_ckm' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| - 1 99 7.2866110E+03 6.44E+02 8.84 0.88* 39.17 + 1 99 7.2865033E+03 6.44E+02 8.84 0.88* 39.17 |-----------------------------------------------------------------------------| - 1 99 7.2866110E+03 6.44E+02 8.84 0.88 39.17 + 1 99 7.2865033E+03 6.44E+02 8.84 0.88 39.17 |=============================================================================| seed = 28 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 28 | Initializing integration for process models_ccww_sm_tt_thresh: | TESTING ONLY: Skip threshold initialization and use tree-level SM. | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_sm_tt_thresh.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_sm_tt_thresh' | Library name = 'model_test_1_lib' | Process index = 30 | Process components: | 1: 'models_ccww_sm_tt_thresh_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_sm_tt_thresh' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 7.1165240E+03 6.76E+02 9.51 0.95* 35.25 |-----------------------------------------------------------------------------| 1 99 7.1165240E+03 6.76E+02 9.51 0.95 35.25 |=============================================================================| seed = 29 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 29 | Initializing integration for process models_ccww_simplest: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_simplest.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_simplest' | Library name = 'model_test_1_lib' | Process index = 31 | Process components: | 1: 'models_ccww_simplest_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_simplest' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 7.7043505E+03 6.46E+02 8.38 0.83* 41.45 |-----------------------------------------------------------------------------| 1 99 7.7043505E+03 6.46E+02 8.38 0.83 41.45 |=============================================================================| seed = 30 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 30 | Initializing integration for process models_ccww_simplest_univ: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_simplest_univ.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_simplest_univ' | Library name = 'model_test_1_lib' | Process index = 32 | Process components: | 1: 'models_ccww_simplest_univ_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_simplest_univ' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.4988268E+03 6.33E+02 9.74 0.97* 34.99 |-----------------------------------------------------------------------------| 1 99 6.4988268E+03 6.33E+02 9.74 0.97 34.99 |=============================================================================| seed = 31 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 31 | Initializing integration for process models_ccww_template: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_template.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_template' | Library name = 'model_test_1_lib' | Process index = 33 | Process components: | 1: 'models_ccww_template_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_template' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.2092827E+03 5.58E+02 10.72 1.07* 28.03 |-----------------------------------------------------------------------------| 1 99 5.2092827E+03 5.58E+02 10.72 1.07 28.03 |=============================================================================| seed = 32 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 32 | Initializing integration for process models_ccww_threeshl: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_threeshl.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_threeshl' | Library name = 'model_test_1_lib' | Process index = 34 | Process components: | 1: 'models_ccww_threeshl_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 5 channels, 2 dimensions | Phase space: found 5 channels, collected in 3 groves. | Phase space: Using 5 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_threeshl' | Integrate: iterations = 1:100 | Integrator: 3 chains, 5 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 100 8.1746005E+03 6.73E+02 8.23 0.82* 42.96 |-----------------------------------------------------------------------------| 1 100 8.1746005E+03 6.73E+02 8.23 0.82 42.96 |=============================================================================| seed = 33 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 33 | Initializing integration for process models_ccww_threeshl_nohf: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_threeshl_nohf.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_threeshl_nohf' | Library name = 'model_test_1_lib' | Process index = 35 | Process components: | 1: 'models_ccww_threeshl_nohf_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_threeshl_nohf' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.5389875E+03 6.31E+02 9.65 0.96* 32.66 |-----------------------------------------------------------------------------| 1 99 6.5389875E+03 6.31E+02 9.65 0.96 32.66 |=============================================================================| seed = 34 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 34 | Initializing integration for process models_ccww_ued: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2500000E+00 GeV) | cbar (mass = 1.2500000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_ued.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_ued' | Library name = 'model_test_1_lib' | Process index = 36 | Process components: | 1: 'models_ccww_ued_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_ued' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.9392959E+08 3.86E+07 6.50 0.65* 41.63 |-----------------------------------------------------------------------------| 1 99 5.9392959E+08 3.86E+07 6.50 0.65 41.63 |=============================================================================| seed = 35 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 35 | Initializing integration for process models_wwww_ssc: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_ssc.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_ssc' | Library name = 'model_test_1_lib' | Process index = 37 | Process components: | 1: 'models_wwww_ssc_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 6 channels, 2 dimensions | Phase space: found 6 channels, collected in 2 groves. | Phase space: Using 10 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_ssc' | Integrate: iterations = 1:100 | Integrator: 2 chains, 6 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 96 1.7800933E+08 1.31E+08 73.83 7.23* 7.03 |-----------------------------------------------------------------------------| 1 96 1.7800933E+08 1.31E+08 73.83 7.23 7.03 |=============================================================================| seed = 36 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 36 | Initializing integration for process models_wwww_ssc_2: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_ssc_2.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_ssc_2' | Library name = 'model_test_1_lib' | Process index = 38 | Process components: | 1: 'models_wwww_ssc_2_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 6 channels, 2 dimensions | Phase space: found 6 channels, collected in 2 groves. | Phase space: Using 10 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_ssc_2' | Integrate: iterations = 1:100 | Integrator: 2 chains, 6 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 96 8.7439298E+06 4.17E+06 47.71 4.68* 10.61 |-----------------------------------------------------------------------------| 1 96 8.7439298E+06 4.17E+06 47.71 4.68 10.61 |=============================================================================| seed = 37 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 37 | Initializing integration for process models_wwww_ssc_altt: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_ssc_altt.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_ssc_altt' | Library name = 'model_test_1_lib' | Process index = 39 | Process components: | 1: 'models_wwww_ssc_altt_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 6 channels, 2 dimensions | Phase space: found 6 channels, collected in 2 groves. | Phase space: Using 10 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_ssc_altt' | Integrate: iterations = 1:100 | Integrator: 2 chains, 6 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 96 9.7035892E+08 6.17E+08 63.60 6.23* 6.93 |-----------------------------------------------------------------------------| 1 96 9.7035892E+08 6.17E+08 63.60 6.23 6.93 |=============================================================================| seed = 38 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 38 | Initializing integration for process models_wwww_alth: | Beam structure: [any particles] | Beam data (collision): | W+ (mass = 8.0419000E+01 GeV) | W- (mass = 8.0419000E+01 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_wwww_alth.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_wwww_alth' | Library name = 'model_test_1_lib' | Process index = 40 | Process components: | 1: 'models_wwww_alth_i1': W+, W- => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 6 channels, 2 dimensions | Phase space: found 6 channels, collected in 2 groves. | Phase space: Using 10 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_wwww_alth' | Integrate: iterations = 1:100 | Integrator: 2 chains, 6 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 96 2.3538217E+08 2.26E+08 96.08 9.41* 6.39 |-----------------------------------------------------------------------------| 1 96 2.3538217E+08 2.26E+08 96.08 9.41 6.39 |=============================================================================| seed = 39 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 39 | Initializing integration for process models_ccww_xdim: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_xdim.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_xdim' | Library name = 'model_test_1_lib' | Process index = 41 | Process components: | 1: 'models_ccww_xdim_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_xdim' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 6.5751622E+08 4.44E+07 6.75 0.67* 43.94 |-----------------------------------------------------------------------------| 1 99 6.5751622E+08 4.44E+07 6.75 0.67 43.94 |=============================================================================| seed = 40 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 40 | Initializing integration for process models_ccww_wzw: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_wzw.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_wzw' | Library name = 'model_test_1_lib' | Process index = 42 | Process components: | 1: 'models_ccww_wzw_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_wzw' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.9222102E+03 6.23E+02 10.52 1.05* 31.63 |-----------------------------------------------------------------------------| 1 99 5.9222102E+03 6.23E+02 10.52 1.05 31.63 |=============================================================================| seed = 41 | RNG: Initializing TAO random-number generator | RNG: Setting seed for random-number generator to 41 | Initializing integration for process models_ccww_zprime: | Beam structure: [any particles] | Beam data (collision): | c (mass = 1.2000000E+00 GeV) | cbar (mass = 1.2000000E+00 GeV) | sqrts = 2.000000000000E+02 GeV | Phase space: generating configuration ... | Phase space: ... success. | Phase space: writing configuration file 'models_ccww_zprime.i1.phs' | ------------------------------------------------------------------------ | Process [scattering]: 'models_ccww_zprime' | Library name = 'model_test_1_lib' | Process index = 43 | Process components: | 1: 'models_ccww_zprime_i1': c, cbar => W+, W- [omega] | ------------------------------------------------------------------------ | Phase space: 3 channels, 2 dimensions | Phase space: found 3 channels, collected in 2 groves. | Phase space: Using 3 equivalences between channels. | Phase space: wood Warning: No cuts have been defined. | Starting integration for process 'models_ccww_zprime' | Integrate: iterations = 1:100 | Integrator: 2 chains, 3 channels, 2 dimensions | Integrator: Using VAMP channel equivalences | Integrator: 100 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 99 5.6983638E+03 6.02E+02 10.57 1.05* 30.61 |-----------------------------------------------------------------------------| 1 99 5.6983638E+03 6.02E+02 10.57 1.05 30.61 |=============================================================================| models_ccww_thdm: 7.5144851E+03 +- 6.76E+02 fb ( 9.0 %) models_ccww_thdm_ckm: 7.3770184E+03 +- 6.73E+02 fb ( 9.1 %) models_ccww_grav_test: 5.5673184E+03 +- 5.70E+02 fb ( 10. %) models_ccww_hsext: 8.3080339E+05 +- 5.36E+04 fb ( 6.5 %) models_ccww_littlest: 6.7081379E+03 +- 5.93E+02 fb ( 8.8 %) models_ccww_littlest_eta: 5.9060694E+03 +- 5.73E+02 fb ( 9.7 %) models_ccww_littlest_tpar: 6.2761051E+03 +- 5.73E+02 fb ( 9.1 %) models_ccww_mssm: 5.9494213E+03 +- 6.28E+02 fb ( 11. %) models_ccww_mssm_ckm: - 6.9673584E+03 +- 6.29E+02 fb ( 9.0 %) + 6.9673757E+03 +- 6.29E+02 fb ( 9.0 %) models_ccww_mssm_grav: 7.0900237E+03 +- 6.59E+02 fb ( 9.3 %) models_ccww_nmssm: 5.3772942E+03 +- 5.41E+02 fb ( 10. %) models_ccww_nmssm_ckm: - 5.6207771E+03 +- 5.42E+02 fb ( 9.6 %) + 5.6207898E+03 +- 5.42E+02 fb ( 9.6 %) models_ccww_nmssm_hgg: 5.8020707E+03 +- 5.14E+02 fb ( 8.9 %) models_ccww_psssm: 6.4182194E+08 +- 5.96E+07 fb ( 9.3 %) models_uudd_qcd: 1.2587507E+05 +- 2.94E+03 fb ( 2.3 %) models_eemm_qed: 2.1772937E+03 +- 4.57E+01 fb ( 2.1 %) models_ccww_sm: 5.4136006E+03 +- 5.84E+02 fb ( 11. %) models_ccww_sm_ckm: - 5.8798557E+03 +- 5.74E+02 fb ( 9.8 %) + 5.8797654E+03 +- 5.74E+02 fb ( 9.8 %) models_ccww_sm_ac: 5.7751152E+03 +- 5.94E+02 fb ( 10. %) models_ccww_sm_ac_ckm: - 5.6718111E+03 +- 6.23E+02 fb ( 11. %) + 5.6717247E+03 +- 6.23E+02 fb ( 11. %) models_ccww_sm_dim6: 6.3233409E+03 +- 5.61E+02 fb ( 8.9 %) models_wwww_sm_rx: 3.0227479E+07 +- 9.65E+06 fb ( 32. %) models_wwww_sm_ul: 9.7080025E+08 +- 7.49E+08 fb ( 77. %) models_wwww_noh_rx: 2.1032248E+10 +- 1.99E+10 fb ( 95. %) models_ccww_sm_top: 6.7123251E+03 +- 6.25E+02 fb ( 9.3 %) models_ccww_sm_top_anom: 5.0274506E+03 +- 5.41E+02 fb ( 11. %) models_ccww_sm_higgs: 6.1901689E+03 +- 5.99E+02 fb ( 9.7 %) models_ccww_sm_higgs_ckm: - 7.2866110E+03 +- 6.44E+02 fb ( 8.8 %) + 7.2865033E+03 +- 6.44E+02 fb ( 8.8 %) models_ccww_sm_tt_thresh: 7.1165240E+03 +- 6.76E+02 fb ( 9.5 %) models_ccww_simplest: 7.7043505E+03 +- 6.46E+02 fb ( 8.4 %) models_ccww_simplest_univ: 6.4988268E+03 +- 6.33E+02 fb ( 9.7 %) models_ccww_template: 5.2092827E+03 +- 5.58E+02 fb ( 11. %) models_ccww_threeshl: 8.1746005E+03 +- 6.73E+02 fb ( 8.2 %) models_ccww_threeshl_nohf: 6.5389875E+03 +- 6.31E+02 fb ( 9.6 %) models_ccww_ued: 5.9392959E+08 +- 3.86E+07 fb ( 6.5 %) models_wwww_ssc: 1.7800933E+08 +- 1.31E+08 fb ( 74. %) models_wwww_ssc_2: 8.7439298E+06 +- 4.17E+06 fb ( 48. %) models_wwww_ssc_altt: 9.7035892E+08 +- 6.17E+08 fb ( 64. %) models_wwww_alth: 2.3538217E+08 +- 2.26E+08 fb ( 96. %) models_ccww_xdim: 6.5751622E+08 +- 4.44E+07 fb ( 6.7 %) models_ccww_wzw: 5.9222102E+03 +- 6.23E+02 fb ( 11. %) models_ccww_zprime: 5.6983638E+03 +- 6.02E+02 fb ( 11. %) | There were no errors and 42 warning(s). | WHIZARD run finished. |=============================================================================| Index: trunk/share/models/SM_Higgs_CKM.mdl =================================================================== --- trunk/share/models/SM_Higgs_CKM.mdl (revision 8479) +++ trunk/share/models/SM_Higgs_CKM.mdl (revision 8480) @@ -1,273 +1,273 @@ ######################################################################## # Standard Model with full CKM matrix, including the effective # Hgg / HAA / HZA / Hmm / Hee couplings model "SM_Higgs_CKM" # Independent parameters ### DO NOT CHANGE THE ORDER OF THESE PARAMETERS parameter GF = 1.16639E-5 # Fermi constant parameter mZ = 91.1882 # Z-boson mass parameter mW = 80.419 # W-boson mass parameter mH = 125 # Higgs mass parameter alphas = 0.1178 # Strong coupling constant (Z point) parameter me = 0.000510997 # electron mass parameter mmu = 0.105658389 # muon mass parameter mtau = 1.77705 # tau-lepton mass parameter ms = 0.095 # s-quark mass parameter mc = 1.2 # c-quark mass parameter mb = 4.2 # b-quark mass parameter mtop = 173.1 # t-quark mass parameter wtop = 1.523 # t-quark width parameter wZ = 2.443 # Z-boson width parameter wW = 2.049 # W-boson width parameter wH = 0.004143 # Higgs width -# Here are the values from PDG 2006 -parameter vckm11 = 0.97383 # Vud -parameter vckm12 = 0.2272 # Vus -parameter vckm13 = 0.00396 # Vub -parameter vckm21 = -0.2271 # Vcd -parameter vckm22 = 0.97296 # Vcs -parameter vckm23 = 0.04221 # Vcb -parameter vckm31 = 0.00814 # Vtd -parameter vckm32 = -0.04161 # Vts -parameter vckm33 = 0.99910 # Vtb +# Here are the values from PDG 2020 +parameter vckm11 = 0.97401 # Vud +parameter vckm12 = 0.22650 # Vus +parameter vckm13 = 0.00361 # Vub +parameter vckm21 = -0.22636 # Vcd +parameter vckm22 = 0.97320 # Vcs +parameter vckm23 = 0.04053 # Vcb +parameter vckm31 = 0.00854 # Vtd +parameter vckm32 = -0.03978 # Vts +parameter vckm33 = 0.999172 # Vtb parameter khgaz = 1.000 # anomaly Higgs couplings K factors parameter khgaga = 1.000 # anomaly Higgs couplings K factors parameter khgg = 1.000 # anomaly Higgs couplings K factors parameter xi0 = 0.000 # R_xi parameter for Z-boson parameter xipm = 0.000 # R_xi parameter for W-boson # Dependent parameters derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev) derived cw = mW / mZ # cos(theta-W) derived sw = sqrt (1-cw**2) # sin(theta-W) derived ee = 2 * sw * mW / v # em-coupling (GF scheme) derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const ######################################################################## # Particle content # The quarks particle D_QUARK 1 parton spin 1/2 charge -1/3 isospin -1/2 color 3 name d down anti dbar D "d~" tex_anti "\bar{d}" particle U_QUARK 2 parton spin 1/2 charge 2/3 isospin 1/2 color 3 name u up anti ubar U "u~" tex_anti "\bar{u}" particle S_QUARK 3 like D_QUARK name s strange anti sbar S "s~" tex_anti "\bar{s}" mass ms particle C_QUARK 4 like U_QUARK name c charm anti cbar C "c~" tex_anti "\bar{c}" mass mc particle B_QUARK 5 like D_QUARK name b bottom anti bbar B "b~" tex_anti "\bar{b}" mass mb particle T_QUARK 6 like U_QUARK name t top anti tbar T "t~" tex_anti "\bar{t}" mass mtop width wtop # The leptons particle E_LEPTON 11 spin 1/2 charge -1 isospin -1/2 name "e-" e1 electron e anti "e+" E1 positron tex_name "e^-" tex_anti "e^+" mass me particle E_NEUTRINO 12 left spin 1/2 isospin 1/2 name nue n1 "nu_e" ve "e-neutrino" anti nuebar N1 "ve~" tex_name "\nu_e" tex_anti "\bar\nu_e" particle MU_LEPTON 13 like E_LEPTON name "mu-" e2 mu muon anti "mu+" E2 tex_name "\mu^-" tex_anti "\mu^+" mass mmu particle MU_NEUTRINO 14 like E_NEUTRINO name numu "nu_mu" n2 vm "mu-neutrino" anti numubar N2 "vm~" tex_name "\nu_\mu" tex_anti "\bar\nu_\mu" particle TAU_LEPTON 15 like E_LEPTON name "tau-" e3 tau "ta-" tauon anti "tau+" E3 "ta+" tex_name "\tau^-" tex_anti "\tau^+" mass mtau particle TAU_NEUTRINO 16 like E_NEUTRINO name nutau "nu_tau" n3 vt "tau_neutrino" anti nutaubar N3 "vt~" tex_name "\nu_\tau" tex_anti "\bar\nu_\tau" # The vector bosons particle GLUON 21 parton gauge spin 1 color 8 name gl g G gluon particle PHOTON 22 gauge spin 1 name A gamma photon tex_name "\gamma" particle Z_BOSON 23 gauge spin 1 name Z mass mZ width wZ particle W_BOSON 24 gauge spin 1 charge 1 name "W+" Wp anti "W-" Wm tex_name "W^+" tex_anti "W^-" mass mW width wW # The Higgs particle HIGGS 25 spin 0 name H h Higgs mass mH width wH # Hadrons particle PROTON 2212 spin 1/2 charge 1 name p "p+" anti pbar "p-" # Beam remnants for proton colliders particle HADRON_REMNANT 90 name hr tex_name "had_r" particle HADRON_REMNANT_SINGLET 91 name hr1 tex_name "had_r^{(1)}" particle HADRON_REMNANT_TRIPLET 92 color 3 name hr3 tex_name "had_r^{(3)}" anti hr3bar tex_anti "had_r^{(\bar 3)}" particle HADRON_REMNANT_OCTET 93 color 8 name hr8 tex_name "had_r^{(8)}" ######################################################################## # Vertices of the Standard model # In graphs with identical structure the first vertex is kept for phase space # therefore, lighter particles come before heavier ones. # QED vertex D d A vertex U u A vertex S s A vertex C c A vertex B b A vertex T t A vertex E1 e1 A vertex E2 e2 A vertex E3 e3 A # QCD vertex G G G vertex G G G G vertex D d G vertex U u G vertex S s G vertex C c G vertex B b G vertex T t G # Neutral currents vertex D d Z vertex U u Z vertex S s Z vertex C c Z vertex B b Z vertex T t Z vertex E1 e1 Z vertex E2 e2 Z vertex E3 e3 Z vertex N1 n1 Z vertex N2 n2 Z vertex N3 n3 Z # Charged currents vertex U d Wp vertex U s Wp vertex U b Wp vertex C d Wp vertex C s Wp vertex C b Wp vertex T d Wp vertex T s Wp vertex T b Wp vertex D u Wm vertex D c Wm vertex D t Wm vertex S u Wm vertex S c Wm vertex S t Wm vertex B u Wm vertex B c Wm vertex B t Wm vertex N1 e1 Wp vertex N2 e2 Wp vertex N3 e3 Wp vertex E1 n1 Wm vertex E2 n2 Wm vertex E3 n3 Wm # Yukawa vertex S s H vertex C c H vertex B b H vertex T t H vertex E1 e1 H vertex E2 e2 H vertex E3 e3 H # Vector-boson self-interactions vertex Wp Wm A vertex Wp Wm Z vertex Wp Wm Z Z vertex Wp Wp Wm Wm vertex Wp Wm Z A vertex Wp Wm A A # Higgs - vector boson vertex H Z A vertex H A A vertex H g g vertex H Wp Wm vertex H Z Z vertex H H Wp Wm vertex H H Z Z # Higgs self-interactions vertex H H H vertex H H H H Index: trunk/share/models/SM_CKM.mdl =================================================================== --- trunk/share/models/SM_CKM.mdl (revision 8479) +++ trunk/share/models/SM_CKM.mdl (revision 8480) @@ -1,267 +1,267 @@ ######################################################################## # Standard Model with nontrivial CKM matrix model "SM_CKM" # Independent parameters ### DO NOT CHANGE THE ORDER OF THESE PARAMETERS parameter GF = 1.16639E-5 # Fermi constant parameter mZ = 91.1882 # Z-boson mass parameter mW = 80.419 # W-boson mass parameter mH = 125 # Higgs mass parameter alphas = 0.1178 # Strong coupling constant (Z point) parameter me = 0.000510997 # electron mass parameter mmu = 0.105658389 # muon mass parameter mtau = 1.77705 # tau-lepton mass parameter ms = 0.095 # s-quark mass parameter mc = 1.2 # c-quark mass parameter mb = 4.2 # b-quark mass parameter mtop = 170.9 # t-quark mass parameter wtop = 1.523 # t-quark width parameter wZ = 2.443 # Z-boson width parameter wW = 2.049 # W-boson width parameter wH = 0.004143 # Higgs width -# Here are the values from PDG 2006 -parameter vckm11 = 0.97383 # Vud -parameter vckm12 = 0.2272 # Vus -parameter vckm13 = 0.00396 # Vub -parameter vckm21 = -0.2271 # Vcd -parameter vckm22 = 0.97296 # Vcs -parameter vckm23 = 0.04221 # Vcb -parameter vckm31 = 0.00814 # Vtd -parameter vckm32 = -0.04161 # Vts -parameter vckm33 = 0.99910 # Vtb +# Here are the values from PDG 2020 +parameter vckm11 = 0.97401 # Vud +parameter vckm12 = 0.22650 # Vus +parameter vckm13 = 0.00361 # Vub +parameter vckm21 = -0.22636 # Vcd +parameter vckm22 = 0.97320 # Vcs +parameter vckm23 = 0.04053 # Vcb +parameter vckm31 = 0.00854 # Vtd +parameter vckm32 = -0.03978 # Vts +parameter vckm33 = 0.999172 # Vtb parameter khgaz = 0.000 # anomaly Higgs couplings K factors parameter khgaga = 0.000 # anomaly Higgs couplings K factors parameter khgg = 0.000 # anomaly Higgs couplings K factors parameter xi0 = 0.000 # R_xi parameter for Z-boson parameter xipm = 0.000 # R_xi parameter for W-boson # Dependent parameters derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev) derived cw = mW / mZ # cos(theta-W) derived sw = sqrt (1-cw**2) # sin(theta-W) derived ee = 2 * sw * mW / v # em-coupling (GF scheme) derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const ######################################################################## # Particle content # The quarks particle D_QUARK 1 parton spin 1/2 charge -1/3 isospin -1/2 color 3 name d down anti dbar D "d~" tex_anti "\bar{d}" particle U_QUARK 2 parton spin 1/2 charge 2/3 isospin 1/2 color 3 name u up anti ubar U "u~" tex_anti "\bar{u}" particle S_QUARK 3 like D_QUARK name s strange anti sbar S "s~" tex_anti "\bar{s}" mass ms particle C_QUARK 4 like U_QUARK name c charm anti cbar C "c~" tex_anti "\bar{c}" mass mc particle B_QUARK 5 like D_QUARK name b bottom anti bbar B "b~" tex_anti "\bar{b}" mass mb particle T_QUARK 6 like U_QUARK name t top anti tbar T "t~" tex_anti "\bar{t}" mass mtop width wtop # The leptons particle E_LEPTON 11 spin 1/2 charge -1 isospin -1/2 name "e-" e1 electron e anti "e+" E1 positron tex_name "e^-" tex_anti "e^+" mass me particle E_NEUTRINO 12 left spin 1/2 isospin 1/2 name nue n1 "nu_e" ve "e-neutrino" anti nuebar N1 "ve~" tex_name "\nu_e" tex_anti "\bar\nu_e" particle MU_LEPTON 13 like E_LEPTON name "mu-" e2 mu muon anti "mu+" E2 tex_name "\mu^-" tex_anti "\mu^+" mass mmu particle MU_NEUTRINO 14 like E_NEUTRINO name numu "nu_mu" n2 vm "mu-neutrino" anti numubar N2 "vm~" tex_name "\nu_\mu" tex_anti "\bar\nu_\mu" particle TAU_LEPTON 15 like E_LEPTON name "tau-" e3 tau "ta-" tauon anti "tau+" E3 "ta+" tex_name "\tau^-" tex_anti "\tau^+" mass mtau particle TAU_NEUTRINO 16 like E_NEUTRINO name nutau "nu_tau" n3 vt "tau_neutrino" anti nutaubar N3 "vt~" tex_name "\nu_\tau" tex_anti "\bar\nu_\tau" # The vector bosons particle GLUON 21 parton gauge spin 1 color 8 name gl g G gluon particle PHOTON 22 gauge spin 1 name A gamma photon tex_name "\gamma" particle Z_BOSON 23 gauge spin 1 name Z mass mZ width wZ particle W_BOSON 24 gauge spin 1 charge 1 name "W+" Wp anti "W-" Wm tex_name "W^+" tex_anti "W^-" mass mW width wW # The Higgs particle HIGGS 25 spin 0 name H h Higgs mass mH width wH # Hadrons particle PROTON 2212 spin 1/2 charge 1 name p "p+" anti pbar "p-" # Beam remnants for proton colliders particle HADRON_REMNANT 90 name hr tex_name "had_r" particle HADRON_REMNANT_SINGLET 91 name hr1 tex_name "had_r^{(1)}" particle HADRON_REMNANT_TRIPLET 92 color 3 name hr3 tex_name "had_r^{(3)}" anti hr3bar tex_anti "had_r^{(\bar 3)}" particle HADRON_REMNANT_OCTET 93 color 8 name hr8 tex_name "had_r^{(8)}" ######################################################################## # Vertices of the Standard model with nontrivial CKM matrix # In graphs with identical structure, the first vertex is kept for phase space, # therefore, lighter particles come before heavier ones. # QED vertex D d A vertex U u A vertex S s A vertex C c A vertex B b A vertex T t A vertex E1 e1 A vertex E2 e2 A vertex E3 e3 A # QCD vertex G G G vertex G G G G vertex D d G vertex U u G vertex S s G vertex C c G vertex B b G vertex T t G # Neutral currents vertex D d Z vertex U u Z vertex S s Z vertex C c Z vertex B b Z vertex T t Z vertex E1 e1 Z vertex E2 e2 Z vertex E3 e3 Z vertex N1 n1 Z vertex N2 n2 Z vertex N3 n3 Z # Charged currents vertex N1 e1 Wp vertex N2 e2 Wp vertex N3 e3 Wp vertex E1 n1 Wm vertex E2 n2 Wm vertex E3 n3 Wm vertex U d Wp vertex U s Wp vertex U b Wp vertex C d Wp vertex C s Wp vertex C b Wp vertex T d Wp vertex T s Wp vertex T b Wp vertex D u Wm vertex S u Wm vertex B u Wm vertex D c Wm vertex S c Wm vertex B c Wm vertex D t Wm vertex S t Wm vertex B t Wm # Yukawa (neutral) ### keeping only 3rd generation for the moment # vertex S s H vertex C c H vertex B b H vertex T t H # vertex E2 e2 H vertex E3 e3 H # Vector-boson self-interactions vertex Wp Wm A vertex Wp Wm Z vertex Wp Wm Z Z vertex Wp Wp Wm Wm vertex Wp Wm Z A vertex Wp Wm A A # Higgs - vector boson #vertex H Z A #vertex H A A #vertex H g g vertex H Wp Wm vertex H Z Z vertex H H Wp Wm vertex H H Z Z # Higgs self-interactions vertex H H H vertex H H H H Index: trunk/share/models/MSSM_CKM.mdl =================================================================== --- trunk/share/models/MSSM_CKM.mdl (revision 8479) +++ trunk/share/models/MSSM_CKM.mdl (revision 8480) @@ -1,1307 +1,1307 @@ ######################################################################## # MSSM with nontrivial CKM matrix model "MSSM_CKM" # Independent parameters ### DO NOT CHANGE THE ORDER OF THESE PARAMETERS parameter GF = 1.16639E-5 # Fermi constant parameter mZ = 91.1882 # Z-boson mass parameter wZ = 2.443 # Z-boson width parameter mW = 80.419 # W-boson mass parameter wW = 2.049 # W-boson width parameter me = 0.000511 # electron mass parameter mmu = 0.1057 # muon mass parameter mtau = 1.777 # tau-lepton mass parameter ms = 0.12 # s-quark mass parameter mc = 1.25 # c-quark mass parameter mb = 4.2 # b-quark mass parameter mtop = 174 # t-quark mass parameter wtop = 1.523 # t-quark width parameter alphas = 0.1178 # Strong coupling constant (Z point) # SUSY: These parameters are not used directly by WHIZARD. # They come into play if an external spectrum generator is invoked. parameter mtype = 1 # MSSM embedding model parameter m_zero = 100 # unified scalar mass (SUGRA/AMSB) parameter m_half = 250 # unified gaugino mass (SUGRA) parameter A0 = -100 # unified trilinear coupling (SUGRA) parameter tanb = 10 # tan(beta) = v2/v1 (MSSM input) parameter sgn_mu = 1 # signum(mu) parameter Lambda = 40000 # SUSY breaking scale (GMSB) parameter M_mes = 80000 # messenger scale (GMSB) parameter N_5 = 3 # messenger index (GMSB) parameter c_grav = 1 # gravitino mass multiplier (GMSB) parameter m_grav = 60000 # gravitino mass (AMSB) # The following parameters are used directly by WHIZARD. # They can be derived from the above, but this is not done by WHIZARD parameter Ae_33 = 0 # Ae soft breaking parameter (3rd gen.) parameter Au_33 = 0 # Au soft breaking parameter (3rd gen.) parameter Ad_33 = 0 # Ad soft breaking parameter (3rd gen.) parameter mh = 125 # light Higgs mass parameter wh = 4.143E-3 # light Higgs width parameter mHH = 1000 # heavy Higgs mass parameter mHA = 1000 # axial Higgs mass parameter mHpm = 1000 # charged Higgs mass parameter wHH = 0 # heavy Higgs width parameter wHpm = 0 # charged Higgs width parameter wHA = 0 # axial Higgs width parameter al_h = 0 # Higgs mixing angle alpha parameter mu_h = 1000 # Higgs mu parameter parameter tanb_h = 10 # Higgs mixing angle tan(beta) parameter msu1 = 1000 # u-squark mass parameter msd1 = 1000 # d-squark mass parameter msc1 = 1000 # c-squark mass parameter mss1 = 1000 # s-squark mass parameter mstop1 = 1000 # t-squark mass parameter msb1 = 1000 # b-squark mass parameter msu2 = 1000 # u-squark mass parameter msd2 = 1000 # d-squark mass parameter msc2 = 1000 # c-squark mass parameter mss2 = 1000 # s-squark mass parameter mstop2 = 1000 # t-squark mass parameter msb2 = 1000 # b-squark mass parameter mse1 = 1000 # selectron1 mass parameter msne = 1000 # electron-sneutrino mass parameter msmu1 = 1000 # smuon1 mass parameter msnmu = 1000 # muon-sneutrino mass parameter mstau1 = 1000 # stau1 mass parameter msntau = 1000 # tau-sneutrino mass parameter mse2 = 1000 # selectron2 mass parameter msmu2 = 1000 # smuon2 mass parameter mstau2 = 1000 # stau2 mass parameter mgg = 1000 # gluino mass parameter mch1 = 1000 # chargino1 mass (signed) parameter mch2 = 1000 # chargino2 mass (signed) parameter mneu1 = 1000 # neutralino1 mass (signed) parameter mneu2 = 1000 # neutralino2 mass (signed) parameter mneu3 = 1000 # neutralino3 mass (signed) parameter mneu4 = 1000 # neutralino4 mass (signed) parameter wsu1 = 0 # u-squark width parameter wsd1 = 0 # d-squark width parameter wsc1 = 0 # c-squark width parameter wss1 = 0 # s-squark width parameter wstop1 = 0 # t-squark width parameter wsb1 = 0 # b-squark width parameter wsu2 = 0 # u-squark width parameter wsd2 = 0 # d-squark width parameter wsc2 = 0 # c-squark width parameter wss2 = 0 # s-squark width parameter wstop2 = 0 # t-squark width parameter wsb2 = 0 # b-squark width parameter wse1 = 0 # selectron1 width parameter wsne = 0 # electron-sneutrino width parameter wsmu1 = 0 # smuon1 width parameter wsnmu = 0 # muon-sneutrino width parameter wstau1 = 0 # stau1 width parameter wsntau = 0 # tau-sneutrino width parameter wse2 = 0 # selectron2 width parameter wsmu2 = 0 # smuon2 width parameter wstau2 = 0 # stau2 width parameter wgg = 0 # gluino width parameter wch1 = 0 # chargino1 width parameter wch2 = 0 # chargino2 width parameter wneu1 = 0 # neutralino1 width parameter wneu2 = 0 # neutralino2 width parameter wneu3 = 0 # neutralino3 width parameter wneu4 = 0 # neutralino4 width parameter mt_11 = 1 # stop mixing matrix parameter mt_12 = 0 # stop mixing matrix parameter mt_21 = 0 # stop mixing matrix parameter mt_22 = 1 # stop mixing matrix parameter mb_11 = 1 # sbottom mixing matrix parameter mb_12 = 0 # sbottom mixing matrix parameter mb_21 = 0 # sbottom mixing matrix parameter mb_22 = 1 # sbottom mixing matrix parameter ml_11 = 1 # stau mixing matrix parameter ml_12 = 0 # stau mixing matrix parameter ml_21 = 0 # stau mixing matrix parameter ml_22 = 1 # stau mixing matrix parameter mn_11 = 1 # neutralino mixing matrix parameter mn_12 = 0 # neutralino mixing matrix parameter mn_13 = 0 # neutralino mixing matrix parameter mn_14 = 0 # neutralino mixing matrix parameter mn_21 = 0 # neutralino mixing matrix parameter mn_22 = 1 # neutralino mixing matrix parameter mn_23 = 0 # neutralino mixing matrix parameter mn_24 = 0 # neutralino mixing matrix parameter mn_31 = 0 # neutralino mixing matrix parameter mn_32 = 0 # neutralino mixing matrix parameter mn_33 = 1 # neutralino mixing matrix parameter mn_34 = 0 # neutralino mixing matrix parameter mn_41 = 0 # neutralino mixing matrix parameter mn_42 = 0 # neutralino mixing matrix parameter mn_43 = 0 # neutralino mixing matrix parameter mn_44 = 1 # neutralino mixing matrix parameter mu_11 = 1 # chargino mixing matrix parameter mu_12 = 0 # chargino mixing matrix parameter mu_21 = 0 # chargino mixing matrix parameter mu_22 = 1 # chargino mixing matrix parameter mv_11 = 1 # chargino mixing matrix parameter mv_12 = 0 # chargino mixing matrix parameter mv_21 = 0 # chargino mixing matrix parameter mv_22 = 1 # chargino mixing matrix -parameter vckm11 = 0.97383 # Vud -parameter vckm12 = 0.2272 # Vus -parameter vckm13 = 0.00396 # Vub -parameter vckm21 = -0.2271 # Vcd -parameter vckm22 = 0.97296 # Vcs -parameter vckm23 = 0.04221 # Vcb -parameter vckm31 = 0.00814 # Vtd -parameter vckm32 = -0.04161 # Vts -parameter vckm33 = 0.99910 # Vtb +parameter vckm11 = 0.97401 # Vud +parameter vckm12 = 0.22650 # Vus +parameter vckm13 = 0.00361 # Vub +parameter vckm21 = -0.22636 # Vcd +parameter vckm22 = 0.97320 # Vcs +parameter vckm23 = 0.04053 # Vcb +parameter vckm31 = 0.00854 # Vtd +parameter vckm32 = -0.03978 # Vts +parameter vckm33 = 0.999172 # Vtb # Dependent parameters derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev) derived cw = mW / mZ # cos(theta-W) derived sw = sqrt (1-cw*cw) # sin(theta-W) derived ee = 2 * sw * mW / v # em-coupling (GF scheme) derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const ######################################################################## # Particle content # The quarks particle D_QUARK 1 parton spin 1/2 charge -1/3 isospin -1/2 color 3 name d down anti dbar D "d~" tex_anti "\bar{d}" particle U_QUARK 2 parton spin 1/2 charge 2/3 isospin 1/2 color 3 name u up anti ubar U "u~" tex_anti "\bar{u}" particle S_QUARK 3 like D_QUARK name s strange anti sbar S "s~" tex_anti "\bar{s}" mass ms particle C_QUARK 4 like U_QUARK name c charm anti cbar C "c~" tex_anti "\bar{c}" mass mc particle B_QUARK 5 like D_QUARK name b bottom anti bbar B "b~" tex_anti "\bar{b}" mass mb particle T_QUARK 6 like U_QUARK name t top anti tbar T "t~" tex_anti "\bar{t}" mass mtop width wtop # The leptons particle E_LEPTON 11 spin 1/2 charge -1 isospin -1/2 name "e-" e1 electron e anti "e+" E1 positron tex_name "e^-" tex_anti "e^+" mass me particle E_NEUTRINO 12 left spin 1/2 isospin 1/2 name nue n1 "nu_e" ve "e-neutrino" anti nuebar N1 "ve~" tex_name "\nu_e" tex_anti "\bar\nu_e" particle MU_LEPTON 13 like E_LEPTON name "mu-" e2 mu muon anti "mu+" E2 tex_name "\mu^-" tex_anti "\mu^+" mass mmu particle MU_NEUTRINO 14 like E_NEUTRINO name numu "nu_mu" n2 vm "mu-neutrino" anti numubar N2 "vm~" tex_name "\nu_\mu" tex_anti "\bar\nu_\mu" particle TAU_LEPTON 15 like E_LEPTON name "tau-" e3 tau "ta-" tauon anti "tau+" E3 "ta+" tex_name "\tau^-" tex_anti "\tau^+" mass mtau particle TAU_NEUTRINO 16 like E_NEUTRINO name nutau "nu_tau" n3 vt "tau_neutrino" anti nutaubar N3 "vt~" tex_name "\nu_\tau" tex_anti "\bar\nu_\tau" # The vector bosons particle GLUON 21 parton gauge spin 1 color 8 name gl g G gluon particle PHOTON 22 gauge spin 1 name A gamma photon tex_name "\gamma" particle Z_BOSON 23 gauge spin 1 name Z mass mZ width wZ particle W_BOSON 24 gauge spin 1 charge 1 name "W+" Wp anti "W-" Wm tex_name "W^+" tex_anti "W^-" mass mW width wW # The Higgses particle LIGHT_HIGGS 25 spin 0 name h h0 Higgs tex_name "h^0" mass mh width wh particle HEAVY_HIGGS 35 spin 0 name H HH HH0 H0 tex_name "H^0" mass mHH width wHH particle AXIAL_HIGGS 36 spin 0 name A0 HA HA0 tex_name "A^0" mass mHA width wHA particle CHARGED_HIGGS 37 spin 0 charge +1 name "H+" Hp anti "H-" Hm tex_name "H^+" tex_anti "H^-" mass mHpm width wHpm # The squarks # Left-handed particle D_SQUARK1 1000001 spin 0 charge -1/3 color 3 name sd1 anti sd1c SD1 "sd1~" tex_name "\tilde{d}_L" tex_anti "\tilde{\bar d}_L" mass msd1 width wsd1 particle U_SQUARK1 1000002 spin 0 charge 2/3 color 3 name su1 anti su1c SU1 "su1~" tex_name "\tilde{u}_L" tex_anti "\tilde{\bar u}_L" mass msu1 width wsu1 particle S_SQUARK1 1000003 like D_SQUARK1 name ss1 anti ss1c SS1 "ss1~" tex_name "\tilde{s}_L" tex_anti "\tilde{\bar s}_L" mass mss1 width wss1 particle C_SQUARK1 1000004 like U_SQUARK1 name sc1 anti sc1c SC1 "sc1~" tex_name "\tilde{c}_L" tex_anti "\tilde{\bar c}_L" mass msc1 width wsc1 particle B_SQUARK1 1000005 like D_SQUARK1 name sb1 anti sb1c SB1 "sb1~" tex_name "\tilde{b}_1" tex_anti "\tilde{\bar b}_1" mass msb1 width wsb1 particle T_SQUARK1 1000006 like U_SQUARK1 name st1 anti st1c ST1 "st1~" tex_name "\tilde{t}_1" tex_anti "\tilde{\bar t}_1" mass mstop1 width wstop1 # Right-handed particle D_SQUARK2 2000001 spin 0 charge -1/3 color 3 name sd2 anti sd2c SD2 "sd2~" tex_name "\tilde{d}_R" tex_anti "\tilde{\bar d}_R" mass msd2 width wsd2 particle U_SQUARK2 2000002 spin 0 charge 2/3 color 3 name su2 anti su2c SU2 "su2~" tex_name "\tilde{u}_R" tex_anti "\tilde{\bar u}_R" mass msu2 width wsu2 particle S_SQUARK2 2000003 like D_SQUARK2 name ss2 anti ss2c SS2 "ss2~" tex_name "\tilde{s}_R" tex_anti "\tilde{\bar s}_R" mass mss2 width wss2 particle C_SQUARK2 2000004 like U_SQUARK2 name sc2 anti sc2c SC2 "sc2~" tex_name "\tilde{c}_R" tex_anti "\tilde{\bar c}_R" mass msc2 width wsc2 particle B_SQUARK2 2000005 like D_SQUARK2 name sb2 anti sb2c SB2 "sb2~" tex_name "\tilde{b}_2" tex_anti "\tilde{\bar b}_2" mass msb2 width wsb2 particle T_SQUARK2 2000006 like U_SQUARK2 name st2 anti st2c ST2 "st2~" tex_name "\tilde{t}_2" tex_anti "\tilde{\bar t}_2" mass mstop2 width wstop2 # The sleptons # Left-handed particle E_SLEPTON1 1000011 spin 0 charge -1 name "se1-" se1 se11 anti "se1+" SE1 SE11 tex_name "\tilde{e}_1^-" tex_anti "\tilde{e}_1^+" mass mse1 width wse1 particle E_SNEUTRINO1 1000012 spin 0 name snue "snu_e1" sn11 sve anti "snue*" SN11 "sve~" tex_name "\tilde\nu_e" tex_anti "\tilde\bar\nu_e" mass msne width wsne particle MU_SLEPTON1 1000013 like E_SLEPTON1 name "smu1-" smu1 se21 anti "smu1+" SMU1 SE21 tex_name "\tilde\mu_1^-" tex_anti "\tilde\mu_1^+" mass msmu1 width wsmu1 particle MU_SNEUTRINO1 1000014 like E_SNEUTRINO1 name snumu "snu_mu1" sn21 svm anti "snumu*" SN21 "svm~" tex_name "\tilde\nu_\mu" tex_anti "\tilde\bar\nu_\mu" mass msnmu width wsnmu particle TAU_SLEPTON1 1000015 like E_SLEPTON1 name "stau1-" stau1 se31 anti "stau1+" STAU1 SE31 tex_name "\tilde\tau_1^-" tex_anti "\tilde\tau_1^+" mass mstau1 width wstau1 particle TAU_SNEUTRINO1 1000016 like E_SNEUTRINO1 name snutau "snu_tau" sn31 svt anti "snutau*" SN31 "svt~" tex_name "\tilde\nu_\tau" tex_anti "\tilde\bar\nu_\tau" mass msntau width wsntau # Right-handed particle E_SLEPTON2 2000011 spin 0 charge -1 name "se2-" se2 se12 anti "se2+" SE2 SE12 tex_name "\tilde{e}_2^-" tex_anti "\tilde{e}_2^+" mass mse2 width wse2 particle MU_SLEPTON2 2000013 like E_SLEPTON2 name "smu2-" smu2 se22 anti "smu2+" SMU2 SE22 tex_name "\tilde\mu_2^-" tex_anti "\tilde\mu_2^+" mass msmu2 width wsmu2 particle TAU_SLEPTON2 2000015 like E_SLEPTON2 name "stau2-" stau2 se32 anti "stau2+" STAU2 SE32 tex_name "\tilde\tau_2^-" tex_anti "\tilde\tau_2^+" mass mstau2 width wstau2 # The gauginos particle GLUINO 1000021 spin 1/2 color 8 name sgl gg GG gluino tex_name "\tilde{g}" mass mgg width wgg particle CHARGINO1 1000024 spin 1/2 charge 1 name "ch1+" "CH1+" anti "ch1-" "CH1-" tex_name "\tilde\chi_1^+" tex_anti "\tilde\chi_1^-" mass mch1 width wch1 particle CHARGINO2 1000037 like CHARGINO1 name "ch2+" "CH2+" anti "ch2-" "CH2-" tex_name "\tilde\chi_2^+" tex_anti "\tilde\chi_2^-" mass mch2 width wch2 particle NEUTRALINO1 1000022 spin 1/2 name neu1 NEU1 tex_name "\tilde\chi_1^0" mass mneu1 width wneu1 particle NEUTRALINO2 1000023 like NEUTRALINO1 name neu2 NEU2 tex_name "\tilde\chi_2^0" mass mneu2 width wneu2 particle NEUTRALINO3 1000025 like NEUTRALINO1 name neu3 NEU3 tex_name "\tilde\chi_3^0" mass mneu3 width wneu3 particle NEUTRALINO4 1000035 like NEUTRALINO1 name neu4 NEU4 tex_name "\tilde\chi_4^0" mass mneu4 width wneu4 # Hadrons particle PROTON 2212 spin 1/2 charge 1 name p "p+" anti pbar "p-" # Beam remnants for proton colliders particle HADRON_REMNANT 90 name hr tex_name "had_r" particle HADRON_REMNANT_SINGLET 91 name hr1 tex_name "had_r^{(1)}" particle HADRON_REMNANT_TRIPLET 92 color 3 name hr3 tex_name "had_r^{(3)}" anti hr3bar tex_anti "had_r^{(\bar 3)}" particle HADRON_REMNANT_OCTET 93 color 8 name hr8 tex_name "had_r^{(8)}" ######################################################################## # Vertices of the MSSM with nontrivial CKM matrix # In graphs with identical structure, the first vertex is kept for phase space, # therefore, lighter particles come before heavier ones. ! QED vertex D d A vertex U u A vertex S s A vertex C c A vertex B b A vertex T t A vertex E1 e1 A vertex E2 e2 A vertex E3 e3 A ! QED/SUSY vertex sd1 SD1 A vertex su1 SU1 A vertex ss1 SS1 A vertex sc1 SC1 A vertex sb1 SB1 A vertex st1 ST1 A vertex se11 SE11 A vertex se21 SE21 A vertex se31 SE31 A vertex sd2 SD2 A vertex su2 SU2 A vertex ss2 SS2 A vertex sc2 SC2 A vertex sb2 SB2 A vertex st2 ST2 A vertex se12 SE12 A vertex se22 SE22 A vertex se32 SE32 A # QCD vertex G G G vertex G G G G ! SUSY QCD vertex g gg gg ! QCD fermion-gluon vertex D d G vertex U u G vertex S s G vertex C c G vertex B b G vertex T t G ! SUSY sfermion-gluon vertex sd1 SD1 g vertex sd2 SD2 g vertex su1 SU1 g vertex su2 SU2 g vertex ss1 SS1 g vertex ss2 SS2 g vertex sc1 SC1 g vertex sc2 SC2 g vertex sb1 SB1 g vertex sb2 SB2 g vertex st1 ST1 g vertex st2 ST2 g ! SUSY gluino vertex d SD1 gg vertex d SD2 gg vertex u SU1 gg vertex u SU2 gg vertex s SS1 gg vertex s SS2 gg vertex c SC1 gg vertex c SC2 gg vertex b SB1 gg vertex b SB2 gg vertex t ST1 gg vertex t ST2 gg vertex sd1 D gg vertex sd2 D gg vertex su1 U gg vertex su2 U gg vertex ss1 S gg vertex ss2 S gg vertex sc1 C gg vertex sc2 C gg vertex sb1 B gg vertex sb2 B gg vertex st1 T gg vertex st2 T gg # Neutral currents vertex D d Z vertex U u Z vertex S s Z vertex C c Z vertex B b Z vertex T t Z vertex E1 e1 Z vertex E2 e2 Z vertex E3 e3 Z vertex N1 n1 Z vertex N2 n2 Z vertex N3 n3 Z ! Neutral currents/SUSY vertex sd1 SD1 Z vertex su1 SU1 Z vertex ss1 SS1 Z vertex sc1 SC1 Z vertex sb1 SB1 Z vertex st1 ST1 Z vertex se11 SE11 Z vertex se21 SE21 Z vertex se31 SE31 Z vertex sn11 SN11 Z vertex sn21 SN21 Z vertex sn31 SN31 Z vertex sd2 SD2 Z vertex su2 SU2 Z vertex ss2 SS2 Z vertex sc2 SC2 Z vertex sb2 SB2 Z vertex st2 ST2 Z vertex se12 SE12 Z vertex se22 SE22 Z vertex se32 SE32 Z ! 3rd gen mixing vertex sb1 SB2 Z vertex sb2 SB1 Z vertex st1 ST2 Z vertex st2 ST1 Z vertex se31 SE32 Z vertex se32 SE31 Z # Charged currents vertex U d Wp vertex U s Wp vertex U b Wp vertex C d Wp vertex C s Wp vertex C b Wp vertex T d Wp vertex T s Wp vertex T b Wp vertex D u Wm vertex S u Wm vertex B u Wm vertex D c Wm vertex S c Wm vertex B c Wm vertex D t Wm vertex S t Wm vertex B t Wm vertex N1 e1 Wp vertex N2 e2 Wp vertex N3 e3 Wp vertex E1 n1 Wm vertex E2 n2 Wm vertex E3 n3 Wm ! SUSY vertex su1 SD1 Wm vertex su1 SS1 Wm vertex su1 SB1 Wm vertex sc1 SD1 Wm vertex sc1 SS1 Wm vertex sc1 SB1 Wm vertex st1 SD1 Wm vertex st1 SS1 Wm vertex st1 SB1 Wm vertex su2 SD1 Wm vertex su2 SS1 Wm vertex su2 SB1 Wm vertex sc2 SD1 Wm vertex sc2 SS1 Wm vertex sc2 SB1 Wm vertex st2 SD1 Wm vertex st2 SS1 Wm vertex st2 SB1 Wm vertex su1 SD2 Wm vertex su1 SS2 Wm vertex su1 SB2 Wm vertex sc1 SD2 Wm vertex sc1 SS2 Wm vertex sc1 SB2 Wm vertex st1 SD2 Wm vertex st1 SS2 Wm vertex st1 SB2 Wm vertex su2 SD2 Wm vertex su2 SS2 Wm vertex su2 SB2 Wm vertex sc2 SD2 Wm vertex sc2 SS2 Wm vertex sc2 SB2 Wm vertex st2 SD2 Wm vertex st2 SS2 Wm vertex st2 SB2 Wm vertex sn11 SE11 Wm vertex sn21 SE21 Wm vertex sn31 SE31 Wm vertex sd1 SU1 Wp vertex sb1 SU1 Wp vertex sd1 SU1 Wp vertex sd1 SC1 Wp vertex ss1 SC1 Wp vertex sb1 SC1 Wp vertex sd1 ST1 Wp vertex ss1 ST1 Wp vertex sb1 ST1 Wp vertex sd1 SU2 Wp vertex sb1 SU2 Wp vertex sd1 SU2 Wp vertex sd1 SC2 Wp vertex ss1 SC2 Wp vertex sb1 SC2 Wp vertex sd1 ST2 Wp vertex ss1 ST2 Wp vertex sb1 ST2 Wp vertex sd2 SU1 Wp vertex sb2 SU1 Wp vertex sd2 SU1 Wp vertex sd2 SC1 Wp vertex ss2 SC1 Wp vertex sb2 SC1 Wp vertex sd2 ST1 Wp vertex ss2 ST1 Wp vertex sb2 ST1 Wp vertex sd2 SU2 Wp vertex sb2 SU2 Wp vertex sd2 SU2 Wp vertex sd2 SC2 Wp vertex ss2 SC2 Wp vertex sb2 SC2 Wp vertex sd2 ST2 Wp vertex ss2 ST2 Wp vertex sb2 ST2 Wp vertex se11 SN11 Wp vertex se21 SN21 Wp vertex se31 SN31 Wp ! 3rd gen mixing (only sleptons here) vertex sn31 SE32 Wm vertex se32 SN31 Wp # Yukawa (neutral) ### keeping only 3rd generation for the moment # vertex S s h # vertex C c h vertex B b h vertex T t h # vertex E2 e2 h vertex E3 e3 h ! SUSY Higgses vertex e3 E3 HH vertex b B HH vertex t T HH vertex e3 E3 HA vertex b B HA vertex t T HA vertex e3 N3 "H+" vertex E3 n3 "H-" vertex b T "H+" vertex B t "H-" ! SUSY-Yukawa vertex se31 SE31 h vertex se31 SE32 h vertex se32 SE31 h vertex se32 SE32 h vertex sb1 SB1 h vertex sb1 SB2 h vertex sb2 SB1 h vertex sb2 SB2 h vertex st1 ST1 h vertex st1 ST2 h vertex st2 ST1 h vertex st2 ST2 h vertex se31 SE31 HH vertex se31 SE32 HH vertex se32 SE31 HH vertex se32 SE32 HH vertex sb1 SB1 HH vertex sb1 SB2 HH vertex sb2 SB1 HH vertex sb2 SB2 HH vertex st1 ST1 HH vertex st1 ST2 HH vertex st2 ST1 HH vertex st2 ST2 HH vertex se31 SE31 HA vertex se31 SE32 HA vertex se32 SE31 HA vertex se32 SE32 HA vertex sb1 SB1 HA vertex sb1 SB2 HA vertex sb2 SB1 HA vertex sb2 SB2 HA vertex st1 ST1 HA vertex st1 ST2 HA vertex st2 ST1 HA vertex st2 ST2 HA vertex se31 SN31 "H+" vertex se32 SN31 "H+" vertex sb1 ST1 "H+" vertex sb1 ST2 "H+" vertex sb2 ST1 "H+" vertex sb2 ST2 "H+" vertex SE31 sn31 "H-" vertex SE32 sn31 "H-" vertex SB1 st1 "H-" vertex SB1 st2 "H-" vertex SB2 st1 "H-" vertex SB2 st2 "H-" # Vector-boson self-interactions vertex Wp Wm A vertex Wp Wm Z vertex Wp Wm Z Z vertex Wp Wp Wm Wm vertex Wp Wm A Z vertex Wp Wm A A # Higgs - vector boson vertex A "H+" "H-" vertex Z "H+" "H-" vertex Z h HA vertex Z HH HA vertex Wp "H-" h vertex Wm "H+" h vertex Wp "H-" HH vertex Wm "H+" HH vertex Wp "H-" HA vertex Wm "H+" HA vertex h Wp Wm vertex h Z Z vertex Wp Wm HH vertex Z Z HH vertex A "H+" "H-" vertex Z Z HH vertex Z h HA vertex Z HH HA vertex Z "H+" "H-" # Higgs self-interactions vertex h h HH vertex h HH HH vertex HH HH HH vertex h HA HA vertex HH HA HA vertex h "H+" "H-" vertex HH "H+" "H-" ! Charginos vertex "ch1+" "ch1-" A vertex "ch2+" "ch2-" A vertex "ch1+" "ch1-" Z vertex "ch1+" "ch2-" Z vertex "ch2+" "ch1-" Z vertex "ch2+" "ch2-" Z ! Neutralinos vertex neu1 neu1 Z vertex neu1 neu2 Z vertex neu1 neu3 Z vertex neu1 neu4 Z vertex neu2 neu2 Z vertex neu2 neu3 Z vertex neu2 neu4 Z vertex neu3 neu3 Z vertex neu3 neu4 Z vertex neu4 neu4 Z ! Charginos + neutralinos vertex "ch1+" neu1 Wm vertex "ch2+" neu1 Wm vertex "ch1+" neu2 Wm vertex "ch2+" neu2 Wm vertex "ch1+" neu3 Wm vertex "ch2+" neu3 Wm vertex "ch1+" neu4 Wm vertex "ch2+" neu4 Wm vertex "ch1-" neu1 Wp vertex "ch2-" neu1 Wp vertex "ch1-" neu2 Wp vertex "ch2-" neu2 Wp vertex "ch1-" neu3 Wp vertex "ch2-" neu3 Wp vertex "ch1-" neu4 Wp vertex "ch2-" neu4 Wp ! Charginos + Higgs vertex "ch1+" "ch1-" h vertex "ch1+" "ch2-" h vertex "ch2+" "ch1-" h vertex "ch2+" "ch2-" h vertex "ch1+" "ch1-" HH vertex "ch1+" "ch2-" HH vertex "ch2+" "ch1-" HH vertex "ch2+" "ch2-" HH vertex "ch1+" "ch1-" HA vertex "ch1+" "ch2-" HA vertex "ch2+" "ch1-" HA vertex "ch2+" "ch2-" HA ! Neutralinos + Higgs vertex neu1 neu1 h vertex neu1 neu2 h vertex neu1 neu3 h vertex neu1 neu4 h vertex neu2 neu2 h vertex neu2 neu3 h vertex neu2 neu4 h vertex neu3 neu3 h vertex neu3 neu4 h vertex neu4 neu4 h vertex neu1 neu1 HH vertex neu1 neu2 HH vertex neu1 neu3 HH vertex neu1 neu4 HH vertex neu2 neu2 HH vertex neu2 neu3 HH vertex neu2 neu4 HH vertex neu3 neu3 HH vertex neu3 neu4 HH vertex neu4 neu4 HH vertex neu1 neu1 HA vertex neu1 neu2 HA vertex neu1 neu3 HA vertex neu1 neu4 HA vertex neu2 neu2 HA vertex neu2 neu3 HA vertex neu2 neu4 HA vertex neu3 neu3 HA vertex neu3 neu4 HA vertex neu4 neu4 HA ! Charginos + neutralinos + Higgs vertex "ch1+" neu1 "H-" vertex "ch2+" neu1 "H-" vertex "ch1+" neu2 "H-" vertex "ch2+" neu2 "H-" vertex "ch1+" neu3 "H-" vertex "ch2+" neu3 "H-" vertex "ch1+" neu4 "H-" vertex "ch2+" neu4 "H-" vertex "ch1-" neu1 "H+" vertex "ch2-" neu1 "H+" vertex "ch1-" neu1 "H+" vertex "ch2-" neu2 "H+" vertex "ch1-" neu3 "H+" vertex "ch2-" neu3 "H+" vertex "ch1-" neu4 "H+" vertex "ch2-" neu4 "H+" ! Lepton-slepton-neutralino vertex e1 SE11 neu1 vertex e1 SE11 neu2 vertex e1 SE11 neu3 vertex e1 SE11 neu4 vertex e1 SE12 neu1 vertex e1 SE12 neu2 vertex e1 SE12 neu3 vertex e1 SE12 neu4 vertex E1 se11 neu1 vertex E1 se11 neu2 vertex E1 se11 neu3 vertex E1 se11 neu4 vertex E1 se12 neu1 vertex E1 se12 neu2 vertex E1 se12 neu3 vertex E1 se12 neu4 vertex e2 SE21 neu1 vertex e2 SE21 neu2 vertex e2 SE21 neu3 vertex e2 SE21 neu4 vertex E2 se21 neu1 vertex E2 se21 neu2 vertex E2 se21 neu3 vertex E2 se21 neu4 vertex e2 SE22 neu1 vertex e2 SE22 neu2 vertex e2 SE22 neu3 vertex e2 SE22 neu4 vertex E2 se22 neu1 vertex E2 se22 neu2 vertex E2 se22 neu3 vertex E2 se22 neu4 vertex e3 SE31 neu1 vertex e3 SE31 neu2 vertex e3 SE31 neu3 vertex e3 SE31 neu4 vertex E3 se31 neu1 vertex E3 se31 neu2 vertex E3 se31 neu3 vertex E3 se31 neu4 vertex e3 SE32 neu1 vertex e3 SE32 neu2 vertex e3 SE32 neu3 vertex e3 SE32 neu4 vertex E3 se32 neu1 vertex E3 se32 neu2 vertex E3 se32 neu3 vertex E3 se32 neu4 ! Neutrino-sneutrino-neutralino vertex n1 SN11 neu1 vertex n1 SN11 neu2 vertex n1 SN11 neu3 vertex n1 SN11 neu4 vertex N1 sn11 neu1 vertex N1 sn11 neu2 vertex N1 sn11 neu3 vertex N1 sn11 neu4 vertex n2 SN21 neu1 vertex n2 SN21 neu2 vertex n2 SN21 neu3 vertex n2 SN21 neu4 vertex N2 sn21 neu1 vertex N2 sn21 neu2 vertex N2 sn21 neu3 vertex N2 sn21 neu4 vertex n3 SN31 neu1 vertex n3 SN31 neu2 vertex n3 SN31 neu3 vertex n3 SN31 neu4 vertex N3 sn31 neu1 vertex N3 sn31 neu2 vertex N3 sn31 neu3 vertex N3 sn31 neu4 ! Quark-squark-neutralino vertex d SD1 neu1 vertex d SD1 neu2 vertex d SD1 neu3 vertex d SD1 neu4 vertex d SD2 neu1 vertex d SD2 neu2 vertex d SD2 neu3 vertex d SD2 neu4 vertex D sd1 neu1 vertex D sd1 neu2 vertex D sd1 neu3 vertex D sd1 neu4 vertex D sd2 neu1 vertex D sd2 neu2 vertex D sd2 neu3 vertex D sd2 neu4 vertex u SU1 neu1 vertex u SU1 neu2 vertex u SU1 neu3 vertex u SU1 neu4 vertex u SU2 neu1 vertex u SU2 neu2 vertex u SU2 neu3 vertex u SU2 neu4 vertex U su1 neu1 vertex U su1 neu2 vertex U su1 neu3 vertex U su1 neu4 vertex U su2 neu1 vertex U su2 neu2 vertex U su2 neu3 vertex U su2 neu4 vertex s SS1 neu1 vertex s SS1 neu2 vertex s SS1 neu3 vertex s SS1 neu4 vertex s SS2 neu1 vertex s SS2 neu2 vertex s SS2 neu3 vertex s SS2 neu4 vertex S ss1 neu1 vertex S ss1 neu2 vertex S ss1 neu3 vertex S ss1 neu4 vertex S ss2 neu1 vertex S ss2 neu2 vertex S ss2 neu3 vertex S ss2 neu4 vertex c SC1 neu1 vertex c SC1 neu2 vertex c SC1 neu3 vertex c SC1 neu4 vertex c SC2 neu1 vertex c SC2 neu2 vertex c SC2 neu3 vertex c SC2 neu4 vertex C sc1 neu1 vertex C sc1 neu2 vertex C sc1 neu3 vertex C sc1 neu4 vertex C sc2 neu1 vertex C sc2 neu2 vertex C sc2 neu3 vertex C sc2 neu4 vertex b SB1 neu1 vertex b SB1 neu2 vertex b SB1 neu3 vertex b SB1 neu4 vertex b SB2 neu1 vertex b SB2 neu2 vertex b SB2 neu3 vertex b SB2 neu4 vertex B sb1 neu1 vertex B sb1 neu2 vertex B sb1 neu3 vertex B sb1 neu4 vertex B sb2 neu1 vertex B sb2 neu2 vertex B sb2 neu3 vertex B sb2 neu4 vertex t ST1 neu1 vertex t ST1 neu2 vertex t ST1 neu3 vertex t ST1 neu4 vertex t ST2 neu1 vertex t ST2 neu2 vertex t ST2 neu3 vertex t ST2 neu4 vertex T st1 neu1 vertex T st1 neu2 vertex T st1 neu3 vertex T st1 neu4 vertex T st2 neu1 vertex T st2 neu2 vertex T st2 neu3 vertex T st2 neu4 ! Lepton-sneutrino-chargino vertex e1 SN11 "ch1+" vertex e1 SN11 "ch2+" vertex E1 sn11 "ch1-" vertex E1 sn11 "ch2-" vertex e2 SN21 "ch1+" vertex e2 SN21 "ch2+" vertex E2 sn21 "ch1-" vertex E2 sn21 "ch2-" vertex e3 SN31 "ch1+" vertex e3 SN31 "ch2+" vertex E3 sn31 "ch1-" vertex E3 sn31 "ch2-" ! Slepton-neutrino-chargino vertex se11 N1 "ch1+" vertex se11 N1 "ch2+" vertex se12 N1 "ch1+" vertex se12 N1 "ch2+" vertex SE11 n1 "ch1-" vertex SE11 n1 "ch2-" vertex SE12 n1 "ch1-" vertex SE12 n1 "ch2-" vertex se21 N2 "ch1+" vertex se21 N2 "ch2+" vertex se22 N2 "ch1+" vertex se22 N2 "ch2+" vertex SE21 n2 "ch1-" vertex SE21 n2 "ch2-" vertex SE22 n2 "ch1-" vertex SE22 n2 "ch2-" vertex se31 N3 "ch1+" vertex se31 N3 "ch2+" vertex se32 N3 "ch1+" vertex se32 N3 "ch2+" vertex SE31 n3 "ch1-" vertex SE31 n3 "ch2-" vertex SE32 n3 "ch1-" vertex SE32 n3 "ch2-" ! Quark-squark-chargino vertex d SU1 "ch1+" vertex d SU1 "ch2+" vertex d SU2 "ch1+" vertex d SU2 "ch2+" vertex s SU1 "ch1+" vertex s SU1 "ch2+" vertex s SU2 "ch1+" vertex s SU2 "ch2+" vertex b SU1 "ch1+" vertex b SU1 "ch2+" vertex b SU2 "ch1+" vertex b SU2 "ch2+" vertex D su1 "ch1-" vertex D su1 "ch2-" vertex D su2 "ch1-" vertex D su2 "ch2-" vertex S su1 "ch1-" vertex S su1 "ch2-" vertex S su2 "ch1-" vertex S su2 "ch2-" vertex B su1 "ch1-" vertex B su1 "ch2-" vertex B su2 "ch1-" vertex B su2 "ch2-" vertex sd1 U "ch1+" vertex sd1 U "ch2+" vertex sd2 U "ch1+" vertex sd2 U "ch2+" vertex ss1 U "ch1+" vertex ss1 U "ch2+" vertex ss2 U "ch1+" vertex ss2 U "ch2+" vertex sb1 U "ch1+" vertex sb1 U "ch2+" vertex sb2 U "ch1+" vertex sb2 U "ch2+" vertex SD1 u "ch1-" vertex SD1 u "ch2-" vertex SD2 u "ch1-" vertex SD2 u "ch2-" vertex SS1 u "ch1-" vertex SS1 u "ch2-" vertex SS2 u "ch1-" vertex SS2 u "ch2-" vertex SB1 u "ch1-" vertex SB1 u "ch2-" vertex SB2 u "ch1-" vertex SB2 u "ch2-" vertex d SC1 "ch1+" vertex d SC1 "ch2+" vertex d SC2 "ch1+" vertex d SC2 "ch2+" vertex s SC1 "ch1+" vertex s SC1 "ch2+" vertex s SC2 "ch1+" vertex s SC2 "ch2+" vertex b SC1 "ch1+" vertex b SC1 "ch2+" vertex b SC2 "ch1+" vertex b SC2 "ch2+" vertex D sc1 "ch1-" vertex D sc1 "ch2-" vertex D sc2 "ch1-" vertex D sc2 "ch2-" vertex S sc1 "ch1-" vertex S sc1 "ch2-" vertex S sc2 "ch1-" vertex S sc2 "ch2-" vertex B sc1 "ch1-" vertex B sc1 "ch2-" vertex B sc2 "ch1-" vertex B sc2 "ch2-" vertex sd1 C "ch1+" vertex sd1 C "ch2+" vertex sd2 C "ch1+" vertex sd2 C "ch2+" vertex ss1 C "ch1+" vertex ss1 C "ch2+" vertex ss2 C "ch1+" vertex ss2 C "ch2+" vertex sb1 C "ch1+" vertex sb1 C "ch2+" vertex sb2 C "ch1+" vertex sb2 C "ch2+" vertex SD1 c "ch1-" vertex SD1 c "ch2-" vertex SD2 c "ch1-" vertex SD2 c "ch2-" vertex SS1 c "ch1-" vertex SS1 c "ch2-" vertex SS2 c "ch1-" vertex SS2 c "ch2-" vertex SB1 c "ch1-" vertex SB1 c "ch2-" vertex SB2 c "ch1-" vertex SB2 c "ch2-" vertex d ST1 "ch1+" vertex d ST1 "ch2+" vertex d ST2 "ch1+" vertex d ST2 "ch2+" vertex s ST1 "ch1+" vertex s ST1 "ch2+" vertex s ST2 "ch1+" vertex s ST2 "ch2+" vertex b ST1 "ch1+" vertex b ST1 "ch2+" vertex b ST2 "ch1+" vertex b ST2 "ch2+" vertex D st1 "ch1-" vertex D st1 "ch2-" vertex D st2 "ch1-" vertex D st2 "ch2-" vertex S st1 "ch1-" vertex S st1 "ch2-" vertex S st2 "ch1-" vertex S st2 "ch2-" vertex B st1 "ch1-" vertex B st1 "ch2-" vertex B st2 "ch1-" vertex B st2 "ch2-" vertex sd1 T "ch1+" vertex sd1 T "ch2+" vertex sd2 T "ch1+" vertex sd2 T "ch2+" vertex ss1 T "ch1+" vertex ss1 T "ch2+" vertex ss2 T "ch1+" vertex ss2 T "ch2+" vertex sb1 T "ch1+" vertex sb1 T "ch2+" vertex sb2 T "ch1+" vertex sb2 T "ch2+" vertex SD1 t "ch1-" vertex SD1 t "ch2-" vertex SD2 t "ch1-" vertex SD2 t "ch2-" vertex SS1 t "ch1-" vertex SS1 t "ch2-" vertex SS2 t "ch1-" vertex SS2 t "ch2-" vertex SB1 t "ch1-" vertex SB1 t "ch2-" vertex SB2 t "ch1-" vertex SB2 t "ch2-" Index: trunk/share/models/SM_ac_CKM.mdl =================================================================== --- trunk/share/models/SM_ac_CKM.mdl (revision 8479) +++ trunk/share/models/SM_ac_CKM.mdl (revision 8480) @@ -1,302 +1,302 @@ ######################################################################## # Standard Model (unitary gauge) with anomalous couplings # and nontrivial CKM matrix # Default Higgs mass is very large model "SM_ac_CKM" # Independent parameters ### DO NOT CHANGE THE ORDER OF THESE PARAMETERS parameter GF = 1.16639E-5 # Fermi constant parameter mZ = 91.1882 # Z-boson mass parameter mW = 80.419 # W-boson mass parameter mH = 125 # Higgs mass parameter alphas = 0.1178 # Strong coupling constant (Z point) parameter me = 0.000511 # electron mass parameter mmu = 0.1057 # muon mass parameter mtau = 1.777 # tau-lepton mass parameter ms = 0.12 # s-quark mass parameter mc = 1.25 # c-quark mass parameter mb = 4.2 # b-quark mass parameter mtop = 174 # t-quark mass parameter wtop = 1.523 # t-quark width parameter wZ = 2.443 # Z-boson width parameter wW = 2.049 # W-boson width parameter wH = 0.004143 # Higgs width parameter xi0 = 0 # R_xi parameter for Z-boson parameter xipm = 0 # R_xi parameter for W-boson -# Here are the values from PDG 2006 -parameter vckm11 = 0.97383 # Vud -parameter vckm12 = 0.2272 # Vus -parameter vckm13 = 0.00396 # Vub -parameter vckm21 = -0.2271 # Vcd -parameter vckm22 = 0.97296 # Vcs -parameter vckm23 = 0.04221 # Vcb -parameter vckm31 = 0.00814 # Vtd -parameter vckm32 = -0.04161 # Vts -parameter vckm33 = 0.99910 # Vtb +# Here are the values from PDG 2020 +parameter vckm11 = 0.97401 # Vud +parameter vckm12 = 0.22650 # Vus +parameter vckm13 = 0.00361 # Vub +parameter vckm21 = -0.22636 # Vcd +parameter vckm22 = 0.97320 # Vcs +parameter vckm23 = 0.04053 # Vcb +parameter vckm31 = 0.00854 # Vtd +parameter vckm32 = -0.03978 # Vts +parameter vckm33 = 0.999172 # Vtb parameter a4 = 0 # Coefficient of L4 parameter a5 = 0 # Coefficient of L5 parameter a6 = 0 # Coefficient of L6 parameter a7 = 0 # Coefficient of L7 parameter a10 = 0 # Coefficient of L10 parameter g1a = 1 # anomalous TGC parameter g1z = 1 # anomalous TGC parameter g4a = 0 # anomalous TGC parameter g4z = 0 # anomalous TGC parameter g5a = 0 # anomalous TGC parameter g5z = 0 # anomalous TGC parameter ka = 1 # anomalous TGC parameter kz = 1 # anomalous TGC parameter la = 0 # anomalous TGC parameter lz = 0 # anomalous TGC parameter k5a = 0 # anomalous TGC parameter k5z = 0 # anomalous TGC parameter l5a = 0 # anomalous TGC parameter l5z = 0 # anomalous TGC parameter fac_gh3 = 1 # gh3 / gh3_SM parameter fghgaga = 1 # ghgaga / ghgaga_SM parameter fghgaz = 1 # ghgaz / ghgaz_SM parameter lambda_h = 5.e3 # scale for anom. Higgs couplings parameter fw = 0 # anomalous Higgs coupling parameter fww = 0 # anomalous Higgs coupling parameter fb = 0 # anomalous Higgs coupling parameter fbb = 0 # anomalous Higgs coupling # Dependent parameters derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev) derived cw = mW / mZ # cos(theta-W) derived sw = sqrt (1-cw**2) # sin(theta-W) derived ee = 2 * sw * mW / v # em-coupling (GF scheme) derived csw = cw**2-sw**2 # Combination of CW and SW derived aZ = a4+a5+(a6+a7+a10)*2 # Combination of Coefficients in 4Z vertex derived aWZ1 = a4+a6 # Combination of Coefficients in WWZZ vertex derived aWZ2 = a5+a7 # Combination of Coefficients in WWZZ vertex derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const ######################################################################## # Particle content # The quarks particle D_QUARK 1 parton spin 1/2 charge -1/3 isospin -1/2 color 3 name d down anti dbar D "d~" tex_anti "\bar{d}" particle U_QUARK 2 parton spin 1/2 charge 2/3 isospin 1/2 color 3 name u up anti ubar U "u~" tex_anti "\bar{u}" particle S_QUARK 3 like D_QUARK name s strange anti sbar S "s~" tex_anti "\bar{s}" mass ms particle C_QUARK 4 like U_QUARK name c charm anti cbar C "c~" tex_anti "\bar{c}" mass mc particle B_QUARK 5 like D_QUARK name b bottom anti bbar B "b~" tex_anti "\bar{b}" mass mb particle T_QUARK 6 like U_QUARK name t top anti tbar T "t~" tex_anti "\bar{t}" mass mtop width wtop # The leptons particle E_LEPTON 11 spin 1/2 charge -1 isospin -1/2 name "e-" e1 electron e anti "e+" E1 positron tex_name "e^-" tex_anti "e^+" mass me particle E_NEUTRINO 12 left spin 1/2 isospin 1/2 name nue n1 "nu_e" ve "e-neutrino" anti nuebar N1 "ve~" tex_name "\nu_e" tex_anti "\bar\nu_e" particle MU_LEPTON 13 like E_LEPTON name "mu-" e2 mu muon anti "mu+" E2 tex_name "\mu^-" tex_anti "\mu^+" mass mmu particle MU_NEUTRINO 14 like E_NEUTRINO name numu "nu_mu" n2 vm "mu-neutrino" anti numubar N2 "vm~" tex_name "\nu_\mu" tex_anti "\bar\nu_\mu" particle TAU_LEPTON 15 like E_LEPTON name "tau-" e3 tau "ta-" tauon anti "tau+" E3 "ta+" tex_name "\tau^-" tex_anti "\tau^+" mass mtau particle TAU_NEUTRINO 16 like E_NEUTRINO name nutau "nu_tau" n3 vt "tau_neutrino" anti nutaubar N3 "vt~" tex_name "\nu_\tau" tex_anti "\bar\nu_\tau" # The vector bosons particle GLUON 21 parton gauge spin 1 color 8 name gl g G gluon particle PHOTON 22 gauge spin 1 name A gamma photon tex_name "\gamma" particle Z_BOSON 23 gauge spin 1 name Z mass mZ width wZ particle W_BOSON 24 gauge spin 1 charge 1 name "W+" Wp anti "W-" Wm tex_name "W^+" tex_anti "W^-" mass mW width wW # The Higgs particle HIGGS 25 spin 0 name H h Higgs mass mH width wH # Hadrons particle PROTON 2212 spin 1/2 charge 1 name p "p+" anti pbar "p-" # Beam remnants for proton colliders particle HADRON_REMNANT 90 name hr tex_name "had_r" particle HADRON_REMNANT_SINGLET 91 name hr1 tex_name "had_r^{(1)}" particle HADRON_REMNANT_TRIPLET 92 color 3 name hr3 tex_name "had_r^{(3)}" anti hr3bar tex_anti "had_r^{(\bar 3)}" particle HADRON_REMNANT_OCTET 93 color 8 name hr8 tex_name "had_r^{(8)}" ######################################################################## # Vertices of the Standard model with anomalous couplings # In graphs with identical structure, the first vertex is kept for phase space, # therefore, lighter particles come before heavier ones. # # Note: The electric charge has opposite sign compared to SM.mdl # QED vertex D d A vertex U u A vertex S s A vertex C c A vertex B b A vertex T t A vertex E1 e1 A vertex E2 e2 A vertex E3 e3 A # QCD vertex G G G vertex D d G vertex U u G vertex S s G vertex C c G vertex B b G vertex T t G # Neutral currents vertex D d Z vertex U u Z vertex S s Z vertex C c Z vertex B b Z vertex T t Z vertex E1 e1 Z vertex E2 e2 Z vertex E3 e3 Z vertex N1 n1 Z vertex N2 n2 Z vertex N3 n3 Z # Charged currents vertex U d Wp vertex U s Wp vertex U b Wp vertex C d Wp vertex C s Wp vertex C b Wp vertex T d Wp vertex T s Wp vertex T b Wp vertex D u Wm vertex S u Wm vertex B u Wm vertex D c Wm vertex S c Wm vertex B c Wm vertex D t Wm vertex S t Wm vertex B t Wm vertex N1 e1 Wp vertex N2 e2 Wp vertex N3 e3 Wp vertex E1 n1 Wm vertex E2 n2 Wm vertex E3 n3 Wm # Yukawa ### keeping only 3rd generation for the moment # vertex S s H # vertex C c H vertex B b H vertex T t H # vertex E2 e2 H vertex E3 e3 H # Vector-boson self-interactions # These have been generated automatically, therefore unsorted vertex A Wm Wp vertex Wm Wp Z vertex A A Wm Wp vertex A Wm Wp Z vertex Wm Wm Wp Wp vertex Wm Wp Z Z vertex Z Z Z Z # Higgs - vector boson vertex H A A vertex H Z A vertex H Wp Wm vertex H Z Z vertex H H Z Z vertex H H Wp Wm # Higgs self-interactions vertex H H H vertex H H H H Index: trunk/share/models/NMSSM_CKM.mdl =================================================================== --- trunk/share/models/NMSSM_CKM.mdl (revision 8479) +++ trunk/share/models/NMSSM_CKM.mdl (revision 8480) @@ -1,1577 +1,1577 @@ ######################################################################## # NMSSM with non-trivial CKM matrix model "NMSSM_CKM" # Independent parameters ### DO NOT CHANGE THE ORDER OF THESE PARAMETERS parameter GF = 1.16639E-5 # Fermi constant parameter mZ = 91.1882 # Z-boson mass parameter wZ = 2.443 # Z-boson width parameter mW = 80.419 # W-boson mass parameter wW = 2.049 # W-boson width parameter me = 0.000511 # electron mass parameter mmu = 0.1057 # muon mass parameter mtau = 1.777 # tau-lepton mass parameter ms = 0.12 # s-quark mass parameter mc = 1.25 # c-quark mass parameter mb = 4.2 # b-quark mass parameter mtop = 174 # t-quark mass parameter wtop = 1.523 # t-quark width parameter alphas = 0.1178 # Strong coupling constant (Z point) # SUSY: These parameters are not used directly by WHIZARD. # They come into play if an external spectrum generator is invoked. parameter mtype = 1 # MSSM embedding model parameter m_zero = 100 # unified scalar mass (SUGRA/AMSB) parameter m_half = 250 # unified gaugino mass (SUGRA) parameter A0 = -100 # unified trilinear coupling (SUGRA) parameter tanb = 10 # tan(beta) = v2/v1 (MSSM input) parameter sgn_mu = 1 # signum(mu) parameter Lambda = 40000 # SUSY breaking scale (GMSB) parameter M_mes = 80000 # messenger scale (GMSB) parameter N_5 = 3 # messenger index (GMSB) parameter c_grav = 1 # gravitino mass multiplier (GMSB) parameter m_grav = 60000 # gravitino mass (AMSB) # The following parameters are used directly by WHIZARD. # They can be derived from the above, but this is not done by WHIZARD parameter Ae_33 = 0 # Ae soft breaking parameter (3rd gen.) parameter Au_33 = 0 # Au soft breaking parameter (3rd gen.) parameter Ad_33 = 0 # Ad soft breaking parameter (3rd gen.) parameter mh0_1 = 125 # first Higgs mass parameter wh0_1 = 4.143E-3 # first Higgs width parameter mh0_2 = 1000 # second Higgs mass parameter mh0_3 = 1000 # third Higgs mass parameter mA0_1 = 1000 # axial Higgs mass 1 parameter mA0_2 = 1000 # axial Higgs mass 2 parameter mHpm = 1000 # charged Higgs mass parameter wh0_2 = 0 # second Higgs width parameter wh0_3 = 0 # third Higgs width parameter wHpm = 0 # charged Higgs width parameter wA0_1 = 0 # first axial Higgs width parameter wA0_2 = 0 # second axial Higgs width parameter tanb_h = 10 # Higgs mixing angle tan(beta) parameter ls = 1 # Higgs - S coupling parameter a_ls = 0 # ssb Higgs - S coupling parameter ks = 1 # trilinear Singlet coupling parameter a_ks = 0 # ssb trilinear Singlet coupling parameter nmu = 1 # Singlet vev (scaled by ls) parameter mixh0_11 = 1 # scalar Higgs mixing matrix parameter mixh0_12 = 0 # scalar Higgs mixing matrix parameter mixh0_13 = 0 # scalar Higgs mixing matrix parameter mixh0_21 = 0 # scalar Higgs mixing matrix parameter mixh0_22 = 1 # scalar Higgs mixing matrix parameter mixh0_23 = 0 # scalar Higgs mixing matrix parameter mixh0_31 = 0 # scalar Higgs mixing matrix parameter mixh0_32 = 0 # scalar Higgs mixing matrix parameter mixh0_33 = 1 # scalar Higgs mixing matrix parameter mixa0_11 = 1 # axial Higgs mixing matrix parameter mixa0_12 = 0 # axial Higgs mixing matrix parameter mixa0_13 = 0 # axial Higgs mixing matrix parameter mixa0_21 = 1 # axial Higgs mixing matrix parameter mixa0_22 = 0 # axial Higgs mixing matrix parameter mixa0_23 = 0 # axial Higgs mixing matrix parameter msu1 = 1000 # u-squark mass parameter msd1 = 1000 # d-squark mass parameter msc1 = 1000 # c-squark mass parameter mss1 = 1000 # s-squark mass parameter mstop1 = 1000 # t-squark mass parameter msb1 = 1000 # b-squark mass parameter msu2 = 1000 # u-squark mass parameter msd2 = 1000 # d-squark mass parameter msc2 = 1000 # c-squark mass parameter mss2 = 1000 # s-squark mass parameter mstop2 = 1000 # t-squark mass parameter msb2 = 1000 # b-squark mass parameter mse1 = 1000 # selectron1 mass parameter msne = 1000 # electron-sneutrino mass parameter msmu1 = 1000 # smuon1 mass parameter msnmu = 1000 # muon-sneutrino mass parameter mstau1 = 1000 # stau1 mass parameter msntau = 1000 # tau-sneutrino mass parameter mse2 = 1000 # selectron2 mass parameter msmu2 = 1000 # smuon2 mass parameter mstau2 = 1000 # stau2 mass parameter mgg = 1000 # gluino mass parameter mch1 = 1000 # chargino1 mass (signed) parameter mch2 = 1000 # chargino2 mass (signed) parameter mneu1 = 1000 # neutralino1 mass (signed) parameter mneu2 = 1000 # neutralino2 mass (signed) parameter mneu3 = 1000 # neutralino3 mass (signed) parameter mneu4 = 1000 # neutralino4 mass (signed) parameter mneu5 = 1000 # neutralino5 mass (signed) parameter wsu1 = 0 # u-squark width parameter wsd1 = 0 # d-squark width parameter wsc1 = 0 # c-squark width parameter wss1 = 0 # s-squark width parameter wstop1 = 0 # t-squark width parameter wsb1 = 0 # b-squark width parameter wsu2 = 0 # u-squark width parameter wsd2 = 0 # d-squark width parameter wsc2 = 0 # c-squark width parameter wss2 = 0 # s-squark width parameter wstop2 = 0 # t-squark width parameter wsb2 = 0 # b-squark width parameter wse1 = 0 # selectron1 width parameter wsne = 0 # electron-sneutrino width parameter wsmu1 = 0 # smuon1 width parameter wsnmu = 0 # muon-sneutrino width parameter wstau1 = 0 # stau1 width parameter wsntau = 0 # tau-sneutrino width parameter wse2 = 0 # selectron2 width parameter wsmu2 = 0 # smuon2 width parameter wstau2 = 0 # stau2 width parameter wgg = 0 # gluino width parameter wch1 = 0 # chargino1 width parameter wch2 = 0 # chargino2 width parameter wneu1 = 0 # neutralino1 width parameter wneu2 = 0 # neutralino2 width parameter wneu3 = 0 # neutralino3 width parameter wneu4 = 0 # neutralino4 width parameter wneu5 = 0 # neutralino5 width parameter mt_11 = 1 # stop mixing matrix parameter mt_12 = 0 # stop mixing matrix parameter mt_21 = 0 # stop mixing matrix parameter mt_22 = 1 # stop mixing matrix parameter mb_11 = 1 # sbottom mixing matrix parameter mb_12 = 0 # sbottom mixing matrix parameter mb_21 = 0 # sbottom mixing matrix parameter mb_22 = 1 # sbottom mixing matrix parameter ml_11 = 1 # stau mixing matrix parameter ml_12 = 0 # stau mixing matrix parameter ml_21 = 0 # stau mixing matrix parameter ml_22 = 1 # stau mixing matrix parameter mixn_11 = 1 # neutralino mixing matrix parameter mixn_12 = 0 # neutralino mixing matrix parameter mixn_13 = 0 # neutralino mixing matrix parameter mixn_14 = 0 # neutralino mixing matrix parameter mixn_15 = 0 # neutralino mixing matrix parameter mixn_21 = 0 # neutralino mixing matrix parameter mixn_22 = 1 # neutralino mixing matrix parameter mixn_23 = 0 # neutralino mixing matrix parameter mixn_24 = 0 # neutralino mixing matrix parameter mixn_25 = 0 # neutralino mixing matrix parameter mixn_31 = 0 # neutralino mixing matrix parameter mixn_32 = 0 # neutralino mixing matrix parameter mixn_33 = 1 # neutralino mixing matrix parameter mixn_34 = 0 # neutralino mixing matrix parameter mixn_35 = 0 # neutralino mixing matrix parameter mixn_41 = 0 # neutralino mixing matrix parameter mixn_42 = 0 # neutralino mixing matrix parameter mixn_43 = 0 # neutralino mixing matrix parameter mixn_44 = 1 # neutralino mixing matrix parameter mixn_45 = 0 # neutralino mixing matrix parameter mixn_51 = 0 # neutralino mixing matrix parameter mixn_52 = 0 # neutralino mixing matrix parameter mixn_53 = 0 # neutralino mixing matrix parameter mixn_54 = 0 # neutralino mixing matrix parameter mixn_55 = 1 # neutralino mixing matrix parameter mu_11 = 1 # chargino mixing matrix parameter mu_12 = 0 # chargino mixing matrix parameter mu_21 = 0 # chargino mixing matrix parameter mu_22 = 1 # chargino mixing matrix parameter mv_11 = 1 # chargino mixing matrix parameter mv_12 = 0 # chargino mixing matrix parameter mv_21 = 0 # chargino mixing matrix parameter mv_22 = 1 # chargino mixing matrix -parameter vckm11 = 0.97383 # Vud -parameter vckm12 = 0.2272 # Vus -parameter vckm13 = 0.00396 # Vub -parameter vckm21 = -0.2271 # Vcd -parameter vckm22 = 0.97296 # Vcs -parameter vckm23 = 0.04221 # Vcb -parameter vckm31 = 0.00814 # Vtd -parameter vckm32 = -0.04161 # Vts -parameter vckm33 = 0.99910 # Vtb +parameter vckm11 = 0.97401 # Vud +parameter vckm12 = 0.22650 # Vus +parameter vckm13 = 0.00361 # Vub +parameter vckm21 = -0.22636 # Vcd +parameter vckm22 = 0.97320 # Vcs +parameter vckm23 = 0.04053 # Vcb +parameter vckm31 = 0.00854 # Vtd +parameter vckm32 = -0.03978 # Vts +parameter vckm33 = 0.999172 # Vtb # Dependent parameters derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev) derived cw = mW / mZ # cos(theta-W) derived sw = sqrt (1-cw*cw) # sin(theta-W) derived ee = 2 * sw * mW / v # em-coupling (GF scheme) derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const ######################################################################## # Particle content # The quarks particle D_QUARK 1 parton spin 1/2 charge -1/3 isospin -1/2 color 3 name d down anti dbar D "d~" tex_anti "\bar{d}" particle U_QUARK 2 parton spin 1/2 charge 2/3 isospin 1/2 color 3 name u up anti ubar U "u~" tex_anti "\bar{u}" particle S_QUARK 3 like D_QUARK name s strange anti sbar S "s~" tex_anti "\bar{s}" mass ms particle C_QUARK 4 like U_QUARK name c charm anti cbar C "c~" tex_anti "\bar{c}" mass mc particle B_QUARK 5 like D_QUARK name b bottom anti bbar B "b~" tex_anti "\bar{b}" mass mb particle T_QUARK 6 like U_QUARK name t top anti tbar T "t~" tex_anti "\bar{t}" mass mtop width wtop # The leptons particle E_LEPTON 11 spin 1/2 charge -1 isospin -1/2 name "e-" e1 electron e anti "e+" E1 positron tex_name "e^-" tex_anti "e^+" mass me particle E_NEUTRINO 12 left spin 1/2 isospin 1/2 name nue n1 "nu_e" ve "e-neutrino" anti nuebar N1 "ve~" tex_name "\nu_e" tex_anti "\bar\nu_e" particle MU_LEPTON 13 like E_LEPTON name "mu-" e2 mu muon anti "mu+" E2 tex_name "\mu^-" tex_anti "\mu^+" mass mmu particle MU_NEUTRINO 14 like E_NEUTRINO name numu "nu_mu" n2 vm "mu-neutrino" anti numubar N2 "vm~" tex_name "\nu_\mu" tex_anti "\bar\nu_\mu" particle TAU_LEPTON 15 like E_LEPTON name "tau-" e3 tau "ta-" tauon anti "tau+" E3 "ta+" tex_name "\tau^-" tex_anti "\tau^+" mass mtau particle TAU_NEUTRINO 16 like E_NEUTRINO name nutau "nu_tau" n3 vt "tau_neutrino" anti nutaubar N3 "vt~" tex_name "\nu_\tau" tex_anti "\bar\nu_\tau" # The vector bosons particle GLUON 21 parton gauge spin 1 color 8 name gl g G gluon particle PHOTON 22 gauge spin 1 name A gamma photon tex_name "\gamma" particle Z_BOSON 23 gauge spin 1 name Z mass mZ width wZ particle W_BOSON 24 gauge spin 1 charge 1 name "W+" Wp anti "W-" Wm tex_name "W^+" tex_anti "W^-" mass mW width wW # The Higgses particle SCALAR_HIGGS_1 25 spin 0 name h01 h h0 h1 "h_1" Higgs tex_name "H^0_1" mass mh0_1 width wh0_1 particle SCALAR_HIGGS_2 35 spin 0 name h02 H HH HH0 H0 h2 "h_2" tex_name "H^0_2" mass mh0_2 width wh0_2 particle SCALAR_HIGGS_3 45 spin 0 name h03 h3 "h_3" tex_name "H^0_3" mass mh0_3 width wh0_3 particle AXIAL_HIGGS_1 36 spin 0 name A01 HA HA0 A0 A1 "A0_1" tex_name "A^0_1" mass mA0_1 width wA0_1 particle AXIAL_HIGGS_2 46 spin 0 name A02 A2 "A0_2" tex_name "A^0_2" mass mA0_2 width wA0_2 particle CHARGED_HIGGS 37 spin 0 charge +1 name "H+" Hp anti "H-" Hm tex_name "H^+" tex_anti "H^-" mass mHpm width wHpm # The squarks # Left-handed particle D_SQUARK1 1000001 spin 0 charge -1/3 color 3 name sd1 anti sd1c SD1 "sd1~" tex_name "\tilde{d}_L" tex_anti "\tilde{\bar d}_L" mass msd1 width wsd1 particle U_SQUARK1 1000002 spin 0 charge 2/3 color 3 name su1 anti su1c SU1 "su1~" tex_name "\tilde{u}_L" tex_anti "\tilde{\bar u}_L" mass msu1 width wsu1 particle S_SQUARK1 1000003 like D_SQUARK1 name ss1 anti ss1c SS1 "ss1~" tex_name "\tilde{s}_L" tex_anti "\tilde{\bar s}_L" mass mss1 width wss1 particle C_SQUARK1 1000004 like U_SQUARK1 name sc1 anti sc1c SC1 "sc1~" tex_name "\tilde{c}_L" tex_anti "\tilde{\bar c}_L" mass msc1 width wsc1 particle B_SQUARK1 1000005 like D_SQUARK1 name sb1 anti sb1c SB1 "sb1~" tex_name "\tilde{b}_1" tex_anti "\tilde{\bar b}_1" mass msb1 width wsb1 particle T_SQUARK1 1000006 like U_SQUARK1 name st1 anti st1c ST1 "st1~" tex_name "\tilde{t}_1" tex_anti "\tilde{\bar t}_1" mass mstop1 width wstop1 # Right-handed particle D_SQUARK2 2000001 spin 0 charge -1/3 color 3 name sd2 anti sd2c SD2 "sd2~" tex_name "\tilde{d}_R" tex_anti "\tilde{\bar d}_R" mass msd2 width wsd2 particle U_SQUARK2 2000002 spin 0 charge 2/3 color 3 name su2 anti su2c SU2 "su2~" tex_name "\tilde{u}_R" tex_anti "\tilde{\bar u}_R" mass msu2 width wsu2 particle S_SQUARK2 2000003 like D_SQUARK2 name ss2 anti ss2c SS2 "ss2~" tex_name "\tilde{s}_R" tex_anti "\tilde{\bar s}_R" mass mss2 width wss2 particle C_SQUARK2 2000004 like U_SQUARK2 name sc2 anti sc2c SC2 "sc2~" tex_name "\tilde{c}_R" tex_anti "\tilde{\bar c}_R" mass msc2 width wsc2 particle B_SQUARK2 2000005 like D_SQUARK2 name sb2 anti sb2c SB2 "sb2~" tex_name "\tilde{b}_2" tex_anti "\tilde{\bar b}_2" mass msb2 width wsb2 particle T_SQUARK2 2000006 like U_SQUARK2 name st2 anti st2c ST2 "st2~" tex_name "\tilde{t}_2" tex_anti "\tilde{\bar t}_2" mass mstop2 width wstop2 # The sleptons # Left-handed particle E_SLEPTON1 1000011 spin 0 charge -1 name "se1-" se1 se11 anti "se1+" SE1 SE11 tex_name "\tilde{e}_1^-" tex_anti "\tilde{e}_1^+" mass mse1 width wse1 particle E_SNEUTRINO1 1000012 spin 0 name snue "snu_e1" sn11 sve anti "snue*" SN11 "sve~" tex_name "\tilde\nu_e" tex_anti "\tilde\bar\nu_e" mass msne width wsne particle MU_SLEPTON1 1000013 like E_SLEPTON1 name "smu1-" smu1 se21 anti "smu1+" SMU1 SE21 tex_name "\tilde\mu_1^-" tex_anti "\tilde\mu_1^+" mass msmu1 width wsmu1 particle MU_SNEUTRINO1 1000014 like E_SNEUTRINO1 name snumu "snu_mu1" sn21 svm anti "snumu*" SN21 "svm~" tex_name "\tilde\nu_\mu" tex_anti "\tilde\bar\nu_\mu" mass msnmu width wsnmu particle TAU_SLEPTON1 1000015 like E_SLEPTON1 name "stau1-" stau1 se31 anti "stau1+" STAU1 SE31 tex_name "\tilde\tau_1^-" tex_anti "\tilde\tau_1^+" mass mstau1 width wstau1 particle TAU_SNEUTRINO1 1000016 like E_SNEUTRINO1 name snutau "snu_tau" sn31 svt anti "snutau*" SN31 "svt~" tex_name "\tilde\nu_\tau" tex_anti "\tilde\bar\nu_\tau" mass msntau width wsntau # Right-handed particle E_SLEPTON2 2000011 spin 0 charge -1 name "se2-" se2 se12 anti "se2+" SE2 SE12 tex_name "\tilde{e}_2^-" tex_anti "\tilde{e}_2^+" mass mse2 width wse2 particle MU_SLEPTON2 2000013 like E_SLEPTON2 name "smu2-" smu2 se22 anti "smu2+" SMU2 SE22 tex_name "\tilde\mu_2^-" tex_anti "\tilde\mu_2^+" mass msmu2 width wsmu2 particle TAU_SLEPTON2 2000015 like E_SLEPTON2 name "stau2-" stau2 se32 anti "stau2+" STAU2 SE32 tex_name "\tilde\tau_2^-" tex_anti "\tilde\tau_2^+" mass mstau2 width wstau2 # The gauginos particle GLUINO 1000021 spin 1/2 color 8 name sgl gg GG gluino tex_name "\tilde{g}" mass mgg width wgg particle CHARGINO1 1000024 spin 1/2 charge 1 name "ch1+" "CH1+" anti "ch1-" "CH1-" tex_name "\tilde\chi_1^+" tex_anti "\tilde\chi_1^-" mass mch1 width wch1 particle CHARGINO2 1000037 like CHARGINO1 name "ch2+" "CH2+" anti "ch2-" "CH2-" tex_name "\tilde\chi_2^+" tex_anti "\tilde\chi_2^-" mass mch2 width wch2 particle NEUTRALINO1 1000022 spin 1/2 name neu1 NEU1 tex_name "\tilde\chi_1^0" mass mneu1 width wneu1 particle NEUTRALINO2 1000023 like NEUTRALINO1 name neu2 NEU2 tex_name "\tilde\chi_2^0" mass mneu2 width wneu2 particle NEUTRALINO3 1000025 like NEUTRALINO1 name neu3 NEU3 tex_name "\tilde\chi_3^0" mass mneu3 width wneu3 particle NEUTRALINO4 1000035 like NEUTRALINO1 name neu4 NEU4 tex_name "\tilde\chi_4^0" mass mneu4 width wneu4 particle NEUTRALINO5 1000045 like NEUTRALINO1 name neu5 NEU5 tex_name "\tilde\chi_5^0" mass mneu5 width wneu5 # Hadrons particle PROTON 2212 spin 1/2 charge 1 name p "p+" anti pbar "p-" # Beam remnants for proton colliders particle HADRON_REMNANT 90 name hr tex_name "had_r" particle HADRON_REMNANT_SINGLET 91 name hr1 tex_name "had_r^{(1)}" particle HADRON_REMNANT_TRIPLET 92 color 3 name hr3 tex_name "had_r^{(3)}" anti hr3bar tex_anti "had_r^{(\bar 3)}" particle HADRON_REMNANT_OCTET 93 color 8 name hr8 tex_name "had_r^{(8)}" ######################################################################## # Vertices of the NMSSM # In graphs with identical structure, the first vertex is kept for phase space, # therefore, lighter particles come before heavier ones. ! QED vertex D d A vertex U u A vertex S s A vertex C c A vertex B b A vertex T t A vertex E1 e1 A vertex E2 e2 A vertex E3 e3 A ! QED/SUSY vertex sd1 SD1 A vertex su1 SU1 A vertex ss1 SS1 A vertex sc1 SC1 A vertex sb1 SB1 A vertex st1 ST1 A vertex se11 SE11 A vertex se21 SE21 A vertex se31 SE31 A vertex sd2 SD2 A vertex su2 SU2 A vertex ss2 SS2 A vertex sc2 SC2 A vertex sb2 SB2 A vertex st2 ST2 A vertex se12 SE12 A vertex se22 SE22 A vertex se32 SE32 A # QCD vertex G G G vertex G G G G ! SUSY QCD vertex g gg gg ! QCD fermion-gluon vertex D d G vertex U u G vertex S s G vertex C c G vertex B b G vertex T t G ! SUSY sfermion-gluon vertex sd1 SD1 g vertex sd2 SD2 g vertex su1 SU1 g vertex su2 SU2 g vertex ss1 SS1 g vertex ss2 SS2 g vertex sc1 SC1 g vertex sc2 SC2 g vertex sb1 SB1 g vertex sb2 SB2 g vertex st1 ST1 g vertex st2 ST2 g ! SUSY gluino vertex d SD1 gg vertex d SD2 gg vertex u SU1 gg vertex u SU2 gg vertex s SS1 gg vertex s SS2 gg vertex c SC1 gg vertex c SC2 gg vertex b SB1 gg vertex b SB2 gg vertex t ST1 gg vertex t ST2 gg vertex sd1 D gg vertex sd2 D gg vertex su1 U gg vertex su2 U gg vertex ss1 S gg vertex ss2 S gg vertex sc1 C gg vertex sc2 C gg vertex sb1 B gg vertex sb2 B gg vertex st1 T gg vertex st2 T gg # Neutral currents vertex D d Z vertex U u Z vertex S s Z vertex C c Z vertex B b Z vertex T t Z vertex E1 e1 Z vertex E2 e2 Z vertex E3 e3 Z vertex N1 n1 Z vertex N2 n2 Z vertex N3 n3 Z ! Neutral currents/SUSY vertex sd1 SD1 Z vertex su1 SU1 Z vertex ss1 SS1 Z vertex sc1 SC1 Z vertex sb1 SB1 Z vertex st1 ST1 Z vertex se11 SE11 Z vertex se21 SE21 Z vertex se31 SE31 Z vertex sn11 SN11 Z vertex sn21 SN21 Z vertex sn31 SN31 Z vertex sd2 SD2 Z vertex su2 SU2 Z vertex ss2 SS2 Z vertex sc2 SC2 Z vertex sb2 SB2 Z vertex st2 ST2 Z vertex se12 SE12 Z vertex se22 SE22 Z vertex se32 SE32 Z ! 3rd gen mixing vertex sb1 SB2 Z vertex sb2 SB1 Z vertex st1 ST2 Z vertex st2 ST1 Z vertex se31 SE32 Z vertex se32 SE31 Z # Charged currents vertex U d Wp vertex U s Wp vertex U b Wp vertex C d Wp vertex C s Wp vertex C b Wp vertex T d Wp vertex T s Wp vertex T b Wp vertex D u Wm vertex D c Wm vertex D t Wm vertex S u Wm vertex S c Wm vertex S t Wm vertex B u Wm vertex B c Wm vertex B t Wm vertex N1 e1 Wp vertex N2 e2 Wp vertex N3 e3 Wp vertex E1 n1 Wm vertex E2 n2 Wm vertex E3 n3 Wm ! SUSY vertex su1 SD1 Wm vertex su1 SS1 Wm vertex su1 SB1 Wm vertex sc1 SD1 Wm vertex sc1 SS1 Wm vertex sc1 SB1 Wm vertex st1 SD1 Wm vertex st1 SS1 Wm vertex st1 SB1 Wm vertex su2 SD1 Wm vertex su2 SS1 Wm vertex su2 SB1 Wm vertex sc2 SD1 Wm vertex sc2 SS1 Wm vertex sc2 SB1 Wm vertex st2 SD1 Wm vertex st2 SS1 Wm vertex st2 SB1 Wm vertex su1 SD2 Wm vertex su1 SS2 Wm vertex su1 SB2 Wm vertex sc1 SD2 Wm vertex sc1 SS2 Wm vertex sc1 SB2 Wm vertex st1 SD2 Wm vertex st1 SS2 Wm vertex st1 SB2 Wm vertex su2 SD2 Wm vertex su2 SS2 Wm vertex su2 SB2 Wm vertex sc2 SD2 Wm vertex sc2 SS2 Wm vertex sc2 SB2 Wm vertex st2 SD2 Wm vertex st2 SS2 Wm vertex st2 SB2 Wm vertex sn11 SE11 Wm vertex sn21 SE21 Wm vertex sn31 SE31 Wm vertex sd1 SU1 Wp vertex sb1 SU1 Wp vertex sd1 SU1 Wp vertex sd1 SC1 Wp vertex ss1 SC1 Wp vertex sb1 SC1 Wp vertex sd1 ST1 Wp vertex ss1 ST1 Wp vertex sb1 ST1 Wp vertex sd1 SU2 Wp vertex sb1 SU2 Wp vertex sd1 SU2 Wp vertex sd1 SC2 Wp vertex ss1 SC2 Wp vertex sb1 SC2 Wp vertex sd1 ST2 Wp vertex ss1 ST2 Wp vertex sb1 ST2 Wp vertex sd2 SU1 Wp vertex sb2 SU1 Wp vertex sd2 SU1 Wp vertex sd2 SC1 Wp vertex ss2 SC1 Wp vertex sb2 SC1 Wp vertex sd2 ST1 Wp vertex ss2 ST1 Wp vertex sb2 ST1 Wp vertex sd2 SU2 Wp vertex sb2 SU2 Wp vertex sd2 SU2 Wp vertex sd2 SC2 Wp vertex ss2 SC2 Wp vertex sb2 SC2 Wp vertex sd2 ST2 Wp vertex ss2 ST2 Wp vertex sb2 ST2 Wp vertex se11 SN11 Wp vertex se21 SN21 Wp vertex se31 SN31 Wp ! 3rd gen mixing (only sleptons here) vertex sn31 SE32 Wm vertex se32 SN31 Wp # Yukawa (neutral) ### keeping only 3rd generation for the moment # vertex S s h1 # vertex C c h1 vertex B b h1 vertex T t h1 # vertex E2 e2 h1 vertex E3 e3 h1 # vertex S s h2 # vertex C c h2 vertex B b h2 vertex T t h2 # vertex E2 e2 h2 vertex E3 e3 h2 # vertex S s h3 # vertex C c h3 vertex B b h3 vertex T t h3 # vertex E2 e2 h3 vertex E3 e3 h3 # vertex S s A01 # vertex C c A01 vertex B b A01 vertex T t A01 # vertex E2 e2 A01 vertex E3 e3 A01 # vertex S s A02 # vertex C c A02 vertex B b A02 vertex T t A02 # vertex E2 e2 A02 vertex E3 e3 A02 # Yukawa (charged) ! SUSY Higgses vertex e3 N3 Hp vertex E3 n3 Hm vertex b T Hp vertex B t Hm ! SUSY-Yukawa vertex se31 SE31 h1 vertex se31 SE32 h1 vertex se32 SE31 h1 vertex se32 SE32 h1 vertex sb1 SB1 h1 vertex sb1 SB2 h1 vertex sb2 SB1 h1 vertex sb2 SB2 h1 vertex st1 ST1 h1 vertex st1 ST2 h1 vertex st2 ST1 h1 vertex st2 ST2 h1 vertex se31 SE31 h2 vertex se31 SE32 h2 vertex se32 SE31 h2 vertex se32 SE32 h2 vertex sb1 SB1 h2 vertex sb1 SB2 h2 vertex sb2 SB1 h2 vertex sb2 SB2 h2 vertex st1 ST1 h2 vertex st1 ST2 h2 vertex st2 ST1 h2 vertex st2 ST2 h2 vertex se31 SE31 h3 vertex se31 SE32 h3 vertex se32 SE31 h3 vertex se32 SE32 h3 vertex sb1 SB1 h3 vertex sb1 SB2 h3 vertex sb2 SB1 h3 vertex sb2 SB2 h3 vertex st1 ST1 h3 vertex st1 ST2 h3 vertex st2 ST1 h3 vertex st2 ST2 h3 vertex se31 SE31 A01 vertex se31 SE32 A01 vertex se32 SE31 A01 vertex se32 SE32 A01 vertex sb1 SB1 A01 vertex sb1 SB2 A01 vertex sb2 SB1 A01 vertex sb2 SB2 A01 vertex st1 ST1 A01 vertex st1 ST2 A01 vertex st2 ST1 A01 vertex st2 ST2 A01 vertex se31 SE31 A02 vertex se31 SE32 A02 vertex se32 SE31 A02 vertex se32 SE32 A02 vertex sb1 SB1 A02 vertex sb1 SB2 A02 vertex sb2 SB1 A02 vertex sb2 SB2 A02 vertex st1 ST1 A02 vertex st1 ST2 A02 vertex st2 ST1 A02 vertex st2 ST2 A02 vertex se31 SN31 Hp vertex se32 SN31 Hp vertex sb1 ST1 Hp vertex sb1 ST2 Hp vertex sb2 ST1 Hp vertex sb2 ST2 Hp vertex SE31 sn31 Hm vertex SE32 sn31 Hm vertex SB1 st1 Hm vertex SB1 st2 Hm vertex SB2 st1 Hm vertex SB2 st2 Hm # Vector-boson self-interactions vertex Wp Wm A vertex Wp Wm Z vertex Wp Wm Z Z vertex Wp Wp Wm Wm vertex Wp Wm A Z vertex Wp Wm A A # Higgs - vector boson vertex A Hp Hm vertex Z Hp Hm vertex Z h1 A01 vertex Z h2 A01 vertex Z h3 A01 vertex Z h1 A02 vertex Z h2 A02 vertex Z h3 A02 vertex Wp Hm h1 vertex Wm Hp h1 vertex Wp Hm h2 vertex Wm Hp h2 vertex Wp Hm h3 vertex Wm Hp h3 vertex Wp Hm A01 vertex Wm Hp A01 vertex Wp Hm A02 vertex Wm Hp A02 vertex h1 Wp Wm vertex h1 Z Z vertex h2 Wp Wm vertex h2 Z Z vertex h3 Wp Wm vertex h3 Z Z # Higgs self-interactions vertex h1 h1 h1 vertex h1 h1 h2 vertex h1 h1 h3 vertex h1 h2 h1 vertex h1 h2 h2 vertex h1 h2 h3 vertex h1 h3 h1 vertex h1 h3 h2 vertex h1 h3 h3 vertex h2 h1 h1 vertex h2 h1 h2 vertex h2 h1 h3 vertex h2 h2 h1 vertex h2 h2 h2 vertex h2 h2 h3 vertex h2 h3 h1 vertex h2 h3 h2 vertex h2 h3 h3 vertex h3 h1 h1 vertex h3 h1 h2 vertex h3 h1 h3 vertex h3 h2 h1 vertex h3 h2 h2 vertex h3 h2 h3 vertex h3 h3 h1 vertex h3 h3 h2 vertex h3 h3 h3 vertex h1 A01 A01 vertex h1 A01 A02 vertex h1 A02 A02 vertex h2 A01 A01 vertex h2 A01 A02 vertex h2 A02 A02 vertex h3 A01 A01 vertex h3 A01 A02 vertex h3 A02 A02 vertex h1 Hp Hm vertex h2 Hp Hm vertex h3 Hp Hm ! SUSY Higgses vertex A Hp Hm vertex Z Hp Hm vertex Z h1 A01 vertex Z h2 A01 vertex Z h3 A01 vertex Z h1 A02 vertex Z h2 A02 vertex Z h3 A02 vertex Wp Hm h1 vertex Wm Hp h1 vertex Wp Hm h2 vertex Wm Hp h2 vertex Wp Hm h3 vertex Wm Hp h3 vertex Wp Hm A01 vertex Wm Hp A01 vertex Wp Hm A02 vertex Wm Hp A02 vertex h1 Wp Wm vertex h1 Z Z vertex h2 Wp Wm vertex h2 Z Z vertex h3 Wp Wm vertex h3 Z Z ! Charginos vertex "ch1+" "ch1-" A vertex "ch2+" "ch2-" A vertex "ch1+" "ch1-" Z vertex "ch1+" "ch2-" Z vertex "ch2+" "ch1-" Z vertex "ch2+" "ch2-" Z ! Neutralinos vertex neu1 neu1 Z vertex neu1 neu2 Z vertex neu1 neu3 Z vertex neu1 neu4 Z vertex neu1 neu5 Z vertex neu2 neu2 Z vertex neu2 neu3 Z vertex neu2 neu4 Z vertex neu2 neu5 Z vertex neu3 neu3 Z vertex neu3 neu4 Z vertex neu3 neu5 Z vertex neu4 neu4 Z vertex neu4 neu5 Z vertex neu5 neu5 Z ! Charginos + neutralinos vertex "ch1+" neu1 Wm vertex "ch2+" neu1 Wm vertex "ch1+" neu2 Wm vertex "ch2+" neu2 Wm vertex "ch1+" neu3 Wm vertex "ch2+" neu3 Wm vertex "ch1+" neu4 Wm vertex "ch2+" neu4 Wm vertex "ch1+" neu5 Wm vertex "ch2+" neu5 Wm vertex "ch1-" neu1 Wp vertex "ch2-" neu1 Wp vertex "ch1-" neu2 Wp vertex "ch2-" neu2 Wp vertex "ch1-" neu3 Wp vertex "ch2-" neu3 Wp vertex "ch1-" neu4 Wp vertex "ch2-" neu4 Wp vertex "ch1-" neu5 Wp vertex "ch2-" neu5 Wp ! Charginos + Higgs vertex "ch1+" "ch1-" h1 vertex "ch1+" "ch2-" h1 vertex "ch2+" "ch1-" h1 vertex "ch2+" "ch2-" h1 vertex "ch1+" "ch1-" h2 vertex "ch1+" "ch2-" h2 vertex "ch2+" "ch1-" h2 vertex "ch2+" "ch2-" h2 vertex "ch1+" "ch1-" h3 vertex "ch1+" "ch2-" h3 vertex "ch2+" "ch1-" h3 vertex "ch2+" "ch2-" h3 vertex "ch1+" "ch1-" A01 vertex "ch1+" "ch2-" A01 vertex "ch2+" "ch1-" A01 vertex "ch2+" "ch2-" A01 vertex "ch1+" "ch1-" A02 vertex "ch1+" "ch2-" A02 vertex "ch2+" "ch1-" A02 vertex "ch2+" "ch2-" A02 ! Neutralinos + Higgs vertex neu1 neu1 h1 vertex neu1 neu2 h1 vertex neu1 neu3 h1 vertex neu1 neu4 h1 vertex neu1 neu5 h1 vertex neu2 neu2 h1 vertex neu2 neu3 h1 vertex neu2 neu4 h1 vertex neu2 neu5 h1 vertex neu3 neu3 h1 vertex neu3 neu4 h1 vertex neu3 neu5 h1 vertex neu4 neu4 h1 vertex neu4 neu5 h1 vertex neu5 neu5 h1 vertex neu1 neu1 h2 vertex neu1 neu2 h2 vertex neu1 neu3 h2 vertex neu1 neu4 h2 vertex neu1 neu5 h2 vertex neu2 neu2 h2 vertex neu2 neu3 h2 vertex neu2 neu4 h2 vertex neu2 neu5 h2 vertex neu3 neu3 h2 vertex neu3 neu4 h2 vertex neu3 neu5 h2 vertex neu4 neu4 h2 vertex neu4 neu5 h2 vertex neu5 neu5 h2 vertex neu1 neu1 h3 vertex neu1 neu2 h3 vertex neu1 neu3 h3 vertex neu1 neu4 h3 vertex neu1 neu5 h3 vertex neu2 neu2 h3 vertex neu2 neu3 h3 vertex neu2 neu4 h3 vertex neu2 neu5 h3 vertex neu3 neu3 h3 vertex neu3 neu4 h3 vertex neu3 neu5 h3 vertex neu4 neu4 h3 vertex neu4 neu5 h3 vertex neu5 neu5 h3 vertex neu1 neu1 A01 vertex neu1 neu2 A01 vertex neu1 neu3 A01 vertex neu1 neu4 A01 vertex neu1 neu5 A01 vertex neu2 neu2 A01 vertex neu2 neu3 A01 vertex neu2 neu4 A01 vertex neu2 neu5 A01 vertex neu3 neu3 A01 vertex neu3 neu4 A01 vertex neu3 neu5 A01 vertex neu4 neu4 A01 vertex neu4 neu5 A01 vertex neu5 neu5 A01 vertex neu1 neu1 A02 vertex neu1 neu2 A02 vertex neu1 neu3 A02 vertex neu1 neu4 A02 vertex neu1 neu5 A02 vertex neu2 neu2 A02 vertex neu2 neu3 A02 vertex neu2 neu4 A02 vertex neu2 neu5 A02 vertex neu3 neu3 A02 vertex neu3 neu4 A02 vertex neu3 neu5 A02 vertex neu4 neu4 A02 vertex neu4 neu5 A02 vertex neu5 neu5 A02 ! Charginos + neutralinos + Higgs vertex "ch1+" neu1 Hm vertex "ch2+" neu1 Hm vertex "ch1+" neu2 Hm vertex "ch2+" neu2 Hm vertex "ch1+" neu3 Hm vertex "ch2+" neu3 Hm vertex "ch1+" neu4 Hm vertex "ch2+" neu4 Hm vertex "ch1+" neu5 Hm vertex "ch2+" neu5 Hm vertex "ch1-" neu1 Hp vertex "ch2-" neu1 Hp vertex "ch1-" neu1 Hp vertex "ch2-" neu2 Hp vertex "ch1-" neu3 Hp vertex "ch2-" neu3 Hp vertex "ch1-" neu4 Hp vertex "ch2-" neu4 Hp vertex "ch1-" neu5 Hp vertex "ch2-" neu5 Hp ! Lepton-slepton-neutralino vertex e1 SE11 neu1 vertex e1 SE11 neu2 vertex e1 SE11 neu3 vertex e1 SE11 neu4 vertex e1 SE11 neu5 vertex e1 SE12 neu1 vertex e1 SE12 neu2 vertex e1 SE12 neu3 vertex e1 SE12 neu4 vertex e1 SE12 neu5 vertex E1 se11 neu1 vertex E1 se11 neu2 vertex E1 se11 neu3 vertex E1 se11 neu4 vertex E1 se11 neu5 vertex E1 se12 neu1 vertex E1 se12 neu2 vertex E1 se12 neu3 vertex E1 se12 neu4 vertex E1 se12 neu5 vertex e2 SE21 neu1 vertex e2 SE21 neu2 vertex e2 SE21 neu3 vertex e2 SE21 neu4 vertex e2 SE21 neu5 vertex E2 se21 neu1 vertex E2 se21 neu2 vertex E2 se21 neu3 vertex E2 se21 neu4 vertex E2 se21 neu5 vertex e2 SE22 neu1 vertex e2 SE22 neu2 vertex e2 SE22 neu3 vertex e2 SE22 neu4 vertex e2 SE22 neu5 vertex E2 se22 neu1 vertex E2 se22 neu2 vertex E2 se22 neu3 vertex E2 se22 neu4 vertex E2 se22 neu5 vertex e3 SE31 neu1 vertex e3 SE31 neu2 vertex e3 SE31 neu3 vertex e3 SE31 neu4 vertex e3 SE31 neu5 vertex E3 se31 neu1 vertex E3 se31 neu2 vertex E3 se31 neu3 vertex E3 se31 neu4 vertex E3 se31 neu5 vertex e3 SE32 neu1 vertex e3 SE32 neu2 vertex e3 SE32 neu3 vertex e3 SE32 neu4 vertex e3 SE32 neu5 vertex E3 se32 neu1 vertex E3 se32 neu2 vertex E3 se32 neu3 vertex E3 se32 neu4 vertex E3 se32 neu5 ! Neutrino-sneutrino-neutralino vertex n1 SN11 neu1 vertex n1 SN11 neu2 vertex n1 SN11 neu3 vertex n1 SN11 neu4 vertex n1 SN11 neu5 vertex N1 sn11 neu1 vertex N1 sn11 neu2 vertex N1 sn11 neu3 vertex N1 sn11 neu4 vertex N1 sn11 neu5 vertex n2 SN21 neu1 vertex n2 SN21 neu2 vertex n2 SN21 neu3 vertex n2 SN21 neu4 vertex n2 SN21 neu5 vertex N2 sn21 neu1 vertex N2 sn21 neu2 vertex N2 sn21 neu3 vertex N2 sn21 neu4 vertex N2 sn21 neu5 vertex n3 SN31 neu1 vertex n3 SN31 neu2 vertex n3 SN31 neu3 vertex n3 SN31 neu4 vertex n3 SN31 neu5 vertex N3 sn31 neu1 vertex N3 sn31 neu2 vertex N3 sn31 neu3 vertex N3 sn31 neu4 vertex N3 sn31 neu5 ! Quark-squark-neutralino vertex d SD1 neu1 vertex d SD1 neu2 vertex d SD1 neu3 vertex d SD1 neu4 vertex d SD1 neu5 vertex d SD2 neu1 vertex d SD2 neu2 vertex d SD2 neu3 vertex d SD2 neu4 vertex d SD2 neu5 vertex D sd1 neu1 vertex D sd1 neu2 vertex D sd1 neu3 vertex D sd1 neu4 vertex D sd1 neu5 vertex D sd2 neu1 vertex D sd2 neu2 vertex D sd2 neu3 vertex D sd2 neu4 vertex D sd2 neu5 vertex u SU1 neu1 vertex u SU1 neu2 vertex u SU1 neu3 vertex u SU1 neu4 vertex u SU1 neu5 vertex u SU2 neu1 vertex u SU2 neu2 vertex u SU2 neu3 vertex u SU2 neu4 vertex u SU2 neu5 vertex U su1 neu1 vertex U su1 neu2 vertex U su1 neu3 vertex U su1 neu4 vertex U su1 neu5 vertex U su2 neu1 vertex U su2 neu2 vertex U su2 neu3 vertex U su2 neu4 vertex U su2 neu5 vertex s SS1 neu1 vertex s SS1 neu2 vertex s SS1 neu3 vertex s SS1 neu4 vertex s SS1 neu5 vertex s SS2 neu1 vertex s SS2 neu2 vertex s SS2 neu3 vertex s SS2 neu4 vertex s SS2 neu5 vertex S ss1 neu1 vertex S ss1 neu2 vertex S ss1 neu3 vertex S ss1 neu4 vertex S ss1 neu5 vertex S ss2 neu1 vertex S ss2 neu2 vertex S ss2 neu3 vertex S ss2 neu4 vertex S ss2 neu5 vertex c SC1 neu1 vertex c SC1 neu2 vertex c SC1 neu3 vertex c SC1 neu4 vertex c SC1 neu5 vertex c SC2 neu1 vertex c SC2 neu2 vertex c SC2 neu3 vertex c SC2 neu4 vertex c SC2 neu5 vertex C sc1 neu1 vertex C sc1 neu2 vertex C sc1 neu3 vertex C sc1 neu4 vertex C sc1 neu5 vertex C sc2 neu1 vertex C sc2 neu2 vertex C sc2 neu3 vertex C sc2 neu4 vertex C sc2 neu5 vertex b SB1 neu1 vertex b SB1 neu2 vertex b SB1 neu3 vertex b SB1 neu4 vertex b SB1 neu5 vertex b SB2 neu1 vertex b SB2 neu2 vertex b SB2 neu3 vertex b SB2 neu4 vertex b SB2 neu5 vertex B sb1 neu1 vertex B sb1 neu2 vertex B sb1 neu3 vertex B sb1 neu4 vertex B sb1 neu5 vertex B sb2 neu1 vertex B sb2 neu2 vertex B sb2 neu3 vertex B sb2 neu4 vertex B sb2 neu5 vertex t ST1 neu1 vertex t ST1 neu2 vertex t ST1 neu3 vertex t ST1 neu4 vertex t ST1 neu5 vertex t ST2 neu1 vertex t ST2 neu2 vertex t ST2 neu3 vertex t ST2 neu4 vertex t ST2 neu5 vertex T st1 neu1 vertex T st1 neu2 vertex T st1 neu3 vertex T st1 neu4 vertex T st1 neu5 vertex T st2 neu1 vertex T st2 neu2 vertex T st2 neu3 vertex T st2 neu4 vertex T st2 neu5 ! Lepton-sneutrino-chargino vertex e1 SN11 "ch1+" vertex e1 SN11 "ch2+" vertex E1 sn11 "ch1-" vertex E1 sn11 "ch2-" vertex e2 SN21 "ch1+" vertex e2 SN21 "ch2+" vertex E2 sn21 "ch1-" vertex E2 sn21 "ch2-" vertex e3 SN31 "ch1+" vertex e3 SN31 "ch2+" vertex E3 sn31 "ch1-" vertex E3 sn31 "ch2-" ! Slepton-neutrino-chargino vertex se11 N1 "ch1+" vertex se11 N1 "ch2+" vertex se12 N1 "ch1+" vertex se12 N1 "ch2+" vertex SE11 n1 "ch1-" vertex SE11 n1 "ch2-" vertex SE12 n1 "ch1-" vertex SE12 n1 "ch2-" vertex se21 N2 "ch1+" vertex se21 N2 "ch2+" vertex se22 N2 "ch1+" vertex se22 N2 "ch2+" vertex SE21 n2 "ch1-" vertex SE21 n2 "ch2-" vertex SE22 n2 "ch1-" vertex SE22 n2 "ch2-" vertex se31 N3 "ch1+" vertex se31 N3 "ch2+" vertex se32 N3 "ch1+" vertex se32 N3 "ch2+" vertex SE31 n3 "ch1-" vertex SE31 n3 "ch2-" vertex SE32 n3 "ch1-" vertex SE32 n3 "ch2-" ! Quark-squark-chargino vertex d SU1 "ch1+" vertex d SU1 "ch2+" vertex d SU2 "ch1+" vertex d SU2 "ch2+" vertex s SU1 "ch1+" vertex s SU1 "ch2+" vertex s SU2 "ch1+" vertex s SU2 "ch2+" vertex b SU1 "ch1+" vertex b SU1 "ch2+" vertex b SU2 "ch1+" vertex b SU2 "ch2+" vertex D su1 "ch1-" vertex D su1 "ch2-" vertex D su2 "ch1-" vertex D su2 "ch2-" vertex S su1 "ch1-" vertex S su1 "ch2-" vertex S su2 "ch1-" vertex S su2 "ch2-" vertex B su1 "ch1-" vertex B su1 "ch2-" vertex B su2 "ch1-" vertex B su2 "ch2-" vertex sd1 U "ch1+" vertex sd1 U "ch2+" vertex sd2 U "ch1+" vertex sd2 U "ch2+" vertex ss1 U "ch1+" vertex ss1 U "ch2+" vertex ss2 U "ch1+" vertex ss2 U "ch2+" vertex sb1 U "ch1+" vertex sb1 U "ch2+" vertex sb2 U "ch1+" vertex sb2 U "ch2+" vertex SD1 u "ch1-" vertex SD1 u "ch2-" vertex SD2 u "ch1-" vertex SD2 u "ch2-" vertex SS1 u "ch1-" vertex SS1 u "ch2-" vertex SS2 u "ch1-" vertex SS2 u "ch2-" vertex SB1 u "ch1-" vertex SB1 u "ch2-" vertex SB2 u "ch1-" vertex SB2 u "ch2-" vertex d SC1 "ch1+" vertex d SC1 "ch2+" vertex d SC2 "ch1+" vertex d SC2 "ch2+" vertex s SC1 "ch1+" vertex s SC1 "ch2+" vertex s SC2 "ch1+" vertex s SC2 "ch2+" vertex b SC1 "ch1+" vertex b SC1 "ch2+" vertex b SC2 "ch1+" vertex b SC2 "ch2+" vertex D sc1 "ch1-" vertex D sc1 "ch2-" vertex D sc2 "ch1-" vertex D sc2 "ch2-" vertex S sc1 "ch1-" vertex S sc1 "ch2-" vertex S sc2 "ch1-" vertex S sc2 "ch2-" vertex B sc1 "ch1-" vertex B sc1 "ch2-" vertex B sc2 "ch1-" vertex B sc2 "ch2-" vertex sd1 C "ch1+" vertex sd1 C "ch2+" vertex sd2 C "ch1+" vertex sd2 C "ch2+" vertex ss1 C "ch1+" vertex ss1 C "ch2+" vertex ss2 C "ch1+" vertex ss2 C "ch2+" vertex sb1 C "ch1+" vertex sb1 C "ch2+" vertex sb2 C "ch1+" vertex sb2 C "ch2+" vertex SD1 c "ch1-" vertex SD1 c "ch2-" vertex SD2 c "ch1-" vertex SD2 c "ch2-" vertex SS1 c "ch1-" vertex SS1 c "ch2-" vertex SS2 c "ch1-" vertex SS2 c "ch2-" vertex SB1 c "ch1-" vertex SB1 c "ch2-" vertex SB2 c "ch1-" vertex SB2 c "ch2-" vertex d ST1 "ch1+" vertex d ST1 "ch2+" vertex d ST2 "ch1+" vertex d ST2 "ch2+" vertex s ST1 "ch1+" vertex s ST1 "ch2+" vertex s ST2 "ch1+" vertex s ST2 "ch2+" vertex b ST1 "ch1+" vertex b ST1 "ch2+" vertex b ST2 "ch1+" vertex b ST2 "ch2+" vertex D st1 "ch1-" vertex D st1 "ch2-" vertex D st2 "ch1-" vertex D st2 "ch2-" vertex S st1 "ch1-" vertex S st1 "ch2-" vertex S st2 "ch1-" vertex S st2 "ch2-" vertex B st1 "ch1-" vertex B st1 "ch2-" vertex B st2 "ch1-" vertex B st2 "ch2-" vertex sd1 T "ch1+" vertex sd1 T "ch2+" vertex sd2 T "ch1+" vertex sd2 T "ch2+" vertex ss1 T "ch1+" vertex ss1 T "ch2+" vertex ss2 T "ch1+" vertex ss2 T "ch2+" vertex sb1 T "ch1+" vertex sb1 T "ch2+" vertex sb2 T "ch1+" vertex sb2 T "ch2+" vertex SD1 t "ch1-" vertex SD1 t "ch2-" vertex SD2 t "ch1-" vertex SD2 t "ch2-" vertex SS1 t "ch1-" vertex SS1 t "ch2-" vertex SS2 t "ch1-" vertex SS2 t "ch2-" vertex SB1 t "ch1-" vertex SB1 t "ch2-" vertex SB2 t "ch1-" vertex SB2 t "ch2-"