Loop gamma-ray process launch fails in v3.2

Could you add the content of the log file that should be present inside the G1_2 directory?
If you do not have any, can you check the log file of the G1/ directory?
and the status in other G1_X directory.

The likely issue is that you have some numerical issue on a particular phase-space point that make one loop routine to crash.
If we do not spot something with those other files (likely) I will need to check how to force to keep the log file for all instance such that we will be able to see the log of G1_2.

Cheers,

Olivier

Hi Olivier,

I replied via email with the requested log on Friday. No rush --- just let me know if I sent it successfully, since I am not experienced with Launchpad.

The actual reason I am writing before you reply is that I have done some further testing. Specifically, I find that MG5_aMC itself cannot compute the loop process that MadDM should be calculating for indirect_spectral_features, namely chi chi > a a [noborn = QED]. I am using the same top-philic scalar model with the same benchmark as in 2107.04598. Compilation fails with the following error message:

Failed the initialization of loop-induced matrix element 'PV0_0_1_chichi_aa' (trying with a maximum of 18 PS points).

So your hypothesis seems correct... but then I don't know how the example in 2107.04598 was performed by the authors.

Hi Taylor,

no the email was likely rejected due to the attachment. Since you create a question and not a bug request, you are not allowed to have attachment (this is annoying but that was the decision of ubuntu not ours). So maybe you can try via email without attachment or directly via the web interface.

> Failed the initialization of loop-induced matrix element 'PV0_0_1_chichi_aa' (trying with a maximum of 18 PS points).

This error is actually a bit random, so typically retry multiple times and it will go through. Maybe I should find the way to increase the maximum here.

Cheers,

Olivier

Hi Olivier,

Ah, thanks for letting me know about attachments.

I will forget about MadGraph for a bit and return to the MadDM problem. Instead of an attachment, I have pasted the log you originally requested in the G1_2 directory for the chichi_aa subprocess.

Obviously it is a bit long, apologies.

----------------------------------------

