VPFIT 2.02DR

 maximum file size      :   100000
 maximum chunk size     :   5200
 maximum no. parameters :   528
 maximum no. chunks     :   16
 exceeding above will corrupt results probably without warning
 

  vacuum wavelengths assumed
  Using data from : /home/rfc/atom.dat
 col den total after    :   v      see vp_setup.dat description

  *****************************************
  * ALL SIGMAS MULTIPLIED BY    1.50000   *
  * SEE VP_SETUP.DAT FILE                 *
  *****************************************

     5.00000E-02< b <   1000.000   see vp_setup.dat description
  Drop system if b <    5.00010E-02 & logN <    14.0000
   or if logN <    8.00000
   or if b >   1001.000

  options:   <CR> for previous value
  I - interactive setup and fit
  F - run from an input file
  D - display profiles from input file
  ? for help
  option (key) (key)...
> i ad il

 Add lines if prob(tot) < 0.0100   or prob(KS) < 0.0100
 Removes ill-conditioned lines
  lower case characters following values are treated as labels.  if the same
  upper case character follows another value then it will be tied to the lower
  case label.  isolated (un-paired) upper case characters fix values.

  setup: ? help, <CR> defaults, n,z,b,cs,sf,il,w,me,p,d,v to change
>

 Column density (n), logN (l) or emission (e), scalefactor
>
  Filename   1 or ? (lam,n) [     ]:
> q2233
  IRAF v2.10 wavelength coefficients
  Data is subset of larger one, wavelengths
   shifted by LTV1 =  -9624 channels
  sigma scale (midpoint value) :    1.50000
 Plot? [y]
>
 at co   co   3
 at er   co   2
 at ti   co   5
 plot parameter? (type he for options list)
> wm 3835
 wm 3835
  Channel  13937
 plot parameter? (type he for options list)
> wx 3860
 wx 3860
  Channel  15076
 plot parameter? (type he for options list)
>

 PGPLOT device? (? for list)
> xw
 Expand plot if needed:
 Cursor ("e" to mark edges, "q" when OK)
 Mark region in which data is to be fitted:
 left limit (space when ready)
 right limit
 Region limits:  3848.955A (channel 14573) -   3852.509A ( 14734)
  Channels   14573  14734
  Vacuum wavelength for start of chunk is     3848.96
 FWHM : 0.074A
  Line   1 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window
CIV 1550
 Wavelength used:    1550.77
 set cursor x- wavelength, y- base of feature
  .. and now half width at half (optical) depth
  Estimated vac wavelength & z are    3850.22    1.48278
     13.3463    9.27845    1.48278
  Line   2 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window
*         * = repeat last line ID
 CIV 1550
 Wavelength used:    1550.77
 set cursor x- wavelength, y- base of feature
  .. and now half width at half (optical) depth
  Estimated vac wavelength & z are    3850.61    1.48303
     13.3308    16.1174    1.48303
  Line   3 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window
*
 CIV 1550
 Wavelength used:    1550.77
 set cursor x- wavelength, y- base of feature
  .. and now half width at half (optical) depth
  Estimated vac wavelength & z are    3850.86    1.48319
     13.0455    5.60617    1.48319
  Line   4 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window
*
 CIV 1550
 Wavelength used:    1550.77
 set cursor x- wavelength, y- base of feature
  .. and now half width at half (optical) depth
  Estimated vac wavelength & z are    3851.17    1.48339
     13.6116    9.25177    1.48339
  Line   5 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window  Four CIV components from 1550 region

  Include other data? y, n (def),
   (or sf for more from the same file)
> y
  Filename   2 or ? (lam,n) [q2233]:
>
  Using filename: q2233               ...
  IRAF v2.10 wavelength coefficients
  Data is subset of larger one, wavelengths
   shifted by LTV1 =  -9624 channels
  sigma scale (midpoint value) :    1.50000
  Line, system number:
