Changeset 4658 for trunk/MagicSoft/Mars/mjobs

Timestamp:

08/17/04 16:26:31 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mjobs

Files:

• MJExtractCalibTest.cc (modified) (9 diffs)
• MJExtractCalibTest.h (modified) (4 diffs)

trunk/MagicSoft/Mars/mjobs/MJExtractCalibTest.cc

@@ -27 +27 @@
 //  MJExtractCalibTest
 //
+// If the flag SetDataCheckDisplay() is set, only the most important distributions
+//  are displayed. 
+// Otherwise, (default: SetNormalDisplay()), a good selection of plots is given
+//
 /////////////////////////////////////////////////////////////////////////////
 #include "MJExtractCalibTest.h"
@@ -53 +57 @@
 #include "MCalibrationRelTimeCam.h"
 #include "MCalibrationQECam.h"
+#include "MCalibrationTestCam.h"
 #include "MCalibrationTestCalc.h"
 #include "MHCamEvent.h"
+#include "MHCalibrationTestCam.h"
 
 #include "MReadMarsFile.h"
@@ -79 +85 @@
 //
 // Sets fUseCosmicsFilter to kTRUE, fRuns to 0, fExtractor to NULL, fTimeExtractor to NULL
+// fDisplay to kNormalDisplay
 //
 MJExtractCalibTest::MJExtractCalibTest(const char *name, const char *title) 
