Changeset 3903 for trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimeCam.cc

Timestamp:

04/30/04 11:05:50 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimeCam.cc (modified) (12 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimeCam.cc

@@ -108 +108 @@
 
 using namespace std;
+
+const UInt_t MHCalibrationRelTimeCam::fgReferencePixel = 1;
 // --------------------------------------------------------------------------
 //
-// Default Constructor. 
+// Default Constructor.
+//
+// Sets: 
+// - fReferencePixel to fgReferencePixel
 //
 MHCalibrationRelTimeCam::MHCalibrationRelTimeCam(const char *name, const char *title) 
@@ -117 +122 @@
   fName  = name  ? name  : "MHCalibrationRelTimeCam";
   fTitle = title ? title : "Histogram class for the relative time calibration of the camera";
-  
+
+  SetReferencePixel();
 }
 
@@ -183 +189 @@
       {
           (*fHiGainArray)[i] = new MHCalibrationRelTimePix("MHCalibrationRelTimePixHiGain",
-                                                          "Relative Arrival Time High Gain Pixel");
+                                                          "Rel. Arr. Time Hi-Gain Pixel ");
           InitHists((*this)[i],(*fBadPixels)[i],i);
       }
   }
-
+  
+  if (fRunHeader)
+    for (Int_t i=0; i<npixels; i++)
+      {
+        TH1F *h = (*this)[i].GetHGausHist();
+          h->SetTitle( Form("%s%i%s", h->GetTitle(),fRunNumbers[fRunNumbers.GetSize()-1]," "));
+      }
+  
   if (fLoGainArray->GetEntries()==0)
   {
@@ -194 +207 @@
       {
           (*fLoGainArray)[i] = new MHCalibrationRelTimePix("MHCalibrationRelTimePixLoGain",
-                                                          "Relative Arrival Time Low Gain Pixel");
+                                                          "Rel. Arr. Time Lo-Gain Pixel ");
           InitHists((*this)(i),(*fBadPixels)[i],i);
       }
   }
+
+  if (fRunHeader)
+    for (Int_t i=0; i<npixels; i++)
+      {
+        TH1F *h = (*this)(i).GetHGausHist();
+        h->SetTitle( Form("%s%i%s", h->GetTitle(),fRunNumbers[fRunNumbers.GetSize()-1]," "));
+      }
 
   if (fAverageHiGainAreas->GetEntries()==0)
@@ -207 +227 @@
         (*fAverageHiGainAreas)[j] = 
           new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageAreaHiGain",
-                                      "Average Relative Arrival Times High Gain Area Idx ");
+                                      "Average Rel. Arr. Times Hi-Gain Area Idx ");
 
         GetAverageHiGainArea(j).GetHGausHist()->SetTitle("Rel. Arr. Times HiGain Area Idx ");
@@ -216 +236 @@
   }
 
+  if (fRunHeader)
+    for (Int_t j=0; j<nareas; j++)
+      {
+        TH1F *h = GetAverageHiGainArea(j).GetHGausHist();
+        h->SetTitle( Form("%s%i%s", h->GetTitle(),fRunNumbers[fRunNumbers.GetSize()-1]," "));
+      }
+  
   if (fAverageLoGainAreas->GetEntries()==0)
   {
@@ -224 +251 @@
         (*fAverageLoGainAreas)[j] = 
           new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageAreaLoGain",
-                                      "Average Relative Arrival Times Low Gain Area Idx ");
+                                      "Average Rel. Arr. Times Lo-Gain Area Idx ");
 
         GetAverageLoGainArea(j).GetHGausHist()->SetTitle("Rel. Arr. Times LoGain Area Idx ");
@@ -232 +259 @@
       }
   }
-
+  
+  if (fRunHeader)
+    for (Int_t j=0; j<nareas; j++)
+        {
+          TH1F *h = GetAverageLoGainArea(j).GetHGausHist();
+          h->SetTitle( Form("%s%i%s", h->GetTitle(),fRunNumbers[fRunNumbers.GetSize()-1]," "));
+        }
+  
   if (fAverageHiGainSectors->GetEntries()==0)
