e+ e- > j j [QCD QED] with eMELA bug with quarkphreco and lepphreco parameters

Asked by Le Van Dung

Dear mg5 team,
I tried to run NLO EW corrections of the process "e+ e- > j j" following this guide https://answers.launchpad.net/mg5amcnlo/+faq/3324:

In mg5, run:
    import model loop_qcd_qed_sm
    generate e+ e- > j j [QCD QED] aW=10 aEW=10
    output ee_jj

exit mg5, the change the following parameter in FKS_params.dat
    #IRPoleCheckThreshold
    1.0d-5
-> #IRPoleCheckThreshold
    -1d0

Launh mg5 again, then run:
    launh ee_jj
    0
    set pdlabel emela
    set lhaid 137002
    set lpp1 -3
    set lpp2 3
    set ebeam 30
    0

The process runs fine. By default, I have the following parameters:
    quarkphreco True

Then I change the following parameter in the 'run_card.dat':
    set quarkphreco False

The program runs fine, then I change the following parameter in the 'run_card.dat':
    set quarkphreco True
Now, the program breaks entirely, I have tried to set everything back to default, it still broke.

The program only works when I "set quarkphreco False", but this is not the default.
And the numbers is not the same as the first time I ran the program with "quarkphreco True". I might have changed the "lepphreco" parameter too. Flipping the two parameters between "True" and "False"

Here's the error message:

    INFO: Compiling source...
    INFO: ...done, continuing with P* directories
    INFO: Compiling directories...
    INFO: Compiling on 16 cores
    INFO: Compiling P0_epem_uux...
    INFO: Compiling P0_epem_ddx...
    INFO: Compiling P0_epem_bbx...
    INFO: Compiling P0_epem_aa...
    Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
    STOP 1
    INFO: P0_epem_ddx done.
    INFO: P0_epem_bbx done.
    INFO: P0_epem_uux done.
    INFO: P0_epem_aa done.
    INFO: Checking test output:
    INFO: P0_epem_uux
    INFO: Result for test_ME:
    INFO: Passed.
    INFO: Result for check_poles:
    INFO: Poles successfully cancel for 20 points over 20 (tolerance=-1.0e+00)
    INFO: P0_epem_ddx
    INFO: Result for test_ME:
    INFO: Passed.
    INFO: Result for check_poles:
    INFO: Poles successfully cancel for 20 points over 20 (tolerance=-1.0e+00)
    INFO: P0_epem_bbx
    INFO: Result for test_ME:
    INFO: Passed.
    INFO: Result for check_poles:
    INFO: Poles successfully cancel for 20 points over 20 (tolerance=-1.0e+00)
    INFO: P0_epem_aa
    INFO: Result for test_ME:
    INFO: Passed.
    INFO: Result for check_poles:
    INFO: Poles successfully cancel for 20 points over 20 (tolerance=-1.0e+00)
    INFO: Starting run
    INFO: Using 16 cores
    INFO: Cleaning previous results
    INFO: Doing fixed order NLO
    INFO: Setting up grids
    WARNING: program /afs/cern.ch/work/c/conex/dung/MG5_aMC_v3_5_6/test/SubProcesses/P0_epem_aa/ajob1 1 all 0 0 launch ends with non zero status: 1. Stop all computation
    INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 12h06 ]
    Command "launch auto " interrupted with error:
    Exception : program /afs/cern.ch/work/c/conex/dung/MG5_aMC_v3_5_6/test/SubProcesses/P0_epem_aa/ajob1 1 all 0 0 launch ends with non zero status: 1. Stop all computation
    Please report this bug on https://bugs.launchpad.net/mg5amcnlo
    More information is found in '/afs/cern.ch/work/c/conex/dung/MG5_aMC_v3_5_6/test/run_11_tag_1_debug.log'.
    Please attach this file to your report.

