Changeset 3038

Timestamp:

02/05/04 17:50:16 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• manalysis/MPedestalCam.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -55 +55 @@
        renormalized to values per slice (this is not to screw up with the 
        pedestal RMS and the derived uncertainties)
+
+   * manalysis/MPedCalcPedRun.cc
+     - added author
+     - split calculation of RMS into two steps to make it more clear to 
+       the user
 
 

trunk/MagicSoft/Mars/manalysis/MPedestalCam.cc

@@ -143 +143 @@
 // Get i-th pixel (pixel number)
 //
-const MHPedestalPixel &MPedestalCam::operator()(UInt_t i) const
+MHPedestalPixel &MPedestalCam::operator()(UInt_t i) const
 {
   return *static_cast<MHPedestalPixel*>(fHArray->UncheckedAt(i));
@@ -155 +155 @@
 
     fArray->ForEach(TObject, Clear)();
-    //    fHArray->ForEach(TObject, Clear)();
 
     fTotalEntries     = 0;
     fExtractSlices    = 0.;
 }
+
+
+// --------------------------------------------------------------------------
+//
+// Our own clone function is necessary since root 3.01/06 or Mars 0.4
+// I don't know the reason
+//
+TObject *MPedestalCam::Clone(const char *) const
+{
+
+  const Int_t n1 = fArray->GetSize();
+  const Int_t n2 = fHArray->GetSize();
+  
+  //
+  // FIXME, this might be done faster and more elegant, by direct copy.
+  //
+  MPedestalCam *cam = new MPedestalCam;
+  
+  cam->fArray->Expand(n1);
+  cam->fHArray->Expand(n2);
+  
+  for (int i=0; i<n1; i++)
+    {
+      delete (*cam->fArray)[i];
+      (*cam->fArray)[i] = (*fArray)[i]->Clone();
+    }
+
+  for (int i=0; i<n2; i++)
+    {
+      delete (*cam->fHArray)[i];
+      (*cam->fHArray)[i] = (*fHArray)[i]->Clone();
+    }
+  return cam;
+}
+
+
 
 // --------------------------------------------------------------------------
@@ -170 +205 @@
 {
 
-  fArray->ForEach(MPedestalPix, InitUseHists)();
   fHArray->Delete();
 
@@ -189 +223 @@
     }
   
+
   Float_t slices = (Float_t)signal->GetNumUsedFADCSlices();
 
@@ -218 +253 @@
           MHPedestalPixel &hist = (*this)(i);
           hist.ChangeHistId(i);
+          MPedestalPix &pix = (*this)[i];
+          pix.InitUseHists();
         }
     }
@@ -229 +266 @@
   for (Int_t i=0; i<n; i++)
     {
+
       const MExtractedSignalPix &pix = (*signal)[i];
       
       const Float_t sig = pix.GetExtractedSignalHiGain();
-      const Float_t sigPerSlice = sig/fExtractSlices;
-      
+
       MHPedestalPixel &hist = (*this)(i);
-
-      hist.FillCharge(sigPerSlice);
-      hist.FillChargevsN(sigPerSlice);
+      //
+      // Don't fill signal per slice, we get completely screwed up 
+      // with the sigma. Better fill like it is and renorm later
+      //
+      //      const Float_t sigPerSlice = sig/fExtractSlices;
+      //      hist.FillCharge(sigPerSlice);
+      //      hist.FillChargevsN(sigPerSlice);
+      const Float_t signal = sig;
+      hist.FillCharge(signal);
+      hist.FillChargevsN(signal);
     }
   
@@ -247 +291 @@
     for (Int_t i=0; i<fHArray->GetSize(); i++)
     {
+
         MHPedestalPixel &hist = (*this)(i);
 
@@ -262 +307 @@
         //
         hist.FitCharge();
+        hist.Renorm(fExtractSlices);
+
     }
     return kTRUE;
@@ -343 +390 @@
 
   const Float_t pederr   = rms/TMath::Sqrt((Float_t)fTotalEntries);
-  const Float_t rmserr   = rms/TMath::Sqrt((Float_t)fTotalEntries/2.);
+  const Float_t rmserr   = rms/TMath::Sqrt((Float_t)fTotalEntries*2.);
 
   const MHPedestalPixel &hist = (*this)(idx);
 
   const Float_t mean     = hist.GetChargeMean();
-  const Float_t meanerr  = hist.GetChargeMeanErr() * TMath::Sqrt((Float_t)fExtractSlices);
-  const Float_t sigma    = hist.GetChargeSigma()   * TMath::Sqrt((Float_t)fExtractSlices);
-  const Float_t sigmaerr = hist.GetChargeSigmaErr()* TMath::Sqrt((Float_t)fExtractSlices);
+  const Float_t meanerr  = hist.GetChargeMeanErr();
+  const Float_t sigma    = hist.GetChargeSigma() ;
+  const Float_t sigmaerr = hist.GetChargeSigmaErr();
   const Float_t prob     = hist.GetChargeProb();
 
@@ -368 +415 @@
       break;
     case 4:
-      //      if ((*this)[idx].IsFitValid())
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = mean;
-        //      else
-        //        return kFALSE;
       break;
     case 5:
-      //      if ((*this)[idx].IsFitValid())
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = meanerr;
-        //      else
-        //        return kFALSE;
       break;
     case 6:
-      //      if ((*this)[idx].IsFitValid())
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = sigma;
-        //      else
-        //        return kFALSE;
       break;
     case 7:
-      //      if ((*this)[idx].IsFitValid())
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = sigmaerr;
-        //      else
-        //        return kFALSE;
       break;
     case 8:
-      //      if ((*this)[idx].IsFitValid())
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = prob;
-        //      else
-        //        return kFALSE;
       break;
     case 9:
-      //      if ((*this)[idx].IsFitValid())
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = 2.*(ped-mean)/(ped+mean);
-        //      else
-        //        return kFALSE;
       break;
     case 10:
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = TMath::Sqrt((pederr*pederr + meanerr*meanerr) * (ped*ped + mean*mean))
             *2./(ped+mean)/(ped+mean);
       break;
     case 11:
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = 2.*(pederr - meanerr)/(pederr + meanerr);
       break;
     case 12:
-      val = 2.*(rms-sigma)/(rms+sigma);
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
+      val = 2.*(sigma-rms)/(sigma+rms);
       break;
     case 13:
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
       val = TMath::Sqrt((rmserr*rmserr + sigmaerr*sigmaerr) * (rms*rms + sigma*sigma))
             *2./(rms+sigma)/(rms+sigma);
       break;
     case 14:
-      //      if ((*this)[idx].IsFitValid())
-      val = 2.*(rmserr - sigmaerr)/(rmserr + sigmaerr);
-              //      else
-        //        return kFALSE;
+      if (!(*this)(idx).IsFitOK())
+        return kFALSE;
+      val = 2.*(sigmaerr - rmserr)/(sigmaerr + rmserr);
       break;
     default:
@@ -431 +479 @@
 void MPedestalCam::DrawPixelContent(Int_t idx) const
 {
-  (*this)[idx].Draw();
-}
+  (*this)(idx).Draw();
+}