generate events with a modified HC model

Asked by Yue Xu on 2019-11-04

Hello,

I want to generate some events with a modified HC model. This the process is(where x0 is a heavy higgs):
import model HC_UFO
generate p p > z x0, z > l- l+, x0 > j j j j
output ./TEDWORK/HC_lljjjj

The model is at:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_UFO

The process and run_01_tag_1_debug.log is at:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_lljjjj

When I generate events with "./bin/madevent", I get some errors:
Error when reading /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/TEDWORK/HC_lljjjj/SubProcesses/P1_qq_zx0_z_ll_x0_ggqq/G6/results.dat
Command "launch " interrupted with error:
Exception : Reported error: End code 4.0
  Full associated log:
  Process in group number 1
  A PDF is used, so alpha_s(MZ) is going to be modified
  Old value of alpha_s from param_card: 0.11839999999999999
   ****************************************

        NNPDFDriver version 1.0.3
    Grid: NNPDF23_lo_as_0130_qed_mem0.grid
   ****************************************
  New value of alpha_s from PDF nn23lo1: 0.13000000000000000
  Define smin to 0.0000000000000000
  Define smin to 0.0000000000000000
  Warning: pt or E min of a jet should in general be >0
  Warning: pt or E min of a jet should in general be >0
  Warning: pt or E min of a jet should in general be >0
  Warning: pt or E min of a jet should in general be >0
  *****************************************************
  * MadGraph/MadEvent *
  * -------------------------------- *
  * http://madgraph.hep.uiuc.edu *
  * http://madgraph.phys.ucl.ac.be *
  * http://madgraph.roma2.infn.it *
  * -------------------------------- *
  * *
  * PARAMETER AND COUPLING VALUES *
  * *
  *****************************************************

   External Params
   ---------------------------------

  mdl_cabi = 0.22773599999999999
  aEWM1 = 127.90000000000001
  mdl_Gf = 1.1663700000000000E-005
  aS = 0.11839999999999999
  mdl_ymb = 4.7000000000000002
  mdl_ymt = 172.00000000000000
  mdl_ymtau = 1.7769999999999999
  mdl_Lambda = 5000.0000000000000
  mdl_ca = 1.0000000000000000
  mdl_kSM = 5.0000000000000003E-002
  mdl_kHtt = 0.0000000000000000
  mdl_kAtt = 0.0000000000000000
  mdl_kHbb = 0.0000000000000000
  mdl_kAbb = 0.0000000000000000
  mdl_kHll = 0.0000000000000000
  mdl_kAll = 0.0000000000000000
  mdl_kHaa = 0.0000000000000000
  mdl_kAaa = 0.0000000000000000
  mdl_kHza = 0.0000000000000000
  mdl_kAza = 0.0000000000000000
  mdl_kHgg = 0.0000000000000000
  mdl_kAgg = 0.0000000000000000
  mdl_kHzz = 7.3405100000000001
  mdl_kAzz = 0.0000000000000000
  mdl_kHww = 9.4465800000000009
  mdl_kAww = 0.0000000000000000
  mdl_kHda = 0.0000000000000000
  mdl_kHdz = 0.0000000000000000
  mdl_kHdwR = 0.0000000000000000
  mdl_kHdwI = 0.0000000000000000
  mdl_kHHgg = 1.0000000000000000
  mdl_kAAgg = 1.0000000000000000
  mdl_kqa = 1.0000000000000000
  mdl_kqb = 1.0000000000000000
  mdl_kla = 1.0000000000000000
  mdl_klb = 1.0000000000000000
  mdl_kw1 = 1.0000000000000000
  mdl_kw2 = 1.0000000000000000
  mdl_kw3 = 0.0000000000000000
  mdl_kw4 = 0.0000000000000000
  mdl_kw5 = 0.0000000000000000
  mdl_kz1 = 0.0000000000000000
  mdl_kz3 = 1.0000000000000000
  mdl_kz5 = 0.0000000000000000
  mdl_kq = 1.0000000000000000
  mdl_kq3 = 1.0000000000000000
  mdl_kl = 1.0000000000000000
  mdl_kg = 1.0000000000000000
  mdl_ka = 1.0000000000000000
  mdl_kz = 1.0000000000000000
  mdl_kw = 1.0000000000000000
  mdl_kza = 0.0000000000000000
  mdl_MZ = 91.187600000000003
  mdl_MTA = 1.7769999999999999
  mdl_MT = 172.00000000000000
  mdl_MB = 4.7000000000000002
  mdl_MX0 = 900.00000000000000
  mdl_MX1 = 125.00000000000000
  mdl_MX2 = 125.00000000000000
  mdl_WZ = 2.4952000000000001
  mdl_WW = 2.0850000000000000
  mdl_WT = 1.5083359999999999
  mdl_WX0 = 33.258000000000003
  mdl_WX1 = 4.0699999999999998E-003
  mdl_WX2 = 4.0699999999999998E-003
   Internal Params
   ---------------------------------

  mdl_cos__cabi = 0.97418004031982097
  mdl_CKM1x1 = ( 0.97418004031982097 , 0.0000000000000000 )
  mdl_sin__cabi = 0.22577256042856930
  mdl_CKM1x2 = ( 0.22577256042856930 , 0.0000000000000000 )
  mdl_CKM2x1 = (-0.22577256042856930 , 0.0000000000000000 )
  mdl_CKM2x2 = ( 0.97418004031982097 , 0.0000000000000000 )
  mdl_ca__exp__2 = 1.0000000000000000
  mdl_sa = 0.0000000000000000
  mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
  mdl_kHdw = ( 0.0000000000000000 , 0.0000000000000000 )
  mdl_MZ__exp__2 = 8315.1783937600012
  mdl_MZ__exp__4 = 69142191.720053151
  mdl_sqrt__2 = 1.4142135623730951
  mdl_nb__2__exp__0_75 = 1.6817928305074290
  mdl_MX0__exp__2 = 810000.00000000000
  mdl_ca__exp__4 = 1.0000000000000000
  mdl_ca__exp__3 = 1.0000000000000000
  mdl_conjg__CKM1x1 = ( 0.97418004031982097 , -0.0000000000000000 )
  mdl_conjg__CKM1x2 = ( 0.22577256042856930 , -0.0000000000000000 )
  mdl_conjg__CKM2x1 = (-0.22577256042856930 , -0.0000000000000000 )
  mdl_conjg__CKM2x2 = ( 0.97418004031982097 , -0.0000000000000000 )
  mdl_conjg__kHdw = ( 0.0000000000000000 , -0.0000000000000000 )
  mdl_aEW = 7.8186082877247844E-003
  mdl_MW = 79.824359746197842
  mdl_sqrt__aEW = 8.8422894590285753E-002
  mdl_ee = 0.31345100004952897
  mdl_MW__exp__2 = 6371.9284088904105
  mdl_sw2 = 0.23369913342182447
  mdl_cw = 0.87538612427783857
  mdl_sqrt__sw2 = 0.48342438232036300
  mdl_sw = 0.48342438232036300
  mdl_ad = -0.12747442148678401
  mdl_al = -1.2074027013887732E-002
  mdl_an = 0.18517461872323218
  mdl_au = 6.9774224250335873E-002
  mdl_bd = -0.18517461872323218
  mdl_bl = -0.18517461872323218
  mdl_bn = 0.18517461872323218
  mdl_bu = 0.18517461872323218
  mdl_gwwz = -0.56759788318358417
  mdl_g1 = 0.35807170271074895
  mdl_gw = 0.64839716719502682
  mdl_vev = 246.22056907348590
  mdl_ee__exp__2 = 9.8251529432049817E-002
  mdl_gAaa = 1.3477023594718659E-005
  mdl_vev__exp__2 = 60624.568634871241
  mdl_cw__exp__2 = 0.76630086657817542
  mdl_gAza = 4.4999802188047602E-006
  mdl_gHaa = 2.6392504539657372E-005
  mdl_gHza = 3.9166441035016922E-005
  mdl_lam = 6.6804599046176811
  mdl_yb = 2.6995322804122722E-002
  mdl_yt = 0.98791394091683138
  mdl_ytau = 1.0206529494239589E-002
  mdl_muH = 636.39610306789280
  mdl_sw__exp__2 = 0.23369913342182450
   Internal Params evaluated point by point
   ----------------------------------------

  mdl_sqrt__aS = 0.34409301068170506
  mdl_G__exp__2 = 1.4878582807401259
  mdl_gAAgg = 3.1083017473614446E-007
  mdl_gAgg = 7.6532782508744533E-005
  mdl_gHgg = -5.1021855005829691E-005
  mdl_gHHgg = 2.0722011649076294E-007
   Couplings of HC_UFO
   ---------------------------------

          GC_6 = -0.12198E+01 0.00000E+00
          GC_7 = 0.00000E+00 0.12198E+01
          GC_8 = 0.00000E+00 0.14879E+01
         GC_10 = -0.00000E+00 0.00000E+00
         GC_11 = 0.00000E+00 0.00000E+00
         GC_12 = 0.00000E+00 -0.00000E+00
         GC_13 = -0.00000E+00 -0.20722E-06
         GC_62 = 0.00000E+00 0.00000E+00
         GC_64 = 0.00000E+00 0.00000E+00
         GC_65 = 0.00000E+00 0.00000E+00
          GC_1 = -0.00000E+00 -0.10448E+00
          GC_2 = 0.00000E+00 0.20897E+00
         GC_16 = -0.00000E+00 -0.00000E+00
         GC_66 = 0.00000E+00 0.00000E+00
         GC_77 = -0.00000E+00 -0.28380E+00
         GC_78 = 0.00000E+00 0.28380E+00
         GC_89 = -0.00000E+00 -0.28850E-01
         GC_90 = 0.00000E+00 0.86550E-01
        GC_100 = 0.00000E+00 0.33771E+01
         GC_38 = -0.00000E+00 -0.00000E+00
         GC_40 = -0.00000E+00 -0.00000E+00
         GC_42 = -0.00000E+00 -0.14681E-02
         GC_68 = 0.00000E+00 0.00000E+00

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

  Running at P P machine @ 13000.000000000000 GeV
  PDF set = nn23lo1
  alpha_s(Mz)= 0.1300 running at 2 loops.
  alpha_s(Mz)= 0.1300 running at 2 loops.
  Renormalization scale set on event-by-event basis
  Factorization scale set on event-by-event basis

  getting user params
 Enter number of events and max and min iterations:
  Number of events and iterations 1000 5 3
 Enter desired fractional accuracy:
  Desired fractional accuracy: 0.10000000000000001
 Enter 0 for fixed, 2 for adjustable grid:
 Suppress amplitude (0 no, 1 yes)?
  Using suppressed amplitude.
 Exact helicity sum (0 yes, n = number/event)?
  Explicitly summing over helicities
 Enter Configuration Number:
 Running Configuration Number: 6
  Not subdividing B.W.
  Attempting mappinvarients 1 8
  Completed mapping 8
  about to integrate 16 1000 5 3 16 1
  Using non-zero grid deformation.
  16 dimensions 1000 events 16 invarients 5 iterations 1 config(s), (0.99)
  Using h-tuple random number sequence.
  Error opening grid
  Using Uniform Grid! 32
  Using uniform alpha 1.0000000000000000
  Grid defined OK
  Set CM energy to 13000.00
  Mapping Graph 6 to config 6
 Setting grid 1 0.59172E-08 1
 Setting grid 2 0.59172E-08 1
  Setting BW -3 1 900.00000000000000
  Setting BW -4 2 91.187600000000003
 Setting grid 5 0.12244E-02 1
  Transforming s_hat 1/s 15 1.2244056105808281E-003 0.0000000000000000 168999999.99999997
    6 1 2 3 4 5 6 7 8 9 10 11 12 13 14
   15 16
  Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
 Using random seed offsets 6 : 2
   with seed 30
  Ranmar initialization seeds 8297 9406

  ********************************************
  * You are using the DiscreteSampler module *
  * part of the MG5_aMC framework *
  * Author: Valentin Hirschi *
  ********************************************

   Particle 3 4 5 6 7 8
       Et > 0.0 0.0 0.0 0.0 0.0 0.0
        E > 0.0 0.0 0.0 0.0 0.0 0.0
      Eta < -1.0 -1.0 -1.0 -1.0 -1.0 -1.0
    xqcut: 0.0 0.0 0.0 0.0 0.0 0.0
 d R # 3 > -0.0 0.0 0.0 0.0 0.0 0.0
 d R # 4 > -0.0 -0.0 0.0 0.0 0.0 0.0
 d R # 5 > -0.0 -0.0 -0.0 0.0 0.0 0.0
 d R # 6 > -0.0 -0.0 -0.0 -0.0 0.0 0.0
 d R # 7 > -0.0 -0.0 -0.0 -0.0 -0.0 0.0
 s min # 3> 0.0 0.0 0.0 0.0 0.0 0.0
 s min # 4> 0.0 0.0 0.0 0.0 0.0 0.0
 s min # 5> 0.0 0.0 0.0 0.0 0.0 0.0
 s min # 6> 0.0 0.0 0.0 0.0 0.0 0.0
 s min # 7> 0.0 0.0 0.0 0.0 0.0 0.0
 xqcutij # 3> 0.0 0.0 0.0 0.0 0.0 0.0
 xqcutij # 4> 0.0 0.0 0.0 0.0 0.0 0.0
 xqcutij # 5> 0.0 0.0 0.0 0.0 0.0 0.0
 xqcutij # 6> 0.0 0.0 0.0 0.0 0.0 0.0
 xqcutij # 7> 0.0 0.0 0.0 0.0 0.0 0.0
  No cut BW -3 1
  No cut BW -4 2
  alpha_s for scale 1103.7531049697764 is 9.2015413872191379E-002
  Error: Failed despite same graph: 6
  Have jets (>0) 0 0 0 0 1 0 1 1
  Should be 5 7
 STOP 4

 ls status:
 input_app.txt
 log.txt
 run1_app.log

Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/TEDWORK/HC_lljjjj/run_01_tag_1_debug.log'.
Please attach this file to your report.

