generate events with HC NLO model

Asked by Yue Xu on 2020-08-06

Dear MadGraph experts,

I'm using MG5 2.6.6 to generate events with HC NLO model(here: https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_NLO_X0_UFO). The precess is:
import model HC_NLO_X0_UFO-VVV
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define v = z w- w+
generate p p > x0 > v v v / a [QCD]

When ran with ./bin/generate_events, I got the error:
Error detected in "launch "
write debug file /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/run_02_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : An error occurred during the collection of results.
 Please check the .log files inside the directories which failed:
 /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zzz_no_a/GF1/log.txt
 /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zwpwm_no_a/GF1/log.txt
 /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zzz_no_a/GF1/log.txt
 /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zwpwm_no_a/GF1/log.txt

It seems the error is caused by add the b-quark in p definition. But when I use HC UFO model to do leading-order generation and add b-quark in p and j definition, it can work.

It can't work in next-to-leading order case. Do you know why ?

the run_02_tag_1_debug.log is:
launch
Traceback (most recent call last):
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 1679, in do_launch
    evt_file = self.run(mode, options)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 1913, in run
    jobs_to_collect,mint_step,mode,mode_dict[mode],fixed_order=False)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 2128, in collect_the_results
    self.append_the_results(jobs_to_run,integration_step)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 2568, in append_the_results
    '\n'.join(error_log)+'\n')
aMCatNLOError: An error occurred during the collection of results.
Please check the .log files inside the directories which failed:
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zzz_no_a/GF1/log.txt
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zwpwm_no_a/GF1/log.txt
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zzz_no_a/GF1/log.txt
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zwpwm_no_a/GF1/log.txt

Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
                  timeout : 60
       cluster_local_path : None
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lhapdf6/bin/lhapdf-config
             cluster_size : 100
           cluster_memory : None
                    pjfry : None
    cluster_status_update : (600, 30)
             cluster_time : None
            f2py_compiler : None
                    ninja : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
               hepmc_path : None
mg5amc_py8_interface_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/MG5aMC_PY8_interface
             pythia8_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/pythia8
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : None
          pythia-pgs_path : None
                 mg5_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : None
         fortran_compiler : None
                  nb_core : 48
                  collier : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                    golem : None
             cpp_compiler : None
      notification_center : True
                 run_mode : 2

the /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zzz_no_a/GF1/log.txt is here:
 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11799999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************
 External Params
  ---------------------------------

 MU_R = 91.188000000000002
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663700000000000E-005
 aS = 0.11799999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 173.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MB = 4.7000000000000002
 mdl_MZ = 91.188000000000002
 mdl_MTA = 1.7769999999999999
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_MX0 = 300.00000000000000
 mdl_WX0 = 1.3131800000000001E-002
 mdl_Lambda = 5000.0000000000000
 mdl_cosa = 1.0000000000000000
 mdl_kSM = 5.0000000000000003E-002
 mdl_kHzz = -0.31488600000000000
 mdl_kHww = -0.31488600000000000
 mdl_kHdz = -0.15744300000000000
 mdl_kHdwR = -0.15744300000000000
  Internal Params
  ---------------------------------

 mdl_CKM11 = 1.0000000000000000
 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_kHdwI = 0.0000000000000000
 mdl_conjg__CKM11 = 1.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_MZ__exp__2 = 8315.2513440000002
 mdl_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
