DST analysis Primer

DST analysis Primer

Y. Unno Physics, KEK

Making symbolic link to the DST data

It is not wise to make a hardcopy of the DST data. Instead, it is wise to make a symbolic link in your working directory. The UNIX command is 
>ln -s other/directory/name yourfilename

E.g., 

>ln -s /afs/cern.ch/user/w/wmurray/public/tb95/paw/lh2035.1.dst lh2035dst

Generating a (function) program with PAW

When using 'ntuple' in the paw, the paw has an built-in capability to generate a (function) program. This is for the program to use the same common block to pass the data from the paw to the user routine. A user can implement his/her code in the program.

In order to get the ntuple structure, you have to get in the paw and set the directory to the ntuple... Here is how. 

>paw
>hist/file 1 filename             <- filename is the ntuple file
>cdir //lun1/aaaa                 <- aaaa is the ntuple name
>nt/uwfunc 10 funcname ""         <- don't miss "", which generate a fuction
                                     with the name 'funcname'
                                     10 is the ntuple id #
>exit

Then, the file named 'funcname' is generated in your directory. Edit the 'funcname' to implement the code whatever you want.

will generatea program, named 'generic.f', for the 95 H8 strip beamtest 
>paw
>exec generic
>ls                   -> you will see 'generic.f' in your directory.

To run the program

After completing the file 'funcname', it can be executed to process the ntuple. 
>paw
>hist/file 1 filename
>cdir //lun1/aaaa
>nt/loop 10 funcname nnnn         <- nnnn is the number of events to be 
                                  processed.

Examples:

I have generated two examples to process the DST. One is to dump the events in the DST. By running this program you can get info on the content of the DST. Also, the program can be a primitive of your program. The other program is to get the correlation 2D map of the telecope hits. This can be another starting point of your analysis program, too. So, look at

There is also a kumac file to run the tele.f program. You will learn how to run the program out of this kumac content. See

which is 
* telescope plane correlation
h/file 1 lh2035dst
cd //lun1/strp
nt/loop 10 tele.f 1000
cd //pawc
exec meta#open                  <- meta.kumac to open paw.ps
zone 2 2
h/plot 102
h/plot 103
h/plot 107
h/plot 108
exec meta#close                 <- meta.kumac to close paw.ps

To execute, you can do 

>paw
>exec example
>exit

Then, a ps file, paw.ps, is generated for the print-out. The meta.kumac file opens and closes the postscript file, paw.ps.

Other information

Other infos such as the DST logs, etc., can be found on WWW. A list of pointers are listed in the beamtest page (of KEK). The page is

You can reach the info via the pointers directly from the ATLAS page of www.cern.ch -> subdetectors and working group -> inner detector -> Barrel sector prototype -> offline analysis -> offline description -> DST's...

Send your comments of this page to Y. Unno