How can I solve the error?

Cheers,
Yue

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2019-11-26
Last reply:
2019-11-26

Hi,

I will take a look. But in order to hope for such process to works, don't you need to have cut on the jets?

Cheers,

Olivier

> On 4 Nov 2019, at 18:02, Yue Xu <email address hidden> wrote:
>
> New question #685616 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/685616
>
> Hello,
>
> I want to generate some events with a modified HC model. This the process is(where x0 is a heavy higgs):
> import model HC_UFO
> generate p p > z x0, z > l- l+, x0 > j j j j
> output ./TEDWORK/HC_lljjjj
>
> The model is at:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_UFO
>
> The process and run_01_tag_1_debug.log is at:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_lljjjj
>
> When I generate events with "./bin/madevent", I get some errors:
> Error when reading /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/TEDWORK/HC_lljjjj/SubProcesses/P1_qq_zx0_z_ll_x0_ggqq/G6/results.dat
> Command "launch " interrupted with error:
> Exception : Reported error: End code 4.0
> Full associated log:
> Process in group number 1
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11839999999999999
> ****************************************
>
> NNPDFDriver version 1.0.3
> Grid: NNPDF23_lo_as_0130_qed_mem0.grid
> ****************************************
> New value of alpha_s from PDF nn23lo1: 0.13000000000000000
> Define smin to 0.0000000000000000
> Define smin to 0.0000000000000000
> Warning: pt or E min of a jet should in general be >0
> Warning: pt or E min of a jet should in general be >0
> Warning: pt or E min of a jet should in general be >0
> Warning: pt or E min of a jet should in general be >0
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> mdl_cabi = 0.22773599999999999
> aEWM1 = 127.90000000000001
> mdl_Gf = 1.1663700000000000E-005
> aS = 0.11839999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 172.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_Lambda = 5000.0000000000000
> mdl_ca = 1.0000000000000000
> mdl_kSM = 5.0000000000000003E-002
> mdl_kHtt = 0.0000000000000000
> mdl_kAtt = 0.0000000000000000
> mdl_kHbb = 0.0000000000000000
> mdl_kAbb = 0.0000000000000000
> mdl_kHll = 0.0000000000000000
> mdl_kAll = 0.0000000000000000
> mdl_kHaa = 0.0000000000000000
> mdl_kAaa = 0.0000000000000000
> mdl_kHza = 0.0000000000000000
> mdl_kAza = 0.0000000000000000
> mdl_kHgg = 0.0000000000000000
> mdl_kAgg = 0.0000000000000000
> mdl_kHzz = 7.3405100000000001
> mdl_kAzz = 0.0000000000000000
> mdl_kHww = 9.4465800000000009
> mdl_kAww = 0.0000000000000000
> mdl_kHda = 0.0000000000000000
> mdl_kHdz = 0.0000000000000000
> mdl_kHdwR = 0.0000000000000000
> mdl_kHdwI = 0.0000000000000000
> mdl_kHHgg = 1.0000000000000000
> mdl_kAAgg = 1.0000000000000000
> mdl_kqa = 1.0000000000000000
> mdl_kqb = 1.0000000000000000
> mdl_kla = 1.0000000000000000
> mdl_klb = 1.0000000000000000
> mdl_kw1 = 1.0000000000000000
> mdl_kw2 = 1.0000000000000000
> mdl_kw3 = 0.0000000000000000
> mdl_kw4 = 0.0000000000000000
> mdl_kw5 = 0.0000000000000000
> mdl_kz1 = 0.0000000000000000
> mdl_kz3 = 1.0000000000000000
> mdl_kz5 = 0.0000000000000000
> mdl_kq = 1.0000000000000000
> mdl_kq3 = 1.0000000000000000
> mdl_kl = 1.0000000000000000
> mdl_kg = 1.0000000000000000
> mdl_ka = 1.0000000000000000
> mdl_kz = 1.0000000000000000
> mdl_kw = 1.0000000000000000
> mdl_kza = 0.0000000000000000
> mdl_MZ = 91.187600000000003
> mdl_MTA = 1.7769999999999999
> mdl_MT = 172.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MX0 = 900.00000000000000
> mdl_MX1 = 125.00000000000000
> mdl_MX2 = 125.00000000000000
> mdl_WZ = 2.4952000000000001
> mdl_WW = 2.0850000000000000
> mdl_WT = 1.5083359999999999
> mdl_WX0 = 33.258000000000003
> mdl_WX1 = 4.0699999999999998E-003
> mdl_WX2 = 4.0699999999999998E-003
> Internal Params
> ---------------------------------
>
> mdl_cos__cabi = 0.97418004031982097
> mdl_CKM1x1 = ( 0.97418004031982097 , 0.0000000000000000 )
> mdl_sin__cabi = 0.22577256042856930
> mdl_CKM1x2 = ( 0.22577256042856930 , 0.0000000000000000 )
> mdl_CKM2x1 = (-0.22577256042856930 , 0.0000000000000000 )
> mdl_CKM2x2 = ( 0.97418004031982097 , 0.0000000000000000 )
> mdl_ca__exp__2 = 1.0000000000000000
> mdl_sa = 0.0000000000000000
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_kHdw = ( 0.0000000000000000 , 0.0000000000000000 )
> mdl_MZ__exp__2 = 8315.1783937600012
> mdl_MZ__exp__4 = 69142191.720053151
> mdl_sqrt__2 = 1.4142135623730951
> mdl_nb__2__exp__0_75 = 1.6817928305074290
> mdl_MX0__exp__2 = 810000.00000000000
> mdl_ca__exp__4 = 1.0000000000000000
> mdl_ca__exp__3 = 1.0000000000000000
> mdl_conjg__CKM1x1 = ( 0.97418004031982097 , -0.0000000000000000 )
> mdl_conjg__CKM1x2 = ( 0.22577256042856930 , -0.0000000000000000 )
> mdl_conjg__CKM2x1 = (-0.22577256042856930 , -0.0000000000000000 )
> mdl_conjg__CKM2x2 = ( 0.97418004031982097 , -0.0000000000000000 )
> mdl_conjg__kHdw = ( 0.0000000000000000 , -0.0000000000000000 )
> mdl_aEW = 7.8186082877247844E-003
> mdl_MW = 79.824359746197842
> mdl_sqrt__aEW = 8.8422894590285753E-002
> mdl_ee = 0.31345100004952897
> mdl_MW__exp__2 = 6371.9284088904105
> mdl_sw2 = 0.23369913342182447
> mdl_cw = 0.87538612427783857
> mdl_sqrt__sw2 = 0.48342438232036300
> mdl_sw = 0.48342438232036300
> mdl_ad = -0.12747442148678401
> mdl_al = -1.2074027013887732E-002
> mdl_an = 0.18517461872323218
> mdl_au = 6.9774224250335873E-002
> mdl_bd = -0.18517461872323218
> mdl_bl = -0.18517461872323218
> mdl_bn = 0.18517461872323218
> mdl_bu = 0.18517461872323218
> mdl_gwwz = -0.56759788318358417
> mdl_g1 = 0.35807170271074895
> mdl_gw = 0.64839716719502682
> mdl_vev = 246.22056907348590
> mdl_ee__exp__2 = 9.8251529432049817E-002
> mdl_gAaa = 1.3477023594718659E-005
> mdl_vev__exp__2 = 60624.568634871241
> mdl_cw__exp__2 = 0.76630086657817542
> mdl_gAza = 4.4999802188047602E-006
> mdl_gHaa = 2.6392504539657372E-005
> mdl_gHza = 3.9166441035016922E-005
> mdl_lam = 6.6804599046176811
> mdl_yb = 2.6995322804122722E-002
> mdl_yt = 0.98791394091683138
> mdl_ytau = 1.0206529494239589E-002
> mdl_muH = 636.39610306789280
> mdl_sw__exp__2 = 0.23369913342182450
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34409301068170506
> mdl_G__exp__2 = 1.4878582807401259
> mdl_gAAgg = 3.1083017473614446E-007
> mdl_gAgg = 7.6532782508744533E-005
> mdl_gHgg = -5.1021855005829691E-005
> mdl_gHHgg = 2.0722011649076294E-007
> Couplings of HC_UFO
> ---------------------------------
>
> GC_6 = -0.12198E+01 0.00000E+00
> GC_7 = 0.00000E+00 0.12198E+01
> GC_8 = 0.00000E+00 0.14879E+01
> GC_10 = -0.00000E+00 0.00000E+00
> GC_11 = 0.00000E+00 0.00000E+00
> GC_12 = 0.00000E+00 -0.00000E+00
> GC_13 = -0.00000E+00 -0.20722E-06
> GC_62 = 0.00000E+00 0.00000E+00
> GC_64 = 0.00000E+00 0.00000E+00
> GC_65 = 0.00000E+00 0.00000E+00
> GC_1 = -0.00000E+00 -0.10448E+00
> GC_2 = 0.00000E+00 0.20897E+00
> GC_16 = -0.00000E+00 -0.00000E+00
> GC_66 = 0.00000E+00 0.00000E+00
> GC_77 = -0.00000E+00 -0.28380E+00
> GC_78 = 0.00000E+00 0.28380E+00
> GC_89 = -0.00000E+00 -0.28850E-01
> GC_90 = 0.00000E+00 0.86550E-01
> GC_100 = 0.00000E+00 0.33771E+01
> GC_38 = -0.00000E+00 -0.00000E+00
> GC_40 = -0.00000E+00 -0.00000E+00
> GC_42 = -0.00000E+00 -0.14681E-02
> GC_68 = 0.00000E+00 0.00000E+00
>
> Collider parameters:
> --------------------
>
> Running at P P machine @ 13000.000000000000 GeV
> PDF set = nn23lo1
> alpha_s(Mz)= 0.1300 running at 2 loops.
> alpha_s(Mz)= 0.1300 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> getting user params
> Enter number of events and max and min iterations:
> Number of events and iterations 1000 5 3
> Enter desired fractional accuracy:
> Desired fractional accuracy: 0.10000000000000001
> Enter 0 for fixed, 2 for adjustable grid:
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Explicitly summing over helicities
> Enter Configuration Number:
> Running Configuration Number: 6
> Not subdividing B.W.
> Attempting mappinvarients 1 8
> Completed mapping 8
> about to integrate 16 1000 5 3 16 1
> Using non-zero grid deformation.
> 16 dimensions 1000 events 16 invarients 5 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> Error opening grid
> Using Uniform Grid! 32
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 13000.00
> Mapping Graph 6 to config 6
> Setting grid 1 0.59172E-08 1
> Setting grid 2 0.59172E-08 1
> Setting BW -3 1 900.00000000000000
> Setting BW -4 2 91.187600000000003
> Setting grid 5 0.12244E-02 1
> Transforming s_hat 1/s 15 1.2244056105808281E-003 0.0000000000000000 168999999.99999997
> 6 1 2 3 4 5 6 7 8 9 10 11 12 13 14
> 15 16
> Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
> Using random seed offsets 6 : 2
> with seed 30
> Ranmar initialization seeds 8297 9406
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Particle 3 4 5 6 7 8
> Et > 0.0 0.0 0.0 0.0 0.0 0.0
> E > 0.0 0.0 0.0 0.0 0.0 0.0
> Eta < -1.0 -1.0 -1.0 -1.0 -1.0 -1.0
> xqcut: 0.0 0.0 0.0 0.0 0.0 0.0
> d R # 3 > -0.0 0.0 0.0 0.0 0.0 0.0
> d R # 4 > -0.0 -0.0 0.0 0.0 0.0 0.0
> d R # 5 > -0.0 -0.0 -0.0 0.0 0.0 0.0
> d R # 6 > -0.0 -0.0 -0.0 -0.0 0.0 0.0
> d R # 7 > -0.0 -0.0 -0.0 -0.0 -0.0 0.0
> s min # 3> 0.0 0.0 0.0 0.0 0.0 0.0
> s min # 4> 0.0 0.0 0.0 0.0 0.0 0.0
> s min # 5> 0.0 0.0 0.0 0.0 0.0 0.0
> s min # 6> 0.0 0.0 0.0 0.0 0.0 0.0
> s min # 7> 0.0 0.0 0.0 0.0 0.0 0.0
> xqcutij # 3> 0.0 0.0 0.0 0.0 0.0 0.0
> xqcutij # 4> 0.0 0.0 0.0 0.0 0.0 0.0
> xqcutij # 5> 0.0 0.0 0.0 0.0 0.0 0.0
> xqcutij # 6> 0.0 0.0 0.0 0.0 0.0 0.0
> xqcutij # 7> 0.0 0.0 0.0 0.0 0.0 0.0
> No cut BW -3 1
> No cut BW -4 2
> alpha_s for scale 1103.7531049697764 is 9.2015413872191379E-002
> Error: Failed despite same graph: 6
> Have jets (>0) 0 0 0 0 1 0 1 1
> Should be 5 7
> STOP 4
>
> ls status:
> input_app.txt
> log.txt
> run1_app.log
>
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/TEDWORK/HC_lljjjj/run_01_tag_1_debug.log'.
> Please attach this file to your report.
>
> How can I solve the error?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