mdl_MX0__exp__2 = 90000.000000000000
 mdl_cosa__exp__2 = 1.0000000000000000
 mdl_sina = 0.0000000000000000
 mdl_kHdw = (-0.15744300000000000 , 0.0000000000000000 )
 mdl_nb__2__exp__0_75 = 1.6817928305074290
 mdl_Ncol__exp__2 = 9.0000000000000000
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29929.000000000000
 mdl_conjg__kHdw = (-0.15744300000000000 , -0.0000000000000000 )
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824863542806241
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6372.0088396276369
 mdl_sw2 = 0.23369618355247201
 mdl_cw = 0.87538780917232795
 mdl_sqrt__sw2 = 0.48342133129649134
 mdl_sw = 0.48342133129649134
 mdl_g1 = 0.35807101351559184
 mdl_gw = 0.64840125943322025
 mdl_vev = 246.22056907348596
 mdl_vev__exp__2 = 60624.568634871270
 mdl_lam = 0.74227332273529756
 mdl_yb = 2.6995322804122715E-002
 mdl_yt = 0.99365762661983603
 mdl_ytau = 1.0206529494239587E-002
 mdl_muH = 212.13203435596427
 mdl_AxialZUp = -0.18517543100305406
 mdl_AxialZDown = 0.18517543100305406
 mdl_VectorZUp = 6.9775986968126480E-002
 mdl_VectorZDown = -0.12747570898559027
 mdl_VectorAUp = 0.20896733336635265
 mdl_VectorADown = -0.10448366668317632
 mdl_VectorWmDxU = 0.22924446373756391
 mdl_AxialWmDxU = -0.22924446373756391
 mdl_VectorWpUxD = 0.22924446373756391
 mdl_AxialWpUxD = -0.22924446373756391
 mdl_I1x33 = ( 2.6995322804122715E-002, 0.0000000000000000 )
 mdl_I2x33 = ( 0.99365762661983603 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99365762661983603 , 0.0000000000000000 )
 mdl_I4x33 = ( 2.6995322804122715E-002, 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.96666230381571328 , 0.0000000000000000 )
 mdl_Axial_tbGp = ( -1.0206529494239587 , -0.0000000000000000 )
mdl_Vector_tbGm = ( 0.96666230381571328 , 0.0000000000000000 )
 mdl_Axial_tbGm = ( -1.0206529494239587 , -0.0000000000000000 )
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_gAaa = 1.3477023594718655E-005
 mdl_cw__exp__2 = 0.76630381644752810
 mdl_gAza = 4.5000939025507963E-006
 mdl_gHaa = 2.6392504539657365E-005
 mdl_gHza = 3.9166796815280137E-005
 mdl_gw__exp__2 = 0.42042419323458619
 mdl_sw__exp__2 = 0.23369618355247204
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_G__exp__4 = 2.1987899468922913
 mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
 mdl_GWcft_UV_b_1EPS_ = -3.1300472141406080E-003
 mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
 mdl_bWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_G__exp__3 = 1.8056676068262196
 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_GWcft_UV_b_FIN_ = -1.8563438626678915E-002
 mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
 mdl_bWcft_UV_FIN_ = -0.13642100947319838
 mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
 mdl_gAgg = 7.6274225811079832E-005
 mdl_gHgg = -5.0849483874053226E-005
  Couplings of HC_NLO_X0_UFO-VVV
  ---------------------------------
 UVWfct_b_0 = -0.13642E+00 -0.00000E+00
    R2_UUZ_V2 = -0.00000E+00 -0.71065E-02
    R2_UUZ_V5 = 0.00000E+00 0.72241E-03
      R2_dxuW = -0.00000E+00 -0.11481E-01
         GC_5 = 0.00000E+00 0.12177E+01
 UVWfct_b_0_1eps -0.18780E-01 0.00000E+00
    GC_3005h2 = 0.00000E+00 0.31489E-04
    GC_3005h3 = 0.00000E+00 0.62977E-04
    GC_3007h2 = 0.00000E+00 0.62977E-04
    GC_3007h3 = 0.00000E+00 0.31489E-04
        GC_22 = 0.00000E+00 0.28380E+00
        GC_23 = -0.00000E+00 -0.28850E-01
        GC_47 = 0.00000E+00 0.45849E+00
    GC_3005h1 = 0.00000E+00 0.25879E+01
    GC_3007h1 = 0.00000E+00 0.33772E+01

 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23nlo
 alpha_s(Mz)= 0.1190 running at 2 loops.
 alpha_s(Mz)= 0.1190 running at 2 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis
Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.9999999999999999E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 1
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the all of this channel
 Normal integration (Sfunction != 1)
 about to integrate 10 -1 12 1
 imode is 0
