/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/lcio/addDifSynchro.cc

Go to the documentation of this file.

#include "lcio/AddDifSynchro.hh" 
#include "lcio/Count.hh" 


#include "Log.hh" 

#include <iostream>
#include <string>

using namespace std;

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

        if ( argc != 4  ) {
                cout  << endl <<"usage:addDifSynchro difSynchro_file slab_file output_file " << endl;
                exit(1);
        }
        // assign file name
        string difSynchroName;
        difSynchroName.assign(argv[1]);

        string slabFileName;
        slabFileName.assign(argv[2]);

        unsigned int synchroNumber = Count(difSynchroName).getNumberOfEvent();
        unsigned int slabNumber = Count(slabFileName).getNumberOfEvent();
        
    if ( synchroNumber != slabNumber ) 
        {
                cout << endl << " !!! ERROR !!! " << endl << endl ;  
                cout << "Number of event differ in both file" << endl ;  
                cout << difSynchroName << " "  << synchroNumber << " number of events." << endl;
                cout << slabFileName << " "  << slabNumber << " number of events." << endl;
                exit ( -1 ) ;
        }
        cout << "--- " << synchroNumber << " event to merge ---" <<endl; 
/*
        // Verify events number 

        Count *countSync = new Count(difSynchroName);
    LCReader* lcReader = LCFactory::getInstance()->createLCReader() ;
    lcReader->open(slabFileName);
    // registere processor to LCReader
    lcReader->registerLCEventListener( evtProc ) ;

    // Start reading slabfile
    lcReader->readStream() ;

        Count *countSlab = new Count(slabFileName);
*/      

        // Create AddDifSynchro processor with difSynchro_file and output file as parameters
        AddDifSynchro *evtProc = new AddDifSynchro(difSynchroName,argv[3]) ;
         
        // Open slabFile to get event information to be coy in new output file
        LCReader* lcReader = LCFactory::getInstance()->createLCReader() ;
        lcReader->open(slabFileName);
        // registere processor to LCReader
        lcReader->registerLCEventListener( evtProc ) ;
        
        // Start reading slabfile
        lcReader->readStream() ;
        
        lcReader->close() ;
        delete lcReader ;
        delete evtProc;
        return 0 ;
}

---