Process in group number 1995
 No PDF is used, alpha_s(MZ) from param_card is used
 Value of alpha_s from param_card: 0.11839999999999999
 The default order of alpha_s running is fixed to 2
 Define smin to 0.0000000000000000
 Warning: pt or E min of a gamma should in general be >0
 Warning: pt or E min of a gamma should in general be >0
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 mdl_lambdachi = 0.84999999999999998
 mdl_sSIe = 1.0000000000000000
 mdl_fSIe = 10000000000.000000
 mdl_vSIe = 1.0000000000000000
 mdl_fSDe = 10000000000.000000
 mdl_vSDe = 1.0000000000000000
 MU_R = 91.188000000000002
 mdl_MD = 5.0400000000000002E-003
 mdl_MU = 2.5500000000000002E-003
 mdl_MS = 0.10100000000000001
 mdl_MC = 1.2700000000000000
 mdl_MB = 4.7000000000000002
 mdl_MT = 172.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_MZ = 91.187600000000003
 mdl_MH = 125.00000000000000
 mdl_Mst = 622.00000000000000
 mdl_Mchi = 450.00000000000000
 mdl_sdmm = 70.000000000000000
 mdl_MFsdm = 10000000000.000000
 mdl_vdmm = 70.000000000000000
 mdl_MFvdm = 10000000000.000000
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663700000000000E-005
 aS = 0.11839999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 172.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_WT = 1.5083359999999999
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WH = 4.0699999999999998E-003
 mdl_Wsdm = 0.0000000000000000
 mdl_WFsdm = 0.0000000000000000
 mdl_Wvdm = 0.0000000000000000
 mdl_WFvdm = 0.0000000000000000
  Internal Params
  ---------------------------------

 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_complexi = (0.0000000000000000,1.0000000000000000)
 mdl_lambdachi__exp__2 = 0.72249999999999992
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29584.000000000000
 mdl_MTA__exp__2 = 3.1577289999999998
 mdl_Mst__exp__2 = 386884.00000000000
 mdl_Mchi__exp__2 = 202500.00000000000
 mdl_lambdachi__exp__4 = 0.52200624999999989
 mdl_MB__exp__3 = 103.82300000000002
 mdl_MT__exp__3 = 5088448.0000000000
 mdl_MTA__exp__3 = 5.6112844329999989
 mdl_sqrt__6 = 2.4494897427831779
 mdl_Wfdm = 0.0000000000000000
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824359746197842
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6371.9284088904105
 mdl_sw2 = 0.23369913342182447
 mdl_cw = 0.87538612427783857
 mdl_sqrt__sw2 = 0.48342438232036300
 mdl_sw = 0.48342438232036300
 mdl_g1 = 0.35807170271074895
 mdl_gw = 0.64839716719502682
 mdl_vev = 246.22056907348590
 mdl_vev__exp__2 = 60624.568634871241
 mdl_lam = 0.12886689630821144
 mdl_yb = 2.6995322804122722E-002
 mdl_yt = 0.98791394091683138
 mdl_ytau = 1.0206529494239589E-002
 mdl_muH = 88.388347648318444
 mdl_I1a33 = (2.69953228041227219E-002,0.0000000000000000)
 mdl_I2a33 = (0.98791394091683138,0.0000000000000000)
 mdl_I3a33 = (0.98791394091683138,0.0000000000000000)
 mdl_I4a33 = (2.69953228041227219E-002,0.0000000000000000)
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_sw__exp__2 = 0.23369913342182450
 mdl_cw__exp__2 = 0.76630086657817542
 mdl_yb__exp__2 = 7.2874745329878823E-004
 mdl_yt__exp__2 = 0.97597395465782455
 mdl_ee__exp__3 = 3.0797040156871746E-002
 mdl_ee__exp__4 = 9.6533630357369515E-003
 mdl_sw__exp__3 = 0.11297585922324961
 mdl_yb__exp__3 = 1.9672772744483135E-005
 mdl_yt__exp__3 = 0.96417827577819637
 mdl_ytau__exp__3 = (1.06324729062986880E-006,0.0000000000000000)
 mdl_cw__exp__3 = 0.67080914562461802
 mdl_cw__exp__4 = 0.58721701811846261
 mdl_sw__exp__4 = 5.4615284962111724E-002
 mdl_ytau__exp__2 = (1.04173244316782646E-004,0.0000000000000000)
 mdl_yb__exp__4 = 5.3107285068946954E-007
 mdl_yt__exp__4 = 0.95252516017043332
 mdl_ytau__exp__4 = (1.08520648314840882E-008,0.0000000000000000)
 mdl_lam__exp__2 = 1.6606676964111319E-002
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34409301068170506
 mdl_G__exp__2 = (1.4878582807401259,0.0000000000000000)
 mdl_G__exp__3 = (1.8148567439626970,0.0000000000000000)
 mdl_G__exp__4 = (2.2137222635669636,0.0000000000000000)
  Couplings of topphilic_NLO_EW_CM_UFO__REAL
  ---------------------------------

         GC_2 = -0.00000E+00 -0.20897E+00
         GC_3 = 0.00000E+00 0.20897E+00
        GC_45 = -0.00000E+00 -0.11540E+00
        GC_46 = 0.00000E+00 0.11540E+00
        GC_48 = -0.00000E+00 -0.48230E-01
        GC_53 = 0.00000E+00 0.25495E+00
         GC_6 = 0.00000E+00 0.87335E-01
        GC_81 = -0.00000E+00 -0.69856E+00
        GC_27 = 0.00000E+00 0.85000E+00

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

 Running at e e machine @ 900.00018983999996 GeV
 PDF set = none
 alpha_s(Mz)= 0.1184 running at 2 loops.
 alpha_s(Mz)= 0.1184 running at 2 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 1000 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
 Not subdividing B.W.
 Attempting mappinvarients 1 4
 Determine nb_t
 T-channel found: 0
 Completed mapping 4
 about to integrate 2 1000 1 1 2 1
 Using unknown grid deformation: -2
  2 dimensions 1000 events 2 invarients 1 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 16
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 900.00
 Mapping Graph 1 to config 1
 Determine nb_t
 T-channel found: 0
Setting grid 1 0.12346E-04 1
   1 1 2
 Masses: 0.450E+03 0.450E+03 0.000E+00 0.000E+00
 Got moffset 2
Using random seed offsets 1 : 1
  with seed 24
 Ranmar initialization seeds 31059 15712
  Particle 3 4
      Et > 0.0 0.0
       E > 0.0 0.0
     Eta < -1.0 -1.0
   xqcut: 0.0 0.0