One work-around for this issue
is
1) allow the cut on the jets:
(set cut_decays True)
2) change the default dynamical to not CKKW mode
(set dynamical_scale_choice 3) #or any other value

Cheers,

Olivier

> On 5 Nov 2019, at 21:15, Olivier Mattelaer <email address hidden> wrote:
>
> Hi,
>
> I will take a look. But in order to hope for such process to works, don't you need to have cut on the jets?
>
> Cheers,
>
> Olivier
>
>> On 4 Nov 2019, at 18:02, Yue Xu <email address hidden> wrote:
>>
>> New question #685616 on MadGraph5_aMC@NLO:
>> https://answers.launchpad.net/mg5amcnlo/+question/685616
>>
>> Hello,
>>
>> I want to generate some events with a modified HC model. This the process is(where x0 is a heavy higgs):
>> import model HC_UFO
>> generate p p > z x0, z > l- l+, x0 > j j j j
>> output ./TEDWORK/HC_lljjjj
>>
>> The model is at:
>> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_UFO
>>
>> The process and run_01_tag_1_debug.log is at:
>> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_lljjjj
>>
>> When I generate events with "./bin/madevent", I get some errors:
>> Error when reading /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/TEDWORK/HC_lljjjj/SubProcesses/P1_qq_zx0_z_ll_x0_ggqq/G6/results.dat
>> Command "launch " interrupted with error:
>> Exception : Reported error: End code 4.0
>> Full associated log:
>> Process in group number 1
>> A PDF is used, so alpha_s(MZ) is going to be modified
>> Old value of alpha_s from param_card: 0.11839999999999999
>> ****************************************
>>
>> NNPDFDriver version 1.0.3
>> Grid: NNPDF23_lo_as_0130_qed_mem0.grid
>> ****************************************
>> New value of alpha_s from PDF nn23lo1: 0.13000000000000000
>> Define smin to 0.0000000000000000
>> Define smin to 0.0000000000000000
>> Warning: pt or E min of a jet should in general be >0
>> Warning: pt or E min of a jet should in general be >0
>> Warning: pt or E min of a jet should in general be >0
>> Warning: pt or E min of a jet should in general be >0
>> *****************************************************
>> * MadGraph/MadEvent *
>> * -------------------------------- *
>> * http://madgraph.hep.uiuc.edu *
>> * http://madgraph.phys.ucl.ac.be *
>> * http://madgraph.roma2.infn.it *
>> * -------------------------------- *
>> * *
>> * PARAMETER AND COUPLING VALUES *
>> * *
>> *****************************************************
>>
>> External Params
>> ---------------------------------
>>
>> mdl_cabi = 0.22773599999999999
>> aEWM1 = 127.90000000000001
>> mdl_Gf = 1.1663700000000000E-005
>> aS = 0.11839999999999999
>> mdl_ymb = 4.7000000000000002
>> mdl_ymt = 172.00000000000000
>> mdl_ymtau = 1.7769999999999999
>> mdl_Lambda = 5000.0000000000000
>> mdl_ca = 1.0000000000000000
>> mdl_kSM = 5.0000000000000003E-002
>> mdl_kHtt = 0.0000000000000000
>> mdl_kAtt = 0.0000000000000000
>> mdl_kHbb = 0.0000000000000000
>> mdl_kAbb = 0.0000000000000000
>> mdl_kHll = 0.0000000000000000
>> mdl_kAll = 0.0000000000000000
>> mdl_kHaa = 0.0000000000000000
>> mdl_kAaa = 0.0000000000000000
>> mdl_kHza = 0.0000000000000000
>> mdl_kAza = 0.0000000000000000
>> mdl_kHgg = 0.0000000000000000
>> mdl_kAgg = 0.0000000000000000
>> mdl_kHzz = 7.3405100000000001
>> mdl_kAzz = 0.0000000000000000
>> mdl_kHww = 9.4465800000000009
>> mdl_kAww = 0.0000000000000000
>> mdl_kHda = 0.0000000000000000
>> mdl_kHdz = 0.0000000000000000
>> mdl_kHdwR = 0.0000000000000000
>> mdl_kHdwI = 0.0000000000000000
>> mdl_kHHgg = 1.0000000000000000
>> mdl_kAAgg = 1.0000000000000000
>> mdl_kqa = 1.0000000000000000
>> mdl_kqb = 1.0000000000000000
>> mdl_kla = 1.0000000000000000
>> mdl_klb = 1.0000000000000000
>> mdl_kw1 = 1.0000000000000000
>> mdl_kw2 = 1.0000000000000000
>> mdl_kw3 = 0.0000000000000000
>> mdl_kw4 = 0.0000000000000000
>> mdl_kw5 = 0.0000000000000000
>> mdl_kz1 = 0.0000000000000000
>> mdl_kz3 = 1.0000000000000000
>> mdl_kz5 = 0.0000000000000000
>> mdl_kq = 1.0000000000000000
>> mdl_kq3 = 1.0000000000000000
>> mdl_kl = 1.0000000000000000
>> mdl_kg = 1.0000000000000000
>> mdl_ka = 1.0000000000000000
>> mdl_kz = 1.0000000000000000
>> mdl_kw = 1.0000000000000000
>> mdl_kza = 0.0000000000000000
>> mdl_MZ = 91.187600000000003
>> mdl_MTA = 1.7769999999999999
>> mdl_MT = 172.00000000000000
>> mdl_MB = 4.7000000000000002
>> mdl_MX0 = 900.00000000000000
>> mdl_MX1 = 125.00000000000000
>> mdl_MX2 = 125.00000000000000
>> mdl_WZ = 2.4952000000000001
>> mdl_WW = 2.0850000000000000
>> mdl_WT = 1.5083359999999999
>> mdl_WX0 = 33.258000000000003
>> mdl_WX1 = 4.0699999999999998E-003
>> mdl_WX2 = 4.0699999999999998E-003
>> Internal Params
>> ---------------------------------
>>
>> mdl_cos__cabi = 0.97418004031982097
>> mdl_CKM1x1 = ( 0.97418004031982097 , 0.0000000000000000 )
>> mdl_sin__cabi = 0.22577256042856930
>> mdl_CKM1x2 = ( 0.22577256042856930 , 0.0000000000000000 )
>> mdl_CKM2x1 = (-0.22577256042856930 , 0.0000000000000000 )
>> mdl_CKM2x2 = ( 0.97418004031982097 , 0.0000000000000000 )
>> mdl_ca__exp__2 = 1.0000000000000000
>> mdl_sa = 0.0000000000000000
>> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
>> mdl_kHdw = ( 0.0000000000000000 , 0.0000000000000000 )
>> mdl_MZ__exp__2 = 8315.1783937600012
>> mdl_MZ__exp__4 = 69142191.720053151
>> mdl_sqrt__2 = 1.4142135623730951
>> mdl_nb__2__exp__0_75 = 1.6817928305074290
>> mdl_MX0__exp__2 = 810000.00000000000
>> mdl_ca__exp__4 = 1.0000000000000000
>> mdl_ca__exp__3 = 1.0000000000000000
>> mdl_conjg__CKM1x1 = ( 0.97418004031982097 , -0.0000000000000000 )
>> mdl_conjg__CKM1x2 = ( 0.22577256042856930 , -0.0000000000000000 )
>> mdl_conjg__CKM2x1 = (-0.22577256042856930 , -0.0000000000000000 )
>> mdl_conjg__CKM2x2 = ( 0.97418004031982097 , -0.0000000000000000 )
>> mdl_conjg__kHdw = ( 0.0000000000000000 , -0.0000000000000000 )
>> mdl_aEW = 7.8186082877247844E-003
>> mdl_MW = 79.824359746197842
>> mdl_sqrt__aEW = 8.8422894590285753E-002
>> mdl_ee = 0.31345100004952897
>> mdl_MW__exp__2 = 6371.9284088904105
>> mdl_sw2 = 0.23369913342182447
>> mdl_cw = 0.87538612427783857
>> mdl_sqrt__sw2 = 0.48342438232036300
>> mdl_sw = 0.48342438232036300
>> mdl_ad = -0.12747442148678401
>> mdl_al = -1.2074027013887732E-002
>> mdl_an = 0.18517461872323218
>> mdl_au = 6.9774224250335873E-002
>> mdl_bd = -0.18517461872323218
>> mdl_bl = -0.18517461872323218
>> mdl_bn = 0.18517461872323218
>> mdl_bu = 0.18517461872323218
>> mdl_gwwz = -0.56759788318358417
>> mdl_g1 = 0.35807170271074895
>> mdl_gw = 0.64839716719502682
>> mdl_vev = 246.22056907348590
>> mdl_ee__exp__2 = 9.8251529432049817E-002
>> mdl_gAaa = 1.3477023594718659E-005
>> mdl_vev__exp__2 = 60624.568634871241
>> mdl_cw__exp__2 = 0.76630086657817542
>> mdl_gAza = 4.4999802188047602E-006
>> mdl_gHaa = 2.6392504539657372E-005
>> mdl_gHza = 3.9166441035016922E-005
>> mdl_lam = 6.6804599046176811
>> mdl_yb = 2.6995322804122722E-002
>> mdl_yt = 0.98791394091683138
>> mdl_ytau = 1.0206529494239589E-002
>> mdl_muH = 636.39610306789280
>> mdl_sw__exp__2 = 0.23369913342182450
>> Internal Params evaluated point by point
>> ----------------------------------------
>>
>> mdl_sqrt__aS = 0.34409301068170506
>> mdl_G__exp__2 = 1.4878582807401259
>> mdl_gAAgg = 3.1083017473614446E-007
>> mdl_gAgg = 7.6532782508744533E-005
>> mdl_gHgg = -5.1021855005829691E-005
>> mdl_gHHgg = 2.0722011649076294E-007
>> Couplings of HC_UFO
>> ---------------------------------
>>
>> GC_6 = -0.12198E+01 0.00000E+00
>> GC_7 = 0.00000E+00 0.12198E+01
>> GC_8 = 0.00000E+00 0.14879E+01
>> GC_10 = -0.00000E+00 0.00000E+00
>> GC_11 = 0.00000E+00 0.00000E+00
>> GC_12 = 0.00000E+00 -0.00000E+00
>> GC_13 = -0.00000E+00 -0.20722E-06
>> GC_62 = 0.00000E+00 0.00000E+00
>> GC_64 = 0.00000E+00 0.00000E+00
>> GC_65 = 0.00000E+00 0.00000E+00
>> GC_1 = -0.00000E+00 -0.10448E+00
>> GC_2 = 0.00000E+00 0.20897E+00
>> GC_16 = -0.00000E+00 -0.00000E+00
>> GC_66 = 0.00000E+00 0.00000E+00
>> GC_77 = -0.00000E+00 -0.28380E+00
>> GC_78 = 0.00000E+00 0.28380E+00
>> GC_89 = -0.00000E+00 -0.28850E-01
>> GC_90 = 0.00000E+00 0.86550E-01
>> GC_100 = 0.00000E+00 0.33771E+01
>> GC_38 = -0.00000E+00 -0.00000E+00
>> GC_40 = -0.00000E+00 -0.00000E+00
>> GC_42 = -0.00000E+00 -0.14681E-02
>> GC_68 = 0.00000E+00 0.00000E+00
>>
>> Collider parameters:
>> --------------------
>>
>> Running at P P machine @ 13000.000000000000 GeV
>> PDF set = nn23lo1
>> alpha_s(Mz)= 0.1300 running at 2 loops.
>> alpha_s(Mz)= 0.1300 running at 2 loops.
>> Renormalization scale set on event-by-event basis
>> Factorization scale set on event-by-event basis
>>
>>
>> getting user params
>> Enter number of events and max and min iterations:
>> Number of events and iterations 1000 5 3
>> Enter desired fractional accuracy:
>> Desired fractional accuracy: 0.10000000000000001
>> Enter 0 for fixed, 2 for adjustable grid:
>> Suppress amplitude (0 no, 1 yes)?
>> Using suppressed amplitude.
>> Exact helicity sum (0 yes, n = number/event)?
>> Explicitly summing over helicities
>> Enter Configuration Number:
>> Running Configuration Number: 6
>> Not subdividing B.W.
>> Attempting mappinvarients 1 8
>> Completed mapping 8
>> about to integrate 16 1000 5 3 16 1
>> Using non-zero grid deformation.
>> 16 dimensions 1000 events 16 invarients 5 iterations 1 config(s), (0.99)
>> Using h-tuple random number sequence.
>> Error opening grid
>> Using Uniform Grid! 32
>> Using uniform alpha 1.0000000000000000
>> Grid defined OK
>> Set CM energy to 13000.00
>> Mapping Graph 6 to config 6
>> Setting grid 1 0.59172E-08 1
>> Setting grid 2 0.59172E-08 1
>> Setting BW -3 1 900.00000000000000
>> Setting BW -4 2 91.187600000000003
>> Setting grid 5 0.12244E-02 1
>> Transforming s_hat 1/s 15 1.2244056105808281E-003 0.0000000000000000 168999999.99999997
>> 6 1 2 3 4 5 6 7 8 9 10 11 12 13 14
>> 15 16
>> Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
>> Using random seed offsets 6 : 2
>> with seed 30
>> Ranmar initialization seeds 8297 9406
>>
>> ********************************************
>> * You are using the DiscreteSampler module *
>> * part of the MG5_aMC framework *
>> * Author: Valentin Hirschi *
>> ********************************************
>>
>> Particle 3 4 5 6 7 8
>> Et > 0.0 0.0 0.0 0.0 0.0 0.0
>> E > 0.0 0.0 0.0 0.0 0.0 0.0
>> Eta < -1.0 -1.0 -1.0 -1.0 -1.0 -1.0
>> xqcut: 0.0 0.0 0.0 0.0 0.0 0.0
>> d R # 3 > -0.0 0.0 0.0 0.0 0.0 0.0
>> d R # 4 > -0.0 -0.0 0.0 0.0 0.0 0.0
>> d R # 5 > -0.0 -0.0 -0.0 0.0 0.0 0.0
>> d R # 6 > -0.0 -0.0 -0.0 -0.0 0.0 0.0
>> d R # 7 > -0.0 -0.0 -0.0 -0.0 -0.0 0.0
>> s min # 3> 0.0 0.0 0.0 0.0 0.0 0.0
>> s min # 4> 0.0 0.0 0.0 0.0 0.0 0.0
>> s min # 5> 0.0 0.0 0.0 0.0 0.0 0.0
>> s min # 6> 0.0 0.0 0.0 0.0 0.0 0.0
>> s min # 7> 0.0 0.0 0.0 0.0 0.0 0.0
>> xqcutij # 3> 0.0 0.0 0.0 0.0 0.0 0.0
>> xqcutij # 4> 0.0 0.0 0.0 0.0 0.0 0.0
>> xqcutij # 5> 0.0 0.0 0.0 0.0 0.0 0.0
>> xqcutij # 6> 0.0 0.0 0.0 0.0 0.0 0.0
>> xqcutij # 7> 0.0 0.0 0.0 0.0 0.0 0.0
>> No cut BW -3 1
>> No cut BW -4 2
>> alpha_s for scale 1103.7531049697764 is 9.2015413872191379E-002
>> Error: Failed despite same graph: 6
>> Have jets (>0) 0 0 0 0 1 0 1 1
>> Should be 5 7
>> STOP 4
>>
>> ls status:
>> input_app.txt
>> log.txt
>> run1_app.log
>>
>> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
>> More information is found in '/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/TEDWORK/HC_lljjjj/run_01_tag_1_debug.log'.
>> Please attach this file to your report.
>>
>> How can I solve the error?
>>
>> Cheers,
>> Yue
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>

