Changeset 3725 for trunk/MagicSoft/Mars/macros

Timestamp:

04/13/04 20:03:32 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/calibration.C (modified) (11 diffs)

trunk/MagicSoft/Mars/macros/calibration.C

@@ -22 +22 @@
 !
 \* ======================================================================== */
-#include "MAGIC.h"
 
 void calibration()
 {
+  
+  const MCalibrationCam::PulserColor_t color = MCalibrationCam::kGREEN;
 
   //  const TString inpath = "/mnt/Data/rootdata/CrabNebula/2004_01_27/";
@@ -41 +42 @@
   gStyle->SetOptFit();
   
+  MCalibrationQECam qecam;
+  MBadPixelsCam     badcam;
+  badcam.AsciiRead("badpixels.dat");
+
+  for (Int_t i=0;i<badcam.GetSize();i++)
+    if (badcam[i].IsBad())
+      cout << "Bad Pixel: " << i << endl;
+
   MStatusDisplay *display = new MStatusDisplay;
   display->SetUpdateTime(3000);
   display->Resize(850,700);
   
-  MBadPixelsCam badcam;
-  badcam.AsciiRead("badpixels.dat");
-
   /************************************/
   /* FIRST LOOP: PEDESTAL COMPUTATION */
   /************************************/
-  
+
   MJPedestal pedloop;
   pedloop.SetInput(&pruns);
@@ -76 +82 @@
   static_cast<MRead&>(read).AddFiles(pruns);
 
-  MGeomApply      geomapl;
-  MExtractSignal  sigcalc;
-  
-  //
-  // Additionally to calculating the pedestals, 
-  // you can fill histograms and look at them
-  //
-  MFillH fill("MHPedestalCam", "MExtractedSignalCam");
-  
-  tlist0.AddToList(&read);
-  tlist0.AddToList(&geomapl);
-  tlist0.AddToList(&sigcalc);
-  tlist0.AddToList(&fill);
-  
-  MGeomCamMagic  geomcam;
-  MHPedestalCam  hpedcam;
-  plist0.AddToList(&geomcam);
-  plist0.AddToList(&pedloop.GetPedestalCam());
-  plist0.AddToList(&badcam);
-  plist0.AddToList(&hpedcam);
-  
-  //
-  // Create and setup the eventloop
-  //
-  MEvtLoop evtloop0;
-  
-  evtloop0.SetParList(&plist0);
-  evtloop0.SetDisplay(display);
-  
-  //
-  // Execute first analysis
-  //
-  if (!evtloop0.Eventloop())
-    return;
-
-  tlist0.PrintStatistics();
-
-  for (Int_t aidx=0;aidx<2;aidx++)
-    {
-      hpedcam.GetAverageHiGainArea(aidx).DrawClone("fourierevents");
-      hpedcam.GetAverageLoGainArea(aidx).DrawClone("fourierevents");
-    }
-  
-  for (Int_t sector=1;sector<7;sector++)
-    {
-      hpedcam.GetAverageHiGainSector(sector).DrawClone("fourierevents");
-      hpedcam.GetAverageLoGainSector(sector).DrawClone("fourierevents");
-    }
-  
-
   // 
   // Now the short version: 
   //
-  /*
   MJCalibration calloop;
+  calloop.SetColor(color);
   calloop.SetInput(&cruns);
   calloop.SetDisplay(display);
-  calloop.SetBadPixels(pedloop.GetBadPixelsCam());
+  calloop.SetQECam(qecam);
+  calloop.SetBadPixels(pedloop.GetBadPixels());
   if (!calloop.Process(pedloop.GetPedestalCam()))
     return;
-    #if 0
-  */
+
+#if 0
   // 
   // The longer version:
@@ -163 +120 @@
   MArrivalTimeCam            timecam;
   MCalibrationRelTimeCam     relcam;
+  MCalibrationQECam          qecam;
   MCalibrationChargeCam      calcam;
   MCalibrationChargePINDiode pindiode;   
@@ -174 +132 @@
   histblind.SetSinglePheCut(600);
   //
-  // As long, as we don't have digital modules,
-  // we have to set the color of the pulser LED by hand
-  //
-  blindpix.SetColor(kCT1);
-  //    pindiode.SetColor(kCT1);
-  //
   // Get the previously created MPedestalCam into the new Parameter List 
   //
@@ -186 +138 @@
   plist.AddToList(&timecam);
   plist.AddToList(&relcam);
+  plist.AddToList(&qecam);
   plist.AddToList(&calcam);
   plist.AddToList(&histtime);
@@ -210 +163 @@
   MFillH fillblind("MHCalibrationChargeBlindPix", "MExtractedSignalBlindPixel");
   MFillH fillcam  ("MHCalibrationChargeCam"     , "MExtractedSignalCam");
+  
+  //
+  // As long, as we don't have digital modules,
+  // we have to set the color of the pulser LED by hand
+  //
+  calcalc.SetPulserColor(MCalibrationCam::kCT1);
   
   //
@@ -300 +259 @@
   MHCamera disp7  (geomcam, "Cal;FFactorConv",          "Conversion Factor to photons (F-Factor Method)");
   MHCamera disp8  (geomcam, "Cal;FFactorFFactor",       "Total F-Factor (F-Factor Method)");
-  MHCamera disp9  (geomcam, "Cal;BlindPixConv",         "Conversion Factor to photons (Blind Pixel Method)");
-  MHCamera disp10 (geomcam, "Cal;BlindPixFFactor",      "Total F-Factor (Blind Pixel Method)");
+  MHCamera disp9  (geomcam, "Cal;CascadesQEFFactor",    "Av. Quantum Efficiency (F-Factor Method)");
+  MHCamera disp10 (geomcam, "Cal;QEFFactor",            "Measured QE (F-Factor Method)");
   MHCamera disp11 (geomcam, "Cal;PINDiodeConv",         "Conversion Factor tp photons (PIN Diode Method)");
   MHCamera disp12 (geomcam, "Cal;PINDiodeFFactor",      "Total F-Factor (PIN Diode Method)");
@@ -351 +310 @@
   disp8.SetCamError(  calcam, 14);
   
-  // Blind Pixel Method
-  disp9.SetCamContent(calcam, 16);
-  disp9.SetCamError(  calcam, 17);
-  disp10.SetCamContent(calcam,18);
-  disp10.SetCamError(  calcam,19);
+  // Quantum efficiency
+  disp9.SetCamContent( qecam, 0 );
+  disp9.SetCamError(   qecam, 1 );
+  disp10.SetCamContent(qecam, 8);
+  disp10.SetCamError(  qecam, 9);
   
   // PIN Diode Method
@@ -411 +370 @@
   disp8.SetYTitle("\\sqrt{N_{Ph}}*\\sigma_{Charge}/\\mu_{Charge} [1] ");
   
-  disp9.SetYTitle("Conversion Factor [Phot/FADC Count]");
-  disp10.SetYTitle("\\sqrt{N_{Ph}}*\\sigma_{Charge}/\\mu_{Charge} [1]");
+  disp9.SetYTitle("Average QE Cascades [1]");
+  disp10.SetYTitle("Measured QE (F-Factor Method) [1]");
   
   disp11.SetYTitle("Conversion Factor [Phot/FADC Count]");
@@ -475 +434 @@
   
 
-  // Blind Pixel Method
-  TCanvas &c5 = display->AddTab("BlindPix");
+  // Quantum Efficiencies
+  TCanvas &c5 = display->AddTab("QE");
   c5.Divide(2, 4);