d R # 3 > -0.0 0.0
s min # 3> 0.0 0.0
xqcutij # 3> 0.0 0.0
 alpha_s for scale 900.00018983999996 is 8.8077955047643222E-002
  ==========================================================================================
 { }
 {   }
 {  ,,  }
 { `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.9.4 (2021-05-30), Ref: arXiv:1103.0621v2, arXiv:1405.0301  }
 {   }
 { }
  ==========================================================================================
 ===============================================================
 INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
 ===============================================================
  > MLReductionLib = 7|6|1|3
  > CTModeRun = -1
  > MLStabThres = 1.0000000000000000E-002
  > 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
 ===============================================================
 ##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.
 INFO:: Attempting to rescue the precision improvement with an alternative method.

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

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

 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Further warnings from the improve_ps routine will now be supressed.
 RESET CUMULATIVE VARIABLE
 RESET CUMULATIVE VARIABLE
 Iteration 1 Mean: 0.1513E-01 Abs mean: 0.1513E-01 Fluctuation: 0.101E-02 0.164E+00 100.0%
  1 0.1513E-01 0.1513E-01 +- 0.1005E-02 2.10
 Writing out events 1.5126275114749196E-005 2.1020811795506842
 Relative summed weights:
  0.1000E+01 0.0000E+00
 Relative number of events:
  0.1000E+01 0.0000E+00
 Events:
        1019 0
 Accuracy: 0.000 0.020 0.066 0.000
 Finished due to accuracy 0.0000000000000000 2.0000000000000000E-002

 -------------------------------------------------------------------------------
 Accumulated results: Integral = 0.1513E-01
                        Std dev = 0.1005E-02
                       Cross sec = 0.1513E-01
             Chi**2 per DoF. = 0.0000
 -------------------------------------------------------------------------------
-------------------------------------------------
---------------------------
 Results Last 1 iters: Integral = 0.1513E-01
                     Abs integral = 0.1513E-01
                          Std dev = 0.1005E-02
                  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
 Process in group number 1995
 No PDF is used, alpha_s(MZ) from param_card is used
 Value of alpha_s from param_card: 0.11839999999999999
 The default order of alpha_s running is fixed to 2
 Define smin to 0.0000000000000000
 Warning: pt or E min of a gamma should in general be >0
 Warning: pt or E min of a gamma should in general be >0
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 mdl_lambdachi = 0.84999999999999998
 mdl_sSIe = 1.0000000000000000
 mdl_fSIe = 10000000000.000000
 mdl_vSIe = 1.0000000000000000
 mdl_fSDe = 10000000000.000000
 mdl_vSDe = 1.0000000000000000
 MU_R = 91.188000000000002
 mdl_MD = 5.0400000000000002E-003
 mdl_MU = 2.5500000000000002E-003
 mdl_MS = 0.10100000000000001
 mdl_MC = 1.2700000000000000
 mdl_MB = 4.7000000000000002
 mdl_MT = 172.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_MZ = 91.187600000000003
 mdl_MH = 125.00000000000000
 mdl_Mst = 622.00000000000000
 mdl_Mchi = 450.00000000000000
 mdl_sdmm = 70.000000000000000
 mdl_MFsdm = 10000000000.000000
 mdl_vdmm = 70.000000000000000
 mdl_MFvdm = 10000000000.000000
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663700000000000E-005
 aS = 0.11839999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 172.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_WT = 1.5083359999999999
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WH = 4.0699999999999998E-003
 mdl_Wsdm = 0.0000000000000000
 mdl_WFsdm = 0.0000000000000000
 mdl_Wvdm = 0.0000000000000000
 mdl_WFvdm = 0.0000000000000000
  Internal Params
  ---------------------------------

 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_complexi = (0.0000000000000000,1.0000000000000000)
 mdl_lambdachi__exp__2 = 0.72249999999999992
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29584.000000000000
 mdl_MTA__exp__2 = 3.1577289999999998
 mdl_Mst__exp__2 = 386884.00000000000
 mdl_Mchi__exp__2 = 202500.00000000000
 mdl_lambdachi__exp__4 = 0.52200624999999989
 mdl_MB__exp__3 = 103.82300000000002
 mdl_MT__exp__3 = 5088448.0000000000
 mdl_MTA__exp__3 = 5.6112844329999989
 mdl_sqrt__6 = 2.4494897427831779
 mdl_Wfdm = 0.0000000000000000
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824359746197842
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6371.9284088904105
 mdl_sw2 = 0.23369913342182447
 mdl_cw = 0.87538612427783857
 mdl_sqrt__sw2 = 0.48342438232036300
 mdl_sw = 0.48342438232036300
 mdl_g1 = 0.35807170271074895
 mdl_gw = 0.64839716719502682
 mdl_vev = 246.22056907348590
 mdl_vev__exp__2 = 60624.568634871241
 mdl_lam = 0.12886689630821144
 mdl_yb = 2.6995322804122722E-002
 mdl_yt = 0.98791394091683138
 mdl_ytau = 1.0206529494239589E-002
 mdl_muH = 88.388347648318444
 mdl_I1a33 = (2.69953228041227219E-002,0.0000000000000000)
 mdl_I2a33 = (0.98791394091683138,0.0000000000000000)
 mdl_I3a33 = (0.98791394091683138,0.0000000000000000)
 mdl_I4a33 = (2.69953228041227219E-002,0.0000000000000000)
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_sw__exp__2 = 0.23369913342182450
 mdl_cw__exp__2 = 0.76630086657817542
 mdl_yb__exp__2 = 7.2874745329878823E-004
 mdl_yt__exp__2 = 0.97597395465782455
 mdl_ee__exp__3 = 3.0797040156871746E-002
 mdl_ee__exp__4 = 9.6533630357369515E-003
 mdl_sw__exp__3 = 0.11297585922324961
 mdl_yb__exp__3 = 1.9672772744483135E-005
 mdl_yt__exp__3 = 0.96417827577819637
 mdl_ytau__exp__3 = (1.06324729062986880E-006,0.0000000000000000)
 mdl_cw__exp__3 = 0.67080914562461802
 mdl_cw__exp__4 = 0.58721701811846261
 mdl_sw__exp__4 = 5.4615284962111724E-002
 mdl_ytau__exp__2 = (1.04173244316782646E-004,0.0000000000000000)
 mdl_yb__exp__4 = 5.3107285068946954E-007
 mdl_yt__exp__4 = 0.95252516017043332
 mdl_ytau__exp__4 = (1.08520648314840882E-008,0.0000000000000000)
 mdl_lam__exp__2 = 1.6606676964111319E-002
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34409301068170506
 mdl_G__exp__2 = (1.4878582807401259,0.0000000000000000)
 mdl_G__exp__3 = (1.8148567439626970,0.0000000000000000)
 mdl_G__exp__4 = (2.2137222635669636,0.0000000000000000)
  Couplings of topphilic_NLO_EW_CM_UFO__REAL
  ---------------------------------

         GC_2 = -0.00000E+00 -0.20897E+00
         GC_3 = 0.00000E+00 0.20897E+00
        GC_45 = -0.00000E+00 -0.11540E+00
        GC_46 = 0.00000E+00 0.11540E+00
        GC_48 = -0.00000E+00 -0.48230E-01
        GC_53 = 0.00000E+00 0.25495E+00
         GC_6 = 0.00000E+00 0.87335E-01
        GC_81 = -0.00000E+00 -0.69856E+00
        GC_27 = 0.00000E+00 0.85000E+00

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

 Running at e e machine @ 900.00018983999996 GeV
 PDF set = none
 alpha_s(Mz)= 0.1184 running at 2 loops.
 alpha_s(Mz)= 0.1184 running at 2 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
 Not subdividing B.W.
 Attempting mappinvarients 1 4
 Determine nb_t
 T-channel found: 0
 Completed mapping 4
 about to integrate 2 2000 1 1 2 1
 Using unknown grid deformation: -2
  2 dimensions 2000 events 2 invarients 1 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Grid read from file
 Grid defined OK
 Set CM energy to 900.00
 Mapping Graph 1 to config 1
 Determine nb_t
 T-channel found: 0
 No modification is setgrid, grid read from file
   1 1 2
 Masses: 0.450E+03 0.450E+03 0.000E+00 0.000E+00
 Got moffset 6
Using random seed offsets 1 : 1
  with seed 24
 Ranmar initialization seeds 31059 28340
  Particle 3 4
      Et > 0.0 0.0
       E > 0.0 0.0
     Eta < -1.0 -1.0
   xqcut: 0.0 0.0
d R # 3 > -0.0 0.0
s min # 3> 0.0 0.0
xqcutij # 3> 0.0 0.0
 alpha_s for scale 900.00018983999996 is 8.8077955047643222E-002
  ==========================================================================================
 { }
 {   }
 {  ,,  }
 { `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.9.4 (2021-05-30), Ref: arXiv:1103.0621v2, arXiv:1405.0301  }
 {   }
 { }
  ==========================================================================================
 ===============================================================
 INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
 ===============================================================
  > MLReductionLib = 7|6|1|3
  > CTModeRun = -1
  > MLStabThres = 1.0000000000000000E-002
  > 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
 ===============================================================
 ##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.
 INFO:: Attempting to rescue the precision improvement with an alternative method.

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

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

 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Attempting to rescue the precision improvement with an alternative method.
 INFO:: Further warnings from the improve_ps routine will now be supressed.
 Iteration 1 Mean: 0.1485E-01 Abs mean: 0.1485E-01 Fluctuation: 0.312E-03 0.676E-01 100.0%
  1 0.1485E-01 0.1485E-01 +- 0.3121E-03 0.94
 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.021 0.000
 Finished due to accuracy 0.0000000000000000 2.0000000000000000E-002

 -------------------------------------------------------------------------------
 Accumulated results: Integral = 0.1485E-01
                        Std dev = 0.3121E-03
                       Cross sec = 0.1485E-01
             Chi**2 per DoF. = 0.0000
 -------------------------------------------------------------------------------
 Found 739 events.
 Wrote 712 events.
 Events wgts > 1: 76
 % Cross section > 1: 1.0643698562543149E-002 0.95987816666105397
