Changeset 3503 for trunk/MagicSoft/Mars/mmain/MEventDisplay.cc

Timestamp:

03/15/04 15:12:34 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mmain/MEventDisplay.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/mmain/MEventDisplay.cc

@@ -32 +32 @@
 // root
 //
+#include <TFile.h>               // TFile
 #include <TList.h>               // TList::Add
 #include <TStyle.h>              // gStyle->SetOptStat
@@ -78 +79 @@
 // Container
 //
-#include "MHillas.h"             // MHillas::Print(const MGeomCam&)
-#include "MHillasExt.h"          // MHillasExt::Print(const MGeomCam&)
-#include "MHillasSrc.h"          // MHillasSrc::Print(const MGeomCam&)
-#include "MHEvent.h"             // MHEvent
-#include "MHCamera.h"            // MHCamera
-#include "MRawEvtData.h"         // MRawEvtData
-#include "MArrivalTimeCam.h"     // MArrivalTimeCam
-#include "MBadPixelsCam.h"       // MBadPixelsCam
+#include "MHillas.h"               // MHillas::Print(const MGeomCam&)
+#include "MHillasExt.h"            // MHillasExt::Print(const MGeomCam&)
+#include "MHillasSrc.h"            // MHillasSrc::Print(const MGeomCam&)
+#include "MHEvent.h"               // MHEvent
+#include "MHCamera.h"              // MHCamera
+#include "MRawEvtData.h"           // MRawEvtData
+#include "MArrivalTimeCam.h"       // MArrivalTimeCam
+#include "MBadPixelsCam.h"         // MBadPixelsCam
+#include "MPedPhotCam.h"           // MPedPhotCam
+#include "MCalibrationChargeCam.h" // MCalibrationChargeCam
 
 ClassImp(MEventDisplay);
@@ -93 +96 @@
 //  Constructor.
 //
