diff --git a/FixedOrderGen/configFO.yml b/FixedOrderGen/configFO.yml
index 2ffcb83..7630ec5 100644
--- a/FixedOrderGen/configFO.yml
+++ b/FixedOrderGen/configFO.yml
@@ -1,95 +1,95 @@
 ## Number of generated events
 events: 200
 
 jets:
   min pt: 20        # minimal jet pt, should be slightly below analysis cut
   peak pt: 30       # peak pt of jets, should be at analysis cut
   algorithm: antikt # jet clustering algorithm
   R: 0.4            # jet R parameter
   max rapidity: 5   # maximum jet rapidity
 
 ## Particle beam
 beam:
   energy: 6500
   particles: [p, p]
 
 ## PDF ID according to LHAPDF
 pdf: 230000
 
 ## Scattering process
-process: p p => h 4j
+process: p p => Wp 4j
 
 ## Particle decays (multiple decays are allowed)
 decays:
   Higgs: {into: [photon, photon], branching ratio: 0.0023568762400521404}
   Wp: {into: [e+, nu_e], branching ratio: 1}
   Wm: {into: [e-, nu_e_bar]} # equivalent to "branching ratio: 1"
 
 ## Fraction of events with two extremal emissions in one direction
 ## that contain an subleading emission e.g. unordered emission
 subleading fraction: 0.05
 
 ## Allow different subleading configurations.
 ## By default all processes are allowed.
 ## This does not check if the processes are implemented in HEJ!
 #
 # subleading channels:
 #   - unordered
 #   - central qqx
 #   - extremal qqx
 
 ## Unweighting parameters
 ## remove to obtain weighted events
 # unweight:
 #   sample size: 200  # should be similar to "events:", but not more than ~10000
 #   max deviation: 0
 
 ## --- The following settings are equivalent to HEJ, --- ##
 ## --- see HEJ for reference.                        --- ##
 
 ## Central scale choice
 scales: max jet pperp
 
 ## Event output files
 event output:
   - HEJFO.lhe
 #  - HEJFO_events.hepmc
 
 ## Analyses
 #
 # analyses:
 ## Rivet analysis
 #   - rivet: MC_XS # rivet analysis name
 #     output: HEJ # name of the yoda files, ".yoda" and scale suffix will be added
 ## Custom analysis
 #   - plugin: /path/to/libmyanalysis.so
 #     my analysis parameter: some value
 
 ## Selection of random number generator and seed
 random generator:
   name: mixmax
 #  seed: 1
 
 ## Vacuum expectation value
 vev: 246.2196508
 
 ## Properties of the weak gauge bosons
 particle properties:
   Higgs:
     mass: 125
     width: 0.004165
   W:
     mass: 80.385
     width: 2.085
   Z:
     mass: 91.187
     width: 2.495
 
 ## Parameters for Higgs-gluon couplings
 ## this requires compilation with QCDloop
 #
 # Higgs coupling:
 #   use impact factors: false
 #   mt: 174
 #   include bottom: true
 #   mb: 4.7
