Event Analysis Package.

Event Analysis Package
Event definition
Event functions
Event functions
Event containers
Event iterators
Event algorithms
Factory and helper classes
Event factory and helper classes
The event analysis package provides means to analysis events generated by the LIGO detectors. Its major components are:
    Event - A set of classes defining a general purpose event
    Containers - A set of containers to store events
    Iterators - Iterators for traversing event containers
    Functions - Function objects acting on events
    Conditions - Conditions acting on events
    Algorithms - A set of algorithms to select and histogram events 

An event is an object consisting of a time, a type, an interferometer set and an arbitrary number of columns.

The layout of the class library is as follows

  Basic event handling:
    Type - Describes an event type
    Name - Describes an event name
    TypeInfo - Implementation details of the event type class
    ColumnType - Describes the supported types of event columns
    ColumnInfo - Describes an event column
    FixedColumns - Interface to fixed event columns
    Layout - Describes the column layout of an event
    LayoutInfo - Implementation details of the event layout class
    Event - The basic event object (inherits from Type and IfoSet)
    EventPtr - A smart pointer for storing events
    Factory - Defines the event factory which manages types and layouts
    IndexList - A template for a vector with a alphabetic index
    List - A list container for storing a series of events
    ListPtr - A smart pointer for the list container
    Chain - A chain of list containers
    ChainPtr - A smart pointer for the chain container
    Set - A set of chain containers
    TimeWindow - Handles a time window (duration and offset)
    Window - A container to manage the events within a time window
    Iterator - An iterator for traversing any event container
    ConstIterator - A const iterator for traversing event containers
    IteratorImp - A base class for implmenting the event iterators
    WindowIterator - An iterator for moving a window thru an event container
    ConstWindowIterator - A const window iterator
    Argument - A utility class to manage event function argument lists
    Value - Describes a column value or a constant event function
    Function - The base class of all event functions
    FunctionPtr - A smart pointer for event functions
    ColumnCache - A utility class to implement fast column access
    Column - An event function which returns a column value
    Info - An event function which returns layout information
    Math - Common math routines for the event functions
    MathOps - Defines the event function operators
    Condition - The basic class for event conditions
    ConditionPtr - A smart pointer for event conditions
    ConditionLogic - Defines the boolean operators between conditions
    Comparison - Defines a class to compare event functions
    ComparisonOps - Defines the relational operators
    Filter - A condition which acts as an event filter
    IfoSet - Describes a condition on the interferometer/detector mask
    Veto - A condition which acts as a veto
    Cluster - A condition which recognizes event clusters
    Algorithm - Conicidence, cluster analysis, etc.

A global event "factory" is used to manage the type and layout information of events. The following utility classes are used internally and generally not accessed directly:

    IndexList - A list with an index
    TypeInfo - Information record for a type
    LayoutInfo - Information record for the column layout
    Factory - Event factory

Header files are in here. Source files are in here.

Most of these classes are utility classes or are transparent to the user. Typically a user will use the Set container to load events from file; the Column class to pick out column values, the Filter and Veto conditions to preselect events and the select, coincidence and histgram algorithms of the Set class to analyze the events.

A typical session could look like this:

    #include "events/Events.hh"
    using namespace std;
    using namespace events;

      // Read some events from file
      Set events ("eventfile.xml");
      // Define an empty event set
      Set selected;
      // Define a second empty event set
      Set selected2;
      // Select a all burst events with amplitude > 5
      selected.Select (events, Filter ("burst::*") && Column("Amplitude") > 5);
      // Define a histogram (20 bins from 0 to 20)
      Histogram1 h1 (20, 0, 20);
      // Fill it with the amplitude of the selected events
      selected.Histogram (h1, Column("Amplitude"));
      // Select another set of events: all Burst which are not vetoed by
      // a PSL glitch
      // Important: use ConditionPtr rather than condition
      ConditionPtr cond_bursts = 
         Filter ("burst::*") && Column("Amplitude") > 5;
      ConditionPtr cond_pslglitch = 
         Filter ("glitch::psl") && (Column("Amplitude") > 10);
      // Set the veto window to -0.5 to +0.5 sec
      selected.Select (events, cond_bursts && !Veto (cond_pslglitch, -0.5, 1.0));
      // Again histogram it
      Histogram1 h2 (20, 0, 20);
      selected.Histogram (h2, Column("Amplitude"));
      // Or do it all in one
      ConditionPtr cond_vetoedbursts = 
         cond_bursts && !Veto (cond_pslglitch, -0.5, 1.0);
      Histogram1 h3 (20, 0, 20);
      events.Histogram (h3, Column("Amplitude"), cond_vetoedbursts);
      // Make sure we only have single H1 and L1 events
      events.Select (IfoSet ("H1L1", IfoSet::kAnyOne));
      // select all non-vetoed H1 burst
      selected.Select (events, IfoSet ("H1") && cond_vetoedbursts);
      // add all non-vetoed L1 burst
      selected.SelectAdd (events, IfoSet ("L1") && cond_vetoedbursts);
      // Now look for coincidences; set the coincidence window is 100 ms
      selected.Coincidence (0.1, IfoSet ("H1[0]") && IfoSet ("L1[1]"));
      // Make a 1D histogram of the time difference
      Histogram1 h4 (100, -0.1, 0.1);
      selected.Histogram (h4, Column("Time[1]") - Column("Time[0]"));
      // Or do it all in one
      ConditionPtr cond_H1 = IfoSet ("H1[0]") && cond_vetoedbursts;
      // Shift the condition to the second event
      ConditionPtr cond_L1 = IfoSet ("L1[1]") && Shift (cond_vetoedbursts);
      events.Histogram (h4, Column("Time[1]") - Column("Time[0]"), 
                       cond_H1 && cond_L1);
      // Check up on clusters
      Set clusters;
      clusters.SetWindow (10.0, -5.0);
      clusters.Clusters (2, 10);
      // Plot cluster number against time
      TSeries t1;
      clusters.TimeSeries (t1, Column("N"));

Written 2001, 2002 by Masahiro Ito and Daniel Sigg

alphabetic index hierarchy of classes

Please send questions and comments to sigg_d@ligo-wa.caltech.edu

generated by doc++