Changeset 7194 for trunk/MagicSoft/Mars/msignal

Timestamp:

07/15/05 13:25:57 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msignal

Files:

• MC_calibration_weights_UV46.dat (added)
• MC_cosmics_weights46.dat (added)
• MExtractTimeAndChargeDigitalFilter.cc (modified) (16 diffs)
• MExtractTimeAndChargeDigitalFilter.h (modified) (4 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc

@@ -70 +70 @@
 #include "MRawRunHeader.h"
 #include "MCalibrationPattern.h"
+#include "MExtractedSignalCam.h"
 
 #include "MPedestalPix.h"
@@ -142 +143 @@
   if (fWindowSizeHiGain > availhirange)
   {
-      // Please simplify this!
-      *fLog << warn << GetDescriptor() 
-            << Form("%s%2i%s%2i%s%2i%s",": Hi Gain window size: ",fWindowSizeHiGain,
-                    " is bigger than available range: [",(int)fHiGainFirst,",",(int)fHiGainLast,"]") << endl;
+      *fLog << warn << GetDescriptor() << ": Hi Gain window size: " << Form("%2i",fWindowSizeHiGain);
+      *fLog << " is bigger than available range: [" << Form("%2i", (int)fHiGainFirst);
+      *fLog << "," << Form("%21", (int)fHiGainLast) << "]" << endl;
+
       fHiGainLast = fHiGainFirst + fWindowSizeHiGain;
-      *fLog << warn << GetDescriptor() 
-            << ": Will set the upper range to: " << (int)fHiGainLast << endl;
+
+      *fLog << warn << GetDescriptor() << ": Will set the upper range to: " << (int)fHiGainLast << endl;
     }
   
@@ -163 +164 @@
       if (fWindowSizeLoGain > availlorange)
         {
-      // Please simplify this!
-          *fLog << warn << GetDescriptor() 
-                << Form("%s%2i%s%2i%s%2i%s",": Lo Gain window size: ",fWindowSizeLoGain,
-                        " is bigger than available range: [",(int)fLoGainFirst,",",(int)fLoGainLast,"]") << endl;
-          fLoGainLast = fLoGainFirst + fWindowSizeLoGain;
-          *fLog << warn << GetDescriptor() 
-            << ": Will set the upper range to: " << (int)fLoGainLast << endl;
+            *fLog << warn << GetDescriptor() << ": Lo Gain window size: " << Form("%2i",fWindowSizeLoGain);
+            *fLog << " is bigger than available range: [" << Form("%2i", (int)fLoGainFirst);
+            *fLog << "," << Form("%21", (int)fLoGainLast) << "]" << endl;
+
+            fLoGainLast = fLoGainFirst + fWindowSizeLoGain;
+
+            *fLog << warn << GetDescriptor() << ": Will set the upper range to: " << (int)fLoGainLast << endl;
         }
       
@@ -205 +206 @@
 // --------------------------------------------------------------------------
 //
-// InitArrays
-//
-// Gets called in the ReInit() and initialized the arrays
-//
-Bool_t MExtractTimeAndChargeDigitalFilter::InitArrays()
-{
-    if (!fRunHeader)
+// The weights are determined using GetAutimaticWeights().
+//
+// kFALSE is returned if it returned an error.
+// kTRUE  is returned if no new weights were set.
+//
+// If new weights are set
+//  fTimeShiftHiGain
+//  fTimeShiftLoGain
+//  fNumHiGainSamples
+//  fNumLoGainSamples
+//  fSqrtHiGainSamples
+//  fSqrtLoGainSamples
+// and
+//  fSignals->SetUsedFADCSlices(...)
+// is updated accordingly.
+//
+Bool_t MExtractTimeAndChargeDigitalFilter::GetWeights()
+{
+    switch (GetAutomaticWeights())
+    {
+    case kERROR: // An error occured
         return kFALSE;
-
-    const Int_t rangehi = (Int_t)(fHiGainLast - fHiGainFirst + 1 + fHiLoLast);
-    const Int_t rangelo = (Int_t)(fLoGainLast - fLoGainFirst + 1);
-
-    fHiGainSignal.Set(rangehi);
-    fLoGainSignal.Set(rangelo);
-
-    // Try to get automatic weights
-    if (!ReadAutomaticWeightsFile())
-        return kFALSE;
-
-    // If still no weights set try conservative way
-    if (fNameWeightsFileSet.IsNull())
-        if (!ReadWeightsFile(fNameWeightsFile))
-            return kFALSE;
-
-    fTimeShiftHiGain = (Float_t)fHiGainFirst + 0.5 + 1./fBinningResolutionHiGain;
+    case kFALSE: // No new weights set
+        return kTRUE;
+    }
+
+    // New weights set
+    fTimeShiftHiGain = 0.5 + 1./fBinningResolutionHiGain + fHiGainFirst;
     fTimeShiftLoGain = 0.5 + 1./fBinningResolutionLoGain;
+
     //
     // We need here the effective number of samples which is about 2.5 in the case of a window
@@ -240 +245 @@
     fSqrtLoGainSamples = TMath::Sqrt(fNumLoGainSamples);
 
+    // From MExtractTimeAndCharge::ReInit
+    if (fSignals)
+        fSignals->SetUsedFADCSlices(fHiGainFirst, fHiGainLast+fHiLoLast, fNumHiGainSamples,
+                                    fLoGainFirst, fLoGainLast, fNumLoGainSamples);
     return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+// InitArrays
+//
+// Gets called in the ReInit() and initialized the arrays
+//
+Bool_t MExtractTimeAndChargeDigitalFilter::InitArrays()
+{
+    if (!fRunHeader)
+        return kFALSE;
+
+    const Int_t rangehi = (Int_t)(fHiGainLast - fHiGainFirst + 1 + fHiLoLast);
+    const Int_t rangelo = (Int_t)(fLoGainLast - fLoGainFirst + 1);
+
+    fHiGainSignal.Set(rangehi);
+    fLoGainSignal.Set(rangelo);
+
+    return GetWeights();
+}
+
+// --------------------------------------------------------------------------
+//
+// Check if reading a new weights file is necessary because the calibration
+// pattern has changed. (Cannot be done in ReInit, because at this time
+// the calibration pattern is not available.
+// Then process the event.
+//
+Int_t MExtractTimeAndChargeDigitalFilter::Process()
+{
+    // Change Weights if the calibration patter changes
+    if (!GetWeights())
+        return kERROR;
+
+    // Process event
+    return MExtractTimeAndCharge::Process();
 }
 
@@ -643 +689 @@
 //----------------------------------------------------------------------------
 //
+// If automatic weights are requested, no default weights (name.IsNull())
+// are requested, fRunHeader is available and fRunHeader->IsMonteCarloRun()
+// is true prepend "MC_" in from of the name.
+//
+// return poth+name;
+//
+TString MExtractTimeAndChargeDigitalFilter::CompileWeightFileName(TString path, const TString &name) const
+{
+    if (fAutomaticWeights && !name.IsNull() && fRunHeader && fRunHeader->IsMonteCarloRun())
+        path += "MC_";
+
+    path += name;
+
+    return path;
+}
+
+//----------------------------------------------------------------------------
+//
 // Read a pre-defined weights file into the class. 
 // This is mandatory for the extraction
@@ -648 +712 @@
 // If filenname is empty, then all weights will be set to 1.
 //
-Bool_t MExtractTimeAndChargeDigitalFilter::ReadWeightsFile(TString filename)
-{
-  
-    // This is a fix for TEnv files edited with windows editors
-    filename.ReplaceAll("\015", "");
-
-    fAmpWeightsHiGain .Set(fBinningResolutionHiGain*fWindowSizeHiGain);
-    fAmpWeightsLoGain .Set(fBinningResolutionLoGain*fWindowSizeLoGain);
-    fTimeWeightsHiGain.Set(fBinningResolutionHiGain*fWindowSizeHiGain);
-    fTimeWeightsLoGain.Set(fBinningResolutionLoGain*fWindowSizeLoGain);
-
-    if (fNameWeightsFile.IsNull() || fNameWeightsFile=="-")
-    {
+// Returns:
+//  kTRUE:  new weights set
+//  kFALSE: no  weights set
+//  kERROR: error
+//
+Int_t MExtractTimeAndChargeDigitalFilter::ReadWeightsFile(TString filename, TString path)
+{
+    if (filename.IsNull())
+    {
+        fAmpWeightsHiGain .Set(fBinningResolutionHiGain*fWindowSizeHiGain);
+        fAmpWeightsLoGain .Set(fBinningResolutionLoGain*fWindowSizeLoGain);
+        fTimeWeightsHiGain.Set(fBinningResolutionHiGain*fWindowSizeHiGain);
+        fTimeWeightsLoGain.Set(fBinningResolutionLoGain*fWindowSizeLoGain);
+
         fAmpWeightsHiGain.Reset(1);
         fTimeWeightsHiGain.Reset(1);
         fAmpWeightsLoGain.Reset(1);
         fTimeWeightsLoGain.Reset(1);
-        fNameWeightsFileSet = "-";
         return kTRUE;
     }
 
+    // Add "MC_" in front of the filename if necessary
+    filename = CompileWeightFileName(path, filename);
+
     //filename = MJob::ExpandPath(filename);
 
     if (fNameWeightsFileSet==filename)
-        return kTRUE;
+        return kFALSE; // No file read
 
     ifstream fin(filename.Data());
@@ -679 +746 @@
         *fLog << err << GetDescriptor() << ": ERROR - Cannot open file " << filename << ": ";
         *fLog << strerror(errno) << endl;
-        return kFALSE;
-    }
-
-    *fLog << all << "DigitalFilter: Reading weights in " << filename << "..." << flush;
+        return kERROR;
+    }
+
+    *fLog << all << GetDescriptor() << ": Reading weights in " << filename << "..." << flush;
 
     Int_t len = 0;
@@ -705 +772 @@
             {
                 *fLog << err << "ERROR - 'High Gain Weights' found twice in line #" << line << "." << endl;
-                return kFALSE;
+                return kERROR;
             }
 
@@ -712 +779 @@
                 *fLog << err << "ERROR - Wrong number of arguments in line #" << line << ":" << endl;
                 *fLog << str << endl;
-                return kFALSE;
+                return kERROR;
             }
 
@@ -727 +794 @@
             {
                 *fLog << err << "ERROR - 'Lo Gain Weights' found twice in line #" << line << "." << endl;
-                return kFALSE;
+                return kERROR;
             }
 
@@ -734 +801 @@
                 *fLog << err << "ERROR - Wrong number of arguments in line #" << line << ":" << endl;
                 *fLog << str << endl;
-                return kFALSE;
+                return kERROR;
             }
 
@@ -758 +825 @@
             *fLog << err << "ERROR - Wrong number of arguments in line #" << line << ":" << endl;
             *fLog << str << endl;
-            return kFALSE;
+            return kERROR;
         }
 
@@ -770 +837 @@
     if (cnt != len)
     {
-        *fLog << err << "Size mismatch in weights file " << filename << endl;
-        return kFALSE;
+        *fLog << err << "ERROR - Size mismatch in weights file " << filename << endl;
+        return kERROR;
     }
 
     if (!hi)
     {
-        *fLog << err << "No correct header found in weights file " << filename << endl;
-        return kFALSE;
+        *fLog << err << "ERROR - No correct header found in weights file " << filename << endl;
+        return kERROR;
     }
 
@@ -824 +891 @@
 //----------------------------------------------------------------------------
 //
-// If automatic weight disabled or run type is pedestal no action is taken
-// and kTRUE is returned.
-// 
-// Otherwise:
-//  - If runtype is kRTData and filename starts with "MC_weights" "MC" is
-//    replaced by "cosmics"
-//  - If runtype is kRTMonteCarlo and filename starts with "cosmics_weights"
-//    "cosmics" is replaced by "MC"
-//  - If name filename begins with "calibration_weights_" for blue and green
-//    "UV" is replaced by "blue"
-//  - If name filename begins with "calibration_weights_" for UV and CT1
-//    "blue" is replaced by "UV"
-//
-// The replacements are done in this order. If no match is found no
-// replacement is done. Afterwards ReadWeightsFile for the resulting
-// path-/filename is called.
-//
-Bool_t MExtractTimeAndChargeDigitalFilter::ReadAutomaticWeightsFile()
-{
-    if (!fAutomaticWeights || fRunHeader->GetRunType()==MRawRunHeader::kRTPedestal)
-        return kTRUE;
-
-    if (fNameWeightsFile.IsNull())
-    {
-        *fLog << err << "ERROR - Cannot get automatic weights without default filename." << endl;
-        return kFALSE;
-    }
-
-    *fLog << dbg << endl;
-    *fLog << underline << GetDescriptor() << endl;
-    *fLog << " Trying to get automatic weight for " << fNameWeightsFile << endl;
-
+// The default (+ prepending possible "MC_") is read for:
+//
+//   - RunType: Pedestal (independant of fAutomaticWeights)
+//   - fAutomaticWeights disabled
+//
+//  if fAutomaticWeights enabled:
+//   - fNameWeightsFile.IsNull()
+//   - !fCalibPattern
+//   - fCalibPattern->GetPulserColor()==MCalibrationCam::kNONE
+//
+// If automatic weights are enabled, the case above didn't take place and
+// fNameWeightsFile starts with "calibration_weights_"
+//   - the color (blue, UV) is replaced by the appropriate one
+//     taken from the calibration pattern
+//
+// In most cases a debug output is printed. Further output about the color
+// determination can be switched on with debug level > 5;
+//
+// Returns:
+//  kFALSE: No new weights set
+//  kTRUE:  New weights set
+//  kERROR: Error
+//
+Int_t MExtractTimeAndChargeDigitalFilter::GetAutomaticWeights()
+{
     const Ssiz_t pos = fNameWeightsFile.Last('/')+1;
     const Ssiz_t len = fNameWeightsFile.Length();
 
+    // Split file name in path and name
     TString path = fNameWeightsFile(0, pos>=0?pos:len);
     TString name = fNameWeightsFile(pos>=0?pos:0, len);
 
-    if (name.BeginsWith("cosmics_weights") && fRunHeader->GetRunType()==MRawRunHeader::kRTMonteCarlo)
-    {
-        name.Remove(0, 7);
-        name.Prepend("MC");
-        *fLog << "Run type: Monte Carlo data" << endl;
-    }
-    if (name.BeginsWith("MC_weights") && fRunHeader->GetRunType()==MRawRunHeader::kRTData)
-    {
-        name.Remove(0, 2);
-        name.Prepend("cosmics");
-        *fLog << "Run type: Cosmics data" << endl;
+    // Remove trailing "MC_" for automatic weights
+    if (fAutomaticWeights && name.BeginsWith("MC_"))
+        name.Remove(0, 3);
+
+    // In case of a pedetsal run no calibration pattern can be available
+    // the default weights are always used.
+    if (fRunHeader->GetRunType()==MRawRunHeader::kRTPedestal)
+    {
+        *fLog << dbg << "Pedestal file... using default weights: " << fNameWeightsFile << endl;
+        return ReadWeightsFile(name, path);
+    }
+
+    // If automatic weights are switched off use default weights
+    if (!fAutomaticWeights)
+    {
+        *fLog << dbg << "Automatic weights switched off... using default weights: " << fNameWeightsFile << endl;
+        return ReadWeightsFile(name, path);
+    }
+
+    // If automatic weights are switched on but no filename is given raise error
+    if (fNameWeightsFile.IsNull())
+    {
+        *fLog << err << "ERROR - Cannot get automatic weights without default filename." << endl;
+        return kERROR;
+    }
+
+    // If this is no pedestal run, automatic weights are requested and a
+    // filename for the weights file is given pedestal-extraction from
+    // cosmics data is assumed.
+    if (!fCalibPattern)
+    {
+        *fLog << dbg << "No decoded calibration pattern available... using default weights: " << fNameWeightsFile << endl;
+        return ReadWeightsFile(name, path);
+    }
+
+    const Bool_t debug = gLog.GetDebugLevel()>5;
+
+    // If no calibration pattern is available do not change the
+    // current weighs or current weights file name.
+    if (fCalibPattern->GetPulserColor()==MCalibrationCam::kNONE)
+    {
+        // If we are extracting data and the calibration pattern is kNONE
+        // we assume that it is a data file without interleaved events
+        // and calibration pattern information available.
+        if ((fRunHeader->GetRunType()!=MRawRunHeader::kRTData && !fRunHeader->IsMonteCarloRun()) || debug)
+            *fLog << dbg << "No calibration color set so far... guessing default: " << fNameWeightsFile << endl;
+
+        return ReadWeightsFile(name, path);
+    }
+
+    if (debug)
+    {
+        *fLog << dbg << endl;
+        *fLog << underline << GetDescriptor() << endl;
+        *fLog << " Trying to get automatic weight for " << fNameWeightsFile << endl;
+        *fLog << " Run type: ";
     }
 
     if (name.BeginsWith("calibration_weights_") && fCalibPattern)
     {
-        *fLog << " Run type: Calibration with color " << (int)fCalibPattern->GetPulserColor() << ", setting ";
+        if (debug)
+            *fLog << " Calibration with color " << fCalibPattern->GetPulserColorStr() << ", setting ";
         switch (fCalibPattern->GetPulserColor())
         {
         case MCalibrationCam::kBLUE:  // 2
         case MCalibrationCam::kGREEN: // 1
-            *fLog << "blue/green" << endl;
+            if (debug)
+                *fLog << "blue/green, ";
             name.ReplaceAll("UV", "blue");
             break;
@@ -888 +996 @@
         case MCalibrationCam::kUV:    // 3
         case MCalibrationCam::kCT1:   // 0
-            *fLog << "UV/CT1" << endl;
+            if (debug)
+                *fLog << "UV/CT1, ";
             name.ReplaceAll("blue", "UV");
             break;
         case MCalibrationCam::kNONE:
-            *fLog << "None" << endl;
             break;
-        default: // kNONE=4
-            *fLog << "???" << endl;
-        }
-    }
-
-    path += name;
-
-    return ReadWeightsFile(path);
+        default: // kNone + etc
+            *fLog << err << "ERROR - Cannot get automatic weights for " << fCalibPattern->GetPulserColorStr() << endl;
+            return kERROR;
+        }
+    }
+
+    return ReadWeightsFile(name, path);
 }
 

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.h

@@ -69 +69 @@
   // MExtractTimeAndChargeDigitalFilter
   void    CalcBinningResArrays();
-  Bool_t  ReadAutomaticWeightsFile();
+  Int_t   GetAutomaticWeights();
+  Bool_t  GetWeights();
+  Int_t   ReadWeightsFile(TString filename, TString path="");
+  TString CompileWeightFileName(TString path, const TString &name) const;
 
   // MExtractTimeAndCharge
@@ -76 +79 @@
   // MTask
   Int_t PreProcess(MParList *pList);
+  Int_t Process();
 
 protected:
@@ -85 +89 @@
   ~MExtractTimeAndChargeDigitalFilter() { }
 
-  Bool_t ReadWeightsFile(TString filename);
-  
   Bool_t WriteWeightsFile(TString filename,
                           TH1F *shapehi, TH2F *autocorrhi,
@@ -92 +94 @@
 
 
-  void SetNameWeightsFile(TString s="") { fNameWeightsFile = s; fNameWeightsFileSet=""; }
+  void SetNameWeightsFile(TString s="") { s.ReplaceAll("\015", ""); // This is a fix for TEnv files edited with windows editors
+      fNameWeightsFile = s; fNameWeightsFileSet=""; }
   void EnableAutomaticWeights(Bool_t b=kTRUE) { fAutomaticWeights = b; }