MadGraph5_aMC@NLO

4j at NLO fails

Asked by Christian Weiss

Dear all,

I am trying to reproduce the total cross section for ee -> 4j as given in table 10 of 1405.0301, also using the parameter and run card provided at the URL given there. I have created the process by

define jbb = j b b~
generate e+ e- > jbb jbb jbb jbb [QCD]

Runnig at LO, I can reproduce the correct result. At NLO, however, the integration starts but then stops without any notification at

INFO: Idle: 86, Running: 0, Completed: 4 [ 3.8s ]

Before that, I get a message

/home/cweiss/Software/Madgraph/MG5_aMC_v2_5_2/4j/SubProcesses/P0_epem_uuxgg/ajob1: line 34: 27839 Segmentation fault (core dumped) ../madevent_mintFO > log.txt < input_app.txt 2>&1
WARNING: program /home/cweiss/Software/Madgraph/MG5_aMC_v2_5_2/4j/SubProcesses/P0_epem_uuxgg/ajob1 4 all 0 0 launch ends with non zero status: 139. Stop all computation

but the program does not stop itself. Do you have an idea what I do wrong? I can provide the file `run_02_tag_1_debug.log` if you want to.

Best,
Christian

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:

This question was reopened

Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi Christian,
did you import the model with a massless b quark before generating the process?
you can do so by typing
import model loop_sm-no_b_mass
Also, can you copy-paste the content of the log.txt inside
/home/cweiss/Software/Madgraph/MG5_aMC_v2_5_2/4j/SubProcesses/P0_epem_uuxgg/all_G4 ?

Let me know.

Best,

Marco

Revision history for this message
Christian Weiss (christian-weiss-g) said : 		#2

Hi Marco,

yes, that apparently solves the problem now. The integration is running.

Thanks,
Christian

Revision history for this message
Christian Weiss (christian-weiss-g) said : 		#3

I'm afraid I was a bit too hasty. The grid setup works fine giving

INFO:
      Results after grid setup:
      Total cross section: 1.131e-01 +- 2.3e-03 pb

However, in the actual integration which follows, I seem to get the same segfault as before:

INFO: Idle: 56, Running: 4, Completed: 0 [ current time: 22h53 ]
INFO: Idle: 55, Running: 4, Completed: 1 [ 3m 53s ]
/home/cweiss/Software/Madgraph/MG5_aMC_v2_5_2/4j_nomb/SubProcesses/P0_epem_uuxgg/ajob1: line 34: 2939 Segmentation fault (core dumped) ../madevent_mintFO > log.txt < input_app.txt 2>&1

