An error occurred during the collection of results (2)

Dear experts,

this question is related to my aMC@NLO bug report a few weeks ago:

https://answers.launchpad.net/madgraph5/+question/228075

The upshot of the error above seemed to be a rarely encountered rounding error. I lowered the threshold as proposed by you, and the problem seemed to have disappeared.

Yesterday, however, I got another error (generating in principle the same process) that to me looks the same!

_____________________________________________________________
_____________________________________________________________

INFO: Idle: 0, Running: 2, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 2, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 3, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 3, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 3, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 3, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 3, Completed: 125 [ 29h 17m ]
INFO: Idle: 0, Running: 2, Completed: 126 [ 29h 17m ]
INFO: Idle: 0, Running: 2, Completed: 126 [ 29h 17m ]
INFO: Idle: 0, Running: 2, Completed: 126 [ 29h 17m ]
INFO: Idle: 0, Running: 1, Completed: 127 [ 29h 18m ]
INFO: Idle: 0, Running: 0, Completed: 128 [ 29h 18m ]
ERROR: not all jobs terminated correctly

ERROR: not all jobs terminated correctly

INFO: Return code of the event collection: None
INFO: Output of the event collection:
N of directories: 128
Determining the number of unweighted events per channel
2 jobs did not terminate correctly
P0_gu_htbxd/GF5
P0_uxg_hbtxdx/GF6
Found 126 correctly terminated jobs
126 128
Integrated abs(cross-section)
P0_uxg_hbtxdx/GF8/res_0: Final result [ABS]: 5.27825985914420846E-004 +/- 1.87335990100897144E-005
2 jobs did not terminate correctly
P0_gu_htbxd/GF5
P0_uxg_hbtxdx/GF6
Found 126 correctly terminated jobs
126 128
Integrated cross-section
P0_uxg_hbtxdx/GF8/res_0: Final result: 9.72692539258686198E-005 +/- 1.91044279895441232E-005

Error detected in "launch auto "
write debug file /portal/ekpcms5/home/bmaier/Programs/MadGraph5_v2_0_0_beta3/aMC_tH_nocuts_43_fixed/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
aMCatNLOError : An error occurred during the collection of results
quit
_______________________________________________________
_______________________________________________________

The log.txt of one of the jobs that didnt terminate correctly is:

_______________________________________________________
_______________________________________________________

