Changeset 5729

Timestamp:

01/07/05 19:56:06 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• manalysis/MEnergyEstParam.h (modified) (1 diff)
• mcalib/MCalibrateData.cc (modified) (11 diffs)
• mcalib/MCalibrateData.h (modified) (2 diffs)

trunk/MagicSoft/Mars/manalysis/MEnergyEstParam.h

@@ -58 +58 @@
     void SetCoeffB(const TArrayD &arr);
 
+    void SetVariables(const TArrayD &arr) { SetCoeff(arr); }
+
     Double_t GetCoeff(Int_t i) { return i<fA.GetSize()? fA[i] : fB[i-fA.GetSize()]; }
 

trunk/MagicSoft/Mars/mcalib/MCalibrateData.cc

@@ -66 +66 @@
 //     SetPedestalFlag(MCalibrateData::kEvent)  (calibration is done for each event)
 //
-//   By calling SetNamePedestalCam() and/or SetNamePedPhotCam() you
-//   can control the name of the MPedestalCam and/or MPedPhotCam container which is used.
+//   By calling AddPedestal() you can control the name of the
+//   MPedestalCam and/or MPedPhotCam container which is used.
 //
 //   Assume you want to calibrate "MPedestalCam" once and "MPedestalFromLoGain" [MPedestalCam]
@@ -76 +76 @@
 //     MCalibrateData cal2;
 //     cal2.SetCalibrationMode(MCalibrateData::kSkip);
-//     cal2.SetNamePedestalCam("MPedestalFromLoGain");
-//     cal2.SetNamePedPhotCam("MPedPhotFromLoGain")
+//     cal2.AddPedestal("FromLoGain");
 //     cal2.SetPedestalFlag(MCalibrateData::kEvent);
 //
@@ -133 +132 @@
 using namespace std;
 
-const TString MCalibrateData::fgNamePedestalCam  = "MPedestalCam";
-const TString MCalibrateData::fgNamePedPhotCam = "MPedPhotCam";
-
 // --------------------------------------------------------------------------
 //
@@ -145 +141 @@
 // - fCalibrationMode to kDefault
 // - fPedestalFlag to kNo
-// - fNamePedADCRunContainer    to "MPedestalCam"  
-// - fNamePedPhotRunContainer   to "MPedPhotCam"    
 //
 MCalibrateData::MCalibrateData(CalibrationMode_t calmode,const char *name, const char *title) 
-    : fGeomCam(NULL),   fPedestalMean(NULL), fPedestalRms(NULL),
-      fBadPixels(NULL), fCalibrations(NULL), fQEs(NULL), fSignals(NULL), 
-      fPedPhot(NULL),   fCerPhotEvt(NULL), fPedestalFlag(kNo)
+    : fGeomCam(NULL),   fBadPixels(NULL), fCalibrations(NULL),
+    fQEs(NULL), fSignals(NULL), fCerPhotEvt(NULL), fPedestalFlag(kNo)
 {
 
@@ -159 +152 @@
   SetCalibrationMode(calmode);
 
-  SetNamePedestalCam();
-  SetNamePedPhotCam();
+  fNamesPedestal.SetOwner();
+}
+
+void MCalibrateData::AddPedestal(const char *name)
+{
+    TString ped(name);
+    TString pho(name);
+    ped.Prepend("MPedestal");
+    pho.Prepend("MPedPhot");
+
+    fNamesPedestal.Add(new TNamed(ped, pho));
+}
+
+void MCalibrateData::AddPedestal(const char *pedestal, const char *pedphot)
+{
+    fNamesPedestal.Add(new TNamed(pedestal, pedphot));
 }
 
@@ -226 +233 @@
     }
 
