test_hdf5.cc
No OneTemporary
Actions

Size

1 KB

Referenced Files

None

Subscribers

None

test_hdf5.cc
View Options

	/**
	* \authors The HEJ collaboration (see AUTHORS for details)
	* \date 2019
	* \copyright GPLv2 or later
	*/
	#include <iostream>

	#include "HEJ/EventReader.hh"

	int main(int argc, char** argv) {
	if(argc != 2) {
	std::cerr << "Usage: " << argv[0] << " file.hdf5\n";
	return EXIT_FAILURE;
	}
	auto reader = HEJ::make_reader(argv[1]);
	if(
	reader->heprup().EBMUP != std::make_pair(7000., 7000.)
	\|\| reader->heprup().PDFSUP != std::make_pair(13000, 13000)
	) {
	std::cerr << "Read incorrect init parameters\n";
	return EXIT_FAILURE;
	}
	int nevent = 0;
	while(reader->read_event()) {
	++nevent;
	if(reader->hepeup().NUP != 13) {
	std::cerr << "Read wrong number of particles: "
	<< reader->hepeup().NUP << " != 13 in event " << nevent;
	return EXIT_FAILURE;
	}
	for(size_t i = 0; i < 2; ++i) {
	for(size_t j = 0; j < 2; ++j) {
	if(reader->hepeup().PUP[i][j] != 0) {
	std::cerr << "Non-vanishing transverse momentum in incoming particle"
	" in event " << nevent;
	return EXIT_FAILURE;
	}
	}
	}
	}
	if(nevent != 51200) {
	std::cerr << "Wrong number of events " << nevent << " != 51200\n";
	return EXIT_FAILURE;
	}
	return EXIT_SUCCESS;
	}