From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- noao/digiphot/daophot/setimpars.cl | 276 +++++++++++++++++++++++++++++++++++++ 1 file changed, 276 insertions(+) create mode 100644 noao/digiphot/daophot/setimpars.cl (limited to 'noao/digiphot/daophot/setimpars.cl') diff --git a/noao/digiphot/daophot/setimpars.cl b/noao/digiphot/daophot/setimpars.cl new file mode 100644 index 00000000..57b1c57b --- /dev/null +++ b/noao/digiphot/daophot/setimpars.cl @@ -0,0 +1,276 @@ +# SETIMPARS -- Initialize the DAOPHOT task and pset parameters. + +procedure setimpars (image, restore, update) + +string image {prompt="Image name"} +bool restore {yes, prompt="Restore the last saved parameter set ?"} +bool update {yes, prompt="Update the last saved parameter set ?"} +bool review {no, + prompt="Review the parameters before saving ?"} +file parfile {"", prompt="Input algorithm parameters file"} +file datapars {"", prompt="The input data dependent parameters file"} +file findpars {"", + prompt="The input object detection parameters file"} +file centerpars {"", prompt="The input centering parameters file"} +file fitskypars {"", prompt="The input sky fitting parameters file"} +file photpars {"", + prompt="The input aperture photometry parameters file"} +file daopars {"", prompt="The input psf fitting parameters file"} +bool unlearn {no, + prompt="Unlearn the current algorithm parameters ?"} + +begin + # Define some temporary variables. + bool trestore, tupdate + string timage, tinparfile, toutparfile, tpars, tfile1, tfile2 + + # Read in the image name. + timage = image + trestore = restore + tupdate = update + print ("Setting parameters for image ", timage, " ...") + + # Set the input image name. + addstar.image = timage + allstar.image = timage + daoedit.image = timage + daofind.image = timage + group.image = timage + nstar.image = timage + peak.image = timage + pexamine.image = timage + phot.image = timage + psf.image = timage + pstselect.image = timage + substar.image = timage + + # Set the input coordinate / sky files back to their defaults. + phot.skyfile = "" + + # Set the input photometry files back to their defaults. + addstar.photfile = "" + addstar.simple_text = no + allstar.photfile = "default" + group.photfile = "default" + grpselect.ingroupfile = "" + nstar.groupfile = "default" + peak.photfile = "default" + phot.coords = "default" + psf.photfile = "default" + psf.pstfile = "" + pstselect.photfile = "default" + substar.photfile = "default" + + # Set the psfimage back to the default. + addstar.psfimage = "default" + allstar.psfimage = "default" + group.psfimage = "default" + nstar.psfimage = "default" + peak.psfimage = "default" + seepsf.psfimage = "" + substar.psfimage = "default" + + # Set the output photometry file names to the default. + allstar.allstarfile = "default" + allstar.rejfile = "default" + daofind.output = "default" + group.groupfile = "default" + grpselect.outgroupfile = "" + nstar.nstarfile = "default" + nstar.rejfile = "default" + peak.peakfile = "default" + peak.rejfile = "default" + phot.output = "default" + psf.groupfile = "default" + psf.opstfile = "default" + pstselect.pstfile = "default" + + # Set the output images back to the default. + addstar.addimage = "default" + allstar.subimage = "default" + daofind.starmap = "" + daofind.skymap = "" + psf.psfimage = "default" + seepsf.image = "" + substar.subimage = "default" + + # Set any output plot files back to the default. + phot.plotfile = "" + pstselect.plotfile = "" + psf.plotfile = "" + + # Get the input parameter file name. + tinparfile = parfile + if (tinparfile == "") { + if (access (timage // ".pars") && trestore) + tinparfile = timage // ".pars" + } else if (! access (tinparfile)) { + print ("File ", tinparfile, " does not exist ...") + return + } + + # Read in the input parameters. + + if (tinparfile != "") { + + print ("Reading algorithm parameters from file ", + tinparfile, " ...") + cl (< tinparfile) + + } else { + + tpars = datapars + if (access (tpars)) { + print ("Reading datapars parameters from file ", tpars, " ...") + tfile1 = mktemp ("tmp$pars") + tfile2 = mktemp ("tmp$pars") + dparam (tpars, > tfile1) + list = tfile1 + while (fscan (list, line) != EOF) { + if (substr (line, 1, 5) != "# EOF") + print ("datapars.", line, >> tfile2) + } + cl (< tfile2) + delete (tfile1 // "," // tfile2, verify-, default_action+, + allversions+, subfiles+, go_ahead+, >& "dev$null") + } else if (unlearn) { + print ("Reading default datapars parameters from disk ...") + unlearn ("daophot.datapars") + } else { + print ("Reading current datapars parameters from disk ...") + } + + tpars = findpars + if (access (tpars)) { + print ("Reading findpars parameters from file ", tpars, " ...") + tfile1 = mktemp ("tmp$pars") + tfile2 = mktemp ("tmp$pars") + dparam (tpars, > tfile1) + list = tfile1 + while (fscan (list, line) != EOF) { + if (substr (line, 1, 5) != "# EOF") + print ("findpars.", line, >> tfile2) + } + cl (< tfile2) + delete (tfile1 // "," // tfile2, verify-, default_action+, + allversions+, subfiles+, go_ahead+, >& "dev$null") + } else if (unlearn) { + print ("Reading default findpars parameters from disk ...") + unlearn ("daophot.findpars") + } else { + print ("Reading current findpars parameters from disk ...") + } + + tpars = centerpars + if (access (tpars)) { + print ("Reading centerpars parameters from file ", tpars, + " ...") + tfile1 = mktemp ("tmp$pars") + tfile2 = mktemp ("tmp$pars") + dparam (tpars, > tfile1) + list = tfile1 + while (fscan (list, line) != EOF) { + if (substr (line, 1, 5) != "# EOF") + print ("centerpars.", line, >> tfile2) + } + cl (< tfile2) + delete (tfile1 // "," // tfile2, verify-, default_action+, + allversions+, subfiles+, go_ahead+, >& "dev$null") + } else if (unlearn) { + unlearn ("daophot.centerpars") + print ("Reading default centerpars parameters from disk ...") + } else { + print ("Reading current centerpars parameters from disk ...") + } + + tpars = fitskypars + if (access (tpars)) { + print ("Reading fitskypars parameters from file ", tpars, + " ...") + tfile1 = mktemp ("tmp$pars") + tfile2 = mktemp ("tmp$pars") + dparam (tpars, > tfile1) + list = tfile1 + while (fscan (list, line) != EOF) { + if (substr (line, 1, 5) != "# EOF") + print ("fitskypars.", line, >> tfile2) + } + cl (< tfile2) + delete (tfile1 // "," // tfile2, verify-, default_action+, + allversions+, subfiles+, go_ahead+, >& "dev$null") + } else if (unlearn) { + unlearn ("daophot.fitskypars") + print ("Reading default fitskypars parameters from disk ...") + } else { + print ("Reading current fitskypars parameters from disk ...") + } + + tpars = photpars + if (access (tpars)) { + print ("Reading photpars parameters from file ", tpars, + " ...") + tfile1 = mktemp ("tmp$pars") + tfile2 = mktemp ("tmp$pars") + dparam (tpars, > tfile1) + list = tfile1 + while (fscan (list, line) != EOF) { + if (substr (line, 1, 5) != "# EOF") + print ("photpars.", line, >> tfile2) + } + cl (< tfile2) + delete (tfile1 // "," // tfile2, verify-, default_action+, + allversions+, subfiles+, go_ahead+, >& "dev$null") + } else if (unlearn) { + unlearn ("daophot.photpars") + print ("Reading default photpars parameters from disk ...") + } else { + print ("Reading current photpars parameters from disk ...") + } + + tpars = daopars + if (access (tpars)) { + print ("Reading psf fitting parameters from file ", tpars, + " ...") + tfile1 = mktemp ("tmp$pars") + tfile2 = mktemp ("tmp$pars") + dparam (tpars, > tfile1) + list = tfile1 + while (fscan (list, line) != EOF) { + if (substr (line, 1, 5) != "# EOF") + print ("daopars.", line, >> tfile2) + } + cl (< tfile2) + delete (tfile1 // "," // tfile2, verify-, default_action+, + allversions+, subfiles+, go_ahead+, >& "dev$null") + } else if (unlearn) { + unlearn ("daophot.daopars") + print ("Reading default daopars parameters from disk ...") + } else { + print ("Reading current daopars parameters from disk ...") + } + } + + # Review the current values for the algorithm parameters. + if (review) { + eparam ("datapars") + eparam ("findpars") + eparam ("centerpars") + eparam ("fitskypars") + eparam ("photpars") + eparam ("daopars") + } + + # Update the output parameter file. + toutparfile = timage // ".pars" + if (tupdate) { + if (access (toutparfile)) { + print ("Updating image parameter file ", toutparfile, " ...") + delete (toutparfile, verify-, default_action+, allversions+, + subfiles+, go_ahead+, >& "dev$null") + } else { + print ("Creating image parameter file ", toutparfile, " ...") + } + dparam ("datapars", "findpars", "centerpars", "fitskypars", + "photpars", "daopars", > toutparfile) + } +end -- cgit