Changeset 6792 for trunk/MagicSoft/Mars/mhcalib

Timestamp:

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

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhcalib

Files:

• MHPedestalCam.cc (modified) (9 diffs)
• MHPedestalCam.h (modified) (2 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)

trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.h

@@ -39 +39 @@
   Bool_t  fRenorm;                        // Flag if the results will be re-normalized
 
+  Axis_t  fFitStart;                      // Charge value to start the fitting procedure
+  
   MArrayD fSum;                           //! sum of values  
   MArrayD fSumSquare;                     //! sum of squared values
@@ -69 +71 @@
   void ResetHists();
 
-  void SetNamePedestalCamOut(const char *name) { fNamePedestalCamOut = name; }
-  void SetPedestalsOut ( MPedestalCam *cam) { fPedestalsOut = cam; }
-  void SetRenorm         (const Bool_t b=kTRUE )  {  fRenorm = b; }
+  void SetFitStart          ( const Axis_t a=-0.5  ) { fFitStart = a;  }
+  void SetNamePedestalCamOut( const char *name     ) { fNamePedestalCamOut = name; }
+  void SetPedestalsOut      ( MPedestalCam *cam    ) { fPedestalsOut = cam; }
+  void SetRenorm            ( const Bool_t b=kTRUE ) { fRenorm = b; }
   
   ClassDef(MHPedestalCam, 1)    // Histogram class for Charge Camera Pedestals