#include <iostream>
#include <iomanip>

#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

#include "JOscProb/JOscParameters.hh"


/**
 * \file
 * Auxiliary program to test JOscParameters.
 *
 * \author bjung
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;

  JOscParameters          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);

  JOscParameters parameters1("sinsqTh12", 0.304,
			     "deltaCP",   1.544,
			     "dM21sq",    7.42e-5);

  NOTICE("parameters1:\n" << parameters1);

  ASSERT(parameters1.is_valid() && parameters1.size() == 3);

  parameters1.set("sinsqTh23", -0.5);

  NOTICE("parameters1.sinsqTh23 = -0.5 -> parameters1:\n" << parameters1);

  ASSERT(!parameters1.is_valid() && parameters1.size() == 4);

  JOscParameters parameters2("dM31sq",    -2.49e-3,
			     "sinsqTh13",  0.02241,
			     "sinsqTh23",  0.570 );

  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));
}