Reported error: End code 1.0

Hi,

Looks like another bug report that I fixed induced a specific issue with the EWA case.
The following patch: https://github.com/mg5amcnlo/mg5amcnlo/commit/c6898224805a1cdcc464daf24f4ca629aa4bc242 <https://github.com/mg5amcnlo/mg5amcnlo/commit/c6898224805a1cdcc464daf24f4ca629aa4bc242>
does fix that issue.

Now when running your script I still have a lot of:

 muf2 too small. setting muf2 to muf2min: 7648.9402763025619 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 7648.9402763025619 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 8216.6374936658685 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 8216.6374936658685 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 8224.5957825656878 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 8224.5957825656878 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 7620.3515912094390 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 7620.3515912094390 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 7543.0216568722153 8315.2513440000002
 muf2 too small. setting muf2 to muf2min: 7543.0216568722153 8315.2513440000002

Which I guess makes sense since in your script, I do not see any of the cuts that we recomend
to be in the area of the validity of the EVA (in particular asking that shat is larger than 1TeV).
I guess that you have a lot of event close to threshold and all those should create the above warning
(and likely indicates that you try to use EVA outside of the theoretical validity regime)

Cheers,

Olivier

> On 26 May 2022, at 10:35, Yuan <email address hidden> wrote:
>
> New question #701968 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/701968
>
> Dear Olivier,
>
> I'm the beginner Yuan, who was helped by you a few ago, I have a problem again and need your help. I am trying to generate SM process: w+ z > mu+ mu+ mu- vm and w+ a > mu+ mu+ mu- vm, in the muon collider. But it is failed.
>
> __________________________________________________________________________
> Here are my command in Madgraph v3.4.0:
>
> set group_subprocesses false
> define vap = w+ z a
> define vam = w- z a
> generate vap vam > mu+ mu+ mu- vm
> output Muon_BSM/TypeII/Xsec_decay/HppHm/scan-vbf-3mu-1.5TeV
> y
> launch
> shower=OFF
> detector=OFF
> analysis=OFF
> madspin=OFF
> done
> set ebeam 1500
> set nevents 1000
> set lpp1 -4
> set lpp2 4
> set pdlabel eva
> set fixed_fac_scale false
> set dynamical_scale_choice 4
> set scalefact 0.5
> set use_syst true
> set nhel 1
> done
>
>
> __________________________________________________________________________
> ################
> The results:
> ################
> __________________________________________________________________________
>
> Error when reading /Users/maxwell/Software/MG5_aMC_v3_4_0_RC/Muon_BSM/TypeII/Xsec_decay/HppHm/scan-vbf-3mu-1.5TeV/SubProcesses/P1_wpa_mupmupmumvm/G39/results.dat
> Command "generate_events run_01" interrupted with error:
> Exception : Reported error: End code 1.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.11799999999999999
> New value of alpha_s from PDF eva : 0.11799999999999999
> Define smin to 900.00000000000000
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MZ = 91.188000000000002
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MH = 125.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WT = 1.4915000000000000
> mdl_WH = 6.3823389999999999E-003
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM1x1 = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_complexi = (0.0000000000000000,1.0000000000000000)
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_I1x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_I2x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I3x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I4x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_cw__exp__2 = 0.77775351421422245
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> Couplings of sm
> ---------------------------------
>
> GC_3 = -0.00000E+00 -0.30795E+00
> GC_5 = 0.00000E+00 0.94836E-01
> GC_36 = 0.00000E+00 0.33188E+00
> GC_50 = -0.00000E+00 -0.28804E+00
> GC_53 = 0.00000E+00 0.57609E+00
> GC_57 = -0.00000E+00 -0.35482E+00
> GC_59 = 0.00000E+00 0.82310E-01
> GC_62 = 0.00000E+00 0.37035E+00
> GC_72 = 0.00000E+00 0.52532E+02
> GC_81 = 0.00000E+00 0.67544E+02
> GC_100 = 0.00000E+00 0.46191E+00
>
> Collider parameters:
> --------------------
>
> Running at m+m- machine @ 3000.0000000000000 GeV
> PDF set = eva
> alpha_s(Mz)= 0.1180 running at 2 loops.
> alpha_s(Mz)= 0.1180 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 4000 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)?
> Monte-Carlo over helicities
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Enter Configuration Number:
> Running Configuration Number: 39
> Not subdividing B.W.
> Attempting mappinvarients 1 6
> Determine nb_t
> T-channel found: 2
> Completed mapping 6
> about to integrate 10 4000 7 3 10 1
> Using non-zero grid deformation.
> 10 dimensions 4000 events 10 invarients 7 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> Error opening grid
> Using Uniform Grid! 24
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 3000.00
> Mapping Graph 39 to config 39
> Determine nb_t
> T-channel found: 2
> Setting grid 1 0.14222E-05 1
> Warning: No cutoff for shat integral found
> Minimum set to 1.0000000000000000E-004
> Setting grid 2 0.10000E-03 1
> Setting grid 3 0.11111E-04 1
> Transforming s_hat 1/s 9 1.0000000000000000E-004 900.00000000000000 9000000.0000000000
> 39 1 2 3 4 5 6 7 8 9 10
> Masses: 0.804E+02 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
> Using random seed offsets 39 : 2
> with seed 24
> Ranmar initialization seeds 31097 9399
> Particle 3 4 5 6
> Et > 10.0 10.0 10.0 0.0
> E > 0.0 0.0 0.0 0.0
> Eta < 2.5 2.5 2.5 -1.0
> xqcut: 0.0 0.0 0.0 0.0
> d R # 3 > -0.0 0.4 0.4 0.0
> d R # 4 > -0.0 -0.0 0.4 0.0
> d R # 5 > -0.0 -0.0 -0.0 0.0
> s min # 3> 0.0 0.0 0.0 0.0
> s min # 4> 0.0 0.0 0.0 0.0
> s min # 5> 0.0 0.0 0.0 0.0
> xqcutij # 3> 0.0 0.0 0.0 0.0
> xqcutij # 4> 0.0 0.0 0.0 0.0
> xqcutij # 5> 0.0 0.0 0.0 0.0
> No cut BW -3 1
> alpha_s for scale 212.62860824938167 is 0.10467187925059432
> muf2 too small. setting muf2 to muf2min: 2052.7133541613543 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2052.7133541613543 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4039.7705598206871 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4039.7705598206871 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6266.8614765283519 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6266.8614765283519 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1666.6154806528189 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1666.6154806528189 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1837.3016859957845 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1837.3016859957845 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4685.1533484090960 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4685.1533484090960 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3574.9355971828295 6467.2159543705357
> ...........
> muf2 too small. setting muf2 to muf2min: 3574.9355971828295 6467.2159543705357
> Iteration 4 Mean: NaN Abs mean: NaN Fluctuation: NaN Infinity 6.5%
> 4 NaN NaN +- NaN NaN
> Found 4405 events.
> Wrote 0 events.
> Actual xsec NaN
> Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
> STOP 1
>
> ls status:
> events.lhe
> input_app.txt
> run1_app.log
>
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/Muon_BSM/TypeII/Xsec_decay/HppHm/scan-vbf-3mu-1.5TeV/run_01_tag_1_debug.log'.
> Please attach this file to your report.
> INFO:
> INFO:
> command not executed:
> command not executed: The results:
>
>
>
>
>
>
>
> __________________________________________________________________________
>
> #######
> Here is the information in the "run_01_tag_1_debug.log":
> #######
> __________________________________________________________________________
>
>
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.4.0.RC 2022-03-18 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #************************************************************
> generate_events run_01
> Traceback (most recent call last):
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/madevent/sum_html.py", line 307, in read_results
> self.axsec, self.xerru, self.xerrc, self.nevents, self.nw,\
> ValueError: not enough values to unpack (expected 10, got 5)
>
> During handling of the above exception, another exception occurred:
>
> Traceback (most recent call last):
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/extended_cmd.py", line 1544, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
> return func(arg, **opt)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/madevent_interface.py", line 2404, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/common_run_interface.py", line 7549, in new_fct
> original_fct(obj, *args, **opts)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/madevent_interface.py", line 2600, in run_generate_events
> self.exec_cmd('refine %s' % nb_event, postcmd=False)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
> return func(arg, **opt)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/interface/madevent_interface.py", line 3599, in do_refine
> x_improve.launch() # create the ajob for the refinment.
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/madevent/gen_ximprove.py", line 1048, in launch
> self.results = sum_html.collect_result(self.cmd,
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/madevent/sum_html.py", line 753, in collect_result
> P_comb.add_results(os.path.basename(G), path, mfactors[G])
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/madevent/sum_html.py", line 434, in add_results
> oneresult.read_results(filepath)
> File "/Users/maxwell/Software/MG5_aMC_v3_4_0_RC/madgraph/madevent/sum_html.py", line 316, in read_results
> raise Exception("Reported error: End code %s \n Full associated log: \n%s"\
> Exception: Reported error: End code 1.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.11799999999999999
> New value of alpha_s from PDF eva : 0.11799999999999999
> Define smin to 900.00000000000000
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MZ = 91.188000000000002
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MH = 125.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WT = 1.4915000000000000
> mdl_WH = 6.3823389999999999E-003
> Internal Params
> ---------------------------------
>
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM1x1 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_complexi = (0.0000000000000000,1.0000000000000000)
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_I1x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_I2x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I3x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I4x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_cw__exp__2 = 0.77775351421422245
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> Couplings of sm
> ---------------------------------
>
> GC_3 = -0.00000E+00 -0.30795E+00
> GC_5 = 0.00000E+00 0.94836E-01
> GC_36 = 0.00000E+00 0.33188E+00
> GC_50 = -0.00000E+00 -0.28804E+00
> GC_53 = 0.00000E+00 0.57609E+00
> GC_57 = -0.00000E+00 -0.35482E+00
> GC_59 = 0.00000E+00 0.82310E-01
> GC_62 = 0.00000E+00 0.37035E+00
> GC_72 = 0.00000E+00 0.52532E+02
> GC_81 = 0.00000E+00 0.67544E+02
> GC_100 = 0.00000E+00 0.46191E+00
>
> Collider parameters:
> --------------------
>
> Running at m+m- machine @ 3000.0000000000000 GeV
> PDF set = eva
> alpha_s(Mz)= 0.1180 running at 2 loops.
> alpha_s(Mz)= 0.1180 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 4000 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)?
> Monte-Carlo over helicities
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Enter Configuration Number:
> Running Configuration Number: 39
> Not subdividing B.W.
> Attempting mappinvarients 1 6
> Determine nb_t
> T-channel found: 2
> Completed mapping 6
> about to integrate 10 4000 7 3 10 1
> Using non-zero grid deformation.
> 10 dimensions 4000 events 10 invarients 7 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> Error opening grid
> Using Uniform Grid! 24
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 3000.00
> Mapping Graph 39 to config 39
> Determine nb_t
> T-channel found: 2
> Setting grid 1 0.14222E-05 1
> Warning: No cutoff for shat integral found
> Minimum set to 1.0000000000000000E-004
> Setting grid 2 0.10000E-03 1
> Setting grid 3 0.11111E-04 1
> Transforming s_hat 1/s 9 1.0000000000000000E-004 900.00000000000000 9000000.0000000000
> 39 1 2 3 4 5 6 7 8 9 10
> Masses: 0.804E+02 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
> Using random seed offsets 39 : 2
> with seed 24
> Ranmar initialization seeds 31097 9399
> Particle 3 4 5 6
> Et > 10.0 10.0 10.0 0.0
> E > 0.0 0.0 0.0 0.0
> Eta < 2.5 2.5 2.5 -1.0
> xqcut: 0.0 0.0 0.0 0.0
> d R # 3 > -0.0 0.4 0.4 0.0
> d R # 4 > -0.0 -0.0 0.4 0.0
> d R # 5 > -0.0 -0.0 -0.0 0.0
> s min # 3> 0.0 0.0 0.0 0.0
> s min # 4> 0.0 0.0 0.0 0.0
> s min # 5> 0.0 0.0 0.0 0.0
> xqcutij # 3> 0.0 0.0 0.0 0.0
> xqcutij # 4> 0.0 0.0 0.0 0.0
> xqcutij # 5> 0.0 0.0 0.0 0.0
> No cut BW -3 1
> alpha_s for scale 212.62860824938167 is 0.10467187925059432
> muf2 too small. setting muf2 to muf2min: 2052.7133541613543 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2052.7133541613543 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4039.7705598206871 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4039.7705598206871 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6266.8614765283519 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6266.8614765283519 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1666.6154806528189 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1666.6154806528189 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1837.3016859957845 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1837.3016859957845 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4685.1533484090960 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4685.1533484090960 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3574.9355971828295 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3574.9355971828295 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6336.4158102708088 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6336.4158102708088 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3278.9611496313842 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3278.9611496313842 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2824.7267118068685 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2824.7267118068685 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6389.3112666897823 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6389.3112666897823 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4222.1433042355884 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4222.1433042355884 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4760.6900585629837 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4760.6900585629837 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2084.7883692021355 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2084.7883692021355 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2415.1621394769891 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2415.1621394769891 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6277.7629267634420 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 6277.7629267634420 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2781.3549254676818 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2781.3549254676818 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2037.9542723782772 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3089.2014783107848 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3089.2014783107848 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1693.1991548463875 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 1693.1991548463875 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2164.3903225686990 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2164.3903225686990 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 5693.5847921116656 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 5693.5847921116656 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4138.3728716798805 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4138.3728716798805 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4011.6631332894258 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4011.6631332894258 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2569.4239366162415 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2569.4239366162415 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2050.1681578268335 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2050.1681578268335 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2171.5833553771072 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2171.5833553771072 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4166.3628817347499 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 4166.3628817347499 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3088.0920517023742 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3088.0920517023742 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2489.4790735472507 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2489.4790735472507 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3976.4949294648441 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3976.4949294648441 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 5592.4656358546290 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 5592.4656358546290 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 2813.4528388238664 6467.2159543705357
> muf2 too small. setting muf2 to muf2min: 3588.9320306506665 6467.2159543705357
> Iteration 4 Mean: NaN Abs mean: NaN Fluctuation: NaN Infinity 6.5%
> 4 NaN NaN +- NaN NaN
> Found 4405 events.
> Wrote 0 events.
> Actual xsec NaN
> Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
> STOP 1
>
> ls status:
> events.lhe
> input_app.txt
> run1_app.log
>
> Run Options
> -----------
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> run_mode : 2
> cluster_queue : None (user set)
> cluster_time : None (user set)
> cluster_size : 100
> cluster_memory : 100 (user set)
> nb_core : 20 (user set)
> cluster_temp_path : None
>
> Configuration Options
> ---------------------
> pythia8_path : None (user set)
> hwpp_path : None (user set)
> thepeg_path : None (user set)
> hepmc_path : None (user set)
> madanalysis_path : None (user set)
> madanalysis5_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> exrootanalysis_path : None (user set)
> syscalc_path : None (user set)
> lhapdf : lhapdf-config
> lhapdf_py2 : None
> lhapdf_py3 : None
> timeout : 60
> f2py_compiler : None
> f2py_compiler_py2 : None
> f2py_compiler_py3 : None
> web_browser : None
> eps_viewer : None
> text_editor : None
> fortran_compiler : None
> cpp_compiler : None
> auto_update : 7 (user set)
> cluster_type : condor
> cluster_status_update : (600, 30)
> cluster_nb_retry : 1
> cluster_local_path : None
> cluster_retry_wait : 300
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.4.0.RC 2022-03-18 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set ignore_six_quark_processes False
> set low_mem_multicore_nlo_generation False
> set complex_mass_scheme False
> set include_lepton_initiated_processes False
> set gauge unitary
> set loop_optimized_output True
> set loop_color_flows False
> set max_npoint_for_channel 0
> set default_unset_couplings 99
> set max_t_for_channel 99
> set zerowidth_tchannel True
> set nlo_mixed_expansion True
> set group_subprocesses false
> define vap = w+ z a
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> define vap = w+ z a
> define vam = w- z a
> generate vap vam > mu+ mu+ mu- vm
> output Muon_BSM/TypeII/Xsec_decay/HppHm/scan-vbf-3mu-1.5TeV
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.700000e+00 # MB
> 6 1.730000e+02 # MT
> 15 1.777000e+00 # MTA
> 23 9.118800e+01 # MZ
> 25 1.250000e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 1 0.000000e+00 # d : 0.0
> 2 0.000000e+00 # u : 0.0
> 3 0.000000e+00 # s : 0.0
> 4 0.000000e+00 # c : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 0.0
> 14 0.000000e+00 # vm : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 1 1.325070e+02 # aEWM1
> 2 1.166390e-05 # Gf
> 3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.491500e+00 # WT
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 6.382339e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> DECAY 1 0.000000e+00 # d : 0.0
> DECAY 2 0.000000e+00 # u : 0.0
> DECAY 3 0.000000e+00 # s : 0.0
> DECAY 4 0.000000e+00 # c : 0.0
> DECAY 5 0.000000e+00 # b : 0.0
> DECAY 11 0.000000e+00 # e- : 0.0
> DECAY 12 0.000000e+00 # ve : 0.0
> DECAY 13 0.000000e+00 # mu- : 0.0
> DECAY 14 0.000000e+00 # vm : 0.0
> DECAY 15 0.000000e+00 # ta- : 0.0
> DECAY 16 0.000000e+00 # vt : 0.0
> DECAY 21 0.000000e+00 # g : 0.0
> DECAY 22 0.000000e+00 # a : 0.0
> #*********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat MadEvent *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a '# ' are info or comments *
> # *
> # mind the format: value = variable ! comment *
> # *
> # To display more options, you can type the command: *
> # update to_full *
> #*********************************************************************
> #
> #*********************************************************************
> # Tag name for the run (one word) *
> #*********************************************************************
> tag_1 = run_tag ! name of the run
> #*********************************************************************
> # Number of events and rnd seed *
> # Warning: Do not generate more than 1M events in a single run *
> #*********************************************************************
> 1000 = nevents ! Number of unweighted events requested
> 0 = iseed ! rnd seed (0=assigned automatically=default))
> #*********************************************************************
> # Collider type and energy *
> # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=elastic photon of proton,*
> # +/-3=PDF of electron/positron beam *
> # +/-4=PDF of muon/antimuon beam *
> #*********************************************************************
> -4 = lpp1 ! beam 1 type
> 4 = lpp2 ! beam 2 type
> 1500.0 = ebeam1 ! beam 1 total energy in GeV
> 1500.0 = ebeam2 ! beam 2 total energy in GeV
> #*********************************************************************
> # Beam polarization from -100 (left-handed) to 100 (right-handed) *
> #*********************************************************************
> 0.0 = polbeam1 ! beam polarization for beam 1
> 0.0 = polbeam2 ! beam polarization for beam 2
>
>
> #*********************************************************************
> # PDF CHOICE: this automatically fixes alpha_s and its evol. *
> # pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
> # iww=Improved Weizsaecker-Williams Approx.[hep-ph/9310350] *
> # eva=Effective W/Z/A Approx. [21yy.zzzzz] *
> # none=No PDF, same as lhapdf with lppx=0 *
> #*********************************************************************
> eva = pdlabel ! PDF set
> 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
> # To see heavy ion options: type "update ion_pdf"
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.188 = scale ! fixed ren scale
> 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> 4 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
> 0.5 = scalefact ! scale factor for event-by-event scales
>
> 0 = ievo_eva ! scale evolution for EW pdfs (eva):
> ! 0 for evo by q^2; 1 for evo by pT^2
>
> #*********************************************************************
> # Type and output format
> #*********************************************************************
> False = gridpack !True = setting up the grid pack
> -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
> average = event_norm ! average/sum. Normalization of the weight in the LHEF
> # To see MLM/CKKW merging options: type "update MLM" or "update CKKW"
>
> #*********************************************************************
> #
> #*********************************************************************
> # Phase-Space Optimization strategy (basic options)
> #*********************************************************************
> 1 = nhel ! using helicities importance sampling or not.
> ! 0: sum over helicity, 1: importance sampling
> 2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
> ! 1 is old strategy (using amp square)
> ! 2 is new strategy (using only the denominator)
> # To see advanced option for Phase-Space optimization: type "update psoptim"
> #*********************************************************************
> # Generation bias, check the wiki page below for more information: *
> # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
> #*********************************************************************
> None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
> {} = bias_parameters ! Specifies the parameters of the module.
> #
> #*******************************
> # Parton level cuts definition *
> #*******************************
> 0.0 = dsqrt_shat ! minimal shat for full process
> #
> #
> #*********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
> #*********************************************************************
> 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 10.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptlmax ! maximum pt for the charged leptons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
> #
> # For display option for energy cut in the partonic center of mass frame type 'update ecut'
> #
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> 2.5 = etal ! max rap for the charged leptons
> 0.0 = etalmin ! main rap for the charged leptons
> {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.4 = drll ! min distance between leptons
> -1.0 = drllmax ! max distance between leptons
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> #*********************************************************************
> 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
> {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
> ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons / neutrino *
> # for pair of lepton includes only same flavor, opposite charge
> #*********************************************************************
> 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0.0 = xptl ! minimum pt for at least one charged lepton
> #*********************************************************************
> # Control the pt's of leptons sorted by pt *
> #*********************************************************************
> 0.0 = ptl1min ! minimum pt for the leading lepton in pt
> 0.0 = ptl2min ! minimum pt for the second lepton in pt
> 0.0 = ptl3min ! minimum pt for the third lepton in pt
> -1.0 = ptl1max ! maximum pt for the leading lepton in pt
> -1.0 = ptl2max ! maximum pt for the second lepton in pt
> -1.0 = ptl3max ! maximum pt for the third lepton in pt
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 4 = maxjetflavor ! Maximum jet pdg code
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: Do not use for interference type of computation *
> #*********************************************************************
> True = use_syst ! Enable systematics studies
> #
> systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
> ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
>
>
>
>
>
>
> I don't know why is it.
> Look forward to your reply, thanks.
> Yuan
>
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thanks so much, that solved my question.