Yue Xu (yuexu) said : #3

Will this cut change the pt distribution of final state? Or change something else ?

Cheers,
Yue

Yue Xu (yuexu) said : #4

And this process can be generated in older MG5 version, such 2.5.5, without settings of cut_decays and dynamical_scale_choice.
But in the new version 2.6.6, I got the above error.

Cheers,
Yue

Hi,

1) setting "cut_decays" on True.
will allow cut like ptj, drjj, ... of the run_card to be applied on the jet comming from x0.

Yes this will change your distribution (if those cut are present, if you do not have any cut within the run_card then it will not change anything).
Now my understanding of your process (but I can be wrong) is that your cross-section should be infinite without any cut preventing two jet to be colinear or one jet to be soft. So those cuts are crucial to have physical results (but if for some reason, that I missed, you do not have such divergencies for this process/model)

2) Changing the dynamical scale choice will change the running of alpha_s.
(more precisely to which scale alpha_s is evaluated for each event)
So it will change the shape of the distribution but ONLY within scale uncertainties.
So it should not change anything to your analysis if you propagate not only the central value but also the theoretical uncertainties.

Cheers,

Olivier

> On 12 Nov 2019, at 10:11, Yue Xu <email address hidden> wrote:
>
> Question #685616 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685616
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Will this cut change the pt distribution of final state? Or change
> something else ?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

HI,

For cut_decays, I still keep my opinion.
I do not see why you do not have soft/collinear singularities in this case.
Now if it works fine in 2.5.5 it might have a reason why you do not have such singularities.
In that case you do not need to set it to True in 2.6.x.
I just hope that you understand why you do not face such singularities.

For the dynamical_scale_choice, this indeed sounds as a new behavior and as I said I will investigate what to do to fix it. By experience this CKKW merging is quite complex in presence of non QCD interaction so it will take a while to fix. This is why I suggest this workaround.
Note that I can not guarantee that the version 2.5.5 works correctly for the default dynamical scale (up to the fact that this is within theoretical uncertainty by definition). Especially since that part of the code faces a lot of bug/fix since 2017.
So for your process, I would recommend to use the dynamical_scale_choice different of -1 even with 2.5.5.

Cheers,

Olivier

> On 12 Nov 2019, at 10:43, Yue Xu <email address hidden> wrote:
>
> Question #685616 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685616
>
> Yue Xu gave more information on the question:
> And this process can be generated in older MG5 version, such 2.5.5, without settings of cut_decays and dynamical_scale_choice.
> But in the new version 2.6.6, I got the above error.
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #7

Hi Olivier,

When I change the default dynamical_scale_choice to 3, it works well !
But I have another question about the generation time. I have 4 processes, A, B, C and D.
Process A:
import model HC_UFO
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-
generate p p > z x0, z > j j, x0 > l- l+ j j