-  {
+    {
       fAverageHiGainSectors->Expand(nsectors);
 
       for (Int_t j=0; j<nsectors; j++)
-      {
+        {
           (*fAverageHiGainSectors)[j] = 
             new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageSectorHiGain",
-                                        "Average Relative Arrival Times High Gain Sector ");
-
+                                        "Average Rel. Arr. Times Hi-Gain Sector ");
+          
           GetAverageHiGainSector(j).GetHGausHist()->SetTitle("Rel. Arr. Times HiGain Sector ");
           GetAverageHiGainSector(j).SetNbins(fAverageNbins);
-
+          
           InitHists(GetAverageHiGainSector(j),fCam->GetAverageBadSector(j),j);
       }
   }
+  
+  if (fRunHeader)
+    for (Int_t j=0; j<nsectors; j++)
+      {
+        TH1F *h = GetAverageHiGainSector(j).GetHGausHist();
+        h->SetTitle( Form("%s%i%s", h->GetTitle(),fRunNumbers[fRunNumbers.GetSize()-1]," "));
+      }
 
   if (fAverageLoGainSectors->GetEntries()==0)
-  {
+    {
       fAverageLoGainSectors->Expand(nsectors);
-
+      
       for (Int_t j=0; j<nsectors; j++)
-      {
+        {
           (*fAverageLoGainSectors)[j] = 
             new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageSectorLoGain",
-                                        "Average Relative Arrival Times Low Gain Sector ");
-
+                                        "Average Rel. Arr. Times Lo-Gain Sector ");
+          
           GetAverageLoGainSector(j).GetHGausHist()->SetTitle("Rel. Arr. Times LoGain Sector ");
           GetAverageLoGainSector(j).SetNbins(fAverageNbins);
-
+          
           InitHists(GetAverageLoGainSector(j),fCam->GetAverageBadSector(j),j);
           
-      }
-  }
-
+        }
+    }
+
+  if (fRunHeader)
+    for (Int_t j=0; j<nsectors; j++)
+      {
+        TH1F *h = GetAverageLoGainSector(j).GetHGausHist();
+        h->SetTitle( Form("%s%i%s", h->GetTitle(),fRunNumbers[fRunNumbers.GetSize()-1]," "));
+      }
+  
   return kTRUE;
 }
@@ -309 +357 @@
   memset(sumsectorhi, 0, nsectors*sizeof(Float_t));
   memset(sumsectorlo, 0, nsectors*sizeof(Float_t));
-  
-  const MArrivalTimePix &refpix = (*arrtime)[1];
+  memset(numareahi,   0, nareas * sizeof(Float_t));
+  memset(numarealo,   0, nareas * sizeof(Float_t));
+  memset(numsectorhi, 0, nsectors*sizeof(Float_t));
+  memset(numsectorlo, 0, nsectors*sizeof(Float_t));
+  
+  const MArrivalTimePix &refpix = (*arrtime)[fReferencePixel];
   const Float_t reftime = refpix.IsLoGainUsed() 
     ? refpix.GetArrivalTimeLoGain() : refpix.GetArrivalTimeHiGain();
@@ -324 +376 @@
 
       const MArrivalTimePix &pix = (*arrtime)[i];
-      const Float_t time    = pix.IsLoGainUsed() 
-        ? pix.GetArrivalTimeLoGain() : pix.GetArrivalTimeHiGain();  
-      const Float_t reltime = time - reftime;
-      
       const Int_t aidx   = (*fGeom)[i].GetAidx();
       const Int_t sector = (*fGeom)[i].GetSector();
@@ -333 +381 @@
       if (pix.IsLoGainUsed())
         { 
-
+          const Float_t reltime = pix.GetArrivalTimeLoGain() - reftime;
           histlo.FillHistAndArray(reltime);
           sumarealo  [aidx]   += reltime;
@@ -342 +390 @@
       else
         {
-
+          const Float_t reltime = pix.GetArrivalTimeHiGain() - reftime;
           histhi.FillHistAndArray(reltime) ;
           sumareahi  [aidx]   += reltime;