Process in group number 0
 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 cteq6_m: 0.11799999999999999
 *****************************************************
 * 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 = 43.750000000000000
 aEWM1 = 132.50700000000001
 Gf = 1.16639000000000003E-005
 aS = 0.11799999999999999
 ymb = 4.7000000000000002
 ymt = -173.00000000000000
 ymtau = 1.7769999999999999
 MT = 173.00000000000000
 MB = 4.7000000000000002
 MZ = 91.188000000000002
 MH = 125.00000000000000
 MTA = 1.7769999999999999
 WT = 0.0000000000000000
 WZ = 2.4414039999999999
 WW = 2.0476000000000001
 WH = 7.41608800000000037E-003
  Internal Params
  ---------------------------------

 CKM3x3 = 1.0000000000000000
 conjg__CKM3x3 = 1.0000000000000000
 CKM22 = 1.0000000000000000
 lhv = 1.0000000000000000
 conjg__CKM22 = 1.0000000000000000
 conjg__CKM33 = 1.0000000000000000
 Ncol = 3.0000000000000000
 CA = 3.0000000000000000
 TF = 0.50000000000000000
 CF = 1.3333333333333333
 complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 MZ__exp__2 = 8315.2513440000002
 MZ__exp__4 = 69143404.913893804
 sqrt__2 = 1.4142135623730951
 MH__exp__2 = 15625.000000000000
 Ncol__exp__2_m_1 = 8.0000000000000000
 MB__exp__2 = 22.090000000000003
 MT__exp__2 = 29929.000000000000
 Ncol__exp__2_m_1_0 = 8.0000000000000000
 aEW = 7.54677111397888348E-003
 MW = 80.419002445756163
 sqrt__aEW = 8.68721538467815552E-002
 ee = 0.30795376724436879
 MW__exp__2 = 6467.2159543705357
 sw2 = 0.22224648578577766
 cw = 0.88190334743339216
 sqrt__sw2 = 0.47143025548407230
 sw = 0.47143025548407230
 g1 = 0.34919219678733299
 gw = 0.65323293034757990
 v = 246.21845810181637
 v__exp__2 = 60623.529110035903
 lam = 0.12886910601690263
 yb = 2.69955542504654901E-002
 yt = -0.99366614581500623
 ytau = 1.02066170006547170E-002
 muH = 88.388347648318430
 AxialZUp = -0.18517701861793787
 AxialZDown = 0.18517701861793787
 VectorZUp = 7.54305075882732989E-002
 VectorZDown = -0.13030376310310560
 VectorAUp = 0.20530251149624587
 VectorADown = -0.10265125574812294
 VectorWmDxU = 0.23095271737156670
 AxialWmDxU = -0.23095271737156670
 VectorWpUxD = 0.23095271737156670
 AxialWpUxD = -0.23095271737156670
 I1x33 = ( 2.69955542504654901E-002, 0.0000000000000000 )
 I2x33 = (-0.99366614581500623 , 0.0000000000000000 )
 I3x33 = (-0.99366614581500623 , 0.0000000000000000 )
 I4x33 = ( 2.69955542504654901E-002, 0.0000000000000000 )
 Vector_tbGp = ( 1.0206617000654716 , 0.0000000000000000 )
 Axial_tbGp = ( 0.96667059156454072 , -0.0000000000000000 )
 Vector_tbGm = ( -1.0206617000654716 , 0.0000000000000000 )
 Axial_tbGm = ( 0.96667059156454072 , -0.0000000000000000 )
 gw__exp__2 = 0.42671326129048615
 cw__exp__2 = 0.77775351421422245
 ee__exp__2 = 9.48355227599988748E-002
 sw__exp__2 = 0.22224648578577769
 yb__exp__2 = 7.28759949289825397E-004
 yt__exp__2 = 0.98737240933884918
  Internal Params evaluated point by point
  ----------------------------------------

 MU_R__exp__2 = 1914.0625000000000
 sqrt__aS = 0.34351128074635334
 G__exp__2 = 1.4828317324943823
 G__exp__3 = 1.8056676068262196
 G__exp__4 = 2.1987899468922913
  Couplings of loop_sm
  ---------------------------------

      UV_GQQb = 0.00000E+00 0.17006E-01
      UV_GQQt = 0.00000E+00 -0.10480E-01
     UV_bMass = 0.00000E+00 0.10231E+01
     UV_tMass = 0.00000E+00 -0.92029E+01
   UVWfct_b_0 = -0.10884E+00 -0.00000E+00
   UVWfct_t_0 = 0.26598E-01 -0.00000E+00
   UVWfct_G_2 = 0.86064E-02 0.00000E+00
   UVWfct_G_1 = -0.13966E-01 0.00000E+00
      R2_sxcW = -0.00000E+00 -0.11566E-01
       R2_bbH = 0.00000E+00 0.95598E-03
       R2_ttH = -0.00000E+00 -0.35188E-01
      R2_GGHb = 0.00000E+00 -0.33698E-02
      R2_GGHt = 0.00000E+00 0.45657E+01
  UV_Htt_1eps = -0.00000E+00 -0.26391E-01
  UV_Hbb_1eps = 0.00000E+00 0.71698E-03
       UV_Htt = 0.00000E+00 0.37377E-01
       UV_Hbb = -0.00000E+00 0.41551E-02
         GC_4 = -0.12177E+01 0.00000E+00
         GC_5 = 0.00000E+00 0.12177E+01
       R2_GQQ = -0.00000E+00 -0.30492E-01
       R2_QQq = 0.00000E+00 0.12520E-01
       R2_QQb = 0.00000E+00 0.11769E+00
       R2_QQt = 0.00000E+00 0.43320E+01
 UV_GQQg_1eps = 0.00000E+00 -0.62890E-01
 UV_GQQq_1eps = 0.00000E+00 0.38115E-02
 UV_bMass_1eps = 0.00000E+00 0.17653E+00
 UV_tMass_1eps = 0.00000E+00 0.64980E+01
 UVWfct_b_0_1eps -0.18780E-01 0.00000E+00
 UVWfct_G_2_1eps -0.31300E-02 0.00000E+00
        GC_31 = 0.00000E+00 0.52532E+02
        GC_33 = -0.00000E+00 -0.19089E-01
        GC_37 = 0.00000E+00 0.70263E+00
        GC_47 = 0.00000E+00 0.46191E+00

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

 Running at P P machine @ 8000.0000000000000 GeV
 PDF set = cteq6_m
 alpha_s(Mz)= 0.1180 running at 2 loops.
 alpha_s(Mz)= 0.1180 running at 2 loops.
 Renormalization scale fixed @ 43.750000000000000
 Factorization scales fixed @ 43.750000000000000 43.750000000000000

 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 5.00000000000000028E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Explicitly summing over helicities
