IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 5143


Ignore:
Timestamp:
Sep 27, 2005, 7:43:11 AM (21 years ago)
Author:
eugene
Message:

fixed millimag / mag issues

Location:
trunk/Ohana/src/relphot
Files:
1 deleted
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/relphot/src/GridOps.v2.c

    r4865 r5143  
    293293void plot_grid (Catalog *catalog) {
    294294
    295   int i, j, m, c, n, N;
     295  int i, j, m, c, n, N, Narea;
    296296  float Msys, Mrel, Mcal, Mmos;
    297297  double *xlist, *Mlist, *dlist, *ylist;
    … …  
    309309  ALLOCATE (dlist, double, N);
    310310
     311  Narea = 0;
    311312  N = 0;
    312313  for (i = 0; i < Ngrid; i++) {
    … …  
    316317      c = clist[i][j];
    317318     
    318       if (catalog[c].measure[m].flags & MEAS_BAD) continue;
     319      if (catalog[c].measure[m].flags & MEAS_BAD) {
     320        Narea ++;
     321        continue;
     322      }
    319323      if ((Mcal  = getMcal  (m, c)) == NO_MAG) continue;
    320324      if ((Mmos  = getMmos  (m, c)) == NO_MAG) continue;
    … …  
    331335    }
    332336  }
     337
     338  fprintf (stderr, "skipped %d meas for area\n", Narea);
    333339
    334340  plot_defaults (&graphdata);

  • trunk/Ohana/src/relphot/src/ImageOps.c

    r5014 r5143  
    313313
    314314  /**** dMcal histgram ****/
    315   for (i = 0; i < NBIN; i++) xlist[i] = 0.0025*i;
     315  for (i = 0; i < NBIN; i++) xlist[i] = 0.00025*i;
    316316  bzero (Mlist, NBIN*sizeof(double));
    317317  for (i = 0; i < Nimage; i++) {
    318     bin = image[i].dMcal_PS / 0.0025;
     318    bin = image[i].dMcal_PS / 0.00025;
    319319    bin = MAX (0, MIN (NBIN - 1, bin));
    320320    Mlist[bin] += 1.0;

  • trunk/Ohana/src/relphot/src/MosaicOps.c

    r5014 r5143  
    630630  plot_list (&graphdata, xlist, Mlist, Nmosaic, "airmass vs Mcal", "airmass.png");
    631631  plot_defaults (&graphdata);
    632   plot_list (&graphdata, Mlist, dlist, Nmosaic, "Mcal vs dMcal", NULL);
     632  graphdata.size = 1.5;
     633  graphdata.ptype = 7;
     634  plot_list (&graphdata, Mlist, dlist, Nmosaic, "Mcal vs dMcal", "MdM.png");
    633635
    634636# define NBIN 200
    … …  
    637639
    638640  /**** dMcal histgram ****/
    639   for (i = 0; i < NBIN; i++) xlist[i] = 0.0025*i;
     641  for (i = 0; i < NBIN; i++) xlist[i] = 0.00005*i;
    640642  bzero (Mlist, NBIN*sizeof(double));
    641643  for (i = 0; i < Nmosaic; i++) {
    642     bin = mosaic[i].dMcal / 0.0025;
     644    bin = mosaic[i].dMcal / 0.00005;
    643645    bin = MAX (0, MIN (NBIN - 1, bin));
    644646    Mlist[bin] += 1.0;

  • trunk/Ohana/src/relphot/src/StarOps.c

    r4865 r5143  
    403403StatType statsStarS (Catalog *catalog, int Ncatalog) {
    404404
    405   int i, j, n, Ntot, dM;
     405  int i, j, n, Ntot;
    406406  double *list, *dlist;
     407  float dM;
    407408  StatType stats;
    408409
    … …  
    446447  ALLOCATE (Mlist, double, NBIN);
    447448
    448   for (i = 0; i < NBIN; i++) xlist[i] = 0.0025*i;
     449  for (i = 0; i < NBIN; i++) xlist[i] = 0.00025*i;
    449450  bzero (Mlist, NBIN*sizeof(double));
    450451  for (i = 0; i < Ncatalog; i++) {
    … …  
    452453      if (catalog[i].average[j].code & STAR_BAD) continue; 
    453454      dMrel = PhotPrimary ? catalog[i].average[j].dM_PS : catalog[i].secfilt[PhotNsec*j+PhotSec].dM_PS;
    454       bin = dMrel / 0.0025;
     455      bin = dMrel / 0.00025;
    455456      bin = MAX (0, MIN (NBIN-1, bin));
    456457      Mlist[bin] += 1.0;

  • trunk/Ohana/src/relphot/src/args.c

    r4865 r5143  
    110110  }
    111111
    112   MIN_ERROR = 1;
     112  MIN_ERROR = 0.001;
    113113  if ((N = get_argument (argc, argv, "-minerror"))) {
    114114    remove_argument (N, &argc, argv);

  • trunk/Ohana/src/relphot/src/global_stats.c

    r4796 r5143  
    2121  msX = statsMosaicX (catalog);
    2222
    23   fprintf (stderr, "STATS           median    mean   sigma     min      max  Nmeas\n");
    24   fprintf (stderr, "meas / image:   %6.0f  %6.1f  %6.1f  %6.0f  %7.0f  %6d\n",   imN.median, imN.mean, imN.sigma, imN.min, imN.max, imN.Nmeas);
    25   fprintf (stderr, "meas / mosaic:  %6.0f  %6.1f  %6.1f  %6.0f  %7.0f  %6d\n",   msN.median, msN.mean, msN.sigma, msN.min, msN.max, msN.Nmeas);
    26   fprintf (stderr, "meas / star:    %6.0f  %6.1f  %6.1f  %6.0f  %7.0f  %6d\n",   stN.median, stN.mean, stN.sigma, stN.min, stN.max, stN.Nmeas);
     23  fprintf (stderr, "STATS            median     mean    sigma      min      max   Nmeas\n");
     24  fprintf (stderr, "meas / image:   %7.0f  %7.0f  %7.0f  %7.0f  %7.0f  %6d\n",   imN.median, imN.mean, imN.sigma, imN.min, imN.max, imN.Nmeas);
     25  fprintf (stderr, "meas / mosaic:  %7.0f  %7.0f  %7.0f  %7.0f  %7.0f  %6d\n",   msN.median, msN.mean, msN.sigma, msN.min, msN.max, msN.Nmeas);
     26  fprintf (stderr, "meas / star:    %7.0f  %7.1f  %7.1f  %7.0f  %7.0f  %6d\n",   stN.median, stN.mean, stN.sigma, stN.min, stN.max, stN.Nmeas);
    2727
    28   fprintf (stderr, "chisq image:    %6.1f  %6.1f  %6.1f  %6.1f  %7.1f  %6d\n",   imX.median, imX.mean, imX.sigma, imX.min, imX.max, imX.Nmeas);
    29   fprintf (stderr, "chisq mosaic:   %6.1f  %6.1f  %6.1f  %6.1f  %7.1f  %6d\n",   msX.median, msX.mean, msX.sigma, msX.min, msX.max, msX.Nmeas);
    30   fprintf (stderr, "chisq star:     %6.1f  %6.1f  %6.1f  %6.1f  %7.1f  %6d\n",   stX.median, stX.mean, stX.sigma, stX.min, stX.max, stX.Nmeas);
     28  fprintf (stderr, "chisq image:    %7.1f  %7.1f  %7.1f  %7.1f  %7.1f  %6d\n",   imX.median, imX.mean, imX.sigma, imX.min, imX.max, imX.Nmeas);
     29  fprintf (stderr, "chisq mosaic:   %7.1f  %7.1f  %7.1f  %7.1f  %7.1f  %6d\n",   msX.median, msX.mean, msX.sigma, msX.min, msX.max, msX.Nmeas);
     30  fprintf (stderr, "chisq star:     %7.1f  %7.1f  %7.1f  %7.1f  %7.1f  %6d\n",   stX.median, stX.mean, stX.sigma, stX.min, stX.max, stX.Nmeas);
    3131
    32   fprintf (stderr, "Mcal image:     %6.3f  %6.3f  %6.3f  %6.3f  %7.3f  %6d\n",   imM.median, imM.mean, imM.sigma, imM.min, imM.max, imM.Nmeas);
    33   fprintf (stderr, "dMcal image:    %6.3f  %6.3f  %6.3f  %6.3f  %7.3f  %6d\n",   imD.median, imD.mean, imD.sigma, imD.min, imD.max, imD.Nmeas);
     32  fprintf (stderr, "Mcal image:     %7.4f  %7.4f  %7.4f  %7.4f  %7.4f  %6d\n",   imM.median, imM.mean, imM.sigma, imM.min, imM.max, imM.Nmeas);
     33  fprintf (stderr, "dMcal image:    %7.4f  %7.4f  %7.4f  %7.4f  %7.4f  %6d\n",   imD.median, imD.mean, imD.sigma, imD.min, imD.max, imD.Nmeas);
    3434
    35   fprintf (stderr, "Mcal mosaic:    %6.3f  %6.3f  %6.3f  %6.3f  %7.3f  %6d\n",   msM.median, msM.mean, msM.sigma, msM.min, msM.max, msM.Nmeas);
    36   fprintf (stderr, "dMcal mosaic:   %6.3f  %6.3f  %6.3f  %6.3f  %7.3f  %6d\n",   msD.median, msD.mean, msD.sigma, msD.min, msD.max, msD.Nmeas);
     35  fprintf (stderr, "Mcal mosaic:    %7.4f  %7.4f  %7.4f  %7.4f  %7.4f  %6d\n",   msM.median, msM.mean, msM.sigma, msM.min, msM.max, msM.Nmeas);
     36  fprintf (stderr, "dMcal mosaic:   %7.4f  %7.4f  %7.4f  %7.4f  %7.4f  %6d\n",   msD.median, msD.mean, msD.sigma, msD.min, msD.max, msD.Nmeas);
    3737
    38   fprintf (stderr, "dMrel star:     %6.3f  %6.3f  %6.3f  %6.3f  %7.3f  %6d\n\n", stS.median, stS.mean, stS.sigma, stS.min, stS.max, stS.Nmeas);
     38  fprintf (stderr, "dMrel star:     %7.4f  %7.4f  %7.4f  %7.4f  %7.4f  %6d\n\n", stS.median, stS.mean, stS.sigma, stS.min, stS.max, stS.Nmeas);
    3939
    4040  initstats (STATMODE);

  • trunk/Ohana/src/relphot/src/plot_scatter.c

    r4865 r5143  
    5151  graphdata.ymin = PlotdMmin;
    5252  graphdata.ymax = PlotdMmax;
    53   plot_list (&graphdata, ilist, ylist, N, "imag vs dmag", NULL);
     53  plot_list (&graphdata, ilist, ylist, N, "imag vs dmag", "iMag.png");
    5454  free (xlist);
    5555  free (ylist);
Note: See TracChangeset for help on using the changeset viewer.