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

Go to the documentation of this file.

#include "lcio.h"

#include "IO/LCWriter.h"
#include "EVENT/LCIO.h"
#include <EVENT/RawCalorimeterHit.h>
#include "DATA/LCFloatVec.h"
#include "DATA/LCIntVec.h"

#include "IMPL/LCEventImpl.h" 
#include "IMPL/LCRunHeaderImpl.h" 
#include "IMPL/LCCollectionVec.h"
#include "IMPL/RawCalorimeterHitImpl.h"
#include "IMPL/SimTrackerHitImpl.h"
#include "IMPL/MCParticleImpl.h" 
#include "IMPL/LCFlagImpl.h" 
#include "IMPL/LCTOOLS.h"

#include "IMPL/TrackerRawDataImpl.h"
#include "IMPL/TrackerDataImpl.h"
#include "IMPL/TrackerPulseImpl.h"

#include "UTIL/LCRelationNavigator.h"
#include "UTIL/LCTime.h"
//#include "UTIL/BitField64.h"
#include "UTIL/CellIDEncoder.h"
#include "UTIL/CellIDDecoder.h"
#include <UTIL/BitField64.h>
#include "UTIL/LCTypedVector.h"

#include "UTIL/LCSplitWriter.h"

#include "mTypes.h"
#include "geometry/Detector.hh"
#include "event/Run.hh"
#include "geometry/Dif.hh"
#include "geometry/Chip.hh"
#include "geometry/Channel.hh"
#include "tools/SteerDesc.hh"
#include "xml/XMLTool.hh"



#include <cstdlib>
#include <iostream>
#include <sstream>


using namespace std ;
using namespace lcio ;



//struct MyTrackLink : public LCObjectLinkTraits< SimTrackerHit, MyTrackLink > {} ;

/** Simple test program to demonstrate writing of data with lcio.
 */



int main(int argc, char** argv )
{
  if( argc !=2 && argc !=3  )
  {
    cout <<  endl <<"usage lcioTest slcioFile [xmlFile]" <<endl;        
        exit(-1);
  }

  bool geometrySet = false;
 
  // get lcio file name and xmlFile file name from parameters
  static string lcioFile = argv[1] ; 

  Detector detector;
  if  ( argc == 3 )
  {
    geometrySet = true;
    string steerName;
    steerName.assign(argv[2]);

    // test if steer file exist
    FILE *file = fopen(steerName.c_str(), "r");
    if(file == NULL){
      cout << "Steer file ["<< steerName.c_str() << "] does not exist" << endl;
      fclose(file);
      exit(0);
    }
  
    // create Run, set detector and steerDesc from steer XML file
    Run run(detector);
    SteerDesc steerDesc;
    XMLTool xml;
    xml.parse(steerDesc, steerName); // parse steer file and fill steerDesc
     // Set detector with xml file content
    if ( ! steerDesc.setRun(run))
    {
      exit(1);
    }
    fclose(file);
    detector.build(steerDesc);
  }

//----------------

        IO::LCReader* lcReader = IOIMPL::LCFactory::getInstance()->createLCReader() ;
        lcReader->open( lcioFile ) ;


        LCEvent* evt = NULL;
        //----------- the event loop -----------
        while( (evt = lcReader->readNextEvent()) != 0 ) 
          {

        cout << endl << endl << endl << " NEW EVENT " << endl;
        UTIL::LCTOOLS::dumpEvent( evt ) ;
        EVENT::long64 timestamp = evt->getTimeStamp();

        LCTime time(timestamp);
        cout << time.getDateString() << endl;

        cout << "Dif synchro: "         << evt->getParameters().getIntVal("Dif_syncro" ) << endl ;
        cout << "Temperature:"          << evt->parameters().getFloatVal("Temperature" ) << endl ;
        cout << "Pressure:"             << evt->parameters().getFloatVal("Pressure")  << endl;
        cout << "globalTriggerCounter:" << evt->parameters().getIntVal("Global trigger counter")  << endl;
        cout << "difTriggerCounter:"    << evt->parameters().getIntVal("Dif trigger counter" )    << endl;
            
            const vector<string>* collNames = evt->getCollectionNames ();
        
            vector<string>::const_iterator it;
        
            for ( it = collNames->begin(); it!=collNames->end(); ++it) 
              {
                cout << *it << endl;
                LCCollection*  collec = evt->getCollection(*it);
                cout << " collection type:" << collec->getTypeName() << endl;
                for ( int i = 0; i < collec->getNumberOfElements(); i++)
                  {
                    try 
                      {
                        RawCalorimeterHit* rawHit = dynamic_cast<RawCalorimeterHit*>(collec->getElementAt(i));  
                        //cout << "RawCalorimeterHit[" << rawHit << "]" <<endl;
        
                    if ( rawHit != NULL )
                      {
                        CellIDDecoder<RawCalorimeterHit> idDecoder( collec ) ;
                        const BitField64& bitField = idDecoder( rawHit );
                        
                        cout << idDecoder( rawHit ).valueString() << endl;
                        lcio::long64  chamberId = bitField["CHAMBER"] ;
                        lcio::long64  difId = bitField["DIF"] ;
                        lcio::long64  chipId = bitField["CHIP"] ;
                        lcio::long64  channelId = bitField["CHANNEL"] ;
                        cout << "Chamber[" << chamberId<< "] " ;
                        cout << "Dif[" <<difId << "] " ;
                        cout << "Chip[" <<chipId << "] "  ;
                        cout << "Channel[" << channelId << "]" << endl;
                        
                        if ( geometrySet )
                          {
                            Channel &channel = detector.getDifById(difId).getChipById(chipId).getChannelById(channelId);
                            cout << " Position[X,Y,Z] = ["  << channel.getX() << ","  << channel.getY() << "," << channel.getZ() << "]" <<  endl;
                          }
                        
                        
                      } 
                    else { cout << " rawHit == NULL "<< endl; } 
                  } catch( Exception& ex)
                  {
                    cout << " an excpetion occured when trying dynamic_cast<RawCalorimeterHit*>: " << endl ;
                    cout << "  " << ex.what() << endl ;
                  }
              }
          }
        const LCParameters& params =  evt->getParameters();
        }
        
        lcReader->close() ;
        
        delete lcReader ;


}

---