Changeset 8888 for trunk/MagicSoft/Mars/mpedestal

Timestamp:

05/14/08 12:03:25 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mpedestal

Files:

• MPedestalSubtract.cc (modified) (5 diffs)
• MPedestalSubtract.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mpedestal/MPedestalSubtract.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MPedestalSubtract.cc,v 1.9 2007-12-19 18:53:03 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MPedestalSubtract.cc,v 1.10 2008-05-14 11:03:24 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -20 +20 @@
 !   Author(s): Thomas Bretz, 10/2006 <mailto:tbretz@astro.uni-wuerzburg.de>
 !
-!   Copyright: MAGIC Software Development, 2000-2006
+!   Copyright: MAGIC Software Development, 2000-2008
 !
 !
@@ -51 +51 @@
 #include "MArrayB.h"
 
+#include "MRawRunHeader.h"
 #include "MRawEvtData.h"
 #include "MRawEvtPixelIter.h"
@@ -141 +142 @@
 }
 
+Bool_t MPedestalSubtract::ReInit(MParList *pList)
+{
+    fRunHeader = (MRawRunHeader*)pList->FindObject(AddSerialNumber("MRawRunHeader"));
+    if (!fRunHeader)
+    {
+        *fLog << err << AddSerialNumber("MRawRunHeader") << " not found... aborting." << endl;
+        return kFALSE;
+    }
+    return kTRUE;
+}
+
 // --------------------------------------------------------------------------
 //
@@ -150 +162 @@
     const Int_t numl = fRawEvt->GetNumLoGainSamples();
 
+    // Check if event is empty (presumably MC event -- sanity check)
+    if (numh+numl==0)
+        return kCONTINUE;
+
+    // Check for consistency (our simulation can do weird things!)
+    if (numh!=fRunHeader->GetNumSamplesHiGain())
+    {
+        *fLog << warn << "WARNING - Number of hi-gain samples (" << numh << ") ";
+        *fLog << " doesn't match run-header (" << fRunHeader->GetNumSamplesHiGain() << ")." << endl;
+    }
+    if (numl!=fRunHeader->GetNumSamplesLoGain())
+    {
+        *fLog << warn << "WARNING - Number of lo-gain samples (" << numl << ") ";
+        *fLog << " doesn't match run-header (" << fRunHeader->GetNumSamplesLoGain() << ")." << endl;
+    }
+
+    // Get scale between FADC units and 256 ;-)
     const UInt_t scale = fRawEvt->GetScale();
 

trunk/MagicSoft/Mars/mpedestal/MPedestalSubtract.h

@@ -6 +6 @@
 #endif
 
+class MRawRunHeader;
 class MRawEvtData;
 class MPedestalCam;
@@ -16 +17 @@
     static const TString fgNamePedestalSubtractedEvt;  //! "MPedestalSubtractedEvt"
 
+    MRawRunHeader          *fRunHeader;      //! Run Header
     MRawEvtData            *fRawEvt;         //! Input Raw data
     MPedestalCam           *fPedestals;      //! Pedestals of all pixels in the camera
@@ -24 +26 @@
 
     Int_t  PreProcess(MParList *pList);
+    Bool_t ReInit(MParList *pList);
     Int_t  Process();