IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 42062


Ignore:
Timestamp:
Feb 24, 2022, 3:50:31 PM (4 years ago)
Author:
eugene
Message:

explicitly define MIN/MAX Error vs MIN/MAX mag

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/eam_branches/ipp-20211108/ippTests/ippCheck/smfcheck.pro

    r42058 r42062  
    676676  book getword psfFields $field MIN_MAG  -var minMag
    677677  book getword psfFields $field MAX_MAG  -var maxMag
     678  book getword psfFields $field MIN_ERR  -var minErr
     679  book getword psfFields $field MAX_ERR  -var maxErr
    678680  # echo $field $useStat
    679681
    680682  # this variable tells the name of the vector which
    681683  # should be used to select the values to test
    682   # the choice of minMag / maxMag depends on this value
     684  # the choice of minErr / maxErr depends on this value
    683685  $keepFieldName = PSF_INST_MAG_SIG
    684686
     
    686688  if ("$useStat" == "VALUE")
    687689    set dtmp = $statScale*(tmp_$v2 - $field\_$v1)
    688     subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minMag) && ($keepFieldName\_$v1 < $maxMag))
     690    subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minErr) && ($keepFieldName\_$v1 < $maxErr))
    689691  end
    690692  if ("$useStat" == "IVALUE")
    691693    set dtmp = tmp_$v2 - $field\_$v1
    692     subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minMag) && ($keepFieldName\_$v1 < $maxMag) && (tmp_$v2 != -4096) && ($field\_$v1 != -4096))
     694    subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minErr) && ($keepFieldName\_$v1 < $maxErr) && (tmp_$v2 != -4096) && ($field\_$v1 != -4096))
    693695  end
    694696  if ("$useStat" == "RATIO")
    695697    set dtmp = $statScale(tmp_$v2 / $field\_$v1)
    696     subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minMag) && ($keepFieldName\_$v1 < $maxMag))
     698    subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minErr) && ($keepFieldName\_$v1 < $maxErr))
    697699  end
    698700  if ("$useStat" == "RVALUE")
    699701    set dtmp = $statScale*(tmp_$v2 - $field\_$v1) / $field\_$v1
    700     subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minMag) && ($keepFieldName\_$v1 < $maxMag))
     702    subset dtmpS = dtmp where (($keepFieldName\_$v1 > $minErr) && ($keepFieldName\_$v1 < $maxErr))
    701703  end
    702704  if (("$useStat" == "NONE") || ("$useStat" == "NULL"))
Note: See TracChangeset for help on using the changeset viewer.