This is the content of the log.txt:

 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 Process in group number 0
 No PDF is used, alpha_s(MZ) from param_card is used
 Value of alpha_s from param_card: 0.11799999999999999
 The default order of alpha_s running is fixed to 2
 *****************************************************
 * 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 = 91.188000000000002
 aEWM1 = 132.50700000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymt = 173.19999999999999
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.19999999999999
 mdl_MZ = 91.188000000000002
 mdl_MH = 125.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_WT = 0.0000000000000000
 mdl_WZ = 0.0000000000000000
 mdl_WW = 0.0000000000000000
 mdl_WH = 0.0000000000000000
  Internal Params
  ---------------------------------

 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_CKM22 = 1.0000000000000000
 mdl_I4x33 = 0.0000000000000000
 mdl_I1x33 = 0.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM22 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 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_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_Ncol__exp__2 = 9.0000000000000000
 mdl_MT__exp__2 = 29998.239999999994
 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_v = 246.21845810181637
 mdl_v__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yt = 0.99481489280438762
 mdl_ytau = 1.0206617000654717E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.5430507588273299E-002
 mdl_VectorZDown = -0.13030376310310560
 mdl_VectorAUp = 0.20530251149624587
 mdl_VectorADown = -0.10265125574812294
 mdl_VectorWmDxU = 0.23095271737156670
 mdl_AxialWmDxU = -0.23095271737156670
 mdl_VectorWpUxD = 0.23095271737156670
 mdl_AxialWpUxD = -0.23095271737156670
 mdl_I2x33 = ( 0.99481489280438762 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99481489280438762 , 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.99481489280438762 , 0.0000000000000000 )
 mdl_Axial_tbGp = (-0.99481489280438762 , -0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.99481489280438762 , 0.0000000000000000 )
 mdl_Axial_tbGm = (-0.99481489280438762 , 0.0000000000000000 )
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_ee__exp__2 = 9.4835522759998875E-002
 mdl_sw__exp__2 = 0.22224648578577769
 mdl_yt__exp__2 = 0.98965667094540521
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__4 = 2.1987899468922913
 mdl_G__exp__2 = 1.4828317324943823
 mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
 mdl_G_UVg_1EPS_ = -5.1645779033320030E-002
 mdl_G_UVq_1EPS_ = 3.1300472141406080E-003
 mdl_G_UVb_1EPS_ = 3.1300472141406080E-003
 mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
 mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_G__exp__3 = 1.8056676068262196
 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_G_UVt_FIN_ = -4.0159988564249073E-003
 mdl_GWcft_UV_t_FIN_ = 4.0159988564249073E-003
 mdl_tWcft_UV_FIN_ = -9.4438457457542298E-004
  Couplings of loop_sm-no_b_mass
  ---------------------------------

       UV_3Gt = 0.48903E-02 0.00000E+00
      UV_GQQt = 0.00000E+00 -0.48903E-02
   UVWfct_G_2 = 0.40160E-02 0.00000E+00
       R2_DDA = 0.00000E+00 0.25704E-02
       R2_UUA = -0.00000E+00 -0.51409E-02
    R2_DDZ_V2 = 0.00000E+00 0.72127E-02
    R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
    R2_UUZ_V5 = 0.00000E+00 0.68702E-03
     R2_GGZup = -0.46369E-02 0.00000E+00
   R2_GGZdown = 0.46369E-02 0.00000E+00
 R2_GGGZvecUp = 0.00000E+00 0.28750E-03
 R2_GGGZvecDown -0.00000E+00 -0.49665E-03
 R2_GGGZaxialUp -0.00000E+00 -0.63522E-02
 R2_GGGZaxialDow 0.00000E+00 0.63522E-02
 R2_GGGAvecUp = 0.00000E+00 0.78251E-03
 R2_GGGAvecDown -0.00000E+00 -0.39126E-03
         GC_4 = -0.12177E+01 0.00000E+00
         GC_5 = 0.00000E+00 0.12177E+01
         GC_6 = 0.00000E+00 0.14828E+01
       R2_3Gq = 0.76230E-02 0.00000E+00
       R2_3Gg = 0.31445E-01 0.00000E+00
       R2_GQQ = -0.00000E+00 -0.30492E-01
       R2_GGq = 0.00000E+00 0.62601E-02
       R2_GGt = -0.00000E+00 -0.11268E+04
     R2_GGg_1 = 0.00000E+00 0.28170E-01
     R2_GGg_2 = -0.00000E+00 -0.18780E-01
       R2_QQq = 0.00000E+00 0.12520E-01
  UV_3Gg_1eps = 0.62890E-01 0.00000E+00
  UV_3Gb_1eps = -0.38115E-02 0.00000E+00
 UV_GQQg_1eps = 0.00000E+00 -0.62890E-01
 UV_GQQq_1eps = 0.00000E+00 0.38115E-02
 UVWfct_G_2_1eps -0.31300E-02 0.00000E+00
         GC_1 = -0.00000E+00 -0.10265E+00
         GC_2 = 0.00000E+00 0.20530E+00
        GC_21 = -0.00000E+00 -0.28804E+00
        GC_22 = 0.00000E+00 0.28804E+00
        GC_23 = -0.00000E+00 -0.27437E-01
        GC_24 = 0.00000E+00 0.82310E-01
      GC_1003 = -0.00000E+00 -0.30795E+00

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

 Running at e e machine @ 1000.0000000000000 GeV
 PDF set = none
 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

 Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 Diagram information for clustering has been set-up for nFKSprocess 5
 Diagram information for clustering has been set-up for nFKSprocess 6
 getting user params
 Number of phase-space points per iteration: 175542
 Maximum number of iterations is: 7
 Desired accuracy is: 5.3061136774599999E-003
 Using adaptive grids: 2
 Using Multi-channel integration
 Do MC over helicities for the virtuals
Running Configuration Number: 4
 Splitting channel: 0
 doing the all of this channel
 Normal integration (Sfunction != 1)
 RESTART: Use old grids, but refil plots
 Not subdividing B.W.
 about to integrate 11 175542 7 4
 Update iterations and points to 7 175542
 imode is -1
#--------------------------------------------------------------------------
# FastJet release 3.1.3 [fjcore]
# M. Cacciari, G.P. Salam and G. Soyez
# A software package for jet finding and analysis at colliders
# http://fastjet.fr
#
# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
# for scientific work and optionally PLB641(2006)57 [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.
#--------------------------------------------------------------------------

  +----------------------------------------------------------------+
  | |
  | Ninja - version 1.1.0 |
  | |
  | Author: Tiziano Peraro |
  | |
  | Based on: |
  | |
  | P. Mastrolia, E. Mirabella and T. Peraro, |
  | "Integrand reduction of one-loop scattering amplitudes |
  | through Laurent series expansion," |
  | JHEP 1206 (2012) 095 [arXiv:1203.0291 [hep-ph]]. |
  | |
  | T. Peraro, |
  | "Ninja: Automated Integrand Reduction via Laurent |
  | Expansion for One-Loop Amplitudes," |
  | Comput.Phys.Commun. 185 (2014) [arXiv:1403.1229 [hep-ph]] |
  | |
  +----------------------------------------------------------------+

 ------- iteration 1
 Update # PS points (even): 175542 --> 174080
