Changeset 8406 for trunk/MagicSoft/Mars/mmovie

Timestamp:

04/15/07 10:38:08 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mmovie

Files:

• MMovieData.h (modified) (2 diffs)
• MMoviePrepare.cc (modified) (7 diffs)
• MMoviePrepare.h (modified) (5 diffs)
• MMovieWrite.cc (modified) (13 diffs)
• MMovieWrite.h (modified) (1 diff)

trunk/MagicSoft/Mars/mmovie/MMovieData.h

@@ -23 +23 @@
 
     Float_t fMax;
+    Float_t fMedianPedestalRms;
     UInt_t  fNumSlices;
 
@@ -41 +42 @@
     Float_t GetMax() const { return fMax; }
 
+    void SetMedianPedestalRms(Float_t m) { fMedianPedestalRms = m; }
+    Float_t GetMedianPedestalRms() const { return fMedianPedestalRms; }
+
     void SetNumSlices(UInt_t n) { fNumSlices=n; }
     UInt_t GetNumSlices() const { return fNumSlices; }

trunk/MagicSoft/Mars/mmovie/MMoviePrepare.cc

@@ -48 +48 @@
 //   MGeomCam
 //   MPedestalSubtractedEvt
+//   MPedestalFundamental [MPedestalCam]
 //   MCalibrateData
 //   MCalibrationRelTimeCam
@@ -70 +71 @@
 #include "MGeomCam.h"
 #include "MGeomPix.h"
+
+#include "MPedestalCam.h"
+#include "MPedestalPix.h"
+
 
 #include "MCalibrationRelTimeCam.h"
@@ -78 +83 @@
 #include "MRawRunHeader.h"
 #include "MCalibrateData.h"
-//#include "MExtractedSignalCam.h"
 #include "MPedestalSubtractedEvt.h"
 
@@ -111 +115 @@
 //   MGeomCam
 //   MPedestalSubtractedEvt
+//   MPedestalFundamental [MPedestalCam]
 //   MCalibrateData
 //   MCalibrationRelTimeCam
@@ -132 +137 @@
         return kFALSE;
     }
+    fPed = (MPedestalCam*)plist->FindObject("MPedestalFundamental", "MPedestalCam");
+    if (!fPed)
+    {
+        *fLog << err << "MPedestalFundamental [MPedestalCam] not found ... aborting." << endl;
+        return kFALSE;
+    }
     fCal = (MCalibrateData*)plist->FindTaskListWithTask("MCalibrateData")->FindTask("MCalibrateData");
     if (!fCal)
@@ -184 +195 @@
     const MArrayF &calco = fCal->GetCalibConsts();
 
-    Int_t n = fCam->GetNumPixWithArea(0);
-
-    MArrayF calco0(n);
+    Int_t n = fCam->GetNumPixWithAidx(0);
+
+    MArrayF arr(n);
 
     for (UInt_t i=0; i<fCam->GetNumPixels(); i++)
         if ((*fCam)[i].GetAidx()==0)
-            calco0[--n] = calco[i];
-
-    return TMath::Median(calco0.GetSize(), calco0.GetArray());
+            arr[--n] = calco[i];
+
+    return TMath::Median(arr.GetSize(), arr.GetArray());
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the median of the pedestal rms of all pixels with area index 0
+//
+Double_t MMoviePrepare::GetMedianPedestalRms() const
+{
+    Int_t n = fCam->GetNumPixWithAidx(0);
+
+    MArrayF rms(n);
+
+    for (UInt_t i=0; i<fCam->GetNumPixels(); i++)
+        if ((*fCam)[i].GetAidx()==0)
+            rms[--n] = (*fPed)[i].GetPedestalRms();
+
+    return TMath::Median(rms.GetSize(), rms.GetArray());
 }
 
@@ -252 +280 @@
     fOut->SetMax(max);
     fOut->SetNumSlices(width);
+    fOut->SetMedianPedestalRms(GetMedianPedestalRms());
 
     return kTRUE;

