generate events with a modified HC model

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.
>

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

Cheers,
Yue

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.

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.

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.