MadGraph5_aMC@NLO

launch ends with non zero status if req_acc_fo is not set to -1

Asked by Prudhvi Bhattiprolu

Dear Madgraph team,

I am trying to compute production cross-sections (for vector-like leptons) at NLO at fixed order. And, in that process, I get the following error:

*******
Command "launch auto " interrupted with error:
Exception : program /nfs/pt3nfs/xdata/prudhvib/VLL/xsec/13tev/ERROR/@@output@@/SubProcesses/P0_uux_nupnupx/ajob1 1 all 0 0 launch ends with non zero status: 1. Stop all computation
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/nfs/pt3nfs/xdata/prudhvib/VLL/xsec/13tev/ERROR/@@output@@/run_01_tag_1_debug.log'.
Please attach this file to your report.
*******

I am getting this error for the default setting of the “req_acc_fo” option in run_card (i.e. req_acc_fo = 0.01). However, this error did not occur if I set “req_acc_fo” to “-1” and instead use the default setting of "npoints_fo_grid=5000, niters_fo_grid=4, npoints_fo=10000, niters_fo=6". I am attaching the log file here (https://www.dropbox.com/s/dz7d0tqx5ldx9nz/run_01_tag_1_debug.log?dl=0). I am surprised that 1% accuracy seems to be very high. Also, this problem does not occur for all processes, but only for one. Am I doing something wrong? Could you please throw some light on this?

Please let me know if anything is needed from my end.

Thank you very much!

Cheers,
Prudhvi

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
marco zaro
Solved:
Last query:
Last reply:
Revision history for this message
Launchpad Janitor (janitor) said : 		#1

This question was expired because it remained in the 'Open' state without activity for the last 15 days.

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

Hi Prudhvi,
what is the process syntax and the model you use?
also, when you get this error

Exception : program /nfs/pt3nfs/xdata/prudhvib/VLL/xsec/13tev/ERROR/@@output@@/SubProcesses/P0_uux_nupnupx/ajob1 1 all 0 0 launch ends with non zero status: 1. Stop all

can you get the content of P0_uux_nupnupx/all_G1/log.txt?

best wishes,

marco

Revision history for this message
Prudhvi Bhattiprolu (prudhvibhattiprolu) said : 		#3

Hi Marco,

Sorry for the late reply.

The Madgraph model file, for vector-like leptons (NLO in QCD), that I've used is attached here (https://www.dropbox.com/sh/r17xmd81io9bp2p/AABqgAlidNb0InSRvBJuP4aCa?dl=0). Using this model file, I was trying to calculate cross-section for the process "generate p p > nup nup~ [QCD]". And, I was getting the error, that I've mentioned in my question, once I do "launch".

However, this error seems to go away if I set “req_acc_fo” to “-1” and instead use the default setting of "npoints_fo_grid=5000, niters_fo_grid=4, npoints_fo=10000, niters_fo=6".

Also, the content of P0_uux_nupnupx/all_G1/log.txt is shown below:

*************
 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11790000000000000
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model. Hence it will be set to: 5
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 MU_R = 91.188000000000002
 mdl_gtaupeRz = 0.0000000000000000
 mdl_gtaupmuRz = 0.0000000000000000
 mdl_gtauptaRz = 0.10000000000000001
 mdl_gtaupeRh = 0.0000000000000000
 mdl_gtaupmuRh = 0.0000000000000000
 mdl_gtauptaRh = 0.10000000000000001
 mdl_gnupeRw = 0.0000000000000000
 mdl_gnupmRw = 0.0000000000000000
 mdl_gnuptRw = 0.10000000000000001
 mdl_Tl = -0.50000000000000000
 aEWM1 = 127.95200000000000
 mdl_Gf = 1.1663789999999999E-005
 aS = 0.11790000000000000
 mdl_ymt = 172.75999999999999
 mdl_MZ = 91.187600000000003
 mdl_MT = 172.75999999999999
 mdl_MH = 125.09999999999999
 mdl_MTAUP = 300.00000000000000
 mdl_MNUP = 300.00000000000000
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WT = 1.4199999999999999
 mdl_WH = 4.1700000000000001E-003
 mdl_WTAUP = 1.0000000000000000
 mdl_WNUP = 1.0000000000000000
  Internal Params
  ---------------------------------

 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15650.009999999998
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_MT__exp__2 = 29846.017599999996
 mdl_aEW = 7.8154307865449541E-003
 mdl_MW = 79.831609916233603
 mdl_sqrt__aEW = 8.8404925126063849E-002
 mdl_ee = 0.31338729995741132
 mdl_MW__exp__2 = 6373.0859418176869
 mdl_sw2 = 0.23355992619469568
 mdl_cw = 0.87546563256663845
 mdl_sqrt__sw2 = 0.48328038051910993
 mdl_sw = 0.48328038051910993
 mdl_g1 = 0.35796642186700228
 mdl_gw = 0.64845856068228980
 mdl_vev = 246.21961912951551
 mdl_gnupb = -0.17898321093350114
 mdl_gtaupb = -0.17898321093350114
 mdl_gtaupnupw3 = -0.32422928034114490
 mdl_gwtaupnup = 0.45852944557691550
 mdl_vev__exp__2 = 60624.100844283676
 mdl_lam = 0.12907416177765593
 mdl_yt = 0.99228297037962632
 mdl_muH = 88.459058326437088
 mdl_I2a33 = ( 0.99228297037962632 , 0.0000000000000000 )
 mdl_I3a33 = ( 0.99228297037962632 , 0.0000000000000000 )
 mdl_ee__exp__2 = 9.8211599774596495E-002
 mdl_sw__exp__2 = 0.23355992619469568
 mdl_cw__exp__2 = 0.76644007380530443
 mdl_yt__exp__2 = 0.98462549330541438
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34336569426778790
 mdl_G__exp__2 = ( 1.4815750954329465 , 0.0000000000000000 )
 mdl_G__exp__4 = ( 2.1950647634071445 , 0.0000000000000000 )
 mdl_G__exp__3 = ( 1.8033727530176447 , 0.0000000000000000 )
  Couplings of VLLD_NLO
  ---------------------------------

   R2GC_100_1 = 0.00000E+00 -0.63803E-02
   R2GC_101_2 = 0.00000E+00 0.78231E-02
  R2GC_125_18 = 0.00000E+00 0.26136E-02
  R2GC_129_20 = -0.00000E+00 -0.52271E-02
  R2GC_168_37 = 0.00000E+00 0.28855E-02
   R2GC_87_45 = -0.00000E+00 -0.14428E-02
 UVGC_106_2_1eps -0.00000E+00 -0.13068E-02
 UVGC_108_3_1eps 0.00000E+00 0.26136E-02
 UVGC_125_21_1ep 0.00000E+00 0.13068E-02
 UVGC_129_28_1ep -0.00000E+00 -0.26136E-02
        GC_12 = 0.00000E+00 0.12172E+01
         GC_1 = -0.00000E+00 -0.10446E+00
         GC_2 = 0.00000E+00 0.20892E+00
        GC_45 = 0.00000E+00 0.57666E-01
        GC_46 = -0.00000E+00 -0.11533E+00
        GC_49 = -0.00000E+00 -0.31268E+00
        GC_50 = 0.00000E+00 0.25502E+00
        GC_56 = 0.00000E+00 0.37035E+00
        GC_58 = -0.00000E+00 0.00000E+00

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

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23nlo
 alpha_s(Mz)= 0.1190 running at 2 loops.
 alpha_s(Mz)= 0.1190 running at 2 loops.
 Renormalization scale fixed @ 300.00000000000000
 Factorization scales fixed @ 300.00000000000000 300.00000000000000

 Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 getting user params
 Number of phase-space points per iteration: -1
 Maximum number of iterations is: 6
 Desired accuracy is: 5.0000000000000003E-002
 Using adaptive grids: 2
 Using Multi-channel integration
 Do MC over helicities for the virtuals
 Number of channels to integrate together: 2
 Running Configuration Number(s): 1 2
 Splitting channel: 0
 Weight multiplier: 1.0000000000000000
 doing the all of this channel
 Normal integration (Sfunction != 1)
 RESTART: Fresh run
 about to integrate 7 -1 6
 imode is 0
channel 1 : 1 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
channel 2 : 2 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
STOP 1
 ------- iteration 1
 Update # PS points (even): 280 --> 256
Using random seed offsets: 0 , 1 , 0
  with seed 1000
 Ranmar initialization seeds 17562 10414
 Total number of FKS directories is 4
 For the Born we use nFKSprocesses # 2 0
tau_min 1 1 : 0.60000E+03 -- 0.60000E+03
tau_min 2 1 : 0.60000E+03 -- 0.60000E+03
tau_min 3 1 : 0.60000E+03 -- 0.60000E+03
tau_min 4 1 : 0.60000E+03 -- 0.60000E+03
 bpower is 0.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.300000D+03 0.300000D+03 1.00
 muF1, muF1_reference: 0.300000D+03 0.300000D+03 1.00
 muF2, muF2_reference: 0.300000D+03 0.300000D+03 1.00
 QES, QES_reference: 0.300000D+03 0.300000D+03 1.00

 muR_reference [functional form]:
    fixed
 muF1_reference [functional form]:
    fixed
 muF2_reference [functional form]:
    fixed
 QES_reference [functional form]:
    fixed scale

 alpha_s= 0.10080414181095018
 ERROR: INTEGRAL APPEARS TO BE ZERO.
 TRIED 100096 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
Time in seconds: 7

*************

Thank you very much!

Best,
Prudhvi

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

Dear Prudhvi,
it looks like one of the two born diagrams (specifically the one with the photon) is zero. This creates problems during the integration. Can you either set the couplings so that it is tiny but non-zero, or apply a restrict card to the model to have it filtered out at generation?

Let me know

cheers,

marco

Revision history for this message
Prudhvi Bhattiprolu (prudhvibhattiprolu) said : 		#5

Thanks marco zaro, that solved my question.

Revision history for this message
Prudhvi Bhattiprolu (prudhvibhattiprolu) said : 		#6

Hi Marco,

Thank you very much for all your help! The issue is that I had previously written the production Lagrangian of "nup" particle (in the FeynRules' .fr file) in terms of Wi[mu,3] and B[mu] instead of writing it in terms of only Z[mu] (as A[mu] piece vanishes for "nup"). So, Madgraph was also dealing with the born diagram with a photon, even though it was zero, and therefore giving rise to this issue that I had reported. Once I recast the relevant piece of the Lagrangian in terms of physical fields (here, Z[mu]) instead of the unphysical, this issue has disappeared.

Thanks again!

Best,
Prudhvi