Using my simulation stuff for multiple lightcurves

From: Robert A. Knop Jr. (robert.a.knop@vanderbilt.edu)
Date: Mon Oct 20 2003 - 09:44:50 PDT

  • Next message: Robert A. Knop Jr.: "Re: Using my simulation stuff for multiple lightcurves"

    See the message I sent before. I haven't looked at Vitaliy's scripts
    yet, so I'm not sure how hard it would be to modify them to get this in
    there; he could probably do it very quickly.

    Run fakeltcv2 as normal to generate your "first" lightcurve. See the
    previous instructions I sent out.

    Then, edit it and recompile it and run it again to generate the points
    in the "second" lightcurve. You may wish, for example, to change the
    "band" variable to "hst/f110w_nic2_shane.dat" to get the HST J-band
    point in there. Put your days and S/N points as usual. Make sure that
    if "obsdayisrest" is "false" that you are giving dates in the *observer*
    frame, and vice versa. Make the second fake lightcurve. Make another
    K-correction with makesnminuitK as usual (choosing the new "rest" and
    "observer" filters properly!). Then, in snminuit, do:

      snminuit,'fake1p2','two','uz','bj',file1='/path/to/lightcurve.1.filename',file2='/path/to/lightcurve.2.filename',kcorr1='/path/to/kcorrectioname.1',kcorr2='/path/to/kcorrectioname.2'

    You will want to make an "Output_zj" directory in your SNMinuit
    directory, as the output from this command will go there. The "1" files
    correspond to the lightcurve and K-corrections for your first (z-band)
    faked lightcurve, the "2" for your second (J or f110W) lightcurve.

    Please ask me to explain this better if this is not clear.

    I should probably write fakeltcv2 so that all of these things can be
    read from parameter files rather than coded at the top of the program,
    but for now editing the variables at the top of the code isn't too
    terrible.

    -Rob

    -- 
    --Prof. Robert Knop
      Department of Physics & Astronomy, Vanderbilt University
      robert.a.knop@vanderbilt.edu
    


    This archive was generated by hypermail 2.1.4 : Mon Oct 20 2003 - 09:45:06 PDT