No OneTemporary
Actions

Size

34 KB

Subscribers

None

View Options

	diff --git a/CMakeLists.txt b/CMakeLists.txt
	index 81df100..6df4bc5 100644
	--- a/CMakeLists.txt
	+++ b/CMakeLists.txt
	@@ -1,134 +1,140 @@
	cmake_minimum_required(VERSION 2.8 FATAL_ERROR)

	set(CMAKE_LEGACY_CYGWIN_WIN32 0)
	set(CMAKE_EXPORT_COMPILE_COMMANDS ON)

	project("Reversed HEJ" C CXX)

	## Flags for the compiler. No warning allowed.
	if ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "GNU")
	set(warnings "-Wall -Wextra -Werror -std=c++1y")
	elseif ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "Clang")
	set(warnings "-Wall -Wextra -Werror -std=c++11")
	elseif ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "MSVC")
	set(warnings "/W4 /WX /EHsc")
	endif()
	if (NOT CONFIGURED_ONCE)
	set(CMAKE_CXX_FLAGS "${warnings}"
	CACHE STRING "Flags used by the compiler during all build types." FORCE)
	set(CMAKE_C_FLAGS "${warnings}"
	CACHE STRING "Flags used by the compiler during all build types." FORCE)
	endif()


	## Add directories and find dependences
	include_directories(${CMAKE_CURRENT_SOURCE_DIR}/include)
	set(CMAKE_MODULE_PATH ${CMAKE_MODULE_PATH} "${CMAKE_SOURCE_DIR}/cmake/Modules/")

	find_package(fastjet REQUIRED)
	find_package(clhep REQUIRED)
	find_package(lhapdf REQUIRED)
	find_package(gsl REQUIRED)
	find_package(Boost REQUIRED COMPONENTS iostreams)
	find_package(HepMC 3)
	if(${HepMC_FOUND})
	set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -DRHEJ_BUILD_WITH_HepMC")
	endif()

	include_directories(SYSTEM ${lhapdf_INCLUDE_PATH})
	include_directories(SYSTEM ${fastjet_INCLUDE_PATH})
	include_directories(SYSTEM ${clhep_INCLUDE_PATH})
	include_directories(SYSTEM ${gsl_INCLUDE_PATH})

	add_subdirectory(src)

	## define executable
	add_executable(rHEJ src/main.cc)

	## link libraries
	target_link_libraries(rHEJ rhej)

	#plugins
	add_subdirectory("analysis-plugins")

	file(GLOB rhej_headers ${CMAKE_CURRENT_SOURCE_DIR}/include/RHEJ/*.hh)
	install(FILES ${rhej_headers} DESTINATION include/RHEJ/)
	install(TARGETS rHEJ DESTINATION bin)

	enable_testing()
	set(tst_dir "${CMAKE_CURRENT_SOURCE_DIR}/t")
	add_executable(test_Matrix ${tst_dir}/test_Matrix.cc)
	add_executable(test_classify ${tst_dir}/test_classify.cc)
	add_executable(test_psp ${tst_dir}/test_psp.cc)
	add_executable(test_ME_h_3j ${tst_dir}/test_ME_h_3j.cc)
	+add_executable(test_ME_hjets_mt174 ${tst_dir}/test_ME_hjets_mt174.cc)
	add_executable(check_res ${tst_dir}/check_res.cc)
	add_executable(check_lhe ${tst_dir}/check_lhe.cc)
	target_link_libraries(test_Matrix rhej)
	target_link_libraries(test_classify rhej)
	target_link_libraries(test_psp rhej)
	target_link_libraries(test_ME_h_3j rhej)
	+target_link_libraries(test_ME_hjets_mt174 rhej)
	target_link_libraries(check_res rhej)
	target_link_libraries(check_lhe rhej)

	add_test(
	NAME t_matrix
	COMMAND test_Matrix
	)
	add_test(
	NAME t_classify
	COMMAND test_classify ${tst_dir}/classify.lhe.gz
	)
	add_test(
	NAME t_psp
	COMMAND test_psp ${tst_dir}/psp_gen.lhe.gz
	)
	add_test(
	NAME t_ME
	COMMAND test_ME_h_3j
	)
	add_test(
	+ NAME t_ME_mt174
	+ COMMAND test_ME_hjets_mt174
	+ )
	+add_test(
	NAME t_2j
	COMMAND check_res ${tst_dir}/2j.lhe.gz 3.382e7 752159
	)
	add_test(
	NAME t_3j
	COMMAND check_res ${tst_dir}/3j.lhe.gz 2.418e6 67247.8
	)
	add_test(
	NAME t_4j
	COMMAND check_res ${tst_dir}/4j.lhe.gz 551096 30572.5
	)
	add_test(
	NAME t_h_3j
	COMMAND check_res ${tst_dir}/h_3j.lhe.gz 0.793107 0.0479054
	)
	add_test(
	NAME t_h_3j_uno
	COMMAND check_res ${tst_dir}/h_3j_uno.lhe.gz 0.0246366 0.001 uno
	)
	if(${HepMC_FOUND})
	file(READ "${tst_dir}/jet_config.yml" config)
	file(WRITE "${tst_dir}/jet_config_withHepMC.yml" "${config} - tst.hepmc3")
	add_executable(check_hepmc ${tst_dir}/check_hepmc.cc)
	target_link_libraries(check_hepmc rhej)
	add_test(
	NAME t_main
	COMMAND rHEJ ${tst_dir}/jet_config_withHepMC.yml ${tst_dir}/2j.lhe.gz
	)
	add_test(
	NAME t_hepmc
	COMMAND check_hepmc tst.hepmc3
	)
	else()
	add_test(
	NAME t_main
	COMMAND rHEJ ${tst_dir}/jet_config.yml ${tst_dir}/2j.lhe.gz
	)
	endif()
	add_test(
	NAME t_lhe
	COMMAND check_lhe tst.lhe
	)

	set(CONFIGURED_ONCE TRUE CACHE INTERNAL
	"A flag showing that CMake has configured at least once.")
	diff --git a/t/ME_hjets_mt174.dat b/t/ME_hjets_mt174.dat
	new file mode 100644
	index 0000000..821fca4
	--- /dev/null
	+++ b/t/ME_hjets_mt174.dat
	@@ -0,0 +1,492 @@
	+Event{
	+ 2.41858843761886e-16,
	+ {PDG(21), {0,0,-680.117145506685,680.117145506685}},
	+ {PDG(21), {0,0,3659.59510295182,3659.59510295182}},
	+ {
	+ {PDG(21), {-279.955685772138,667.735101894511,-51.1713137794792,725.853880372723}},
	+ {PDG(21), {27.5142656992431,-14.4104401037,0.942189551118422,31.0738366170177}},
	+ {PDG(21), {10.7340626328525,-14.7795338976409,42.7206058838912,46.4618648993374}},
	+ {PDG(21), {12.0512256423466,-12.4464813880144,41.2141700349669,44.7074350650549}},
	+ {PDG(21), {-1.77612744061724,4.90817213613428,32.9089496311042,33.3203203499875}},
	+ {PDG(21), {-0.863102454281029,3.33521150291264,29.9172043145537,30.1149081953523}},
	+ {PDG(21), {-0.396353150549452,1.59373177696446,28.503203619967,28.5504762376769}},
	+ {PDG(21), {-16.18472886774,-26.2763398759877,580.223321038203,581.043452559619}},
	+ {PDG(21), {-1.46164323748031,6.30413835855595,129.284251838758,129.446113633839}},
	+ {PDG(21), {47.424158708006,7.9664387825176,2006.462374337,2007.0385583257}},
	+ {PDG(25), {202.913928240358,-623.929999186253,138.473000975052,682.101402202195}}
	+ }
	+},
	+Event{
	+ 1.70239154781339e-13,
	+ {PDG(2), {0,0,-101.039878086234,101.039878086234}},
	+ {PDG(21), {0,0,1403.98431968767,1403.98431968767}},
	+ {
	+ {PDG(2), {-62.7762937466985,-18.9590916883968,-11.4064793094033,66.5613850851727}},
	+ {PDG(21), {-8.53082502152277,29.1311519730581,-2.37778588750478,30.4475427013751}},
	+ {PDG(21), {0.863585608169176,1.21045173247596,0.3574540024056,1.52929619862431}},
	+ {PDG(21), {17.522259001433,4.19651954710253,37.564621609456,41.662226698583}},
	+ {PDG(21), {-1.22291064453758,2.58971166953696,6.11708072464106,6.75431666178263}},
	+ {PDG(21), {2.74348966983159,1.35891821530508,8.38636257662644,8.92773608206476}},
	+ {PDG(21), {-0.543324027882566,4.1965913160823,18.3489976711284,18.8306212114356}},
	+ {PDG(21), {-0.305992521695085,1.41867374611182,7.62092062349878,7.75787972134844}},
	+ {PDG(21), {-1.32053088020405,0.0425194257961066,8.11680904767516,8.22363659359798}},
	+ {PDG(21), {39.1268647288542,-14.0265730315661,276.622654219559,279.727991309403}},
	+ {PDG(21), {-51.3819022287447,74.4162458671198,812.07659491703,817.096244966061}},
	+ {PDG(25), {65.8255800629972,-85.5751187726256,141.517211406325,217.505320544457}}
	+ }
	+},
	+Event{
	+ 5.34901661444232e-12,
	+ {PDG(21), {0,0,-405.402177074155,405.402177074155}},
	+ {PDG(21), {0,0,38.0986580811438,38.0986580811438}},
	+ {
	+ {PDG(21), {-27.1802918847255,-60.8133439796394,-119.707616143782,136.992497735285}},
	+ {PDG(21), {35.7325549408384,-7.63574700373474,-36.3364648702955,51.5311439245753}},
	+ {PDG(25), {-8.55226305611281,68.4490909833741,-211.259437978933,254.977193495439}}
	+ }
	+},
	+Event{
	+ 1.33550696279992e-16,
	+ {PDG(21), {0,0,-87.0631830426542,87.0631830426542}},
	+ {PDG(2), {0,0,854.329839061515,854.329839061515}},
	+ {
	+ {PDG(21), {-56.1870623050712,-10.9768805506147,-11.2280135542579,58.3399191417936}},
	+ {PDG(21), {26.701957312687,-9.22570603044948,20.1179330302083,34.6819752248452}},
	+ {PDG(21), {3.84322135690833,11.5116152689347,10.6966152220222,16.1773055143977}},
	+ {PDG(21), {4.58729617315599,-2.1744250951473,4.8183229428243,6.99911756266496}},
	+ {PDG(21), {1.35079846694462,-2.53949074197776,5.15002700988233,5.89885140763931}},
	+ {PDG(2), {-138.882682639874,54.0800596019451,274.550722594371,312.395825293631}},
	+ {PDG(25), {158.586471635249,-40.6751724526906,463.161048773811,506.900027959197}}
	+ }
	+},
	+Event{
	+ 2.00092597157308e-10,
	+ {PDG(21), {0,0,-617.807743703745,617.807743703745}},
	+ {PDG(21), {0,0,1394.0427641665,1394.0427641665}},
	+ {
	+ {PDG(21), {8.95282214356842,-95.1553510617665,540.577773072874,548.961767885934}},
	+ {PDG(21), {-36.9270394494355,5.10424180853012,238.144291096428,241.04431731273}},
	+ {PDG(21), {-35.5918634797112,-33.9337278674189,594.330347955196,596.361334371606}},
	+ {PDG(25), {63.5660807855783,123.984837120655,-596.81739166174,625.483088299977}}
	+ }
	+},
	+Event{
	+ 6.10485794239085e-14,
	+ {PDG(21), {0,0,-210.479775924028,210.479775924028}},
	+ {PDG(1), {0,0,3347.52410659223,3347.52410659223}},
	+ {
	+ {PDG(21), {-29.345169132615,-45.7983095873898,-113.441918354287,125.80815932357}},
	+ {PDG(21), {-11.8767796022229,-16.9997140738415,-34.0704615149647,39.8853923154183}},
	+ {PDG(21), {3.25082473702128,-7.43553770383388,-11.9535524091472,14.4479236437003}},
	+ {PDG(21), {-0.821259007904344,1.97651688216404,-0.110391870427621,2.14319194394716}},
	+ {PDG(21), {-24.7270186011526,0.598731114982081,3.24786002763389,24.9465934068942}},
	+ {PDG(21), {-1.76825657201983,-0.393173714152914,2.04530355399262,2.73213899755612}},
	+ {PDG(21), {-3.23817575061384,-3.87772915771791,6.06637241507212,7.89451961114969}},
	+ {PDG(21), {-13.6044985077785,25.2407046948561,114.01662066349,117.56685477065}},
	+ {PDG(21), {-6.49553455764411,-4.47435843874359,36.7643883276061,37.600958787389}},
	+ {PDG(21), {0.833550157095199,-2.08222782502452,13.1871371661159,13.3765116984083}},
	+ {PDG(1), {-32.8290045969523,-196.488789098994,1995.94454324984,2005.86146269208}},
	+ {PDG(25), {120.621321434787,249.733886907696,1125.34842941327,1165.74017532549}}
	+ }
	+},
	+Event{
	+ 1.15327935186291e-13,
	+ {PDG(4), {0,0,-340.664687276466,340.664687276466}},
	+ {PDG(2), {0,0,2038.60644774981,2038.60644774981}},
	+ {
	+ {PDG(4), {93.227687614025,-0.237385800108499,-280.862849721683,295.931408342644}},
	+ {PDG(21), {-52.4413575317427,-34.7494646378399,10.5765956997302,63.7925203216185}},
	+ {PDG(21), {-8.36540885973258,16.4405420131842,9.85530995529677,20.9140771107103}},
	+ {PDG(21), {19.9408369042266,-8.22668516423455,42.2028366901075,47.3961469945004}},
	+ {PDG(21), {0.23737419735702,-0.594573967142463,2.48010555658275,2.56140357690214}},
	+ {PDG(21), {-1.60930349133105,0.417687612200225,10.8539945867435,10.9805973952993}},
	+ {PDG(21), {-24.3804397377401,-3.95406400311759,193.620411008942,195.189405509669}},
	+ {PDG(21), {-1.8492854063318,2.33009137487207,24.7111672925605,24.889575554642}},
	+ {PDG(2), {-117.320411538935,8.14866075962713,1252.5334385282,1258.04232610318}},
	+ {PDG(25), {92.5603078502044,20.4251918125594,431.970750876869,459.573674117109}}
	+ }
	+},
	+Event{
	+ 1.23013302494519e-13,
	+ {PDG(2), {0,0,-3052.69985945628,3052.69985945628}},
	+ {PDG(21), {0,0,948.211487135056,948.211487135056}},
	+ {
	+ {PDG(2), {-21.8643801848913,28.6777167771716,-52.5141789276578,63.7040151706876}},
	+ {PDG(21), {-2.39761243200898,3.60510082430021,-2.35633486347936,4.92926072716207}},
	+ {PDG(21), {-2.71615068832628,-7.86392074072438,-3.82343187607456,9.15627409425217}},
	+ {PDG(21), {3.71481143654212,-0.593826015910158,-0.571807737874626,3.80518297001986}},
	+ {PDG(21), {-1409.73537044466,-164.832918114848,-5.40793033852321,1419.34948173012}},
	+ {PDG(25), {1432.99870231335,141.007847270011,-2039.81468857762,2499.9671318991}}
	+ }
	+},
	+Event{
	+ 8.84444580897612e-11,
	+ {PDG(2), {0,0,-1355.85516458651,1355.85516458651}},
	+ {PDG(21), {0,0,762.973068470164,762.973068470164}},
	+ {
	+ {PDG(2), {31.6920661491847,-20.0650356752864,-88.0817739086889,95.7360517692066}},
	+ {PDG(21), {10.7374500960332,5.21496639984996,-23.5222614151694,26.3777461357094}},
	+ {PDG(21), {11.3191440860765,-1.0583249820641,-21.0506660864576,23.9243310730427}},
	+ {PDG(21), {3.55346622033097,-10.5928500980898,-14.7616526205119,18.5132920753793}},
	+ {PDG(21), {-1.90387732688831,0.857758277367415,-2.69850669970789,3.41210148538991}},
	+ {PDG(21), {7.35005095379305,6.70445492129264,-11.3926092176962,15.1249631272984}},
	+ {PDG(21), {11.6288204797665,0.660707344701538,0.00398251555856957,11.6475755334065}},
	+ {PDG(21), {-0.597325538999893,0.0140750129896864,0.868140965166798,1.05388075271015}},
	+ {PDG(21), {0.181875741278548,-1.15024801659397,3.16189939830898,3.36953366060805}},
	+ {PDG(21), {1.89639812168121,-9.37996926770805,29.1000022474337,30.6331565480887}},
	+ {PDG(21), {2.20135893465102,1.9395929828194,12.3445618061109,12.6884281249909}},
	+ {PDG(21), {-81.5865436002047,27.0279219747078,693.765293761538,699.068777368102}},
	+ {PDG(25), {3.52711568329717,-0.173048873986144,-1170.61850686223,1177.27839540274}}
	+ }
	+},
	+Event{
	+ 2.82897015974001e-12,
	+ {PDG(2), {0,0,-875.812694091832,875.812694091832}},
	+ {PDG(21), {0,0,57.0924203359919,57.0924203359919}},
	+ {
	+ {PDG(2), {-50.099164226706,-17.95306742076,-12.3571288507816,54.6345817176597}},
	+ {PDG(21), {-28.3876522119705,-27.1230636296487,15.3765217853456,42.1658250361709}},
	+ {PDG(25), {78.4868164386766,45.0761310504086,-821.739666690405,836.104707673994}}
	+ }
	+},
	+Event{
	+ 3.88530103257739e-12,
	+ {PDG(1), {0,0,-692.199788415069,692.199788415069}},
	+ {PDG(21), {0,0,133.544525386842,133.544525386842}},
	+ {
	+ {PDG(1), {92.3817704071039,-159.70188645106,-422.314076911156,460.856011783398}},
	+ {PDG(21), {-18.9438002627979,-66.3822494147389,-17.6432749454592,71.2513561734684}},
	+ {PDG(25), {-73.4379701443061,226.084135865799,-118.697911171612,293.636945845044}}
	+ }
	+},
	+Event{
	+ 1.89512471026535e-16,
	+ {PDG(21), {0,0,-9.70721396206708,9.70721396206708}},
	+ {PDG(2), {0,0,4972.11152387058,4972.11152387058}},
	+ {
	+ {PDG(21), {59.0285496032886,24.0995479860389,243.772685243208,251.9727762122}},
	+ {PDG(21), {-41.8958472734926,34.7258546036119,541.854498089259,544.580062154572}},
	+ {PDG(21), {3.57059106888938,2.52571646985078,62.0914041344665,62.2452474624124}},
	+ {PDG(21), {23.6493805334707,11.3814642200049,911.158676703928,911.536595568309}},
	+ {PDG(2), {14.1476188418534,-46.7808803752718,1859.29459989825,1859.93683094303}},
	+ {PDG(25), {-58.5002927740095,-25.9517029042347,1344.23244583939,1351.54722549212}}
	+ }
	+},
	+Event{
	+ 1.06230333729696e-12,
	+ {PDG(1), {0,0,-673.267061150677,673.267061150677}},
	+ {PDG(2), {0,0,3753.4734253635,3753.4734253635}},
	+ {
	+ {PDG(1), {-51.2138885211792,29.9443154124237,-16.45848609113,61.5662745940536}},
	+ {PDG(21), {-14.96107972804,7.9809570994548,-3.51009205672024,17.3161869099048}},
	+ {PDG(21), {-1.51150764857182,0.795036276751962,0.160919751702901,1.71541051049912}},
	+ {PDG(21), {-3.43192788833337,-4.1044326662858,0.982667470863932,5.4396812315694}},
	+ {PDG(21), {-0.869686723042832,-0.544413629987343,0.508247367872377,1.14501379192662}},
	+ {PDG(21), {-0.383436377997686,3.68633256199393,5.04487446673279,6.2599384660459}},
	+ {PDG(2), {1169.51285589063,-196.825208114573,2985.52041775983,3212.44963345787}},
	+ {PDG(25), {-1097.14132900347,159.067413060222,107.957815543675,1120.84834755231}}
	+ }
	+},
	+Event{
	+ 2.43248258895893e-11,
	+ {PDG(1), {0,0,-3609.43952382001,3609.43952382001}},
	+ {PDG(4), {0,0,328.844648852218,328.844648852218}},
	+ {
	+ {PDG(1), {21.7193670553475,-75.5800911872304,-3123.09807448816,3124.08797314685}},
	+ {PDG(21), {-0.744850169611073,-0.816712982616582,-23.7967948014986,23.8224529529555}},
	+ {PDG(21), {1.33787759560138,-2.92642914999928,-66.1568241949618,66.2350306982179}},
	+ {PDG(21), {-13.9922568724697,-6.99798415971765,-238.679704284382,239.191881701626}},
	+ {PDG(21), {-0.0730667391275246,-0.219917787096389,-2.372118027567,2.38341069019934}},
	+ {PDG(21), {2.51086763666141,3.05877939884658,-29.2362056757196,29.502818679113}},
	+ {PDG(21), {-0.681958002353936,6.08642606895845,-14.3728250017642,15.6233078296887}},
	+ {PDG(21), {12.915385818678,10.7387096403812,-8.22435034444319,18.7020590891449}},
	+ {PDG(21), {-7.87667931706262,-0.396784235522235,-0.425892319795344,7.89815795369159}},
	+ {PDG(21), {17.7870103239941,-2.46823517076154,-0.879687189993476,17.9789813581386}},
	+ {PDG(21), {0.115334200308525,-0.592637803160529,0.0317327885691422,0.604589541231171}},
	+ {PDG(21), {-0.0151307646454415,-19.7198295865988,5.89589951058375,20.5823598963637}},
	+ {PDG(21), {-10.5374066407399,7.13293776452366,6.79661471339502,14.4260081598207}},
	+ {PDG(21), {-2.36827797535072,-4.24959031195396,4.61431801579641,6.70519866513197}},
	+ {PDG(21), {1.20825309197649,-3.46469020292517,3.61125199451062,5.14830988814517}},
	+ {PDG(21), {-7.52809896447883,-1.63174532073958,11.512786252748,13.8520436799847}},
	+ {PDG(4), {-15.5789385981766,110.227142121701,169.534018134277,202.816442856384}},
	+ {PDG(25), {1.80256832144892,-18.1793470960883,24.6509799506125,128.723145885539}}
	+ }
	+},
	+Event{
	+ 3.78086571500299e-12,
	+ {PDG(1), {0,0,-515.75612918479,515.75612918479}},
	+ {PDG(21), {0,0,45.255197092535,45.255197092535}},
	+ {
	+ {PDG(1), {-44.9232695261218,-78.4343251675974,-184.938188471832,205.84503167348}},
	+ {PDG(21), {-0.789896787696013,-3.49425776005511,-2.86038855949182,4.58427714478406}},
	+ {PDG(21), {39.2341856799301,-5.96018061260582,-15.2932404589689,42.5291462721503}},
	+ {PDG(25), {6.47898063388767,87.8887635402583,-267.409114601962,308.05287118691}}
	+ }
	+},
	+Event{
	+ 1.98244708559733e-10,
	+ {PDG(-3), {0,0,-224.330383335872,224.330383335872}},
	+ {PDG(-2), {0,0,170.771028822514,170.771028822514}},
	+ {
	+ {PDG(-3), {32.8202531191564,13.1037487329351,-119.069331793241,124.202991185198}},
	+ {PDG(21), {0.37564059390238,1.41038391799038,-3.40407009672132,3.70377940424447}},
	+ {PDG(21), {0.63521472398526,0.441271590606225,-1.10804403157733,1.35128825132093}},
	+ {PDG(21), {8.87538401212426,2.63951102873116,-9.34563787541138,13.156002703397}},
	+ {PDG(21), {-3.99498503716358,4.33898221800619,-2.91212769442599,6.57777772838043}},
	+ {PDG(21), {1.75917803287605,2.21910282041749,-0.317694605292067,2.84957093983817}},
	+ {PDG(21), {6.12933969889579,9.27724272383803,-0.192024547062053,11.1208323037505}},
	+ {PDG(-2), {53.7605785018433,-8.14883691362613,25.3268608158271,59.9837746618368}},
	+ {PDG(25), {-100.36060364562,-25.2814061188984,57.4627153145456,172.15539498042}}
	+ }
	+},
	+Event{
	+ 3.36135186973605e-14,
	+ {PDG(2), {0,0,-83.6308005825699,83.6308005825699}},
	+ {PDG(2), {0,0,1594.9580341022,1594.9580341022}},
	+ {
	+ {PDG(2), {69.2700931296573,-24.537962284628,19.6216085762928,76.0622437112821}},
	+ {PDG(21), {-0.640411544448317,-0.121536864300778,0.203901139034718,0.682988894599612}},
	+ {PDG(21), {20.387552719669,23.8741715428291,15.990208091389,35.2323023312362}},
	+ {PDG(21), {10.4247649761852,-16.4889062247163,19.2708744312976,27.4212755838058}},
	+ {PDG(2), {114.447395613312,-132.326163608661,300.055338986615,347.335034789796}},
	+ {PDG(25), {-213.889394894376,149.600397439477,1156.185302295,1191.85498937405}}
	+ }
	+},
	+Event{
	+ 1.24164842721612e-16,
	+ {PDG(21), {0,0,-4163.3968380791,4163.3968380791}},
	+ {PDG(21), {0,0,427.036915807944,427.036915807944}},
	+ {
	+ {PDG(21), {106.732103151261,-343.085152631255,-1392.48829180645,1438.09693922807}},
	+ {PDG(21), {1.92934827926003,0.64112295647063,-4.93196610137752,5.33457712036617}},
	+ {PDG(21), {-102.542946664284,-26.6613393533565,-172.654373495821,202.572001061212}},
	+ {PDG(21), {14.2652496363979,-0.43151854879601,7.04812099346319,15.9172725359969}},
	+ {PDG(21), {9.87774103876447,0.17423949094779,13.357608189581,16.614024917986}},
	+ {PDG(21), {-58.6774462908419,-1.89814743385483,97.5905399225204,113.888362659511}},
	+ {PDG(21), {-60.6174182878834,96.7859027541654,233.889248084702,260.280930422215}},
	+ {PDG(25), {89.0333691373258,274.474892765678,-2518.17080805777,2537.72964594168}}
	+ }
	+},
	+Event{
	+ 4.37143657120759e-10,
	+ {PDG(2), {0,0,-2136.85584933779,2136.85584933779}},
	+ {PDG(21), {0,0,1535.0933385543,1535.0933385543}},
	+ {
	+ {PDG(2), {104.503496367527,-32.3872132573183,82.4597838109833,137.001928022078}},
	+ {PDG(21), {2.68327599459302,-18.174655889621,15.8982219172697,24.2955046644446}},
	+ {PDG(21), {-0.854731047516613,-0.031774996308675,3.40171957095986,3.50760186643676}},
	+ {PDG(21), {0.372235882566706,2.23324951123747,12.081335879082,12.2916491796222}},
	+ {PDG(21), {4.38990156130304,-3.10557717749447,31.3124101666806,31.7707865148118}},
	+ {PDG(21), {1.28104907274355,4.89268663824963,65.3503567790372,65.5457748479209}},
	+ {PDG(21), {-28.2445918369738,-71.2536590425608,1287.14707312343,1289.42717077981}},
	+ {PDG(25), {-84.1306359942426,117.826944213816,-2099.41341203094,2108.10877201698}}
	+ }
	+},
	+Event{
	+ 9.31368951066563e-16,
	+ {PDG(21), {0,0,-698.439847271152,698.439847271152}},
	+ {PDG(21), {0,0,90.9552739539947,90.9552739539947}},
	+ {
	+ {PDG(21), {-37.5373123292732,52.3552562658943,-200.837990844143,210.917095660558}},
	+ {PDG(21), {27.9237075287326,-13.3157432192522,-70.3264379158699,76.8300093034354}},
	+ {PDG(21), {1.36880790725906,-2.95181851908393,-4.17344113102872,5.29192578658696}},
	+ {PDG(21), {6.17710679155582,2.71343481048643,-0.967530110525438,6.8158265309365}},
	+ {PDG(21), {-2.82366428666703,6.80703234929341,1.12115792335197,7.45424473016687}},
	+ {PDG(21), {12.7669535786477,23.2663366046878,15.7199663244537,30.8453377988301}},
	+ {PDG(21), {-32.4487752025597,16.4243437934552,27.5226017914867,45.6089430985263}},
	+ {PDG(25), {24.5731760123048,-85.298842085481,-375.542899354883,405.631738316107}}
	+ }
	+},
	+Event{
	+ 2.36452558089734e-10,
	+ {PDG(21), {0,0,-664.049185150277,664.049185150277}},
	+ {PDG(21), {0,0,146.509751294701,146.509751294701}},
	+ {
	+ {PDG(21), {-180.981929813298,-39.6342396414616,-472.959038678737,507.952344358028}},
	+ {PDG(21), {0.57880540710976,-0.0675151313985046,-1.26903657751655,1.39643396812817}},
	+ {PDG(21), {2.86132453451057,3.15954361837591,-6.79137957538417,8.01827478359013}},
	+ {PDG(21), {89.5755490131266,-51.6291981844296,24.9090752907339,106.347614538407}},
	+ {PDG(25), {87.9662508585508,88.1714093389138,-61.4290543146726,186.844268796825}}
	+ }
	+},
	+Event{
	+ 1.20964859638042e-13,
	+ {PDG(21), {0,0,-8.1373403430798,8.1373403430798}},
	+ {PDG(1), {0,0,3295.61434986056,3295.61434986056}},
	+ {
	+ {PDG(21), {-16.5854415353967,52.6781851700513,275.050222907465,280.54000282523}},
	+ {PDG(21), {0.465988539586193,2.52758693128769,12.8172644674686,13.0724178881825}},
	+ {PDG(21), {-23.7362754046221,8.03907391298122,217.556106447963,218.994741791349}},
	+ {PDG(21), {-5.19074973569163,-2.99613400850315,87.7851089089995,87.9894655511692}},
	+ {PDG(21), {2.41443481716938,-1.80587486542991,47.3139626381618,47.4099329259315}},
	+ {PDG(1), {-9.81510053951171,42.2382563798738,768.683460256626,769.905623159561}},
	+ {PDG(25), {52.4471438584665,-100.681093520261,1878.2708838908,1885.83950606222}}
	+ }
	+},
	+Event{
	+ 1.46852610716719e-18,
	+ {PDG(2), {0,0,-5714.75573186164,5714.75573186164}},
	+ {PDG(21), {0,0,77.6740669371232,77.6740669371232}},
	+ {
	+ {PDG(2), {-18.8305225798133,-52.6505424832289,-2363.15500932613,2363.81646206037}},
	+ {PDG(21), {-1.62439035792315,-0.269354344074759,-14.1396125203896,14.2351620301402}},
	+ {PDG(21), {2.0007959876589,1.77946988407915,-13.8437148468598,14.1002885932786}},
	+ {PDG(21), {30.0490753374022,14.2799607788029,-114.222430529362,118.969020524302}},
	+ {PDG(21), {20.6411627836056,-8.99132803899525,-56.7006517033978,61.0070937232233}},
	+ {PDG(21), {3.92866444093065,-1.08877623852802,-9.81966066742811,10.6322891989661}},
	+ {PDG(21), {-42.5786890965395,66.5196561269822,-99.9443242306918,127.383976081615}},
	+ {PDG(21), {5.3784867294076,-1.17310107831838,-6.39552783455823,8.43842769246585}},
	+ {PDG(21), {69.3178595092519,-35.7774407699638,-79.688278907507,111.513284904684}},
	+ {PDG(21), {-123.887209236263,-57.5035034857678,-78.9249842760433,157.746146287195}},
	+ {PDG(25), {55.6047664822826,74.8749596490128,-2800.24747008214,2804.58764770252}}
	+ }
	+},
	+Event{
	+ 6.69354348520286e-12,
	+ {PDG(21), {0,0,-363.55197745489,363.55197745489}},
	+ {PDG(21), {0,0,558.719769887286,558.719769887286}},
	+ {
	+ {PDG(21), {-142.284916508267,43.963017997484,-115.883016401246,188.697159245875}},
	+ {PDG(21), {-0.93450438496766,-2.21313206858824,-0.523804457924861,2.4587848845882}},
	+ {PDG(21), {-4.29462885949269,-0.836460930276429,4.50510309400993,6.28008422048018}},
	+ {PDG(21), {15.1772362339943,0.866443360588954,44.9789198091798,47.4784419605383}},
	+ {PDG(21), {-12.744502710776,-2.27910996056292,40.1995842043464,42.2329641602345}},
	+ {PDG(21), {-46.9587515102193,-18.4240424401213,362.146630288179,365.642928978965}},
	+ {PDG(25), {192.040067739729,-21.076715958524,-140.255624104148,269.481383891494}}
	+ }
	+},
	+Event{
	+ 8.3174288820234e-12,
	+ {PDG(2), {0,0,-646.963789143026,646.963789143026}},
	+ {PDG(21), {0,0,113.424378870348,113.424378870348}},
	+ {
	+ {PDG(2), {21.5504990596639,-116.39970001518,-411.977621677728,428.647728253145}},
	+ {PDG(21), {-42.8281857584533,-60.5461885723979,13.7513387178484,75.4267443453443}},
	+ {PDG(25), {21.2776866987894,176.945888587578,-135.313127312798,256.313695414885}}
	+ }
	+},
	+Event{
	+ 2.94783281607854e-13,
	+ {PDG(2), {0,0,-2816.05204722295,2816.05204722295}},
	+ {PDG(21), {0,0,1423.09979850593,1423.09979850593}},
	+ {
	+ {PDG(2), {-33.3470005399893,22.6755832627906,-932.741855747253,933.613182203502}},
	+ {PDG(21), {5.5405939109513,-4.66908983780988,-19.7483299589629,21.0355679021909}},
	+ {PDG(21), {30.3810352506951,-93.8195958005309,-228.742509383718,249.09488043362}},
	+ {PDG(21), {-1.46431117893921,-2.49975044956676,-4.64082450471947,5.47085108758045}},
	+ {PDG(21), {0.286951244640676,-2.47793418066104,-3.61748074003404,4.39416268758979}},
	+ {PDG(21), {-7.82558934948559,-10.7464611198216,-10.6310266736306,17.0218977614737}},
	+ {PDG(21), {399.183480515315,334.175375092356,1321.55855569806,1420.40052399779}},
	+ {PDG(25), {-392.755159853188,-242.638126966757,-1514.38877740677,1588.12077965514}}
	+ }
	+},
	+Event{
	+ 8.88417087090488e-16,
	+ {PDG(-2), {0,0,-1851.29062696483,1851.29062696483}},
	+ {PDG(21), {0,0,186.280215398404,186.280215398404}},
	+ {
	+ {PDG(-2), {28.8723070834901,-56.1218938140198,-823.234524059411,825.65026414634}},
	+ {PDG(21), {0.384283890138024,0.830090317808507,-4.07704475717608,4.17839897520043}},
	+ {PDG(21), {5.8754348302866,1.41366020070531,-17.4804367841423,18.49553567141}},
	+ {PDG(21), {3.82454019498235,-2.95888721338262,-12.4954846268925,13.3984796639567}},
	+ {PDG(21), {-9.02267114154693,-10.1599797482432,-17.3779502711594,22.0596223594173}},
	+ {PDG(21), {-7.88986624251134,-11.5202620391288,-6.4291628246172,15.3720708234236}},
	+ {PDG(21), {-7.09200177508962,-24.3334528415271,4.61903043410503,25.7633239027864}},
	+ {PDG(21), {-42.6636865852933,-23.0491377478917,50.4772658471066,69.9957660960211}},
	+ {PDG(21), {-1.80081232075941,15.5565268006063,16.5311957613373,22.7712293127147}},
	+ {PDG(21), {-39.1600132047787,-10.8237568918042,67.1090858551584,78.4492813973866}},
	+ {PDG(25), {68.6724852710823,121.167092976877,-922.652386140737,941.436870014579}}
	+ }
	+},
	+Event{
	+ 2.27899854271868e-13,
	+ {PDG(21), {0,0,-1258.96355189766,1258.96355189766}},
	+ {PDG(21), {0,0,81.176311893459,81.176311893459}},
	+ {
	+ {PDG(21), {-41.9932914021352,-17.5759899460602,-719.625550036912,721.063994532591}},
	+ {PDG(21), {17.8538286650283,-1.33872983518764,-186.734122603116,187.59046868092}},
	+ {PDG(21), {8.12488322996903,-3.54242517420957,-74.0820567544673,74.6104123871971}},
	+ {PDG(21), {5.02326173422592,4.53278078162111,-26.1549372807581,27.0159213099454}},
	+ {PDG(21), {2.010676181176,0.0895780065556372,-6.4834859594208,6.78869892621654}},
	+ {PDG(21), {-1.51519968702226,-6.96114066829676,-5.5449537476347,9.02772515968311}},
	+ {PDG(21), {-81.0198955165549,-11.2215177464934,9.34693505445439,82.3256407503794}},
	+ {PDG(25), {91.5157367953132,36.0174445820708,-168.509068676352,231.717002044191}}
	+ }
	+},
	+Event{
	+ 1.62897655848059e-12,
	+ {PDG(21), {0,0,-728.028572277565,728.028572277565}},
	+ {PDG(2), {0,0,5471.40544830623,5471.40544830623}},
	+ {
	+ {PDG(21), {-349.486945492137,244.211353481189,-564.788495891036,707.648468753658}},
	+ {PDG(21), {-0.998171526577376,-2.87842558819181,-2.9957497751261,4.27272711255882}},
	+ {PDG(21), {-2.18567107401271,0.123664193098539,-1.97315179196501,2.94716454758897}},
	+ {PDG(21), {0.210441507117186,-1.76348235906386,-1.41608475680581,2.27144264666862}},
	+ {PDG(21), {-39.0408466663652,-4.00706680002269,-21.516865128205,44.757343282548}},
	+ {PDG(21), {-5.74951807667289,-20.2509774201203,-2.04223453304746,21.1501717835147}},
	+ {PDG(21), {-0.288339610595367,-4.34190624651721,-0.16319191991985,4.35452881346658}},
	+ {PDG(21), {-0.673662507441454,-0.788923270724528,-0.0141101935888307,1.0375067221884}},
	+ {PDG(21), {-24.8456498835464,-13.8150835485464,0.2209610268658,28.4290639199478}},
	+ {PDG(2), {-63.5398575176866,-21.6044476238811,120.314641580918,137.766754441141}},
	+ {PDG(25), {486.598220847918,-174.884704817219,5217.75115741058,5244.79884856052}}
	+ }
	+},
	+Event{
	+ 8.20436158138027e-07,
	+ {PDG(21), {0,0,-209.379297853835,209.379297853835}},
	+ {PDG(21), {0,0,846.606619135812,846.606619135812}},
	+ {
	+ {PDG(21), {52.6432900514396,43.6321210819995,-169.464498912454,182.738322114447}},
	+ {PDG(21), {2.99782511715041,0.0198972711790634,-6.46510403565276,7.12635401353538}},
	+ {PDG(21), {0.537242823558143,0.0616236447709689,-1.06040989951642,1.19033460844084}},
	+ {PDG(21), {5.44352610355536,9.37833798368214,-5.65913864968275,12.2315595911804}},
	+ {PDG(21), {-3.27352348464913,-0.30889202684913,-0.0284978601558312,3.28818831833486}},
	+ {PDG(21), {-33.75658978502,-11.8419390620898,25.7704101413361,44.0891473439558}},
	+ {PDG(25), {-24.5917708260343,-40.9411488926927,794.134560498101,805.322010999753}}
	+ }
	+},
	+Event{
	+ 8.65827810734935e-15,
	+ {PDG(2), {0,0,-3952.66731044239,3952.66731044239}},
	+ {PDG(21), {0,0,12.1104637453713,12.1104637453713}},
	+ {
	+ {PDG(2), {-12.1495546786762,50.4831663564115,-2500.38371270986,2500.92280419628}},
	+ {PDG(21), {-7.32621475208922,0.0757155910718654,-255.164934976177,255.270098516926}},
	+ {PDG(21), {0.110760540399791,-13.5727211001368,-351.352792300759,351.614868405049}},
	+ {PDG(21), {-10.4566354158336,0.246553264119039,-73.1353297869597,73.8794861634748}},
	+ {PDG(21), {14.9367180687427,-0.246124492175714,-73.5151099290482,75.0175813513791}},
	+ {PDG(21), {-27.862967306244,-26.5044093700876,-133.842504516773,139.257476203213}},
	+ {PDG(25), {42.7478935437005,-10.4821802492022,-553.162462477451,568.815459351445}}
	+ }
	+},
	+Event{
	+ 7.60567299762933e-15,
	+ {PDG(-2), {0,0,-2284.43915680122,2284.43915680122}},
	+ {PDG(21), {0,0,21.9134732742032,21.9134732742032}},
	+ {
	+ {PDG(-2), {-41.9672742572463,17.5358016673363,-206.838639435157,211.780497714228}},
	+ {PDG(21), {0.497690219900095,-0.653121114279033,-3.06325520816853,3.17140272044925}},
	+ {PDG(21), {2.61154217477597,-3.77261039786484,-13.6322154739765,14.3836727046202}},
	+ {PDG(21), {0.63866619657849,-10.1281419798583,-28.2064184787692,29.9764774093667}},
	+ {PDG(21), {3.09140279142949,-2.79088869674741,-11.5721117062287,12.2987641768564}},
	+ {PDG(21), {13.1982468647104,-34.0640014053258,-99.9785870782433,106.443730610148}},
	+ {PDG(21), {4.65063703385033,5.37170869834737,-9.26614121782567,11.6766884102133}},
	+ {PDG(21), {-37.1443101714842,-7.86605243333928,-23.1304929011714,44.4589053042313}},
	+ {PDG(25), {54.4233991474857,36.3673056617309,-1866.83782202747,1872.16249102531}}
	+ }
	+},
	+Event{
	+ 2.79469396816861e-14,
	+ {PDG(21), {0,0,-234.137906817215,234.137906817215}},
	+ {PDG(1), {0,0,5464.65616882501,5464.65616882501}},
	+ {
	+ {PDG(21), {-30.9409068565341,-44.1316916488309,-92.5170176968088,107.07167920792}},
	+ {PDG(21), {-9.49202136819425,3.66757658764439,-5.85864536592669,11.7419467382623}},
	+ {PDG(21), {8.19685511086654,14.5833371724506,57.3193947721863,59.7107626298949}},
	+ {PDG(21), {8.7438973655622,-18.6047438202084,121.455765716908,123.183177664953}},
	+ {PDG(21), {7.02790232157018,-4.41019320059009,122.499743518413,122.780407138863}},
	+ {PDG(21), {-9.23386216732118,-10.2352103312775,298.511078947806,298.829195353522}},
	+ {PDG(1), {19.8782407329069,160.704600743324,4797.6163673083,4800.3483124685}},
	+ {PDG(25), {5.81989486114369,-101.573675502512,-68.5084251930821,175.128594440306}}
	+ }
	+}
	diff --git a/t/test_ME_hjets_mt174.cc b/t/test_ME_hjets_mt174.cc
	new file mode 100644
	index 0000000..0aaa81b
	--- /dev/null
	+++ b/t/test_ME_hjets_mt174.cc
	@@ -0,0 +1,78 @@
	+#include "LHEF/LHEF.h"
	+#include "RHEJ/MatrixElement.hh"
	+#include "RHEJ/Event.hh"
	+
	+struct Event{
	+ double weight;
	+ std::array<RHEJ::Sparticle, 2> incoming;
	+ std::vector<RHEJ::Sparticle> outgoing;
	+
	+ Event(
	+ double wt,
	+ RHEJ::Sparticle in1, RHEJ::Sparticle in2,
	+ std::initializer_list<RHEJ::Sparticle> out
	+ ):
	+ weight{wt},
	+ incoming{std::move(in1), std::move(in2)},
	+ outgoing{out}
	+ {}
	+};
	+
	+constexpr
	+RHEJ::pid::ParticleID PDG(int id){
	+ return static_cast<RHEJ::pid::ParticleID>(id);
	+}
	+
	+constexpr double alpha_s = 0.112654;
	+constexpr double mu = 125.;
	+constexpr double R = 0.4;
	+constexpr double min_jet_pt = 35.;
	+constexpr auto jet_def = fastjet::antikt_algorithm;
	+constexpr double ep = 1e-4;
	+constexpr double alpha_s_mH_RHEJ = 0.113559;
	+
	+void dump(Event const & ev){
	+ LHEF::Writer writer{std::cerr};
	+ std::cerr << std::setprecision(15);
	+ RHEJ::UnclusteredEvent tmp;
	+ tmp.incoming = ev.incoming;
	+ tmp.outgoing = ev.outgoing;
	+ tmp.central = {mu, mu, 0.};
	+ RHEJ::Event out_ev{std::move(tmp), {jet_def, R}, min_jet_pt};
	+ writer.hepeup = to_HEPEUP(std::move(out_ev), nullptr);
	+ writer.writeEvent();
	+}
	+
	+int main(){
	+ std::vector<Event> events = {
	+#include "ME_hjets_mt174.dat"
	+ };
	+
	+ RHEJ::HiggsCouplingSettings settings;
	+ settings.mt = 174.;
	+ settings.use_impact_factors = false;
	+ settings.mb = 5.;
	+ settings.include_bottom = true;
	+
	+ RHEJ::MatrixElement ME{
	+ {jet_def, R}, min_jet_pt, false, settings
	+ };
	+ for(size_t i = 0; i < events.size(); ++i){
	+ std::sort(
	+ begin(events[i].outgoing), end(events[i].outgoing),
	+ RHEJ::rapidity_less{}
	+ );
	+ const double our_ME = ME.tree(
	+ alpha_s, mu, events[i].incoming, events[i].outgoing, true
	+ )alpha_salpha_s/(alpha_s_mH_RHEJ*alpha_s_mH_RHEJ);
	+ const double deviation = our_ME/events[i].weight - 1.;
	+ if(std::abs(deviation) > ep){
	+ std::cerr
	+ << "Matrix element deviates by factor of " << 1+deviation << ":\n"
	+ "Is " << our_ME << ", should be " << events[i].weight << "\n"
	+ "Event:\n";
	+ dump(events[i]);
	+ return EXIT_FAILURE;
	+ }
	+ }
	+}

File Metadata

Mime Type: text/x-diff
Expires: Tue, Jan 21, 2:16 AM (1 d, 21 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4225678
Default Alt Text: (34 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions