#include #include #include "Jeep/JPrint.hh" #include "Jeep/JParser.hh" #include "Jeep/JMessage.hh" #include "JOscProb/JOscParametersGrid.hh" /** * \file * Auxiliary program to test JOscParametersGrid. * * \author bjung */ int main(int argc, char **argv) { using namespace std; using namespace JPP; JOscParametersGrid parameters; int debug; try { JParser<> zap; zap['@'] = make_field(parameters) = JPARSER::initialised(); zap['d'] = make_field(debug) = 3; zap(argc, argv); } catch (const exception& error) { FATAL(error.what() << endl); } NOTICE("User-defined parameters:\n" << parameters); ASSERT(parameters.size() > 0); JOscParametersGrid parameters1("sinsqTh12", make_grid(5, 0.3, 0.4), "deltaCP", make_grid(6, 1.5, 1.6), "dM21sq", make_grid(7, 7.4e-5, 7.5e-5)); NOTICE("parameters1:\n" << parameters1); ASSERT(parameters1.is_valid() && parameters1.size() == 3); parameters1.set("sinsqTh23", make_grid(3, -0.6, -0.4)); NOTICE("parameters1.sinsqTh23 = -0.5 -> parameters1:\n" << parameters1); ASSERT(!parameters1.is_valid() && parameters1.size() == 4); JOscParametersGrid parameters2("dM31sq", make_grid(5, -2.6e-3, -2.0e-3), "sinsqTh13", make_grid(2, 0.01, 0.03 ), "sinsqTh23", make_grid(3, 0.57, 0.59 )); NOTICE("parameters2:\n" << parameters2); ASSERT(parameters2.is_valid() && parameters2.size() == 3); ASSERT(!parameters1.contains(parameters2)); NOTICE("Joining parameters2 with parameters1" << endl); parameters1.join(parameters2); ASSERT(parameters1.is_valid() && parameters1.size() == 6); ASSERT(parameters1.contains(parameters2)); }