diff --git a/FixedOrderGen/t/CMakeLists.txt b/FixedOrderGen/t/CMakeLists.txt
index ef2d76f..387dfe9 100644
--- a/FixedOrderGen/t/CMakeLists.txt
+++ b/FixedOrderGen/t/CMakeLists.txt
@@ -1,205 +1,209 @@
 set(tst_dir "${CMAKE_CURRENT_SOURCE_DIR}")
 set(runcard_dir ${tst_dir}/runcards)
 
 function(add_long_test)
   set(options NONE)
   set(oneValueArgs NAME)
   set(multiValueArgs COMMAND)
   cmake_parse_arguments(TEST "${options}" "${oneValueArgs}" "${multiValueArgs}" ${ARGN})
   add_test(
     NAME ${TEST_NAME}_short
     COMMAND ${TEST_COMMAND} short
   )
   if(TEST_ALL)
     add_test(
       NAME ${TEST_NAME}
       COMMAND ${TEST_COMMAND}
     )
   endif()
 endfunction()
 
 foreach(tst W_reconstruct_enu W_2j_classify W_nj_classify)
   add_executable(test_${tst} ${tst_dir}/${tst}.cc)
   target_link_libraries(test_${tst} hejfog_lib)
   add_test(NAME ${tst}
     COMMAND test_${tst}
     WORKING_DIRECTORY ${runcard_dir}
   )
 endforeach()
 # this only tests if the runcard actually works, not if the result is correct
 add_test(
+  NAME main_example
+  COMMAND HEJFOG ${PROJECT_SOURCE_DIR}/configFO.yml
+)
+add_test(
   NAME main_2j
   COMMAND HEJFOG ${runcard_dir}/2j.yml
 )
 add_test(
   NAME main_h2j
   COMMAND HEJFOG ${runcard_dir}/h2j.yml
 )
 add_test(
   NAME main_h2j_decay
   COMMAND HEJFOG ${runcard_dir}/h2j_decay.yml
 )
 add_test(
   NAME peakpt
   COMMAND HEJFOG ${runcard_dir}/2j_peak.yml
 )
 # special case where coupling not trivial
 add_test(
   NAME main_W4j_uno+eqqx
   COMMAND HEJFOG ${runcard_dir}/Wp4j_uno+eqqx.yml
 )
 
 # check that uno emission doesn't change FKL xs
 add_executable(FKL_uno FKL_uno.cc)
 target_link_libraries(FKL_uno hejfog_lib)
 add_test(
   NAME FKL_uno
   # calculated with HEJ revision 9570e3809613272ac4b8bf3236279ba23cf64d20
   COMMAND FKL_uno ${runcard_dir}/h3j_uno.yml 0.0243548 0.000119862
 )
 # xs tests
 add_executable(xs_gen xs_gen.cc)
 target_link_libraries(xs_gen hejfog_lib)
 ## Higgs
 add_test(
   NAME xs_hqQ
   # calculated with Sherpa see #132
   COMMAND xs_gen ${runcard_dir}/hqQ.yml 1.612e-02 1.26303e-04
 )
 add_test(
   NAME xs_h2j
   # calculated with HEJ revision 9570e3809613272ac4b8bf3236279ba23cf64d20
   COMMAND xs_gen ${runcard_dir}/h2j.yml 2.04928 0.00956022
 )
 add_test(
   NAME xs_h3j
   # calculated with HEJ revision bd4388fe55cbc3f5a7b6139096456c551294aa31
   COMMAND xs_gen ${runcard_dir}/h3j.yml 1.07807 0.0132409
 )
 add_long_test(
   NAME xs_h3j_uno
   # calculated with HEJ revision 9570e3809613272ac4b8bf3236279ba23cf64d20
   COMMAND xs_gen ${runcard_dir}/h3j_uno.yml 0.00347538 3.85875e-05
 )
 add_test(
   NAME xs_h5j
   # calculated with HEJ revision dbde2ffbb3b383ae6709b2424d8f0f9d5658270b
   COMMAND xs_gen ${runcard_dir}/h5j.yml 0.0112504 0.000199633
 )
 add_test(
   NAME xs_h2j_decay
   # calculated with HEJ revision 9570e3809613272ac4b8bf3236279ba23cf64d20
   COMMAND xs_gen ${runcard_dir}/h2j_decay.yml 0.00466994 2.20995e-05
 )
 ## pure jets
 add_test(
   NAME xs_qQ
   # calculated with Sherpa see #132
   COMMAND xs_gen ${runcard_dir}/qQ.yml 7.354e+05 1.905e+03
 )
 add_test(
   NAME xs_2j
   # calculated with "combined" HEJ svn r3480
   COMMAND xs_gen ${runcard_dir}/2j.yml 86.42031848e06 590570
 )
 add_test(
   NAME xs_3j_uno
   # calculated with HEJ revision 9401196fba75b5d16bc33f2a309175fecaca00f1
   COMMAND xs_gen ${runcard_dir}/3j_uno.yml 900 14.3131
 )
 add_test(
   NAME xs_3j_qqx
   # calculated with HEJ revision 9401196fba75b5d16bc33f2a309175fecaca00f1
   COMMAND xs_gen ${runcard_dir}/3j_qqx.yml 62040 1005
 )
 add_test(
   NAME xs_4j_qqx
   # calculated with HEJ revision 9401196fba75b5d16bc33f2a309175fecaca00f1
   COMMAND xs_gen ${runcard_dir}/4j_qqx.yml 28936 550
 )
 add_test(
   NAME xs_4j
   # calculated with HEJ revision 13207b5f67a5f40a2141aa7ee515b022bd4efb65
   COMMAND xs_gen ${runcard_dir}/4j.yml 915072 15402.4
 )
 ## W
 add_test(
   NAME xs_WqQ
   # calculated with Sherpa see #132
   COMMAND xs_gen ${runcard_dir}/WpqQ.yml 3.086e+00 4.511e-02
 )
 add_test(
   NAME xs_W2j
   # calculated with HEJ revision 449f2f6b597020e9c9e35699568edc05c827fc11+1
   COMMAND xs_gen ${runcard_dir}/Wm2j.yml 4.177443e+02 7.446928e+00
 )
 add_long_test(
   NAME xs_W3j_eqqx
   # calculated with HEJ revision 449f2f6b597020e9c9e35699568edc05c827fc11+1
   COMMAND xs_gen ${runcard_dir}/Wp3j_qqx.yml 8.710751e+00 1.245725e-01
 )
 add_long_test(
   NAME xs_W3j_uno
   # calculated with HEJ revision 449f2f6b597020e9c9e35699568edc05c827fc11+1
   COMMAND xs_gen ${runcard_dir}/Wp3j_uno.yml 3.000955e-01 5.831799e-03
 )
 add_long_test(
   NAME xs_W4j_qqx
   # calculated with HEJ revision be065dc9a21e5965ce57583f6c0a3d953664b82b
   COMMAND xs_gen ${runcard_dir}/Wp4j_qqx.yml 9.274718e-02 4.875742e-03
 )
 
 # Test that sum of partons == proton
 add_executable(PSP_channel PSP_channel.cc)
 target_link_libraries(PSP_channel hejfog_lib)
 # pure jets
 add_test(
   NAME channel_2j
   COMMAND PSP_channel ${runcard_dir}/2j.yml
 )
 add_test(
   NAME channel_3j_qqx
   COMMAND PSP_channel ${runcard_dir}/3j_qqx.yml
 )
 add_test(
   NAME channel_3j_uno
   COMMAND PSP_channel ${runcard_dir}/3j_uno.yml
 )
 add_long_test(
   NAME channel_4j_qqx
   COMMAND PSP_channel ${runcard_dir}/4j_qqx.yml
 )
 # W+jets
 # pure jets
 add_test(
   NAME channel_W2j
   COMMAND PSP_channel ${runcard_dir}/Wm2j.yml
 )
 add_test(
   NAME channel_W3j_uno
   COMMAND PSP_channel ${runcard_dir}/Wp3j_uno.yml
 )
 add_test(
   NAME channel_W3j_eqqx
   COMMAND PSP_channel ${runcard_dir}/Wp3j_qqx.yml
 )
 add_long_test(
   NAME channel_W4j_qqx
   COMMAND PSP_channel ${runcard_dir}/Wp4j_qqx.yml
 )
 
 # Test that each subleading channel is generated consistently
 add_executable(PSP_subl PSP_subl.cc)
 target_link_libraries(PSP_subl hejfog_lib)
 add_long_test(
   NAME subl_5j
   COMMAND PSP_subl ${runcard_dir}/5j.yml
 )
 add_long_test(
   NAME subl_h5j
   COMMAND PSP_subl ${runcard_dir}/h5j.yml
 )
 add_long_test(
   NAME subl_W5j
   COMMAND PSP_subl ${runcard_dir}/Wm5j.yml
 )
