MadGraph event generation produces errors while running on the cluster, but not while running on my PC

Asked by Thomas Hugle on 2019-06-13

Hi,

I'm trying to simulate the process g g > l+ l- xd xd~ [QCD] in a two Higgs doublet + pseudoscalar model, where xd is dark matter. Doing so on my PC works fine (using MG 2.6.4), but when trying to simulate the same process in the same way on our cluster (using MG 2.6.5) the event generation produces errors that look like the following:
...
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
INFO: Using LHAPDF v6.1.6 interface for PDFs
initMadLoop -r -f
Initializing MadLoop loop-induced matrix elements (this can take some time)...
MadLoop initialization finished.
compile Source Directory
Using random number seed offset = 11
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: P0_gg_llxdxdx
INFO: Idle: 0, Running: 14, Completed: 0 [ current time: 18h13 ]
STOP 1
ERROR DETECTED
STOP 1
ERROR DETECTED
STOP 1
ERROR DETECTED
STOP 1
ERROR DETECTED
INFO: Idle: 0, Running: 13, Completed: 1 [ 13m 4s ]
INFO: Idle: 0, Running: 12, Completed: 2 [ 15m 25s ]
INFO: Idle: 0, Running: 11, Completed: 3 [ 16m 16s ]
...

Each "STOP 1; ERROR DETECTED" pair corresponds to an error "rm: cannot remove ‘results.dat’: No such file or directory".

Interestingly the errors do not appear when forcing an s-channel propagator via "generate g g > h2 > l+ l- xd xd~ [QCD]" or some other simpler processes that I ran for testing purposes.
My suspicion is that it could be related to the fact that on the cluster we have gcc version 4.8.5, whereas on my PC there is gcc version 5.4.0. For Python the cluster is running version 2.7.5 and my PC uses 2.7.12.

Thank you in advance for any insights on how to solve the problem.
For completeness I will post underneath the full debug log for one of the runs.

Best,
Thomas

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 2469, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 6869, in new_fct
    original_fct(obj, *args, **opts)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 2511, in run_generate_events
    self.exec_cmd('refine %s' % nb_event, postcmd=False)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 3421, in do_refine
    x_improve.launch() # create the ajob for the refinment.
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/madevent/gen_ximprove.py", line 860, in launch
    main_dir=pjoin(self.cmd.me_dir,'SubProcesses')) #main_dir is for gridpack readonly mode
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/madevent/sum_html.py", line 747, in collect_result
    P_comb.add_results(os.path.basename(G), path, mfactors[G])
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/madevent/sum_html.py", line 425, in add_results
    oneresult.read_results(filepath)
  File "/lfs/l8/lin/thugle/MG5_aMC_v2_6_5/madgraph/madevent/sum_html.py", line 315, in read_results
    % (error_code, open(log).read())
Exception: Reported error: End code 1.0
 Full associated log:
==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.1.6 loading /lfs/l8/lin/thugle/MG5_aMC_v2_6_5/HEPTools/lhapdf6/share/LHAPDF/NNPDF30_lo_as_0130/NNPDF30_lo_as_0130_0000.dat
NNPDF30_lo_as_0130 PDF set, member #0, version 2; LHAPDF ID = 263000
 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 lhapdf : 0.13000256957493794
 using LHAPDF
 Define smin to 400.00000000000000
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 MU_R = 91.188000000000002
 mdl_gPXd = 1.0000000000000000
 mdl_lam3 = 0.0000000000000000
 mdl_laP1 = 0.0000000000000000
 mdl_laP2 = 0.0000000000000000
 mdl_sinp = 0.29999999999999999
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 172.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_tanbeta = 1.0000000000000000
 mdl_sinbma = 1.0000000000000000
 mdl_MZ = 91.187600000000003
 mdl_MTA = 1.7769999999999999
 mdl_MT = 172.00000000000000
 mdl_MB = 4.7000000000000002
 mdl_mhc = 100.00000000000000
 mdl_mh1 = 125.00000000000000
 mdl_mh2 = 100.00000000000000
 mdl_mh3 = 100.00000000000000
 mdl_mh4 = 100.00000000000000
 mdl_MXd = 10.000000000000000
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WT = 1.5083359999999999
 mdl_whc = 2.0852590000000000E-004
 mdl_Wh1 = 6.5336320000000002E-003
 mdl_Wh2 = 4.4988019999999997E-003
 mdl_Wh3 = 0.35499249999999999
 mdl_Wh4 = 3.5480290000000001
  Internal Params
  ---------------------------------

 mdl_mh1__exp__2 = 15625.000000000000
 mdl_mh2__exp__2 = 10000.000000000000
 mdl_dm2 = -5625.0000000000000
 mdl_mh3__exp__2 = 10000.000000000000
 mdl_mh4__exp__2 = 10000.000000000000
 mdl_asin__sinp = 0.30469265401539752
 mdl_sqrt__2 = 1.4142135623730951
 mdl_mA = 100.00000000000000
 mdl_mP = 100.00000000000000
 mdl_asin__sinbma = 1.5707963267948966
 mdl_mu1 = -7812.5000000000000
 mdl_mu3 = 0.0000000000000000
 mdl_TH1x1 = 1.0000000000000000
 mdl_sinbma__exp__2 = 1.0000000000000000
 mdl_TH1x2 = 0.0000000000000000
 mdl_TH2x1 = -0.0000000000000000
 mdl_TH2x2 = 1.0000000000000000
 mdl_sinp__exp__2 = 8.9999999999999997E-002
 mdl_TH3x3 = 0.95393920141694566
 mdl_TH3x4 = 0.29999999999999999
 mdl_TH4x3 = -0.29999999999999999
 mdl_TH4x4 = 0.95393920141694566
 mdl_tmc = 0.0000000000000000
 mdl_tanbeta__exp__2 = 1.0000000000000000
 mdl_lP = 0.0000000000000000
 mdl_lP1 = 0.0000000000000000
 mdl_lP2 = 0.0000000000000000
 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_mhc__exp__2 = 10000.000000000000
 mdl_tanbeta__exp__3 = 1.0000000000000000
 mdl_mA__exp__2 = 10000.000000000000
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_TH1x1__exp__3 = 1.0000000000000000
 mdl_TH1x1__exp__2 = 1.0000000000000000
 mdl_TH2x1__exp__2 = 0.0000000000000000
 mdl_TH2x1__exp__3 = -0.0000000000000000
 mdl_TH1x2__exp__2 = 0.0000000000000000
 mdl_TH2x2__exp__2 = 1.0000000000000000
 mdl_TH1x2__exp__3 = 0.0000000000000000
 mdl_TH2x2__exp__3 = 1.0000000000000000
 mdl_TH1x2__exp__4 = 0.0000000000000000
 mdl_TH2x2__exp__4 = 1.0000000000000000
 mdl_TH3x3__exp__2 = 0.91000000000000003
 mdl_TH3x4__exp__2 = 8.9999999999999997E-002
 mdl_TH4x3__exp__2 = 8.9999999999999997E-002
 mdl_TH3x3__exp__3 = 0.86808467328942063
 mdl_TH3x3__exp__4 = 0.82810000000000006
 mdl_TH4x4__exp__2 = 0.91000000000000003
 mdl_TH3x4__exp__3 = 2.7000000000000000E-002
 mdl_TH3x4__exp__4 = 8.0999999999999996E-003
 mdl_TH1x1__exp__4 = 1.0000000000000000
 mdl_TH2x1__exp__4 = 0.0000000000000000
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29584.000000000000
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824660036055974
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6371.9763498719121
 mdl_sw2 = 0.23369336794341478
 mdl_cw = 0.87538941737753784
 mdl_sqrt__sw2 = 0.48341841911889827
 mdl_sw = 0.48341841911889827
 mdl_g1 = 0.35807035569216145
 mdl_gw = 0.64840516548963911
 mdl_vev = 246.21845810181625
 mdl_vev__exp__2 = 60623.529110035844
 mdl_mu2 = 17812.500000000000
 mdl_l1 = 0.12886910601690277
 mdl_l2 = 0.12886910601690277
 mdl_l3 = -0.25773821203380554
 mdl_l4 = 0.0000000000000000
 mdl_l5 = 0.0000000000000000
 mdl_l6 = 0.0000000000000000
 mdl_l7 = 0.0000000000000000
 mdl_muP = 0.0000000000000000
 mdl_yb = 2.6995554250465501E-002
 mdl_yt = 0.98792241086809907
 mdl_ytau = 1.0206617000654720E-002
 mdl_I1a33 = ( 2.6995554250465501E-002, 0.0000000000000000 )
 mdl_I2a33 = (-0.98792241086809907 , 0.0000000000000000 )
 mdl_I3a33 = (-0.98792241086809907 , 0.0000000000000000 )
 mdl_I4a33 = ( 2.6995554250465501E-002, 0.0000000000000000 )
 mdl_I5a33 = ( 2.6995554250465501E-002, 0.0000000000000000 )
 mdl_I6a33 = ( 0.98792241086809907 , 0.0000000000000000 )
 mdl_I7a33 = ( 0.98792241086809907 , 0.0000000000000000 )
 mdl_I8a33 = ( 2.6995554250465501E-002, 0.0000000000000000 )
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_sw__exp__2 = 0.23369336794341478
 mdl_cw__exp__2 = 0.76630663205658511
 mdl_yb__exp__2 = 7.2875994928982594E-004
 mdl_yt__exp__2 = 0.97599068989543714
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_G__exp__3 = 1.8056676068262196
 mdl_G__exp__4 = 2.1987899468922913
  Couplings of Pseudoscalar_2HDM
  ---------------------------------

  R2GC_319_21 = -0.00000E+00 -0.16849E-02
  R2GC_319_25 = 0.00000E+00 -0.22565E+01
  R2GC_320_27 = -0.00000E+00 -0.16849E-02
  R2GC_320_31 = 0.00000E+00 0.22565E+01
        GC_10 = 0.00000E+00 0.12177E+01
       GC_156 = 0.30000E+00 0.00000E+00
       GC_172 = -0.95394E+00 0.00000E+00
         GC_1 = -0.00000E+00 -0.10448E+00
       GC_128 = -0.00000E+00 0.00000E+00
       GC_130 = 0.35329E+00 0.00000E+00
       GC_149 = -0.00000E+00 0.00000E+00
       GC_151 = 0.11111E+00 0.00000E+00
         GC_2 = 0.00000E+00 0.20897E+00
       GC_256 = -0.18209E-01 0.00000E+00
       GC_257 = -0.57266E-02 0.00000E+00
       GC_258 = -0.00000E+00 -0.19089E-01
       GC_259 = -0.00000E+00 -0.19089E-01
       GC_283 = 0.66639E+00 0.00000E+00
       GC_285 = 0.20957E+00 0.00000E+00
       GC_299 = 0.00000E+00 -0.69857E+00
         GC_3 = -0.00000E+00 -0.31345E+00
       GC_300 = 0.00000E+00 0.69857E+00
        GC_53 = 0.00000E+00 0.28380E+00
        GC_60 = -0.00000E+00 -0.28850E-01
        GC_62 = 0.00000E+00 0.17310E+00
        GC_64 = 0.00000E+00 -0.19725E+00

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

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = lhapdf
 alpha_s(Mz)= 0.1300 running at 1 loops.
 alpha_s(Mz)= 0.1300 running at 1 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 4000 5 3
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.0000000000000000E-002
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Monte-Carlo over helicities

 ********************************************
 * You are using the DiscreteSampler module *
 * part of the MG5_aMC framework *
 * Author: Valentin Hirschi *
 ********************************************

Enter Configuration Number:
Running Configuration Number: 1
 Using dconfig= 14
 BW Setting 2 1 1 0
 Attempting mappinvarients 1 6
 Completed mapping 6
 about to integrate 10 4000 5 3 10 1
 Using non-zero grid deformation.
 10 dimensions 4000 events 10 invarients 5 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 24
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 13000.00
 Mapping Graph 1 to config 1
Setting grid 1 0.23669E-05 1
 Setting BW -2 2 91.187600000000003
 Setting PDF BW 9 3 125.00000000000000
 Transforming s_hat BW 9.2455621301775159E-005 4.8325680473372789E-009
   1 1 2 3 4 5 6 7 8 9 10
 Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.100E+02 0.100E+02
Using random seed offsets 1 : 1
  with seed 14
 Ranmar initialization seeds 18869 9388
  Particle 3 4 5 6
      Et > 10.0 10.0 0.0 0.0
       E > 0.0 0.0 0.0 0.0
     Eta < 2.5 2.5 -1.0 -1.0
   xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.4 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 0.0 0.0
s min # 4> 0.0 0.0 0.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
 Excluding BW -1 1
 Requiring BW -2 2
 Requiring BW -3 3
 alpha_s for scale 125.00588138446317 is 0.12380727518388480
  ==========================================================================================
 { }
 {   }
 {  ,,  }
 { `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.6.5 (2018-02-03), 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-004
  > NRotations_DP = 1
  > NRotations_QP = 0
  > CTStabThres = 1.0000000000000000E-002
  > CTLoopLibrary = 2
  > CTModeInit = 1
  > CheckCycle = 4
  > MaxAttempts = 10
  > UseLoopFilter = F
  > HelicityFilterLevel = 1
  > ImprovePSPoint = 2
  > DoubleCheckHelicityFilter = F
  > LoopInitStartOver = F
  > HelInitStartOver = F
  > ZeroThres = 1.0000000000000001E-009
  > OSThres = 1.0000000000000000E-008
  > WriteOutFilters = F
  > 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
 ===============================================================
 ##INFO: For loop-induced processes it is preferable to always set the parameter LoopInitStartOver to True, so it is hard-set here to True.
 ##INFO: For loop-induced processes it is preferable to always set the parameter HelInitStartOver to True, so it is hard-set here to True.
 ##INFO: Due to the dynamic setting of the reference scale for contributions comparisons, it is preferable to set the parameter CheckCycle to a value larger than 4, so it is hard-set here to 5.

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

  +----------------------------------------------------------------+
  | |
  | Ninja - version 1.2.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]] |
  | |
  +----------------------------------------------------------------+

STOP 1
Thanks for using LHAPDF 6.1.6. Please make sure to cite the paper:
  Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
# #
# 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. #
# #
########################################################################
 RESET CUMULATIVE VARIABLE
 RESET CUMULATIVE VARIABLE

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

 Iteration 1 Mean: 0.3922E-16 Abs mean: 0.3922E-16 Fluctuation: 0.283E-17 0.398E-14 2.7%
  1 0.3922E-16 0.3922E-16 +- 0.2834E-17 4.57
 Writing out events 7.2648745727604389E-020 4.5698840605824955
 Relative summed weights:
  0.1000E+01 0.0000E+00
 Relative number of events:
  0.1000E+01 0.0000E+00
 Events:
        4000 0
 DiscreteSampler:: Error, no point could be picked with random variable 0.7042701840400696 using upper bound found of 0.0000000000000000.

ls status:
input_app.txt
run1_app.log

                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
          cluster_temp_path : None
             cluster_memory : None (user set)
               cluster_size : 100
              cluster_queue : None (user set)
                    nb_core : 1 (user set)
               cluster_time : 1 (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : /lfs/l8/lin/thugle/MG5_aMC_v2_6_5/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /lfs/l8/lin/thugle/MG5_aMC_v2_6_5/Delphes (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : /lfs/l8/lin/thugle/MG5_aMC_v2_6_5/HEPTools/madanalysis5/madanalysis5 (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : /lfs/l8/lin/thugle/MG5_aMC_v2_6_5/HEPTools/lhapdf6/bin/lhapdf-config (user set)
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : None (user set)
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* 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 default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set low_mem_multicore_nlo_generation False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model Pseudoscalar_2HDM
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~
generate g g > l+ l- xd xd~ [QCD]
output cluster_runs/pseudoscalar/job_100_100
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR DMINPUTS
###################################
BLOCK DMINPUTS #
      1 1.000000e+00 # gpxd
###################################
## INFORMATION FOR FRBLOCK
###################################
BLOCK FRBLOCK #
      2 1.000000e+00 # tanbeta
      3 1.000000e+00 # sinbma
###################################
## INFORMATION FOR HIGGS
###################################
BLOCK HIGGS #
      1 0.000000e+00 # lam3
      2 0.000000e+00 # lap1
      3 0.000000e+00 # lap2
      5 3.000000e-01 # sinp
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      5 4.700000e+00 # mb
      6 1.720000e+02 # mt
      15 1.777000e+00 # mta
      23 9.118760e+01 # mz
      25 1.250000e+02 # mh1
      35 1.000000e+02 # mh2
      36 1.000000e+02 # mh3
      37 1.000000e+02 # mhc
      55 1.000000e+02 # mh4
      42000 1.000000e+01 # mxd
      1 0.000000e+00 # d : 0.0
      2 0.000000e+00 # u : 0.0
      3 0.000000e+00 # s : 0.0
      4 0.000000e+00 # c : 0.0
      11 0.000000e+00 # e- : 0.0
      12 0.000000e+00 # ve : 0.0
      13 0.000000e+00 # mu- : 0.0
      14 0.000000e+00 # vm : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 7.982466e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
      9000002 9.118760e+01 # ghz : mz
      9000003 7.982466e+01 # ghwp : mw
      9000004 7.982466e+01 # ghwm : mw
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.279000e+02 # aewm1
      2 1.166390e-05 # gf
      3 1.180000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      5 4.700000e+00 # ymb
      6 1.720000e+02 # ymt
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
BLOCK QNUMBERS 9000001 # gha
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000002
###################################
BLOCK QNUMBERS 9000002 # ghz
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000003
###################################
BLOCK QNUMBERS 9000003 # ghwp
      1 3 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000004
###################################
BLOCK QNUMBERS 9000004 # ghwm
      1 -3 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 82
###################################
BLOCK QNUMBERS 82 # ghg
      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)
###################################
## INFORMATION FOR QNUMBERS 37
###################################
BLOCK QNUMBERS 37 # h+
      1 3 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 35
###################################
BLOCK QNUMBERS 35 # h2
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 36
###################################
BLOCK QNUMBERS 36 # h3
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 55
###################################
BLOCK QNUMBERS 55 # h4
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 42000
###################################
BLOCK QNUMBERS 42000 # xd
      1 0 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.508336e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.495200e+00
#
# PDG Width
DECAY 24 2.085000e+00
#
# PDG Width
DECAY 25 6.533632e-03
# BR NDA ID1 ID2 ...
   8.250815e-01 2 5 -5 # 0.00539077895804
   3.960247e-02 2 15 -15 # 0.000258747981073
   2.126229e-02 3 -2 1 24 # 0.000138920006097
   2.126229e-02 3 -4 3 24 # 0.000138920006097
   2.111690e-02 3 -24 -1 2 # 0.00013797003549
   2.111690e-02 3 -24 -3 4 # 0.00013797003549
   7.083273e-03 3 -12 11 24 # 4.62795001586e-05
   7.083273e-03 3 -14 13 24 # 4.62795001586e-05
   7.018148e-03 3 -24 -15 16 # 4.58539982005e-05
   7.000012e-03 3 -24 -11 12 # 4.57355016341e-05
   7.000012e-03 3 -24 -13 14 # 4.57355016341e-05
   6.997640e-03 3 -16 15 24 # 4.5720003277e-05
   1.921749e-03 3 -1 1 23 # 1.25560033233e-05
   1.921749e-03 3 -3 3 23 # 1.25560033233e-05
   1.539878e-03 3 -5 5 23 # 1.00609976904e-05
   1.497514e-03 3 -2 2 23 # 9.78420215959e-06
   1.494361e-03 3 -4 4 23 # 9.76360251775e-06
#
# PDG Width
DECAY 35 4.498802e-03
# BR NDA ID1 ID2 ...
   9.540195e-01 2 5 -5 # 0.00429194527529
   4.598048e-02 2 15 -15 # 0.000206857099433
#
# PDG Width
DECAY 36 3.549925e-01
# BR NDA ID1 ID2 ...
   9.883689e-01 2 42000 -42000 # 0.350863560555
   1.110020e-02 2 5 -5 # 0.00394048835563
   5.309352e-04 2 15 -15 # 0.000188478025495
#
# PDG Width
DECAY 37 2.085259e-04
# BR NDA ID1 ID2 ...
   9.932522e-01 2 16 -15 # 0.000207118771357
   6.747844e-03 3 -5 5 24 # 1.4071e-06
#
# PDG Width
DECAY 55 3.548029e+00
# BR NDA ID1 ID2 ...
   9.998849e-01 2 42000 -42000 # 3.54762044561
   1.098409e-04 2 5 -5 # 0.000389718628578
   5.253814e-06 2 15 -15 # 1.86406838401e-05
#
# PDG Width
DECAY 42000 0.000000e+00
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#*********************************************************************
#
#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
  10000 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     6500.0 = ebeam1 ! beam 1 total energy in GeV
     6500.0 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     lhapdf = pdlabel ! PDF set
     263000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.188 = scale ! fixed ren scale
 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  3.0 = lhe_version ! Change the way clustering information pass to shower.
  True = clusinfo ! include clustering tag in output
  average = event_norm ! average/sum. Normalization of the weight in the LHEF

#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
 0 = ickkw ! 0 no matching, 1 MLM
 1.0 = alpsfact ! scale factor for QCD emission vx
 False = chcluster ! cluster only according to channel diag
 4 = asrwgtflavor ! highest quark flavor for a_s reweight
 False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
                                   ! (turn off for VBF and single top processes)
 0.0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
   1 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
   False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20.0 = ptj ! minimum pt for the jets
 0.0 = ptb ! minimum pt for the b
 10.0 = pta ! minimum pt for the photons
 10.0 = ptl ! minimum pt for the charged leptons
 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
 -1.0 = ptjmax ! maximum pt for the jets
 -1.0 = ptbmax ! maximum pt for the b
 -1.0 = ptamax ! maximum pt for the photons
 -1.0 = ptlmax ! maximum pt for the charged leptons
 -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
 {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0.0 = ej ! minimum E for the jets
  0.0 = eb ! minimum E for the b
  0.0 = ea ! minimum E for the photons
  0.0 = el ! minimum E for the charged leptons
  -1.0 = ejmax ! maximum E for the jets
 -1.0 = ebmax ! maximum E for the b
 -1.0 = eamax ! maximum E for the photons
 -1.0 = elmax ! maximum E for the charged leptons
 {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
  5.0 = etaj ! max rap for the jets
  -1.0 = etab ! max rap for the b
 2.5 = etaa ! max rap for the photons
 2.5 = etal ! max rap for the charged leptons
 0.0 = etajmin ! min rap for the jets
 0.0 = etabmin ! min rap for the b
 0.0 = etaamin ! min rap for the photons
 0.0 = etalmin ! main rap for the charged leptons
 {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
 0.4 = drjj ! min distance between jets
 0.0 = drbb ! min distance between b's
 0.4 = drll ! min distance between leptons
 0.4 = draa ! min distance between gammas
 0.0 = drbj ! min distance between b and jet
 0.4 = draj ! min distance between gamma and jet
 0.4 = drjl ! min distance between jet and lepton
 0.0 = drab ! min distance between gamma and b
 0.0 = drbl ! min distance between b and lepton
 0.4 = dral ! min distance between gamma and lepton
 -1.0 = drjjmax ! max distance between jets
 -1.0 = drbbmax ! max distance between b's
 -1.0 = drllmax ! max distance between leptons
 -1.0 = draamax ! max distance between gammas
 -1.0 = drbjmax ! max distance between b and jet
 -1.0 = drajmax ! max distance between gamma and jet
 -1.0 = drjlmax ! max distance between jet and lepton
 -1.0 = drabmax ! max distance between gamma and b
 -1.0 = drblmax ! max distance between b and lepton
 -1.0 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#*********************************************************************
 0.0 = mmjj ! min invariant mass of a jet pair
 0.0 = mmbb ! min invariant mass of a b pair
 0.0 = mmaa ! min invariant mass of gamma gamma pair
 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1.0 = mmjjmax ! max invariant mass of a jet pair
 -1.0 = mmbbmax ! max invariant mass of a b pair
 -1.0 = mmaamax ! max invariant mass of gamma gamma pair
 -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
 {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
 {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = ptheavy ! minimum pt for at least one heavy final state
 0.0 = xptj ! minimum pt for at least one jet
 0.0 = xptb ! minimum pt for at least one b
 0.0 = xpta ! minimum pt for at least one photon
 0.0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0.0 = ptj1min ! minimum pt for the leading jet in pt
 0.0 = ptj2min ! minimum pt for the second jet in pt
 0.0 = ptj3min ! minimum pt for the third jet in pt
 0.0 = ptj4min ! minimum pt for the fourth jet in pt
 -1.0 = ptj1max ! maximum pt for the leading jet in pt
 -1.0 = ptj2max ! maximum pt for the second jet in pt
 -1.0 = ptj3max ! maximum pt for the third jet in pt
 -1.0 = ptj4max ! maximum pt for the fourth jet in pt
 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
 0.0 = ptl1min ! minimum pt for the leading lepton in pt
 0.0 = ptl2min ! minimum pt for the second lepton in pt
 0.0 = ptl3min ! minimum pt for the third lepton in pt
 0.0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1.0 = ptl1max ! maximum pt for the leading lepton in pt
 -1.0 = ptl2max ! maximum pt for the second lepton in pt
 -1.0 = ptl3max ! maximum pt for the third lepton in pt
 -1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
 0.0 = ihtmin !inclusive Ht for all partons (including b)
 -1.0 = ihtmax !inclusive Ht for all partons (including b)
 0.0 = ht2min ! minimum Ht for the two leading jets
 0.0 = ht3min ! minimum Ht for the three leading jets
 0.0 = ht4min ! minimum Ht for the four leading jets
 -1.0 = ht2max ! maximum Ht for the two leading jets
 -1.0 = ht3max ! maximum Ht for the three leading jets
 -1.0 = ht4max ! maximum Ht for the four leading jets
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#***********************************************************************
 0.0 = ptgmin ! Min photon transverse momentum
 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)
#*********************************************************************
# WBF cuts *
#*********************************************************************
 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#***********************************************************************
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#***********************************************************************
 -1.0 = ktdurham
 0.4 = dparameter
 -1.0 = ptlund
 1, 2, 3, 4, 5, 6, 21, 82 = pdgs_for_merging_cut ! PDGs for two cuts above
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Valentin Hirschi Edit question
Solved by:
Thomas Hugle
Solved:
2019-07-05
Last query:
2019-07-05
Last reply:
2019-07-05

Hi,

Not sure what I would be able to do for you here.
Would you be able to put a link to your model, such that I can try locally to see what the issue can be?
I'm more worried about some IO issue in this case actually. But it would be nice to run this process with full compiler security flag to be sure that noting special occurs.

Thanks,

Olivier

Thomas Hugle (thugle3) said : #2

Hi,

thanks for your answer!

The model can be found under
https://github.com/LHC-DMWG/model-repository/tree/master/models/Pseudoscalar_2HDM

Is there an easy way to run the process with full compiler security flag for me?

Best,
Thomas

Thanks I will take a look.

> Is there an easy way to run the process with full compiler security flag
> for me?

Sure you can go to the Source directory.
You have a file make_opts
and inside a line
#FFLAGS+= -g -fbounds-check -ffpe-trap=invalid,zero,overflow,underflow,denormal -Wall

That you can uncomment or replace it by
FFLAGS+= -g -fcheck=all -ffpe-trap=invalid,zero,overflow,underflow,denormal -Wall
.
then in that directory you can do a "make clean" to be sure that everyting will be recompile.

Cheers,

Olivier

> On 19 Jun 2019, at 10:12, Thomas Hugle <email address hidden> wrote:
>
> Question #681396 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/681396
>
> Status: Answered => Open
>
> Thomas Hugle is still having a problem:
> Hi,
>
> thanks for your answer!
>
> The model can be found under
> https://github.com/LHC-DMWG/model-repository/tree/master/models/Pseudoscalar_2HDM
>
> Is there an easy way to run the process with full compiler security flag
> for me?
>
> Best,
> Thomas
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thomas Hugle (thugle3) said : #4

As an update:

I ran the process with full compiler security flag and the errors appear to all be of the same type and happen at the same place in the code.
The error is a "Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation" and the full backtrace is:
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
#0 0x2AE04EF56697
#1 0x2AE04EF56CDE
#2 0x2AE04F9E927F
#3 0xF75FAC in __reductiond_MOD_calcdpv1
#4 0xF7F971 in __reductiond_MOD_calcdred
#5 0xF90C63 in __reductiond_MOD_calcd
#6 0xE4C8BB in __reductiontn_MOD_calctn
#7 0xEA6452 in __collier_coefs_MOD_tn_main_checked_cll
#8 0xEB17F0 in __collier_coefs_MOD_tn_main_cll
#9 0xD6F6AC in __collier_tensors_MOD_tnten_main_checked_cll
#10 0xD72AA1 in __collier_tensors_MOD_tnten_main_cll
#11 0x467064 in ml5_0_0_1_collierloop_ at COLLIER_interface.f:362
#12 0x463F33 in ml5_0_0_1_tirloop_ at TIR_interface.f:204 (discriminator 1)
#13 0x485C67 in ml5_0_0_1_loop_4_ at CT_interface.f:1013 (discriminator 8)
#14 0x45B919 in ml5_0_0_1_loop_ct_calls_1_ at loop_CT_calls_1.f:114
#15 0x449385 in ml5_0_0_1_sloopmatrix_ at loop_matrix.f:1132
#16 0x44F9F4 in ml5_0_0_1_sloopmatrixhel_ at loop_matrix.f:2939 (discriminator 2)
#17 0x45025B in ml5_0_0_1_sloopmatrixhel_thres_ at loop_matrix.f:2981 (discriminator 2)
#18 0x44090D in smatrix1_ at matrix1.f:309 (discriminator 2)
#19 0x43F3FB in dsig1_ at auto_dsig1.f:125
#20 0x43C7C8 in dsigproc_ at auto_dsig.f:474 (discriminator 1)
#21 0x43E38E in dsig_ at auto_dsig.f:266
#22 0x10EBB8F in sample_full_ at dsample.f:137
#23 0x40BA0A in MAIN__ at driver.f:177
ERROR DETECTED

Maybe this helps in pinpointing the source of the error.

Best,
Thomas

Thanks,

However this points to an issue in collier where the only proposal that I can give is to unlink collier to see if this improves the situation. (I will put Valentin in cc if he has something smarter to tell)

One other idea is to reduce the flag to
 -g -fcheck=all -Wall

Since not allowing floating point exception might be too trash for external code like Collier.

Cheers,

Olivier

Thomas Hugle (thugle3) said : #6

Thanks for the ideas on how to address the problem.
Unluckily neither of them freed me of the errors.

Reducing the flag to
-g -fcheck=all -Wall
I get the (assumingly) same error as in the beginning "STOP 1; ERROR DETECTED".

When disabling collier and running with full compiler security flag I also still get errors, but the backtrace changes a bit:
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
#0 0x2AE86402A697
#1 0x2AE86402ACDE
#2 0x2AE864ABD27F
#3 0xD97885 in __discretesampler_MOD_ds_get_point_with_binid at DiscreteSampler.f:2049
#4 0xD9A05D in __discretesampler_MOD_ds_get_point_with_integerbinid at DiscreteSampler.f:1877
#5 0x43E2F9 in dsig_ at auto_dsig.f:330
#6 0xD8A34F in sample_full_ at dsample.f:137
#7 0x40B10A in MAIN__ at driver.f:177
ERROR DETECTED

Best,
Thomas

Thomas Hugle (thugle3) said : #7

Here a bit more information on things that I tried:
Changing the compiler from gcc-4.8.5 to gcc-5.4.0 or gcc-7.2.0 did not help for the two Higgs doublet + pseudoscalar model, however it seems to have helped with another model (two Higgs doublet + scalar).

Digging a bit deeper, the error that triggers the breakdown of the simulation seems to be happening in the DiscreteSampler module (in principle this is already apparent from the full debug log in my very first post, but I did not fully realize it back then):
DiscreteSampler:: Error, no point could be picked with random variable 0.8378990292549133 using upper bound found of 0.0000000000000000.
(An error like this appears in the "run1_app.log" files of the SubProcesses that fail. Interestingly, not all SubProcesses fail.)

Also, before the error there are floating-point exceptions (cf. "run1_app.log" files of the SubProcesses):
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP 1

For completeness and since it might help understanding the problem, I will add the log files of a SubProcess that fails on the cluster but works on my PC.

Best,
Thomas

I leave out the model parameters in both files (they are identical to the full debug log posted above)

CLUSTER (giving an error):
=====================

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

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = lhapdf
 alpha_s(Mz)= 0.1300 running at 1 loops.
 alpha_s(Mz)= 0.1300 running at 1 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 4000 5 3
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.0000000000000000E-002
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Monte-Carlo over helicities

 ********************************************
 * You are using the DiscreteSampler module *
 * part of the MG5_aMC framework *
 * Author: Valentin Hirschi *
 ********************************************

Enter Configuration Number:
Running Configuration Number: 1
 Using dconfig= 14
 BW Setting 2 1 1 0
 Attempting mappinvarients 1 6
 Completed mapping 6
 about to integrate 10 4000 5 3 10 1
 Using non-zero grid deformation.
 10 dimensions 4000 events 10 invarients 5 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 24
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 13000.00
 Mapping Graph 1 to config 1
Setting grid 1 0.23669E-05 1
 Setting BW -2 2 91.187600000000003
 Setting PDF BW 9 3 125.00000000000000
 Transforming s_hat BW 9.2455621301775159E-005 4.8325680473372789E-009
   1 1 2 3 4 5 6 7 8 9 10
 Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.100E+02 0.100E+02
Using random seed offsets 1 : 1
  with seed 14
 Ranmar initialization seeds 18869 9388
  Particle 3 4 5 6
      Et > 10.0 10.0 0.0 0.0
       E > 0.0 0.0 0.0 0.0
     Eta < 2.5 2.5 -1.0 -1.0
   xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.4 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 100.0 0.0
s min # 4> 0.0 0.0 100.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
 Excluding BW -1 1
 Requiring BW -2 2
 Requiring BW -3 3
 alpha_s for scale 125.00588138446317 is 0.12380727518388480
  ==========================================================================================
 { }
 { ^[[32m ^[[0m }
 { ^[[32m ,, ^[[0m }
 { ^[[32m`7MMM. ,MMF' `7MM `7MMF' ^[[0m }
 { ^[[32m MMMb dPMM MM MM ^[[0m }
 { ^[[32m M YM ,M MM ,6"Yb. ,M""bMM MM ,pW"Wq. ,pW"Wq.`7MMpdMAo. ^[[0m }
 { ^[[32m M Mb M' MM 8) MM ,AP MM MM 6W' `Wb 6W' `Wb MM `Wb ^[[0m }
 { ^[[32m M YM.P' MM ,pm9MM 8MI MM MM , 8M M8 8M M8 MM M8 ^[[0m }
 { ^[[32m M `YM' MM 8M MM `Mb MM MM ,M YA. ,A9 YA. ,A9 MM ,AP ^[[0m }
 { ^[[32m.JML. `' .JMML.`Moo9^Yo.`Wbmd"MML..JMMmmmmMMM `Ybmd9' `Ybmd9' MMbmmd' ^[[0m }
 { ^[[32m MM ^[[0m }
 { ^[[32m .JMML. ^[[0m }
 { ^[[32m^[[0mv2.6.6 (2018-06-28), Ref: arXiv:1103.0621v2, arXiv:1405.0301^[[32m ^[[0m }
 { ^[[32m ^[[0m }
 { }
  ==========================================================================================
 ===============================================================
 INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
 ===============================================================
  > MLReductionLib = 6|7|1
  > CTModeRun = -1
  > MLStabThres = 1.0000000000000000E-004
  > NRotations_DP = 1
  > NRotations_QP = 0
  > CTStabThres = 1.0000000000000000E-002
  > CTLoopLibrary = 2
  > CTModeInit = 1
  > CheckCycle = 4
  > MaxAttempts = 10
  > UseLoopFilter = F
  > HelicityFilterLevel = 1
  > ImprovePSPoint = 2
  > DoubleCheckHelicityFilter = F
  > LoopInitStartOver = F
  > HelInitStartOver = F
  > ZeroThres = 1.0000000000000001E-009
  > OSThres = 1.0000000000000000E-008
  > WriteOutFilters = F
  > 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
 ===============================================================
 ##INFO: For loop-induced processes it is preferable to always set the parameter LoopInitStartOver to True, so it is hard-set here to True.
 ##INFO: For loop-induced processes it is preferable to always set the parameter HelInitStartOver to True, so it is hard-set here to True.
 ##INFO: Due to the dynamic setting of the reference scale for contributions comparisons, it is preferable to set the parameter CheckCycle to a value larger than 4, so it is hard-set here to 5.
------------------------------------------------------------------------
| 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, #

  +----------------------------------------------------------------+
  | |
  | Ninja - version 1.2.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]] |
  | |
  +----------------------------------------------------------------+

Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP 1
Thanks for using LHAPDF 6.1.6. Please make sure to cite the paper:
  Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
# 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. #
# #
########################################################################
 RESET CUMULATIVE VARIABLE
 RESET CUMULATIVE VARIABLE

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

 Iteration 1 Mean: 0.3922E-16 Abs mean: 0.3922E-16 Fluctuation: 0.283E-17 0.398E-14 2.7%
  1 0.3922E-16 0.3922E-16 +- 0.2834E-17 4.57
 Writing out events 7.2648745727604389E-020 4.5698840605824955
 Relative summed weights:
  0.1000E+01 0.0000E+00
 Relative number of events:
  0.1000E+01 0.0000E+00
 Events:
        4000 0
 DiscreteSampler:: Error, no point could be picked with random variable 0.7042701840400696 using upper bound found of 0.0000000000000000.

ls status:
input_app.txt
run1_app.log

PC (works):
=========
 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = lhapdf
 alpha_s(Mz)= 0.1300 running at 1 loops.
 alpha_s(Mz)= 0.1300 running at 1 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 2000 1 1
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.0000000000000000E-002
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Monte-Carlo over helicities

 ********************************************
 * You are using the DiscreteSampler module *
 * part of the MG5_aMC framework *
 * Author: Valentin Hirschi *
 ********************************************

Enter Configuration Number:
Running Configuration Number: 1
 Using dconfig= 14
 BW Setting 2 1 1 0
 Attempting mappinvarients 1 6
 Completed mapping 6
 about to integrate 10 2000 1 1 10 1
 Using unknown grid deformation: -2
 10 dimensions 2000 events 10 invarients 1 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 24
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 13000.00
 Mapping Graph 1 to config 1
Setting grid 1 0.23669E-05 1
 Setting BW -2 2 91.187600000000003
 Setting PDF BW 9 3 125.00000000000000
 Transforming s_hat BW 9.2455621301775159E-005 4.8325680473372789E-009
   1 1 2 3 4 5 6 7 8 9 10
 Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.100E+02 0.100E+02
 Got moffset 1
Using random seed offsets 1 : 1
  with seed 14
 Ranmar initialization seeds 18869 12545
  Particle 3 4 5 6
      Et > 10.0 10.0 0.0 0.0
       E > 0.0 0.0 0.0 0.0
     Eta < 2.5 2.5 -1.0 -1.0
   xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.4 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 0.0 0.0
s min # 4> 0.0 0.0 0.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
 Excluding BW -1 1
 Requiring BW -2 2
 Requiring BW -3 3
 alpha_s for scale 124.99908171426037 is 0.12380829279978629
  ==========================================================================================
 { }
 {   }
 {  ,,  }
 { `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.6.5 (2018-02-03), 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-004
  > NRotations_DP = 1
  > NRotations_QP = 0
  > CTStabThres = 1.0000000000000000E-002
  > CTLoopLibrary = 2
  > CTModeInit = 1
  > CheckCycle = 4
  > MaxAttempts = 10
  > UseLoopFilter = F
  > HelicityFilterLevel = 1
  > ImprovePSPoint = 2
  > DoubleCheckHelicityFilter = F
  > LoopInitStartOver = F
  > HelInitStartOver = F
  > ZeroThres = 1.0000000000000001E-009
  > OSThres = 1.0000000000000000E-008
  > WriteOutFilters = F
  > 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
 ===============================================================
 ##INFO: For loop-induced processes it is preferable to always set the parameter LoopInitStartOver to True, so it is hard-set here to True.
 ##INFO: For loop-induced processes it is preferable to always set the parameter HelInitStartOver to True, so it is hard-set here to True.
 ##INFO: Due to the dynamic setting of the reference scale for contributions comparisons, it is preferable to set the parameter CheckCycle to a value larger than 4, so it is hard-set here to 5.

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

  +----------------------------------------------------------------+
  | |
  | Ninja - version 1.2.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]] |
  | |
  +----------------------------------------------------------------+

Thanks for using LHAPDF 6.1.6. Please make sure to cite the paper:
  Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
# 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. #
# #
########################################################################
 RESET CUMULATIVE VARIABLE
 RESET CUMULATIVE VARIABLE
 Iteration 1 Mean: 0.3474E-16 Abs mean: 0.3474E-16 Fluctuation: 0.301E-17 0.266E-14 2.2%
  1 0.3474E-16 0.3474E-16 +- 0.3009E-17 3.87
 Writing out events 7.9538527730590813E-019 3.8733910142143664
 Relative summed weights:
  0.1000E+01 0.0000E+00
 Relative number of events:
  0.1000E+01 0.0000E+00
 Events:
        2000 0
 Accuracy: 0.000 0.020 0.087 0.000
 Finished due to accuracy 0.0000000000000000 2.0000000000000000E-002

 -------------------------------------------------------------------------------
 Accumulated results: Integral = 0.3474E-16
                        Std dev = 0.3009E-17
                       Cross sec = 0.3474E-16
             Chi**2 per DoF. = 0.0000
 -------------------------------------------------------------------------------
-------------------------------------------------
---------------------------
 Results Last 1 iters: Integral = 0.3474E-16
                     Abs integral = 0.3474E-16
                          Std dev = 0.3009E-17
                  Chi**2 per DoF. = 0.0000
-------------------------------------------------
---------------------------
 Status 2.0000000000000000E-002 2 1

ls status:
ftn26
grid_information
input_app.txt
moffset.dat
results.dat
run1_app.log

Hi,

DiscreteSampler:: Error, no point could be picked with random variable 0.8378990292549133 using upper bound found of 0.0000000000000000.

This message is due to a consistency check of the code.
The code realises that something is wrong and stops.

One thing that you need to check is if the IEEE_DIVIDE_BY_ZERO is coming from your model.
One method is to put your param_card within the UFO model directory and renamed it restrict_xxx.dat
then do import model 2HDM_pseudoscalar-xxx
If you have a division by zero, it should crash and point you to the issue in your code.

Doing that will also remove from the model any coupling which are zero (or less than 1e-10).
Which can also help fo the other flag (or hurt you if you need one of such suppressed coupling)

Cheers,

Olivier

On 5 Jul 2019, at 11:47, Thomas Hugle <<email address hidden><mailto:<email address hidden>>> wrote:

Question #681396 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/681396

Thomas Hugle gave more information on the question:
Here a bit more information on things that I tried:
Changing the compiler from gcc-4.8.5 to gcc-5.4.0 or gcc-7.2.0 did not help for the two Higgs doublet + pseudoscalar model, however it seems to have helped with another model (two Higgs doublet + scalar).

Digging a bit deeper, the error that triggers the breakdown of the simulation seems to be happening in the DiscreteSampler module (in principle this is already apparent from the full debug log in my very first post, but I did not fully realize it back then):
DiscreteSampler:: Error, no point could be picked with random variable 0.8378990292549133 using upper bound found of 0.0000000000000000.
(An error like this appears in the "run1_app.log" files of the SubProcesses that fail. Interestingly, not all SubProcesses fail.)

Also, before the error there are floating-point exceptions (cf. "run1_app.log" files of the SubProcesses):
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP 1

For completeness and since it might help understanding the problem, I
will add the log files of a SubProcess that fails on the cluster but
works on my PC.

Best,
Thomas

I leave out the model parameters in both files (they are identical to
the full debug log posted above)

CLUSTER (giving an error):
=====================

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

Running at P P machine @ 13000.000000000000 GeV
PDF set = lhapdf
alpha_s(Mz)= 0.1300 running at 1 loops.
alpha_s(Mz)= 0.1300 running at 1 loops.
Renormalization scale set on event-by-event basis
Factorization scale set on event-by-event basis

getting user params
Enter number of events and max and min iterations:
Number of events and iterations 4000 5 3
Enter desired fractional accuracy:
Desired fractional accuracy: 2.0000000000000000E-002
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Monte-Carlo over helicities

********************************************
* You are using the DiscreteSampler module *
* part of the MG5_aMC framework *
* Author: Valentin Hirschi *
********************************************

Enter Configuration Number:
Running Configuration Number: 1
Using dconfig= 14
BW Setting 2 1 1 0
Attempting mappinvarients 1 6
Completed mapping 6
about to integrate 10 4000 5 3 10 1
Using non-zero grid deformation.
10 dimensions 4000 events 10 invarients 5 iterations 1 config(s), (0.99)
Using h-tuple random number sequence.
Error opening grid
Using Uniform Grid! 24
Using uniform alpha 1.0000000000000000
Grid defined OK
Set CM energy to 13000.00
Mapping Graph 1 to config 1
Setting grid 1 0.23669E-05 1
Setting BW -2 2 91.187600000000003
Setting PDF BW 9 3 125.00000000000000
Transforming s_hat BW 9.2455621301775159E-005 4.8325680473372789E-009
  1 1 2 3 4 5 6 7 8 9 10
Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.100E+02 0.100E+02
Using random seed offsets 1 : 1
 with seed 14
Ranmar initialization seeds 18869 9388
 Particle 3 4 5 6
     Et > 10.0 10.0 0.0 0.0
      E > 0.0 0.0 0.0 0.0
    Eta < 2.5 2.5 -1.0 -1.0
  xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.4 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 100.0 0.0
s min # 4> 0.0 0.0 100.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
Excluding BW -1 1
Requiring BW -2 2
Requiring BW -3 3
alpha_s for scale 125.00588138446317 is 0.12380727518388480
 ==========================================================================================
{ }
{ ^[[32m ^[[0m }
{ ^[[32m ,, ^[[0m }
{ ^[[32m`7MMM. ,MMF' `7MM `7MMF' ^[[0m }
{ ^[[32m MMMb dPMM MM MM ^[[0m }
{ ^[[32m M YM ,M MM ,6"Yb. ,M""bMM MM ,pW"Wq. ,pW"Wq.`7MMpdMAo. ^[[0m }
{ ^[[32m M Mb M' MM 8) MM ,AP MM MM 6W' `Wb 6W' `Wb MM `Wb ^[[0m }
{ ^[[32m M YM.P' MM ,pm9MM 8MI MM MM , 8M M8 8M M8 MM M8 ^[[0m }
{ ^[[32m M `YM' MM 8M MM `Mb MM MM ,M YA. ,A9 YA. ,A9 MM ,AP ^[[0m }
{ ^[[32m.JML. `' .JMML.`Moo9^Yo.`Wbmd"MML..JMMmmmmMMM `Ybmd9' `Ybmd9' MMbmmd' ^[[0m }
{ ^[[32m MM ^[[0m }
{ ^[[32m .JMML. ^[[0m }
{ ^[[32m^[[0mv2.6.6 (2018-06-28), Ref: arXiv:1103.0621v2, arXiv:1405.0301^[[32m ^[[0m }
{ ^[[32m ^[[0m }
{ }
 ==========================================================================================
===============================================================
INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
===============================================================
MLReductionLib = 6|7|1
CTModeRun = -1
MLStabThres = 1.0000000000000000E-004
NRotations_DP = 1
NRotations_QP = 0
CTStabThres = 1.0000000000000000E-002
CTLoopLibrary = 2
CTModeInit = 1
CheckCycle = 4
MaxAttempts = 10
UseLoopFilter = F
HelicityFilterLevel = 1
ImprovePSPoint = 2
DoubleCheckHelicityFilter = F
LoopInitStartOver = F
HelInitStartOver = F
ZeroThres = 1.0000000000000001E-009
OSThres = 1.0000000000000000E-008
WriteOutFilters = F
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
===============================================================
##INFO: For loop-induced processes it is preferable to always set the parameter LoopInitStartOver to True, so it is hard-set here to True.
##INFO: For loop-induced processes it is preferable to always set the parameter HelInitStartOver to True, so it is hard-set here to True.
##INFO: Due to the dynamic setting of the reference scale for contributions comparisons, it is preferable to set the parameter CheckCycle to a value larger than 4, so it is hard-set here to 5.
------------------------------------------------------------------------
| 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, #

 +----------------------------------------------------------------+
 | |
 | Ninja - version 1.2.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]] |
 | |
 +----------------------------------------------------------------+

Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP 1
Thanks for using LHAPDF 6.1.6. Please make sure to cite the paper:
 Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
# 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. #
# #
########################################################################
RESET CUMULATIVE VARIABLE
RESET CUMULATIVE VARIABLE

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

Iteration 1 Mean: 0.3922E-16 Abs mean: 0.3922E-16 Fluctuation: 0.283E-17 0.398E-14 2.7%
 1 0.3922E-16 0.3922E-16 +- 0.2834E-17 4.57
Writing out events 7.2648745727604389E-020 4.5698840605824955
Relative summed weights:
 0.1000E+01 0.0000E+00
Relative number of events:
 0.1000E+01 0.0000E+00
Events:
       4000 0
DiscreteSampler:: Error, no point could be picked with random variable 0.7042701840400696 using upper bound found of 0.0000000000000000.

ls status:
input_app.txt
run1_app.log

PC (works):
=========
Collider parameters:
--------------------

Running at P P machine @ 13000.000000000000 GeV
PDF set = lhapdf
alpha_s(Mz)= 0.1300 running at 1 loops.
alpha_s(Mz)= 0.1300 running at 1 loops.
Renormalization scale set on event-by-event basis
Factorization scale set on event-by-event basis

getting user params
Enter number of events and max and min iterations:
Number of events and iterations 2000 1 1
Enter desired fractional accuracy:
Desired fractional accuracy: 2.0000000000000000E-002
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Monte-Carlo over helicities

********************************************
* You are using the DiscreteSampler module *
* part of the MG5_aMC framework *
* Author: Valentin Hirschi *
********************************************

Enter Configuration Number:
Running Configuration Number: 1
Using dconfig= 14
BW Setting 2 1 1 0
Attempting mappinvarients 1 6
Completed mapping 6
about to integrate 10 2000 1 1 10 1
Using unknown grid deformation: -2
10 dimensions 2000 events 10 invarients 1 iterations 1 config(s), (0.99)
Using h-tuple random number sequence.
Error opening grid
Using Uniform Grid! 24
Using uniform alpha 1.0000000000000000
Grid defined OK
Set CM energy to 13000.00
Mapping Graph 1 to config 1
Setting grid 1 0.23669E-05 1
Setting BW -2 2 91.187600000000003
Setting PDF BW 9 3 125.00000000000000
Transforming s_hat BW 9.2455621301775159E-005 4.8325680473372789E-009
  1 1 2 3 4 5 6 7 8 9 10
Masses: 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.100E+02 0.100E+02
Got moffset 1
Using random seed offsets 1 : 1
 with seed 14
Ranmar initialization seeds 18869 12545
 Particle 3 4 5 6
     Et > 10.0 10.0 0.0 0.0
      E > 0.0 0.0 0.0 0.0
    Eta < 2.5 2.5 -1.0 -1.0
  xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.4 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 0.0 0.0
s min # 4> 0.0 0.0 0.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
Excluding BW -1 1
Requiring BW -2 2
Requiring BW -3 3
alpha_s for scale 124.99908171426037 is 0.12380829279978629
 ==========================================================================================
{ }
{   }
{  ,,  }
{ `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.6.5 (2018-02-03), 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-004
NRotations_DP = 1
NRotations_QP = 0
CTStabThres = 1.0000000000000000E-002
CTLoopLibrary = 2
CTModeInit = 1
CheckCycle = 4
MaxAttempts = 10
UseLoopFilter = F
HelicityFilterLevel = 1
ImprovePSPoint = 2
DoubleCheckHelicityFilter = F
LoopInitStartOver = F
HelInitStartOver = F
ZeroThres = 1.0000000000000001E-009
OSThres = 1.0000000000000000E-008
WriteOutFilters = F
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
===============================================================
##INFO: For loop-induced processes it is preferable to always set the parameter LoopInitStartOver to True, so it is hard-set here to True.
##INFO: For loop-induced processes it is preferable to always set the parameter HelInitStartOver to True, so it is hard-set here to True.
##INFO: Due to the dynamic setting of the reference scale for contributions comparisons, it is preferable to set the parameter CheckCycle to a value larger than 4, so it is hard-set here to 5.

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

 +----------------------------------------------------------------+
 | |
 | Ninja - version 1.2.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]] |
 | |
 +----------------------------------------------------------------+

Thanks for using LHAPDF 6.1.6. Please make sure to cite the paper:
 Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
# 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. #
# #
########################################################################
RESET CUMULATIVE VARIABLE
RESET CUMULATIVE VARIABLE
Iteration 1 Mean: 0.3474E-16 Abs mean: 0.3474E-16 Fluctuation: 0.301E-17 0.266E-14 2.2%
 1 0.3474E-16 0.3474E-16 +- 0.3009E-17 3.87
Writing out events 7.9538527730590813E-019 3.8733910142143664
Relative summed weights:
 0.1000E+01 0.0000E+00
Relative number of events:
 0.1000E+01 0.0000E+00
Events:
       2000 0
Accuracy: 0.000 0.020 0.087 0.000
Finished due to accuracy 0.0000000000000000 2.0000000000000000E-002

-------------------------------------------------------------------------------
Accumulated results: Integral = 0.3474E-16
                       Std dev = 0.3009E-17
                      Cross sec = 0.3474E-16
            Chi**2 per DoF. = 0.0000
-------------------------------------------------------------------------------
-------------------------------------------------
---------------------------
Results Last 1 iters: Integral = 0.3474E-16
                    Abs integral = 0.3474E-16
                         Std dev = 0.3009E-17
                 Chi**2 per DoF. = 0.0000
-------------------------------------------------
---------------------------
Status 2.0000000000000000E-002 2 1

ls status:
ftn26
grid_information
input_app.txt
moffset.dat
results.dat
run1_app.log

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

Thomas Hugle (thugle3) said : #9

Hey,

the problem is solved :)

The issue was that I had set (under input/mg5_configuration.txt) run_mode = 0 (single machine), which is automatically switched to run_mode=2 (multicore) with one core for loop induced processes as the one I am dealing with. However, this seems to have broken the discrete sampling routine.

Directly working with run_mode=2 and nb_core=4 got rid of the errors and the simulation appears to work as intended.

Best,
Thomas

Thomas Hugle (thugle3) said : #10

Thanks, for your support!

You are not using run_mode=1 on the cluster?
Do I understand correctly that you want to run on 4 core of a single machine of the cluster (via manual submision of the job?)

On your PC, what was the option used? run_mode=0 or run_mode=2?

Cheers,

Olivier

Thomas Hugle (thugle3) said : #12

Yes, I am not using run_mode=1 on the cluster. Since I am running a parameter scan, it seemed simpler to me to manually submit my jobs because didn't manage to figure out how to change my run parameters when using the cluster mode. Is there a simple way to do so?
My impression was that the cluster is well suited to calculate many events for one parameter set but not that well suited to calculate less events for many different parameter points. Of course this impression might be wrong.

Yes, I assume the cluster gives me 4 cores of a single node (via manual submission).

On my PC I did not set the run_mode option explicitly, but I assume it defaults to run_mode=2 since my computer used all of its cores.

Best,
Thomas

Thanks for the info, I will check on my laptop if I can reproduce issue with run_mode=0 then.

> Since I am running a parameter scan, it seemed simpler to me to manually submit my jobs because didn't manage to figure out how to change my run parameters when using the cluster mode. Is there a simple way to do so?

in the param_card you can set any value to
scan:[val1, val2, val3]

You have more advance syntax explained here:
https://answers.launchpad.net/mg5amcnlo/+faq/2735

> My impression was that the cluster is well suited to calculate many events for one parameter set but not that well suited to calculate less events for many different parameter points. Of course this impression might be wrong.

Not completely wrong, since you can not run within the same directory two benchmark simultaneously.
The various benchmark are done sequentially. This being said the computation of loop-induced process are very slow and therefore can easily take advantages of the cluster mode.

Now, they are nothing wrong to run like you did. I was just surprised.

Cheers,

Olivier

> On 5 Jul 2019, at 15:13, Thomas Hugle <email address hidden> wrote:
>
> Question #681396 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/681396
>
> Thomas Hugle posted a new comment:
> Yes, I am not using run_mode=1 on the cluster. Since I am running a parameter scan, it seemed simpler to me to manually submit my jobs because didn't manage to figure out how to change my run parameters when using the cluster mode. Is there a simple way to do so?
> My impression was that the cluster is well suited to calculate many events for one parameter set but not that well suited to calculate less events for many different parameter points. Of course this impression might be wrong.
>
> Yes, I assume the cluster gives me 4 cores of a single node (via manual
> submission).
>
> On my PC I did not set the run_mode option explicitly, but I assume it
> defaults to run_mode=2 since my computer used all of its cores.
>
> Best,
> Thomas
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thomas Hugle (thugle3) said : #14

Thanks for the additional information and the link.

Best,
Thomas