Changeset 4949 for trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc

Timestamp:

09/11/04 21:44:12 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc

@@ -50 +50 @@
 /////////////////////////////////////////////////////////////////////////////
 #include "MHCalibrationCam.h"
+#include "MHCalibrationPix.h"
 
 #include <TVirtualPad.h>
@@ -60 +61 @@
 #include "MLogManip.h"
 
+#include "MCalibrationIntensityCam.h"
+#include "MCalibrationCam.h"
 #include "MCalibrationPix.h"
-#include "MCalibrationCam.h"
-
-#include "MHCalibrationPix.h"
 
 #include "MBadPixelsPix.h"
@@ -80 +80 @@
 
 const Int_t   MHCalibrationCam::fgPulserFrequency = 500;
+const TString MHCalibrationCam::gsHistName   = "Hist";
+const TString MHCalibrationCam::gsHistTitle  = "";
+const TString MHCalibrationCam::gsHistXTitle = "";
+const TString MHCalibrationCam::gsHistYTitle = "Nr. events";
 // --------------------------------------------------------------------------
 //
@@ -96 +100 @@
 //
 MHCalibrationCam::MHCalibrationCam(const char *name, const char *title)
-    :  fColor(MCalibrationCam::kNONE), 
+    :  fHistName(gsHistName),fHistTitle(gsHistTitle),
+       fHistXTitle(gsHistXTitle),fHistYTitle(gsHistYTitle), 
+       fColor(MCalibrationCam::kNONE), 
        fBadPixels(NULL), fIntensCam(NULL), fCam(NULL), fGeom(NULL), fRunHeader(NULL)
 {
@@ -300 +306 @@
 
 
-const TArrayI &MHCalibrationCam::GetRunNumbers() const 
-{
-  return fRunNumbers;
-}
-
 // --------------------------------------------------------------------------
 //
@@ -427 +428 @@
 }
 
-
-Bool_t MHCalibrationCam::SetupHists(const MParList *pList)
-{
-  return kTRUE;
-}
 
 // --------------------------------------------------------------------------
@@ -570 +566 @@
 }
 
