Index: /trunk/Mars/mreport/MReportFitsRead.cc =================================================================== --- /trunk/Mars/mreport/MReportFitsRead.cc (revision 12741) +++ /trunk/Mars/mreport/MReportFitsRead.cc (revision 12741) @@ -0,0 +1,204 @@ +/* ======================================================================== *\ +! +! * +! * This file is part of MARS, the MAGIC Analysis and Reconstruction +! * Software. It is distributed to you in the hope that it can be a useful +! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. +! * It is distributed WITHOUT ANY WARRANTY. +! * +! * Permission to use, copy, modify and distribute this software and its +! * documentation for any purpose is hereby granted without fee, +! * provided that the above copyright notice appear in all copies and +! * that both that copyright notice and this permission notice appear +! * in supporting documentation. It is provided "as is" without express +! * or implied warranty. +! * +! +! +! Author(s): Thomas Bretz, 12/2011 +! +! Copyright: MAGIC Software Development, 2000-2011 +! +! +\* ======================================================================== */ + +////////////////////////////////////////////////////////////////////////////// +// +// MReportFitsRead +// +////////////////////////////////////////////////////////////////////////////// +#include "MReportFitsRead.h" + +#include +#include + +#include +#include +#include + +#include "fits.h" + +#include "MLog.h" +#include "MLogManip.h" + +#include "MParList.h" +#include "MReport.h" + +ClassImp(MReportFitsRead); + +using namespace std; + +// -------------------------------------------------------------------------- +// +// Default constructor. It tries to open the given file and creates a +// THashTable which allows faster access to the MReport* objects. +// +MReportFitsRead::MReportFitsRead(const char *fname, const char *name, const char *title) + : fFileName(fname), fReport(NULL), fIn(NULL) +{ + fName = name ? name : "MReportFitsRead"; + fTitle = title ? title : "Read task to read general report files"; +} + +// -------------------------------------------------------------------------- +// +// Destructor. Delete input stream and hash table. +// +MReportFitsRead::~MReportFitsRead() +{ + delete fIn; +} + +// -------------------------------------------------------------------------- +// +// Call SetupReading for all MReportHelp objects scheduled. +// Try to open the file and check the file header. +// +Int_t MReportFitsRead::PreProcess(MParList *pList) +{ + // Add the MReport instances first to the paramter list + // so that SetupReading can find them if needed + //fList->R__FOR_EACH(MReportHelp, AddToList)(*pList); + + fReport = (MReport*)pList->FindCreateObj(fReportName); + if (!fReport) + return kFALSE; + + if (!fReport->InheritsFrom(MReport::Class())) + { + *fLog << err << fReportName << " does not derive from MReport." << endl; + return kFALSE; + } + + + // + // open the input stream + // first of all check if opening the file in the constructor was + // successfull + // + if (fIn) + delete fIn; + fIn = new fits(fFileName.Data()); + if (!(*fIn)) + { + *fLog << err << "Cannot open file " << fFileName << ": "; + *fLog << strerror(errno) << endl; + return kFALSE; + } + + if (fIn->GetStr("TELESCOP")!="FACT") + { + *fLog << err << "Not a valid FACT FITS file (key TELESCOP not 'FACT')." << endl; + return kFALSE; + } + + fMjdRef = fIn->HasKey("MJDREF") ? fIn->GetUInt("MJDREF") : 0; + + return + fIn->SetRefAddress("QoS", fBufQos) && + fIn->SetRefAddress("Time", fBufTime) && + fReport->SetupReading(*pList) && + fReport->SetupReadingFits(*fIn); +} + +// -------------------------------------------------------------------------- +// +// Read the file line by line as long as a matching MReport* class is found. +// In this case call its interpreter (Interprete()) and remove the identifier +// first (XYZ-REPORT) +// +Int_t MReportFitsRead::Process() +{ + if (!fIn->GetNextRow()) + { + *fLog << inf << "End-of-file detected." << endl; + return kFALSE; + } + + if (fBufTime==0) + return kCONTINUE; + + MTime &time = *fReport->GetTime(); + + time.SetMjd(fBufTime+fMjdRef); + + // return -1: This is the special case: out of time limit + if (fStart && timefStop) + return kCONTINUE; + + const Int_t rc = fReport->InterpreteFits(*fIn); + + /* + switch (rc) + { + case kTRUE: fNumReports++; break; + case kCONTINUE: fNumSkipped++; break; + }*/ + + switch (rc) + { + case kFALSE: + *fLog << err << "ERROR - Interpreting '" << fReport->GetName() << "' failed (l." << fIn->GetRow() << ")... abort." << endl; + break; + case kCONTINUE: + *fLog << warn << "WARNING - Interpreting '" << fReport->GetName() << "' failed (l." << fIn->GetRow() << ")... skipped." << endl; + break; + case -1: // This is the special case: out of time limit + return kCONTINUE; + } + + return rc; +} + +// -------------------------------------------------------------------------- +// +// Close the file and print some execution statistics +// +Int_t MReportFitsRead::PostProcess() +{ + fIn->close(); + + if (!GetNumExecutions()) + return kTRUE; + + *fLog << inf << endl; + *fLog << GetDescriptor() << " statistics:" << endl; + *fLog << dec << setfill(' '); + /* + *fLog << inf; + *fLog << " " << setw(7) << rep->GetNumReports() << " ("; + *fLog << setw(3) << (int)(100.*rep->GetNumReports()/GetNumExecutions()); + *fLog << "%): " << rep->GetName() << endl; + + if (rep->GetNumSkipped()==0) + continue; + + *fLog << warn; + *fLog << " " << setw(7) << rep->GetNumSkipped() << " ("; + *fLog << setw(3) << (int)(100.*rep->GetNumSkipped()/GetNumExecutions()); + *fLog << "%): " << rep->GetName() << " skipped!" << endl; + */ + return kTRUE; +} Index: /trunk/Mars/mreport/MReportFitsRead.h =================================================================== --- /trunk/Mars/mreport/MReportFitsRead.h (revision 12741) +++ /trunk/Mars/mreport/MReportFitsRead.h (revision 12741) @@ -0,0 +1,64 @@ +#ifndef MARS_MReportFitsRead +#define MARS_MReportFitsRead + +#ifndef MARS_MTask +#include "MTask.h" +#endif +#ifndef MARS_MTime +#include "MTime.h" +#endif + +#include +#include + +namespace std +{ + class fits; +} +class THashTable; + +class MTime; +class MReport; +class MReportHelp; + +class MReportFitsRead : public MTask +{ +private: + TString fFileName; // Name of the input file + TString fReportName; // Name of the report class + + MReport *fReport; // List of possible reports to be interpreted + + MTime fStart; // Time range which should be read from file + MTime fStop; // Time range which should be read from file + + int32_t fBufQos; //! Buffer for Quality of Service + Double_t fBufTime; //! Buffer for Time + + UInt_t fMjdRef; + + enum { kHasNoHeader = BIT(14) }; + + Int_t PreProcess(MParList *pList); + Int_t Process(); + Int_t PostProcess(); + + virtual Int_t CheckFileHeader() { return 0; } + +protected: + std::fits *fIn; //! buffered input stream (file to read from) + +public: + MReportFitsRead(const char *filename, const char *name=NULL, const char *title=NULL); + ~MReportFitsRead(); + + void SetHasNoHeader() { SetBit(kHasNoHeader); } + void SetTimeStart(const MTime &tm) { fStart = tm; } + void SetTimeStop(const MTime &tm) { fStop = tm; } + + void SetReportName(const char *name) { fReportName = name; } + + ClassDef(MReportFitsRead, 0)// Task to read general report file +}; + +#endif