Changeset 6792 for trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.cc

Timestamp:

03/08/05 15:14:57 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.cc (modified) (9 diffs)

trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.cc

@@ -119 +119 @@
 using namespace std;
 
-const Int_t   MHPedestalCam::fgNbins      = 200;
-const Axis_t  MHPedestalCam::fgFirst      = -49.5;
-const Axis_t  MHPedestalCam::fgLast       = 150.5;
+const Int_t   MHPedestalCam::fgNbins      =  50;
+const Axis_t  MHPedestalCam::fgFirst      = -57.5;
+const Axis_t  MHPedestalCam::fgLast       = 192.5;
 const TString MHPedestalCam::gsHistName   = "Pedestal";
 const TString MHPedestalCam::gsHistTitle  = "Pedestal";
@@ -166 +166 @@
   SetHistYTitle(gsHistYTitle.Data());
 
+  SetFitStart();
 }
 
@@ -280 +281 @@
   
 
-  fPedestalsOut = (MPedestalCam*)pList->FindObject(AddSerialNumber(fNamePedestalCamOut),"MPedestalCam");
+  if (!fPedestalsOut)
+    fPedestalsOut = (MPedestalCam*)pList->FindObject(AddSerialNumber(fNamePedestalCamOut),"MPedestalCam");
 
   if (!fPedestalsOut)
@@ -558 +560 @@
       //      else
       TH1F *gaush = hist.GetHGausHist();
-      hist.FitGaus("RQ0",-0.5,gaush->GetBinCenter(gaush->GetXaxis()->GetLast()));
-      //
-      // 3) Take histogram means and RMS
-      //
-      hist.BypassFit();
+      hist.FitGaus("RQ0",fFitStart,gaush->GetBinCenter(gaush->GetXaxis()->GetLast()));
       //
       // 4) Check for oscillations
@@ -596 +594 @@
       // 2) Fit the Hi Gain histograms with a Gaussian
       //
-      hist.FitGaus();
-      //
-      // 3) Take histogram means and RMS
-      //
-      hist.BypassFit();
+      TH1F *gaush = hist.GetHGausHist();
+      hist.FitGaus("RQ0",fFitStart,gaush->GetBinCenter(gaush->GetXaxis()->GetLast()));
       //
       // 4) Check for oscillations
@@ -640 +635 @@
       // 2) Fit the Hi Gain histograms with a Gaussian
       //
-      hist.FitGaus();
-      //
-      // 3) Take histogram means and RMS
-      //
-      hist.BypassFit();
+      TH1F *gaush = hist.GetHGausHist();
+      hist.FitGaus("RQ0",fFitStart,gaush->GetBinCenter(gaush->GetXaxis()->GetLast()));
       //
       // 4) Check for oscillations
@@ -698 +690 @@
       
       MCalibrationPix &pix = (*fCam)[i];
+      MPedestalPix    &ped = (*fPedestalsOut)[i];
       pix.SetHiGainMean     ( pix.GetHiGainMean()     / fExtractHiGainSlices  );
       pix.SetLoGainMean     ( pix.GetLoGainMean()     / fExtractHiGainSlices  );
+
+      ped.SetPedestal(pix.GetHiGainMean());
       //
       // Mean error goes with PedestalRMS/Sqrt(entries) -> scale with sqrt(slices)
@@ -705 +700 @@
       pix.SetHiGainMeanVar  ( pix.GetHiGainMeanVar()  / fExtractHiGainSlices );
       pix.SetLoGainMeanVar  ( pix.GetHiGainMeanVar()  / fExtractHiGainSlices );
+
       //
       // Sigma goes like PedestalRMS -> scale with sqrt(slices)    
@@ -710 +706 @@
       pix.SetHiGainSigma    ( pix.GetHiGainSigma()    / sqslices  );
       pix.SetLoGainSigma    ( pix.GetLoGainSigma()    / sqslices  );
+
+      ped.SetPedestalRms(pix.GetHiGainSigma());
       //
       // Sigma error goes like PedestalRMS/2.(entries) -> scale with sqrt(slices)