diff --git a/Changes-API.md b/Changes-API.md index 50d28ee..9f1ad56 100644 --- a/Changes-API.md +++ b/Changes-API.md @@ -1,132 +1,133 @@ # Changelog for HEJ API This log lists only changes on the HEJ API. These are primarily code changes relevant for calling HEJ as an API. This file should only be read as an addition to [`Changes.md`](Changes.md), where the main features are documented. ## Version 2.2 ### 2.2.0 #### New and updated functions * New functions to decide whether particles are massive, charged, charged leptons, antiparticles. * Added `to_string` function for `Event`. * Updated function `Event::EventData.reconstruct_intermediate()` - Now requires an `EWConstants` argument - Added support for WpWp and WmWm events. - In the case of WW same-flavour the reconstruction will minimise the total off-shell momentum. -* Renamed `no_2_jets` error flag to `not_enough_jets`. +* Replaced `no_2_jets` and `bad_final_state` enumerators in `EventType` with + `invalid` and `unknown`. * Added helper functions `is_backward_g_to_h`, `is_forward_g_to_h` to detect incoming gluon to outgoing Higgs transitions. * Updated function `implemented_types()` to now be in HEJ namespace. ## Version 2.1 ### 2.1.0 #### Changes to Event class * Restructured `Event` class - `Event` can now only be build from a (new) `Event::EventData` class - Removed default constructor for `Event` - `Event::EventData` replaces the old `UnclusteredEvent` struct. - `UnclusteredEvent` is now **deprecated**, and will be removed in version 2.2.0 - Removed `Event.unclustered()` function - Added new member function `Events.parameters()`, to directly access (underlying) `Parameters` - New member functions `begin_partons`, `end_partons` (`rbegin_partons`, `rend_partons`) with aliases `cbegin_partons`, `cend_partons` (`crbegin_partons`, `crend_partons`) for constant (reversed) iterators over outgoing partons. * New function `Event::EventData.reconstruct_intermediate()` to reconstruct bosons from decays, e.g. `positron + nu_e => Wp` * Added optional Colour charges to particles (`Particle.colour`) - Colour connection in the HEJ limit can be generated via `Event.generate_colours` (automatically done in the resummation) - Colour configuration of input events can be checked with `Event.is_leading_colour` * Added function `Event.valid_hej_state` to check if event _could have_ been produced by `HEJ` according to the `soft pt regulator` cut on the jets #### New and updated functions * Renamed `EventType::nonHEJ` to `EventType::non_resummable` and `is_HEJ()` to `is_resummable()` such that run card is consistent with internal workings * Made `MatrixElement.tree_kin(...)` and `MatrixElement.tree_param(...)` public * New `EventReweighter` member function `treatment` to query the treatment with respect to resummation for the various event types. * Added auxiliary functions to obtain a `std::string` from `EventDescription` (`to_string` for human readable, and `to_simple_string` for easy parsable string) * New `get_analyses` function to read in multiple `HEJ::Analysis` at once, similar to `get_analysis` * New `get_ew_parameters` function to extract electroweak parameters from YAML configuration. #### New classes * New class `Unweighter` to do unweighting * New class `CrossSectionAccumulator` to keep track of Cross Section of the different subprocess * New template struct `Parameters` similar to old `Weights` - `Weights` are now an alias for `Parameters`. Calling `Weights` did not change - `Weights.hh` was replaced by `Parameters.hh`. The old `Weights.hh` header will be removed in version 2.2.0 * Function to multiplication and division of `EventParameters.weight` by double - This can be combined with `Parameters`, e.g. `Parameters*Weights`, see also `Events.parameters()` - Moved `EventParameters` to `Parameters.hh` header * new class `EWConstants` replaces previously hard coded `vev` - `EWConstants` have to be set in the general `Config` and the `MatrixElementConfig` #### Input/Output * New abstract `EventReader` class, as base for reading events from files - Moved LHE file reader to `HEJ::LesHouchesReader` * New (optional) function `finish()` in abstract class `EventWriter`. `finish()` is called _after_ all events are written. * Support reading (`HDF5Reader`) and writing (`HDF5Writer`) `hdf5` files * New `BufferedEventReader` class that allows to put events back into the reader. * New `SherpaLHEReader` to read Sherpa LHE files with correct weights * `get_analysis` now requires `YAML::Node` and `LHEF::HEPRUP` as arguments * Replaced `HepMCInterface` and `HepMCWriter` by `HepMCInterfaceX` and `HepMCWriterX` respectively, with `X` being the major version of HepMC (2 or 3) - Renamed `HepMCInterfaceX::init_kinematics` to `HepMCInterfaceX::init_event` and protected it, use `HepMCInterfaceX::operator()` instead - Removed redundant `HepMCInterfaceX::add_variation` function #### Linking * Export cmake target `HEJ::HEJ` to link directly against `libHEJ` * Preprocessor flags (`HEJ_BUILD_WITH_XYZ`) for enabled/disabled dependencies are now written to `ConfigFlags.hh` * Provide links to version specific object files, e.g. `libHEJ.so -> libHEJ.so.2.1 (soname) -> libHEJ.so.2.1.0` * Removed `LHAPDF` from public interface #### Miscellaneous * Capitalisation of `Config.hh` now matches class `Config` (was `config.hh`) * Renamed `Config::max_ext_soft_pt_fraction` to `Config::soft_pt_regulator`. The new `Config::soft_pt_regulator` parameter is optional and the old `Config::max_ext_soft_pt_fraction` is **deprecated**. * Replaced redundant member `EventReweighterConfig::jet_param` with getter function `EventReweighter.jet_param()` (`JetParameters` are already in `PhaseSpacePointConfig`) * Avoid storing reference to the Random Number Generator inside classes - Constructors of `EventReweighter` now expect `std::shared_ptr` (was reference) - Moved reference to `HEJ::RNG` from constructor of `JetSplitter` to `JetSplitter.split` ## Version 2.0 ### 2.0.4 * Fixed wrong path of `HEJ_INCLUDE_DIR` in `hej-config.cmake` ### 2.0.0 * First release diff --git a/include/HEJ/event_types.hh b/include/HEJ/event_types.hh index 165a552..c62633a 100644 --- a/include/HEJ/event_types.hh +++ b/include/HEJ/event_types.hh @@ -1,119 +1,113 @@ /** \file * \brief Define different types of events. * * \authors The HEJ collaboration (see AUTHORS for details) * \date 2019-2020 * \copyright GPLv2 or later */ #pragma once #include #include "HEJ/exceptions.hh" namespace HEJ { //! Namespace for event types namespace event_type { //! Possible event types enum EventType: std::size_t { non_resummable = 0, //!< event configuration not covered by All Order resummation - bad_final_state = 1, //!< event with an unsupported final state - not_enough_jets = 2, //!< event with less than two jets + invalid = 1, //!< unphysical event, e.g. violating conservation laws + unknown = 2, //!< configuration not considered by HEJ, e.g di-Higgs FKL = 4, //!< FKL-type event unordered_backward = 8, //!< event with unordered backward emission unordered_forward = 16, //!< event with unordered forward emission extremal_qqbar_backward = 32, //!< event with a backward extremal qqbar extremal_qqbar_forward = 64, //!< event with a forward extremal qqbar central_qqbar = 128, //!< event with a central qqbar unob = unordered_backward, //!< alias for unordered_backward unof = unordered_forward, //!< alias for unordered_forward qqbar_exb = extremal_qqbar_backward, //!< alias for extremal_qqbar_backward qqbar_exf = extremal_qqbar_forward, //!< alias for extremal_qqbar_forward qqbar_mid = central_qqbar, //!< alias for central_qqbar - invalid = 256, //!< unphysical event, e.g. violating conservation laws - unknown = 512, //!< configuration not considered by HEJ, e.g di-Higgs first_type = non_resummable, //!< alias for numerically smallest enumerator - last_type = unknown //!< alias for numerically largest enumerator + last_type = central_qqbar //!< alias for numerically largest enumerator }; constexpr std::size_t UNO = unordered_backward | unordered_forward; constexpr std::size_t EXTREMAL_QQBAR = extremal_qqbar_backward | extremal_qqbar_forward; constexpr std::size_t QQBAR = EXTREMAL_QQBAR | central_qqbar; constexpr auto NLL = UNO | QQBAR; constexpr auto RESUMMABLE = FKL | NLL; constexpr auto VALID = RESUMMABLE | non_resummable; //! Event type names /** * For example, name(FKL) is the string "FKL" */ inline std::string name(EventType type) { switch(type) { case FKL: return "FKL"; case unordered_backward: return "unordered backward"; case unordered_forward: return "unordered forward"; case extremal_qqbar_backward: return "extremal qqbar backward"; case extremal_qqbar_forward: return "extremal qqbar forward"; case central_qqbar: return "central qqbar"; case non_resummable: return "non-resummable"; - case not_enough_jets: - return "not enough jets"; - case bad_final_state: - return "bad final state"; case invalid: return "invalid"; case unknown: return "unknown"; default: throw std::logic_error{"Unreachable"}; } } //! Returns True for a HEJ \ref event_type::EventType "EventType" inline constexpr bool is_resummable(EventType type) { return type & RESUMMABLE; } //! Returns True for an unordered \ref event_type::EventType "EventType" inline constexpr bool is_uno(EventType type) { return type & UNO; } //! Returns True for an extremal_qqbar \ref event_type::EventType "EventType" inline constexpr bool is_ex_qqbar(EventType type) { return type & EXTREMAL_QQBAR; } //! Returns True for an central_qqbar \ref event_type::EventType "EventType" inline constexpr bool is_mid_qqbar(EventType type) { return type == central_qqbar; } //! Returns True for any qqbar event \ref event_type::EventType "EventType" inline constexpr bool is_qqbar(EventType type) { return type & QQBAR; } } // namespace event_type } // namespace HEJ diff --git a/src/YAMLreader.cc b/src/YAMLreader.cc index 5d301af..99f99b3 100644 --- a/src/YAMLreader.cc +++ b/src/YAMLreader.cc @@ -1,608 +1,606 @@ /** * \authors The HEJ collaboration (see AUTHORS for details) * \date 2019-2020 * \copyright GPLv2 or later */ #include "HEJ/YAMLreader.hh" #include #include #include #include #include #include #include #include #include "HEJ/ConfigFlags.hh" #include "HEJ/Constants.hh" #include "HEJ/ScaleFunction.hh" #include "HEJ/event_types.hh" #include "HEJ/output_formats.hh" namespace HEJ { class Event; namespace { //! Get YAML tree of supported options /** * The configuration file is checked against this tree of options * in assert_all_options_known. */ YAML::Node const & get_supported_options(){ const static YAML::Node supported = [](){ YAML::Node supported; static const auto opts = { "trials", "min extparton pt", "max ext soft pt fraction", "soft pt regulator", "scales", "scale factors", "max scale ratio", "import scales", "log correction", "event output", "analysis", "analyses", "vev", "regulator parameter", "max events", "off-shell tolerance", "require low pt jet" }; // add subnodes to "supported" - the assigned value is irrelevant for(auto && opt: opts) supported[opt] = ""; for(auto && jet_opt: {"min pt", "algorithm", "R"}){ supported["resummation jets"][jet_opt] = ""; supported["fixed order jets"][jet_opt] = ""; } for(auto && opt: {"mt", "use impact factors", "include bottom", "mb"}){ supported["Higgs coupling"][opt] = ""; } for(auto && opt: {"name", "seed"}){ supported["random generator"][opt] = ""; } for(auto && opt: {"enabled", "nlo order"}){ supported["NLO truncation"][opt] = ""; } for( auto && opt: { "FKL", "unordered", "extremal qqbar", "central qqbar", "non-resummable", "unknown", "invalid" }){ supported["event treatment"][opt] = ""; } for(auto && particle_type: {"Higgs", "W", "Z"}){ for(auto && particle_opt: {"mass", "width"}){ supported["particle properties"][particle_type][particle_opt] = ""; } } for(auto && opt: {"type", "trials", "max deviation"}){ supported["unweight"][opt] = ""; } return supported; }(); return supported; } fastjet::JetAlgorithm to_JetAlgorithm(std::string const & algo){ using namespace fastjet; static const std::map known = { {"kt", kt_algorithm}, {"cambridge", cambridge_algorithm}, {"antikt", antikt_algorithm}, {"cambridge for passive", cambridge_for_passive_algorithm}, {"plugin", plugin_algorithm} }; const auto res = known.find(algo); if(res == known.end()){ throw std::invalid_argument("Unknown jet algorithm \"" + algo + "\""); } return res->second; } EventTreatment to_EventTreatment(std::string const & name){ static const std::map known = { {"reweight", EventTreatment::reweight}, {"keep", EventTreatment::keep}, {"discard", EventTreatment::discard}, {"abort", EventTreatment::abort} }; const auto res = known.find(name); if(res == known.end()){ throw std::invalid_argument("Unknown event treatment \"" + name + "\""); } return res->second; } WeightType to_weight_type(std::string const & setting){ if(setting == "weighted") return WeightType::weighted; if(setting =="resummation") return WeightType::unweighted_resum; if(setting =="partial") return WeightType::partially_unweighted; throw std::invalid_argument{"Unknown weight type \"" + setting + "\""}; } } // namespace namespace detail{ void set_from_yaml(fastjet::JetAlgorithm & setting, YAML::Node const & yaml){ setting = to_JetAlgorithm(yaml.as()); } void set_from_yaml(EventTreatment & setting, YAML::Node const & yaml){ setting = to_EventTreatment(yaml.as()); } void set_from_yaml(ParticleID & setting, YAML::Node const & yaml){ setting = to_ParticleID(yaml.as()); } void set_from_yaml(WeightType & setting, YAML::Node const & yaml){ setting = to_weight_type(yaml.as()); } } // namespace detail JetParameters get_jet_parameters( YAML::Node const & node, std::string const & entry ){ assert(node); JetParameters result; fastjet::JetAlgorithm jet_algo = fastjet::antikt_algorithm; double R = NAN; set_from_yaml_if_defined(jet_algo, node, entry, "algorithm"); set_from_yaml(R, node, entry, "R"); result.def = fastjet::JetDefinition{jet_algo, R}; set_from_yaml(result.min_pt, node, entry, "min pt"); return result; } RNGConfig to_RNGConfig( YAML::Node const & node, std::string const & entry ){ assert(node); RNGConfig result; set_from_yaml(result.name, node, entry, "name"); set_from_yaml_if_defined(result.seed, node, entry, "seed"); return result; } NLOConfig to_NLOConfig( YAML::Node const & node, std::string const & entry ){ assert(node); NLOConfig result; set_from_yaml_if_defined(result.enabled, node, entry, "enabled"); set_from_yaml_if_defined(result.nj, node, entry, "nlo order"); return result; } ParticleProperties get_particle_properties( YAML::Node const & node, std::string const & entry, std::string const & boson ){ ParticleProperties result{}; set_from_yaml(result.mass, node, entry, boson, "mass"); set_from_yaml(result.width, node, entry, boson, "width"); return result; } EWConstants get_ew_parameters(YAML::Node const & node){ EWConstants result; double vev = NAN; set_from_yaml(vev, node, "vev"); result.set_vevWZH(vev, get_particle_properties(node, "particle properties", "W"), get_particle_properties(node, "particle properties", "Z"), get_particle_properties(node, "particle properties", "Higgs") ); return result; } HiggsCouplingSettings get_Higgs_coupling( YAML::Node const & node, std::string const & entry ){ assert(node); static constexpr double mt_max = 2e4; #ifndef HEJ_BUILD_WITH_QCDLOOP if(node[entry].IsDefined()){ throw std::invalid_argument{ "Higgs coupling settings require building HEJ 2 " "with QCDloop support" }; } #endif HiggsCouplingSettings settings; set_from_yaml_if_defined(settings.mt, node, entry, "mt"); set_from_yaml_if_defined(settings.mb, node, entry, "mb"); set_from_yaml_if_defined(settings.include_bottom, node, entry, "include bottom"); set_from_yaml_if_defined(settings.use_impact_factors, node, entry, "use impact factors"); if(settings.use_impact_factors){ if(settings.mt != std::numeric_limits::infinity()){ throw std::invalid_argument{ "Conflicting settings: " "impact factors may only be used in the infinite top mass limit" }; } } else{ // huge values of the top mass are numerically unstable settings.mt = std::min(settings.mt, mt_max); } return settings; } FileFormat to_FileFormat(std::string const & name){ static const std::map known = { {"Les Houches", FileFormat::Les_Houches}, {"HepMC", FileFormat::HepMC}, {"HepMC2", FileFormat::HepMC2}, {"HepMC3", FileFormat::HepMC3}, {"HDF5", FileFormat::HDF5} }; const auto res = known.find(name); if(res == known.end()){ throw std::invalid_argument("Unknown file format \"" + name + "\""); } return res->second; } std::string extract_suffix(std::string const & filename){ size_t separator = filename.rfind('.'); if(separator == std::string::npos) return {}; return filename.substr(separator + 1); } FileFormat format_from_suffix(std::string const & filename){ const std::string suffix = extract_suffix(filename); if(suffix == "lhe") return FileFormat::Les_Houches; if(suffix == "hepmc") return FileFormat::HepMC; if(suffix == "hepmc3") return FileFormat::HepMC3; if(suffix == "hepmc2") return FileFormat::HepMC2; if(suffix == "hdf5") return FileFormat::HDF5; throw std::invalid_argument{ "Can't determine format for output file \"" + filename + "\"" }; } void assert_all_options_known( YAML::Node const & conf, YAML::Node const & supported ){ if(!conf.IsMap()) return; if(!supported.IsMap()) throw invalid_type{"must not have sub-entries"}; for(auto const & entry: conf){ const auto name = entry.first.as(); if(! supported[name]) throw unknown_option{name}; /* check sub-options, e.g. 'resummation jets: min pt' * we don't check analyses sub-options * those depend on the analysis being used and should be checked there * similar for "import scales" */ if(name != "analyses" && name != "analysis" && name != "import scales"){ try{ assert_all_options_known(conf[name], supported[name]); } catch(unknown_option const & ex){ throw unknown_option{name + ": " + ex.what()}; } catch(invalid_type const & ex){ throw invalid_type{name + ": " + ex.what()}; } } } } } // namespace HEJ namespace YAML { Node convert::encode(HEJ::OutputFile const & outfile) { Node node; node[to_string(outfile.format)] = outfile.name; return node; } bool convert::decode(Node const & node, HEJ::OutputFile & out) { switch(node.Type()){ case NodeType::Map: { YAML::const_iterator it = node.begin(); out.format = HEJ::to_FileFormat(it->first.as()); out.name = it->second.as(); return true; } case NodeType::Scalar: out.name = node.as(); out.format = HEJ::format_from_suffix(out.name); return true; default: return false; } } } // namespace YAML namespace HEJ { namespace detail{ void set_from_yaml(OutputFile & setting, YAML::Node const & yaml){ setting = yaml.as(); } } namespace { void update_fixed_order_jet_parameters( JetParameters & fixed_order_jets, YAML::Node const & yaml ){ if(!yaml["fixed order jets"]) return; set_from_yaml_if_defined( fixed_order_jets.min_pt, yaml, "fixed order jets", "min pt" ); fastjet::JetAlgorithm algo = fixed_order_jets.def.jet_algorithm(); set_from_yaml_if_defined(algo, yaml, "fixed order jets", "algorithm"); double R = fixed_order_jets.def.R(); set_from_yaml_if_defined(R, yaml, "fixed order jets", "R"); fixed_order_jets.def = fastjet::JetDefinition{algo, R}; } // like std::stod, but throw if not the whole string can be converted double to_double(std::string const & str){ std::size_t pos = 0; const double result = std::stod(str, &pos); if(pos < str.size()){ throw std::invalid_argument(str + " is not a valid double value"); } return result; } using EventScale = double (*)(Event const &); void import_scale_functions( std::string const & file, std::vector const & scale_names, std::unordered_map & known ) { void * handle = dlopen(file.c_str(), RTLD_NOW); char * error = dlerror(); if(error != nullptr) throw std::runtime_error{error}; for(auto const & scale: scale_names) { void * sym = dlsym(handle, scale.c_str()); error = dlerror(); if(error != nullptr) throw std::runtime_error{error}; known.emplace(scale, reinterpret_cast(sym)); // NOLINT } } auto get_scale_map( YAML::Node const & yaml ) { std::unordered_map scale_map; scale_map.emplace("H_T", H_T); scale_map.emplace("max jet pperp", max_jet_pt); scale_map.emplace("jet invariant mass", jet_invariant_mass); scale_map.emplace("m_j1j2", m_j1j2); if(yaml["import scales"].IsDefined()) { if(! yaml["import scales"].IsMap()) { throw invalid_type{"Entry 'import scales' is not a map"}; } for(auto const & import: yaml["import scales"]) { const auto file = import.first.as(); const auto scale_names = import.second.IsSequence() ?import.second.as>() :std::vector{import.second.as()}; import_scale_functions(file, scale_names, scale_map); } } return scale_map; } // simple (as in non-composite) scale functions /** * An example for a simple scale function would be H_T, * H_T/2 is then composite (take H_T and then divide by 2) */ ScaleFunction parse_simple_ScaleFunction( std::string const & scale_fun, std::unordered_map const & known ) { assert( scale_fun.empty() || (!std::isspace(scale_fun.front()) && !std::isspace(scale_fun.back())) ); const auto it = known.find(scale_fun); if(it != end(known)) return {it->first, it->second}; try{ const double scale = to_double(scale_fun); return {scale_fun, FixedScale{scale}}; } catch(std::invalid_argument const &){} throw std::invalid_argument{"Unknown scale choice: \"" + scale_fun + "\""}; } std::string trim_front(std::string const & str){ const auto new_begin = std::find_if( begin(str), end(str), [](char c){ return std::isspace(c) == 0; } ); return std::string(new_begin, end(str)); } std::string trim_back(std::string str){ size_t pos = str.size() - 1; // use guaranteed wrap-around behaviour to check whether we have // traversed the whole string for(; pos < str.size() && std::isspace(str[pos]); --pos) {} str.resize(pos + 1); // note that pos + 1 can be 0 return str; } ScaleFunction parse_ScaleFunction( std::string const & scale_fun, std::unordered_map const & known ){ assert( scale_fun.empty() || (!std::isspace(scale_fun.front()) && !std::isspace(scale_fun.back())) ); // parse from right to left => a/b/c gives (a/b)/c const size_t delim = scale_fun.find_last_of("*/"); if(delim == std::string::npos){ return parse_simple_ScaleFunction(scale_fun, known); } const std::string first = trim_back(std::string{scale_fun, 0, delim}); const std::string second = trim_front(std::string{scale_fun, delim+1}); if(scale_fun[delim] == '/'){ return parse_ScaleFunction(first, known) / parse_ScaleFunction(second, known); } assert(scale_fun[delim] == '*'); return parse_ScaleFunction(first, known) * parse_ScaleFunction(second, known); } EventTreatMap get_event_treatment( YAML::Node const & node, std::string const & entry ){ using namespace event_type; EventTreatMap treat { - {not_enough_jets, EventTreatment::discard}, - {bad_final_state, EventTreatment::discard}, {FKL, EventTreatment::discard}, {unob, EventTreatment::discard}, {unof, EventTreatment::discard}, {qqbar_exb, EventTreatment::discard}, {qqbar_exf, EventTreatment::discard}, {qqbar_mid, EventTreatment::discard}, {non_resummable, EventTreatment::discard}, {unknown, EventTreatment::abort}, {invalid, EventTreatment::abort} }; set_from_yaml(treat.at(FKL), node, entry, "FKL"); set_from_yaml(treat.at(unob), node, entry, "unordered"); treat.at(unof) = treat.at(unob); set_from_yaml(treat.at(qqbar_exb), node, entry, "extremal qqbar"); treat.at(qqbar_exf) = treat.at(qqbar_exb); set_from_yaml(treat.at(qqbar_mid), node, entry, "central qqbar"); set_from_yaml(treat.at(non_resummable), node, entry, "non-resummable"); set_from_yaml(treat.at(unknown), node, entry, "unknown"); set_from_yaml(treat.at(invalid), node, entry, "invalid"); if(treat[non_resummable] == EventTreatment::reweight){ throw std::invalid_argument{"Cannot reweight non-resummable events"}; } return treat; } Config to_Config(YAML::Node const & yaml){ try{ assert_all_options_known(yaml, get_supported_options()); } catch(unknown_option const & ex){ throw unknown_option{std::string{"Unknown option '"} + ex.what() + "'"}; } Config config; config.resummation_jets = get_jet_parameters(yaml, "resummation jets"); config.fixed_order_jets = config.resummation_jets; update_fixed_order_jet_parameters(config.fixed_order_jets, yaml); set_from_yaml_if_defined(config.min_extparton_pt, yaml, "min extparton pt"); if(config.min_extparton_pt!=0) std::cerr << "WARNING: \"min extparton pt\" is deprecated." << " Please remove this entry or set \"soft pt regulator\" instead.\n"; set_from_yaml_if_defined( config.max_ext_soft_pt_fraction, yaml, "max ext soft pt fraction" ); if(config.max_ext_soft_pt_fraction){ std::cerr << "WARNING: \"max ext soft pt fraction\" is deprecated." << " Please remove this entry or set \"soft pt regulator\" instead.\n"; config.soft_pt_regulator = *config.max_ext_soft_pt_fraction; } else { set_from_yaml_if_defined( config.soft_pt_regulator, yaml, "soft pt regulator" ); } // Sets the standard value, then changes this if defined config.regulator_lambda=CLAMBDA; set_from_yaml_if_defined(config.regulator_lambda, yaml, "regulator parameter"); set_from_yaml_if_defined(config.max_events, yaml, "max events"); set_from_yaml(config.trials, yaml, "trials"); config.weight_type = WeightType::weighted; set_from_yaml_if_defined(config.weight_type, yaml, "unweight", "type"); if(config.weight_type == WeightType::partially_unweighted) { config.unweight_config = PartialUnweightConfig{}; set_from_yaml( config.unweight_config->trials, yaml, "unweight", "trials" ); set_from_yaml( config.unweight_config->max_dev, yaml, "unweight", "max deviation" ); } else if(yaml["unweight"].IsDefined()) { for(auto && opt: {"trials", "max deviation"}) { if(yaml["unweight"][opt].IsDefined()) { throw std::invalid_argument{ "'unweight: " + std::string{opt} + "' " "is only supported if 'unweight: type' is set to 'partial'" }; } } } set_from_yaml(config.log_correction, yaml, "log correction"); config.treat = get_event_treatment(yaml, "event treatment"); set_from_yaml_if_defined(config.output, yaml, "event output"); config.rng = to_RNGConfig(yaml, "random generator"); set_from_yaml_if_defined(config.lowpt, yaml, "require low pt jet"); set_from_yaml_if_defined(config.analyses_parameters, yaml, "analyses"); if(yaml["analysis"].IsDefined()){ std::cerr << "WARNING: Configuration entry 'analysis' is deprecated. " " Use 'analyses' instead.\n"; set_from_yaml(config.analysis_parameters, yaml, "analysis"); if(!config.analysis_parameters.IsNull()){ config.analyses_parameters.push_back(config.analysis_parameters); } } config.scales = to_ScaleConfig(yaml); config.ew_parameters = get_ew_parameters(yaml); config.Higgs_coupling = get_Higgs_coupling(yaml, "Higgs coupling"); //HEJ@NLO Truncation config.nlo = to_NLOConfig(yaml, "NLO truncation"); set_from_yaml_if_defined( config.off_shell_tolerance, yaml, "off-shell tolerance" ); return config; } } // namespace ScaleConfig to_ScaleConfig(YAML::Node const & yaml){ ScaleConfig config; auto scale_funs = get_scale_map(yaml); std::vector scales; set_from_yaml(scales, yaml, "scales"); config.base.reserve(scales.size()); std::transform( begin(scales), end(scales), std::back_inserter(config.base), [scale_funs](auto const & entry){ return parse_ScaleFunction(entry, scale_funs); } ); set_from_yaml_if_defined(config.factors, yaml, "scale factors"); config.max_ratio = std::numeric_limits::infinity(); set_from_yaml_if_defined(config.max_ratio, yaml, "max scale ratio"); return config; } Config load_config(std::string const & config_file){ try{ return to_Config(YAML::LoadFile(config_file)); } catch(...){ std::cerr << "Error reading " << config_file << ":\n "; throw; } } } // namespace HEJ diff --git a/t/check_res.cc b/t/check_res.cc index b01e9db..6445d36 100644 --- a/t/check_res.cc +++ b/t/check_res.cc @@ -1,153 +1,151 @@ /** * \authors The HEJ collaboration (see AUTHORS for details) * \date 2019-2020 * \copyright GPLv2 or later */ #include "hej_test.hh" #include #include #include #include #include #include #include #include #include "HEJ/Config.hh" #include "HEJ/CrossSectionAccumulator.hh" #include "HEJ/Event.hh" #include "HEJ/event_types.hh" #include "HEJ/EventReweighter.hh" #include "HEJ/EWConstants.hh" #include "HEJ/Fraction.hh" #include "HEJ/HiggsCouplingSettings.hh" #include "HEJ/make_RNG.hh" #include "HEJ/Mixmax.hh" #include "HEJ/Parameters.hh" #include "HEJ/ScaleFunction.hh" #include "HEJ/EventReader.hh" #include "HEJ/YAMLreader.hh" #include "fastjet/JetDefinition.hh" #include "LHEF/LHEF.h" namespace HEJ { struct RNG; } namespace { using EventTreatment = HEJ::EventTreatment; using namespace HEJ::event_type; HEJ::EventTreatMap TREAT{ - {not_enough_jets, EventTreatment::discard}, - {bad_final_state, EventTreatment::discard}, {non_resummable, EventTreatment::discard}, {unof, EventTreatment::discard}, {unob, EventTreatment::discard}, {qqbar_exb, EventTreatment::discard}, {qqbar_exf, EventTreatment::discard}, {qqbar_mid, EventTreatment::discard}, {FKL, EventTreatment::reweight} }; bool correct_colour(HEJ::Event const & ev){ if(!HEJ::event_type::is_resummable(ev.type())) return true; return ev.is_leading_colour(); } } // namespace int main(int argn, char** argv) { if(argn == 6 && std::string(argv[5]) == "unof"){ --argn; TREAT[unof] = EventTreatment::reweight; TREAT[unob] = EventTreatment::discard; TREAT[FKL] = EventTreatment::discard; } else if(argn == 6 && std::string(argv[5]) == "unob"){ --argn; TREAT[unof] = EventTreatment::discard; TREAT[unob] = EventTreatment::reweight; TREAT[FKL] = EventTreatment::discard; } else if(argn == 6 && std::string(argv[5]) == "splitf"){ --argn; TREAT[qqbar_exb] = EventTreatment::discard; TREAT[qqbar_exf] = EventTreatment::reweight; TREAT[FKL] = EventTreatment::discard; } else if(argn == 6 && std::string(argv[5]) == "splitb"){ --argn; TREAT[qqbar_exb] = EventTreatment::reweight; TREAT[qqbar_exf] = EventTreatment::discard; TREAT[FKL] = EventTreatment::discard; } else if(argn == 6 && std::string(argv[5]) == "qqbar_mid"){ --argn; TREAT[qqbar_mid] = EventTreatment::reweight; TREAT[FKL] = EventTreatment::discard; } if(argn != 5){ std::cerr << "Usage: check_res yaml eventfile xsection tolerance [uno] \n"; return EXIT_FAILURE; } HEJ::Config config = HEJ::load_config(argv[1]); auto reader = HEJ::make_reader(argv[2]); const double xsec_ref = std::stod(argv[3]); const double tolerance = std::stod(argv[4]); reader->read_event(); HEJ::EventReweighterConfig conf = HEJ::to_EventReweighterConfig(config); conf.treat = TREAT; HEJ::MatrixElementConfig ME_conf = HEJ::to_MatrixElementConfig(config); HEJ::PhaseSpacePointConfig psp_conf = HEJ::to_PhaseSpacePointConfig(config); const fastjet::JetDefinition BORN_JET_DEF = config.fixed_order_jets.def; const double BORN_JETPTMIN = config.fixed_order_jets.min_pt; const int NUM_TRIES = config.trials; HEJ::ScaleGenerator scale_gen{ config.scales.base, config.scales.factors, config.scales.max_ratio }; std::shared_ptr ran{ HEJ::make_RNG(config.rng.name,config.rng.seed) }; HEJ::EventReweighter hej{reader->heprup(), std::move(scale_gen), conf, ran}; HEJ::CrossSectionAccumulator xs; do{ auto ev_data = HEJ::Event::EventData{reader->hepeup()}; shuffle_particles(ev_data); ev_data.reconstruct_intermediate(ME_conf.ew_parameters); HEJ::Event ev{ ev_data.cluster( BORN_JET_DEF, BORN_JETPTMIN ) }; auto resummed_events = hej.reweight(ev, NUM_TRIES); for(auto const & res_ev: resummed_events) { ASSERT(correct_colour(res_ev)); ASSERT(std::isfinite(res_ev.central().weight)); // we fill the xs uncorrelated since we only want to test the uncertainty // of the resummation xs.fill(res_ev); } } while(reader->read_event()); const double xsec = xs.total().value; const double xsec_err = std::sqrt(xs.total().error); const double significance = std::abs(xsec - xsec_ref) / std::sqrt( xsec_err*xsec_err + tolerance*tolerance ); std::cout << xsec_ref << " +/- " << tolerance << " ~ " << xsec << " +- " << xsec_err << " => " << significance << " sigma\n"; if(significance > 3.){ std::cerr << "Cross section is off by over 3 sigma!\n"; return EXIT_FAILURE; } return EXIT_SUCCESS; } diff --git a/t/classify_ref_4j b/t/classify_ref_4j index 98687cd..7f5a5f2 100644 --- a/t/classify_ref_4j +++ b/t/classify_ref_4j @@ -1,3255 +1,3255 @@ 4 4 0 4 4 8 4 4 4 16 4 32 0 4 4 0 4 4 128 4 4 4 4 4 4 4 0 4 64 4 4 8 4 64 4 4 -512 +2 8 4 -512 +2 4 4 4 4 4 4 4 8 4 4 4 4 4 -512 +2 4 4 4 16 128 4 4 4 4 4 4 64 4 4 64 4 4 4 4 4 4 0 8 32 0 4 4 64 -512 +2 4 4 4 64 4 0 4 4 4 128 4 4 0 16 8 4 4 4 4 4 4 4 4 4 -512 +2 4 32 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 64 4 4 4 32 4 4 4 -512 +2 4 4 4 4 0 4 8 4 16 4 4 32 8 4 4 -512 +2 4 -512 +2 4 4 64 4 64 4 4 4 4 4 8 4 128 4 128 4 4 -512 +2 4 64 -512 -512 +2 +2 8 4 4 128 4 4 8 4 4 16 8 0 4 128 4 4 4 4 0 4 4 4 4 0 4 4 4 -512 +2 0 4 4 4 4 4 4 4 4 4 4 4 4 -512 +2 8 4 4 4 4 4 4 4 128 4 64 4 8 4 4 64 0 128 4 16 4 4 4 4 32 4 4 4 4 4 4 4 4 4 4 32 4 4 128 4 64 4 8 4 4 4 4 0 4 8 4 8 4 4 4 4 4 4 4 4 128 0 4 4 32 4 4 4 4 4 16 4 4 4 -512 +2 0 64 4 4 4 4 0 64 8 4 4 4 4 4 32 4 4 4 16 64 4 4 16 4 4 32 4 4 16 4 -512 +2 4 4 64 4 4 4 0 0 4 4 4 4 4 4 -512 +2 4 4 0 0 4 4 4 64 4 4 32 4 4 4 4 4 0 4 16 0 4 4 4 4 0 4 64 4 4 4 -512 +2 16 -512 +2 4 128 32 4 0 4 4 4 4 0 4 4 4 4 4 -512 +2 0 4 4 16 4 4 64 4 4 4 4 4 4 4 4 4 4 4 4 0 4 4 16 0 4 4 128 4 4 64 -512 +2 4 64 4 4 32 4 4 4 4 -512 +2 0 4 16 8 8 4 4 0 4 4 -512 +2 4 4 4 4 4 4 4 4 4 32 4 4 4 -512 +2 128 4 4 4 4 4 4 8 4 4 4 4 128 4 4 32 8 4 4 4 4 4 -512 +2 16 0 4 4 -512 +2 4 4 4 4 4 -512 +2 0 4 4 8 4 64 4 4 -512 +2 4 4 4 4 4 4 4 0 4 4 0 4 8 4 4 4 0 64 16 4 4 4 0 4 4 4 4 4 4 4 4 4 4 4 4 16 4 4 4 4 4 4 4 4 4 -512 +2 4 4 4 16 64 4 128 16 -512 +2 4 4 4 64 4 0 4 4 64 0 4 4 4 4 4 4 4 4 4 4 4 4 4 4 32 4 0 64 4 0 4 4 4 4 4 4 -512 +2 4 4 4 4 4 4 4 4 32 8 4 0 4 32 128 4 -512 +2 4 4 64 4 4 4 4 128 16 4 4 4 4 128 -512 +2 8 4 4 4 64 4 4 4 0 4 4 4 4 0 4 4 4 4 4 8 64 4 32 4 128 -512 +2 4 4 0 4 4 -512 +2 0 0 4 4 4 64 4 4 4 128 -512 +2 4 4 128 4 8 0 4 4 4 -512 +2 4 32 0 4 4 4 4 4 -512 +2 4 4 0 4 0 64 32 0 4 4 4 4 -512 -512 +2 +2 4 4 4 4 32 4 0 4 4 0 4 4 4 4 4 8 4 -512 +2 4 4 4 8 4 -512 +2 4 4 4 0 0 4 4 4 4 4 4 4 4 4 4 4 64 4 4 4 4 4 4 -512 -512 -512 +2 +2 +2 4 4 4 4 4 4 4 4 0 0 4 4 4 4 8 4 8 4 0 8 64 4 4 4 -512 -512 +2 +2 4 -512 +2 4 32 4 4 4 4 0 4 -512 +2 0 4 0 4 0 4 4 4 4 0 4 4 4 0 4 4 4 64 128 4 4 4 4 4 4 4 4 64 4 8 4 0 4 -512 +2 4 4 16 4 -512 +2 4 128 -512 +2 8 4 4 0 4 4 4 4 4 4 4 4 4 4 32 4 4 4 4 4 0 4 4 0 32 4 -512 +2 4 4 4 4 -512 +2 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 8 4 4 4 4 4 4 4 128 4 64 4 128 -512 +2 128 16 4 128 0 4 4 4 4 4 4 16 4 0 -512 +2 4 4 4 4 4 4 4 4 4 64 4 16 4 4 4 32 -512 -512 +2 +2 4 32 8 4 4 4 4 4 4 4 4 4 64 4 64 32 4 4 4 32 4 0 4 4 8 4 4 4 4 4 -512 +2 4 0 -512 +2 4 4 -512 +2 4 4 0 4 128 4 4 4 4 4 16 64 4 0 4 4 4 4 -512 +2 4 4 4 4 4 0 4 4 128 4 -512 +2 0 4 4 4 0 4 4 4 4 4 4 4 4 32 4 4 4 4 4 128 128 0 4 4 4 4 128 32 4 4 4 4 64 4 4 4 4 128 0 4 4 4 8 4 4 128 4 4 64 4 4 4 0 4 -512 +2 4 4 4 4 4 4 4 4 4 32 4 -512 +2 32 4 4 16 4 4 4 128 128 4 4 4 4 4 4 4 128 4 32 4 64 4 4 4 0 4 16 4 4 -512 +2 4 4 4 4 4 4 128 4 4 4 4 4 4 4 64 32 4 4 -512 +2 4 -512 +2 16 4 -512 +2 4 4 8 16 4 4 4 4 4 128 4 4 16 0 4 4 0 4 128 4 4 4 4 4 4 4 0 4 4 4 4 64 4 4 4 4 4 4 4 4 4 4 4 8 4 -512 +2 4 4 4 4 4 4 4 4 4 -512 +2 4 16 8 4 -512 +2 4 0 4 4 8 4 4 4 4 0 4 16 32 4 16 4 0 4 4 4 4 4 4 128 4 4 4 4 0 4 0 4 4 128 -512 +2 0 4 4 4 4 4 64 4 4 4 0 4 4 4 4 4 4 4 8 4 4 4 -512 +2 8 4 4 4 4 4 4 4 4 128 4 4 4 4 -512 +2 16 32 -512 +2 4 4 -512 +2 32 32 4 4 4 -512 +2 4 64 4 128 64 0 4 4 4 -512 +2 4 4 4 32 4 4 0 4 4 4 4 32 4 4 4 4 4 4 0 64 4 4 4 4 64 0 4 4 0 4 4 4 4 4 -512 +2 4 4 0 128 4 4 0 4 16 0 4 16 4 4 0 4 4 128 -512 +2 4 4 4 4 4 4 128 4 4 8 4 4 0 4 4 4 4 4 8 8 64 16 4 4 4 4 4 16 16 -512 +2 128 4 4 4 4 4 16 16 4 64 4 4 4 64 128 -512 +2 8 4 4 -512 +2 4 4 0 4 4 4 4 4 4 4 -512 +2 32 0 4 4 4 -512 +2 128 4 4 4 4 4 4 4 4 4 4 4 4 64 128 4 4 4 -512 +2 4 4 4 4 8 4 4 4 4 0 4 8 4 4 4 4 4 32 4 0 4 4 4 4 4 64 4 4 4 -512 +2 0 4 0 4 128 4 4 4 4 4 4 16 4 4 4 8 4 -512 +2 4 4 4 4 -512 +2 4 4 0 4 4 4 32 0 4 4 4 4 4 4 4 4 8 4 64 4 4 4 4 4 0 4 32 32 4 4 4 8 4 4 4 4 4 4 64 4 4 0 4 -512 +2 32 -512 +2 8 4 4 4 4 32 -512 +2 4 4 4 -512 +2 4 4 -512 -512 +2 +2 4 4 4 4 4 128 4 4 4 0 32 4 4 4 4 -512 +2 4 -512 +2 4 -512 +2 4 4 4 -512 +2 4 4 -512 +2 4 4 8 4 4 8 4 -512 +2 4 4 16 128 4 4 4 4 -512 +2 4 -512 +2 4 4 4 4 4 4 16 4 4 4 4 -512 +2 64 4 4 4 -512 +2 0 4 16 4 4 4 4 4 4 4 4 4 4 4 4 4 64 4 4 4 4 4 4 0 -512 +2 4 0 4 4 4 4 4 4 4 0 4 -512 +2 64 4 4 8 -512 +2 32 4 4 4 16 4 4 4 4 4 4 0 64 64 4 4 4 16 4 4 0 4 4 4 -512 +2 8 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 0 4 16 32 4 4 4 4 128 4 4 4 4 128 4 128 32 4 64 0 4 16 16 4 4 0 4 4 4 4 -512 +2 4 128 0 4 4 4 4 8 4 128 4 4 4 4 4 4 4 -512 +2 4 128 4 64 4 4 16 8 4 64 4 128 4 4 128 4 32 4 8 4 4 4 4 4 4 4 4 0 4 -512 +2 4 4 4 4 128 4 4 4 4 4 4 0 4 4 -512 +2 4 64 4 4 64 64 4 4 4 64 4 4 4 4 4 16 4 8 4 4 4 8 4 4 -512 +2 4 4 4 128 4 4 32 4 0 4 4 4 4 0 4 16 4 4 4 4 4 64 4 -512 +2 4 32 4 -512 +2 32 8 4 128 4 4 4 4 4 -512 +2 64 4 4 4 4 4 4 4 -512 +2 4 4 4 32 4 4 64 4 4 4 4 0 4 4 4 4 4 4 4 4 64 4 4 4 4 4 4 4 0 -512 +2 4 4 4 4 4 4 0 4 0 4 32 4 0 4 4 4 128 8 64 128 4 4 64 8 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 4 4 32 4 4 4 4 4 8 4 4 4 4 4 4 8 4 4 64 0 -512 +2 32 0 4 4 4 4 0 4 4 4 4 4 4 4 4 4 4 4 4 4 0 0 0 4 4 0 4 4 32 4 4 4 4 4 4 16 4 0 128 4 4 4 4 4 -512 +2 16 32 4 4 4 4 0 -512 +2 4 64 4 4 4 4 4 4 4 4 4 -512 +2 4 16 4 4 4 4 -512 +2 0 4 4 4 4 4 4 4 4 4 4 4 4 4 4 16 4 8 4 8 32 4 4 4 4 4 128 64 4 4 4 4 4 4 -512 +2 4 4 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 0 -512 +2 4 4 8 0 4 4 4 4 0 4 4 4 4 0 4 64 4 128 -512 +2 4 4 4 4 16 4 4 4 4 4 4 4 4 4 4 4 4 64 4 4 4 -512 +2 4 4 4 4 4 4 -512 +2 4 4 4 4 4 4 -512 +2 0 4 4 4 4 4 4 4 4 4 4 4 -512 +2 4 32 4 128 4 4 4 64 4 -512 +2 4 4 4 0 4 128 32 4 128 4 4 8 0 4 -512 +2 4 4 0 4 4 4 4 4 4 4 4 4 4 -512 +2 4 4 128 4 4 4 4 0 4 -512 +2 32 4 4 4 4 -512 +2 4 8 4 4 4 8 4 4 128 4 4 4 4 0 0 4 4 -512 +2 4 4 4 4 32 4 -512 +2 4 4 128 8 16 128 4 4 4 32 128 4 32 0 4 4 -512 +2 4 32 4 4 32 8 4 4 4 4 4 4 8 4 4 4 4 4 4 4 64 4 4 4 4 4 4 8 4 4 -512 -512 +2 +2 4 4 4 4 4 32 4 4 4 4 32 4 4 4 4 4 64 4 4 4 4 8 4 -512 +2 4 4 64 -512 +2 4 4 4 4 4 4 0 4 128 4 4 4 4 -512 +2 4 4 4 64 4 4 4 4 -512 +2 128 4 4 4 4 4 8 4 4 128 4 4 4 4 128 4 4 4 0 4 4 16 64 4 4 32 4 32 64 16 4 4 64 4 4 4 0 4 4 4 4 0 4 64 4 64 16 4 4 4 4 64 4 0 128 4 4 8 4 -512 +2 4 0 4 -512 +2 0 4 128 4 4 4 4 4 4 4 4 4 4 -512 +2 4 16 64 64 4 4 4 4 4 4 4 4 4 -512 +2 8 4 4 4 -512 +2 32 16 4 4 128 32 4 4 -512 +2 32 4 4 4 4 64 4 4 -512 +2 4 32 4 4 0 4 128 64 4 4 4 16 4 4 4 -512 +2 4 4 4 4 4 4 128 4 4 4 0 4 -512 +2 4 4 4 4 0 0 32 4 32 4 4 4 4 4 4 128 4 4 4 4 4 4 4 4 16 4 32 4 4 4 4 4 -512 +2 0 4 4 4 64 16 4 4 4 4 4 4 4 4 4 4 -512 +2 8 -512 +2 4 4 4 0 4 4 128 4 4 4 0 4 4 4 4 4 64 4 4 4 4 4 0 4 4 4 4 4 4 4 4 4 4 4 4 4 -512 +2 4 4 32 4 64 64 128 0 8 4 0 4 0 16 4 4 4 -512 +2 16 4 4 4 4 4 0 8 0 4 4 4 4 4 4 4 -512 +2 16 4 4 4 0 4 128 4 4 32 16 4 4 4 4 4 4 4 4 4 4 0 4 4 16 128 4 4 4 4 4 4 4 4 16 4 4 4 -512 +2 4 -512 +2 4 -512 +2 4 4 4 4 4 4 4 32 32 4 4 16 4 -512 +2 0 4 4 -512 +2 4 4 4 4 16 4 128 4 4 4 4 64 0 4 128 4 4 4 4 4 8 4 4 4 4 4 4 4 4 4 16 4 4 4 4 4 4 16 32 64 0 4 0 8 4 4 4 4 4 4 4 4 32 4 4 4 4 0 4 4 128 -512 +2 16 -512 +2 4 32 4 16 0 4 4 4 4 4 4 4 64 8 4 0 4 4 64 4 4 4 4 4 64 -512 +2 4 32 4 0 4 128 4 4 16 4 0 4 4 4 64 4 4 4 4 128 4 -512 +2 4 4 4 128 -512 +2 4 4 4 16 32 4 4 4 4 64 4 4 -512 +2 4 128 -512 +2 4 0 4 -512 +2 4 4 0 4 64 4 4 4 0 128 8 16 4 -512 +2 4 4 4 4 4 8 4 4 4 4 4 4 4 4 4 16 4 4 4 4 32 4 128 4 4 4 4 16 4 4 4 64 4 4 4 4 4 4 4 4 4 4 0 4 4 4 4 4 64 4 4 64 4 8 4 4 4 4 4 4 4 4 4 16 4 4 4 8 16 -512 +2 4 4 32 4 128 4 4 4 4 4 4 4 4 0 8 -512 +2 4 4 128 4 4 4 4 4 4 64 8 4 4 4 4 0 16 4 4 4 -512 +2 4 4 4 4 4 4 4 32 0 4 4 4 4 4 4 4 0 32 4 64 4 4 4 4 4 4 4 4 0 4 4 4 4 4 -512 +2 4 128 4 4 4 64 4 128 0 4 4 8 16 4 8 4 32 4 4 16 -512 +2 4 4 4 4 4 4 4 16 4 4 0 4 4 4 4 4 4 8 4 0 4 4 4 -512 +2 4 -512 +2 4 4 128 4 4 0 4 4 4 4 4 -512 +2 16 4 16 4 4 4 4 8 4 4 4 4 4 4 4 4 4 4 4 4 4 0 4 4 4 4 4 -512 +2 4 4 4 0 16 4 4 4 4 4 4 -512 +2 4 4 16 4 4 64 64 -512 +2 4 4 16 4 128 4 4 4 4 4 128 -512 +2 4 4 4 4 4 128 4 4 8 4 32 128 4 16 -512 +2 4 4 4 -512 +2 4 4 32 4 -512 +2 4 4 4 4 16 4 4 4 4 8 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 128 16 32 4 8 4 4 4 4 4 4 4 32 128 4 0 32 4 4 4 4 0 4 128 4 4 4 128 4 4 4 -512 +2 4 4 4 4 4 128 4 0 4 4 4 64 0 -512 +2 0 4 32 16 4 4 4 4 -512 +2 4 64 4 64 4 4 32 128 4 0 4 4 32 4 4 4 4 128 0 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 16 8 4 4 16 4 4 4 4 4 4 4 4 -512 +2 4 4 128 16 4 4 16 4 4 4 4 4 4 4 4 4 0 4 4 8 4 4 -512 +2 4 4 4 4 128 4 0 -512 +2 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 64 8 4 0 4 -512 +2 4 4 4 4 8 4 diff --git a/t/classify_ref_W4j b/t/classify_ref_W4j index 3c3640a..982ef2b 100644 --- a/t/classify_ref_W4j +++ b/t/classify_ref_W4j @@ -1,3000 +1,3000 @@ 0 0 0 4 4 64 32 0 4 0 32 4 128 4 0 4 4 16 0 0 4 4 0 4 128 4 4 4 4 16 4 4 0 0 8 -512 +2 32 64 4 -512 +2 0 4 4 16 32 0 0 8 4 4 128 4 8 -512 +2 4 0 4 4 4 4 0 0 128 4 4 128 0 8 16 0 16 4 0 32 128 128 4 4 16 32 16 4 4 64 0 0 128 0 0 0 16 16 8 0 16 4 4 0 4 4 16 8 16 0 8 4 4 16 0 -512 +2 0 -512 +2 128 16 4 128 32 8 16 0 0 8 4 4 4 4 -512 +2 4 0 4 64 4 4 4 8 4 0 4 64 16 4 4 4 128 0 16 4 0 0 0 4 0 4 128 0 16 32 8 -512 +2 0 4 4 4 16 32 8 0 4 128 4 4 4 4 0 4 4 4 16 4 128 16 4 16 0 0 0 4 4 16 4 4 0 4 4 16 0 -512 +2 0 128 0 4 4 0 16 32 0 8 128 128 16 16 4 4 0 4 0 8 0 0 4 64 4 0 4 0 8 16 8 4 4 32 8 0 16 16 0 4 128 0 4 128 16 4 4 128 0 0 8 4 8 0 0 4 0 0 0 0 4 0 0 0 4 16 0 4 0 8 4 0 0 4 4 0 4 16 0 4 4 8 4 0 8 128 4 16 16 4 4 4 4 4 4 4 0 4 8 4 8 4 64 0 8 16 0 64 0 8 8 4 0 4 8 4 0 4 4 8 4 4 4 16 8 4 0 -512 +2 0 4 -512 +2 16 16 4 64 0 4 8 16 4 4 32 4 32 4 64 4 4 0 16 4 16 16 0 0 4 16 4 0 4 -512 +2 0 16 0 0 4 0 0 16 128 32 8 4 4 8 128 8 8 64 4 4 8 8 16 0 4 4 128 32 4 0 4 0 0 4 16 16 4 0 32 16 0 8 4 8 16 4 32 4 4 4 128 0 128 4 8 4 128 4 8 128 0 0 4 4 4 16 8 32 4 4 4 0 4 4 128 16 16 16 4 0 4 16 16 0 4 0 64 8 8 4 4 16 16 4 0 8 0 8 64 4 4 8 0 128 4 -512 +2 4 4 0 4 4 4 0 4 4 0 8 0 16 64 0 0 16 0 4 16 16 0 4 4 4 16 4 0 4 4 0 8 4 8 4 128 4 0 16 0 0 4 4 0 16 32 4 0 0 8 4 0 16 4 4 0 16 128 8 0 4 4 64 4 128 0 0 4 4 0 4 16 0 8 0 0 16 4 4 4 4 0 4 16 16 0 8 4 4 4 4 4 4 0 16 8 16 128 4 4 0 4 64 4 4 128 0 16 4 0 8 0 128 4 0 4 4 8 64 0 4 8 128 128 4 8 4 32 128 128 4 4 16 0 0 128 4 4 0 8 8 8 16 8 16 128 0 8 8 4 4 0 16 4 0 4 0 4 4 0 0 4 4 4 16 4 -512 +2 4 4 8 0 4 0 4 32 4 0 4 16 0 4 0 8 4 4 0 4 4 4 4 4 0 128 4 16 4 16 0 4 0 8 4 8 4 4 4 8 0 16 0 8 128 64 0 8 4 128 4 4 0 4 4 16 8 4 8 16 0 16 0 4 4 32 0 4 4 0 0 32 4 4 4 16 4 -512 +2 128 4 0 4 32 4 16 0 8 0 4 16 4 4 16 4 8 4 8 8 64 128 0 8 0 4 0 0 -512 +2 0 32 4 128 4 0 -512 +2 8 4 0 16 4 4 0 8 8 8 8 4 8 16 4 4 0 0 0 0 4 0 4 8 16 4 0 0 0 4 16 8 128 32 16 0 4 4 0 0 32 32 128 4 8 4 16 16 4 0 0 16 4 8 0 0 16 4 4 16 4 64 -512 -512 +2 +2 0 0 0 0 128 4 4 8 8 0 4 16 8 8 4 32 16 4 8 8 0 64 0 16 4 64 0 8 4 16 4 0 4 32 4 64 0 32 0 4 0 0 0 4 8 8 128 8 0 4 4 4 0 4 4 8 4 -512 +2 4 4 4 8 32 4 4 16 -512 +2 0 4 4 4 0 4 4 4 4 4 4 4 0 32 4 4 4 64 8 4 16 -512 +2 4 64 4 4 0 128 0 4 -512 -512 +2 +2 16 8 32 8 4 16 4 4 128 4 4 0 0 4 0 0 0 64 0 16 8 0 16 8 4 4 4 4 16 4 4 16 0 16 0 4 16 -512 +2 4 4 0 4 8 8 4 4 8 0 4 4 4 0 4 16 16 0 4 4 128 4 4 4 4 8 0 64 4 128 16 64 0 8 4 4 4 32 4 4 0 4 8 0 0 0 4 0 128 4 4 4 -512 +2 64 4 4 4 64 4 4 8 64 4 32 4 0 4 32 0 4 4 4 4 4 4 0 0 64 8 0 -512 +2 4 4 4 4 0 0 4 64 16 16 4 16 0 8 4 0 4 8 16 8 32 8 4 0 4 16 4 16 4 0 4 16 64 0 16 4 16 0 4 4 32 8 0 4 128 16 4 0 4 4 4 4 16 4 32 128 16 32 0 4 0 4 0 4 8 -512 +2 4 16 16 -512 +2 64 4 8 4 4 4 4 0 0 8 4 0 0 -512 +2 16 4 0 8 4 0 8 16 0 0 16 4 0 0 32 0 8 0 0 4 64 4 8 64 4 0 4 64 4 8 -512 +2 0 8 4 4 16 0 4 64 4 4 4 -512 +2 16 16 0 0 4 -512 +2 0 4 8 4 4 4 0 4 4 0 0 4 32 4 4 4 16 64 0 4 32 4 4 4 8 8 4 128 4 8 4 64 4 4 16 4 4 4 128 -512 +2 8 0 0 -512 +2 0 0 0 -512 +2 4 0 4 64 0 -512 +2 4 4 0 16 4 0 16 8 0 0 4 16 0 4 4 0 4 8 0 32 0 0 0 16 0 4 8 0 0 4 8 4 0 -512 +2 8 8 0 0 16 8 0 16 0 32 4 16 0 4 128 4 4 0 0 4 128 128 4 4 16 128 16 4 8 4 0 0 -512 +2 4 64 4 0 32 4 4 0 4 8 4 32 4 16 128 4 4 0 4 0 4 16 0 4 4 4 4 16 16 8 16 8 4 16 0 4 0 0 0 8 4 16 4 4 0 0 4 4 4 0 4 4 8 64 4 8 4 0 4 4 64 4 4 0 4 16 4 4 16 4 16 4 4 0 4 0 4 4 0 32 4 128 0 16 0 0 4 0 8 64 4 0 4 4 4 4 4 4 4 64 8 4 4 4 4 16 4 4 16 16 4 128 0 4 8 32 4 4 8 0 0 8 8 128 0 16 4 8 4 4 4 16 4 0 4 128 0 4 32 4 4 4 4 4 4 0 8 8 -512 +2 4 4 16 0 4 4 4 16 4 8 16 0 8 128 4 4 4 4 16 8 4 4 4 0 4 0 0 0 16 8 16 8 0 4 64 128 16 4 32 0 4 0 32 0 4 128 4 4 4 4 32 4 64 4 8 64 4 8 0 4 -512 +2 0 4 4 4 0 4 8 0 16 4 -512 +2 16 0 4 0 4 0 4 128 0 0 4 64 16 4 4 4 64 8 8 8 64 0 4 4 4 0 32 4 4 4 4 0 4 16 -512 +2 4 4 8 0 16 8 4 8 8 4 4 -512 +2 64 16 4 4 4 16 4 4 8 0 16 32 4 4 4 4 4 4 4 0 4 16 8 8 4 0 -512 +2 4 4 8 16 16 4 4 16 0 128 0 4 128 32 4 0 0 128 16 0 4 4 8 8 0 0 0 4 4 4 4 0 -512 +2 32 4 8 -512 +2 4 4 0 4 4 4 0 0 16 4 -512 +2 64 4 16 0 0 4 4 0 16 64 16 16 128 8 -512 +2 0 64 4 128 4 4 16 8 4 4 0 4 16 4 128 4 128 64 0 4 128 8 16 32 16 16 4 4 32 4 4 0 4 128 0 64 16 128 0 4 0 0 4 0 128 128 8 0 8 4 8 8 4 4 128 0 16 4 0 8 0 8 4 4 4 0 4 8 4 8 4 16 32 8 32 8 32 4 128 0 16 0 4 4 32 16 64 0 128 0 8 16 8 4 -512 +2 4 16 16 -512 +2 128 4 4 -512 +2 16 0 128 0 4 4 128 4 0 64 4 4 0 8 0 32 4 4 0 4 8 4 4 0 4 -512 +2 4 8 -512 +2 0 0 16 -512 +2 4 0 0 8 128 4 16 4 0 4 8 0 4 0 32 0 4 8 64 16 32 32 8 0 4 4 4 4 16 0 4 0 64 4 32 4 0 4 -512 +2 4 -512 +2 16 16 4 0 16 32 32 0 4 0 128 4 4 8 4 0 0 0 0 4 4 16 4 0 16 0 0 16 8 4 128 8 0 0 -512 +2 4 4 0 4 4 4 8 8 8 4 4 4 0 32 0 4 4 32 0 32 0 0 0 128 4 8 0 8 4 128 4 4 0 8 0 0 16 128 8 64 4 4 8 0 4 16 8 4 128 32 128 4 4 32 0 4 4 4 16 8 4 4 16 16 16 -512 +2 4 0 4 32 4 4 4 4 4 128 0 4 0 8 16 32 16 4 8 64 8 0 4 16 4 4 4 0 -512 +2 0 0 16 0 4 -512 +2 4 4 4 4 4 64 8 0 128 4 4 4 4 4 32 4 0 4 16 4 4 0 4 0 4 4 4 4 0 4 4 8 16 4 0 128 16 0 128 16 4 0 128 0 4 128 4 32 4 8 4 4 0 4 0 8 16 0 16 4 4 4 8 4 0 64 16 0 0 4 4 0 4 4 0 16 8 16 4 4 4 4 0 0 16 -512 +2 0 0 128 0 0 8 0 0 8 0 32 0 0 16 4 8 -512 +2 0 0 4 0 -512 +2 0 32 0 16 128 0 4 4 -512 +2 4 16 0 16 4 0 0 128 -512 +2 128 16 4 8 128 4 4 0 0 4 8 0 0 4 16 64 4 16 0 0 8 8 4 128 4 32 4 0 4 32 128 8 16 4 64 4 4 0 0 4 4 16 0 4 4 16 4 16 128 4 0 4 4 32 4 0 4 4 4 32 -512 +2 4 8 4 64 8 8 16 16 0 64 -512 +2 4 8 8 4 16 0 4 4 0 0 32 0 4 4 -512 +2 16 16 16 128 -512 +2 0 0 0 8 4 128 4 16 4 4 4 -512 +2 8 4 4 4 128 4 4 0 0 128 0 16 128 8 0 8 0 4 4 -512 +2 16 8 4 4 0 128 4 0 32 -512 +2 16 4 4 4 8 8 0 4 4 0 8 0 16 128 0 8 4 4 4 64 -512 +2 4 4 4 4 8 0 128 0 4 4 128 4 32 0 128 4 0 0 0 4 4 4 4 4 0 8 4 0 0 0 0 16 4 4 4 0 4 16 4 4 16 4 4 0 0 4 16 4 4 0 16 4 4 16 4 4 4 4 0 4 0 16 32 4 32 8 8 0 32 0 4 4 16 64 0 0 4 8 128 4 16 4 4 0 16 32 0 0 0 16 0 4 0 16 8 4 128 4 4 4 8 0 4 8 0 4 0 8 128 4 128 0 0 4 -512 +2 4 16 0 0 0 64 4 8 8 0 0 0 4 16 4 4 4 16 4 4 0 0 0 -512 +2 128 4 4 4 8 4 0 0 0 0 0 0 128 8 4 16 4 4 16 4 8 16 0 64 4 8 4 16 4 4 0 -512 +2 0 8 4 0 4 64 0 128 -512 +2 0 0 -512 +2 0 16 4 0 0 4 4 0 0 16 8 16 4 0 4 4 0 8 32 0 4 4 16 4 4 4 4 4 16 0 4 16 8 4 32 16 4 4 16 0 8 8 4 8 4 0 4 8 -512 +2 0 32 4 0 4 0 0 4 128 8 -512 +2 0 0 0 -512 +2 4 0 4 8 32 128 64 16 4 8 4 128 4 4 16 16 8 0 4 4 16 16 16 8 0 -512 +2 64 0 16 4 8 128 0 4 0 -512 +2 8 64 8 4 -512 +2 0 8 16 16 4 8 0 128 8 8 0 16 -512 +2 4 4 64 64 64 16 4 0 4 4 16 16 4 0 0 4 8 16 8 16 8 4 16 4 0 0 0 4 0 4 4 4 0 4 8 0 0 0 4 32 4 4 4 0 4 16 16 4 4 8 128 4 16 4 4 4 0 8 0 0 0 4 0 -512 +2 64 4 4 0 0 0 128 0 0 4 4 4 4 16 -512 +2 32 64 0 8 4 128 0 8 16 128 0 4 0 128 0 64 4 4 0 0 8 0 4 0 0 8 16 0 64 8 0 4 4 4 4 4 128 4 0 16 128 4 8 4 0 16 4 4 8 16 4 16 4 16 4 32 8 4 0 32 8 0 4 0 4 4 4 4 4 16 128 4 0 128 16 128 16 4 32 32 4 16 4 4 4 4 8 4 8 4 0 0 4 4 32 16 64 32 4 16 4 0 4 4 16 0 0 4 0 0 32 16 8 4 4 0 16 16 4 0 -512 +2 16 16 0 4 128 4 4 0 0 0 8 4 -512 +2 8 16 4 0 8 8 64 16 8 4 4 4 4 0 0 8 8 0 16 4 16 0 4 4 0 8 4 4 128 4 8 0 64 32 64 0 4 4 4 0 8 8 0 16 4 16 0 64 4 4 4 4 0 4 4 -512 +2 64 64 4 64 0 4 0 0 4 0 4 4 16 4 0 4 0 8 16 8 16 16 0 4 32 16 16 4 4 4 16 0 0 4 16 4 16 8 8 4 4 4 4 4 0 16 8 8 8 4 16 4 4 8 16 -512 +2 4 0 16 4 0 4 32 4 0 16 4 4 0 4 16 0 4 16 0 8 16 4 4 4 4 0 16 16 4 8 4 8 4 4 4 0 4 4 0 8 4 4 -512 +2 4 128 4 0 4 16 8 32 0 64 4 0 4 64 4 4 8 0 4 0 0 8 0 4 0 4 8 0 0 4 32 16 0 4 64 16 16 0 4 4 16 128 32 4 16 4 16 64 16 4 4 32 8 4 0 0 4 0 0 8 4 0 4 8 4 4 4 4 4 0 0 0 16 -512 +2 16 0 4 4 0 0 8 64 0 16 4 0 0 0 128 0 0 4 16 8 16 4 128 4 diff --git a/t/classify_ref_Z4j b/t/classify_ref_Z4j index ec65988..5653d40 100644 --- a/t/classify_ref_Z4j +++ b/t/classify_ref_Z4j @@ -1,4000 +1,4000 @@ 4 0 0 4 0 4 0 16 0 16 0 16 0 0 0 0 0 16 0 4 8 8 8 0 4 0 0 0 16 0 8 0 4 0 4 0 0 0 0 0 0 4 0 0 0 4 4 0 0 0 4 8 0 0 4 0 0 16 0 0 0 4 16 0 0 4 0 0 16 0 4 0 0 0 4 8 0 16 0 8 0 0 8 4 0 0 0 0 0 4 0 0 0 4 0 0 0 0 0 0 0 0 16 4 0 0 0 0 0 -512 +2 0 0 0 0 16 8 4 0 0 0 16 0 0 8 0 0 4 0 0 0 0 0 8 0 0 0 0 8 0 0 16 0 4 4 16 4 0 0 0 0 0 8 0 0 0 0 0 16 4 4 0 0 0 4 0 0 0 0 0 4 0 4 16 0 0 16 0 0 0 0 4 0 0 4 0 16 0 8 16 16 0 0 0 0 8 0 0 0 16 0 8 0 8 0 0 4 0 0 0 0 0 4 0 0 16 4 16 4 4 0 0 0 0 8 0 0 16 0 0 0 8 0 0 16 0 16 0 8 0 0 0 0 0 4 8 0 0 0 0 -512 +2 0 0 0 4 0 8 16 16 0 -512 +2 8 0 0 0 0 0 16 0 0 0 0 4 4 0 0 0 0 0 0 0 8 8 4 8 0 0 16 0 4 0 0 0 0 0 0 4 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 4 0 0 0 0 4 16 4 0 0 0 0 8 0 0 0 16 0 8 0 8 0 4 -512 +2 16 0 0 4 4 0 8 4 4 0 0 0 0 4 0 0 0 0 0 0 0 16 0 0 16 0 8 4 4 8 16 16 0 0 0 8 16 0 0 8 4 0 16 0 8 0 0 16 0 0 16 0 0 0 4 4 0 4 4 0 0 4 0 0 4 0 0 4 4 0 0 16 0 0 0 4 16 0 16 0 0 0 4 16 8 8 0 0 0 0 0 0 4 0 4 8 16 0 4 0 8 0 0 16 4 0 0 0 0 0 8 0 0 8 0 16 4 0 0 16 0 8 0 0 0 0 0 0 0 0 0 0 8 0 16 0 0 0 0 0 0 0 0 4 8 16 0 0 8 8 16 16 0 0 4 4 0 0 0 0 0 0 4 0 0 0 0 0 16 4 0 8 16 0 0 0 0 0 0 8 0 0 4 4 0 0 4 4 0 0 0 0 4 0 4 0 0 16 0 0 8 0 0 0 16 0 4 8 0 0 0 16 16 0 8 0 0 0 8 0 4 4 0 0 0 0 0 0 0 0 4 0 0 0 0 0 16 4 16 0 0 8 8 0 4 0 4 0 4 0 8 0 0 8 4 0 0 8 0 0 0 0 16 4 0 0 8 0 16 0 0 16 8 0 8 0 0 0 0 4 0 0 0 4 0 0 0 0 0 0 8 4 0 4 0 0 0 0 0 0 4 0 0 0 16 4 16 0 0 0 0 0 0 0 8 0 8 0 16 0 8 4 0 0 0 8 16 8 0 4 0 4 0 0 8 0 0 4 4 0 0 4 8 0 0 0 0 0 0 0 16 0 8 0 16 16 8 0 4 0 0 4 0 8 0 0 4 0 0 8 0 0 4 0 4 4 8 16 4 0 0 0 0 0 0 0 4 0 0 0 0 0 4 4 0 0 0 0 16 0 0 8 4 4 4 0 0 0 0 4 0 0 0 0 8 0 8 4 4 0 0 0 0 16 0 4 0 16 0 0 4 0 4 0 0 0 0 0 0 0 0 0 0 0 0 8 0 0 0 4 0 0 8 0 0 0 0 0 0 4 0 0 0 16 0 0 0 0 4 4 8 0 0 0 0 0 0 16 4 0 16 0 0 0 4 0 8 0 0 16 4 4 0 16 0 8 0 16 0 0 8 0 4 0 8 16 0 0 -512 +2 0 4 4 0 0 0 16 0 0 0 0 8 0 0 0 0 0 0 -512 +2 4 8 4 0 0 16 0 0 0 8 0 0 0 4 0 0 0 0 0 16 0 0 0 0 0 4 4 16 0 4 16 0 0 0 8 0 0 0 4 0 0 0 16 16 4 4 -512 +2 0 4 0 0 0 0 4 8 0 4 8 16 4 4 8 0 0 0 0 16 16 8 0 0 0 0 16 4 0 0 4 4 16 4 0 0 0 0 0 0 0 8 8 0 4 16 0 8 16 4 8 0 0 0 4 0 4 0 0 4 0 0 0 4 8 0 8 0 16 0 8 0 0 0 0 0 0 0 0 0 0 0 8 0 4 0 0 0 0 0 16 16 0 0 0 0 0 0 16 4 4 0 0 16 4 0 16 0 4 16 0 -512 +2 0 0 4 8 0 0 4 8 4 8 4 0 0 0 8 8 0 0 0 16 4 0 0 0 0 0 8 0 0 16 0 0 0 0 4 0 8 16 4 0 4 0 8 0 8 16 0 0 0 0 16 8 4 0 0 8 0 16 4 0 4 0 0 0 0 0 0 0 0 0 0 16 8 0 16 0 0 4 0 8 0 0 16 0 4 0 8 0 0 0 0 16 4 0 0 -512 +2 8 8 0 0 0 0 4 0 0 4 0 0 0 4 16 0 0 4 0 8 16 0 0 0 4 0 0 0 16 0 16 0 0 0 0 0 0 0 0 0 8 0 0 0 0 8 0 4 0 0 16 16 0 0 0 4 4 4 0 0 0 0 0 16 0 0 8 0 16 0 0 4 0 0 4 0 0 0 0 0 16 0 0 0 4 16 8 16 16 0 0 0 16 16 16 0 0 0 0 0 0 0 0 8 4 0 0 0 4 4 0 16 0 0 0 0 0 4 0 0 0 0 0 0 4 4 0 0 16 16 0 0 4 0 0 4 4 0 0 0 0 16 8 4 16 0 0 4 0 0 4 0 0 8 8 0 16 0 0 0 4 4 16 0 0 4 0 0 0 8 0 0 0 0 0 0 8 0 0 0 0 4 0 8 0 0 4 0 8 0 0 0 16 4 4 8 0 8 0 0 8 0 0 16 0 0 0 16 0 0 0 0 8 0 8 16 4 0 0 8 4 0 8 0 4 0 16 0 4 0 8 8 0 8 4 0 8 0 4 0 0 0 0 0 4 0 0 0 4 0 16 0 4 0 0 4 0 0 0 0 0 0 16 0 0 0 0 4 -512 +2 0 0 8 0 8 0 0 0 0 0 8 0 4 0 4 16 0 0 8 0 0 0 4 0 0 0 16 0 0 8 0 4 0 0 16 0 4 0 0 0 0 0 4 0 8 0 0 0 8 -512 +2 0 0 0 0 4 0 0 0 8 0 0 0 0 0 4 16 4 0 0 16 0 4 0 16 0 0 0 16 0 0 0 0 0 0 8 0 0 0 0 0 4 0 0 0 0 16 0 0 8 0 0 0 16 16 0 0 0 0 0 4 0 4 0 0 8 0 0 4 0 16 0 0 0 0 0 0 8 0 0 4 0 0 16 0 0 0 0 0 16 0 0 0 4 0 0 0 16 0 0 0 0 0 0 0 0 0 0 0 4 0 4 0 0 0 0 0 4 0 16 16 0 8 0 0 4 0 0 8 4 0 0 0 0 0 0 0 0 -512 +2 16 8 0 0 4 0 0 4 0 0 8 8 0 0 0 4 16 16 8 0 8 4 4 0 16 0 16 0 16 0 0 8 0 16 0 8 0 0 0 0 8 4 0 16 0 0 0 8 0 0 4 0 0 4 4 0 16 16 0 8 16 0 0 0 0 16 8 8 0 4 4 4 0 0 0 0 0 -512 +2 0 0 0 0 4 0 16 0 8 0 0 0 4 0 0 4 0 0 4 4 0 16 0 16 0 0 0 0 8 0 0 16 8 0 0 0 0 4 0 8 4 0 0 0 0 4 0 0 0 0 0 8 0 8 16 4 4 4 0 0 4 0 0 0 8 0 0 0 4 -512 +2 4 0 4 0 0 0 0 0 8 0 0 0 0 4 0 0 0 4 0 0 8 0 0 0 4 8 0 0 0 0 4 0 0 16 0 0 -512 +2 0 4 0 4 0 16 0 8 0 4 0 0 0 0 0 0 0 0 4 8 0 0 0 0 4 0 0 4 16 0 8 0 0 16 8 0 0 0 0 0 4 16 0 8 16 0 8 4 -512 +2 4 0 0 16 8 16 8 4 4 0 0 0 0 0 0 0 0 8 0 0 0 0 8 4 0 0 0 0 0 4 0 0 0 0 0 0 0 0 16 0 0 4 8 0 8 4 16 0 16 16 4 0 4 8 4 0 0 0 0 0 4 0 0 0 0 0 8 0 4 0 0 0 0 4 0 8 0 8 16 4 16 0 0 16 4 0 0 16 8 0 0 0 0 8 4 16 0 0 0 0 0 16 0 -512 +2 4 0 0 0 4 0 4 4 4 0 0 0 0 0 0 0 0 0 8 0 4 0 0 0 0 4 0 0 4 0 4 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 4 0 0 16 4 0 4 0 4 4 0 4 0 0 0 0 0 0 0 16 0 0 0 0 0 0 4 8 0 0 0 0 0 8 4 16 0 0 0 0 0 0 4 0 0 0 0 0 0 4 0 4 0 0 4 0 0 0 4 0 0 0 0 0 0 0 0 0 0 -512 +2 0 16 0 0 0 0 0 4 8 4 0 16 0 0 0 0 4 0 0 8 0 0 8 0 0 4 0 -512 +2 0 0 16 4 8 0 4 8 0 8 0 0 4 0 0 16 0 0 0 0 8 0 8 0 4 0 0 0 0 8 0 0 16 4 8 0 0 0 0 0 16 0 0 0 16 0 0 16 0 16 16 4 16 0 4 0 0 16 16 16 16 0 0 0 16 0 0 0 0 4 4 0 0 4 0 16 0 16 0 4 0 0 0 0 4 16 16 4 -512 +2 0 0 0 16 0 0 16 8 0 8 0 0 0 0 0 8 0 0 0 16 4 4 4 0 0 0 0 0 0 0 4 0 0 8 0 16 16 0 0 4 0 4 0 4 0 0 4 8 0 4 0 0 4 0 0 0 0 0 0 0 4 4 4 4 0 0 0 0 0 -512 +2 4 8 0 0 4 8 0 0 0 0 0 0 0 0 0 0 0 0 16 4 8 4 8 0 0 0 0 4 0 0 0 4 0 0 0 16 8 4 0 0 0 0 0 0 0 0 0 0 16 16 16 0 0 4 0 0 8 0 8 0 0 4 16 0 0 8 0 0 4 0 0 8 8 0 0 8 4 0 0 0 8 0 0 0 4 0 0 0 8 -512 +2 0 0 8 0 16 0 0 4 4 0 0 0 0 0 0 0 4 0 -512 +2 0 0 0 0 0 4 8 0 4 0 0 16 4 0 0 0 0 4 0 16 4 0 0 0 4 0 0 0 4 0 0 4 16 4 4 0 0 0 0 16 16 0 0 0 0 0 0 0 4 0 0 8 0 0 0 0 16 0 0 0 4 0 8 0 0 0 16 16 0 0 16 0 0 0 16 0 0 8 0 8 0 0 8 0 0 0 0 16 0 0 0 8 0 4 0 4 8 16 4 0 0 0 0 8 8 0 0 0 0 0 0 8 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 -512 +2 8 0 0 8 0 0 8 0 4 0 0 0 0 0 8 0 0 0 0 0 0 0 4 0 0 0 8 -512 +2 4 0 8 0 0 0 0 0 0 8 0 0 0 0 0 16 16 0 0 0 0 0 16 0 0 0 0 4 16 8 0 0 0 0 0 0 0 8 4 0 0 0 0 16 0 0 4 0 0 4 0 0 8 0 0 0 16 0 0 4 16 16 8 4 0 0 0 0 4 8 0 16 4 0 0 0 4 0 0 4 0 0 0 0 0 0 4 4 8 0 0 0 0 0 0 0 4 0 4 0 0 0 0 0 16 4 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 16 0 0 0 0 0 0 8 0 4 0 0 0 4 0 0 4 0 0 0 0 0 4 16 0 4 0 0 0 0 8 0 0 4 0 4 0 0 8 8 0 0 0 16 0 4 0 0 8 0 16 0 0 16 8 0 0 0 0 0 0 8 0 4 0 0 0 0 0 8 0 0 4 0 0 0 0 0 0 0 0 8 0 0 8 0 0 4 4 0 16 4 0 0 0 0 16 0 4 0 0 16 0 -512 +2 0 0 0 0 0 0 0 0 0 0 8 8 4 4 8 0 0 4 -512 +2 0 8 4 0 4 0 4 -512 +2 0 4 16 0 0 -512 +2 0 4 0 8 0 0 16 0 0 4 0 8 8 4 8 0 16 0 0 8 8 4 0 0 4 4 0 8 8 0 0 0 4 0 16 0 0 0 0 0 16 0 4 -512 +2 0 0 0 4 0 4 0 0 0 0 4 8 0 0 0 0 0 8 16 16 4 0 0 8 0 8 0 0 0 16 0 0 0 4 0 4 0 16 0 0 0 0 0 0 4 4 4 4 -512 +2 0 8 0 0 4 0 0 0 0 0 16 0 0 0 0 4 0 0 8 0 0 0 0 -512 +2 0 8 0 0 0 8 0 0 8 0 0 0 0 0 0 8 0 0 0 0 4 4 0 4 4 16 0 0 0 0 8 16 4 0 0 0 0 0 0 0 0 -512 +2 0 0 0 4 0 0 0 0 16 4 0 0 8 0 8 0 0 8 4 0 0 0 0 0 16 8 0 0 0 0 0 0 0 0 8 0 0 0 0 0 0 0 4 0 0 0 0 8 0 8 0 0 0 4 16 0 4 0 0 0 16 4 0 8 4 16 4 4 0 0 4 0 0 0 0 4 0 0 16 0 0 0 0 8 8 16 0 8 4 0 16 0 0 0 4 0 4 0 0 0 0 0 16 16 0 0 0 0 0 0 0 0 8 4 4 0 0 4 0 0 4 0 8 0 0 0 0 0 0 0 0 0 0 0 4 8 0 0 0 4 0 4 0 0 16 0 0 0 0 0 8 4 4 0 0 0 0 4 16 8 4 4 4 0 0 0 16 8 0 4 0 0 0 0 4 4 0 0 4 0 0 0 8 16 -512 +2 4 0 4 0 0 4 0 8 0 0 4 -512 +2 4 4 16 0 0 4 0 0 0 0 0 0 4 4 0 0 0 0 0 4 0 4 0 16 4 0 0 0 0 16 16 8 0 0 0 0 0 0 0 16 4 16 4 8 8 0 0 0 0 16 8 4 0 4 0 0 16 4 0 8 4 0 0 0 0 0 0 0 4 0 0 0 0 4 0 0 4 0 -512 +2 0 0 4 0 0 0 0 0 0 0 0 0 0 0 16 0 4 4 0 0 0 0 0 4 4 8 0 16 0 0 0 16 0 16 0 0 0 0 0 8 8 0 0 0 0 16 0 16 0 0 0 4 4 4 4 8 0 0 0 0 0 4 0 4 8 4 0 8 0 16 16 0 0 0 0 8 0 4 16 0 8 16 0 16 4 0 0 0 0 0 4 0 4 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 16 0 0 0 0 0 0 0 0 0 4 0 8 8 0 4 0 0 0 4 0 0 0 0 0 4 4 0 0 4 0 4 0 0 0 0 0 0 0 4 8 0 0 0 0 0 4 8 -512 +2 0 0 0 0 4 0 0 0 0 0 8 0 0 0 0 0 0 0 8 4 4 0 16 0 0 4 0 0 4 0 0 0 0 0 0 0 0 0 8 0 0 8 0 8 8 0 0 8 0 4 0 0 4 0 16 16 0 16 0 0 0 0 0 0 0 4 0 4 0 4 0 0 16 0 0 4 0 4 0 4 0 16 8 16 0 0 16 0 0 0 0 0 8 0 8 0 0 8 8 0 0 0 0 4 16 4 4 0 8 0 8 0 4 0 16 0 8 16 0 0 4 4 4 4 0 0 4 0 4 0 0 0 0 16 0 0 0 16 8 0 16 0 0 0 0 0 0 0 -512 +2 4 -512 +2 0 16 0 8 8 8 0 0 0 0 0 0 0 0 16 4 16 4 0 0 0 0 0 8 8 0 0 0 4 0 0 0 0 0 0 0 4 0 0 0 16 0 0 0 4 0 0 8 8 0 16 0 0 0 0 8 8 0 0 4 0 16 0 16 4 -512 +2 0 0 0 16 4 0 0 4 0 0 0 0 0 0 8 0 0 16 0 0 0 0 0 0 0 0 0 0 -512 +2 0 0 8 0 0 0 0 0 0 4 0 0 0 0 0 0 0 8 4 8 0 4 0 0 4 0 0 8 0 0 16 0 4 0 0 0 0 0 16 0 0 0 8 4 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -512 +2 0 0 0 0 0 0 0 8 0 16 0 0 0 16 16 0 0 16 0 8 0 0 4 0 0 0 0 -512 +2 0 0 8 0 16 16 0 8 4 16 8 0 0 0 0 0 0 0 0 4 0 0 16 0 4 0 0 8 0 4 0 8 0 4 16 8 0 16 8 4 0 8 16 0 0 8 4 0 -512 +2 0 0 0 0 0 0 0 0 4 0 0 0 8 0 4 4 0 0 4 0 0 16 16 0 -512 +2 0 0 0 0 4 0 0 4 0 16 16 0 0 0 8 0 0 16 0 0 0 0 0 0 0 16 16 0 0 0 4 4 0 8 0 0 0 0 8 0 16 0 0 8 0 4 4 0 8 16 4 0 0 0 0 0 16 0 4 0 0 0 0 0 0 4 0 0 -512 +2 0 4 0 4 0 0 0 0 0 0 0 0 16 0 0 0 16 0 8 0 0 0 0 0 0 0 4 4 16 8 4 0 0 0 8 0 0 8 0 0 8 0 16 0 4 8 0 0 4 8 0 0 4 8 0 0 4 0 0 4 0 16 0 0 0 16 0 0 8 8 4 0 0 0 0 0 0 -512 +2 4 0 -512 +2 0 0 0 8 0 0 8 4 0 0 4 16 0 0 0 0 0 0 4 4 0 0 4 8 4 0 16 8 4 0 0 0 8 0 0 4 0 0 0 0 0 4 4 0 0 0 0 0 4 8 0 8 8 0 0 0 0 0 16 0 0 16 0 0 4 0 4 8 8 0 0 0 16 0 16 0 0 16 4 4 16 0 0 0 16 8 0 0 0 4 0 4 0 4 0 4 0