-------------------------------------------------
---------------------------
 Results Last 1 iters: Integral = 0.1485E-01
                     Abs integral = 0.1485E-01
                          Std dev = 0.3121E-03
                  Chi**2 per DoF. = 0.0000
-------------------------------------------------
---------------------------
 Status 2.0000000000000000E-002 2 1

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

From the log, I do not see any issue at all.
 @danielle can you try to reproduce the issue?

Hi Taylor Murphy,

I tested again the complete example of A.4 of 2107.04598 with two versions of MG5_aMC:
- v2.9.3
- v2.9.6 (latest version)
and everything is working fine as expected.

So, my suggestion is to update to the latest MG5_aMC version and see if the problem persists.

Cheers.
--
Daniele

Hi all,

I have tried again using quite a few configurations. I continue to get the same error as in the original post, which is really bizarre.

Specifically, I have used

MG5 v2.9.3 in Fedora 31 (admittedly outdated OS)
MG5 v2.9.6 in Fedora 31
MG5 v2.9.6 in Fedora 35 (fully up to date)
MG5 v2.9.6 in Ubuntu 20.04 (virtual machine within Windows 10... completely independent of the previous instances).

I am installing MadDM v3.2 each time via "install maddm" within MG5_aMC.

My results are exactly the same every time.

