MadGraph5_aMC@NLO

launch ends with non zero status: 1. Stop all computation

Asked by David Vannerom

Dear all,

I am trying to generate events for a NLO VBF process (ZToNuNu+Jets). I first did:

./bin/mg5_aMC process_card.dat

This created a folder name after the process, and now I am trying to generate events with:

cd process_folder
./bin/generate_events

But at some point it fails witht he following warning:

INFO: Idle: 2873, Running: 48, Completed: 103 [ 46m 59s ]
INFO: Idle: 2872, Running: 48, Completed: 104 [ 47m 24s ]
INFO: Idle: 2871, Running: 48, Completed: 105 [ 47m 26s ]
INFO: Idle: 2870, Running: 48, Completed: 106 [ 47m 32s ]
WARNING: program /scratch/9121102.cream02.iihe.ac.be/ZToNuNu_VBF_5f_NLO_FXFX_pT250To400/SubProcesses/P0_ub_vevexbc/ajob1 1 F 0 0 launch ends with non zero status: 1. Stop all computation
INFO: Idle: 2870, Running: 47, Completed: 107 [ 47m 35s ]
INFO: Idle: 2870, Running: 46, Completed: 108 [ 47m 35s ]
INFO: Idle: 2870, Running: 45, Completed: 109 [ 47m 35s ]
INFO: Idle: 2870, Running: 44, Completed: 110 [ 47m 35s ]

Since I did not have any log file, I followed a proposition on a similar thread so I have done the following thing:

cd SubProcesses/P0_ub_vevexbc/
chmod +x ./ajob1
./ajob1 1 F 0 0
vim GF1/log.txt

The log then reads:

==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.1.6 loading /cvmfs/cms.cern.ch/slc6_amd64_gcc530/external/lhapdf/6.1.6-giojec/share/LHAPDF/NNPDF30_nlo_nf_5_pdfas/NNPDF30_nlo_nf_5_pdfas_0000.dat
NNPDF30_nlo_nf_5_pdfas PDF set, member #0, version 2; LHAPDF ID = 292200
 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = -1.0000000000000000
  > PrecisionVirtualAtRunTime = -1.0000000000000000
  > 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.11799999999999999
 New value of alpha_s from PDF lhapdf : 0.11800222660863767
 using LHAPDF
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 MU_R = 91.188000000000002
 aEWM1 = 132.50700000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_lamWS = 0.22530000000000000
 mdl_AWS = 0.80800000000000005
 mdl_rhoWS = 0.13200000000000001
 mdl_etaWS = 0.34100000000000003
 mdl_ymt = 173.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MZ = 91.188000000000002
 mdl_MH = 125.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 2.0476000000000001
 mdl_WH = 6.3823389999999999E-003
  Internal Params
  ---------------------------------

 mdl_CKM33 = 1.0000000000000000
 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_I4x33 = 0.0000000000000000
 mdl_I1x33 = 0.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 mdl_lamWS__exp__2 = 5.0760090000000001E-002
 mdl_CKM11 = ( 0.97461995499999998 , 0.0000000000000000 )
 mdl_CKM12 = ( 0.22530000000000000 , 0.0000000000000000 )
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_lamWS__exp__3 = 1.1436248277000001E-002
 mdl_CKM13 = ( 1.2197444962317122E-003, -3.1510066152652563E-003)
 mdl_CKM21 = (-0.22530000000000000 , 0.0000000000000000 )
 mdl_CKM22 = ( 0.97461995499999998 , 0.0000000000000000 )
 mdl_CKM23 = ( 4.1014152720000002E-002, 0.0000000000000000 )
 mdl_CKM31 = ( 8.0207441115842885E-003, -3.1510066152652563E-003)
 mdl_CKM32 = ( -4.1014152720000002E-002, 0.0000000000000000 )
 mdl_MZ__exp__2 = 8315.2513440000002
 mdl_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_CKM1x1 = ( 0.97461995499999998 , 0.0000000000000000 )
 mdl_CKM1x2 = ( 0.22530000000000000 , 0.0000000000000000 )
 mdl_CKM1x3 = ( 1.2197444962317122E-003, -3.1510066152652563E-003)
 mdl_CKM2x1 = (-0.22530000000000000 , 0.0000000000000000 )
 mdl_CKM2x2 = ( 0.97461995499999998 , 0.0000000000000000 )
 mdl_CKM2x3 = ( 4.1014152720000002E-002, 0.0000000000000000 )
 mdl_CKM3x1 = ( 8.0207441115842885E-003, -3.1510066152652563E-003)
 mdl_CKM3x2 = ( -4.1014152720000002E-002, 0.0000000000000000 )
 mdl_conjg__CKM1x3 = ( 1.2197444962317122E-003, 3.1510066152652563E-003)
 mdl_conjg__CKM2x3 = ( 4.1014152720000002E-002, -0.0000000000000000 )
 mdl_conjg__CKM2x1 = (-0.22530000000000000 , -0.0000000000000000 )
 mdl_conjg__CKM3x1 = ( 8.0207441115842885E-003, 3.1510066152652563E-003)
 mdl_conjg__CKM2x2 = ( 0.97461995499999998 , -0.0000000000000000 )
 mdl_conjg__CKM3x2 = ( -4.1014152720000002E-002, -0.0000000000000000 )
 mdl_Ncol__exp__2 = 9.0000000000000000
 mdl_MT__exp__2 = 29929.000000000000
 mdl_conjg__CKM11 = ( 0.97461995499999998 , -0.0000000000000000 )
 mdl_conjg__CKM12 = ( 0.22530000000000000 , -0.0000000000000000 )
 mdl_conjg__CKM13 = ( 1.2197444962317122E-003, 3.1510066152652563E-003)
 mdl_conjg__CKM21 = (-0.22530000000000000 , -0.0000000000000000 )
 mdl_conjg__CKM22 = ( 0.97461995499999998 , -0.0000000000000000 )
 mdl_conjg__CKM23 = ( 4.1014152720000002E-002, -0.0000000000000000 )
 mdl_conjg__CKM31 = ( 8.0207441115842885E-003, 3.1510066152652563E-003)
 mdl_conjg__CKM32 = ( -4.1014152720000002E-002, -0.0000000000000000 )
 mdl_aEW = 7.5467711139788835E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.6872153846781555E-002
 mdl_ee = 0.30795376724436879
 mdl_MW__exp__2 = 6467.2159543705357
 mdl_sw2 = 0.22224648578577766
 mdl_cw = 0.88190334743339216
 mdl_sqrt__sw2 = 0.47143025548407230
 mdl_sw = 0.47143025548407230
 mdl_g1 = 0.34919219678733299
 mdl_gw = 0.65323293034757990
 mdl_v = 246.21845810181637
 mdl_v__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yt = 0.99366614581500623
 mdl_ytau = 1.0206617000654717E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.5430507588273299E-002
 mdl_VectorZDown = -0.13030376310310560
 mdl_VectorAUp = 0.20530251149624587
 mdl_VectorADown = -0.10265125574812294
 mdl_VectorWmDxU = 0.23095271737156670
 mdl_AxialWmDxU = -0.23095271737156670
 mdl_VectorWpUxD = 0.23095271737156670
 mdl_AxialWpUxD = -0.23095271737156670
 mdl_I2x13 = ( 7.9699418879263656E-003, 3.1310485988282153E-003)
 mdl_I2x23 = ( -4.0754375057150453E-002, -0.0000000000000000 )
 mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I3x31 = ( 7.9699418879263656E-003, -3.1310485988282153E-003)
 mdl_I3x32 = ( -4.0754375057150453E-002, 0.0000000000000000 )
 mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_Vector_tdGp = ( -7.9699418879263656E-003, -3.1310485988282153E-003)
 mdl_Vector_tsGp = ( 4.0754375057150453E-002, 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.99366614581500623 , 0.0000000000000000 )
 mdl_Axial_tdGp = ( -7.9699418879263656E-003, -3.1310485988282153E-003)
 mdl_Axial_tsGp = ( 4.0754375057150453E-002, 0.0000000000000000 )
 mdl_Axial_tbGp = (-0.99366614581500623 , -0.0000000000000000 )
 mdl_Vector_tdGm = ( 7.9699418879263656E-003, -3.1310485988282153E-003)
 mdl_Vector_tsGm = ( -4.0754375057150453E-002, 0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_Axial_tdGm = ( -7.9699418879263656E-003, 3.1310485988282153E-003)
 mdl_Axial_tsGm = ( 4.0754375057150453E-002, -0.0000000000000000 )
 mdl_Axial_tbGm = (-0.99366614581500623 , 0.0000000000000000 )
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_ee__exp__2 = 9.4835522759998875E-002
 mdl_sw__exp__2 = 0.22224648578577769
 mdl_yt__exp__2 = 0.98737240933884918
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__4 = 2.1987899468922913
 mdl_G__exp__2 = 1.4828317324943823
 mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
 mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
 mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_G__exp__3 = 1.8056676068262196
 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
 mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
  Couplings of loop_sm-ckm_no_b_mass
  ---------------------------------

       R2_DDA = 0.00000E+00 0.25704E-02
       R2_UUA = -0.00000E+00 -0.51409E-02
    R2_DDZ_V2 = 0.00000E+00 0.72127E-02
    R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
    R2_UUZ_V5 = 0.00000E+00 0.68702E-03
      R2_dxcW = 0.00000E+00 0.26059E-02
      R2_sxcW = -0.00000E+00 -0.11273E-01
      R2_bxuW = -0.36445E-04 -0.14108E-04
      R2_uxsW = -0.00000E+00 -0.26059E-02
      R2_uxbW = 0.36445E-04 -0.14108E-04
      R2_cxbW = -0.00000E+00 -0.47438E-03
         GC_5 = 0.00000E+00 0.12177E+01
       R2_QQq = 0.00000E+00 0.12520E-01
         GC_1 = -0.00000E+00 -0.10265E+00
         GC_2 = 0.00000E+00 0.20530E+00
         GC_7 = 0.00000E+00 0.57609E+00
        GC_14 = 0.14555E-02 0.56341E-03
        GC_15 = -0.00000E+00 -0.10407E+00
        GC_17 = 0.00000E+00 0.18945E-01
        GC_21 = -0.00000E+00 -0.28804E+00
        GC_22 = 0.00000E+00 0.28804E+00
        GC_23 = -0.00000E+00 -0.27437E-01
        GC_28 = 0.00000E+00 0.37035E+00
        GC_40 = 0.00000E+00 0.10407E+00
        GC_41 = -0.14555E-02 0.56341E-03
        GC_43 = 0.00000E+00 0.45018E+00
        GC_47 = 0.00000E+00 0.46191E+00

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

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

 Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 Diagram information for clustering has been set-up for nFKSprocess 5
 Diagram information for clustering has been set-up for nFKSprocess 6
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.9999999999999999E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 1
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the all of this channel
 Normal integration (Sfunction != 1)
 Not subdividing B.W.
 about to integrate 13 -1 12 1
 imode is 0
#--------------------------------------------------------------------------
# FastJet release 3.1.3 [fjcore]
# M. Cacciari, G.P. Salam and G. Soyez
# A software package for jet finding and analysis at colliders
# http://fastjet.fr
#
# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].
#
# FastJet is provided without warranty under the terms of the GNU GPLv2.
# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
# and 3rd party plugin jet algorithms. See COPYING file for details.
#--------------------------------------------------------------------------
 ------- iteration 1
 Update # PS points (even): 1040 --> 1040
Using random seed offsets: 1 , 20 , 0
  with seed 33
 Ranmar initialization seeds 11949 9427
 Total number of FKS directories is 6
 FKS process map (sum= 3 ) :
           1 --> 2 : 1 5
           2 --> 2 : 2 6
           3 --> 1 : 3
           4 --> 1 : 4
 ================================
 process combination map (specified per FKS dir):
  1 map 1
  1 inv. map 1
  2 map 1
  2 inv. map 1
  3 map 1
  3 inv. map 1
  4 map 1
  4 inv. map 1
  5 map 1
  5 inv. map 1
  6 map 1
  6 inv. map 1
 ================================
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.00000E+00 0.00000E+00 0.13000E+05
nFKSprocess: 1. Lower bound for tau is 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.00000E+00 0.00000E+00 0.13000E+05
nFKSprocess: 2. Lower bound for tau is 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 5. Absolute lower bound for tau at the Born is 0.00000E+00 0.00000E+00 0.13000E+05
nFKSprocess: 5. Lower bound for tau is 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 5. Lower bound for tau is (taking resonances into account) 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 6. Absolute lower bound for tau at the Born is 0.00000E+00 0.00000E+00 0.13000E+05
nFKSprocess: 6. Lower bound for tau is 0.23669E-05 0.20000E+02 0.13000E+05
nFKSprocess: 6. Lower bound for tau is (taking resonances into account) 0.23669E-05 0.20000E+02 0.13000E+05
 bpower is 0.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.172338D+03 0.172338D+03 1.00
 muF1, muF1_reference: 0.242916D+03 0.242916D+03 1.00
 muF2, muF2_reference: 0.242916D+03 0.242916D+03 1.00
 QES, QES_reference: 0.462910D+03 0.462910D+03 1.00

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

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

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

########################################################################
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #

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

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)
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
########################################################################
########################################################################
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
########################################################################
 ---- POLES CANCELLED ----
Only 1 independent helicities: switching to explicitly summing over them
 ERROR: INTEGRAL APPEARS TO BE ZERO.
 TRIED 100880 PS POINTS AND ONLY 23 GAVE A NON-ZERO INTEGRAND.
Time in seconds: 5

Thanks a lot,
David

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi David,
can you copy-paste the content of your process card?
Ant can you also tell us if you change any cut/parameter in the run/param cards?
Thanks a lot!
Best,

Marco

Revision history for this message
David Vannerom (david.vannerom) said : 		#2

Dear Marco,

Here is the content of the process card:

import model loop_sm-ckm_no_b_mass
#switch to diagonal ckm matrix if needed for speed
#import model loop_sm-no_b_mass

define p = p b b~
define j = j b b~

generate p p > ve ve~ j j QED=4 QCD=0 [QCD]

output $TMPDIR/ZToNuNu_VBF_5f_NLO_FXFX_pT100To250 -nojpeg

As a matter of fact, I have a "_cuts.f" card to implement pT cut for the Z boson.

There is also a FKS card where I have put the "#IRPoleCheckThreshold" and "#PrecisionVirtualAtRunTime" both to "-1.0d0" value as was indicated in some of the threads in this forum for VBF-like processes.

Cheers,
David

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

Hi David,
thanks,
can you try the simpler process where, instead of ve ve~, you put an undecayed Z boson?
Why do you actually need to decay it?
I see you have FxFx in the name of your directory, are you doing any merging?

