diff.cxx

Go to the documentation of this file.

//Dear emacs, this is -*- c++ -*-
 
/**
 * @file diff.cxx
 * @author <a href="mailto:[email protected]>Andr� Rabello dos
 * ANJOS</a>
 * $Author: zhangy $
 * $Revision: 1.1.1.1 $
 * $Date: 2009/06/19 07:35:41 $
 *
 * This source code implements a comparison application that can be used to
 * show differences in event-formatted files, at the ROB level
 */
 
#include "eformat/eformat.h"
#include <algorithm>
 
/**
 * The maximum event size, in words
 */
const size_t MAX_EVENT_SIZE = 2500000;
 
/**
 * The maximum number of ROBS
 */
const size_t MAX_ROBS = 2048;
 
/**
 * Compares two ROB fragments for equality
 *
 * @param p1 The first ROB fragment
 * @param p2 The second ROB fragment
 */
int cmp_source_id (const uint32_t* p1, const uint32_t* p2)
{
  eformat::ROBFragment<const uint32_t*> rob1(p1);
  eformat::ROBFragment<const uint32_t*> rob2(p2);
  return rob1.rod_source_id() < rob2.rod_source_id();
}
 
/**
 * Shows differences between to event-formatted files at the ROB level.
 */
int main (int argc, char** argv)
{
  using namespace eformat;
  
  if ( argc != 3 ) {
    std::cerr << "usage: " << argv[0] << " <file1> <file2>"
          << std::endl;
    std::cerr << "OBS: The event order should be the same on both files" 
          << std::endl;
    std::exit(1);
  }
 
  //open normally a file
  size_t event_counter = 0;
  std::fstream in1(argv[1], std::ios::in|std::ios::binary);
  std::fstream in2(argv[2], std::ios::in|std::ios::binary);
  const uint32_t* robp1[MAX_ROBS];
  const uint32_t* robp2[MAX_ROBS];
 
  uint32_t* event1 = new uint32_t[MAX_EVENT_SIZE];
  uint32_t* event2 = new uint32_t[MAX_EVENT_SIZE];
 
  while (true) {
    if (!(next_fragment(in1, event1, MAX_EVENT_SIZE*4)) || 
    !(next_fragment(in2, event2, MAX_EVENT_SIZE*4))) break;
    ++event_counter;
    try {
      FullEventFragment<const uint32_t*> fe1(event1);
      FullEventFragment<const uint32_t*> fe2(event2);
      fe1.check_tree(); // no checks, so we dump all
      fe2.check_tree(); // no checks, so we dump all
      size_t s1 = get_robs(event1, robp1, MAX_ROBS);
      size_t s2 = get_robs(event2, robp2, MAX_ROBS);
      if (s1 != s2) {
    std::cerr << "The number of ROB's of event " << event_counter
          << " is different in file `" << argv[1]
          << "' (" << s1 << ") and in file `" << argv[2]
          << "' (" << s2 << "). Skipping comparison." << std::endl;
    continue;
      }
      //if we get here, the number of ROD's is the same
      std::sort(robp1, robp1+s1, cmp_source_id);
      std::sort(robp2, robp2+s2, cmp_source_id);
      bool mark = false;
      for (size_t i=0; i<s1; ++i) {
    eformat::ROBFragment<const uint32_t*> rob1(robp1[i]);
    eformat::ROBFragment<const uint32_t*> rob2(robp2[i]);
    if (rob1.rod_source_id() != rob2.rod_source_id()) {
      std::cerr << "! ROB[" << i << "]" << std::endl
            << "- " << HEX(rob1.rod_source_id()) << std::endl
            << "+ " << HEX(rob2.rod_source_id()) << std::endl;
      mark = true;
    }
      }
      if (mark) {
    std::cerr << "Event #" << fe1.lvl1_id() << " differ." << std::endl;
      }
    }
    catch (eformat::Issue& ex) {
      std::cerr << std::endl
        << "Uncaught eformat issue: " << ex.what() << std::endl;
      std::cout << "Trying to continue..." << std::endl;
      continue;
    }
    catch (ers::Issue& ex) {
      std::cerr << std::endl
        << "Uncaught ERS issue: " << ex.what() << std::endl;
      delete[] event1;
      delete[] event2;
      std::exit(1);
    }
    catch (std::exception& ex) {
      std::cerr << std::endl
        << "Uncaught std exception: " << ex.what() << std::endl;
      delete[] event1;
      delete[] event2;
      std::exit(1);
    }
    catch (...) {
      std::cerr << std::endl << "Uncaught unknown exception" << std::endl;
      delete[] event1;
      delete[] event2;
      std::exit(1);
    }
 
  }
 
  std::cerr << "In the absence of remarks, files do NOT differ." << std::endl;
 
  delete[] event1;
  delete[] event2;
  return 0;
}