Process B:
import model HC_UFO
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-
generate p p > z x0, z > l- l+, x0 > j j j j
add process p p > z x0, z > l- l+, x0 > l- vl~ j j
add process p p > z x0, z > l- l+, x0 > l+ vl j j

Process C:
import model HC_UFO
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-
generate p p > w+ x0, w+ > j j, x0 > l- l+ j j
add process p p > w+ x0, w+ > l+ vl, x0 > l- l+ j j
add process p p > w+ x0 > w+ w+ w-, w+ > l+ vl, w- > j j

Process D:
import model HC_UFO
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-
generate p p > w- x0, w- > j j, x0 > l- l+ j j
add process p p > w- x0, w- > l- vl~, x0 > l- l+ j j
add process p p > w- x0 > w- w- w+, w- > l- vl~, w+ > j j

I run the 4 processes in cern lxplus with 8 cores with same setting of parameters and event number.

Process A has 3135 diagrams. Process B has 2662 diagrams. Process C and process D both have 5832 diagrams.

In principle, the generation time is proportional to the number of diagrams.
But process A takes more than 10 hours, while process B takes about 10 minutes and process C/D takes about 40 minutes.
It's strange that process A takes such a long time.

Why does process A take such a long time?

Cheers,
Yue

P.S:
Four output log files can be found here:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/TimingMeasurement/logfilesForTiming

