Table of Contents
pvtrace - PV diagram envelope tracing
pvtrace in=image [parameter=value]
pvtrace tries to recover the rotation curve from a single
position-velocity (PV) diagram, by tracing the outer edges of the signal
that, after correction for a velocity dispersion factor, is the material
with the maximum amount of tangential motion. This is particularly appropriate
for high-inclination galaxies, where significant amounts of gas is observed
below the rotation speed, due to line-of-sight integrations and simple moments
in velocity or gaussian fits fail to recover the appropriate rotation speed.
The rotation velocity, V, is determined from a terminal velocity V_t in
the usual way:
V = (V_t - V_sys)/sin(i) - sigma
but corrected for a dispersion factor sigma. The terminal velocity is defined
by the velocity at which the intensity becomes equal to
I_t = sqrt( (eta*I_max)^2 + (I_lc)^2 )
where I_max is the maximum intensity in the map, and I_lc the lowest counter
(usually 3 times the noise in the PV diagram). Both eta and I_lc can be
set via parameters to the program. Typical values for eta = 0.2 - 0.5.
In addition,
this program will also perform the usual peak fit, gaussian fit, and simple
first moment analysis. See OUTPUT below for the output format. For lower
inclinations, the envelope tracing method should give simular results to
the other methods.
- in=
- Input image file - must be an PV image(5NEMO)
with position along the 1st axis, and velocity along the 2nd. Units are
not important.
- eta=
- [0..1] weight factor for I_max adding to I_lc. [Default:
0.2]
- ilc=
- Lowest Contour value. Normally taken as something like 3 times
the noise level. [Default: 0]
- sign=-1|1
- Designates rotation of the galaxy.
This is only important for tracing the 'right' part of the PV diagram. sign=1
is meant for PV diagrams where the approaching side is for positive positions.
By using the "wrong" value of sign you can find out something about the
radial extent of the tracer, as for large radii this would cause small
values of dV/dR. For a known rotation curve using sign=1 one could determine
the radial extent from sign=-1. [Default: 1]
- sigma=
- Velocity dispersion correction
factor. This includes both the observational resolution as well as the velocity
dispersion of the tracer component that is being observed. [Default: 0]
- clip=
- Clipping values. If value(s) given, signal will be clipped as follows:
if one value given, clipping occurs between -clip and clip. If two numbers
given, clipping occurs between clip(1)
and clip(2)
. Default: None
- vsys=
- Systemic velocity of the galaxy. Since the velocity axis (assumed to be
the 2nd or Y-axis), this number is simply subtraced from each velocity determination.
[Default: 0]
- inc=
- Inclination of the galaxy. Each rotation speed is corrected
by sin(inc). [Default: 90]
- center=
- Center of the galaxy along the Position
axis. [Default: 0]
Output of the program will be an ASCII table,
with various derived velocities as function of position along the P axis
in the PV-diagram. The following columns are returned:
1 position from the P-axis of the PV diagram
2 envelope tracing (corrected with vsys, sini, and sigma)
3 first moment (sum(I*v)/sum(I)) (corrected with vsys, sini)
4 velocity at peak of signal
5 velocity at peak of parabolic fit of 3 points at peak signal
Position axis must be carefully created such that the galaxy center
is at P=0, whereas for velocity it can be corrected with VSYS=.
rotcur(1NEMO)
,
ccdvel(1NEMO)
, snapgrid(1NEMO)
, ccdmom(1NEMO)
Sofue, Y. 1996 ApJ 458, 120
Sofue et al 1997 PASJ 49, 17
Sofue et al 523, 136
src/image/rotcur sources
Peter Teuben
5-May-01 V1.0 Created PJT
6-may-01 V1.1 clip=, interpolated to find trace PJT
Table of Contents