trunk/MagicSoft/Mars/mmovie/MMoviePrepare.h

@@ -9 +9 @@
 class MExtractor;
 class MMovieData;
-class MSignalCam;
-class MRawEvtHeader;
+class MPedestalCam;
 class MRawRunHeader;
 class MCalibrateData;
@@ -24 +23 @@
     MRawRunHeader          *fRun;    //! Run Header (sampling frequency)
     MGeomCam               *fCam;    //! Pixel size
-//    MCalibConstCam         *fCal;
+    MPedestalCam           *fPed;    //! Fundamental pedestal (slice by slice)
 
     MMovieData             *fOut;    //! Prepared output container
@@ -31 +30 @@
     UInt_t fLastSlice;
 
-    //Float_t fSlope;
-
     // MTask
     Int_t  PreProcess(MParList *plist);
-    //Bool_t ReInit(MParList *plist);
     Int_t  Process();
 
@@ -43 +39 @@
     // MMoviePrepare
     Double_t GetMedianCalibConst() const;
+    Double_t GetMedianPedestalRms() const;
 
 public:
@@ -53 +50 @@
 
 #endif
-

trunk/MagicSoft/Mars/mmovie/MMovieWrite.cc

@@ -60 +60 @@
 //
 //  The interpolation of the frames is done using a TSpline3. If the spline
-//  would extrapolate the contents is set to zero. Unsuitable pixels are
-//  interpolated frame by frame using the surrounding suitable pixels.
-//
-//
-// Input:
-//   - MGeomCam
-//   - MRawRunHeader
-//   - MRawEvtHeader
-//   - MSignalCam
-//   - MBadPixelsCam
-//   - MPedestalFundamental [MPedestalCam]
-//   - MMovieData
+//  would extrapolate due to the shift by the relative time calibration the
+//  contents is set to zero. Unsuitable pixels are interpolated frame by
+//  frame using the surrounding suitable pixels.
+//
+//  A few words about file size: MPEG is a motion compensation compression,
+//  which means that if a region of a past frame is shown again at the same
+//  place or somewhere else this region is referenced instead of encoded again.
+//  This means that in our case (almost all frames are identical!) the
+//  increase of file size is far from linear with the number of encoded events!
+//
+//
+//  Input:
+//   MGeomCam
+//   MRawRunHeader
+//   MRawEvtHeader
+//   MSignalCam
+//   MBadPixelsCam
+//   MMovieData
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -91 +97 @@
 #include "MGeomCam.h"
 #include "MGeomPix.h"
-
-#include "MPedestalCam.h"
-#include "MPedestalPix.h"
 
 #include "MH.h"
@@ -117 +120 @@
 //
 MMovieWrite::MMovieWrite(const char *name, const char *title)