-    //fPedestalMean = 0;
-    //fPedestalRms  = 0;
-    fPedPhot  = 0;
+    if (fNamesPedestal.GetSize()>0 && fPedestalFlag==kNo)
+    {
+        *fLog << warn << "Pedestal list contains entries, but mode is set to kNo... setting to kEvent." << endl;
+        fPedestalFlag = kEvent;
+    }
+
     if (fPedestalFlag)
     {
-        if (!fPedestalMean)
-        {
-            fPedestalMean = (MPedestalCam*)pList->FindObject(AddSerialNumber(fNamePedestalCam), "MPedestalCam");
-            if (!fPedestalMean)
+        if (fNamesPedestal.GetSize()==0)
+        {
+            *fLog << inf << "No container names specified... using default: MPedestalCam and MPedPhotCam." << endl;
+            AddPedestal();
+        }
+
+        fPedestalCams.Clear();
+        fPedPhotCams.Clear();
+
+        TIter Next(&fNamesPedestal);
+        TObject *o=0;
+        while ((o=Next()))
+        {
+            TObject *pedcam = pList->FindObject(AddSerialNumber(o->GetName()), "MPedestalCam");
+            if (!pedcam)
             {
-                *fLog << err << AddSerialNumber(fNamePedestalCam) << " [MPedestalCam] not found ... aborting" << endl;
+                *fLog << err << AddSerialNumber(o->GetName()) << " [MPedestalCam] not found ... aborting" << endl;
                 return kFALSE;
             }
-        }
-
-        if (!fPedestalRms)
-        {
-            fPedestalRms = (MPedestalCam*)pList->FindObject(AddSerialNumber(fNamePedestalCam), "MPedestalCam");
-            if (!fPedestalRms)
-            {
-                *fLog << err << AddSerialNumber(fNamePedestalCam) << " [MPedestalCam] not found ... aborting" << endl;
+            TObject *pedphot = pList->FindCreateObj("MPedPhotCam", AddSerialNumber(o->GetTitle()));
+            if (!pedphot)
                 return kFALSE;
-            }
-        }
-
-        fPedPhot = (MPedPhotCam*)pList->FindCreateObj("MPedPhotCam", AddSerialNumber(fNamePedPhotCam));
-        if (!fPedPhot)
-            return kFALSE;
+
+            fPedestalCams.Add(pedcam);
+            fPedPhotCams.Add(pedphot);
+        }
     }
 
@@ -281 +294 @@
 
     // Sizes might have changed
-    if (fPedestalFlag && (Int_t)fPedestalMean->GetSize() != fSignals->GetSize())
-    {
-        *fLog << err << "Size mismatch of MPedestalCam for Mean and MCalibrationCam... abort." << endl;
-        return kFALSE;
-    }
-
-    if (fPedestalFlag && (Int_t)fPedestalRms->GetSize() != fSignals->GetSize())
-    {
-        *fLog << err << "Size mismatch of MPedestalCam for RMS and MCalibrationCam... abort." << endl;
-        return kFALSE;
+    if (fPedestalFlag)
+    {
+        TIter Next(&fPedestalCams);
+        MPedestalCam *cam=0;
+        while ((cam=(MPedestalCam*)Next()))
+        if ((Int_t)cam->GetSize() != fSignals->GetSize())
+        {
+            *fLog << err << "Size mismatch of " << cam->GetDescriptor() << " and MCalibrationCam... abort." << endl;
+            return kFALSE;
+        }
     }
 
@@ -523 +536 @@
         if (pedestal)
         {
+            /*
             // pedestals/(used FADC slices)   in [ADC] counts
             const Float_t pedes  = (*fPedestalMean)[pixidx].GetPedestal()   * slices;
@@ -534 +548 @@
 
             (*fPedPhot)[pixidx].Set(pedphot, pedphotrms);
+            */
+            TIter NextPed(&fPedestalCams);
+            TIter NextPhot(&fPedPhotCams);
+
+            MPedestalCam *pedestal = 0;
+            MPedPhotCam  *pedphot  = 0;
+
+            while ((pedestal=(MPedestalCam*)NextPed()) &&
+                   (pedphot =(MPedPhotCam*)NextPhot()))
+            {
+                // pedestals/(used FADC slices)   in [ADC] counts
+                const Float_t pedes  = (*pedestal)[pixidx].GetPedestal()    * slices;
+                const Float_t pedrms = (*pedestal)[pixidx].GetPedestalRms() * sqrtslices;
+
+                // pedestals/(used FADC slices)   in [number of photons]
+                const Float_t mean = pedes  * calibConv;
+                const Float_t rms  = pedrms * calibConv;
+
+                (*pedphot)[pixidx].Set(mean, rms);
+                pedphot->SetReadyToSave();
+            }
         }
     }
@@ -542 +577 @@
         return kCONTINUE;
     }
