Документ взят из кэша поисковой машины. Адрес оригинального документа : http://zebu.uoregon.edu/~uochep/talks/talks02/G020237.pdf Дата изменения: Thu Sep 26 00:48:33 2002 Дата индексирования: Tue Oct 2 10:35:59 2012 Кодировка: Поисковые слова: п п п п п п п п п п п п п п п п п п п п п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п

Event Analysis Tool
Masahiro Ito (U of Oregon), Daniel Sigg (LHO) 3/22/2002
LIGO-G020237-00-D


Purpose
- Provide various analysis algorithms and make the analysis easier.

Overview
- C++ Package.
- A part of the gds tree. - gds/Triggers/events

- Run under ROOT.
- ROOT command line environment. - Stand alone executable.

- Flexible column layouts.
- Each event contains layout information.

- Import event data from LIGO database.
- Import/Export XML event tables.

- Do analysis in a single line of code.
- Select, Sort, Coincidence and Cluster. - Plot histogram and time series.


Flow Chart
GUILD & LIGO Tools Local File (XML) Restore Events Import XML event data into an event container

LDAS Meta Database

Save Events

Process events (Sort, Coincidence, Cluster)

Create Histogram and/or Time Series

Plot Histogram/ Time Series


Event Containers Hierarchy
"Chain" container
- A set of List containers. Example :

"Set" container
- A set of Chain containers. - Most of operations can be done by handling Set containers.

Chain - 1

Chain-1 : H0 events Chain-2 : L0 events .........

Chain - 2

Chain - 3

"List" container
- A set of Event containers. Example:
List-1 : Jan. events List-2 : Feb. events .........

Set
......

"Event" container
- The smallest element of event containers. - Holds one event.
List - 1

Chain - N

List - 2

......

List - M

Event - 1

Event - 2

......

Event - X


Event
- Column layout information.
Type. Fixed columns (Time , Name and IFO). Variable columns. Can be extended.

- Column data.
- Int, Real, Complex, Time, String and Event.
Event

Layout Data

List of Layouts

Layout Information

Column Information

Layout - 1 Layout - 2

Type Column - 0 Column - 1 ...... Column - N

Name Var. Type Offset

Value - 0
Data Block

Layout - 3 Layout - 4 ...... Layout - M

Value - 1 ...... Value - N


Restore & Save Events
- XML Files
- Support read/write LIGO-LW format. - Built in layout types(see T990101) : GDS_Trigger, Sngl_Inspiral, Sngl_Burst, Sngl_Ringdown, Sngl_Unmodeled and Sngl_DPeriodic.

- ASCII Files
- Support only importing GDS Trigger format generated by GUILD.


Processing Events
- Select
- Select events by specifying conditions.

- Sort
- Sort events by a column or given function of columns.

- Coincidence
- Find coincident events which fulfill conditions. - Capable of 2 to N multi-coincidences.

- Cluster
- Find cluster events which satisfy certain conditions, time window and threshold.
Event

Sliding Time Window - Events are not independent.
Time Window Time

Cluster of more than 2 events = 3 Coincidence of 2 events = 4


Make Histogram & Time Series
- 1-D and 2-D Histogram
- Fill a DMT histogram container with a specified column or value of an event function. - Plot histogram with DMT graphics or ROOT graphics libraries.

- Time Series
- Make a Time Series (DMT TSeries container) of the event rate, values of specified column or event function. - Plot time series with DMT graphics libraries.


Example 1
using namespace events; Set eventset("pemLVEA.xml"); // Import events from pemLVEA.xml. TSeries tRate; // Make a time series tRate.setName("Event Rate"); // Give eventset.TimeSeries(tRate, true); // // Plot(tRate); // Plot time series. container. it a name. Generate time series of event rate.

Histogram1 hSize("SIZE", 30, 5.0, 20.0, "Size", "# of Events"); // Make a histogram container. eventset.Histogram(hSize, Column("SIZE"), Column("DURATION") < 5.0); // Fill hSize. PlotHistogram(hSize); // Plot histogram.


Example 2
eventset.AddChain("pemEX.xml"); // Import another set of events // from pemEX.xml. eventset.Sort(); // Sort events in chronological order.

eventset.Coincidence(TimeWindow(20.0, -10.0)); // Find coincidence events and store the result.

Histogram2 h2Size("Size", 20, 0.0, 20.0, 20, 0.0, 20.0 "Size-1", "Size-2"); eventset.Histogram(hDiff, Column("SIZE(0)"), Column("SIZE(1)")); PlotRootHistogram(h2Size, "col");


Event Function & Condition
eventset.Histogram(hSize, Column("SIZE") , Column("DURATION") < 5.0 ); eventset.Histogram(h2Size, Column("SIZE(0)") , Column("SIZE(1)") );

-

Event Function
Value : Represents a constant value. Column : An access to column value of

- Operators : ! , && , || , == , != , < , <= , > , >= , WithIn. - Filter : An event filter by Type and Name information. - Veto/Coincidence : Check for nearby events. - IFO Set : Select events by IFO. - Cluster : Count nearby events.

Event Condition

an event. Info : An access to layout information

of an event. - Math Operation : + , * ,& , << , sqrt , sin , conj and more.


Conclusion
- It's available. Use it. - Add String operations.
- Support regular expression.

- User feedback. - Bug fixes. - Support for time segments. - Event look up through Lidax. - Event classification by DMT Trigger environment.