Project

General

Profile

Tests Programming » History » Version 2

Ole Hansen, 03/05/2018 03:17 PM

1 1 Ole Hansen
h1. Dynamically defined logical expressions (tests/cuts)
2
3 2 Ole Hansen
15 April 2004
4
5 1 Ole Hansen
h2. Overview
6
7
A useful feature of Podd is the ability to define arithmetic and logical expressions dynamically (on the fly) without the need to recompile any code. Underlying both types of expressions are "global variables" that are represented by the "@THaVar@":http://hallaweb.jlab.org/podd/html/THaVar.html class and collected in the global variable list "@THaVarList@":http://hallaweb.jlab.org/podd/html/THaVarList.html. Each analysis module (detector, apparatus, physics module) normally adds internal variables of interest to the global variable list @gHaVars@ as part of its [[Writing Detector Code#Initialization|initialization]] and removes those variables from the list in its destructor. For an example, see "@THaScintillator@":http://hallaweb.jlab.org/podd/html/THaScintillator.html. In this way, analysis results are conveniently available for use in expressions and by "@THaOutput@":http://hallaweb.jlab.org/podd/html/THaVarList.html.
8
9
Arithmetic expression parsing is implemented via the "@THaFormula@":http://hallaweb.jlab.org/podd/html/THaFormula.html class, while tests/cuts are special cases of @THaFormula@s which evaluate to either 1 or 0 (true or false) and are represented by the "@THaCut@":http://hallaweb.jlab.org/podd/html/THaCut.html class. In the C++ analyzer, @THaFormulas@ are primarily used in the "@THaOutput@":http://hallaweb.jlab.org/podd/html/THaOutout.html system, while @THaCuts@ are used for two purposes: (1) in @THaOutput@ as conditions on histograms; and (2) to control the flow of the data analysis in the standard analysis algorithm. The latter is supported by a special global list of tests, "@THaCutList@":http://hallaweb.jlab.org/podd/html/THaCutList.html, @THaCutList@ supports the concept of "blocks" of cuts which can be evaluated as a unit. Blocks are evaluated at the end of each stage of the analysis (Decode/Reconstruct/Physics).
10
11
The interactive interface to the analyzer, "@THaInterface@":http://hallaweb.jlab.org/podd/html/THaInterface.html, automatically creates an instance of "@THaVarList@":http://hallaweb.jlab.org/podd/html/THaVarList.html and "@THaCutList@":http://hallaweb.jlab.org/podd/html/THaCutList.html upon startup. These are called the "global variable list" and "global cut list", respectively. They are accessible via the global program variables @gHaVars@ and @gHaCuts@ from anywhere in the analyzer if you  include @THaGlobals.h@.
12
13
h2. Examples
14
15
A short tutorial on using the test/cut classes along with the global variable system follows. As mentioned, the relevant classes are:
16
17
* "@THaCutList@":http://hallaweb.jlab.org/podd/html/THaCutList.html -- Define and evaluate cuts. Automatically manages cuts with the help of internal lists.
18
* "@THaVarList@":http://hallaweb.jlab.org/podd/html/THaVarList.html -- Manage "global" symbolic variables.
19
* "@THaCut@":http://hallaweb.jlab.org/podd/html/THaCut.html -- Definition of a single cut. This is a low-level class which should normally only be used internally by @THaCutList@.
20
* "@THaVar@":http://hallaweb.jlab.org/podd/html/THaVar.html -- Definition of a single symbolic variable. This is a low-level class which should normally only be used internally by @THaVarList@.
21
22
Here is a sample session to demonstrate the use of these classes:
23
24
<pre>
25
analyzer [0] double xvar=10             // Define a variable xvar
26
analyzer [1] gHaVars->Define("x",xvar)  // Add xvar to global var list, name it "x"
27
analyzer [2] gHaVars->PrintFull()       // Show all global vars defined including their current values
28
  OBJ: THaVar     x     x
29
  (Double_t)[1]  10
30
31
analyzer [3] gHaCuts->Define("cut1","x>0")      // Define a cut named "cut1" that is true if x>0.
32
analyzer [4] gHaCuts->Define("cut2","abs(x)>5") // dto., but |x|>5
33
analyzer [5] gHaCuts->Print()                   // List all defined cuts
34
  Name  Def       T  Block    Called     Passed
35
  -------------------------------------------------
36
  cut1  x>0       0  Default  0          0 (0.0%)
37
  cut2  abs(x)>5  0  Default  0          0 (0.0%)
38
39
analyzer [6] gHaCuts->Eval()                   // Evaluate all defined cuts
40
analyzer [7] gHaCuts->Print()
41
  Name  Def       T  Block    Called     Passed
42
  -------------------------------------------------
43
  cut1  x>0       1  Default  1          1 (100%)
44
  cut2  abs(x)>5  1  Default  1          1 (100%)
45
46
(Note the current value ("T") of each of the cuts. Both are true since both
47
conditions are true for x=10).
48
49
analyzer [8] xvar=2                           // Give xvar a new value
50
analyzer [9] gHaCuts->Eval()                  // Evaluate the cuts again
51
analyzer [10] gHaCuts->Print()
52
  Name  Def       T  Block    Called     Passed
53
  -------------------------------------------------
54
  cut1  x>0       1  Default  2          2 (100%)
55
  cut2  abs(x)>5  0  Default  2          1 (50%)
56
57
(Note that cut2 is now false since |x| is less than 5.)
58
59
analyzer [11] gHaCuts->Result("cut2")        // Retrieve result of cut2
60
                                             // NB: This does not re-evaluate the cut
61
  (Int_t)0
62
analyzer [11] gHaCuts->Result("cut1")
63
  (Int_t)1
64
</pre>
65
Tests may refer to other tests already defined in gHaCuts. When evaluating cuts containing other tests, the referenced tests are not re-evaluated, but the result of their last evaluation is used (as in the call to "@THaCutList::Result()@":http://hallaweb.jlab.org/podd/html/THaCutList.html#THaCutList:Result above). This ensures that the test statistics (i.e. number of calls/number of passes) remain consistent when evaluating all the tests for an event. Here's an example, continuing from above:
66
<pre>
67
analyzer [12] gHaCuts->Define("cut3","cut1&&!cut2") // Define cut based on two previously-defined cuts
68
analyzer [13] gHaCuts->Eval()                       // Evaluate defined cuts again
69
analyzer [14] gHaCuts->Print()
70
  Name  Def          T  Block    Called     Passed
71
  ----------------------------------------------------
72
  cut1  x>0          1  Default  3          3 (100%)
73
  cut2  abs(x)>5     0  Default  3          1 (33.3%)
74
  cut3  cut1&&!cut2  1  Default  1          1 (100%)
75
</pre>