-    : fUseCosmicsFilter(kTRUE), fRuns(NULL), fExtractor(NULL), fTimeExtractor(NULL)
+    : fUseCosmicsFilter(kTRUE), fRuns(NULL), fExtractor(NULL), fTimeExtractor(NULL),
+      fDisplayType(kNormalDisplay)
 {
     fName  = name  ? name  : "MJExtractCalibTest";
@@ -105 +113 @@
   // Get container from list
   //
-  MGeomCam &geomcam = *(MGeomCam*)plist.FindObject("MGeomCam");
+  MGeomCam  &geomcam = *(MGeomCam*) plist.FindObject("MGeomCam");
+  MHCalibrationTestCam &testcam = *(MHCalibrationTestCam*)plist.FindObject("MHCalibrationTestCam");
   
   // Create histograms to display
@@ -112 +121 @@
   MHCamera disp3 (geomcam, "Test;PhotonsPerArea",    "Equiv. Cherenkov Photons per Area");
   MHCamera disp4 (geomcam, "Test;SigmaPhotPerArea",  "Sigma equiv. Cher. Photons per Area");
+  MHCamera disp5 (geomcam, "Test;Phot",              "Calibrated Photons");
+  MHCamera disp6 (geomcam, "Test;PhotPerArea",       "Calibrated Photons per Area");
+  MHCamera disp7 (geomcam, "Test;NotInterpolate",    "Not interpolated pixels");
+  MHCamera disp8 (geomcam, "Test;DeviatingPhots",    "Deviating Number Photons");
 
   // Fitted charge means and sigmas
-  disp1.SetCamContent(fTestCam,  0);
-  disp1.SetCamError(  fTestCam,  1);
-  disp2.SetCamContent(fTestCam,  2);
-  disp2.SetCamError(  fTestCam,  3);
-  disp3.SetCamContent(fTestCam,  7);
-  disp3.SetCamError(  fTestCam,  8);
-  disp4.SetCamContent(fTestCam,  9);
-  disp4.SetCamError(  fTestCam,  10);
+  disp1.SetCamContent(testcam,  0);
+  disp1.SetCamError(  testcam,  1);
+  disp2.SetCamContent(testcam,  2);
+  disp2.SetCamError(  testcam,  3);
+  disp3.SetCamContent(testcam,  7);
+  disp3.SetCamError(  testcam,  8);
+  disp4.SetCamContent(testcam,  9);
+  disp4.SetCamError(  testcam,  10);
+
+  disp5.SetCamContent(fTestCam,  0);
+  disp5.SetCamError(  fTestCam,  1);
+  disp6.SetCamContent(fTestCam,  2);
+  disp6.SetCamError(  fTestCam,  3);
+  disp7.SetCamError(  fTestCam,  4);
+
+  disp8.SetCamError(  fBadPixels, 22);
+
 
   disp1.SetYTitle("Photons");
@@ -127 +149 @@
   disp3.SetYTitle("Photons per Area [mm^{-2}]");
   disp4.SetYTitle("\\sigma_{phot} per Area [mm^{-2}]");
+
+  disp5.SetYTitle("Photons");
+  disp6.SetYTitle("Photons per Area [mm^{-2}]");
+  disp7.SetYTitle("[1]");
+  disp8.SetYTitle("[1]");
   
   gStyle->SetOptStat(1111);
   gStyle->SetOptFit();
 
-  TCanvas &c = fDisplay->AddTab("TestCharges");
-  c.Divide(4,4);
-
-  disp1.CamDraw(c, 1, 4, 2, 1);
-  disp2.CamDraw(c, 2, 4, 2, 1);        
-  disp3.CamDraw(c, 3, 4, 1, 1);        
-  disp4.CamDraw(c, 4, 4, 2, 1);        
+  if (fDisplayType == kNormalDisplay)
+    {
+
+      TCanvas &c = fDisplay->AddTab("TestCharges");
+      c.Divide(4,4);
+      
+      disp1.CamDraw(c, 1, 4, 2, 1);
+      disp2.CamDraw(c, 2, 4, 2, 1);        
+      disp3.CamDraw(c, 3, 4, 1, 1);        
+      disp4.CamDraw(c, 4, 4, 2, 1);        
+    }
+  
+  TCanvas &c2 = fDisplay->AddTab("TestResult");
+  c2.Divide(2,4);
+
+  disp5.CamDraw(c2, 1, 2, 2, 1);
+  disp6.CamDraw(c2, 2, 2, 2, 1);        
+
+  TCanvas &c3 = fDisplay->AddTab("TestDefects");
+  c3.Divide(2,2);
+
+  disp7.CamDraw(c3, 1, 2, 0);
+  disp8.CamDraw(c3, 2, 2, 0);        
 
   return;
@@ -234 +277 @@
   *fLog << endl;
   
-  MCerPhotEvt    cerphot;
-  MPedPhotCam    pedphot;
+  MCerPhotEvt          cerphot;
+  MPedPhotCam          pedphot;
+  MHCalibrationTestCam testcam;
 
   // Setup Lists
@@ -244 +288 @@
   plist.AddToList(&cerphot);
   plist.AddToList(&pedphot);
+  plist.AddToList(&testcam);
   plist.AddToList(&fTestCam);
   plist.AddToList(&fBadPixels);
@@ -493 +538 @@
     if (fTestCam.Write()<=0)
     {
-        *fLog << err << "Unable to write MHCalibrationTestCam to " << oname << endl;
+        *fLog << err << "Unable to write MCalibrationTestCam to " << oname << endl;
         return kFALSE;
     }

trunk/MagicSoft/Mars/mjobs/MJExtractCalibTest.h

@@ -2 +2 @@
 #define MARS_MJExtractCalibTest
 
-#ifndef MARS_MHCalibrationTestCam
-#include "MHCalibrationTestCam.h"
+#ifndef MARS_MCalibrationTestCam
+#include "MCalibrationTestCam.h"
 #endif
 #ifndef MARS_MHCalibrationTestTimeCam
@@ -35 +35 @@
   
   MBadPixelsCam            fBadPixels;
-  MHCalibrationTestCam     fTestCam;
+  MCalibrationTestCam      fTestCam;
   MHCalibrationTestTimeCam fTestTimeCam;
   MPedPhotCam              fPedPhotCam;
 
+  enum  Display_t                                      // Possible Display types
+    {
+      kDataCheckDisplay,
+      kNormalDisplay
+    };
+  
+  Display_t fDisplayType;                              // Chosen Display type
+  
   Bool_t ReadPedPhotCam();
   
@@ -58 +66 @@
   const char* GetOutputFile() const;
   
-  MHCalibrationTestCam     &GetTestCam()      { return fTestCam; }
+  MCalibrationTestCam      &GetTestCam()      { return fTestCam; }
   MHCalibrationTestTimeCam &GetTestTimeCam()  { return fTestTimeCam; }  
   MPedPhotCam          &GetPedPhotCam() { return fPedPhotCam; }
@@ -67 +75 @@
   void SetBadPixels(const MBadPixelsCam &bad) { bad.Copy(fBadPixels); }
   void SetUseCosmicsFilter( const Bool_t b )    { fUseCosmicsFilter = b;  }
-     
+
+  // Displays
+  void SetDataCheckDisplay() { fDisplayType = kDataCheckDisplay; }
+  void SetNormalDisplay()    { fDisplayType = kNormalDisplay;    }
+  
   Bool_t ProcessD(MPedestalCam &pedcam, MCalibrationChargeCam &calcam, MCalibrationQECam &qecam);
   Bool_t ProcessT(MPedestalCam &pedcam, MCalibrationRelTimeCam &relcam);