Changeset 1664 for trunk/MagicSoft/Mars/mfileio/MCT1ReadPreProc.cc

Timestamp:

11/21/02 16:20:21 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mfileio/MCT1ReadPreProc.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/mfileio/MCT1ReadPreProc.cc

@@ -112 +112 @@
 void MCT1ReadPreProc::AddFile(const char *txt)
 {
-    ifstream *store = fIn;
-
-    fIn = new ifstream(gSystem->ExpandPathName(txt));
-
-    if (!(*fIn))
-    {
-        *fLog << warn << "Cannot open file '" << txt << "'... ignored." << endl;
-        fIn = store;
+    const TString fname = gSystem->ExpandPathName(txt);
+
+    if (!CheckHeader(fname))
+    {
+        *fLog << warn << "WARNING - Problem reading header... ignored." << endl;
         return;
     }
 
-    fEntries += GetNumEvents();
-
-    delete fIn;
-
-    fIn = store;
+    const Int_t n = GetNumEvents(fname);
+    if (n==0)
+    {
+        *fLog << warn << "WARNING - File contains no data... ignored." << endl;
+        return;
+    }
+
+    fEntries += n;
+
+    *fLog << inf << "File " << txt << " contains " << n << " events (Total=" << fEntries << ")" << endl;
 
     fFileNames->AddLast(new TNamed(txt, ""));
@@ -390 +392 @@
            &dummy/*&outpars.eruntype*/,   &outpars.iproc_mjdate);
 
-    if (STRUCT_VERSION != structversion)
+    if (fpreprocversion<0.6)
+    {
+        *fLog << err << "Sorry only files from PreProc V0.6 and newer are supported." << endl;
+        return kFALSE;
+    }
+
+    //
+    // This is a stupid way of getting rid of numerical uncertanties when
+    // comparing floating point numbers (Argh...)
+    //
+    TString s1 = Form("%.2f", structversion);
+    TString s2 = Form("%.2f", STRUCT_VERSION);
+
+    if (s1 != s2)
     {
         *fLog << warn << "WARNING: Version of C-structures of file (V";
-        *fLog << structversion << ") not identical with current structures (V";
-        *fLog << STRUCT_VERSION << ")" << endl;
+        *fLog << s1 << ") not identical with current structures (V";
+        *fLog << s2 << ")" << endl;
     }
 
@@ -404 +419 @@
 }
 
-Bool_t MCT1ReadPreProc::ReadRunFooter()
+Int_t MCT1ReadPreProc::ReadRunFooter()
 {
     char cheadertitle[iHEADERTITLELENGTH];
@@ -420 +435 @@
     {
         fIn->seekg(-iHEADERTITLELENGTH, ios::cur);
-        return kFALSE;
+        return 0;
     }
 
@@ -444 +459 @@
      */
 
-    fNumFilterEvts += filterres.ifilter_passed_evts;
+    if (fNumEventsInRun!=(UInt_t)filterres.ifilter_passed_evts)
+    {
+        *fLog << err << "ERROR! Number of events in run (" << (UInt_t)filterres.ifilter_passed_evts;
+        *fLog << ") doesn't match number of read events (";
+        *fLog << fNumEventsInRun << ")" << endl;
+        *fLog << "       File corrupted." << endl;
+        return -1;
+    }
+
+    fNumFilterEvts += fNumEventsInRun;
     fNumRuns++;
 
-    *fLog << inf << "Read " << fNumEventsInRun << " events from run." << endl;
-
-    if (fNumEventsInRun!=(UInt_t)filterres.ifilter_passed_evts)
-    {
-        *fLog << warn << "WARNING! Number of events in run doesn't match number of read events." << endl;
-        *fLog << "         File might be corrupt." << endl;
-    }
-
-    return kTRUE;
+    *fLog << inf << "Read " << fNumEventsInRun << " events from run (Total=";
+    *fLog << fNumFilterEvts << "/" << fEntries << " [";
+    *fLog << 100*fNumFilterEvts/fEntries << "%], Runs=" << fNumRuns << ")";
+    *fLog << endl;
+
+    return 1;
 }
 
@@ -481 +502 @@
     // open the file which is the first one in the chain
     //
-    const char *name = file->GetName();
-
-    fIn = new ifstream(gSystem->ExpandPathName(name));
-
-    if (!(*fIn))
-    {
-        *fLog << dbginf << "Cannot open file '" << name << "'" << endl;
-        return kFALSE;
-    }
-
-    *fLog << "Open file: '" << name << "'" << endl;
+    const TString name  = file->GetName();
+    const TString fname = gSystem->ExpandPathName(name);
 
     //
@@ -498 +510 @@
     fFileNames->Remove(file);
 
+    *fLog << inf << "Open file: '" << name << "'" << endl;
+
+    if (!CheckHeader(fname))
+        return kFALSE;
+
+    fIn = new ifstream(fname);
+
     *fLog << inf << "-----------------------------------------------------------------------" << endl;
 
