ttZ@NLO aMCatNLOError : An error occurred during the collection of results.

Hi all,

I've run into a problem running ttbar+Z production at NLO. Maybe I'm just doing something stupid and I'm sure this worked in the past but I've run out of ideas.

I get this error:

ERROR: not all jobs terminated correctly

INFO: Return code of the event collection: None
INFO: Output of the event collection:
N of directories: 24
8 jobs did not terminate correctly:
P0_gg_zttx/all_G1
P0_gg_zttx/all_G2
P0_gg_zttx/all_G3
P0_gg_zttx/all_G4
P0_gg_zttx/all_G5
P0_gg_zttx/all_G6
P0_gg_zttx/all_G7
P0_gg_zttx/all_G8
The results are probably not correct. Please check the relevant log files corresponding to the above jobs for more information.
Found 16 correctly terminated jobs
16 24
P0_uxu_zttx/all_G4/res_0: Final result: 2.24602897282302309E-002 +/- 5.51327402591432524E-004

Error detected in "calculate_xsect "
write debug file /afs/cern.ch/work/m/mcfayden/mcgen/MadGraph5_aMCatNLO/ttV_xs/13TeV/ttZ-NLO-14TeV/run_03_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
aMCatNLOError : An error occurred during the collection of results.
 Please check the .log files inside the directories which failed.
quit
INFO:

An example log file is:

 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.00000000000000008E-005
  > PrecisionVirtualAtRunTime = 1.00000000000000002E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.00000000000000010E-003
 ===============================================================
 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
 ==== STAND-ALONE LHAGLUE MODE TO USE LHAPDF ====
 *************************************
 * LHAPDF Version 5.8.8 *
 * Configured for the following: *
 * All PDFs *
 * LOW MEMORY option *
 * Maximum 3 concurrent set(s) *
 *************************************

 >>>>>> PDF description: <<<<<<
 CTEQ6.6
 Reference:
 P.M.Nadolsky, H.-L.Lai, Q.-H.Cao, J. Huston,
 J.Pumplin, D.Stump, W.-K.Tung and C.-P.Yuan.
 arXiv:0802.0007
 This set has 44 member PDFs.
     mem=0 --> central value
     mem=1-44 --> eigenvector sets
         use these for standard uncertainty calculations
         all have alpha_s(mz) = 0.118
 >>>>>> <<<<<<

 ==============================================
 PDFset name ../../../lib/PDFsets/cteq66.LHgrid
 with 44 members
 ==== initialized. ===========================
 Strong coupling at Mz for PDF is: 0.11798
 New value of alpha_s from PDF lhapdf : 0.11798150967350456
 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.16639000000000003E-005
 aS = 0.11799999999999999
 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 = 0.0000000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 2.0476000000000001
 mdl_WH = 6.38233899999999989E-003
  Internal Params
  ---------------------------------

 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_CKM22 = 1.0000000000000000
 mdl_I4x33 = 0.0000000000000000
 mdl_I1x33 = 0.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM22 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 mdl_complexi = ( 0.0000000000000000 , 1.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_Ncol__exp__2_m_1 = 8.0000000000000000
 mdl_MT__exp__2 = 29929.000000000000
 mdl_Ncol__exp__2_m_1_0 = 8.0000000000000000
 mdl_aEW = 7.54677111397888348E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.68721538467815552E-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.02066170006547170E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.54305075882732989E-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_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.99366614581500623 , 0.0000000000000000 )
 mdl_Axial_tbGp = (-0.99366614581500623 , -0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.99366614581500623 , 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.48355227599988748E-002
 mdl_sw__exp__2 = 0.22224648578577769
 mdl_yt__exp__2 = 0.98737240933884918
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_MU_R__exp__2 = 8315.2513440000002
 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 loop_sm-no_b_mass
  ---------------------------------

       UV_3Gt = 0.48815E-02 0.00000E+00
      UV_GQQt = 0.00000E+00 -0.48815E-02
     UV_tMass = 0.00000E+00 0.34177E+00
   UVWfct_t_0 = -0.98778E-03 -0.00000E+00
   UVWfct_G_2 = 0.40088E-02 0.00000E+00
    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_GGZup = -0.11592E-02 0.00000E+00
   R2_GGZdown = 0.11592E-02 0.00000E+00
 R2_GGGZvecUp = 0.00000E+00 0.28750E-03
 R2_GGGZvecDown -0.00000E+00 -0.49665E-03
 R2_GGGZaxialUp 0.00000E+00 0.63522E-02
 R2_GGGZaxialDow -0.00000E+00 -0.63522E-02
         GC_4 = -0.12177E+01 0.00000E+00
         GC_5 = 0.00000E+00 0.12177E+01
         GC_6 = 0.00000E+00 0.14828E+01
       R2_3Gq = 0.76230E-02 0.00000E+00
       R2_3Gg = 0.31445E-01 0.00000E+00
       R2_GQQ = -0.00000E+00 -0.30492E-01
       R2_GGq = 0.00000E+00 0.62601E-02
       R2_GGt = -0.00000E+00 -0.11242E+04
     R2_GGg_1 = 0.00000E+00 0.28170E-01
     R2_GGg_2 = -0.00000E+00 -0.18780E-01
       R2_QQq = 0.00000E+00 0.12520E-01
       R2_QQt = 0.00000E+00 0.43320E+01
  UV_3Gg_1eps = 0.62890E-01 0.00000E+00
  UV_3Gb_1eps = -0.38115E-02 0.00000E+00
 UV_GQQg_1eps = 0.00000E+00 -0.62890E-01
 UV_GQQq_1eps = 0.00000E+00 0.38115E-02
 UV_tMass_1eps = 0.00000E+00 0.64980E+01
 UVWfct_t_0_1eps -0.18780E-01 0.00000E+00
 UVWfct_G_2_1eps -0.31300E-02 0.00000E+00
        GC_21 = -0.00000E+00 -0.28804E+00
        GC_22 = 0.00000E+00 0.28804E+00
        GC_23 = -0.00000E+00 -0.27437E-01

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

 Running at P P machine @ 14000.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 fixed @ 212.69999999999999
 Factorization scales fixed @ 212.69999999999999 212.69999999999999

 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
 Diagram information for clustering has been set-up for nFKSprocess 7
 Diagram information for clustering has been set-up for nFKSprocess 8
 Diagram information for clustering has been set-up for nFKSprocess 9
 Diagram information for clustering has been set-up for nFKSprocess 10
 Diagram information for clustering has been set-up for nFKSprocess 11
 Diagram information for clustering has been set-up for nFKSprocess 12
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 6
Enter desired accuracy:
 Desired absolute accuracy: 0.10000000000000001
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)?
 Summing over 1 helicities/event for virt
