Changeset 2765

Timestamp:

01/09/04 22:57:10 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mcalib/MCalibrationCalc.cc (modified) (2 diffs)
• mcalib/MCalibrationCam.cc (modified) (2 diffs)
• mcalib/MCalibrationCam.h (modified) (3 diffs)
• mcalib/MCalibrationPix.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -5 +5 @@
                                                  -*-*- END OF LINE -*-*-
  2004/01/09: Markus Gaug
+
+   * mcalib/MCalibrationCam.[h,cc]
+     - take out 1-dim histos because we can use 
+       the projections in MHCamera (s.b.)
 
    * mcalib/MHCalibrationBlindPixel.h

trunk/MagicSoft/Mars/mcalib/MCalibrationCalc.cc

@@ -239 +239 @@
     for (UInt_t i=0;i<cam->GetNumPixels();i++)
       {
+
         MCalibrationPix &pix = (*fCalibrations)[i];
         pix.DefinePixId(i);
+        MHCalibrationPixel *hist = pix.GetHist();
+
+        hist->SetTimeFitRangesHiGain(fSignals->GetFirstUsedSliceHiGain(),
+                                    fSignals->GetLastUsedSliceHiGain());
+        hist->SetTimeFitRangesLoGain(fSignals->GetFirstUsedSliceLoGain(),
+                                    fSignals->GetLastUsedSliceLoGain());
+        
       }
     
@@ -453 +461 @@
       MCalibrationPix &pix = (*fCalibrations)[pixid];
 
+      //
+      // get the pedestals
+      //
       const Float_t ped    = (*fPedestals)[pixid].GetPedestal() * fNumHiGainSamples;
-      const Float_t prms   = (*fPedestals)[pixid].GetPedestalRms() * TMath::Sqrt((float)fNumHiGainSamples);
-
+      const Float_t prms   = (*fPedestals)[pixid].GetPedestalRms() * fSqrtHiGainSamples;
+
+      //
+      // set them in the calibration camera
+      //
       pix.SetPedestal(ped,prms);
 
+      //
+      // perform the Gauss fits to the charges
+      //
       pix.FitCharge();
 
+      //
+      // Fill 1-dim histograms with the results
+      //
+      fCalibrations->Fill1DimCharges(pix.GetCharge(),pix.GetErrCharge(),
+                                     pix.GetSigmaCharge(),pix.GetChargeProb());
+
+      //
+      // Perform the Gauss fits to the arrival times
+      //
       if (TESTBIT(fFlags,kUseTimeFits))
         pix.FitTime();

trunk/MagicSoft/Mars/mcalib/MCalibrationCam.cc

@@ -83 +83 @@
     fBlindPixel = new MCalibrationBlindPix();
     fPINDiode   = new MCalibrationPINDiode();
+
 }
 
@@ -337 +338 @@
       break;
     case 7:
-      val = (*this)[idx].GetTimeChiSquare();
+      val = (*this)[idx].GetTimeProb();
       break;
     case 8:

trunk/MagicSoft/Mars/mcalib/MCalibrationCam.h

@@ -57 +57 @@
   Float_t fMeanPhotErrOutsidePlexiglass;      //  The uncertainty about the number of photons in an INNER PIXEL  
 
-  TH1F* fHChargeAllPixels;                     //! 
-  TH1F* fHChargeErrAllPixels;                  //! 
-  TH1F* fHChargeSigmaAllPixels;                //! 
-  TH1F* fHChargeProbAllPixels;                 //! 
-
-  TH1F* fHTimeAllPixels;                       //! 
-  TH1F* fHTimeSigmaAllPixels;                  //! 
-  TH1F* fHTimeProbAllPixels;                   //! 
-
-  TH1F* fHPedAllPixels;                        //! 
-  TH1F* fHPedRMSAllPixels;                     //! 
-  TH1F* fHSigmaPerChargeAllPixels;             //! 
-  TH1F* fHPhEAllPixels;                        //! 
-
-  TH1F* fHConvPhEperFADCAllPixels;             //! 
-  TH1F* fHConvPhperFADCAllPixels;              //! 
-
-  TH1D* fOffsets;                              //! 
-  TH1D* fSlopes;                               //! 
+  TH1D* fOffsets;                             //! 
+  TH1D* fSlopes;                              //! 
   
-  TH2D* fOffvsSlope;                           //! 
+  TH2D* fOffvsSlope;                          //! 
 
 public:
@@ -95 +78 @@
   
   void InitSize(const Int_t i);
-  Int_t GetSize() const;
+  Int_t GetSize()                    const;
 
-  UInt_t GetNumPixels() const { return fNumPixels; }
+  UInt_t GetNumPixels()              const { return fNumPixels; }
 
   Bool_t IsPixelUsed(Int_t idx)      const;
@@ -117 +100 @@
   MCalibrationPINDiode *GetPINDiode()                const { return fPINDiode;    }
 
-  void SetColor(CalibrationColor_t color)    { fColor = color; }
+  void SetColor(CalibrationColor_t color)                     { fColor = color; }
 
   Bool_t GetConversionFactorFFactor(Int_t ipx, Float_t &mean, Float_t &err, Float_t &sigma);

trunk/MagicSoft/Mars/mcalib/MCalibrationPix.cc

@@ -389 +389 @@
 {
 
-  if (fTime < fHist->GetTimeLowerFitRange()+1.)
+  Float_t lowerrange;
+  Float_t upperrange;
+
+  if (fHiGainSaturation)
+    {
+      lowerrange = (Float_t)fHist->GetTimeLowerFitRangeLoGain()+1.;
+      upperrange = (Float_t)fHist->GetTimeUpperFitRangeLoGain()+1.;
+    }
+  else
+    {
+      lowerrange = (Float_t)fHist->GetTimeLowerFitRangeHiGain()+1.;
+      upperrange = (Float_t)fHist->GetTimeUpperFitRangeHiGain()+1.;
+    }
+
+
+  if (fTime < lowerrange)
     {
       *fLog << warn 
             << "WARNING: Mean Fitted Time inside or smaller than first used FADC slice in Pixel " 
-            << fPixId << endl;
-      return kFALSE;
-    }
-
-  if (fTime < fHist->GetTimeUpperFitRange()-1.)
+            << fPixId << " time: " << fTime << " Range: " << lowerrange << endl;
+      return kFALSE;
+    }
+
+  if (fTime > upperrange)
     {
       *fLog << warn 
             << "WARNING: Mean Fitted Time inside or greater than last used FADC slice in Pixel " 
-            << fPixId << endl;
+            << fPixId << " time: " << fTime << " Range: " << upperrange << endl;
       return kFALSE;
     }
@@ -484 +499 @@
   fSigmaTime     = fHist->GetTimeSigma();
   fTimeChiSquare = fHist->GetTimeChiSquare();
+  fTimeProb      = fHist->GetTimeProb();
 
   if (!CheckTimeFitValidity())