Segmentation fault in p p > j j w+ at fixed order

Dear Madgraph experts,

We are studying the following process p p > j j w+ [QCD] using the loop_sm model and the zero_mass_ckm restriction file on MG5_v3_4_0. After launching madevent, a segmentation fault seems to arise from the interface between aMC@NLO and Fixed order such that all jobs crash.

./ajob1 : line 38 : 1355624 Segmentation fault ../madevent_mintFO > log.txt < input_app.txt 2>&1

The error looks a bit similar to the question #675872 but I do not use an external version of FastJet and the jets have non-zero minimum pT.

What is the issue here ? Thank you in advance.

Best,

Julien

Here is the proc_card of the process generated :

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.0 2022-05-06 *
#* *
#* 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 group_subprocesses Auto
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
import model sm
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~
set complex_mass_scheme
import model loop_sm-zeromass_ckm
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate p p > w+ j [QCD]
add process p p > w+ j j [QCD]
output SM_NLO_FIXEDORDER_CMS_jjw+/

For completeness here is the banner of my run:

<LesHouchesEvents version="3.0">
<header>
<!--
#*********************************************************************
# *
# MadGraph5_aMC@NLO *
# *
# Going Beyond *
# *
# http://madgraph.hep.uiuc.edu *
# http://madgraph.phys.ucl.ac.be *
# http://amcatnlo.cern.ch *
# *
# The MadGraph5_aMC@NLO team *
# *
#....................................................................*
# *
# This file contains all the information necessary to reproduce *
# the events generated: *
# *
# 1. software version *
# 2. proc_card : code generation info including model *
# 3. param_card : model primary parameters in the LH format *
# 4. run_card : running parameters (collider and cuts) *
# 5. pythia_card : present only if pythia has been run *
# 6. pgs_card : present only if pgs has been run *
# 7. delphes_cards : present only if delphes has been run *
# *
# *
#*********************************************************************
-->
<MGVersion>
3.4.0
</MGVersion>
<MGRunCard>
<![CDATA[
#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# Some of the values of variables can be list. These can either be *
# comma or space separated. *
# *
# To display additional parameter, you can use the command: *
# update to_full *
#***********************************************************************
#
#*******************
# Running parameters
#*******************
#
#***********************************************************************
# Tag name for the run (one word) *
#***********************************************************************
  tag_1 = run_tag ! name of the run
#***********************************************************************
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#***********************************************************************
 10000 = nevents ! Number of unweighted events requested
 -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
 -1 = nevt_job! Max number of events per job in event generation.
                 ! (-1= no split).
#***********************************************************************
# Output format
#***********************************************************************
  -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/bias. Normalization of the weight in the LHEF
#***********************************************************************
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#***********************************************************************
 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
                     ! number of points and iter. below)
# These numbers are ignored except if req_acc_FO is equal to -1
 5000 = npoints_FO_grid ! number of points to setup grids
 4 = niters_FO_grid ! number of iter. to setup grids
 10000 = npoints_FO ! number of points to compute Xsec
 6 = niters_FO ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
 0 = iseed ! rnd seed (0=assigned automatically=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
 1 = lpp1 ! beam 1 type (0 = no PDF)
 1 = lpp2 ! beam 2 type (0 = no PDF)
 6500.0 = ebeam1 ! beam 1 energy in GeV
 6500.0 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
 nn23nlo = pdlabel ! PDF set
 244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
              ! numbers for central PDF sets are allowed. Can be a list;
              ! PDF sets beyond the first are included via reweighting.
#***********************************************************************
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#***********************************************************************
  PYTHIA8 = parton_shower
  1.0 = shower_scale_factor ! multiply default shower starting
                                  ! scale by this factor
#***********************************************************************
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses divided by two of all final state particles *
# and partons. This can be changed in SubProcesses/set_scales.f or via *
# dynamical_scale_choice option) *
#***********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.118 = muR_ref_fixed ! fixed ren reference scale
 91.118 = muF_ref_fixed ! fixed fact reference scale
 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
           ! dynamical choices. Can be a list; scale choices beyond the
           ! first are included via reweighting
 1.0 = muR_over_ref ! ratio of current muR over reference muR
 1.0 = muF_over_ref ! ratio of current muF over reference muF

#***********************************************************************
# Reweight variables for scale dependence and PDF uncertainty *
#***********************************************************************
 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
 True = reweight_scale ! Reweight to get scale variation using the
            ! rw_rscale and rw_fscale factors. Should be a list of
            ! booleans of equal length to dynamical_scale_choice to
            ! specify for which choice to include scale dependence.
 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
            ! list booleans of equal length to lhaid to specify for
            ! which PDF set to include the uncertainties.
#***********************************************************************
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#***********************************************************************
 False = store_rwgt_info ! Store info for reweighting in LHE file
#***********************************************************************
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#***********************************************************************
 3 = ickkw
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
# written in the LHE event file *
#***********************************************************************
 15.0 = bwcutoff
#***********************************************************************
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - If gamma_is_j, photons are also clustered with jets. *
# Otherwise, they will be treated as tagged particles and photon *
# isolation will be applied. Note that photons in the real emission *
# will always be clustered with QCD partons. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  -1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
  0.4 = jetradius ! The radius parameter for the jet algorithm
 30.0 = ptj ! Min jet transverse momentum
 4.5 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
 False = gamma_is_j! Wether to cluster photons as jets or not
#***********************************************************************
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  30.0 = ptl ! Min lepton transverse momentum
  2.5 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
  0.0 = drll ! Min distance between opposite sign lepton pairs
  0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
  0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
  30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#***********************************************************************
# Fermion-photon recombination parameters *
# If Rphreco=0, no recombination is performed *
#***********************************************************************
 0.1 = Rphreco ! Minimum fermion-photon distance for recombination
 -1.0 = etaphreco ! Maximum abs(pseudo-rap) for photons to be recombined (a value .lt.0 means no cut)
 False = lepphreco ! Recombine photons and leptons together
 False = quarkphreco ! Recombine photons and quarks together
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# Not applied if gamma_is_j *
# When ptgmin=0, all the other parameters are ignored *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  20.0 = ptgmin ! Min photon transverse momentum
  -1.0 = etagamma ! Max photon abs(pseudo-rap)
  0.4 = R0gamma ! Radius of isolation code
  1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
  1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
 True = isoEM ! isolate photons from EM energy (photons and leptons)
#***********************************************************************
# Cuts associated to MASSIVE particles identified by their PDG codes. *
# All cuts are applied to both particles and anti-particles, so use *
# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
# {6:100, 25:200} for multiple particles *
#***********************************************************************
  {} = pt_min_pdg ! Min pT for a massive particle
  {} = pt_max_pdg ! Max pT for a massive particle
  {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
#***********************************************************************
# Use PineAPPL to generate PDF-independent fast-interpolation grid *
# (https://zenodo.org/record/3992765#.X2EWy5MzbVo) *
#***********************************************************************
 False = pineappl ! PineAPPL switch
#***********************************************************************
]]>
</MGRunCard>
<slha>
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      6 1.730000e+02 # mt
      23 9.118800e+01 # mz
      24 8.041900e+01 # mw
      25 1.250000e+02 # mh
      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
      5 0.000000e+00 # b : 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
      15 0.000000e+00 # ta- : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.325070e+02 # aewm1
      3 1.190000e-01 # as (note that parameter not used if you use a pdf set)
###################################
## INFORMATION FOR WOLFENSTEIN
###################################
BLOCK WOLFENSTEIN #
      1 2.253000e-01 # lamws
      2 8.080000e-01 # aws
      3 1.320000e-01 # rhows
      4 3.410000e-01 # etaws
###################################
## 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
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
###################################
## INFORMATION FOR QNUMBERS 82
###################################
BLOCK QNUMBERS 82 # gh
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
</slha>
<run_settings>
order = NLO
fixed_order = ON
shower = OFF
madspin = OFF
reweight = OFF
madanalysis = OFF
runshower = False
</run_settings>
<foanalyse>
<![CDATA[
#######################################################################
#
# This file contains the settings for analyses to be linked to fixed
# order runs. Analysis files are meant to be put (or linked) inside
# <PROCDIR>/FixedOrderAnalysis/ (<PROCDIR> is the name of the exported
# process directory). See the
# <PROCDIR>/FixedOrderAnalysis/analysis_*_template.f file for details
# on how to write your own analysis.
#
#######################################################################
#
# Analysis format.
# Can either be 'topdrawer', 'root', 'HwU', 'LHE' or 'none'.
# When choosing HwU, it comes with a GnuPlot wrapper. When choosing
# topdrawer, the histogramming package 'dbook.f' is included in the
# code, while when choosing root the 'rbook_fe8.f' and 'rbook_be8.cc'
# are included. If 'none' is chosen, all the other entries below have
# to be set empty.
FO_ANALYSIS_FORMAT = HwU
#
#
# Needed extra-libraries (FastJet is already linked):
FO_EXTRALIBS =
#
# (Absolute) path to the extra libraries. Directory names should be
# separated by white spaces.
FO_EXTRAPATHS =
#
# (Absolute) path to the dirs containing header files needed by the
# libraries (e.g. C++ header files):
FO_INCLUDEPATHS =
#
# User's analysis (to be put in the <PROCDIR>/FixedOrderAnalysis/
# directory). Please use .o as extension and white spaces to separate
# files.
FO_ANALYSE = analysis_HwU_general.o
#
#
## When linking with root, the following settings are a working
## example on lxplus (CERN) as of July 2014. When using this, comment
## out the lines above and replace <PATH_TO_ROOT> with the physical
## path to root,
## e.g. /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.11/x86_64-slc6-gcc46-dbg/root/
#FO_ANALYSIS_FORMAT = root
#FO_EXTRALIBS = Core Cint Hist Matrix MathCore RIO dl Thread
#FO_EXTRAPATHS = <PATH_TO_ROOT>/lib
#FO_INCLUDEPATHS = <PATH_TO_ROOT>/include
#FO_ANALYSE = analysis_root_template.o
]]>
</foanalyse>
</header>
</LesHouchesEvents>

Here is the log file from the run :

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.0 2022-05-06 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch auto
Traceback (most recent call last):
  File "/home/jtoucheque/MG5_aMC_v3_4_0/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/jtoucheque/MG5_aMC_v3_4_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/jtoucheque/MG5_aMC_v3_4_0/madgraph/interface/amcatnlo_run_interface.py", line 1784, in do_launch
    evt_file = self.run(mode, options)
  File "/home/jtoucheque/MG5_aMC_v3_4_0/madgraph/interface/amcatnlo_run_interface.py", line 1940, in run
    jobs_to_run,jobs_to_collect=self.collect_the_results(options,req_acc,jobs_to_run, \
  File "/home/jtoucheque/MG5_aMC_v3_4_0/madgraph/interface/amcatnlo_run_interface.py", line 2259, in collect_the_results
    self.append_the_results(jobs_to_run,integration_step)
  File "/home/jtoucheque/MG5_aMC_v3_4_0/madgraph/interface/amcatnlo_run_interface.py", line 2750, in append_the_results
    raise aMCatNLOError('An error occurred during the collection of results.\n' +
madgraph.interface.amcatnlo_run_interface.aMCatNLOError: An error occurred during the collection of results.
Please check the .log files inside the directories which failed:
/home/jtoucheque/MG5_aMC_v3_4_0/SM_NLO_FIXEDORDER_CMS_jjw+/SubProcesses/P0_gu_wpd/all_G1/log.txt
.... (list of all directories associated to each process)

The log file in one of the SubProcesses directories is the following:

 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > 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: T F
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11900000000000002
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
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 = 91.188000000000002
 aEWM1 = 132.50700000000001
 aS = 0.11900000000000002
 mdl_lamWS = 0.22530000000000000
 mdl_AWS = 0.80800000000000005
 mdl_rhoWS = 0.13200000000000001
 mdl_etaWS = 0.34100000000000003
 mdl_MT = 173.00000000000000
 mdl_MZ = 91.188000000000002
 mdl_MH = 125.00000000000000
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 0.0000000000000000
 mdl_WH = 6.3823389999999999E-003
 mdl_MW = 80.418999999999997
  Internal Params
  ---------------------------------

 mdl_lhv = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_I1x33 = 0.0000000000000000
 mdl_I4x33 = 0.0000000000000000
 mdl_CKM33 = 1.0000000000000000
 mdl_conjg__CKM3x3 = 1.0000000000000000
 CMASS_mdl_MH = (125.00000004073425,-3.19116949896008093E-003)
 CMASS_mdl_MW = (80.418999999999997,0.0000000000000000)
 CMASS_mdl_MZ = (91.196168726742613,-1.2205926578948283)
 CMASS_mdl_MT = (173.00160731255821,-0.74574307143234730)
 mdl_ymt = (173.00160731255821,-0.74574307143234730)
 mdl_Ncol = (3.0000000000000000,0.0000000000000000)
 mdl_CA = (3.0000000000000000,0.0000000000000000)
 mdl_TF = (0.50000000000000000,0.0000000000000000)
 mdl_CF = (1.3333333333333333,0.0000000000000000)
 mdl_lamWS__exp__2 = (5.07600900000000008E-002,0.0000000000000000)
 mdl_CKM11 = (0.97461995499999998,0.0000000000000000)
 mdl_CKM12 = (0.22530000000000000,0.0000000000000000)
 mdl_complexi = (0.0000000000000000,1.0000000000000000)
 mdl_lamWS__exp__3 = (1.14362482770000006E-002,0.0000000000000000)
 mdl_CKM13 = (1.21974449623171223E-003,-3.15100661526525631E-003)
 mdl_CKM21 = (-0.22530000000000000,0.0000000000000000)
 mdl_CKM22 = (0.97461995499999998,0.0000000000000000)
 mdl_CKM23 = (4.10141527200000022E-002,0.0000000000000000)
 mdl_CKM31 = (8.02074411158428847E-003,-3.15100661526525631E-003)
 mdl_CKM32 = (-4.10141527200000022E-002,-0.0000000000000000)
 mdl_MZ__exp__2 = (8315.2513440000002,-222.62674795199996)
 mdl_MZ__exp__4 = (69093842.244990125,-3702394.7302364339)
 mdl_sqrt__2 = (1.4142135623730951,0.0000000000000000)
 mdl_MH__exp__2 = (15625.000000000000,-0.79779237500000000)
 mdl_CKM1x1 = (0.97461995499999998,0.0000000000000000)
 mdl_CKM1x2 = (0.22530000000000000,0.0000000000000000)
 mdl_CKM1x3 = (1.21974449623171223E-003,-3.15100661526525631E-003)
 mdl_CKM2x1 = (-0.22530000000000000,0.0000000000000000)
 mdl_CKM2x2 = (0.97461995499999998,0.0000000000000000)
 mdl_CKM2x3 = (4.10141527200000022E-002,0.0000000000000000)
 mdl_CKM3x1 = (8.02074411158428847E-003,-3.15100661526525631E-003)
 mdl_CKM3x2 = (-4.10141527200000022E-002,-0.0000000000000000)
 mdl_conjg__CKM1x3 = (1.21974449623171223E-003,3.15100661526525631E-003)
 mdl_conjg__CKM2x3 = (4.10141527200000022E-002,-0.0000000000000000)
 mdl_conjg__CKM2x1 = (-0.22530000000000000,-0.0000000000000000)
 mdl_conjg__CKM3x1 = (8.02074411158428847E-003,3.15100661526525631E-003)
 mdl_conjg__CKM2x2 = (0.97461995499999998,-0.0000000000000000)
 mdl_conjg__CKM3x2 = (-4.10141527200000022E-002,0.0000000000000000)
 mdl_Ncol__exp__2 = (9.0000000000000000,0.0000000000000000)
 mdl_MT__exp__2 = (29929.000000000004,-258.02949999999998)
 mdl_conjg__CKM11 = (0.97461995499999998,-0.0000000000000000)
 mdl_conjg__CKM12 = (0.22530000000000000,-0.0000000000000000)
 mdl_conjg__CKM13 = (1.21974449623171223E-003,3.15100661526525631E-003)
 mdl_conjg__CKM21 = (-0.22530000000000000,-0.0000000000000000)
 mdl_conjg__CKM22 = (0.97461995499999998,-0.0000000000000000)
 mdl_conjg__CKM23 = (4.10141527200000022E-002,-0.0000000000000000)
 mdl_conjg__CKM31 = (8.02074411158428847E-003,3.15100661526525631E-003)
 mdl_conjg__CKM32 = (-4.10141527200000022E-002,0.0000000000000000)
 mdl_aEW = (7.54677111397888348E-003,0.0000000000000000)
 Gf = (1.15341318059218847E-005,1.07719755703855034E-006)
 mdl_sqrt__aEW = (8.68721538467815552E-002,0.0000000000000000)
 mdl_ee = (0.30795376724436879,0.0000000000000000)
 mdl_MW__exp__2 = (6467.2155609999991,0.0000000000000000)
 mdl_sw2 = (0.22280363504870571,-2.08081141617049332E-002)
 mdl_cw = (0.88166638559643473,1.18004465757354134E-002)
 mdl_sqrt__sw2 = (0.47253402945974582,-2.20175827183230755E-002)
 mdl_sw = (0.47253402945974582,-2.20175827183230755E-002)
 mdl_g1 = (0.34922348844991208,-4.67409576430357840E-003)
 mdl_gw = (0.65029523469870720,3.03003132656498343E-002)
 mdl_v = (246.79492935034500,-11.499336413181682)
 mdl_v__exp__2 = (60775.502415098257,-5675.9558353340462)
 mdl_lam = (0.12743596415884256,1.18949573345149291E-002)
 mdl_yt = (0.98940537059448697,4.18277039442539814E-002)
 mdl_muH = (88.388347677121914,-2.25649759263027048E-003)
 mdl_AxialZUp = (-0.18447574045214427,-6.12270387629205645E-003)
 mdl_AxialZDown = (0.18447574045214427,6.12270387629205645E-003)
 mdl_VectorZUp = (7.45310271927428908E-002,1.27211881094349798E-002)
 mdl_VectorZDown = (-0.12950338382244359,-9.42194599286351811E-003)
 mdl_VectorAUp = (0.20530251149624587,0.0000000000000000)
 mdl_VectorADown = (-0.10265125574812294,-0.0000000000000000)
 mdl_VectorWmDxU = (0.22991408511437664,1.07127784911088490E-002)
 mdl_AxialWmDxU = (-0.22991408511437664,-1.07127784911088490E-002)
 mdl_VectorWpUxD = (0.22991408511437664,1.07127784911088490E-002)
 mdl_AxialWpUxD = (-0.22991408511437664,-1.07127784911088490E-002)
 mdl_I2x13 = (7.80396792833590140E-003,3.45311217803416720E-003)
 mdl_I2x23 = (-4.05796229715504886E-002,-1.71552783749657936E-003)
 mdl_I2x33 = (0.98940537059448697,4.18277039442539814E-002)
 mdl_I3x31 = (8.06756667199530313E-003,-2.78213355781023507E-003)
 mdl_I3x32 = (-4.05796229715504886E-002,-1.71552783749657936E-003)
 mdl_I3x33 = (0.98940537059448697,4.18277039442539814E-002)
 mdl_Vector_tdGp = (-7.80396792833590140E-003,-3.45311217803416720E-003)
 mdl_Vector_tsGp = (4.05796229715504886E-002,1.71552783749657936E-003)
 mdl_Vector_tbGp = (-0.98940537059448697,-4.18277039442539814E-002)
 mdl_Axial_tdGp = (-7.80396792833590140E-003,-3.45311217803416720E-003)
 mdl_Axial_tsGp = (4.05796229715504886E-002,1.71552783749657936E-003)
 mdl_Axial_tbGp = (-0.98940537059448697,-4.18277039442539814E-002)
 mdl_Vector_tdGm = (8.06756667199530313E-003,-2.78213355781023507E-003)
 mdl_Vector_tsGm = (-4.05796229715504886E-002,-1.71552783749657936E-003)
 mdl_Vector_tbGm = (0.98940537059448697,4.18277039442539814E-002)
 mdl_Axial_tdGm = (-8.06756667199530313E-003,2.78213355781023507E-003)
 mdl_Axial_tsGm = (4.05796229715504886E-002,1.71552783749657936E-003)
 mdl_Axial_tbGm = (-0.98940537059448697,-4.18277039442539814E-002)
 mdl_gw__exp__2 = (0.42196578328785017,3.94082986530602183E-002)
 mdl_cw__exp__2 = (0.77719636495129429,2.08081141617049332E-002)
 mdl_ee__exp__2 = (9.48355227599988748E-002,0.0000000000000000)
 mdl_sw__exp__2 = (0.22280363504870573,-2.08081141617049332E-002)
 mdl_yt__exp__2 = (0.97717343054396588,8.27691098441621920E-002)
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = (0.34496376621320685,0.0000000000000000)
 mdl_G__exp__4 = (2.2362154867812229,0.0000000000000000)
 mdl_G__exp__2 = (1.4953981031087418,0.0000000000000000)
 mdl_R2MixedFactor_FIN_ = (-2.52525843039140656E-002,-0.0000000000000000)
 mdl_G_UVg_1EPS_ = (-5.20834551268227613E-002,-0.0000000000000000)
 mdl_G_UVb_1EPS_ = (3.15657303798925820E-003,0.0000000000000000)
 mdl_G__exp__3 = (1.8286695702272273,0.0000000000000000)
 mdl_MU_R__exp__2 = (8315.2513440000002,0.0000000000000000)
 mdl_G_UVt_FIN_ = (-4.04285583279235909E-003,2.72133644204851571E-005)
  Couplings of loop_sm-zeromass_ckm
  ---------------------------------

         GC_4 = -0.12229E+01 0.00000E+00
         GC_5 = 0.00000E+00 0.12229E+01
         GC_6 = 0.00000E+00 0.14954E+01
       R2_3Gq = 0.77201E-02 0.00000E+00
       R2_3Gg = 0.31846E-01 0.00000E+00
       R2_GQQ = -0.00000E+00 -0.30880E-01
       R2_GGq = 0.00000E+00 0.63131E-02
       R2_GGt = -0.97739E+01 -0.11337E+04
     R2_GGg_1 = 0.00000E+00 0.28409E-01
     R2_GGg_2 = -0.00000E+00 -0.18939E-01
       R2_QQq = 0.00000E+00 0.12626E-01
  UV_3Gg_1eps = 0.63691E-01 0.00000E+00
  UV_3Gb_1eps = -0.38601E-02 -0.00000E+00
 UV_GQQg_1eps = 0.00000E+00 -0.63691E-01
 UV_GQQb_1eps = 0.00000E+00 0.38601E-02
 UVWfct_G_2_1eps -0.31566E-02 -0.00000E+00
      R2_bxcW = 0.22191E-04 -0.47625E-03
      R2_cxdW = -0.12190E-03 0.26161E-02
      R2_cxsW = 0.52732E-03 -0.11317E-01
      R2_uxbW = 0.37249E-04 -0.12459E-04
       UV_3Gt = 0.49439E-02 -0.33278E-04
      UV_GQQt = -0.33278E-04 -0.49439E-02
   UVWfct_G_2 = 0.40429E-02 -0.27213E-04
        GC_12 = -0.20882E-01 0.44816E+00
        GC_13 = -0.48272E-02 0.10360E+00
        GC_17 = -0.87875E-03 0.18859E-01
        GC_41 = -0.14751E-02 0.49336E-03

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

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23nlo
 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= 202 for QCD,QED, = 2 , 2 ,
 AMP_SPLIT: 1 correspond to S.O. 2 2
 orders_tag_plot= 204 for QCD,QED, = 4 , 2 ,
 AMP_SPLIT: 2 correspond to S.O. 4 2
 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
#--------------------------------------------------------------------------
# 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.
#--------------------------------------------------------------------------
 ------- iteration 1
 Update # PS points (even_rn): 280 --> 256
Using random seed offsets: 0 , 1 , 0
  with seed 53
 Ranmar initialization seeds 6247 9429
 initial-final FKS maps:
           0 : 10 1 2 3 4 5 6 7 8 9 10
           1 : 1 3 0 0 0 0 0 0 0 0 0
           2 : 9 1 2 4 5 6 7 8 9 10 0
 Total number of FKS directories is 10
 For the Born we use nFKSprocesses:
           1 2 3 1 1 1 1 1 1 2
tau_min 1 1 : 0.80419E+02 -- 0.11042E+03
tau_min 2 1 : 0.80419E+02 -- 0.11042E+03
tau_min 3 1 : 0.11042E+03 0.11042E+03 0.11042E+03
tau_min 4 1 : 0.80419E+02 -- 0.11042E+03
tau_min 5 1 : 0.80419E+02 -- 0.11042E+03
tau_min 6 1 : 0.80419E+02 -- 0.11042E+03
tau_min 7 1 : 0.80419E+02 -- 0.11042E+03
tau_min 8 1 : 0.80419E+02 -- 0.11042E+03
tau_min 9 1 : 0.80419E+02 -- 0.11042E+03
tau_min 10 1 : 0.80419E+02 -- 0.11042E+03
 Scale values (may change event by event):
 muR, muR_reference: 0.153497D+03 0.153497D+03 1.00
 muF1, muF1_reference: 0.153497D+03 0.153497D+03 1.00
 muF2, muF2_reference: 0.153497D+03 0.153497D+03 1.00
 QES, QES_reference: 0.750524D+02 0.750524D+02 1.00

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

 alpha_s= 0.11027800010843979
 BORN: keeping split order 1
 counterterm S.O 1 QCD
 BORN: keeping split order 1
 counterterm S.O 2 QED
 BORN: not keeping split order 1
 Charge-linked born are not used
 Color-linked born are used
 alpha_s value used for the virtuals is (for the first PS point): 0.11027800010843979
  ==========================================================================================
 { }
 { }
 { ,, }
 { `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. }
 { v3.4.0 (2022-05-06), Ref: arXiv:1103.0621v2, arXiv:1405.0301 }
 { }
 { }
  ==========================================================================================

  +----------------------------------------------------------------+
  | |
  | 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]] |
  | |
  +----------------------------------------------------------------+

 ===============================================================
 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 #

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0 0x7f84c440dd21 in ???
#1 0x7f84c440cef5 in ???
#2 0x7f84c408a08f in ???
 at /build/glibc-SzIz7B/glibc-2.31/signal/../sysdeps/unix/sysv/linux/x86_64/sigaction.c:0
#3 0x5575f6b37509 in ???
#4 0x5575f6b2069a in ???
#5 0x5575f6bd521d in ???
#6 0x5575f6b8e350 in ???
#7 0x5575f6bd66cb in ???
#8 0x5575f6bbcbf3 in ???
#9 0x5575f6bbeba6 in ???
#10 0x5575f6bd81cf in ???
#11 0x5575f6bd9580 in ???
#12 0x7f84c406b082 in __libc_start_main
 at ../csu/libc-start.c:308
#13 0x5575f6b1e5ad in ???
#14 0xffffffffffffffff in ???
Time in seconds: 1