Using random seed offsets: 4 , 1 , 0
  with seed 35
 Ranmar initialization seeds 14390 9410
 Total number of FKS directories is 6
 For the Born we use nFKSprocesses # 0 3
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 5. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 5. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 6. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
nFKSprocess: 6. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
 bpower is 2.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.474445D+03 0.474445D+03 1.00
 muF1, muF1_reference: 0.474445D+03 0.474445D+03 1.00
 muF2, muF2_reference: 0.474445D+03 0.474445D+03 1.00
 QES, QES_reference: 0.474445D+03 0.474445D+03 1.00

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

 alpha_s= 9.4588782287462023E-002
 alpha_s value used for the virtuals is (for the first PS point): 0.10561757838397419
  ==========================================================================================
 { }
 {   }
 {  ,,  }
 { `7MMM. ,MMF' `7MM `7MMF'  }
 {  MMMb dPMM MM MM  }
 {  M YM ,M MM ,6"Yb. ,M""bMM MM ,pW"Wq. ,pW"Wq.`7MMpdMAo.  }
 {  M Mb M' MM 8) MM ,AP MM MM 6W' `Wb 6W' `Wb MM `Wb  }
 {  M YM.P' MM ,pm9MM 8MI MM MM , 8M M8 8M M8 MM M8  }
 {  M `YM' MM 8M MM `Mb MM MM ,M YA. ,A9 YA. ,A9 MM ,AP  }
 { .JML. `' .JMML.`Moo9^Yo.`Wbmd"MML..JMMmmmmMMM `Ybmd9' `Ybmd9' MMbmmd'  }
 {  MM  }
 {  .JMML.  }
 { v2.5.2 (2016-12-10), Ref: arXiv:1103.0621v2, arXiv:1405.0301  }
 {   }
 { }
  ==========================================================================================
 ===============================================================
 INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
 ===============================================================
  > MLReductionLib = 6|7|1
  > CTModeRun = -1
  > MLStabThres = 1.0000000000000000E-003
  > NRotations_DP = 0
  > NRotations_QP = 0
  > CTStabThres = 1.0000000000000000E-002
  > CTLoopLibrary = 2
  > CTModeInit = 1
  > CheckCycle = 3
  > MaxAttempts = 10
  > UseLoopFilter = F
  > HelicityFilterLevel = 2
  > ImprovePSPoint = 2
  > DoubleCheckHelicityFilter = T
  > LoopInitStartOver = F
  > HelInitStartOver = F
  > ZeroThres = 1.0000000000000001E-009
  > OSThres = 1.0000000000000000E-008
  > WriteOutFilters = T
  > UseQPIntegrandForNinja = T
  > UseQPIntegrandForCutTools = T
  > IREGIMODE = 2
  > IREGIRECY = T
  > COLLIERMode = 1
  > COLLIERRequiredAccuracy = 1.0000000000000000E-008
  > COLLIERCanOutput = F
  > COLLIERComputeUVpoles = T
  > COLLIERComputeIRpoles = T
  > COLLIERGlobalCache = -1
  > COLLIERUseCacheForPoles = F
  > COLLIERUseInternalStabilityTest = T
 ===============================================================

------------------------------------------------------------------------
| You are using CutTools - Version 1.9.3 |
| Authors: G. Ossola, C. Papadopoulos, R. Pittau |
| Published in JHEP 0803:042,2008 |
| http://www.ugr.es/~pittau/CutTools |
| |
| Compiler with 34 significant digits detetected |
  ----------------------------------------------------------------------

########################################################################
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
########################################################################
########################################################################
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
########################################################################
 ---- POLES CANCELLED ----
ABS integral = 0.1690E-01 +/- 0.2688E-03 ( 1.590 %)
Integral = 0.7682E-02 +/- 0.2712E-03 ( 3.531 %)
Virtual = -.6913E-04 +/- 0.1072E-03 ( 155.011 %)
Virtual ratio = 0.1108E+02 +/- 0.1857E+00 ( 1.676 %)
ABS virtual = 0.1801E-02 +/- 0.1071E-03 ( 5.946 %)
Born*ao2pi = 0.1089E-03 +/- 0.3131E-05 ( 2.875 %)
Chi^2 per d.o.f. 0.0000E+00
accumulated results ABS integral = 0.1690E-01 +/- 0.2688E-03 ( 1.590 %)
accumulated results Integral = 0.7682E-02 +/- 0.2712E-03 ( 3.531 %)
accumulated results Virtual = -.6913E-04 +/- 0.1072E-03 ( 155.011 %)
accumulated results Virtual ratio = 0.1108E+02 +/- 0.1857E+00 ( 1.676 %)
accumulated results ABS virtual = 0.1801E-02 +/- 0.1071E-03 ( 5.946 %)
accumulated results Born*ao2pi = 0.1089E-03 +/- 0.3131E-05 ( 2.875 %)
accumulated result Chi^2 per DoF = 0.0000E+00
update virtual fraction to: 0.070 -5.627
  1: 0 1 2 3 45 6
 ------- iteration 2
 Update # PS points (even): 175542 --> 174080
ABS integral = 0.1668E-01 +/- 0.2370E-03 ( 1.421 %)
Integral = 0.7622E-02 +/- 0.2397E-03 ( 3.144 %)
Virtual = 0.1370E-04 +/- 0.9630E-04 ( 702.801 %)
Virtual ratio = 0.1134E+02 +/- 0.2008E+00 ( 1.770 %)
ABS virtual = 0.1679E-02 +/- 0.9621E-04 ( 5.730 %)
Born*ao2pi = 0.1062E-03 +/- 0.3639E-05 ( 3.427 %)
Chi^2= 0.2022E+00
accumulated results ABS integral = 0.1678E-01 +/- 0.1778E-03 ( 1.059 %)
accumulated results Integral = 0.7650E-02 +/- 0.1796E-03 ( 2.348 %)
accumulated results Virtual = -.2550E-04 +/- 0.7162E-04 ( 280.851 %)
accumulated results Virtual ratio = 0.1121E+02 +/- 0.1363E+00 ( 1.216 %)
accumulated results ABS virtual = 0.1737E-02 +/- 0.7156E-04 ( 4.121 %)
accumulated results Born*ao2pi = 0.1076E-03 +/- 0.2373E-05 ( 2.205 %)
accumulated result Chi^2 per DoF = 0.2022E+00
update virtual fraction to: 0.056 -5.927
  1: 0 1 2 3 4 5 6
 ------- iteration 3
 Update # PS points (even): 175542 --> 174080

         *******************************************
         * C O L L I E R *
         * *
         * Complex One-Loop Library *
         * In Extended Regularizations *
         * *
         * by A.Denner, S.Dittmaier, L.Hofer *
         * *
         * version 1.1 *
         * *
         *******************************************

ABS integral = 0.2739E+03 +/- 0.2739E+03 ( 99.993 %)
Integral = -.2738E+03 +/- 0.2739E+03 ( 100.003 %)
Virtual = -.2739E+03 +/- 0.2739E+03 ( 100.000 %)
Virtual ratio = -.1458E+07 +/- 0.1457E+07 ( 99.991 %)
ABS virtual = 0.2739E+03 +/- 0.2739E+03 ( 99.999 %)
Born*ao2pi = 0.1088E-03 +/- 0.3770E-05 ( 3.464 %)
Chi^2= 0.1000E+01
accumulated results ABS integral = 0.1696E-01 +/- 0.1778E-03 ( 1.048 %)
accumulated results Integral = 0.7470E-02 +/- 0.1796E-03 ( 2.404 %)
accumulated results Virtual = -.9713E-04 +/- 0.7162E-04 ( 73.743 %)
accumulated results Virtual ratio = 0.1107E+02 +/- 0.1363E+00 ( 1.231 %)
accumulated results ABS virtual = 0.1808E-02 +/- 0.7156E-04 ( 3.958 %)
accumulated results Born*ao2pi = 0.1081E-03 +/- 0.2008E-05 ( 1.858 %)
accumulated result Chi^2 per DoF = 0.6011E+00
update virtual fraction to: 0.113 *******
  1: 0 1 2 3 4 5 6

Revision history for this message
marco zaro (marco-zaro) said : 		#4

Hi Christian,
can you try to disable collier?
to do so open
Cards/MadLoopParams.dat and change this lines (should be close to the top of the file)
#MLReductionLib
!6|7|1

to
#MLReductionLib
6|1

(remove also the exclamation mark at the beginning of the line…

Let me know if this helps
Cheers,

Marco
On 02 Mar 2017, at 18:22, Christian Weiss <email address hidden> wrote:

> Question #503117 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/503117
>
> Status: Solved => Open
>
> Christian Weiss is still having a problem:
> I'm afraid I was a bit too hasty. The grid setup works fine giving
>
> INFO:
> Results after grid setup:
> Total cross section: 1.131e-01 +- 2.3e-03 pb
>
> However, in the actual integration which follows, I seem to get the same
> segfault as before:
>
> INFO: Idle: 56, Running: 4, Completed: 0 [ current time: 22h53 ]
> INFO: Idle: 55, Running: 4, Completed: 1 [ 3m 53s ]
> /home/cweiss/Software/Madgraph/MG5_aMC_v2_5_2/4j_nomb/SubProcesses/P0_epem_uuxgg/ajob1: line 34: 2939 Segmentation fault (core dumped) ../madevent_mintFO > log.txt < input_app.txt 2>&1
>
> This is the content of the log.txt:
>
> ===============================================================
> INFO: MadFKS read these parameters from FKS_params.dat
> ===============================================================
>> IRPoleCheckThreshold = 1.0000000000000001E-005
>> PrecisionVirtualAtRunTime = 1.0000000000000000E-003
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000001E-003
> ===============================================================
> Process in group number 0
> No PDF is used, alpha_s(MZ) from param_card is used
> Value of alpha_s from param_card: 0.11799999999999999
> The default order of alpha_s running is fixed to 2
> *****************************************************
> * 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 = 91.188000000000002
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymt = 173.19999999999999
> mdl_ymtau = 1.7769999999999999
> mdl_MT = 173.19999999999999
> mdl_MZ = 91.188000000000002
> mdl_MH = 125.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WT = 0.0000000000000000
> mdl_WZ = 0.0000000000000000
> mdl_WW = 0.0000000000000000
> mdl_WH = 0.0000000000000000
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_CKM22 = 1.0000000000000000
> mdl_I4x33 = 0.0000000000000000
> mdl_I1x33 = 0.0000000000000000
> mdl_lhv = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM22 = 1.0000000000000000
> mdl_conjg__CKM33 = 1.0000000000000000
> 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_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_Ncol__exp__2 = 9.0000000000000000
> mdl_MT__exp__2 = 29998.239999999994
> 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_v = 246.21845810181637
> mdl_v__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yt = 0.99481489280438762
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_AxialZUp = -0.18517701861793787
> mdl_AxialZDown = 0.18517701861793787
> mdl_VectorZUp = 7.5430507588273299E-002
> mdl_VectorZDown = -0.13030376310310560
> mdl_VectorAUp = 0.20530251149624587
> mdl_VectorADown = -0.10265125574812294
> mdl_VectorWmDxU = 0.23095271737156670
> mdl_AxialWmDxU = -0.23095271737156670
> mdl_VectorWpUxD = 0.23095271737156670
> mdl_AxialWpUxD = -0.23095271737156670
> mdl_I2x33 = ( 0.99481489280438762 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99481489280438762 , 0.0000000000000000 )
> mdl_Vector_tbGp = (-0.99481489280438762 , 0.0000000000000000 )
> mdl_Axial_tbGp = (-0.99481489280438762 , -0.0000000000000000 )
> mdl_Vector_tbGm = ( 0.99481489280438762 , 0.0000000000000000 )
> mdl_Axial_tbGm = (-0.99481489280438762 , 0.0000000000000000 )
> mdl_gw__exp__2 = 0.42671326129048615
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_yt__exp__2 = 0.98965667094540521
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__4 = 2.1987899468922913
> mdl_G__exp__2 = 1.4828317324943823
> mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
> mdl_G_UVg_1EPS_ = -5.1645779033320030E-002
> mdl_G_UVq_1EPS_ = 3.1300472141406080E-003
> mdl_G_UVb_1EPS_ = 3.1300472141406080E-003
> mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
> mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_G__exp__3 = 1.8056676068262196
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_G_UVt_FIN_ = -4.0159988564249073E-003
> mdl_GWcft_UV_t_FIN_ = 4.0159988564249073E-003
> mdl_tWcft_UV_FIN_ = -9.4438457457542298E-004
> Couplings of loop_sm-no_b_mass
> ---------------------------------
>
> UV_3Gt = 0.48903E-02 0.00000E+00
> UV_GQQt = 0.00000E+00 -0.48903E-02
> UVWfct_G_2 = 0.40160E-02 0.00000E+00
> R2_DDA = 0.00000E+00 0.25704E-02
> R2_UUA = -0.00000E+00 -0.51409E-02
> R2_DDZ_V2 = 0.00000E+00 0.72127E-02
> R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
> R2_UUZ_V5 = 0.00000E+00 0.68702E-03
> R2_GGZup = -0.46369E-02 0.00000E+00
> R2_GGZdown = 0.46369E-02 0.00000E+00
> R2_GGGZvecUp = 0.00000E+00 0.28750E-03
> R2_GGGZvecDown -0.00000E+00 -0.49665E-03
> R2_GGGZaxialUp -0.00000E+00 -0.63522E-02
> R2_GGGZaxialDow 0.00000E+00 0.63522E-02
> R2_GGGAvecUp = 0.00000E+00 0.78251E-03
> R2_GGGAvecDown -0.00000E+00 -0.39126E-03
> GC_4 = -0.12177E+01 0.00000E+00
> GC_5 = 0.00000E+00 0.12177E+01
> GC_6 = 0.00000E+00 0.14828E+01
> R2_3Gq = 0.76230E-02 0.00000E+00
> R2_3Gg = 0.31445E-01 0.00000E+00
> R2_GQQ = -0.00000E+00 -0.30492E-01
> R2_GGq = 0.00000E+00 0.62601E-02
> R2_GGt = -0.00000E+00 -0.11268E+04
> R2_GGg_1 = 0.00000E+00 0.28170E-01
> R2_GGg_2 = -0.00000E+00 -0.18780E-01
> R2_QQq = 0.00000E+00 0.12520E-01
> UV_3Gg_1eps = 0.62890E-01 0.00000E+00
> UV_3Gb_1eps = -0.38115E-02 0.00000E+00
> UV_GQQg_1eps = 0.00000E+00 -0.62890E-01
> UV_GQQq_1eps = 0.00000E+00 0.38115E-02
> UVWfct_G_2_1eps -0.31300E-02 0.00000E+00
> GC_1 = -0.00000E+00 -0.10265E+00
> GC_2 = 0.00000E+00 0.20530E+00
> GC_21 = -0.00000E+00 -0.28804E+00
> GC_22 = 0.00000E+00 0.28804E+00
> GC_23 = -0.00000E+00 -0.27437E-01
> GC_24 = 0.00000E+00 0.82310E-01
> GC_1003 = -0.00000E+00 -0.30795E+00
>
> Collider parameters:
> --------------------
>
> Running at e e machine @ 1000.0000000000000 GeV
> PDF set = none
> 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
>
>
> Diagram information for clustering has been set-up for nFKSprocess 1
> Diagram information for clustering has been set-up for nFKSprocess 2
> Diagram information for clustering has been set-up for nFKSprocess 3
> Diagram information for clustering has been set-up for nFKSprocess 4
> Diagram information for clustering has been set-up for nFKSprocess 5
> Diagram information for clustering has been set-up for nFKSprocess 6
> getting user params
> Number of phase-space points per iteration: 175542
> Maximum number of iterations is: 7
> Desired accuracy is: 5.3061136774599999E-003
> Using adaptive grids: 2
> Using Multi-channel integration
> Do MC over helicities for the virtuals
> Running Configuration Number: 4
> Splitting channel: 0
> doing the all of this channel
> Normal integration (Sfunction != 1)
> RESTART: Use old grids, but refil plots
> Not subdividing B.W.
> about to integrate 11 175542 7 4
> Update iterations and points to 7 175542
> imode is -1
> #--------------------------------------------------------------------------
> # FastJet release 3.1.3 [fjcore]
> # M. Cacciari, G.P. Salam and G. Soyez
> # A software package for jet finding and analysis at colliders
> # http://fastjet.fr
> #
> # Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
> # for scientific work and optionally PLB641(2006)57 [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.
> #--------------------------------------------------------------------------
>
> +----------------------------------------------------------------+
> | |
> | Ninja - version 1.1.0 |
> | |
> | Author: Tiziano Peraro |
> | |
> | Based on: |
> | |
> | P. Mastrolia, E. Mirabella and T. Peraro, |
> | "Integrand reduction of one-loop scattering amplitudes |
> | through Laurent series expansion," |
> | JHEP 1206 (2012) 095 [arXiv:1203.0291 [hep-ph]]. |
> | |
> | T. Peraro, |
> | "Ninja: Automated Integrand Reduction via Laurent |
> | Expansion for One-Loop Amplitudes," |
> | Comput.Phys.Commun. 185 (2014) [arXiv:1403.1229 [hep-ph]] |
> | |
> +----------------------------------------------------------------+
>
>
> ------- iteration 1
> Update # PS points (even): 175542 --> 174080
> Using random seed offsets: 4 , 1 , 0
> with seed 35
> Ranmar initialization seeds 14390 9410
> Total number of FKS directories is 6
> For the Born we use nFKSprocesses # 0 3
> nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 5. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 5. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 6. Absolute lower bound for tau at the Born is 0.14400E-01 0.12000E+03 0.10000E+04
> nFKSprocess: 6. Lower bound for tau is (taking resonances into account) 0.14400E-01 0.12000E+03 0.10000E+04
> bpower is 2.0000000000000000
> Scale values (may change event by event):
> muR, muR_reference: 0.474445D+03 0.474445D+03 1.00
> muF1, muF1_reference: 0.474445D+03 0.474445D+03 1.00
> muF2, muF2_reference: 0.474445D+03 0.474445D+03 1.00
> QES, QES_reference: 0.474445D+03 0.474445D+03 1.00
>
> muR_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF1_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF2_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> QES_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
>
> alpha_s= 9.4588782287462023E-002
> alpha_s value used for the virtuals is (for the first PS point): 0.10561757838397419
> ==========================================================================================
> { }
> {   }
> {  ,,  }
> { `7MMM. ,MMF' `7MM `7MMF'  }
> {  MMMb dPMM MM MM  }
> {  M YM ,M MM ,6"Yb. ,M""bMM MM ,pW"Wq. ,pW"Wq.`7MMpdMAo.  }
> {  M Mb M' MM 8) MM ,AP MM MM 6W' `Wb 6W' `Wb MM `Wb  }
> {  M YM.P' MM ,pm9MM 8MI MM MM , 8M M8 8M M8 MM M8  }
> {  M `YM' MM 8M MM `Mb MM MM ,M YA. ,A9 YA. ,A9 MM ,AP  }
> { .JML. `' .JMML.`Moo9^Yo.`Wbmd"MML..JMMmmmmMMM `Ybmd9' `Ybmd9' MMbmmd'  }
> {  MM  }
> {  .JMML.  }
> { v2.5.2 (2016-12-10), Ref: arXiv:1103.0621v2, arXiv:1405.0301  }
> {   }
> { }
> ==========================================================================================
> ===============================================================
> INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
> ===============================================================
>> MLReductionLib = 6|7|1
>> CTModeRun = -1
>> MLStabThres = 1.0000000000000000E-003
>> NRotations_DP = 0
>> NRotations_QP = 0
>> CTStabThres = 1.0000000000000000E-002
>> CTLoopLibrary = 2
>> CTModeInit = 1
>> CheckCycle = 3
>> MaxAttempts = 10
>> UseLoopFilter = F
>> HelicityFilterLevel = 2
>> ImprovePSPoint = 2
>> DoubleCheckHelicityFilter = T
>> LoopInitStartOver = F
>> HelInitStartOver = F
>> ZeroThres = 1.0000000000000001E-009
>> OSThres = 1.0000000000000000E-008
>> WriteOutFilters = T
>> UseQPIntegrandForNinja = T
>> UseQPIntegrandForCutTools = T
>> IREGIMODE = 2
>> IREGIRECY = T
>> COLLIERMode = 1
>> COLLIERRequiredAccuracy = 1.0000000000000000E-008
>> COLLIERCanOutput = F
>> COLLIERComputeUVpoles = T
>> COLLIERComputeIRpoles = T
>> COLLIERGlobalCache = -1
>> COLLIERUseCacheForPoles = F
>> COLLIERUseInternalStabilityTest = T
> ===============================================================
>
> ------------------------------------------------------------------------
> | You are using CutTools - Version 1.9.3 |
> | Authors: G. Ossola, C. Papadopoulos, R. Pittau |
> | Published in JHEP 0803:042,2008 |
> | http://www.ugr.es/~pittau/CutTools |
> | |
> | Compiler with 34 significant digits detetected |
> ----------------------------------------------------------------------
>
> ########################################################################
> # #
> # You are using OneLOop-3.6 #
> # #
> # for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
> # #
> # author: Andreas van Hameren <email address hidden> #
> # date: 18-02-2015 #
> # #
> # Please cite #
> # A. van Hameren, #
> # Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
> # A. van Hameren, C.G. Papadopoulos and R. Pittau, #
> # JHEP 0909:106,2009, arXiv:0903.4665 #
> # in publications with results obtained with the help of this program. #
> # #
> ########################################################################
> ########################################################################
> # #
> # You are using OneLOop-3.6 #
> # #
> # for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
> # #
> # author: Andreas van Hameren <email address hidden> #
> # date: 18-02-2015 #
> # #
> # Please cite #
> # A. van Hameren, #
> # Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
> # A. van Hameren, C.G. Papadopoulos and R. Pittau, #
> # JHEP 0909:106,2009, arXiv:0903.4665 #
> # in publications with results obtained with the help of this program. #
> # #
> ########################################################################
> ---- POLES CANCELLED ----
> ABS integral = 0.1690E-01 +/- 0.2688E-03 ( 1.590 %)
> Integral = 0.7682E-02 +/- 0.2712E-03 ( 3.531 %)
> Virtual = -.6913E-04 +/- 0.1072E-03 ( 155.011 %)
> Virtual ratio = 0.1108E+02 +/- 0.1857E+00 ( 1.676 %)
> ABS virtual = 0.1801E-02 +/- 0.1071E-03 ( 5.946 %)
> Born*ao2pi = 0.1089E-03 +/- 0.3131E-05 ( 2.875 %)
> Chi^2 per d.o.f. 0.0000E+00
> accumulated results ABS integral = 0.1690E-01 +/- 0.2688E-03 ( 1.590 %)
> accumulated results Integral = 0.7682E-02 +/- 0.2712E-03 ( 3.531 %)
> accumulated results Virtual = -.6913E-04 +/- 0.1072E-03 ( 155.011 %)
> accumulated results Virtual ratio = 0.1108E+02 +/- 0.1857E+00 ( 1.676 %)
> accumulated results ABS virtual = 0.1801E-02 +/- 0.1071E-03 ( 5.946 %)
> accumulated results Born*ao2pi = 0.1089E-03 +/- 0.3131E-05 ( 2.875 %)
> accumulated result Chi^2 per DoF = 0.0000E+00
> update virtual fraction to: 0.070 -5.627
> 1: 0 1 2 3 45 6
> ------- iteration 2
> Update # PS points (even): 175542 --> 174080
> ABS integral = 0.1668E-01 +/- 0.2370E-03 ( 1.421 %)
> Integral = 0.7622E-02 +/- 0.2397E-03 ( 3.144 %)
> Virtual = 0.1370E-04 +/- 0.9630E-04 ( 702.801 %)
> Virtual ratio = 0.1134E+02 +/- 0.2008E+00 ( 1.770 %)
> ABS virtual = 0.1679E-02 +/- 0.9621E-04 ( 5.730 %)
> Born*ao2pi = 0.1062E-03 +/- 0.3639E-05 ( 3.427 %)
> Chi^2= 0.2022E+00
> accumulated results ABS integral = 0.1678E-01 +/- 0.1778E-03 ( 1.059 %)
> accumulated results Integral = 0.7650E-02 +/- 0.1796E-03 ( 2.348 %)
> accumulated results Virtual = -.2550E-04 +/- 0.7162E-04 ( 280.851 %)
> accumulated results Virtual ratio = 0.1121E+02 +/- 0.1363E+00 ( 1.216 %)
> accumulated results ABS virtual = 0.1737E-02 +/- 0.7156E-04 ( 4.121 %)
> accumulated results Born*ao2pi = 0.1076E-03 +/- 0.2373E-05 ( 2.205 %)
> accumulated result Chi^2 per DoF = 0.2022E+00
> update virtual fraction to: 0.056 -5.927
> 1: 0 1 2 3 4 5 6
> ------- iteration 3
> Update # PS points (even): 175542 --> 174080
>
> *******************************************
> * C O L L I E R *
> * *
> * Complex One-Loop Library *
> * In Extended Regularizations *
> * *
> * by A.Denner, S.Dittmaier, L.Hofer *
> * *
> * version 1.1 *
> * *
> *******************************************
>
> ABS integral = 0.2739E+03 +/- 0.2739E+03 ( 99.993 %)
> Integral = -.2738E+03 +/- 0.2739E+03 ( 100.003 %)
> Virtual = -.2739E+03 +/- 0.2739E+03 ( 100.000 %)
> Virtual ratio = -.1458E+07 +/- 0.1457E+07 ( 99.991 %)
> ABS virtual = 0.2739E+03 +/- 0.2739E+03 ( 99.999 %)
> Born*ao2pi = 0.1088E-03 +/- 0.3770E-05 ( 3.464 %)
> Chi^2= 0.1000E+01
> accumulated results ABS integral = 0.1696E-01 +/- 0.1778E-03 ( 1.048 %)
> accumulated results Integral = 0.7470E-02 +/- 0.1796E-03 ( 2.404 %)
> accumulated results Virtual = -.9713E-04 +/- 0.7162E-04 ( 73.743 %)
> accumulated results Virtual ratio = 0.1107E+02 +/- 0.1363E+00 ( 1.231 %)
> accumulated results ABS virtual = 0.1808E-02 +/- 0.7156E-04 ( 3.958 %)
> accumulated results Born*ao2pi = 0.1081E-03 +/- 0.2008E-05 ( 1.858 %)
> accumulated result Chi^2 per DoF = 0.6011E+00
> update virtual fraction to: 0.113 *******
> 1: 0 1 2 3 4 5 6
>
> --
> You received this question notification because you are subscribed to
> the question.

Can you help with this problem?

Provide an answer of your own, or ask Christian Weiss for more information if necessary.

To post a message you must log in.