Thanks

Cheers,

Marco

On 19 Dec 2016, at 09:58, David Vannerom <email address hidden> wrote:

> Question #406808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/406808
>
> Status: Needs information => Open
>
> David Vannerom gave more information on the question:
> Dear Marco,
>
> Here is the content of the process card:
>
> import model loop_sm-ckm_no_b_mass
> #switch to diagonal ckm matrix if needed for speed
> #import model loop_sm-no_b_mass
>
> define p = p b b~
> define j = j b b~
>
> generate p p > ve ve~ j j QED=4 QCD=0 [QCD]
>
> output $TMPDIR/ZToNuNu_VBF_5f_NLO_FXFX_pT100To250 -nojpeg
>
> As a matter of fact, I have a "_cuts.f" card to implement pT cut for the
> Z boson.
>
> There is also a FKS card where I have put the "#IRPoleCheckThreshold"
> and "#PrecisionVirtualAtRunTime" both to "-1.0d0" value as was indicated
> in some of the threads in this forum for VBF-like processes.
>
> Cheers,
> David
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
David Vannerom (david.vannerom) said : 		#4

Dear Marco,

So, we inforce the decay of the Z boson to neutrinos to allow for off-shell Z's, does that make sense to you?

Concerning the FxFx algorithm, we do use it for the jet merging at NLO.

