Changeset 5911 for trunk/MagicSoft/Mars/mfbase

Timestamp:

01/20/05 14:27:59 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mfbase

Files:

• MF.cc (modified) (5 diffs)
• MF.h (modified) (3 diffs)
• MFilterList.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mfbase/MF.cc

@@ -350 +350 @@
     if (!fF)
     {
+        if (IsAllowEmpty())
+            return kTRUE;
+
         *fLog << err << dbginf << "No filter rule available." << endl;
         return kFALSE;
@@ -370 +373 @@
 Int_t MF::Process()
 {
-    return fF->CallProcess();
+    return fF ? fF->CallProcess() : kTRUE;
 }
 
@@ -379 +382 @@
 Int_t MF::PostProcess()
 {
-    return fF->CallPostProcess();
-}
-
-// --------------------------------------------------------------------------
-//
-// Return the result of the filter rule.
+    return fF ? fF->CallPostProcess() : kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the result of the filter rule. If no filter available the
+// condition is always true
 //
 Bool_t MF::IsExpressionTrue() const
 {
-    return fF->IsConditionTrue();
+    return fF ? fF->IsConditionTrue() : kTRUE;
 }
 
@@ -398 +402 @@
 TString MF::GetDataMember() const
 {
-    return fF->GetDataMember();
+    return fF ? fF->GetDataMember() : "";
 }
 
@@ -422 +426 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Print the rule.
+// if now filter is set up '<n/a>' is printed.
+//
 void MF::Print(Option_t *opt) const
 {
     *fLog << all << underline << GetDescriptor() << endl;
-    fF->Print();
+    if (fF)
+        fF->Print();
+    else
+        *fLog << "<n/a>";
     *fLog << endl << endl;
 }
+
+// --------------------------------------------------------------------------
+//
+// Check for corresponding entries in resource file and setup filters.
+//
+// Assuming your MF-filter is called (Set/GetName): MyFilter
+//
+// First you can setup whether the filter is inverted or not:
+//     MyFilter.Inverted: yes, no  <default=no>
+//
+// Now setup the condition, eg:
+//     MyFilter.Condition: MHillas.fSize<1000
+// or
+//     MyFilter.Condition: MHillas.fSize>500 && MHillas.fSize<10000
+//
+// If you want to use more difficult filters you can split the
+// condition into subcondistions. Subcondistions are identified
+// by {}-brackets. Avoid trailing 0's! For example:
+//
+//     MyFilter.Condition: MHillas.fSize>500 && {0} && {1}
+//     MyFilter.0: MHillas.fSize>1000
+//     MyFilter.1: MHillas.fSize<10000
+//
+// The numbering must be continous and start with 0. You can use
+// a subcondition more than once. All {}-brackets are simply replaced
+// by the correspodning conditions. The rules how conditions can
+// be written can be found in the class description of MF and MDataChain.
+//
+Int_t MF::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Bool_t rc = kFALSE;
+    if (IsEnvDefined(env, prefix, "Inverted", print))
+    {
+        rc = kTRUE;
+        SetInverted(GetEnvValue(env, prefix, "Inverted", IsInverted()));
+    }
+
+    if (!IsEnvDefined(env, prefix, "Condition", print))
+        return rc;
+
+    TString rule = GetEnvValue(env, prefix, "Condition", "");
+    rule.ReplaceAll(" ", "");
+
+    Int_t idx=0;
+    while (1)
+    {
+        TString cond;
+        if (IsEnvDefined(env, prefix, Form("%d", idx), print))
+        {
+            cond += "(";
+            cond += GetEnvValue(env, prefix, Form("%d", idx), "");
+            cond += ")";
+        }
+
+        if (cond.IsNull())
+            break;
+
+        rule.ReplaceAll(Form("{%d}", idx), cond);
+        idx++;
+    }
+
+    if (fF)
+    {
+        delete fF;
+        fF = 0;
+    }
+
+    if (rule.IsNull())
+    {
+        *fLog << warn << "MF::ReadEnv - WARNING: Empty condition found." << endl;
+        SetAllowEmpty();
+        return kTRUE;
+    }
+
+    SetAllowEmpty(kFALSE);
+
+    if (!(fF=ParseString(rule, 1)))
+    {
+        *fLog << err << "MF::ReadEnv - ERROR: Parsing '" << rule << "' failed." << endl;
+        return kERROR;
+    }
+
+    if (print)
+    {
+        *fLog << inf << "found: ";
+        fF->Print();
+        *fLog << endl;
+    }
+
+    return kTRUE;
+}

trunk/MagicSoft/Mars/mfbase/MF.h

@@ -32 +32 @@
     void StreamPrimitive(ofstream &out) const;
 
+    enum { kAllowEmpty = BIT(14) };
+
 public:
     MF();
@@ -37 +39 @@
     ~MF();
 
+    // MF
+    void SetAllowEmpty(Bool_t b=kTRUE) { b ? SetBit(kAllowEmpty) : ResetBit(kAllowEmpty); }
+    Bool_t IsAllowEmpty() const { return TestBit(kAllowEmpty); }
+
+    // MFilter
     TString GetRule() const { return fF ? fF->GetRule() : MFilter::GetRule(); }
     TString GetDataMember() const;
@@ -42 +49 @@
     Bool_t IsExpressionTrue() const;
 
+    // MTask
     Int_t PreProcess(MParList *pList);
     Int_t Process();
     Int_t PostProcess();
 
+    // TObject
     void Print(Option_t *opt="") const;
 
+    // MParContainer
     void SetVariables(const TArrayD &arr) { if (fF) fF->SetVariables(arr); }
+    Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
 
     ClassDef(MF, 0) // A Filter for cuts in any data member

trunk/MagicSoft/Mars/mfbase/MFilterList.cc

@@ -439 +439 @@
 }
 
-
+// --------------------------------------------------------------------------
+//
+// Loop over all set filters and distribute arr to their SetVariables.
+// Make sure, that the values are unique (not two filters using the
+// same index in the array with different meanings)
+//
 void MFilterList::SetVariables(const TArrayD &arr)
 {