Process A, B, C and D are here:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/TimingMeasurement

Hi,

The timing is indeed linked to the speed of the matrix-element.
But this is not the only factor, a quite determinant factor is the number of times that we need to
evaluate the matrix-element in order to reach the requested accuracy.
This is linked to the variance of the function that you integrating (after the change of variable done analytically and numerically --via machine learning-- in order to reduce such variance)

Cheers,

Olivier

In top of that lxplus is an heterogeous cluster and therefore you can have huge variation of speed depending of the cpu generation that you are using.

> On 22 Nov 2019, at 14:46, Yue Xu <email address hidden> wrote:
>
> Question #685616 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685616
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi Olivier,
>
> When I change the default dynamical_scale_choice to 3, it works well !
> But I have another question about the generation time. I have 4 processes, A, B, C and D.
> Process A:
> import model HC_UFO
> 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-
> generate p p > z x0, z > j j, x0 > l- l+ j j
>
> Process B:
> import model HC_UFO
> 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-
> generate p p > z x0, z > l- l+, x0 > j j j j
> add process p p > z x0, z > l- l+, x0 > l- vl~ j j
> add process p p > z x0, z > l- l+, x0 > l+ vl j j
>
> Process C:
> import model HC_UFO
> 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-
> generate p p > w+ x0, w+ > j j, x0 > l- l+ j j
> add process p p > w+ x0, w+ > l+ vl, x0 > l- l+ j j
> add process p p > w+ x0 > w+ w+ w-, w+ > l+ vl, w- > j j
>
> Process D:
> import model HC_UFO
> 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-
> generate p p > w- x0, w- > j j, x0 > l- l+ j j
> add process p p > w- x0, w- > l- vl~, x0 > l- l+ j j
> add process p p > w- x0 > w- w- w+, w- > l- vl~, w+ > j j
>
> I run the 4 processes in cern lxplus with 8 cores with same setting of
> parameters and event number.
>
> Process A has 3135 diagrams. Process B has 2662 diagrams. Process C and
> process D both have 5832 diagrams.
>
> In principle, the generation time is proportional to the number of diagrams.
> But process A takes more than 10 hours, while process B takes about 10 minutes and process C/D takes about 40 minutes.
> It's strange that process A takes such a long time.
>
> Why does process A take such a long time?
>
> Cheers,
> Yue
>
> P.S:
> Four output log files can be found here:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/TimingMeasurement/logfilesForTiming
>
> Process A, B, C and D are here:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/TimingMeasurement
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #9

Hi Olivier,

I don't understand why the number of time that we need to evaluate the matrix-element in order to reach the requested accuracy in case A is much larger.

I have tried several times. Each time, case A takes a much longer time. I also tested the timing in a standalone MadGraph, and also case A takes a much longer time.

Cheers,
Yue

Hi,

You can read the slides of my lecture on this topic:
https://cp3.irmp.ucl.ac.be/projects/madgraph/raw-attachment/wiki/Madrid2015_2%20MC/15_09_07_Madrid_MC.pdf

from page ~60 to ~120

Cheers,

Olivier

> On 26 Nov 2019, at 03:27, Yue Xu <email address hidden> wrote:
>
> Question #685616 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685616
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi Olivier,
>
> I don't understand why the number of time that we need to evaluate the
> matrix-element in order to reach the requested accuracy in case A is
> much larger.
>
> I have tried several times. Each time, case A takes a much longer time.
> I also tested the timing in a standalone MadGraph, and also case A takes
> a much longer time.
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

Provide an answer of your own, or ask Yue Xu for more information if necessary.

To post a message you must log in.