Enter Configuration Number:
Running Configuration Number: 5
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the novB of this channel
 Normal integration (Sfunction != 1)
 Not subdividing B.W.
 about to integrate 13 -1 12 1
 imode is 0
#--------------------------------------------------------------------------
# FastJet release 3.0.1
# M. Cacciari, G.P. Salam and G. Soyez
# A software package for jet finding and analysis at colliders
# http://fastjet.fr
#
# Please cite arXiv:1111.6097 if you use this package for scientific
# work and optionally also Phys. Lett. B641 (2006) [hep-ph/0512210].
#
# FastJet is provided without warranty under the terms of the GNU GPLv2.
# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
# and 3rd party plugin jet algorithms. See COPYING file for details.
#--------------------------------------------------------------------------
 ------- iteration 1
 Update # PS points (even): 1040 --> 1040
Using random seed offsets: 5 , 1 , 0
  with seed 33
 Ranmar initialization seeds 11953 9408
 Total number of FKS directories is 10
 For the Born we use nFKSprocesses # 2 5
 FKS process map (sum= 3 ) :
           1 --> 5 : 1 6 7 8 9
           2 --> 2 : 2 10
           3 --> 1 : 3
           4 --> 1 : 4
           5 --> 1 : 5
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 1. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 2. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 5. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 5. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 6. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 6. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 6. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 7. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 7. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 7. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 8. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 8. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 8. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 9. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 9. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 9. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 10. Absolute lower bound for tau at the Born is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 10. Lower bound for tau is 0.14317E-02 0.30270E+03 0.80000E+04
nFKSprocess: 10. Lower bound for tau is (taking resonances into account) 0.14317E-02 0.30270E+03 0.80000E+04
 bpower is 1.0000000000000000
 =================================================
 From cuts.f
 Jets are defined with the kT algorithm
 with a mimumal pT of 0.0000000000000000 GeV
 and no maximal pseudo-rapidity.
 Charged leptons are required to have at least 0.0000000000000000 GeV of transverse momentum and
 pseudo rapidity of maximum 100.00000000000000 .
 Opposite charged lepton pairs need to be separated by at least 0.0000000000000000
 and have an invariant mass of 0.0000000000000000 GeV
 =================================================
 Scale values (may change event by event):
 muR, muR_reference: 0.437500D+02 0.437500D+02 1.00
 muF1, muF1_reference: 0.437500D+02 0.437500D+02 1.00
 muF2, muF2_reference: 0.437500D+02 0.437500D+02 1.00
 QES, QES_reference: 0.370422D+03 0.370422D+03 1.00

 muR_reference [functional form]:
    fixed
 muF1_reference [functional form]:
    fixed
 muF2_reference [functional form]:
    fixed
 QES_reference [functional form]:
    sum_i mT(i), i=final state

 alpha_s= 0.13274136713136939
 ================================
 process combination map (specified per FKS dir):
  1 map 1 2
  1 inv. map 1 2
  2 map 1 2
  2 inv. map 1 2
  3 map 1 2
  3 inv. map 1 2
  4 map 1 2
  4 inv. map 1 2
  5 map 1 2
  5 inv. map 1 2
  6 map 1 1 1 1 2 2 2 2
  6 inv. map 4 8
  7 map 1 2
  7 inv. map 1 2
  8 map 1 2
  8 inv. map 1 2
  9 map 1 1 2 2
  9 inv. map 2 4
 10 map 1 2
 10 inv. map 1 2
 ================================
 |int|= 0.17330050572292940 +/- 9.64194474107150645E-002 ( 55.637141396961198 %)
  int = 0.10056116102495566 +/- 9.65187245221813955E-002
 Chi^2 per d.o.f. 0.0000000000000000
 accumulated result |int|= 0.17330050572292940 +/- 9.64194474107150645E-002 ( 55.637141396961198 %)
 accumulated result int = 0.10056116102495566 +/- 9.65187245221813955E-002
 accumulated result Chi^2 per DoF = 0.0000000000000000
  1: 01 2 3 45
 ------- iteration 2
 Update # PS points (even): 2080 --> 2080
 |int|= 0.50615727454507264 +/- 0.44841867710478045 ( 88.592755583294363 %)
  int = -0.43389984395007847 +/- 0.44845508890078872
 Chi^2= 0.37323280384913621
 accumulated result |int|= 0.23220582874026977 +/- 9.42649440631464081E-002 ( 40.595425435502314 %)
 accumulated result int = 5.90432938193407447E-003 +/- 9.43580465806350965E-002
 accumulated result Chi^2 per DoF = 0.37323280384913621
  1: 1 3 45
 ------- iteration 3
 Update # PS points (even): 4160 --> 4160
 |int|= 0.48496992829649122 +/- 0.46106504607586379 ( 95.070852680578383 %)
  int = 0.45447206631336418 +/- 0.46107251486712303
 Chi^2= 0.20717079156430676
 accumulated result |int|= 0.27511147786589046 +/- 9.23544960093967832E-002 ( 33.569844750140575 %)
 accumulated result int = 8.21082302001049924E-002 +/- 9.24421067966431642E-002
 accumulated result Chi^2 per DoF = 0.29020179770672150
  1: 0 1 3 45
 ------- iteration 4
 Update # PS points (even): 8320 --> 8320
 |int|= 3.57942243257004854E-004 +/- 1.85045017149119040E-004 ( 51.696892623052463 %)
  int = 2.48618581346270260E-004 +/- 1.85066551537499190E-004
 Large fluctuation ( >30 % ).Not including iteration in results.
 ------- iteration 5
 Update # PS points (even): 16640 --> 16384
 |int|= 1.00968096051904672E-002 +/- 9.66982481674736813E-003 ( 95.771092006889006 %)
  int = -7.44217359815812097E-003 +/- 9.66997175590711931E-003
 Large fluctuation ( >30 % ).Not including iteration in results.
 2nd bad iteration in a row. Resetting grids and starting from scratch...
 ------- iteration 1
 Update # PS points (even): 4160 --> 4160
 error C in xmcsubt_PY6Q
  2.27308778773548327E-005 2.25872662428771262E-005
Time in seconds: 26

__________________________________________________
__________________________________________________

My impression is that I get this error always when generating events with aMC@NLO with Pythia6q subtraction terms at fixed renormalization/hadronization scale. The old generation was at mu=91 GeV, this one is with mu=43.75 GeV. I did not detect this with Herwig6 or with dynamic scale so far.

So maybe, with decreasing mu (with pythia), one has to further lower the threshold?

What do you propose?

Thank you very much!
Benedikt