STAR OO data model and 3D visualization

1
STAR OO data model and 3D visualization

• by Valery Fine
• STAR collaboration meeting
• August, 1999

2
OO model of the STAR simulation / reconstruction
chain
St_DataSet object the "named" list of
St_DataSet objects
St_Table St_TableSorter
Data definition
base container class
St_ObjectSet
abstract TObject
St_DataSet St_DataSetIter
file system description
St_FileSet
GEANT Geometry structure
St_Node/ St_NodeView
Flow control
StMaker
3
OO model of STAR geometry

• Detector geometry supplied by St_geant_Maker
  (GEANT 3.21)
• Event geometry supplied by bfc.C Makers of
  the reconstruction chain.

4
Detector and Event data model
STAR 3D geometry Object Model

St_DataSet Legend
"classA

"classB" is derived from
v v the
"classA
St_Node "classB

"classQ

"classQ" has a pointer
TShape to
"classT
"classT
v -------------------
------------------------------------------


. . .

v
v v v v
St_PolyLineShape TBRIK TTUBE
TPCON TTRD1
R O O T S H A
P E S v
(see GEANT 3.21 shapes as well) St_PolyLine3D
TPoints3DABC
v
--------------------------------------------------
------------



V
v v
v StHits3DPoints StHelix3DPoints
St_Table3Points St_Points3D




StObjArray StTrack / StHelixD
St_TableSorter flat floating

point array
(see St_PolyLine3D as well)

St_Table S T A R S H A P
E S
Detector geometry
Event geometry
5
StEventDisplayMaker
Detector geometry St_Node / St_NodeView (St_gea
nt_Maker/StDBMaker)
Event geometry StEvent / St_Tables ( bfc
Makers)
StEventDisplayMaker
StVirtualEventFilter
3D Viewers TBrowser TPad X3D OpenGL
6
TBrowser viewer
Full GEANT tree (20000 nodes)
sectors only tree (24 nodes)
St_NodeView
St_Node
7
X3D / TPad viewers
Run 3002 27th-July-99
8
OpenGL viewer
9
How to use this (with bfc.c) ?
Add display parameter to turn graphics on
root4star 'bfc.C(2,"tdaq display","/disk1/star/daq
/990729.3104.01.daq")'where the 3d parameter,
the file name, can be either daq file found on
Web site
http//duvall.star.bnl.gov/cgi-bin/prod/dbrunlog.p
l?showyesrun
The command above reads first 2 events from the
DAQ file with no stop. To get any next event (if
present) one should type root4star1
chain-gtMakeEvent() from the ROOT4STAR command
prompt.
10
Typical StMaker list of bfc.C
QAInfo 7 ER99
is ON Turn on Year 1a geometry and 1999
engineering run real QAInfo 10
Y2a is ON Turn on Year 2a
geometry (and corresponding Makers)
QAInfo 13 XIN
is ON Read XDF input file with g2t
QAInfo 16 TDAQ
is ON TPC DAQ chain QAInfo
19 FieldOn is ON Use nominal STAR
field QAInfo 22 TPC
is ON TPC in Chain
(St_tcl_tpt_Maker) QAInfo
27 TCL is ON St_tcl_Maker
in Chain QAInfo 28 TPT
is ON St_tpt_Maker in
Chain QAInfo 29 SVT
is ON SVT in Chain
(St_srsstk_Maker) QAInfo
30 SRS is ON St_srs_Maker
in Chain QAInfo 31 STK
is ON St_stk_Maker in
Chain QAInfo 32 FTPC
is ON FTPC in Chain
(St_fclfpt_Maker) QAInfo
34 FCL is ON St_fcl_Maker
in Chain QAInfo 35 FTP
is ON St_fpt_Maker in
Chain QAInfo 36 EMS
is ON St_ems_Maker in Chain
QAInfo 38 TRG
is ON Trigger in Chain
(St_ctfmwctrg_Maker) QAInfo
39 CTF is ON St_ctf_Maker
in Chain QAInfo 40
MWC is ON St_mwc_Maker
in Chain QAInfo 43 GLOBAL
is ON GLOBAL in
Chain QAInfo 44 MATCH
is ON StMatchMaker in Chain
QAInfo 45 PRIMARY
is ON StPrimaryMaker in Chain
QAInfo 46 V0
is ON StV0Maker in Chain
QAInfo 47 XI
is ON StXiMaker in Chain
QAInfo 48 KINK
is ON StKinkMaker in Chain
QAInfo 49 DST
is ON St_dst_Maker in Chain
QAInfo 50 EVENT
is ON StEventMaker in Chain
QAInfo 51 ANALYSIS
is ON StAnalysisMaker in Chain
QAInfo 52 QA
is ON St_QA_Maker in Chain
QAInfo 53 TREE
is ON StTreeMaker in Chain
QAInfo 55 DISPLAY
is ON StEventDisplayMaker in Chain
QAInfo 57 DEFAULT
is ON Default has been set
QAInfoInput file
name /disk1/star/daq/990727.3002.01.daq with
No. of Events to process 1 QAInfoOutput root
file name 990727.3002.01.root

11
How to see that (with doEvents.C )
Tracks
StHelix
vertices
StVertex
.x drawEvent.C(1,"/disk00000/star/test/dev/tfs_Lin
ux/Fri/year_1b/psc0050_01_40evts.dst.root","")
12
Pad control panel
.x drawEvent.C(1,"/disk00000/star/test/dev/tfs_Lin
ux/Fri/year_1b/psc0050_01_40evts.dst.root","")
Try to recognize such case without visualization
.x PadControlPanel.C
13
From TPad to OpenGL view
.x drawEvent.C(1,"/disk00000/star/test/dev/tfs_Lin
ux/Fri/year_1b/psc0050_01_40evts.dst.root","")
Select View with OpenGL
Note Free MesaGL viewer works quite slow over
TCP/IP, one has to open GL window as small as
possible to get the reasonable response time.
14
Filterslist of the default filters provided
Filter for Primary Vertex is OFF
Filter for Tpc Hit is OFF Filter
for Svt Hit is OFF Filter for
Ftpc Hit is OFF Filter for
EmcTower Hit is OFF Filter for
EmcPreShower Hit is OFF Filter for
SmdPhi Hit is OFF Filter for SmdEta
Hit is OFF Filter for Vertices is
ON Filter for Global Tracks is ON Filter for
Track is ON Filter for Track Tpc
Hits is OFF Filter for Track Svt Hits is
OFF Filter for Track Ftpc Hits is
OFF Filter for St_Table generic object is
OFF
15
Example of the user-defined filter
root4star session with the simple user-defined
filter changes the previous picture
significantly, leaving only 15 tracks of 7039
ones
Int_t St_TLA_EventFilterFilter(StGlobalTrack
globTrack,Width_t ,Style_t ) // This
filters the short tracks (lt 250 cm) off Int_t
color StVirtualEventFilterFilter(globTrack,siz
e,style) if (GetFlag() 2
globTrack-gtlength() lt 250) color 0 return
color
root4star 1 .x drawEvent.C(1,"/disk00000/star/te
st/dev/tfs_Linux/Fri/year_1b/psc0050_01_40evts.dst
.root","") root4star 2 disp-gtSetTrackFlag(2) roo
t4star 3 disp-gtReDraw
16
Filters - TBrowser and command line interfaces
Click here to redraw the picture
disp-gtReDraw()
Click here to turn filter on/off
disp-gtSetSmdPhiHitFlag(1)