diff --git a/t/CMakeLists.txt b/t/CMakeLists.txt
index af9db00..f0375ef 100644
--- a/t/CMakeLists.txt
+++ b/t/CMakeLists.txt
@@ -1,493 +1,497 @@
 set(tst_dir "${CMAKE_CURRENT_SOURCE_DIR}")
 set(tst_ME_data_dir "${tst_dir}/ME_data")
 
 # small library for common test functions
 add_library(hej_test SHARED hej_test.cc)
 target_include_directories(hej_test PUBLIC ${tst_dir})
 target_link_libraries(hej_test HEJ)
 # test event classification
 # test explicit configurations
 add_executable(test_classify ${tst_dir}/test_classify.cc)
 target_compile_options(test_classify PRIVATE "-O0") # avoid compiler optimisation
 target_link_libraries(test_classify HEJ hej_test)
 add_test(
   NAME classify
   COMMAND test_classify
   )
 # test against reference data
 add_executable(test_classify_ref ${tst_dir}/test_classify_ref.cc)
 target_link_libraries(test_classify_ref HEJ hej_test)
 add_test(
   NAME classify_ref
   COMMAND test_classify_ref ${tst_dir}/classify_ref ${tst_dir}/classify.lhe.gz
   )
 add_test(
   NAME classify_ref_4j
   COMMAND test_classify_ref ${tst_dir}/classify_ref_4j ${tst_dir}/4j.lhe.gz
   )
 add_test(
   NAME classify_ref_W4j
   COMMAND test_classify_ref
           ${tst_dir}/classify_ref_W4j
           ${tst_dir}/MGemnubar4j_unweighted.lhe.gz
   )
 add_test(
   NAME t_classify_ref_Z4j
   COMMAND test_classify_ref
           ${tst_dir}/classify_ref_Z4j
           ${tst_dir}/Z_4j.lhe.gz
   )
 
 # test for valid W decays
 add_executable(test_decay ${tst_dir}/test_decay.cc)
 target_link_libraries(test_decay HEJ hej_test)
 add_test(
   NAME valid_decay
   COMMAND test_decay
   )
 
 # test valid jet cuts on tagging jets
 add_executable(test_jet_cuts ${tst_dir}/test_jet_cuts.cc)
 target_link_libraries(test_jet_cuts HEJ hej_test)
 add_test(
   NAME jet_cuts
   COMMAND test_jet_cuts
   )
 
 # test phase space point
 add_executable(test_psp ${tst_dir}/test_psp.cc)
 target_link_libraries(test_psp HEJ hej_test)
 add_test(
   NAME PhaseSpace
   COMMAND test_psp ${tst_dir}/psp_gen.lhe.gz
   )
 
 # test importing analyses
 file(COPY "${tst_dir}/analysis_config.yml" DESTINATION "${CMAKE_CURRENT_BINARY_DIR}")
 set(test_config "${CMAKE_CURRENT_BINARY_DIR}/analysis_config.yml")
 
 get_target_property(ANALYSIS_PATH AnalysisTemplate_lib BINARY_DIR)
 get_target_property(ANALYSIS_LIB AnalysisTemplate_lib OUTPUT_NAME)
 set(ANALYSIS_PARAMETERS "")
 file(READ ${test_config} config)
 file(WRITE ${test_config} "${config}  - plugin: ${ANALYSIS_PATH}/lib${ANALYSIS_LIB}.so\n  ${ANALYSIS_PARAMETERS}\n")
 configure_file( ${tst_dir}/analysis_config_single.yml.in
                 analysis_config_simple.yml  @ONLY )
 add_test(
   NAME analysis_simple
   COMMAND $<TARGET_FILE:HEJ_main>
           analysis_config_simple.yml
           ${tst_dir}/2j.lhe.gz
 )
 
 get_target_property(ANALYSIS_PATH AnalysisPrint_lib BINARY_DIR)
 get_target_property(ANALYSIS_LIB AnalysisPrint_lib OUTPUT_NAME)
 set(ANALYSIS_PARAMETERS "  output: ana_output")
 file(READ ${test_config} config)
 file(WRITE ${test_config} "${config}  - plugin: ${ANALYSIS_PATH}/lib${ANALYSIS_LIB}.so\n  ${ANALYSIS_PARAMETERS}\n")
 configure_file( ${tst_dir}/analysis_config_single.yml.in
                 analysis_config_print.yml  @ONLY )
 add_test(
   NAME analysis_print
   COMMAND $<TARGET_FILE:HEJ_main>
           analysis_config_print.yml
           ${tst_dir}/2j.lhe.gz
 )
 
 if(ROOT_FOUND)
   get_target_property(ANALYSIS_PATH AnalysisROOT_lib BINARY_DIR)
   get_target_property(ANALYSIS_LIB AnalysisROOT_lib OUTPUT_NAME)
   set(ANALYSIS_PARAMETERS "")
   file(READ ${test_config} config)
   file(WRITE ${test_config} "${config}  - plugin: ${ANALYSIS_PATH}/lib${ANALYSIS_LIB}.so\n  ${ANALYSIS_PARAMETERS}\n")
   configure_file( ${tst_dir}/analysis_config_single.yml.in
                   analysis_config_root.yml  @ONLY )
   add_test(
     NAME analysis_root
     COMMAND $<TARGET_FILE:HEJ_main>
             analysis_config_root.yml
             ${tst_dir}/2j.lhe.gz
   )
 endif()
 if(RIVET_FOUND)
   file(READ ${test_config} config)
   file(WRITE ${test_config} "${config}  - rivet: MC_XS\n    output: ana_rivet\n")
 endif()
 
 add_test(
   NAME analysis_all
   COMMAND $<TARGET_FILE:HEJ_main>
           ${test_config}
           ${tst_dir}/2j.lhe.gz
 )
 
 # test importing scales (from examples/softestptScale)
 add_executable(test_scale_import ${tst_dir}/test_scale_import)
 target_link_libraries(test_scale_import HEJ)
 
 get_target_property(SCALE_PATH softestptScale_lib BINARY_DIR)
 get_target_property(SCALE_LIB softestptScale_lib OUTPUT_NAME)
 set(SCALE_NAME "softest_jet_pt")
 configure_file( ${tst_dir}/jet_config_with_import.yml.in
                 jet_config_with_import.yml  @ONLY )
 
 add_test(
   NAME scale_import
   COMMAND test_scale_import jet_config_with_import.yml
 )
 
 # test scale arithmetic (e.g. 2*H_T/4)
 add_executable(test_scale_arithmetics ${tst_dir}/test_scale_arithmetics)
 target_link_libraries(test_scale_arithmetics HEJ hej_test)
 add_test(
   NAME scale_arithmetics
   COMMAND test_scale_arithmetics ${tst_dir}/jet_config.yml ${tst_dir}/2j.lhe.gz
   )
 
 # test "ParameterDescription"
 add_executable(test_descriptions ${tst_dir}/test_descriptions)
 target_link_libraries(test_descriptions HEJ hej_test)
 add_test(
   NAME descriptions
   COMMAND test_descriptions
   )
 
 # test "EventParameters*Weight"
 add_executable(test_parameters ${tst_dir}/test_parameters)
 target_link_libraries(test_parameters HEJ hej_test)
 add_test(
   NAME parameters
   COMMAND test_parameters
   )
 
 # test Particle IDs
 add_executable(test_pid ${tst_dir}/test_pid)
 target_link_libraries(test_pid HEJ hej_test)
 add_test(
   NAME pid
   COMMAND test_pid
   )
 
 # test unweighting
 add_executable(test_unweighter ${tst_dir}/test_unweighter)
 target_link_libraries(test_unweighter HEJ hej_test)
 add_test(
   NAME unweighter
   COMMAND test_unweighter ${tst_dir}/4j.lhe.gz
   )
 
 # test colour generation
 add_executable(test_colours ${tst_dir}/test_colours)
 target_link_libraries(test_colours HEJ hej_test)
 add_test(
   NAME colour_flow
   COMMAND test_colours
   )
 add_executable(test_colours2 ${tst_dir}/test_colours2)
 target_link_libraries(test_colours2 HEJ hej_test)
 add_test(
   NAME colour_flow2
   COMMAND test_colours2
   )
 
 # test matrix elements
 add_executable(test_ME_generic ${tst_dir}/test_ME_generic.cc)
 target_link_libraries(test_ME_generic HEJ hej_test)
 add_test(
   NAME ME_j
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_mtinf.yml
           ${tst_ME_data_dir}/ME_jets_tree.dat ${tst_ME_data_dir}/PSP_jets.lhe.gz
   )
 add_test(
   NAME ME_j_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_mtinf.yml
           ${tst_ME_data_dir}/ME_jets_tree_virt.dat ${tst_ME_data_dir}/PSP_jets.lhe.gz
   )
 add_test(
   NAME ME_h
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_mtinf.yml
           ${tst_ME_data_dir}/ME_h_mtinf_tree.dat ${tst_ME_data_dir}/PSP_h.lhe.gz
   )
 add_test(
   NAME ME_h_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_mtinf.yml
           ${tst_ME_data_dir}/ME_h_mtinf_virt.dat ${tst_ME_data_dir}/PSP_h.lhe.gz
   )
 if(QCDloop_FOUND)
   add_test(
     NAME ME_h_mt
     COMMAND test_ME_generic ${tst_ME_data_dir}/config_mt.yml
             ${tst_ME_data_dir}/ME_h_mt_tree.dat ${tst_ME_data_dir}/PSP_h.lhe.gz
     )
   add_test(
     NAME ME_h_mtmb
     COMMAND test_ME_generic ${tst_ME_data_dir}/config_mtmb.yml
             ${tst_ME_data_dir}/ME_h_mtmb_tree.dat ${tst_ME_data_dir}/PSP_h.lhe.gz
     )
 endif()
 add_test(
   NAME ME_j_subl
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_pure.yml
           ${tst_ME_data_dir}/ME_jets_subl_tree.dat ${tst_ME_data_dir}/PSP_jets_subl.lhe.gz
   )
 add_test(
   NAME ME_j_subl_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_pure.yml
           ${tst_ME_data_dir}/ME_jets_subl_tree_virt.dat ${tst_ME_data_dir}/PSP_jets_subl.lhe.gz
   )
 add_test(
   NAME ME_j_subl_4j
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_pure.yml
           ${tst_ME_data_dir}/ME_jets_subl_tree_new.dat ${tst_dir}/4j.lhe.gz
   )
 add_test(
   NAME ME_j_subl_4j_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_pure.yml
           ${tst_ME_data_dir}/ME_jets_subl_tree_new_virt.dat ${tst_dir}/4j.lhe.gz
   )
 add_test(
   NAME ME_w_FKL
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_w_ME.yml
           ${tst_ME_data_dir}/ME_w_FKL_tree.dat ${tst_ME_data_dir}/PSP_w_FKL.lhe.gz
   )
 add_test(
   NAME ME_w_FKL_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_w_ME.yml
           ${tst_ME_data_dir}/ME_w_FKL_virt.dat ${tst_ME_data_dir}/PSP_w_FKL.lhe.gz
   )
 add_test(
   NAME ME_Wp
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_w_ME.yml
           ${tst_ME_data_dir}/ME_Wp.dat ${tst_ME_data_dir}/PSP_Wp.lhe.gz
   )
 add_test(
   NAME ME_Wp_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_w_ME.yml
           ${tst_ME_data_dir}/ME_Wp_virt.dat ${tst_ME_data_dir}/PSP_Wp.lhe.gz
   )
 add_test(
   NAME ME_Wm
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_w_ME.yml
           ${tst_ME_data_dir}/ME_Wm.dat ${tst_ME_data_dir}/PSP_Wm.lhe.gz
   )
 add_test(
   NAME ME_Wm_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_w_ME.yml
           ${tst_ME_data_dir}/ME_Wm_virt.dat ${tst_ME_data_dir}/PSP_Wm.lhe.gz
   )
 add_test(
   NAME t_ME_Z_FKL
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_Z_ME.yml
           ${tst_ME_data_dir}/ME_Z_FKL_tree.dat ${tst_ME_data_dir}/PSP_Z_FKL.lhe.gz
   )
 add_test(
   NAME t_ME_Z_FKL_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_Z_ME.yml
           ${tst_ME_data_dir}/ME_Z_FKL_virt.dat ${tst_ME_data_dir}/PSP_Z_FKL.lhe.gz
   )
 add_test(
   NAME t_ME_Z_uno
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_Z_ME.yml
           ${tst_ME_data_dir}/ME_Z_uno_tree.dat ${tst_ME_data_dir}/PSP_Z_uno.lhe.gz
   )
 add_test(
   NAME t_ME_Z_uno_virt
   COMMAND test_ME_generic ${tst_ME_data_dir}/config_Z_ME.yml
           ${tst_ME_data_dir}/ME_Z_uno_virt.dat ${tst_ME_data_dir}/PSP_Z_uno.lhe.gz
   )
 
 # test main executable
 file(COPY "${tst_dir}/jet_config.yml" DESTINATION "${CMAKE_CURRENT_BINARY_DIR}")
 set(test_config "${CMAKE_CURRENT_BINARY_DIR}/jet_config.yml")
 if(HighFive_FOUND)
   file(READ ${test_config} config)
   file(WRITE ${test_config} "${config}  - tst.hdf5\n")
 endif()
 if(HepMC3_FOUND)
   file(READ ${test_config} config)
   file(WRITE ${test_config} "${config}  - tst.hepmc\n")
 endif()
 if(HepMC_FOUND)
   file(READ ${test_config} config)
   file(WRITE ${test_config} "${config}  - tst.hepmc2\n")
 endif()
 if(rivet_FOUND)
   file(READ ${test_config} config)
   file(WRITE ${test_config} "${config}\nanalysis:\n  rivet: MC_XS\n  output: tst\n")
 endif()
 set(test_cmd_main "$<TARGET_FILE:HEJ_main>\\\;${test_config}\\\;${tst_dir}/2j.lhe.gz")
 # check that HepMC3 output is correct
 if(HepMC3_FOUND)
   add_executable(check_hepmc ${tst_dir}/check_hepmc.cc)
   target_link_libraries(check_hepmc HEJ ${HEPMC3_LIBRARIES})
   target_include_directories(check_hepmc PRIVATE ${HEPMC3_INCLUDE_DIR})
   set(test_cmd_hepmc "$<TARGET_FILE:check_hepmc>\\\;tst.hepmc")
 else()
   set(test_cmd_hepmc "")
 endif()
 # check that LHEF output is correct
 add_executable(check_lhe ${tst_dir}/check_lhe.cc)
 target_link_libraries(check_lhe HEJ hej_test)
 set(test_cmd_lhe "$<TARGET_FILE:check_lhe>\\\;tst.lhe")
 # check that rivet interface is consistent with naive rivet
 if(rivet_FOUND)
   # this assumes "rivet" and "yodadiff" are found in PATH
   if(rivet_USE_HEPMC3)
     set(hepmc_file "tst.hepmc")
   else()
     set(hepmc_file "tst.hepmc2")
   endif()
   if(rivet_USE_HEPMC3 OR (rivet_VERSION VERSION_LESS 3))
     set(histo_exclude "")
   else()
     # rivet 3 with HepMC 2 is inconsistent in order of weights
     # -> interface != direct call (by permutation)
     # REQUIRES Yoda 1.7.5
     set(histo_exclude "-M\\\;\\\\d")
   endif()
   set(test_cmd_rivet "rivet\\\;-a\\\;MC_XS\\\;${hepmc_file}\\\;-o\\\;tst_direct.yoda\
 \;yodadiff\\\;${histo_exclude}\\\;tst.yoda\\\;tst_direct.yoda")
 else()
   set(test_cmd_rivet "")
 endif()
 # Run dependent tests in one command to ensure correct execution order
 # Note: The commands are concatenated with "\;" to escape CMake lists.
 #       Thus arguments have to be escaped twice "\\\;".
 #       e.g. "cmd1\;cmd2\\\;arg1\\\;arg2" is executed like "cmd1 && cmd2 arg1 arg2"
 add_test(
   NAME main
   COMMAND ${CMAKE_COMMAND}
     -DCMDS=${test_cmd_main}\;${test_cmd_hepmc}\;${test_cmd_lhe}\;${test_cmd_rivet}
     -P ${PROJECT_SOURCE_DIR}/cmake/run_multiple_tests.cmake
   )