Cheers,
David

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

Dear David,
it makes sense to force the Z to decay into neutrinos, I was wondering if a simpler process (with the Z onshell) was working or not.

For the FxFx, which samples do you want to merge exactly?
Thanks
Cheers,

MArco
On 19 Dec 2016, at 15:33, David Vannerom <email address hidden> wrote:

> Question #406808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/406808
>
> Status: Answered => Open
>
> David Vannerom is still having a problem:
> Dear Marco,
>
> So, we inforce the decay of the Z boson to neutrinos to allow for off-
> shell Z's, does that make sense to you?
>
> Concerning the FxFx algorithm, we do use it for the jet merging at NLO.
>
> Cheers,
> David
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
David Vannerom (david.vannerom) said : 		#6

Hi Marco,

So, I would like to produce parton-level samples that I would shower with a PYTHIA8 script later on, so I guess I do not actually use the FxFx algorithm in the event generation, do I?

Thanks,
David

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

Hi David,
no, if you just use the sample with two jets, you don’t do FxFx.
I was just asking because I saw the name in the output folder.
Anyway, note that FxFx does not support processes which have jets at the lowest multiplicity. (e.g. you cannot do z+ 2 jets and z + 3 jets but you have to start from Z+0 jets)
Cheers,

Marco
On 19 Dec 2016, at 15:47, David Vannerom <email address hidden> wrote:

> Question #406808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/406808
>
> Status: Answered => Open
>
> David Vannerom is still having a problem:
> Hi Marco,
>
> So, I would like to produce parton-level samples that I would shower
> with a PYTHIA8 script later on, so I guess I do not actually use the
> FxFx algorithm in the event generation, do I?
>
> Thanks,
> David
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
David Vannerom (david.vannerom) said : 		#8

Hi Marco,

Ah but this is interesting as I am asking for VBF processes which have at least two jets in the final state:

generate p p > ve ve~ j j QED=4 QCD=0 [QCD]

So what would you suggest as an alternative solution?

Cheers,
David

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

Hi David,
You can do
pp > v v j j QCD=0 QED=4 [QCD], just the un-merged one (that is, without turning FxFx merging on, which may be relevant if you want to combine this process with
pp > v v j j j QCD=0 QED=4 [QCD])

As I said, try to start with

pp > z j j QCD=0 QED=3 [QCD].
If this works, then we may have a look at the more complex one
pp > v v j j QCD=0 QED=4 [QCD]

Cheers,

Marco
On 19 Dec 2016, at 16:02, David Vannerom <email address hidden> wrote:

> Question #406808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/406808
>
> Status: Answered => Open
>
> David Vannerom is still having a problem:
> Hi Marco,
>
> Ah but this is interesting as I am asking for VBF processes which have
> at least two jets in the final state:
>
> generate p p > ve ve~ j j QED=4 QCD=0 [QCD]
>
> So what would you suggest as an alternative solution?
>
> Cheers,
> David
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
David Vannerom (david.vannerom) said : 		#10

Dear Marco,

So just to be sure, if I only generate one process, I should set the ickkw parameter to 0 in the run card?

0 = ickkw ! 0 no merging, 3 FxFx merging

Thanks,
David

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

Hi David,
yes, exactly, leave it to 0
Cheers,

Marco
On 19 Dec 2016, at 16:38, David Vannerom <email address hidden> wrote:

> Question #406808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/406808
>
> Status: Answered => Open
>
> David Vannerom is still having a problem:
> Dear Marco,
>
> So just to be sure, if I only generate one process, I should set the
> ickkw parameter to 0 in the run card?
>
> 0 = ickkw ! 0 no merging, 3 FxFx merging
>
> Thanks,
> David
>
> --
> You received this question notification because you are subscribed to
> the question.

Can you help with this problem?

Provide an answer of your own, or ask David Vannerom for more information if necessary.

To post a message you must log in.