I guess, there is something to do with the "P0_epem_aa", previously, when the program was running fine, I had this message when the program tried to run "P0_epem_aa":

    INFO: P0_epem_bbx
    INFO: Result for test_ME:
    INFO: Passed.
    INFO: Result for check_poles:
    INFO: Poles successfully cancel for 20 points over 20 (tolerance=-1.0e+00)
    INFO: P0_epem_aa
    INFO: Result for test_ME:
    INFO: Passed.
    INFO: Result for check_poles:
    INFO: Poles successfully cancel for 20 points over 20 (tolerance=-1.0e+00)
    INFO: Starting run
    INFO: Using 16 cores
    INFO: Cleaning previous results
    INFO: Doing fixed order NLO
    WARNING: No integration channels found for contribution P0_epem_aa
    INFO: Setting up grids
    INFO: Idle: 0, Running: 3, Completed: 0 [ current time: 10h54 ]
    INFO: Idle: 0, Running: 2, Completed: 1 [ 3.7s ]
    INFO: Idle: 0, Running: 1, Completed: 2 [ 3.9s ]
    INFO: Idle: 0, Running: 0, Completed: 3 [ 4s ]
    sum of cpu time of last step: 0 second
    INFO:
        Results after grid setup:
        Total cross section: 1.267e+02 +- 1.2e+00 pb

    INFO: Refining results, step 1
    INFO: Idle: 0, Running: 3, Completed: 0 [ current time: 10h54 ]
    INFO: Idle: 0, Running: 2, Completed: 1 [ 5.9s ]
    INFO: Idle: 0, Running: 1, Completed: 2 [ 6.3s ]
    INFO: Idle: 0, Running: 0, Completed: 3 [ 6.4s ]
    sum of cpu time of last step: 0 second
    INFO:
    --------------------------------------------------------------
        Final results and run summary:
        Process e+ e- > j j [QCD QED] aS=10 aEW=10
        Run at e+-e- collider (30.0 + 30.0 GeV)
        Total cross section: 1.279e+02 +- 6.8e-01 pb
    --------------------------------------------------------------
        Scale variation (computed from histogram information):
            Dynamical_scale_choice -1 (envelope of 9 values):
                1.279e+02 pb +0.6% -0.9%