-    : fPipe(0), fTargetLength(5), fThreshold(2), fNumEvents(1000), fFilename("movie.mpg")
+    : fPipe(0), fTargetLength(5), fThreshold(2), fNumEvents(25000), fFilename("movie.mpg")
 {
     fName  = name  ? name  : "MMovieWrite";
@@ -172 +175 @@
 
     fPipe = gSystem->OpenPipe(name, "w");
-    if (fPipe)
-        return kTRUE;
-
-    *fLog << err;
-    *fLog << "Pipe: " << name << endl;
-    *fLog << "Couldn't open pipe... aborting." << endl;
-    CheckPipe();
-    return kFALSE;
+    if (!fPipe)
+    {
+        *fLog << err;
+        *fLog << "Pipe: " << name << endl;
+        *fLog << "Couldn't open pipe... aborting." << endl;
+        CheckPipe();
+        return kFALSE;
+    }
+
+    *fLog << inf << "Setup pipe to ppmtoy4m and mpeg2enc to encode " << fFilename << "." << endl;
+
+    return kTRUE;
 
     //  1: 37M name += "-f 9 -E 40 -H -4 1 -2 1 --dualprime-mpeg2";
@@ -229 +236 @@
 //   - MSignalCam
 //   - MBadPixelsCam
-//   - MPedestalFundamental [MPedestalCam]
 //   - MMovieData
 //
@@ -265 +271 @@
     {
         *fLog << err << "MBadPixelsCam not found ... aborting." << endl;
-        return kFALSE;
-    }
-    fPed = (MPedestalCam*)plist->FindObject("MPedestalFundamental", "MPedestalCam");
-    if (!fPed)
-    {
-        *fLog << err << "MPedestalFundamental [MPedestalCam] not found ... aborting." << endl;
         return kFALSE;
     }
@@ -317 +317 @@
 // and everything above gets colors.
 //
-void MMovieWrite::SetPalette(Double_t rms, const TH1 &h) const
+Int_t MMovieWrite::SetPalette(Double_t rms, const TH1 &h) const
 {
     const Double_t min = h.GetMinimum();
@@ -335 +335 @@
 
     TArrayI col(99);
-    gStyle->CreateGradientColorTable(6, s, r, g, b, col.GetSize());
+
+    const Int_t rc = gStyle->CreateGradientColorTable(6, s, r, g, b, col.GetSize());
 
     // --- Overwrite the 'underflow' bin with white ---
@@ -347 +348 @@
 
     gStyle->SetPalette(col.GetSize(), col.GetArray());
+
+    return rc;
+}
+
+// --------------------------------------------------------------------------
+//
+// The created colors are not overwritten and must be deleted manually
+// because having more than 32768 color in a palette will crash
+// gPad->PaintBox
+//
+void MMovieWrite::DeletePalette(Int_t colidx) const
+{
+    for (int i=0; i<99; i++)
+    {
+        TColor *col = gROOT->GetColor(colidx+i);
+        if (col)
+            delete col;
+    }
 }
 
@@ -563 +582 @@
 // --------------------------------------------------------------------------
 //
-// Return the median of the pedestal rms of all pixels with area index 0
-//
-Double_t MMovieWrite::GetMedianPedestalRms() const
-{
-    Int_t n = fCam->GetNumPixWithAidx(0);
-
-    MArrayF rms(n);
-
-    for (UInt_t i=0; i<fCam->GetNumPixels(); i++)
-        if ((*fCam)[i].GetAidx()==0)
-            rms[--n] = (*fPed)[i].GetPedestalRms();
-
-    return TMath::Median(rms.GetSize(), rms.GetArray());
-}
-
-// --------------------------------------------------------------------------
-//
 Bool_t MMovieWrite::Process(TH1 &h, TVirtualPad &c)
 {
@@ -589 +591 @@
     //const Float_t len    = (slices-2)/freq;               // [ns]  length of data stream in data-time
 
-    const Double_t rms = GetMedianPedestalRms();
+    const Double_t rms = fIn->GetMedianPedestalRms();
 
     h.SetMinimum(fIn->GetMax()-(fIn->GetMax()-rms)*99/98); // rms0
@@ -601 +603 @@
 
     // Set new adapted palette for further rendering
-    SetPalette(rms, h);
+    const Int_t colidx = SetPalette(rms, h);
 
     // Get the pad containing the camera with the movie
@@ -638 +640 @@
             return kFALSE;
     }
+
+    DeletePalette(colidx);
 
     cout << setw(3) << GetNumExecutions() << ": " << Form("%6.2f", (float)numframes/(slices-2)) << " f/sl " << slices << " " << numframes+1 << endl;

trunk/MagicSoft/Mars/mmovie/MMovieWrite.h

@@ -54 +54 @@
     void TreatBadPixels(TH1 &h) const;
 
-    void SetPalette(Double_t rms, const TH1 &h) const;
+    Int_t SetPalette(Double_t rms, const TH1 &h) const;
+    void DeletePalette(Int_t colidx) const;
 
     void UpdateImage(TASImage &img, TVirtualPad &pad);