-
-
-Bool_t MHCalibrationCam::ReInitHists(MParList *pList)
-{
-  return kTRUE;
-}
+//--------------------------------------------------------------------------------------
+//
+// Initializes the High Gain Arrays:
+//
+// - Expand fHiGainArrays to npixels
+// - Expand fAverageHiGainAreas to nareas
+// - Expand fAverageHiGainSectors to nsectors
+// 
+// - For every entry in the expanded arrays: 
+//   * call new MHCalibrationPix(fHistName,fHistTitle)
+//   * SetNbins(fNbins), SetFirst(fFirst), SetLast(fLast)
+//   * Set fHistName and fHistTitle in the Histogram of MHCalibrationPix
+//   * Set X-axis and Y-axis titles with fHistXTitle and fHistYTitle
+//   * Call InitHists
+//
+void MHCalibrationCam::InitHiGainArrays(const Int_t npixels, const Int_t nareas, const Int_t nsectors)
+{
+
+  if (fHiGainArray->GetEntries()==0)
+  {
+      fHiGainArray->Expand(npixels);
+      for (Int_t i=0; i<npixels; i++)
+      {
+        (*fHiGainArray)[i] = new MHCalibrationPix(Form("%s%s",fHistName.Data(),"HiGainPix"),
+                                                  Form("%s%s",fHistTitle.Data()," High Gain Pixel"));
+        
+        MHCalibrationPix &pix = (*this)[i];          
+        pix.SetNbins(fNbins);
+        pix.SetFirst(fFirst);
+        pix.SetLast (fLast);
+        
+        TH1F *h =  pix.GetHGausHist();
+        
+        h->SetName (Form("%s%s%s","H",fHistName.Data(),"HiGainPix"));
+        h->SetTitle(Form("%s%s",fHistTitle.Data()," High Gain Pixel "));  
+        h->SetXTitle(fHistXTitle.Data());
+        h->SetYTitle(fHistYTitle.Data());
+        
+        InitHists(pix,(*fBadPixels)[i],i);
+      }
+  }
+
+  if (fAverageHiGainAreas->GetEntries()==0)
+  {
+    fAverageHiGainAreas->Expand(nareas);
+    
+    for (Int_t j=0; j<nareas; j++)
+      {
+        (*fAverageHiGainAreas)[j] = new MHCalibrationPix(Form("%s%s",fHistName.Data(),"HiGainArea"),
+                                                         Form("%s%s",fHistTitle.Data()," High Gain Area Idx "));
+        
+        MHCalibrationPix &pix = GetAverageHiGainArea(j);
+        
+        pix.SetNbins(fNbins*(Int_t)TMath::Sqrt((Float_t)npixels/nareas));
+        pix.SetFirst(fFirst);
+        pix.SetLast (fLast);
+        
+        TH1F *h =  pix.GetHGausHist();
+
+        h->SetName (Form("%s%s%s","H",fHistName.Data(),"HiGainArea"));
+        h->SetXTitle(fHistXTitle.Data());
+        h->SetYTitle(fHistYTitle.Data());
+
+        if (fGeom && fGeom->InheritsFrom("MGeomCamMagic"))
+          {
+            h->SetTitle(Form("%s%s%s%s",fHistTitle.Data()," averaged on event-by-event basis ",
+                             j==0 ? "Inner Pixels " : "Outer Pixels ","High Gain Runs: "));
+            pix.InitBins();
+            pix.SetEventFrequency(fPulserFrequency);
+          }
+        else
+          {
+            h->SetTitle(Form("%s%s",fHistTitle.Data()," averaged on event-by-event basis High Gain Area Idx ")); 
+            InitHists(pix,fIntensCam ? fIntensCam->GetAverageBadArea(j) : fCam->GetAverageBadArea(j),j);
+          }
+      }
+  }
+
+  if (fAverageHiGainSectors->GetEntries()==0)
+    {
+      fAverageHiGainSectors->Expand(nsectors);
+
+      for (Int_t j=0; j<nsectors; j++)
+        {
+          (*fAverageHiGainSectors)[j] = new MHCalibrationPix(Form("%s%s",fHistName.Data(),"HiGainSector"),
+                                                             Form("%s%s",fHistTitle.Data()," High Gain Sector "));
+          MHCalibrationPix &pix = GetAverageHiGainSector(j);
+
+          pix.SetNbins(fNbins*(Int_t)TMath::Sqrt((Float_t)npixels/nsectors));
+          pix.SetFirst(fFirst);
+          pix.SetLast (fLast);
+          
+          TH1F *h =  pix.GetHGausHist();
+
+          h->SetName (Form("%s%s%s","H",fHistName.Data(),"HiGainSector"));
+          h->SetTitle(Form("%s%s",fHistTitle.Data()," High Gain Sector "));  
+          h->SetXTitle(fHistXTitle.Data());
+          h->SetYTitle(fHistYTitle.Data());
+
+          InitHists(pix,fIntensCam ? fIntensCam->GetAverageBadSector(j) : fCam->GetAverageBadSector(j),j);
+      }
+  }
+}
+
+//--------------------------------------------------------------------------------------
+//
+// Return, if IsLoGain() is kFALSE 
+//
+// Initializes the Low Gain Arrays:
+//
+// - Expand fLoGainArrays to npixels
+// - Expand fAverageLoGainAreas to nareas
+// - Expand fAverageLoGainSectors to nsectors
+// 
+// - For every entry in the expanded arrays: 
+//   * call new MHCalibrationPix(fHistName,fHistTitle)
+//   * SetBins(fNbins), SetFirst(fFirst), SetLast(fLast)
+//   * Set fHistName and fHistTitle in the Histogram of MHCalibrationPix
+//   * Set X-axis and Y-axis titles with fHistXTitle and fHistYTitle
+//   * Call InitHists
+//
+void MHCalibrationCam::InitLoGainArrays(const Int_t npixels, const Int_t nareas, const Int_t nsectors)
+{
+
+  if (!IsLoGain())
+    return;
+
+  if (fLoGainArray->GetEntries()==0)
+  {
+      fLoGainArray->Expand(npixels);
+      for (Int_t i=0; i<npixels; i++)
+      {
+        (*fLoGainArray)[i] = new MHCalibrationPix(Form("%s%s",fHistName.Data(),"LoGainPix"),
+                                                  Form("%s%s",fHistTitle.Data()," Low Gain Pixel"));
+
+          MHCalibrationPix &pix = (*this)[i];          
+          pix.SetNbins(fNbins);
+          pix.SetFirst(fFirst);
+          pix.SetLast (fLast);
+
+          TH1F *h =  pix.GetHGausHist();
+          
+          h->SetName (Form("%s%s%s","H",fHistName.Data(),"LoGainPix"));
+          h->SetTitle(Form("%s%s",fHistTitle.Data()," Low Gain Pixel "));  
+          h->SetXTitle(fHistXTitle.Data());
+          h->SetYTitle(fHistYTitle.Data());
+
+          InitHists(pix,(*fBadPixels)[i],i);
+      }
+  }
+
+  if (fAverageLoGainAreas->GetEntries()==0)
+  {
+    fAverageLoGainAreas->Expand(nareas);
+    
+    for (Int_t j=0; j<nareas; j++)
+      {
+        (*fAverageLoGainAreas)[j] = new MHCalibrationPix(Form("%s%s",fHistName.Data(),"LoGainArea"),
+                                                         Form("%s%s",fHistTitle.Data()," Low Gain Area Idx "));
+        
+        MHCalibrationPix &pix = GetAverageLoGainArea(j);
+        
+        pix.SetNbins(fNbins*(Int_t)TMath::Sqrt((Float_t)npixels/nareas));
+        pix.SetFirst(fFirst);
+        pix.SetLast (fLast);
+        
+        TH1F *h =  pix.GetHGausHist();
+
+        h->SetName (Form("%s%s%s","H",fHistName.Data(),"LoGainArea"));
+        h->SetXTitle(fHistXTitle.Data());
+        h->SetYTitle(fHistYTitle.Data());
+
+        if (fGeom && fGeom->InheritsFrom("MGeomCamMagic"))
+          {
+            h->SetTitle(Form("%s%s%s%s",fHistTitle.Data()," averaged on event-by-event basis ",
+                             j==0 ? "Inner Pixels " : "Outer Pixels ","Low Gain Runs: "));
+            pix.InitBins();
+            pix.SetEventFrequency(fPulserFrequency);
+          }
+        else
+          {
+            h->SetTitle(Form("%s%s",fHistTitle.Data()," averaged on event-by-event basis Low Gain Area Idx ")); 
+            InitHists(pix,fIntensCam ? fIntensCam->GetAverageBadArea(j) : fCam->GetAverageBadArea(j),j);
+          }
+      }
+  }
+
+  if (fAverageLoGainSectors->GetEntries()==0)
+    {
+      fAverageLoGainSectors->Expand(nsectors);
+
+      for (Int_t j=0; j<nsectors; j++)
+        {
+          (*fAverageLoGainSectors)[j] = new MHCalibrationPix(Form("%s%s",fHistName.Data(),"LoGainSector"),
+                                                             Form("%s%s",fHistTitle.Data()," Low Gain Sector "));
+          MHCalibrationPix &pix = GetAverageLoGainSector(j);
+
+          pix.SetNbins(fNbins*(Int_t)TMath::Sqrt((Float_t)npixels/nsectors));
+          pix.SetFirst(fFirst);
+          pix.SetLast (fLast);
+          
+          TH1F *h =  pix.GetHGausHist();
+
+          h->SetName (Form("%s%s%s","H",fHistName.Data(),"LoGainSector"));
+          h->SetTitle(Form("%s%s",fHistTitle.Data()," Low Gain Sector "));  
+          h->SetXTitle(fHistXTitle.Data());
+          h->SetYTitle(fHistYTitle.Data());
+
+          InitHists(pix,fIntensCam ? fIntensCam->GetAverageBadSector(j) : fCam->GetAverageBadSector(j),j);
+      }
+  }
+}
+
+
 
 //--------------------------------------------------------------------------------
@@ -1099 +1303 @@
     if (IsEnvDefined(env, prefix, "Oscillations", print))
     {
-        SetDebug(GetEnvValue(env, prefix, "Oscillations", IsOscillations()));
+        SetOscillations(GetEnvValue(env, prefix, "Oscillations", IsOscillations()));
         rc = kTRUE;
     }
 
+    if (IsEnvDefined(env, prefix, "Nbins", print))
+    {
+        SetNbins(GetEnvValue(env, prefix, "Nbins", fNbins));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "First", print))
+    {
+        SetFirst(GetEnvValue(env, prefix, "First", fFirst));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "Last", print))
+    {
+        SetLast(GetEnvValue(env, prefix, "Last", fLast));
+        rc = kTRUE;
+    }
 
     return rc;