JPrintAAnet.cc File Reference

Example program to print track fit results from Evt formatted data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
#include <cmath>
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/Evt.hh"
#include "JAAnet/JAAnetToolkit.hh"
#include "JReconstruction/JEvtToolkit.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JSupport.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Example program to print track fit results from Evt formatted data.

Author

mdejong

Definition in file JPrintAAnet.cc.

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 83 of file JPrintAAnet.cc.

{
  using namespace std;
  using namespace JPP;
 
  JMultipleFileScanner<Evt> inputFile;
  JLimit_t&      numberOfEvents = inputFile.getLimit();
  vector<string> keys;
  int            debug;
 
  try {
 
    JParser<> zap("Example program to print track fit results from Evt formatted data.");
 
    zap['f'] = make_field(inputFile);
    zap['n'] = make_field(numberOfEvents)      = JLimit::max();
    zap['k'] = make_field(keys,   "print optional weights: " << get_keys(getWeight))         = JPARSER::initialised();
    zap['d'] = make_field(debug)               = 2;
 
    zap(argc, argv);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
    
  while (inputFile.hasNext()) {
 
    cout << "event " << setw(10) << inputFile.getCounter() << endl;
    
    const Evt* evt = inputFile.next();
 
    if (has_reconstructed_jppmuon(*evt)) {
 
      const Trk trk = get_best_reconstructed_jppmuon(*evt);
 
      print(cout, "", trk);
      print(cout, "->", *evt, trk.mother_id);
 
      for (const auto& key : keys) {
        cout << ' ' << SCIENTIFIC(12,3) << getWeight(trk, key, 0.0);
      }
 
      cout << endl;
    }
 
    if (debug >= debug_t) {
      for (const Trk& i : evt->trks) {
        print(cout, "", i);
        cout << endl;
      }
    }
  }
}