channel 1 : 1 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
 ------- iteration 1
 Update # PS points (even): 800 --> 800
Using random seed offsets: 1 , 17 , 0
  with seed 34
 Ranmar initialization seeds 13168 9425
 Total number of FKS directories is 4
 FKS process map (sum= 3 ) :
           1 --> 2 : 1 3
           2 --> 2 : 2 4
 ================================
 process combination map (specified per FKS dir):
  1 map 1
  1 inv. map 1
  2 map 1
  2 inv. map 1
  3 map 1
  3 inv. map 1
  4 map 1
  4 inv. map 1
 ================================
tau_min 1 1 : 0.27356E+03 -- 0.39119E+03
tau_min 2 1 : 0.27356E+03 -- 0.39119E+03
tau_min 3 1 : 0.27356E+03 -- 0.39119E+03
tau_min 4 1 : 0.27356E+03 -- 0.39119E+03
 Mass shell violation [nocms]
j= 1
 mass= 4.7000000000000002
 mass computed= 0.0000000000000000
0.42197153D+03 0.00000000D+00 0.00000000D+00 0.42197153D+03
 Mass shell violation [nocms]
 j= 1
 mass= 4.7000000000000002
 mass computed= 0.0000000000000000
0.42148321D+03 0.00000000D+00 0.00000000D+00 0.42148321D+03
 bpower is 0.0000000000000000
 Mass shell violation [nocms]
 j= 1
 mass= 4.7000000000000002
 mass computed= 0.0000000000000000
0.42197153D+03 0.00000000D+00 0.00000000D+00 0.42197153D+03
 Mass shell violation [nocms]
 j= 1
 mass= 4.7000000000000002
 mass computed= 0.0000000000000000
0.42148321D+03 0.00000000D+00 0.00000000D+00 0.42148321D+03
 Fatal error in set_tau_min
Time in seconds: 0

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Yue Xu
Solved:
2020-08-06
Last query:
2020-08-06
Last reply:
2020-08-06

Hi,

The issue is the presence of the top-quark resonances at the real diagram contribution which breaks NLO accuracy. You have already the same issue at LO when you try to do matched/merged sample (with either MLM or CKKW(-L))

This issue is not present with 4 flavor computation since you do not face onshell top contribution.

If you want to do such computation with 5 flavor you need to choose a scheme of substraction of the resonances. Such substraction scheme is not part of the default aMC@NLO but are available via the use of a plugin (MadSTR).

Cheers,

Olivier

