Whipple data analysis at Purdue is performed with a set of programs designed to apply step-wise image parameterization and Supercuts. These steps are outlined below. The programs make use of several external libraries, including:
CERNLIB
HDF
STARLINK
MATLAB
fz2hdf
(CERN,HDF,[MATLAB]) - conversion of ZEBRA data
format to a reduced HDF format. Reduced data format example:
gt009630.hdf .
Several diagnostic distributions are written to a MATLAB
file yyxxxxxxd.mat. See mdiagnose below for explanation.
USAGE:
fz2hdf [options] yyxxxxxx.fz
where [options] are
-h help
-o %f oscillator freq
[MHz] pre FZ version 76 [ns/s] post version 76
outputs yyxxxxxx.hdf.
Example: fz2hdf gt008899.fz
fz2hdf -o 1.0 gt008899.fz
cpeds (HDF)
- performs ADC pedestals calibration. Pedestals and pedestal
variances are written to the ASCII database file hrcyy.cpeds
(yy = year).
USAGE:
cpeds [options] yyxxxxxx.hdf
where [options] are
-e 11m data
-g %s path to database files [/usr/local/db/]
-h help
-l %f lower pedestal variance threshold [2.00]
-u %f upper pedestal variance threshold [7.00]
Example: cpeds -g ./ gt008899.hdf
cn2gains (HDF)
- performs gain calibration for nitrogen pulser runs. Gain
factors for each PMT are written to the ASCII database file
hrcyy.cn2gains.
USAGE:
cn2gains [options] yyxxxxxx.hdf
where [options] are
-e 11m data
-g %s path to database files [/usr/local/db/]
-h help
cparam (HDF,STARLINK)
- performs image parameterization. Hillas parameters are written to
an HDF file. Example:
gt009630p.hdf.
USAGE:
cparam [options] yyxxxxxx.hdf
where [options] are
-b %f,%f picture,boundary threshold [4.25,2.25]
-d %f declination of source [ddmmss.s]
-e 11m data
-f fixed two dimensional analysis, no derotation
-g %s path to database files [/usr/local/db/]
-h help
-i %s padding pair id
-n %s nitrogen id
-p software padding
-r right ascension of source [hhmmss.s]
-s %s source name [UNKNOWN]
-t %d number of trigger tubes [91]
-v verbose text output of all events
-w %f elongation factor [1.85]
Example: cparam -g ./ -p yyxxxxxx.hdf
mcut (HDF,MATLAB)
- performs Supercuts and outputs MATLAB binary files of events
passing cuts and alpha, radial and 2d distributions.
USAGE:
mcut [options] gxyyyyyyp.hdf
where gxyyyyyyp.hdf -> parameterized file and
[options]
%d is an integer %f is a float (real)
-a %f -> alpha cut [15.00]
-b %f,%f,%f-> bin size alpha,radial,image [5.00,0.05,0.10]
-c %f,%f -> aperture center [0.00,0.00]
-d %f,%f -> distance cuts [0.51,1.20]
-f %f -> frac3 cut [0.975]
-g -> Gaussian smoothing [TRASH CAN]
-h -> help
-k %f -> run duration [28.00]
-l %f,%f -> length cuts [0.16,0.33]
-n %d -> number of phase bins [20]
-p %f,%f,%f,%f,%f -> periodic analysis ra,dec,freq,freq',epoch
-r %f -> aperture size [0.32]
-s %d,%d -> size cut [0,999999]
-t %d,%d,%d-> software trigger [80,40,0]
-w %f,%f -> width cuts [0.07,0.16]
Example: mcut -s 400,999999 -t 100,80,0 -g yyxxxxxxp.hdf
mtabout (MATLAB)
- MATLAB program which tabulates
results and produces an output summary page. Output summary page
example: J0751 .
USAGE:
mtabout(file_list,'PropertyName','ProperyValue')
The first argument is the file name of the files list and is
required. The next arguments are optional Properties as list below.
option 'Grid', values ['RaDec','Camera'] - default 'RaDec'
option 'Epoch', - default 2000.0
option 'Mode', values ['Pair','Tracking'] - default 'Pair'
option 'Plot', values ['Image','Spectrum','LightCurve','OPhase','GPhase']
- default 'Image'
option 'AlphaRatio', - default 0.334
option 'RadialRatio', - default 0.162
Example: >>mtabout('files','Mode','Tracking');
canalyze - is a BASH program which runs
all of the above software automatically given a pair directive
input file in the form mode on off nitrogen date.
USAGE:
canalyze [options] list
[options]
-c #tubes camera configuration for cuts [331]
-g path path to database files [/usr/local/db]
-h help
-l low elevation cuts
-m manual mode - cdb required
-n nice mode - remove fz files
-s output command script instead of executing
where list is a file containing the pair directives. Example:
pr gt008898 gt008899 gt008895 971201
pr gt008900 gt008901 gt008895 971201
tr gt008904 none gt008903 971202
Options may now be included for canalyze at any location in the file.
# - comment tag ignores rest of line
(# followed by a least one space)
source %s - source name
pulsar %s - pulsar catalog name - will epoch fold
camera %d - camera configuration [#pixels]
position %d %d - source position ra [hhmmss] dec [ddmmss]
onoff - onoff mode for off position calculation
offon - offon mode for off position calculation
option mcut %s - additional cuts to be passed to mcut
option cparam %s - additional cuts to be passed to cparam
Example:
source Crab
camera 109
position 53421 220048
#
# Crab 1998 - full 331 trigger
#
option cparam -t 331
option mcut -t 54,52,0 -w 0.07,0.18 -l 0.16,0.37 -d 0.51,1.18
# option mcut -t 80,45,0 -w 0.073,0.16 -l 0.16,0.50 -d 0.51,1.3 -a 10.0
# CFD (85 mV)
pr gt009930 gt009931 gt009932 980130
pr gt009967 gt009968 gt009969 980202
# CFD (95 - 110mV) optimizing trigger
pr gt010019 gt010020 gt010022 980222
pr gt010036 gt010037 gt010039 980223
# CFD 100mV
pr gt010057 gt010058 gt010063 980228
pr gt010060 gt010061 gt010063 980228
offon
pr gt010078 gt010077 gt010083 980301
onoff
pr gt010079 gt010080 gt010083 980301
pr gt010091 gt010092 gt010093 980302
# adjusted summed output levels
pr gt010157 gt010158 gt010163 980322
You should find the scipt option, -s, useful for doing none
standard analysis. With this option, canalyze will output
to the screen all the commands it would have executed. This
output may then be re-directed to a file and edited for custom
use.
SUN Workshop Pro (C,FORTRAN,etc)
HDF(Hierarchical Data Format)
lessard@physics.purdue.edu
