#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include "Quaero_utils/QuaeroEvent.hh"
#include "Quaero_utils/QuaeroItemFile.hh"
using namespace std;


int main(int argc, char* argv[])
{
  if(argc==1)
    {
      cout << "Usage:  calculateVariable <files>" << endl;
      exit(1);
    }
  assert(argc>=1);
  QuaeroEvent e;
  double w=0;
  for(int i=1; i<argc; i++)
    {
      QuaeroItemFile<QuaeroEvent> ef(argv[i]); // open Quaero file
      while(ef.nextEvent(e)) // read next Quaero event e
	{
	  HepLorentzVector l = e.getThisObject("e+")->getFourVector();
	  HepLorentzVector v = e.getPmiss();
	  HepLorentzVector j = e.getThisObject("j")->getFourVector();
	  HepLorentzVector b = e.getThisObject("b")->getFourVector();
	  QuaeroRecoObject::ChiSqdConstrainNeutrino(l,v);
	  HepLorentzVector w = l+v;
	  cout << (w+j+b).m() << " " << (j+b).m() << endl;
	} // loop over all events in file
    } // loop over all files
  return(0);
}