Here is the content in the all_G1/log.txt . Thank you in advance.
Best, Dung.

    LHAPDF 6.5.4 loading ../../../lib/PDFsets/LL_ALPMZ/LL_ALPMZ_0000.dat
    LL_ALPMZ PDF set, member #0, version 1
    Reading INFO from grid
        2508966
    ===============================================================
    INFO: MadFKS read these parameters from FKS_params.dat
    ===============================================================
    > IRPoleCheckThreshold = -1.0000000000000000
    > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
    > SelectedContributionTypes = All
    > VetoedContributionTypes = None
    > QCD_squared_selected = All
    > QED_squared_selected = All
    > SelectedCouplingOrders = All
    > NHelForMCoverHels = 4
    > VirtualFraction = 1.0000000000000000
    > MinVirtualFraction = 5.0000000000000001E-003
    > SeparateFlavourConfigs = F
    > UsePolyVirtual = F
    ===============================================================
    SPLIT TYPE USED: F T
    A PDF is used, so alpha_s(MZ) is going to be modified
    Old value of alpha_s from param_card: 0.11900000000000002

    ======= Evolution parameters =======

    Perturbative order: LO
    Factorisation scheme: Delta
    Renormalisation scheme: alpha(MZ)
    method: magnus
    Evolution scheme: VFNS
    NLMAX = 3, NUMAX = 2, NDMAX = 3
    Mass thresholds:
    me = 0.000510999 GeV
    mu = 0.002160000 GeV
    md = 0.004670000 GeV
    ms = 0.093000000 GeV
    mm = 0.105660000 GeV
    mc = 1.270000000 GeV
    mt = 1.776860000 GeV
    mb = 4.180000000 GeV
    MW = 80.379000000 GeV
    MZ = 91.187600000 GeV
    Alpha FIXED with reference value:
    Alpha(Qref) = 0.007755544
    W effects included in evolution
    New value of alpha_s from PDF emela : 0.11900000000000002
    WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model.
    Hence it will be set to: 5
    *****************************************************
    * MadGraph/MadEvent *
    * -------------------------------- *
    * http://madgraph.hep.uiuc.edu *
    * http://madgraph.phys.ucl.ac.be *
    * http://madgraph.roma2.infn.it *
    * -------------------------------- *
    * *
    * PARAMETER AND COUPLING VALUES *
    * *
    *****************************************************

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

    MU_R = 1000.0000000000000
    aEWM1 = 128.93000000000001
    mdl_MW = 80.418999999999997
    aS = 0.11900000000000002
    mdl_ymt = 173.30000000000001
    mdl_MZ = 91.188000000000002
    mdl_MT = 173.30000000000001
    mdl_MH = 125.00000000000000
    mdl_WZ = 2.4888300000000001
    mdl_WW = 2.0846499999999999
    mdl_WT = 1.3672800000000001
    mdl_WH = 4.1871800000000004E-003
    Internal Params
    ---------------------------------

    mdl_lhv = 1.0000000000000000
    mdl_SCKM22 = 1.0000000000000000
    mdl_epsIR = 1.0000000000000000
    mdl_epsUV = 1.0000000000000000
    mdl_CKM33 = 1.0000000000000000
    mdl_MC__exp__2 = 0.0000000000000000
    mdl_MS__exp__2 = 0.0000000000000000
    mdl_LogSwitch = 1.0000000000000000
    mdl_MB__exp__4 = 0.0000000000000000
    mdl_Ql__exp__4 = 1.0000000000000000
    mdl_vep = 0.0000000000000000
    mdl_CMSParam = 0.0000000000000000
    mdl_Me__exp__2 = 0.0000000000000000
    mdl_conjg__CKM22 = 1.0000000000000000
    mdl_Ql__exp__2 = 1.0000000000000000
    mdl_CKM22 = 1.0000000000000000
    mdl_yb = 0.0000000000000000
    mdl_MU__exp__2 = 0.0000000000000000
    mdl_MB__exp__2 = 0.0000000000000000
    mdl_conjg__CKM11 = 1.0000000000000000
    mdl_SCKM22__exp__2 = 1.0000000000000000
    mdl_conjg__CKM33 = 1.0000000000000000
    mdl_SCKM33__exp__2 = 1.0000000000000000
    mdl_SCKM11 = 1.0000000000000000
    mdl_SCKM33 = 1.0000000000000000
    mdl_ep = 1.0000000000000000E-010
    mdl_Ncol = 3.0000000000000000
    mdl_CA = 3.0000000000000000
    mdl_TF = 0.50000000000000000
    mdl_CF = 1.3333333333333333
    mdl_complexi = (0.0000000000000000,1.0000000000000000)
    mdl_MZ__exp__2 = 8315.2513440000002
    mdl_MW__exp__2 = 6467.2155609999991
    mdl_sqrt__2 = 1.4142135623730951
    mdl_sw2 = 0.22224653309289089
    mdl_cw = 0.88190332061236121
    mdl_sqrt__sw2 = 0.47143030565810135
    mdl_sw = 0.47143030565810135
    mdl_MH__exp__2 = 15625.000000000000
    mdl_I3d = -0.50000000000000000
    mdl_I3u = 0.50000000000000000
    mdl_I3l = -0.50000000000000000
    mdl_I3v = 0.50000000000000000
    mdl_Qd = -0.33333333333333331
    mdl_Qu = 0.66666666666666663
    mdl_Ql = -1.0000000000000000
    mdl_Ncol__exp__2 = 9.0000000000000000
    mdl_sw__exp__2 = 0.22224653309289086
    mdl_MT__exp__2 = 30032.890000000003
    mdl_cw__exp__2 = 0.77775346690710911
    mdl_MT__exp__4 = 901974481.75210023
    mdl_sw__exp__4 = 4.9393521471809435E-002
    mdl_MW__exp__4 = 41824877.112440534
    mdl_MT__exp__3 = 5204699.8370000012
    mdl_MH__exp__4 = 244140625.00000000
    mdl_MZ__exp__4 = 69143404.913893804
    mdl_MT__exp__6 = 27088900393267.840
    mdl_MW__exp__6 = 270490496098.48813
    mdl_MT__exp__8 = 8.1355796573196979E+017
    mdl_MW__exp__8 = 1749320345470752.0
    mdl_MT__exp__10 = 2.4433496893452020E+022
    mdl_MW__exp__10 = 1.1313231759402340E+019
    mdl_MZ__exp__6 = 574944790638.99182
    mdl_sw__exp__3 = 0.10477375102743487
    mdl_MZ__exp__3 = 758251.13955667207
    mdl_cw__exp__4 = 0.60490045528602765
    mdl_cw__exp__6 = 0.47046342623239684
    mdl_MH__exp__6 = 3814697265625.0000
    mdl_cw__exp__3 = 0.68590336508315575
    mdl_Qd__exp__3 = -3.7037037037037028E-002
    mdl_Qd__exp__2 = 0.11111111111111110
    mdl_I3d__exp__2 = 0.25000000000000000
    mdl_I3d__exp__3 = -0.12500000000000000
    mdl_Qu__exp__3 = 0.29629629629629622
    mdl_Qu__exp__2 = 0.44444444444444442
    mdl_I3u__exp__2 = 0.25000000000000000
    mdl_I3u__exp__3 = 0.12500000000000000
    mdl_Ql__exp__3 = -1.0000000000000000
    mdl_I3l__exp__2 = 0.25000000000000000
    mdl_I3l__exp__3 = -0.12500000000000000
    mdl_I3v__exp__3 = 0.12500000000000000
    mdl_MW__exp__3 = 520087.00820005889
    mdl_sw__exp__6 = 1.0977538904358909E-002
    mdl_sw__exp__8 = 2.4397199633860995E-003
    mdl_Qu__exp__4 = 0.19753086419753083
    mdl_Qd__exp__4 = 1.2345679012345677E-002
    mdl_I3l__exp__4 = 6.2500000000000000E-002
    mdl_I3u__exp__4 = 6.2500000000000000E-002
    mdl_I3d__exp__4 = 6.2500000000000000E-002
    mdl_I3v__exp__4 = 6.2500000000000000E-002
    mdl_I3v__exp__2 = 0.25000000000000000
    mdl_aEW = 7.7561467462964398E-003
    mdl_Gf = (1.19874983504616246E-005,0.0000000000000000)
    mdl_sqrt__aEW = 8.8068988561788533E-002
    mdl_ee = 0.31219643584339185
    mdl_g1 = 0.35400301659666522
    mdl_gw = 0.66223242777652103
    mdl_vev = 242.87243157213206
    mdl_vev__exp__2 = 58987.018017759969
    mdl_lam = 0.13244439645428069
    mdl_yt = 1.0091026337275697
    mdl_muH = 88.388347648318430
    mdl_I233 = (1.0091026337275697,0.0000000000000000)
    mdl_I333 = (1.0091026337275697,0.0000000000000000)
    mdl_AxialZUp = -0.18772818184783885
    mdl_AxialZDown = 0.18772818184783885
    mdl_VectorZUp = 7.6469681635801659E-002
    mdl_VectorZDown = -0.13209893174182025
    mdl_AxialG0Up3 = -0.71354331522196934
    mdl_VectorHUp3 = (-0.0000000000000000,-0.71354331522196934)
    mdl_VectorAUp = 0.20813095722892791
    mdl_VectorADown = -0.10406547861446395
    mdl_VectorWmDxU = 0.23413452020120429
    mdl_AxialWmDxU = -0.23413452020120429
    mdl_VectorWpUxD = 0.23413452020120429
    mdl_AxialWpUxD = -0.23413452020120429
    mdl_VectorGpUx3D3 = (-0.50455131686378485,0.0000000000000000)
    mdl_AxialGpUx3D3 = (-0.50455131686378485,0.0000000000000000)
    mdl_VectorGmDx3U3 = (0.50455131686378485,0.0000000000000000)
    mdl_AxialGmDx3U3 = (-0.50455131686378485,0.0000000000000000)
    mdl_ee__exp__2 = 9.7466614553317080E-002
    mdl_R2VV_FIN_ = 9.8754327521535911E-003
    mdl_ee__exp__3 = 3.0428729677267258E-002
    mdl_R2SFF_FIN_ = 3.0830749076334492E-003
    mdl_ee__exp__4 = 9.4997409524848805E-003
    mdl_ZMass2_UV_EW_1EPS_ = (-131.83588536647306,0.0000000000000000)
    mdl_cWcft_UV_EW_R_1EPS_ = (-7.83874748570571384E-005,0.0000000000000000)
    mdl_bWcft_UV_EW_R_1EPS_ = (-1.95968687142642846E-005,0.0000000000000000)
    mdl_eWcft_UV_EW_R_1EPS_ = (-1.76371818428378545E-004,0.0000000000000000)
    mdl_cWcft_UV_EW_L_1EPS_ = (-1.83059772155613233E-003,0.0000000000000000)
    mdl_bWcft_UV_EW_L_1EPS_ = (-5.26052834620677982E-003,0.0000000000000000)
    mdl_dWcft_UV_EW_L_1EPS_ = (-2.03633590389266537E-003,0.0000000000000000)
    mdl_eWcft_UV_EW_L_1EPS_ = (-1.66405327664812821E-003,0.0000000000000000)
    mdl_cWcft_BAR_UV_EW_R_1EPS_ = (-7.83874748570571384E-005,0.0000000000000000)
    mdl_bWcft_BAR_UV_EW_R_1EPS_ = (-1.95968687142642846E-005,0.0000000000000000)
    mdl_eWcft_BAR_UV_EW_R_1EPS_ = (-1.76371818428378545E-004,0.0000000000000000)
    mdl_cWcft_BAR_UV_EW_L_1EPS_ = (-1.83059772155613233E-003,0.0000000000000000)
    mdl_bWcft_BAR_UV_EW_L_1EPS_ = (-5.26052834620677982E-003,0.0000000000000000)
    mdl_dWcft_BAR_UV_EW_L_1EPS_ = (-2.03633590389266537E-003,0.0000000000000000)
    mdl_eWcft_BAR_UV_EW_L_1EPS_ = (-1.66405327664812821E-003,0.0000000000000000)
    mdl_ZZWcft_UV_EW_1EPS_ = (-3.00516333930618204E-003,0.0000000000000000)
    mdl_AZWcft_UV_EW_1EPS_ = (-7.20327294279212577E-003,0.0000000000000000)
    mdl_ZAWcft_UV_EW_1EPS_ = (4.61848593104902116E-003,0.0000000000000000)
    mdl_AAWcft_UV_EW_1EPS_ = (7.54373335233954985E-004,0.0000000000000000)
    mdl_eCoup_UV_EW_1EPS_ = (1.13156000285093193E-003,0.0000000000000000)
    mdl_SWCoup_UV_EW_1EPS_ = (2.60641210157978431E-003,0.0000000000000000)
    Internal Params evaluated point by point
    ----------------------------------------

    mdl_sqrt__aS = 0.34496376621320685
    mdl_G__exp__4 = 2.2362154867812229
    mdl_G__exp__2 = 1.4953981031087418
    mdl_R2MixedFactor_FIN_ = -2.5252584303914066E-002
    mdl_G__exp__3 = 1.8286695702272273
    mdl_MU_R__exp__2 = 1000000.0000000000
    mdl_ZMass2_UV_EW_FIN_ = (-385.12955353367431,0.0000000000000000)
    mdl_cWcft_UV_EW_R_FIN_ = (-3.36255924164604474E-004,-0.0000000000000000)
    mdl_bWcft_UV_EW_R_FIN_ = (-8.40639810411511186E-005,0.0000000000000000)
    mdl_eWcft_UV_EW_R_FIN_ = (-7.56575829370360040E-004,-0.0000000000000000)
    mdl_cWcft_UV_EW_L_FIN_ = (-8.20166283160201660E-003,-0.0000000000000000)
    mdl_bWcft_UV_EW_L_FIN_ = (-2.13238746930339181E-002,0.0000000000000000)
    mdl_dWcft_UV_EW_L_FIN_ = (-9.08421049324912092E-003,-0.0000000000000000)
    mdl_eWcft_UV_EW_L_FIN_ = (-7.48724313203721556E-003,-0.0000000000000000)
    mdl_cWcft_BAR_UV_EW_R_FIN_ = (-3.36255924164604474E-004,-0.0000000000000000)
    mdl_bWcft_BAR_UV_EW_R_FIN_ = (-8.40639810411511186E-005,0.0000000000000000)
    mdl_eWcft_BAR_UV_EW_R_FIN_ = (-7.56575829370360040E-004,-0.0000000000000000)
    mdl_cWcft_BAR_UV_EW_L_FIN_ = (-8.20166283160201660E-003,-0.0000000000000000)
    mdl_bWcft_BAR_UV_EW_L_FIN_ = (-2.13238746930339181E-002,0.0000000000000000)
    mdl_dWcft_BAR_UV_EW_L_FIN_ = (-9.08421049324912092E-003,-0.0000000000000000)
    mdl_eWcft_BAR_UV_EW_L_FIN_ = (-7.48724313203721556E-003,-0.0000000000000000)
    mdl_ZZWcft_UV_EW_FIN_ = (-1.01950909161341863E-002,0.0000000000000000)
    mdl_AZWcft_UV_EW_FIN_ = (-3.65197478696778557E-002,0.0000000000000000)
    mdl_ZAWcft_UV_EW_FIN_ = (2.32818320295721908E-002,-0.0000000000000000)
    mdl_AAWcft_UV_EW_FIN_ = (5.89919023924206833E-003,0.0000000000000000)
    mdl_eCoup_UV_EW_FIN_ = (8.53848558153882074E-003,0.0000000000000000)
    mdl_SWCoup_UV_EW_FIN_ = (2.27007143175545385E-002,0.0000000000000000)
    Couplings of loop_qcd_qed_sm
    ---------------------------------

            GC_11 = 0.00000E+00 0.12229E+01
        R2_DDA = 0.00000E+00 0.26279E-02
        R2_UUA = -0.00000E+00 -0.52558E-02
        R2_DDZ_V2 = 0.00000E+00 0.73741E-02
        R2_DDZ_V3 = 0.00000E+00 0.70239E-03
    c_UVZMass1EW = -0.00000E+00 -0.10195E-01
    c_UVZMass2EW = -0.00000E+00 -0.46990E+03
    c_UVAMass1EW = 0.00000E+00 0.58992E-02
    c_UVAZMass1EW = 0.00000E+00 0.66190E-02
    c_UVAZMass2EW = 0.00000E+00 -0.96797E+02
    c_UVemepMass1EW -0.00000E+00 0.74872E-02
    c_UVemepMass2EW 0.00000E+00 -0.75658E-03
    c_UVepemA1EW = 0.00000E+00 0.11789E-02
    c_UVepemA2EW = 0.00000E+00 -0.52931E-02
    c_UVcxcA1EW = -0.00000E+00 -0.23915E-02
    c_UVcxcA2EW = 0.00000E+00 0.36162E-02
    c_UVdxdA1EW = 0.00000E+00 0.34729E-02
    c_UVbxbA1EW = 0.00000E+00 0.47466E-02
    c_UVbxbA2EW = 0.00000E+00 -0.18343E-02
    c_UVepemZ1EW = 0.00000E+00 0.18389E-01
    c_UVepemZ2EW = 0.00000E+00 0.50798E-02
    c_UVcxcZ1EW = 0.00000E+00 -0.17259E-01
    c_UVcxcZ2EW = 0.00000E+00 -0.34333E-02
    c_UVdxdZ1EW = 0.00000E+00 0.16262E-01
    c_UVbxbZ1EW = 0.00000E+00 0.20176E-01
    c_UVbxbZ2EW = 0.00000E+00 0.17307E-02
            GC_1 = -0.00000E+00 -0.10407E+00
            GC_2 = 0.00000E+00 0.20813E+00
            GC_3 = -0.00000E+00 -0.31220E+00
            GC_5 = 0.00000E+00 0.97467E-01
            GC_6 = 0.00000E+00 0.19493E+00
            GC_30 = -0.10091E+01 -0.00000E+00
            GC_54 = 0.00000E+00 0.34109E+00
            GC_68 = -0.00000E+00 -0.29201E+00
            GC_70 = -0.00000E+00 -0.58402E+00
            GC_75 = -0.00000E+00 -0.36466E+00
            GC_76 = -0.00000E+00 -0.27815E-01
            GC_77 = 0.00000E+00 0.83444E-01
            GC_78 = -0.37546E+00 0.00000E+00
            GC_79 = 0.00000E+00 -0.20857E+00
            GC_80 = 0.00000E+00 0.37546E+00
            GC_81 = 0.00000E+00 0.13023E+00
            GC_82 = 0.00000E+00 0.87002E-01
            GC_83 = 0.00000E+00 0.28193E+00
            GC_84 = -0.13421E+02 0.00000E+00
            GC_92 = -0.25107E+02 0.00000E+00
            GC_99 = 0.00000E+00 0.68474E+02
        GC_124 = 0.00000E+00 0.46827E+00
    R2_AAboson1 = 0.00000E+00 0.61721E-03
    R2_AAboson2 = -0.00000E+00 -0.41148E-03
    R2_AAboson3 = -0.00000E+00 -0.79833E+01
        R2_AAU = 0.00000E+00 0.54864E-03
        R2_AAD = 0.00000E+00 0.13716E-03
        R2_AAt3 = -0.00000E+00 -0.98863E+02
    R2_AZboson1 = 0.00000E+00 0.11546E-02
    R2_AZboson2 = -0.00000E+00 -0.76975E-03
    R2_AZboson3 = -0.00000E+00 -0.14934E+02
        R2_AZl = 0.00000E+00 0.27468E-04
        R2_AZU = 0.00000E+00 0.20157E-03
        R2_AZD = 0.00000E+00 0.17411E-03
        R2_AZt3 = -0.00000E+00 -0.36323E+02
    R2_ZZboson1 = 0.00000E+00 0.21599E-02
    R2_ZZboson2 = -0.00000E+00 -0.14400E-02
    R2_ZZboson3 = -0.00000E+00 -0.27938E+02
        R2_ZZl = 0.00000E+00 0.15061E-03
        R2_ZZU = 0.00000E+00 0.52040E-03
        R2_ZZD = 0.00000E+00 0.66735E-03
        R2_ZZt3 = -0.00000E+00 -0.93775E+02
        R2_ZZv = 0.00000E+00 0.29756E-03
        R2_LLCm = 0.00000E+00 0.79359E-03
        R2_LLCp0 = 0.00000E+00 0.89269E-03
        R2_QQCpcs = 0.00000E+00 0.13886E-02
        R2_bbA2Cp = 0.00000E+00 0.14909E-03
    R2_bbA2Cp_t = 0.00000E+00 -0.38205E-03
        R2_bbA2Cm = 0.00000E+00 0.18352E-04
        R2_ccA2Cp = 0.00000E+00 -0.29818E-03
    R2_ccA2Cp_s = 0.00000E+00 -0.57801E-03
        R2_ccA2Cm = 0.00000E+00 -0.14682E-03
        R2_eeACp = 0.00000E+00 0.55739E-03
        R2_eeACm = 0.00000E+00 0.49551E-03
    R2_ddZ2Cp_u = -0.00000E+00 0.88821E-03
        R2_bbZ2Cp = 0.00000E+00 0.45821E-03
    R2_bbZ2Cp_t = -0.00000E+00 0.12469E-02
        R2_bbZ2Cm = 0.00000E+00 -0.98104E-05
        R2_ccZ2Cp = 0.00000E+00 -0.37851E-03
    R2_ccZ2Cp_s = -0.00000E+00 -0.73372E-03
        R2_ccZ2Cm = 0.00000E+00 0.78483E-04
        R2_eeZCp = 0.00000E+00 0.37237E-03
        R2_llZCp = -0.00000E+00 0.57923E-03
        R2_eeZCm = 0.00000E+00 -0.26488E-03
        R2_AWWlv = 0.00000E+00 -0.28901E-03
        R2_AWWcs = 0.00000E+00 -0.86702E-03
    c_UVZMass1EW_1e -0.00000E+00 -0.30052E-02
    c_UVZMass2EW_1e -0.00000E+00 -0.15682E+03
    c_UVAMass1EW_1e 0.00000E+00 0.75437E-03
    c_UVAZMass1EW_1 0.00000E+00 0.12924E-02
    c_UVAZMass2EW_1 0.00000E+00 -0.19202E+02
    c_UVemepMass1EW 0.00000E+00 0.16641E-02
    c_UVemepMass2EW -0.00000E+00 -0.17637E-03
    c_UVepemA1EW_1e 0.00000E+00 0.53012E-03
    c_UVepemA2EW_1e 0.00000E+00 -0.80135E-03
    c_UVcxcA1EW_1ep -0.00000E+00 -0.67708E-03
    c_UVcxcA2EW_1ep 0.00000E+00 0.55462E-03
    c_UVdxdA1EW_1ep 0.00000E+00 0.79346E-03
    c_UVbxbA1EW_1ep 0.00000E+00 0.11290E-02
    c_UVbxbA2EW_1ep 0.00000E+00 -0.28343E-03
    c_UVepemZ1EW_1e 0.00000E+00 0.19690E-02
    c_UVepemZ2EW_1e -0.00000E+00 -0.29434E-04
    c_UVcxcZ1EW_1ep 0.00000E+00 -0.21056E-02
    c_UVcxcZ2EW_1ep 0.00000E+00 0.87213E-05
    c_UVdxdZ1EW_1ep 0.00000E+00 0.22732E-02
    c_UVbxbZ1EW_1ep 0.00000E+00 0.33044E-02
    c_UVbxbZ2EW_1ep -0.00000E+00 -0.10902E-05

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

    Running at e+e- machine @ 60.000000000000000 GeV
    PDF set = emela
    alpha_s(Mz)= 0.1190 running at 2 loops.
    alpha_s(Mz)= 0.1190 running at 2 loops.
    Renormalization scale set on event-by-event basis
    Factorization scale set on event-by-event basis

    INFO: orders_tag_plot is computed as: + QCD * 1 + QED * 100
    orders_tag_plot= 400 for QCD,QED, = 0 , 4 ,
    AMP_SPLIT: 1 correspond to S.O. 0 4
    orders_tag_plot= 600 for QCD,QED, = 0 , 6 ,
    AMP_SPLIT: 2 correspond to S.O. 0 6
    getting user params
    Number of phase-space points per iteration: -1
    Maximum number of iterations is: 6
    Desired accuracy is: 5.0000000000000003E-002
    Using adaptive grids: 2
    Using Multi-channel integration
    Do MC over helicities for the virtuals
    Number of channels to integrate together: 4
    Running Configuration Number(s): 1 1 2 2
    initial-or-final 1 2 1 2
    Splitting channel: 0
    Weight multiplier: 1.0000000000000000
    doing the all of this channel
    Normal integration (Sfunction != 1)
    RESTART: Fresh run
    about to integrate 7 -1 6
    imode is 0
    channel 1 : 1 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
    channel 2 : 1 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
    channel 3 : 2 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
    channel 4 : 2 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
    Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
    STOP 1
    Thanks for using LHAPDF 6.5.4. Please make sure to cite the paper:
    Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
    ------- iteration 1
    Update # PS points (even_rn): 280 --> 256
    Using random seed offsets: 0 , 4 , 0
    with seed 42
    Ranmar initialization seeds 22919 9420
    initial-final FKS maps:
            0 : 8 1 2 3 4 5 6 7 8
            1 : 4 5 6 7 8 0 0 0 0
            2 : 4 1 2 3 4 0 0 0 0
    Total number of FKS directories is 8
    For the Born we use nFKSprocesses:
            1 2 1 2 2 2 2 2
    tau_min 1 1 : 0.00000E+00 -- 0.20000E+02
    tau_min 2 1 : 0.00000E+00 -- 0.20000E+02
    tau_min 3 1 : 0.00000E+00 -- 0.20000E+02
    tau_min 4 1 : 0.00000E+00 -- 0.20000E+02
    tau_min 5 1 : 0.20000E+02 0.20000E+02 0.20000E+02
    tau_min 6 1 : 0.20000E+02 0.20000E+02 0.20000E+02
    tau_min 7 1 : 0.20000E+02 0.20000E+02 0.20000E+02
    tau_min 8 1 : 0.20000E+02 0.20000E+02 0.20000E+02
    ERROR, too many photons 2
    Time in seconds: 0

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said :
#1