Enter Configuration Number:
Running Configuration Number: 1
 enter id string for this run
 enter 1 if you want restart files
 enter 0 to exclude, 1 for new run, 2 to restart
 '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 10 -1 6 1
 imode is 0
 ------- iteration 1
 Update # PS points (even): 800 --> 800
Using random seed offsets: 1 , 1 , 0
  with seed 34
 Ranmar initialization seeds 13168 9409
 Total number of FKS directories is 12
 For the Born we use nFKSprocesses # 2 4
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 1. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 2. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 5. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 5. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 5. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 6. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 6. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 6. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 7. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 7. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 7. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 8. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 8. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 8. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 9. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 9. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 9. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 10. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 10. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 10. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 11. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 11. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 11. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 12. Absolute lower bound for tau at the Born is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 12. Lower bound for tau is 0.97517E-03 0.43719E+03 0.14000E+05
nFKSprocess: 12. Lower bound for tau is (taking resonances into account) 0.97517E-03 0.43719E+03 0.14000E+05
 bpower is 2.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.212700D+03 0.212700D+03 1.00
 muF1, muF1_reference: 0.212700D+03 0.212700D+03 1.00
 muF2, muF2_reference: 0.212700D+03 0.212700D+03 1.00
 QES, QES_reference: 0.212700D+03 0.212700D+03 1.00

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

 alpha_s= 0.10464527456295308
 alpha_s value used for the virtuals is (for the first PS point): 0.10464527456295308
  ==========================================================================================
 { }
 { }
 { ,, }
 { `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.2.1 (2014-09-25), Ref: arXiv:1103.0621v2, arXiv:1405.0301 }
 { }
 { }
  ==========================================================================================
 ===============================================================
 INFO: MadLoop read these parameters from ../MadLoop5_resources/MadLoopParams.dat
 ===============================================================
  > MLReductionLib = 1|4|3|2
  > IREGIMODE = 2
  > IREGIRECY = T
  > CTModeRun = -1
  > MLStabThres = 1.00000000000000002E-003
  > NRotations_DP = 1
  > NRotations_QP = 0
  > CTStabThres = 1.00000000000000002E-002
  > CTLoopLibrary = 2
  > CTModeInit = 1
  > CheckCycle = 3
  > MaxAttempts = 10
  > UseLoopFilter = F
  > ImprovePSPoint = 2
  > DoubleCheckHelicityFilter = T
  > LoopInitStartOver = F
  > HelInitStartOver = F
  > ZeroThres = 1.00000000000000006E-009
  > OSThres = 1.00000000000000002E-008
  > WriteOutFilters = 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.4 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 02-01-2014 #
# #
# 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. #
# #
########################################################################
 #####################################################################################
 # #
 # IREGI-alpha-1.0.0 #
 # package for one-loop tensor Integral REduction with General propagator Indices #
 # Author: Hua-Sheng Shao (<email address hidden><email address hidden>) #
 # a) Physics School, Peking University, Beijing, China #
 # b) PH Department, TH Unit, CERN, Geneva, Switzerland #
 # #
 #####################################################################################
 alpha_s value used for the virtuals is (for the first PS point): 0.10464527456295308
 ---- POLES CANCELLED ----
 ERROR #51 in dsig: NaN skipping event
 ERROR #51 in dsig: NaN skipping event
 ERROR #51 in dsig: NaN skipping event
 ERROR #51 in dsig: NaN skipping event
 ERROR #51 in dsig: NaN skipping event
 ERROR #51 in dsig: NaN skipping event
 ERROR #51 in dsig: NaN skipping event
ABS integral = NaN +/- NaN ( NaN %)
Integral = NaN +/- NaN ( NaN %)
Virtual = NaN +/- NaN ( NaN %)
Virtual ratio = NaN +/- NaN ( NaN %)
ABS virtual = NaN +/- NaN ( NaN %)
Born*ao2pi = 0.1764E-03 +/- 0.1627E-04 ( 9.223 %)
Chi^2 per d.o.f. 0.0000E+00
accumulated results ABS integral = NaN +/- NaN ( NaN %)
accumulated results Integral = NaN +/- NaN ( NaN %)
accumulated results Virtual = NaN +/- NaN ( NaN %)
accumulated results Virtual ratio = NaN +/- NaN ( NaN %)
accumulated results ABS virtual = NaN +/- NaN ( NaN %)
accumulated results Born*ao2pi = 0.1764E-03 +/- 0.1627E-04 ( 9.223 %)
accumulated result Chi^2 per DoF = 0.0000E+00
update virtual fraction to: 1.000 NaN
  1: 0 1 2 3 5 67 9 1 2
 ------- iteration 2
 Update # PS points (even): 1600 --> 1600

Program received signal 11 (SIGSEGV): Segmentation fault.

Backtrace for this error:
  + /lib64/libc.so.6(+0x326a0) [0x7f4a6c0a26a0]
  + function get_ave_virt_ (0x48F380)
  + function mint_ (0x490A32)
  + in the main program
    from file driver_mintFO.f
  + /lib64/libc.so.6(__libc_start_main+0xfd) [0x7f4a6c08ed5d]
Time in seconds: 1116

The process card is:

import model loop_sm-no_b_mass
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > t t~ z [QCD]
output ttZ-NLO

Run card:

#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
#***********************************************************************
#
#*******************
# Running parameters
#*******************
#
#***********************************************************************
# Tag name for the run (one word) *
#***********************************************************************
  tag_1 = run_tag ! name of the run
#***********************************************************************
# Number of events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#***********************************************************************
 50000 = nevents ! Number of unweighted events requested
    -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
    -1 = nevt_job! Max number of events per job in event generation.
                 ! (-1= no split).
average = event_norm ! Normalize events to sum or average to the X sect.
#***********************************************************************
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#***********************************************************************
 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
                     ! number of points and iter. below)
# These numbers are ignored except if req_acc_FO is equal to -1
 5000 = npoints_FO_grid ! number of points to setup grids
 4 = niters_FO_grid ! number of iter. to setup grids
 10000 = npoints_FO ! number of points to compute Xsec
 6 = niters_FO ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
     0 = iseed ! rnd seed (0=assigned automatically=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
    1 = lpp1 ! beam 1 type (0 = no PDF)
    1 = lpp2 ! beam 2 type (0 = no PDF)
 6500 = ebeam1 ! beam 1 energy in GeV
 6500 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
 lhapdf = pdlabel ! PDF set
  10550 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# 21141 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#***********************************************************************
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#***********************************************************************
  HERWIG6 = parton_shower
#***********************************************************************
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses of all final state particles and partons. This *
# can be changed in SubProcesses/set_scales.f) *
#***********************************************************************
 T = fixed_ren_scale ! if .true. use fixed ren scale
 T = fixed_fac_scale ! if .true. use fixed fac scale
 212.7 = muR_ref_fixed ! fixed ren reference scale
 212.7 = muF1_ref_fixed ! fixed fact reference scale for pdf1
 212.7 = muF2_ref_fixed ! fixed fact reference scale for pdf2
# 172.5 = muR_ref_fixed ! fixed ren reference scale
# 172.5 = muF1_ref_fixed ! fixed fact reference scale for pdf1
# 172.5 = muF2_ref_fixed ! fixed fact reference scale for pdf2
#***********************************************************************
# Renormalization and factorization scales (advanced and NLO options) *
#***********************************************************************
 T = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
 212.7 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
# 172.5 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
 1 = muR_over_ref ! ratio of current muR over reference muR
 1 = muF1_over_ref ! ratio of current muF1 over reference muF1
 1 = muF2_over_ref ! ratio of current muF2 over reference muF2
 1 = QES_over_ref ! ratio of current QES over reference QES
#***********************************************************************
# Reweight flags to get scale dependence and PDF uncertainty *
# For scale dependence: factor rw_scale_up/down around central scale *
# For PDF uncertainty: use LHAPDF with supported set *
#***********************************************************************
# .true. = reweight_scale ! reweight to get scale dependence
 .false. = reweight_scale ! reweight to get scale dependence
  0.5 = rw_Rscale_down ! lower bound for ren scale variations
  2.0 = rw_Rscale_up ! upper bound for ren scale variations
  0.5 = rw_Fscale_down ! lower bound for fact scale variations
  2.0 = rw_Fscale_up ! upper bound for fact scale variations
# .true. = reweight_PDF ! reweight to get PDF uncertainty
 .false. = reweight_PDF ! reweight to get PDF uncertainty
# 21142 = PDF_set_min ! First of the error PDF sets
# 21181 = PDF_set_max ! Last of the error PDF sets
  10551 = PDF_set_min ! First of the error PDF sets
  10594 = PDF_set_max ! Last of the error PDF sets
#***********************************************************************
# Merging - WARNING! Applies merging only at the hard-event level. *
# After showering an MLM-type merging should be applied as well. *
# See http://amcatnlo.cern.ch/FxFx_merging.htm for more details. *
#***********************************************************************
 0 = ickkw ! 0 no merging, 3 FxFx merging
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) *
#***********************************************************************
 15 = bwcutoff
#***********************************************************************
# Cuts on the jets *
# When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# (more specific cuts can be specified in SubProcesses/cuts.f) *
#***********************************************************************
   1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
 0.7 = jetradius ! The radius parameter for the jet algorithm
  10 = ptj ! Min jet transverse momentum
  -1 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
#***********************************************************************
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# (more specific gen cuts can be specified in SubProcesses/cuts.f) *
#***********************************************************************
   0 = ptl ! Min lepton transverse momentum
  -1 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
   0 = drll ! Min distance between opposite sign lepton pairs
  30 = mll ! Min inv. mass of all oppositely charged lepton pairs
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
#***********************************************************************
  20 = ptgmin ! Min photon transverse momentum
  -1 = etagamma ! Max photon abs(pseudo-rap)
 0.4 = R0gamma ! Radius of isolation code
 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
 .true. = isoEM ! isolate photons from EM energy (photons and leptons)
#***********************************************************************
# maximal pdg code for quark to be considered as a jet *
#***********************************************************************
 5 = maxjetflavor
#***********************************************************************

Cheers,

Josh.