-
-    if (pedestal)
-        fPedPhot->SetReadyToSave();
 
     if (data)
@@ -610 +642 @@
     }
 
-    if (fNamePedestalCam != fgNamePedestalCam)
-    {
-        out << "   " << GetUniqueName() << ".SetNamePedestalCam(";
-        out << fNamePedestalCam << ");" << endl;
-    }
-
-    if (fNamePedPhotCam  != fgNamePedPhotCam)
-    {
-        out << "   " << GetUniqueName() << ".SetNamePedPhotCam(";
-        out << fNamePedPhotCam << ");" << endl;
+    TIter Next(&fNamesPedestal);
+    TObject *o=0;
+    while ((o=Next()))
+    {
+        out << "   " << GetUniqueName() << ".AddPedestal(\"";
+        out << o->GetName() << "\", \"" << o->GetTitle() << "\");" << endl;
     }
 }

trunk/MagicSoft/Mars/mcalib/MCalibrateData.h

@@ -34 +34 @@
 {
 private:
-  static const TString fgNamePedestalCam;    //! "MPedestalCam"
-  static const TString fgNamePedPhotCam;     //! "MPedPhotCam"
+  MGeomCam              *fGeomCam;       //! Camera geometry container
+  MBadPixelsCam         *fBadPixels;     //! Bad Pixels information
+  MCalibrationChargeCam *fCalibrations;  //! Calibration constants
+  MCalibrationQECam     *fQEs;           //! Quantum efficiencies
+  MExtractedSignalCam   *fSignals;       //! Integrated charge in FADCs counts
+  MCerPhotEvt           *fCerPhotEvt;    //! Cerenkov Photon Event used for calculation
   
-  MGeomCam              *fGeomCam;          //! Camera geometry container
-  MPedestalCam          *fPedestalMean;     //! Pedestals/slice [ADC counts]
-  MPedestalCam          *fPedestalRms;      //! Pedestals/slice [ADC counts]
-  MBadPixelsCam         *fBadPixels;        //! Bad Pixels information
-  MCalibrationChargeCam *fCalibrations;     //! Calibration constants
-  MCalibrationQECam     *fQEs;              //! Quantum efficiencies
-  MExtractedSignalCam   *fSignals;          //! Integrated charge in FADCs counts
+  UShort_t fCalibrationMode;             // Flag defining the calibration mode (CalibrationMode_t)
+  Byte_t   fPedestalFlag;                // Flags defining to calibrate the pedestal each event or each run
   
-  MPedPhotCam           *fPedPhot;          //! Pedestals/(used slices) [photons]
-  MCerPhotEvt           *fCerPhotEvt;       //! Cerenkov Photon Event used for calculation
-  
-  UShort_t fCalibrationMode;                // Flag defining the calibration mode (CalibrationMode_t)
-  Byte_t   fPedestalFlag;                   // Flags defining to calibrate the pedestal each event or each run
-  
-  TString  fNamePedestalCam;            // name of fPedestal
-  TString  fNamePedPhotCam;           // name of fPedPhot
+  TList fNamesPedestal;                  // Names of input and output pedestal conatainer
+  TList fPedestalCams;                   //! List of pointers to input MPedestalCam
+  TList fPedPhotCams;                    //! List of pointers to corresponding output MPedPhotCam
 
   Int_t Calibrate(Bool_t data, Bool_t pedestal) const;
@@ -95 +89 @@
   
   void   SetCalibrationMode ( CalibrationMode_t calmode=kDefault ) { fCalibrationMode=calmode; }
-  
-  void   SetNamePedestalCam(const char *name=fgNamePedestalCam.Data())
-  {
-      fNamePedestalCam  = name;
-  }
-  void SetPedestalCamMean(MPedestalCam *cam) { fPedestalMean = cam; }
-  void SetPedestalCamRms(MPedestalCam *cam)  { fPedestalRms  = cam; }
 
-  void   SetNamePedPhotCam(const char *name=fgNamePedPhotCam.Data())
-  {
-      fNamePedPhotCam   = name;
-  }
-  
+  void   AddPedestal(const char *name="Cam");
+  void   AddPedestal(const char *pedestal, const char *pedphot);
+
   ClassDef(MCalibrateData, 0)   // Task to calibrate FADC counts into Cherenkov photons 
 };