Hi,
thanks for your message.
The problem is that the definition of EW corrections to jets in e+e- collisions is ill-defined in a model with massless partons and leptons (such as those that we employ in mg5_aMC)
The reason is that, you start from a process like

e+ e- > u u~

and you add real radiation

e+ e- > u u~ gamma

IR safety tells you that you must cluster the photon with the quarks in the jet algorithm or, at least, do some kind of recombination.

However, the same real emission has also another underlying Born, namely
e+ e-> gamma gamma

which, in turn, can give you also processes like
e+ e- > gamma mu+ mu-
for which, yet another underlying Born is needed
e+ e- > mu+ mu-

All these happen at the same order in alpha (alpha^2 at the born), in a way similar to the case of the LO_3 at pp colliders (see also 1612.06548)

So, the thing which is most correct, at least formally, is to generate

e+ e- > j j
with j = quarks, photon, leptons, gluon

and to edit cuts.f so that the jet clustering will run on any kind of particles.
However, in this case, you will have also contributions where jets are leptons or photons.

An alternative would be to have a model with massive quarks, in this case one would not have collinear singularities but only photon emissions...

So this is the situation...

Best,

Marco

Can you help with this problem?

Provide an answer of your own, or ask Le Van Dung for more information if necessary.

To post a message you must log in.