> 1548,3
 Line   1  wavelength   1548.19
  Wavelength    3844.46
     3317.6001107618    2.1959550678730D-02
  lo =  13798   hi =  14938
 FWHM : 0.074A
 Adding the previous lines to the displayfor the 1548 region display
   1 C IV      13.3463    1.48278    9.27845
   2 C IV      13.3308    1.48303    16.1174
   3 C IV      13.0455    1.48319    5.60617
   4 C IV      13.6116    1.48339    9.25177
 Convolved:  0.    3.13352E-02
 Plot? [y]
>
 at co   co   3
 at er   co   2
 at ti   co   5
 plot parameter? (type he for options list)
>

 Expand plot if needed:
 Cursor ("e" to mark edges, "q" when OK)
 Mark region in which data is to be fitted:
 left limit (space when ready)
 right limit
 Region limits:  3841.618A (channel 14239) -   3845.950A ( 14436)
  Channels   14239  14436
  Vacuum wavelength for start of chunk is     3841.62
 FWHM : 0.074A
  Line   5 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window
HI    Extra component present which is not seen in 1550
 Wavelength used:    1215.67
 set cursor x- wavelength, y- base of feature
  .. and now half width at half (optical) depth
  Estimated vac wavelength & z are    3843.10    2.16130
     13.0640    26.2072    2.16130
  Line   6 :ion, lamda0, N,b,z? <CR> to end
  .. in the graphics window

  Include other data? y, n (def),
   (or sf for more from the same file)
> n

 fort.26 summary to: f26.w3851 see vp_setup.dat description
 

 Summary to: f13.w3851

  Continue (c), fit (f), or stop (s)? [f]
> f

  no. of ions for fitting is   5
 

   ion          N           z           b         bturb   temp

 iteration   :   0 ( 1 )
 chi-squared :        4.575 (  1578.3702,  345 )

   C IV      13.3463    1.482781       9.28       0.00   0.00E+00  0    !  1
   C IV      13.3308    1.483030      16.12       0.00   0.00E+00  0    !  2
   C IV      13.0455    1.483189       5.61       0.00   0.00E+00  0    !  3
   C IV      13.6116    1.483393       9.25       0.00   0.00E+00  0    !  4
   H I       13.0640    2.161300      26.21       0.00   0.00E+00  0    !  5
 

 ......

 iteration   :   7 ( 1 )
 chi-squared :        1.765 (   608.9744,  345 )

   C IV      13.1366    1.482771       5.63       0.00   0.00E+00  0    !  1
   C IV      13.1503    1.482954      15.71       0.00   0.00E+00  0    !  2
   C IV      12.7645    1.483175       3.10       0.00   0.00E+00  0    !  3
   C IV      13.7834    1.483392      18.79       0.00   0.00E+00  0    !  4
   H I       13.0702    2.161266      26.66       0.00   0.00E+00  0    !  5

 parameter errors:
   C IV       0.0770    0.000005       1.01                             !  1
   C IV       0.0939    0.000019       4.66                             !  2
   C IV       0.1201    0.000007       1.77                             !  3
   C IV       0.0204    0.000006       1.11                             !  4
   H I        0.0446    0.000024       3.18                             !  5

 statistics for whole fit:
  Runs test  K-S test Chi-squared  Chans ndf   APr   Xp(.68) Xp(.95) Xp(.99)
   0.00000   0.00465    608.97     360  345   0.000  356.95  389.03  408.37

 statistics for each region :
    Start      End    Chi-squared  Chans df?
   3848.74   3852.71    335.69     162  150   0.000 < Prob < 0.000 g= 0.000   1
   3848.71   maxdev       2.9378              0.000   0.956   1.358   1.627

   3841.40   3846.17    273.28     198  183   0.000 < Prob < 0.000 g= 0.000   2
   3841.38   maxdev       1.7242              0.005   0.956   1.358   1.627
 

  Worst region is:  1   Fit is not good enough, so tries adding a line

   WLM:    3851.51 EW:    5.18753E-02

 Try adding C IV    1550.77 at z= 1.4836128 b= 12.3 logN=  13.13

  All old lines lost

.....           In this case, has four attempts

 iteration   :   2 ( 4 )
 chi-squared :        1.625 (   545.8841,  336 )

   C IV      13.1168    1.482770       5.41       0.00   0.00E+00  0    !  1
   C IV      13.1817    1.482953      17.31       0.00   0.00E+00  0    !  2
   C IV      13.0164    1.483189       5.57       0.00   0.00E+00  0    !  3
   C IV      13.6509    1.483372      11.01       0.00   0.00E+00  0    !  4
   H I       13.0695    2.161270      26.48       0.00   0.00E+00  0    !  5
   C IV      12.4353    1.483500       1.92       0.00   0.00E+00  0    !  6
   C IV      12.5093    1.483545       8.82       0.00   0.00E+00  0    !  7
   C IV      12.8438    1.483572      13.02       0.00   0.00E+00  0    !  8

 iteration   :   3 ( 4 )
 chi-squared :        1.623 (   545.3870,  336 )

   C IV      13.1166    1.482770       5.41       0.00   0.00E+00  0    !  1
   C IV      13.1813    1.482953      17.28       0.00   0.00E+00  0    !  2
   C IV      13.0196    1.483189       5.63       0.00   0.00E+00  0    !  3
   C IV      13.6511    1.483372      11.00       0.00   0.00E+00  0    !  4
   H I       13.0695    2.161270      26.48       0.00   0.00E+00  0    !  5
   C IV      12.4497    1.483501       1.55       0.00   0.00E+00  0    !  6
   C IV      12.5164    1.483548       8.83       0.00   0.00E+00  0    !  7
   C IV      12.8504    1.483573      12.40       0.00   0.00E+00  0    !  8

 parameter errors:
   C IV       0.0930    0.000005       1.14                             !  1
   C IV       0.1142    0.000022       6.07                             !  2
   C IV       0.1290    0.000009       1.74                             !  3
   C IV       0.1557    0.000023       3.28                             !  4
   H I        0.0427    0.000023       3.04                             !  5
   C IV       1.1997    0.000017      11.31                             !  6
   C IV      11.6805    0.000295      76.73                             !  7
   C IV       4.9466    0.000357      22.74                             !  8

 statistics for whole fit:
  Runs test  K-S test Chi-squared  Chans ndf   APr   Xp(.68) Xp(.95) Xp(.99)
   0.00000   0.02942    545.39     360  336   0.000  347.79  379.46  398.57

 statistics for each region :
    Start      End    Chi-squared  Chans df?
   3848.74   3852.71    292.98     162  141   0.000 < Prob < 0.000 g= 0.000   1
   3848.71   maxdev       2.8366              0.000   0.956   1.358   1.627

   3841.40   3846.17    252.41     198  174   0.000 < Prob < 0.005 g= 0.000   2
   3841.38   maxdev       1.9641              0.001   0.956   1.358   1.627

  Plot? y,n, or c=change device [y]
>
  Identify which dataset:
  # ID   Filename       data cont    lamcoeffs      nchans  KS prob
  1  1    q2233            1   0    3528.96  0.0221  37879    0.000  CIV1550 region
  2  2    q2233            1   0    3528.96  0.0221  37879    0.001  CIV1548 region
 Which # ? [ 1 ]
> 1
  Line, system number: (* for internal guesses)
> *
 at co   co   3
 at er   co   2
 at ti   co   5
 plot parameter? (type he for options list)
>

  Plot? y,n, or c=change device [y]
>
  Identify which dataset:
  # ID   Filename       data cont    lamcoeffs      nchans  KS prob
  1  1    q2233            1   0    3528.96  0.0221  37879    0.000
  2  2    q2233            1   0    3528.96  0.0221  37879    0.001
 Which # ? [ 1 ]
> 2
  Line, system number: (* for internal guesses)
> *
 at co   co   3
 at er   co   2
 at ti   co   5
 plot parameter? (type he for options list)
>

  Plot? y,n, or c=change device [y]
> n
  Restore continuum? [n]
> n
  Fit more lines? [n]
> n
[exits]

In this case the adding lines iterations terminated because there was no improvement in the normalized chi-squared. A look at the spectrum reveals why -- in the CIV 1550 region the red wing of the complex does not match the CIV 1548 very well. Since both lines are in the Ly-a forest of the quasar, the way to get a good fit would be to add a Ly-a in the red wing of the first region, rather than trying with more and more CIV components.