Run Pluto strips through pipeline

26 May 2002, asb

Initialize

<< asbGroup`

<< asbGroup`stripPipeline`

Set up directories and files, and auto function

datadir = "Macintosh HD:Remote Files:From \
aspen.lowell.edu:tocha:data2:boetiker:planet:agraph:";

outdirRef = "Macintosh HD:Local \
Files:Projects:Pluto:Strip Astrometry 200205:RUNS:20020522:";

outdir = "Macintosh HD:Local Files:Projects:Pluto:Strip Astrometry \
200205:RUNS:";

plutoEph = "Macintosh HD:Local \
Files:Reference:Ephemerides:plutoDE406_020501.020512.eph";

astromRefFile = 
        outdirRef <> "20020522.020.use";
targetRefFile = outdirRef <> "20020522.020.stars";

\(run[datestring_, 
        filenumber_] := 
Module[
{str}, 
str = 
          ToString[
            NumberForm[filenumber, 2, 
              NumberPadding -> {"0", ""}]]; 
Print["Working on file: ", 
          datestring <> "." <> str]; 
Print[
          Date[]]; 
res[
            number] = 
spAnalyzeStrip[
datadir <> 
              datestring <> ":" <> datestring <> "." <> 
              str <> ".fits",  (*FITS\ data\ file*) 
plutoEph, \  (*planet\ ephemeris, \ 
              for\ extraction\ purposes\ only*) 
astromRefFile, \  (*file\ of\ astrometric\ reference\ stars, \ 
              can\ be\ any\ *\(\(.\)\(astrom\)\(.\)\(dat\)\)\ file*) 
targetRefFile, \  (*file\ of\ target\ stars\ and\ planet, \ 
              can\ be\ any\ *\(\(.\)\(targ\)\(.\)\(dat\)\)\ file*) 
outdir <> datestring <> ":" <> 
              datestring <> "." <> 
              str, \  (*output\ directory\ and\ base\ name*) 
4, \ \  (*the\ centering\ star, \ 
              index\ number\ in\ astromRefFile*) 
{0.864653566288418`, \(-0.8497163417033069`\)}, \  (*row\ and\ col\ plate\ scale, \ 
              arcsec\ per\ pixel*) 
opPrintDiagnostics -> True, \  (*show\ status*) 
opPlanetName -> "Pluto", \  (*just\ to\ put\ in\ the\ *\(\(.\)\(targ\)\(.\)\(dat\)\)\ file*) 
opPlanetMagnitude -> 
              13.8, \  (*just\ to\ put\ in\ the\ *\(\(.\)\(targ\)\(.\)\(dat\)\)\ file*) 
opDoAstrometryCentroiding -> 
              True\  (*get\ centers\ for\ all\ astrometric\ reference\ stars*) 
]; \ 
Print[
          Date[]]
]\)

Sample runs: 20020522

Sample "normal" analysis of a strip

\(\(\( (*perform\ analysis\ on\ strip\ number\ 21\ from\ 20020522*) \)\(
\)\(run["20020522", 21]\)\)\)

If you want to re-do just the Pluto/target stars portion

Not quite as easy...  When doing it this way, the 3rd argument is changed from the standard astrometric reference file to the *.astrom.dat file for the strip.  It must have been already created in order for this to work.  This allows you to (for example) re-do the Pluto fit by changing the near-ness parameter or plate scale if necessary.

Date[]

datestring = "20020522";
str = "021";

\(\(spAnalyzeStrip[
datadir <> datestring <> ":" <> 
          datestring <> "." <> 
          str <> ".fits",  (*FITS\ data\ file*) 
plutoEph, \  (*planet\ ephemeris, \ 
          for\ extraction\ purposes\ only*) 
outdir <> 
          datestring <> ":" <> datestring <> "." <> 
          str <> ".astrom.dat", \  (*file\ of\ astrometric\ reference\ stars, \ can\ be\ any\ *\(\(.\)\(astrom\)\(.\)\(dat\)\)\ file*) 
targetRefFile, \  (*file\ of\ target\ stars\ and\ planet, \ 
          can\ be\ any\ *\(\(.\)\(targ\)\(.\)\(dat\)\)\ file*) 
outdir <> datestring <> ":" <> 
          datestring <> "." <> 
          str, \  (*output\ directory\ and\ base\ name*) 
4, \ \  (*the\ centering\ star, \ 
          index\ number\ in\ astromRefFile*) 
{0.864653566288418`, \(-0.8497163417033069`\)}, \  (*row\ and\ col\ plate\ scale, \ 
          arcsec\ per\ pixel*) 
opPrintDiagnostics -> 
          True, \  (*show\ status*) 
opPlanetName -> "Pluto", \  (*just\ to\ put\ in\ the\ *\(\(.\)\(targ\)\(.\)\(dat\)\)\ file*) 
opPlanetMagnitude -> 
          13.8, \  (*just\ to\ put\ in\ the\ *\(\(.\)\(targ\)\(.\)\(dat\)\)\ file*) 
opDoAstrometryCentroiding -> 
          False\  (*do\ not\ refit\ for\ centers\ of\ astrometric\ reference\ stars*) 
];\)\)

Date[]

Setting up to do several in a row

I'm leaving this running overnight.  It may crash Mathematica if it runs out of memory.

Do[
run["20020522", i], 
{i, 
        22, 34}]