> On 6 Aug 2020, at 15:55, Yue Xu <email address hidden> wrote:
>
> New question #692244 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/692244
>
> Dear MadGraph experts,
>
> I'm using MG5 2.6.6 to generate events with HC NLO model(here: https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_NLO_X0_UFO). The precess is:
> import model HC_NLO_X0_UFO-VVV
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = g u c d s b u~ c~ d~ s~ b~
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> define v = z w- w+
> generate p p > x0 > v v v / a [QCD]
>
> When ran with ./bin/generate_events, I got the error:
> Error detected in "launch "
> write debug file /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/run_02_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
> aMCatNLOError : An error occurred during the collection of results.
> Please check the .log files inside the directories which failed:
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zzz_no_a/GF1/log.txt
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zwpwm_no_a/GF1/log.txt
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zzz_no_a/GF1/log.txt
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zwpwm_no_a/GF1/log.txt
>
> It seems the error is caused by add the b-quark in p definition. But when I use HC UFO model to do leading-order generation and add b-quark in p and j definition, it can work.
>
> It can't work in next-to-leading order case. Do you know why ?
>
> the run_02_tag_1_debug.log is:
> launch
> Traceback (most recent call last):
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/extended_cmd.py", line 1514, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
> return func(arg, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 1679, in do_launch
> evt_file = self.run(mode, options)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 1913, in run
> jobs_to_collect,mint_step,mode,mode_dict[mode],fixed_order=False)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 2128, in collect_the_results
> self.append_the_results(jobs_to_run,integration_step)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/bin/internal/amcatnlo_run_interface.py", line 2568, in append_the_results
> '\n'.join(error_log)+'\n')
> aMCatNLOError: An error occurred during the collection of results.
> Please check the .log files inside the directories which failed:
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zzz_no_a/GF1/log.txt
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zwpwm_no_a/GF1/log.txt
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zzz_no_a/GF1/log.txt
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bxb_x0_zwpwm_no_a/GF1/log.txt
>
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_local_path : None
> cluster_queue : None
> madanalysis_path : None
> lhapdf : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lhapdf6/bin/lhapdf-config
> cluster_size : 100
> cluster_memory : None
> pjfry : None
> cluster_status_update : (600, 30)
> cluster_time : None
> f2py_compiler : None
> ninja : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
> hepmc_path : None
> mg5amc_py8_interface_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/MG5aMC_PY8_interface
> pythia8_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/pythia8
> hwpp_path : None
> automatic_html_opening : False
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> mg5_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6
> td_path : None
> delphes_path : None
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : None
> fortran_compiler : None
> nb_core : 48
> collier : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> notification_center : True
> run_mode : 2
>
> the /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/WorkArea/tmp300_vvv_withb_error/HC_NLO_LHC_VHH_VV/SubProcesses/P0_bbx_x0_zzz_no_a/GF1/log.txt is here:
> ===============================================================
> INFO: MadFKS read these parameters from FKS_params.dat
> ===============================================================
>> IRPoleCheckThreshold = 1.0000000000000001E-005
>> PrecisionVirtualAtRunTime = 1.0000000000000000E-003
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000001E-003
> ===============================================================
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> ****************************************
>
> NNPDFDriver version 1.0.3
> Grid: NNPDF23nlo_as_0119_qed_mem0.grid
> ****************************************
> New value of alpha_s from PDF nn23nlo: 0.11899999999999999
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
> External Params
> ---------------------------------
>
> MU_R = 91.188000000000002
> aEWM1 = 127.90000000000001
> mdl_Gf = 1.1663700000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MZ = 91.188000000000002
> mdl_MTA = 1.7769999999999999
> mdl_WT = 1.4915000000000000
> mdl_WZ = 2.4952000000000001
> mdl_WW = 2.0850000000000000
> mdl_MX0 = 300.00000000000000
> mdl_WX0 = 1.3131800000000001E-002
> mdl_Lambda = 5000.0000000000000
> mdl_cosa = 1.0000000000000000
> mdl_kSM = 5.0000000000000003E-002
> mdl_kHzz = -0.31488600000000000
> mdl_kHww = -0.31488600000000000
> mdl_kHdz = -0.15744300000000000
> mdl_kHdwR = -0.15744300000000000
> Internal Params
> ---------------------------------
>
> mdl_CKM11 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_kHdwI = 0.0000000000000000
> mdl_conjg__CKM11 = 1.0000000000000000
> mdl_lhv = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM33 = 1.0000000000000000
> mdl_Ncol = 3.0000000000000000
> mdl_CA = 3.0000000000000000
> mdl_TF = 0.50000000000000000
> mdl_CF = 1.3333333333333333
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MX0__exp__2 = 90000.000000000000
> mdl_cosa__exp__2 = 1.0000000000000000
> mdl_sina = 0.0000000000000000
> mdl_kHdw = (-0.15744300000000000 , 0.0000000000000000 )
> mdl_nb__2__exp__0_75 = 1.6817928305074290
> mdl_Ncol__exp__2 = 9.0000000000000000
> mdl_MB__exp__2 = 22.090000000000003
> mdl_MT__exp__2 = 29929.000000000000
> mdl_conjg__kHdw = (-0.15744300000000000 , -0.0000000000000000 )
> mdl_aEW = 7.8186082877247844E-003
> mdl_MW = 79.824863542806241
> mdl_sqrt__aEW = 8.8422894590285753E-002
> mdl_ee = 0.31345100004952897
> mdl_MW__exp__2 = 6372.0088396276369
> mdl_sw2 = 0.23369618355247201
> mdl_cw = 0.87538780917232795
> mdl_sqrt__sw2 = 0.48342133129649134
> mdl_sw = 0.48342133129649134
> mdl_g1 = 0.35807101351559184
> mdl_gw = 0.64840125943322025
> mdl_vev = 246.22056907348596
> mdl_vev__exp__2 = 60624.568634871270
> mdl_lam = 0.74227332273529756
> mdl_yb = 2.6995322804122715E-002
> mdl_yt = 0.99365762661983603
> mdl_ytau = 1.0206529494239587E-002
> mdl_muH = 212.13203435596427
> mdl_AxialZUp = -0.18517543100305406
> mdl_AxialZDown = 0.18517543100305406
> mdl_VectorZUp = 6.9775986968126480E-002
> mdl_VectorZDown = -0.12747570898559027
> mdl_VectorAUp = 0.20896733336635265
> mdl_VectorADown = -0.10448366668317632
> mdl_VectorWmDxU = 0.22924446373756391
> mdl_AxialWmDxU = -0.22924446373756391
> mdl_VectorWpUxD = 0.22924446373756391
> mdl_AxialWpUxD = -0.22924446373756391
> mdl_I1x33 = ( 2.6995322804122715E-002, 0.0000000000000000 )
> mdl_I2x33 = ( 0.99365762661983603 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99365762661983603 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995322804122715E-002, 0.0000000000000000 )
> mdl_Vector_tbGp = (-0.96666230381571328 , 0.0000000000000000 )
> mdl_Axial_tbGp = ( -1.0206529494239587 , -0.0000000000000000 )
> mdl_Vector_tbGm = ( 0.96666230381571328 , 0.0000000000000000 )
> mdl_Axial_tbGm = ( -1.0206529494239587 , -0.0000000000000000 )
> mdl_ee__exp__2 = 9.8251529432049817E-002
> mdl_gAaa = 1.3477023594718655E-005
> mdl_cw__exp__2 = 0.76630381644752810
> mdl_gAza = 4.5000939025507963E-006
> mdl_gHaa = 2.6392504539657365E-005
> mdl_gHza = 3.9166796815280137E-005
> mdl_gw__exp__2 = 0.42042419323458619
> mdl_sw__exp__2 = 0.23369618355247204
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> mdl_G__exp__4 = 2.1987899468922913
> mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
> mdl_GWcft_UV_b_1EPS_ = -3.1300472141406080E-003
> mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
> mdl_bWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_G__exp__3 = 1.8056676068262196
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_GWcft_UV_b_FIN_ = -1.8563438626678915E-002
> mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
> mdl_bWcft_UV_FIN_ = -0.13642100947319838
> mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
> mdl_gAgg = 7.6274225811079832E-005
> mdl_gHgg = -5.0849483874053226E-005
> Couplings of HC_NLO_X0_UFO-VVV
> ---------------------------------
> UVWfct_b_0 = -0.13642E+00 -0.00000E+00
> R2_UUZ_V2 = -0.00000E+00 -0.71065E-02
> R2_UUZ_V5 = 0.00000E+00 0.72241E-03
> R2_dxuW = -0.00000E+00 -0.11481E-01
> GC_5 = 0.00000E+00 0.12177E+01
> UVWfct_b_0_1eps -0.18780E-01 0.00000E+00
> GC_3005h2 = 0.00000E+00 0.31489E-04
> GC_3005h3 = 0.00000E+00 0.62977E-04
> GC_3007h2 = 0.00000E+00 0.62977E-04
> GC_3007h3 = 0.00000E+00 0.31489E-04
> GC_22 = 0.00000E+00 0.28380E+00
> GC_23 = -0.00000E+00 -0.28850E-01
> GC_47 = 0.00000E+00 0.45849E+00
> GC_3005h1 = 0.00000E+00 0.25879E+01
> GC_3007h1 = 0.00000E+00 0.33772E+01
>
> Collider parameters:
> --------------------
>
> Running at P P machine @ 13000.000000000000 GeV
> PDF set = nn23nlo
> alpha_s(Mz)= 0.1190 running at 2 loops.
> alpha_s(Mz)= 0.1190 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
> Diagram information for clustering has been set-up for nFKSprocess 1
> Diagram information for clustering has been set-up for nFKSprocess 2
> Diagram information for clustering has been set-up for nFKSprocess 3
> Diagram information for clustering has been set-up for nFKSprocess 4
> getting user params
> Enter number of events and iterations:
> Number of events and iterations -1 12
> Enter desired fractional accuracy:
> Desired fractional accuracy: 2.9999999999999999E-002
> Enter alpha, beta for G_soft
> Enter alpha<0 to set G_soft=1 (no ME soft)
> for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
> Enter alpha, beta for G_azi
> Enter alpha>0 to set G_azi=0 (no azi corr)
> for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
> Doing the S and H events together
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Do MC over helicities for the virtuals
> Enter Configuration Number:
> Running Configuration Number: 1
> Enter running mode for MINT:
> 0 to set-up grids, 1 to integrate, 2 to generate events
> MINT running mode: 0
> Set the three folding parameters for MINT
> xi_i, phi_i, y_ij
> 1 1 1
> 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
> Enter 'born0' or 'virt0' to perform
> a pure n-body integration (no S functions)
> doing the all of this channel
> Normal integration (Sfunction != 1)
> about to integrate 10 -1 12 1
> imode is 0
> channel 1 : 1 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
> ------- iteration 1
> Update # PS points (even): 800 --> 800
> Using random seed offsets: 1 , 17 , 0
> with seed 34
> Ranmar initialization seeds 13168 9425
> Total number of FKS directories is 4
> FKS process map (sum= 3 ) :
> 1 --> 2 : 1 3
> 2 --> 2 : 2 4
> ================================
> process combination map (specified per FKS dir):
> 1 map 1
> 1 inv. map 1
> 2 map 1
> 2 inv. map 1
> 3 map 1
> 3 inv. map 1
> 4 map 1
> 4 inv. map 1
> ================================
> tau_min 1 1 : 0.27356E+03 -- 0.39119E+03
> tau_min 2 1 : 0.27356E+03 -- 0.39119E+03
> tau_min 3 1 : 0.27356E+03 -- 0.39119E+03
> tau_min 4 1 : 0.27356E+03 -- 0.39119E+03
> Mass shell violation [nocms]
> j= 1
> mass= 4.7000000000000002
> mass computed= 0.0000000000000000
> 0.42197153D+03 0.00000000D+00 0.00000000D+00 0.42197153D+03
> Mass shell violation [nocms]
> j= 1
> mass= 4.7000000000000002
> mass computed= 0.0000000000000000
> 0.42148321D+03 0.00000000D+00 0.00000000D+00 0.42148321D+03
> bpower is 0.0000000000000000
> Mass shell violation [nocms]
> j= 1
> mass= 4.7000000000000002
> mass computed= 0.0000000000000000
> 0.42197153D+03 0.00000000D+00 0.00000000D+00 0.42197153D+03
> Mass shell violation [nocms]
> j= 1
> mass= 4.7000000000000002
> mass computed= 0.0000000000000000
> 0.42148321D+03 0.00000000D+00 0.00000000D+00 0.42148321D+03
> Fatal error in set_tau_min
> Time in seconds: 0
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #2

Hi Olivier,

Thanks for the quick action.

We may try to use MadSTR or just not add b-quark in p definition.

Cheers,
Yue