if not keyword_set(lc) then begin
    lcPDS = pds_lcparam_getcatalog(/dump)
    val = where( strCmp(lcPDS.p.qual,'2',1) or strCmp(lcPDS.p.qual,'3',1) )
    lc = lcPDS[val]
    ok =where( lc.p.val gt 0)
    lc=lc[ok]
  endif

  help, lc.p,/str

  p=lc.p.val
  
  frac = 1/p

  cgHistoPlot, frac, minInput=0.0001, maxInput=1, xRange=[0.001,1], loc=x, histData=y

 

  cgPlot, x, 1.*y/max(y), /xlog, xRange=[0.0001, 1], psym=10, yRange=[0,1.1]

  cum = total(y,/Cum)
  tot = 1.*max(cum)
  cgPlot, /OverPlot, x, cum/tot, color='Cornflower Blue'




  root = '/data/euclid/'
  dir = root+'lc/pds/'
  readcol, dir+'statpds.csv', delimiter=',', /Silent, $
           LOW,HIGH,NEA,MC,MBA,Trojan,Centaur,KBO, comet
  x=(low+high)/2.

;  forprint, high-low

  readcol, root+'pop', delimiter=',', $
           format='(A,I1,F4.1,F4.1,F4.1,F4.1,F3.1,F3,A,A)', $
           popName, popMain, popHmin, popHmax, xMin, xMax, yMin, yMax, popC, popS, /Silent
  popName = strTrim(popName,2)
  nbPop=n_elements(popName)

  wdelete

  cgPS_open, xSize=24, ySize=13, Filename=dir+'euclid-periodstat.eps', $
             /Metric, /Decomposed, /Landscape, /Encapsulated, Language_Level=2, /Quiet


  cgPlot, x, NEA, /xLog, /NoData, $
          xRange=[0.01,1], $
          yRange=[0,1.05], $
          xTitle='Fraction of rotation period', $
          yTitle='Cumulative distribution', $
          position=[0.12,0.15,0.96,0.99]


  xKey = [0.013,0.02]
  yKey = 0.9
  yShi = -0.075


  kD=0
  ind = [7,6,5,4,2,1,0]
;  for kPop=0, nbPop-1 do begin
  for kI=0,6 do begin

    kPop = ind[kI]
    case popName[kPop] of
      'NEA': y=NEA
      'MC':  y=MC
      'MB':  y=MBA
      'Trojan': y=Trojan
      'Centaur': y=Centaur
      'KBO': y=KBO
      'Comet': y=comet
      else: goto, j2Next
    endcase

    cgPlot, /OverPlot, x, y, color=popC[kPop]
    cgPlot, /OverPlot, x, y, color=popC[kPop], psym=popS[kPop]

    cgPlot, /OverPlot, xKey, [yKey,yKey]+kD*yShi, color=popC[kPop]
    cgPlot, /OverPlot, mean(xKey), yKey+kD*yShi, color=popC[kPop], psym=popS[kPop]
    cgText, 0.022, yKey+kD*yShi-0.015, popName[kPop], color=popC[kPop]
    

    tttt = min( abs(x-0.20),pp )
    print, popName[kPop], 1-y[pp]


    kD++

    j2Next:
  endfor


  cgPS_close, /png, Delete_PS=0




end