I am sorry to prolong our exchange, but I wonder if you can share some details about your setup (maybe OS, compilers, anything relevant).

Hi Taylor Murphy,

my setup is:
- Ubuntu 20.04
- gcc 9.3.0
- python 2.7.18

Could you please check that you have installed all the looptools, including Collier and Ninja?
Unfortunately I am not able to tell you much more. A test you could do is to increase the value of 'vave_indirect_line' from 7.5e-4 to 1e-3 or 1e-2 for instance and see if something changes.
I think @olivier is better suited to help you here as it seems to me this is more MG5_aMC related.

Cheers.
--
Daniele

Hi Daniele,

You found the issue. My run proceeds as normal with larger values of vave_indirect_line, including 1e-3 and 9e-4.

It fails for 8e-4 and smaller, hence the failure when I tried to run according to the settings in 2107.04598.

Any idea what might cause this behavior or how to mitigate it? And how serious an issue is it? Certainly not all choices for mean velocity at the galactic center are valid.

But thank you for suggesting something that reveals the problem.

Hi Taylor Murphy,

the mean velocity at the galactic centre was taken from 1810.09466, fig. 3. We estimated a range [200, 250] km/s, which corresponds to [5e-4, 1e-3]c, so I think you are 'safe' (from a phenomenological point of view) if you stay below 1e-3 c.
Remember that if you are dealing with p-wave processes, their <sigma v> will depend on the velocity.

This certainly does not explain why you obtain this issue (and why I don't). This is likely due to some instabilities during the computation, when the two incoming particles have very similar energies. We will try to inspect this.

Cheers.
--
Daniele