-MEventDisplay::MEventDisplay(const char *filename) : MStatusDisplay()
+MEventDisplay::MEventDisplay(const char *fname, const char *pname, const char *cname) : MStatusDisplay()
 {
     //
     // Setup Task list for hillas calculation
     //
-    SetupTaskList("Events", filename);
+    SetupTaskList("Events", fname, pname, cname);
 
     //
     // Add MEventDisplay GUI elements to the display
     //
-    AddUserFrame(filename);
+    AddUserFrame(fname);
 
     //
@@ -134 +137 @@
 //  preprocess tasks and read in first event (process)
 //
-void MEventDisplay::SetupTaskList(const char *tname, const char *fname)
-{
+void MEventDisplay::SetupTaskList(const char *tname, const char *fname,
+                                  const char *pname, const char *cname)
+{
+    MCalibrationChargeCam *ccam=0;
+    MPedPhotCam           *pcam=0;
+
+    MBadPixelsCam *badpix = new MBadPixelsCam;
+
+    TFile file(pname, "READ");
+    if (!file.IsZombie())
+        pcam = new MPedPhotCam;
+    if (pcam)
+    {
+        if (pcam->Read()<=0)
+        {
+            delete pcam;
+            pcam = NULL;
+        }
+
+        if (file.FindKey("MBadPixelsCam"))
+        {
+            MBadPixelsCam bad;
+            if (bad.Read()>0)
+                badpix->Merge(bad);
+        }
+    }
+    file.Close();
+    file.Open(cname, "READ");
+    if (!file.IsZombie())
+        ccam = new MCalibrationChargeCam;
+    if (ccam)
+    {
+        if (ccam->Read()<=0)
+        {
+            delete ccam;
+            ccam = NULL;
+        }
+
+        if (file.FindKey("MBadPixelsCam"))
+        {
+            MBadPixelsCam bad;
+            if (bad.Read()>0)
+                badpix->Merge(bad);
+        }
+    }
+    file.Close();
+
     //
     // Setup an empty job, with a reader task only.
@@ -154 +202 @@
     plist->SetOwner();
     plist->AddToList(tlist);
-
-    MBadPixelsCam *badpix = new MBadPixelsCam;
     plist->AddToList(badpix);
+    if (pcam)
+        plist->AddToList(pcam);
+    if (ccam)
+        plist->AddToList(ccam);
 
     fEvtLoop = new MEvtLoop;
@@ -180 +230 @@
     evt8->SetName("Arrival Time");
 
+    // This makes sure, that the containers are deleted...
     plist->AddToList(evt1);
     plist->AddToList(evt2);
@@ -189 +240 @@
     plist->AddToList(evt8);
 
-
-    MMcPedestalCopy   *pcopy = new MMcPedestalCopy;
-    MMcPedestalNSBAdd *pdnsb = new MMcPedestalNSBAdd;
     MCerPhotAnal2     *nanal = new MCerPhotAnal2;
     MFillH            *fill1 = new MFillH(evt1, "MCerPhotEvt", "MFillH1");
@@ -203 +251 @@
     MHillasCalc       *hcalc = new MHillasCalc;
     MHillasSrcCalc    *scalc = new MHillasSrcCalc;
-    MArrivalTimeCalc  *tcalc = new MArrivalTimeCalc;
-    MFillH            *fill7 = new MFillH(evt7, "MRawEvtData",  "MFillH7");
-    MFillH            *fill8 = new MFillH(evt8, "MArrivalTimeCam", "MFillH8");
-
-    MFilter *f1=new MFDataMember("MRawRunHeader.fRunType", '>', 255.5);
-    MFilter *f2=new MFDataMember("MRawRunHeader.fRunType", '<', 255.5);
-    f1->SetName("MFMonteCarlo");
-    f2->SetName("MFRealData");
-
-    MExtractSignal* sigextract = new MExtractSignal();
-    sigextract->SetRange(5, 9, 5, 9);
-
-    MMcCalibrationUpdate* mccalibupdate = new MMcCalibrationUpdate;
-    mccalibupdate->SetOuterPixelsGainScaling(kFALSE);
-    MCalibrate* mccalib = new MCalibrate;
-
-    // MC
-    sigextract->SetFilter(f1);
-    mccalibupdate->SetFilter(f1);
-    mccalib->SetFilter(f1);
-
-    // Data
-    nanal->SetFilter(f2);
-
-    tlist->AddToList(f1);
-    tlist->AddToList(f2);
-    tlist->AddToList(pcopy);
-    tlist->AddToList(pdnsb);
-
-    tlist->AddToList(sigextract);
-    tlist->AddToList(mccalibupdate);
-    tlist->AddToList(mccalib);
-
-    tlist->AddToList(nanal);
+
+    // If no pedestal or no calibration file is availble
+    if (!pcam || !ccam)
+    {
+        MFilter *f1=new MFDataMember("MRawRunHeader.fRunType", '>', 255.5);
+        MFilter *f2=new MFDataMember("MRawRunHeader.fRunType", '<', 255.5);
+        f1->SetName("MFMonteCarlo");
+        f2->SetName("MFRealData");
+
+        MMcPedestalCopy   *pcopy = new MMcPedestalCopy;
+        MMcPedestalNSBAdd *pdnsb = new MMcPedestalNSBAdd;
+
+        MExtractSignal* extra = new MExtractSignal();
+        extra->SetRange(5, 9, 5, 9);
+
+        MMcCalibrationUpdate* mcupd = new MMcCalibrationUpdate;
+        mcupd->SetOuterPixelsGainScaling(kFALSE);
+
+        MCalibrate* mccal = new MCalibrate;
+
+        // MC
+        extra->SetFilter(f1);
+        mcupd->SetFilter(f1);
+        mccal->SetFilter(f1);
+
+        // Data
+        nanal->SetFilter(f2);
+
+        // TaskList
+        tlist->AddToList(f1);
+        tlist->AddToList(f2);
+        tlist->AddToList(pcopy);
+        tlist->AddToList(pdnsb);
+
+        tlist->AddToList(extra);
+        tlist->AddToList(mcupd);
+        tlist->AddToList(mccal);
+
+        tlist->AddToList(nanal);
+    }
+    else
+    {
+        MCalibrate* calib = new MCalibrate;
+        tlist->AddToList(calib);
+    }
 
     tlist->AddToList(fill1);
@@ -248 +307 @@
     tlist->AddToList(hcalc);
     tlist->AddToList(scalc);
-    tlist->AddToList(tcalc);
-    tlist->AddToList(fill7);
-    tlist->AddToList(fill8);
+    if (!pcam || !ccam)
+    {
+        MArrivalTimeCalc  *tcalc = new MArrivalTimeCalc;
+        MFillH            *fill7 = new MFillH(evt7, "MRawEvtData",     "MFillH7");
+        MFillH            *fill8 = new MFillH(evt8, "MArrivalTimeCam", "MFillH8");
+        tlist->AddToList(tcalc);
+        tlist->AddToList(fill7);
+        tlist->AddToList(fill8);
+    }
 
     //
@@ -408 +473 @@
     MGeomCam    *geom  = (MGeomCam*)   plist->FindObject("MGeomCam");
     MRawEvtData *raw   = (MRawEvtData*)plist->FindObject("MRawEvtData");
-    if (!raw)
-        return;
 
     //
@@ -436 +499 @@
         reader->DecEventNum();
 
-    } while (raw->GetNumPixels()<1 && dir!=0);
+    } while (raw && raw->GetNumPixels()<1 && dir!=0);
 
     //