Changeset 1342 for trunk/MagicSoft/Mars/mbase/MReadMarsFile.cc

Timestamp:

06/03/02 15:43:45 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mbase/MReadMarsFile.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mbase/MReadMarsFile.cc

@@ -41 +41 @@
 #include "MParList.h"
 #include "MTaskList.h"
+
+#include "MRawRunHeader.h"
+
+#include "MMcRunHeader.hxx"
 
 ClassImp(MReadMarsFile);
@@ -109 +113 @@
 {
     Int_t runtype = -1;
-/*
-    const MRawRunHeader *runheader = (MRawRunHeader*)fParList->FindObject("MRawRunHeader");
-    if (runheader)
-        runtype = runheader->GetRunType();
-*/
+
+    const MRawRunHeader *rawheader = (MRawRunHeader*)fParList->FindObject("MRawRunHeader");
+    if (rawheader)
+        runtype = rawheader->GetRunType();
+
     //
     // Try to read the new run headers. If reading the new run header
@@ -120 +124 @@
     const Int_t idx = GetFileIndex();
     fRun->SetEventNum(idx<0?0:idx); // Assumption: One Entry per File!
-    if (fRun->Process())
+    if (!fRun->Process())
     {
-        *fLog << inf << "MReadMarsFile: Switching to #" << GetFileIndex();
-        *fLog << " '" << GetFileName() << "' (before event #";
-        *fLog << GetEventNum()-1 << ")" << endl;
-/*
-        if (runheader)
-            if (runtype != runheader->GetRunType())
-                *fLog << warn << "Warning - You are mixing files with different run types!" << endl;
-*/
-        if (fTaskList->ReInit())
-            return kTRUE;
+        *fLog << err << "ERROR - Cannot read new runheaders #" << idx;
+        *fLog << " after reading event #" << GetEventNum() << endl;
+        return kFALSE;
+    }
 
-        //MReadTree::Notify();
+    *fLog << inf << "MReadMarsFile: Switching to #" << GetFileIndex();
+    *fLog << " '" << GetFileName() << "' (before event #";
+    *fLog << GetEventNum()-1 << ")" << endl;
 
-        *fLog << err << "ERROR - ReInit of '" << fTaskList->GetName() << "' failed." << endl;
+    if (rawheader)
+    {
+        if (runtype != rawheader->GetRunType())
+            *fLog << warn << "Warning - You are mixing files with different run types!" << endl;
     }
-    else
-        *fLog << err << "ERROR - Cannot read new runheaders after reading event #" << GetEventNum() << endl;
 
+    const MMcRunHeader *mcheader = (MMcRunHeader*)fParList->FindObject("MMcRunHeader");
+    if (mcheader)
+    {
+        if (mcheader->GetReflVersion()<=40 && mcheader->GetTelesTheta()>15)
+        {
+            *fLog << warn << "Warning - You may use Monte Carlo data produced with a version" << endl;
+            *fLog << "of the reflector program < 0.4 and a zenith angle > 15," << endl;
+            *fLog << "in this case you may get less photons than you would expect." << endl;
+        }
+    }
+
+    MTaskList *tlist = (MTaskList*)fParList->FindObject("MTaskList");
+    if (!tlist)
+    {
+        *fLog << err << dbginf << "ERROR - Task List not found in Parameter List." << endl;
+        return kFALSE;
+    }
+
+    if (tlist->ReInit())
+        return kTRUE;
+
+    //MReadTree::Notify();
+
+    *fLog << err << "ERROR - ReInit of '" << tlist->GetName() << "' failed." << endl;
     return kFALSE;
 }
@@ -151 +176 @@
 {
     fParList  = pList;
-    fTaskList = (MTaskList*)pList->FindObject("MTaskList");
 
     if (!fRun->PreProcess(pList))