-    *fLog << "File contains " << GetNumEvents() << " events." << endl;
-
-    // WORKAROUND for not working seekg(0) in GetNumEvents
-    fIn->close();
-    fIn->open(gSystem->ExpandPathName(name));
-
-    Bool_t rc = ReadRunHeader();
-
-    if (!rc)
+    if (!ReadRunHeader())
+    {
         *fLog << warn << "Unable to read first run header... skipping file." << endl;
-
-    return rc;
-}
-
-Int_t MCT1ReadPreProc::GetNumEvents()
-{
-    *fLog << inf << "Scanning file for size" << flush;
+        return kFALSE;
+    }
+
+    return kTRUE;
+}
+
+Bool_t MCT1ReadPreProc::CheckHeader(const TString fname) const
+{
+    ifstream fin(fname);
+    if (!fin)
+    {
+        *fLog << dbginf << err << "ERROR - Cannot open file '" << fname << "'" << endl;
+        return kFALSE;
+    }
+
+    char cheadertitle[iHEADERTITLELENGTH];
+    fin.read(cheadertitle, iHEADERTITLELENGTH);
+
+    Float_t fpreprocversion, structversion, dummy;
+
+    sscanf(cheadertitle, cTITLE_TEMPLATE,
+           &fpreprocversion, &structversion,
+           &dummy, &dummy, &dummy, &dummy);
+
+    if (fpreprocversion < 0.6)
+    {
+        *fLog << dbginf << err << "ERROR - You must use PreProc V0.6 or higher." << endl;
+        return kFALSE;
+    }
+
+    if (STRUCT_VERSION > structversion)
+    {
+        *fLog << warn << "WARNING: Version of C-structures of file (V";
+        *fLog << structversion << ") newer than current structures (V";
+        *fLog << STRUCT_VERSION << ")" << endl;
+    }
+
+    *fLog << "Current structures: " << STRUCT_VERSION << "  ";
+    *fLog << "Structures in file: " << structversion << "  ";
+    *fLog << "Used preproc version: " << fpreprocversion << endl;
+
+    return kTRUE;
+}
+
+
+Int_t MCT1ReadPreProc::GetNumEvents(const TString fname) const
+{
+    *fLog << inf << "Scanning file " << fname << " for size" << flush;
+
+    ifstream fin(fname);
+    if (!fin)
+    {
+        *fLog << dbginf << err << "ERROR - Opening file." << endl;
+        return 0;
+    }
 
     const TString m(cEND_EVENTS_TEMPLATE);
@@ -525 +585 @@
     Int_t nruns = 0;
 
-    while (!fIn->eof())
-    {
-        fIn->seekg(iHEADERTITLELENGTH, ios::cur);
-        fIn->seekg(sizeof(struct outputpars), ios::cur);
+    while (!fin.eof() && fin.peek()!=EOF)
+    {
+        fin.seekg(iHEADERTITLELENGTH, ios::cur);
+        fin.seekg(sizeof(struct outputpars), ios::cur);
 
         while (1)
         {
-            if (fIn->peek()==cEND_EVENTS_TEMPLATE[0])
+            if (fin.peek()==cEND_EVENTS_TEMPLATE[0])
             {
                 char cheadertitle[iHEADERTITLELENGTH];
-                fIn->read(cheadertitle, iHEADERTITLELENGTH);
+                fin.read(cheadertitle, iHEADERTITLELENGTH);
 
                 const TString s = cheadertitle;
                 if (s.BeginsWith(test))
                 {
-                    fIn->seekg(sizeof(struct filterresults), ios::cur);
+                    fin.seekg(sizeof(struct filterresults), ios::cur);
                     nruns++;
                     break;
                 }
 
-                fIn->seekg(-iHEADERTITLELENGTH, ios::cur);
+                fin.seekg(-iHEADERTITLELENGTH, ios::cur);
             }
 
-            fIn->seekg(sizeof(struct eventrecord), ios::cur);
-            if (fIn->eof())
+            fin.seekg(sizeof(struct eventrecord), ios::cur);
+            if (fin.eof())
                 break;
+
             nevts++;
         }
@@ -746 +807 @@
         // must be an event
         //
-        if (!ReadRunFooter())
+        switch (ReadRunFooter())
+        {
+        case -1:
+            return kFALSE;
+        case 0:
             return kTRUE;
+        }
 
         *fLog << inf << "Footer found." << endl;
+
+        const char c = fIn->peek();
 
         //
         // No after reading the footer check if we reached the end of the file
         //
-        if (fIn->eof())
+        if (fIn->eof() || c==EOF)
         {
             *fLog << "End of file." << endl;
@@ -763 +831 @@
         // If the eof isn't reached a new header must follow. Check for it.
         //
-        if (fIn->peek()!=cTITLE_TEMPLATE[0])
+        if (c!=cTITLE_TEMPLATE[0])
         {
             *fLog << inf << "Error finding new run header in file (possible EOF)... skipping rest of file." << endl;
@@ -813 +881 @@
     *fLog << "Number of Events read from file: " << fNumEvents << endl;
     *fLog << "Number of Runs read from file:   " << fNumRuns << endl;
+    *fLog << "Number of events detected first: " << fEntries << endl;
 
     if (fNumEvents!=fNumFilterEvts)