diff.cxx

Go to the documentation of this file.

// Dear emacs, this is -*- c++ -*-
 
/**
 * @file diff.cxx
 * @author <a href="mailto:Andre.dos.Anjos@cern.ch>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;
}