+add_test(
+  NAME main_example
+  COMMAND $<TARGET_FILE:HEJ_main> ${PROJECT_SOURCE_DIR}/config.yml ${tst_dir}/h_3j.lhe.gz
+  )
 # check that Sherpas LHE input can be read
 add_executable(check_lhe_sherpa ${tst_dir}/check_lhe_sherpa.cc)
 target_link_libraries(check_lhe_sherpa HEJ hej_test)
 add_test(
   NAME sherpa_reader
   COMMAND check_lhe_sherpa ${tst_dir}/SherpaLHE.lhe 1.62624e+08
   )
 # check HDF5 reader & writer
 if(HighFive_FOUND)
   add_executable(test_hdf5 ${tst_dir}/test_hdf5.cc)
   target_link_libraries(test_hdf5 HEJ)
   add_test(
     NAME hdf5_read
     COMMAND test_hdf5 ${tst_dir}/Wm9-g4-repack.hdf5
     )
   add_executable(test_hdf5_write ${tst_dir}/test_hdf5_write.cc)
   target_link_libraries(test_hdf5_write HEJ hej_test)
   add_test(
     NAME hdf5_write
     COMMAND test_hdf5_write ${tst_dir}/Wm9-g4-repack.hdf5 output.hdf5
     )
 endif()
 # check rivet interface
 if(RIVET_FOUND)
   add_executable(check_rivet ${tst_dir}/check_rivet.cc)
   target_link_libraries(check_rivet HEJ rivet::rivet)
   add_test(
     NAME rivet
     COMMAND check_rivet
     )
 endif()
 
 # test boson reconstruction
 add_executable(cmp_events ${tst_dir}/cmp_events.cc)
 target_link_libraries(cmp_events HEJ)
 add_test(
   NAME reconstruct_W
   COMMAND cmp_events ${tst_dir}/epnu2jLOFKL_unweight.lhe.tar.gz
           ${tst_dir}/epnu2jLOFKL_unweight_noW.lhe.gz
   )
 
 # test resummed result
 add_executable(check_res ${tst_dir}/check_res.cc)
 target_link_libraries(check_res HEJ hej_test)
 
 if(TEST_ALL) # deactivate long tests by default
   add_test(
     NAME xs_2j
     COMMAND check_res ${tst_dir}/2j.lhe.gz 3.49391e+07 419684
     )
   add_test(
     NAME xs_3j
     COMMAND check_res ${tst_dir}/3j.lhe.gz  2.37902e+06 25746.6
     )
   add_test(
     NAME xs_3j_unof
     COMMAND check_res ${tst_dir}/3j.lhe.gz  133399 4688.83 unof
     )
   add_test(
     NAME xs_3j_unob
     COMMAND check_res ${tst_dir}/3j.lhe.gz  105247 3449.45 unob
     )
   add_test(
     NAME xs_3j_splitf
     COMMAND check_res ${tst_dir}/3j.lhe.gz  97659.9 2748.65 splitf
     )
   add_test(
     NAME xs_3j_splitb
     COMMAND check_res ${tst_dir}/3j.lhe.gz  107150 2799.8 splitb
     )
   add_test(
     NAME xs_4j
     COMMAND check_res ${tst_dir}/4j.lhe.gz  603713 72822.6
     )
   add_test(
     NAME xs_4j_qqxmid
     COMMAND check_res ${tst_dir}/4j.lhe.gz  21866.7 1716.96 qqxmid
     )
   add_test(
     NAME xs_h_3j
     COMMAND check_res ${tst_dir}/h_3j.lhe.gz 0.821622 0.0220334
     )
   add_test(
     NAME xs_h_3j_unof
     COMMAND check_res ${tst_dir}/h_3j_uno.lhe.gz 0.0127362 0.000271555 unof
     )
   add_test(
     NAME xs_h_3j_unob
     COMMAND check_res ${tst_dir}/h_3j_uno.lhe.gz 0.0130615 0.000224793 unob
     )
   add_test(
     NAME xs_epnu_2j
     COMMAND check_res ${tst_dir}/epnu2jLOFKL_unweight.lhe.tar.gz 262.7 3
     )
   add_test(
     NAME xs_epnu_3j
     COMMAND check_res ${tst_dir}/MGepnu3j_unweighted.lhe.gz 38.9512 1
     )
   add_test(
     NAME xs_emnubar_3j
     COMMAND check_res ${tst_dir}/MGemnubar3j_unweighted.lhe.gz 24.1575 1
     )
   add_test(
     NAME xs_epnu_3j_unof
     COMMAND check_res ${tst_dir}/MGepnu3j_unweighted.lhe.gz 9.63702 0.128355 unof
     )
   add_test(
     NAME xs_epnu_3j_unob
     COMMAND check_res ${tst_dir}/MGepnu3j_unweighted.lhe.gz 9.70119 0.108436 unob
     )
   add_test(
     NAME xs_epnu_3j_splitf
     COMMAND check_res ${tst_dir}/MGepnu3j_unweighted.lhe.gz 2.91995 0.0463182 splitf
     )
   add_test(
     NAME xs_epnu_3j_splitb
     COMMAND check_res ${tst_dir}/MGepnu3j_unweighted.lhe.gz 3.40708 0.0550975 splitb
     )
   add_test(
     NAME xs_epnu_4j
     COMMAND check_res ${tst_dir}/MGepnu4j_unweighted.lhe.gz 10.2542 0.135106
     )
   add_test(
     NAME xs_emnubar_4j
     COMMAND check_res ${tst_dir}/MGemnubar4j_unweighted.lhe.gz 5.73282 0.0583738
     )
   add_test(
     NAME xs_epnu_4j_qqxmid
     COMMAND check_res ${tst_dir}/MGepnu4j_unweighted.lhe.gz 0.802603 0.0126908 qqxmid
     )
 endif()