/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/eformat/eformat-00-00-04/test/diff.cxx

Go to the documentation of this file.

//Dear emacs, this is -*- c++ -*-

#include "eformat/eformat.h"
#include <algorithm>

const size_t MAX_EVENT_SIZE = 2500000;

const size_t MAX_ROBS = 2048;

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

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

---