Changeset 3855

Timestamp:

04/27/04 19:24:37 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• macros/calibration.C (modified) (4 diffs)
• mcalib/MHCalibrationTestCam.cc (modified) (2 diffs)
• mjobs/JobsLinkDef.h (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -63 +63 @@
        was not correct
 
+   * macros/calibration.C
+     - added the test class MJExtractCalibTest
 
  2004/04/27: Abelardo Moralejo

trunk/MagicSoft/Mars/macros/calibration.C

@@ -51 +51 @@
 //  uncommented as well. 
 // 
+//  Last, a third loop is performed over the calibration file again in order to 
+//  "calibrate" it and test the resulting outcome.
+//
 /////////////////////////////////////////////////////////////////////////////
-const TString inpath = "/home/rootdata/Calib/2004_04_16/";
-const Int_t pedrun  = 22265;
-const Int_t calrun1 = 22299;
-const Int_t calrun2 = 22300;
+static const TString inpath = "/home/rootdata/Calib/2004_04_16/";
+static const Int_t pedrun  = 22265;
+static const Int_t calrun1 = 22300;
+static const Int_t calrun2 = 0;
 
 void calibration(const Int_t prun=pedrun, const Int_t crun1=calrun1, const Int_t crun2=calrun2)
 {
+
   
   MRunIter pruns;
@@ -72 +76 @@
   MCalibrationQECam qecam;
   MBadPixelsCam     badcam;
-
+  MGeomCamMagic     geomcam;
+  MGeomApply        geomapl;
   //
   // If you want to exclude pixels from the beginning, read 
@@ -95 +100 @@
     return;
 
-  //
-  // Now setup the tasks and for the calibration:
-  // ---------------------------------------------------
-  //
+  /****************************************/
+  /* SECOND LOOP: CALIBRATION COMPUTATION */
+  /****************************************/
+
   MJCalibration calloop;
-  //  calloop.SetColor(color);
 
   //
   // If you want to run the data-check on RAW DATA!!!, choose:
   // calloop.SetDataCheck();
+  // 
+  // If you want to see the data-check plots only, choose:
+  // calloop.SetDataCheckDisplay();
   // 
   // For everything, you ever dreamed of, choose:
@@ -119 +126 @@
   if (!calloop.Process(pedloop.GetPedestalCam()))
     return;
-  
+
+
+  /********************************************************************/
+  /* THIRD LOOP: APPLY CALIBRATION TO THE CALIBRATION FILES FOR TESTS */
+  /********************************************************************/
+
+  MJExtractCalibTest testloop;
+
+  testloop.SetInput(&cruns);
+  testloop.SetDisplay(display);
+  testloop.SetBadPixels(calloop.GetBadPixels());
+  if (!testloop.ProcessD(pedloop.GetPedestalCam(),calloop.GetCalibrationCam(),calloop.GetQECam()))
+    return;
+
+
 }
 

trunk/MagicSoft/Mars/mcalib/MHCalibrationTestCam.cc

@@ -453 +453 @@
 // 6: Fourier spectrum not OK                       (calls: MHGausEvents::IsFourierSpectrumOK())
 //
+// Converted values:
+// =================
+//
+// 7: Fitted Mean Test Calibration (MHGausEvents::GetMean()) by MGeomPix::GetA()
+// 8: Fitted Mean Error Calibration (MHGausEvents::GetSigma()) by MGeomPix::GetA()
+// 9: Fitted Sigma Test Calibration (MHGausEvents::GetMean()) by MGeomPix::GetA()
+// 8: Fitted Sigma Error Calibration (MHGausEvents::GetSigma()) by MGeomPix::GetA()
+//
 Bool_t MHCalibrationTestCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
 {
@@ -488 +496 @@
       if (!pix.IsFourierSpectrumOK())
         val = 1.;
+      break;
+    case 7:
+      val = pix.GetMean()/cam[idx].GetA();
+      break;
+    case 8:
+      val = pix.GetMeanErr()/cam[idx].GetA();
+      break;
+    case 9:
+      val = pix.GetSigma()/cam[idx].GetA();
+      break;
+    case 10:
+      val = pix.GetSigmaErr()/cam[idx].GetA();
       break;
     default:

trunk/MagicSoft/Mars/mjobs/JobsLinkDef.h

@@ -8 +8 @@
 #pragma link C++ class MJCalibration+;
 #pragma link C++ class MJExtractSignal+;
+#pragma link C++ class MJExtractCalibTest+;
 
 #pragma link C++ class MGCamDisplays+;