testFDasEventFile.cc

Go to the documentation of this file.

#include <io/FDasEventFile.h>
#include <evt/Event.h>

#include <fwk/CentralConfig.h>
#include <det/Detector.h>

#include <iostream>
#include <utl/Trace.h>
#include <utl/Photon.h>
#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/ErrorLogger.h>
#include <fevt/FEvent.h>
#include <fevt/Eye.h>
#include <fevt/Telescope.h>
#include <fevt/Channel.h>
#include <fevt/Pixel.h>
#include <fevt/ChannelSimData.h>
#include <fevt/PixelSimData.h>
#include <fevt/PixelRecData.h>

#include <AugerEvent.h>
#include <EyeEvent.hh>

#include <boost/filesystem/operations.hpp>
#include <boost/filesystem/path.hpp>

namespace fs = boost::filesystem;
using namespace std;
using namespace io;
using namespace utl;
using namespace evt;
using namespace fevt;
using namespace det;


#define PRINTFADCTRACES 0


void Dump(const FEvent& fevent);


int
main(int /*argc*/, char* argv[])
{
  ErrorLogger::GetInstance().SetVerbosity(Verbosity::eVerbose);

  // Set up the detector description
  fwk::CentralConfig::GetInstance(BOOTSTRAPFILE);

  TimeStamp theTime = UTCDateTime(2002,4,10).GetTimeStamp();
  det::Detector::GetInstance().Update(theTime);
  cout << endl << "Starting " << argv[0] << endl << endl;

  string theInputFilename = TESTFILE;

  FDasEventFile file(theInputFilename, io::eRead);

  evt::Event event;

  file.Read(event);

  AugerEvent& rawevent = event.GetRawEvent();

  for (AugerEvent::EyeIterator it = rawevent.EyesBegin();
       it != rawevent.EyesEnd(); ++it)
    it->Dump();

  cout << " done" << endl;
}


/* This test has to be rewritten , sa
void Dump(const FEvent& fevent){
  cout << endl;
  cout << " ******* Now starting readback *********" << endl;
  cout << endl;

// iterate over eyes

  FEvent::ConstEyeIterator iter;
  for (iter = fevent.EyesBegin(); iter != fevent.EyesEnd();iter++){

    Eye& eye = *iter;

// iterate over telescopes

    Eye::ConstTelescopeIterator teliter;

    for (teliter = eye.TelescopesBegin();
         teliter!=eye.TelescopesEnd(); teliter++){

      Telescope& tel = *teliter;

            if (tel.HasTriggerData()) {

        TelescopeTriggerData& trig = tel.GetTriggerData();

//#define PRINTSLT 1


        for (int a=1; a<= 440 ; a++){

          for (unsigned int i=0; i< TelescopeTriggerData::GetNTriggerBins(); i++){
            int dummy = trig.HasFLT(i,a);
#if PRINSLT
            cout << dummy;
#endif
          }
#if PRINSLT
          cout << endl;
#endif
        }

      }

      Telescope::ConstChannelIterator chiter;

// iterate over channels

      for (chiter = tel.ChannelsBegin();
           chiter != tel.ChannelsEnd(); chiter++){
        Channel& ch = *chiter;

#if PRINTFADCTRACES
                cout << endl;
                cout << "Examining channel "<< ch.GetId()
                     << " Tel " << ch.GetTelescopeId()
                     << " Eye " << ch.GetEyeId()<< endl;

#else
                ;
#endif

// retrieve raw fadc trace

        if (!ch.HasRawData()){ERROR("Raw Data not found");exit(1);}
           ChannelRawData& rawdata= ch.GetRawData();

        if (!rawdata.HasFADCTrace()){ERROR("Raw FADC not found");exit(1);}


        TraceI& rawtrace = rawdata.GetFADCTrace();

#if PRINTFADCTRACES
        cout << "    raw fadc values: " ;
#endif

       for (unsigned int i = rawtrace.GetStart(); i<= rawtrace.GetStop(); i++ )
#if PRINTFADCTRACES
             cout << "[" << i<< "]" << rawtrace[i]<< " " ;

        cout << endl;

#else
             int  a=rawtrace[i] ;
#endif


      }// loop on channels

      cout << endl;

// iterate over pixel

      Telescope::ConstPixelIterator pixiter;

      for (pixiter = tel->PixelsBegin(); pixiter != tel->PixelsEnd(); pixiter++){

        Pixel * pix = *pixiter;
        cout << endl;
        cout << "Examining pixel "<< pix->GetId() << " Tel " << pix->GetTelescopeId() << " Eye " \
            << pix->GetEyeId()<< endl;

// retrieve reconstructed photon trace

        TraceD * recph = pix->GetRecData()->GetPhotonTrace();

        cout << "    reconstructed photons: " ;

        for (int i = recph->GetStart(); i<= recph->GetStop(); i++ )
           cout << "[" << i<< "]" << (*recph)[i]<< " " ;

        cout <<endl;

// retrieve reconstructed Fluorescence pulse information

        cout << "    Pulse Start "<< pix->GetRecData()->GetPulseStart();
        cout <<     " Pulse Stop  "<< pix->GetRecData()->GetPulseStop();
        cout <<     " Charge " << pix->GetRecData()->GetTotalCharge();
        cout << endl;

// retrieve Mie scattered component of reconstructed photon trace

        TraceD * recphScat = pix->GetRecData()->GetPhotonTrace(Telescope::eMieScattered);

        cout << "    reconstructed photons Mie Scattered: " ;

        for (int i = recphScat->GetStart(); i<= recphScat->GetStop(); i++ )
           cout << "[" << i<< "]" << (*recphScat)[i]<< " " ;

        cout << endl;

// retrieve simuated photon trace

        TraceD * simph = pix->GetSimData()->GetPhotonTrace();

        cout << "    Simulated photons: " ;

        for (int i = simph->GetStart(); i<= simph->GetStop(); i++ )
           cout << "[" << i<< "]" << (*simph)[i]<< " " ;
        cout << endl;

// retrieve direct component of simulated photon trace

        TraceD * simphDir = pix->GetSimData()->GetPhotonTrace(Telescope::eDirect);
        cout << "    Simulated photons Direct component: " ;


        for (int i = simphDir->GetStart(); i<= simphDir->GetStop(); i++ )
           cout << "[" << i<< "]" << (*simphDir)[i]<< " " ;
           cout << "" <<endl;

      }// loop on pixels

    }// loop on telescopes

  }// loop on eyes


}

*/

// Local Variables:
// mode:c++
// compile-command: "make -C .. testFDasEventFile  -k"
// End: