KBO={mag:1.5, dMag:2.0, sep:0.43 , dSep:0.60 , f:5.5, dF:4. }
  MBA={mag:5.4, dMag:2.7, sep:0.30 , dSep:0.25 , f:2.5, dF:1. }
  sMB={mag:2.5, dMag:0.9, sep:0.01 , dSep:0.005, f:17., dF:3. }
  NEA={mag:1.8, dMag:2.0, sep:0.002, dSep:0.006, f:17., dF:3. }

  binO= 0.010
  binS= 0.010


 ;--- Orbit max Distance projected in one direction
  cghistoplot, abs(sin(randomu(1,1000000)*2*!PI)),loc=oD, hist=oPDF, binSize=binD
  nbO = n_elements(oD)
  wdelete


  help, oD, oPDF
  

  for k=0, 2 do begin
    case k of
      0: pop=KBO
      1: pop=MBA
      2: pop=sMB
      3: pop=NEA
    endcase

    nbS = 200
    sD = findgen(nbS) * binS

    sPDF = evalGauss( sD, [1., pop.sep, pop.dSep] )
    sPDF /= total(sPDF)

    help, sD, sPDF


    distArr = fltarr( nbO, nbS )
    pdfArr  = fltarr( nbO, nbS )

    for kO=0, nbO-1 do for kS=0, nbS-1 do begin
      distArr[kO,kS] = oD[kO] * sD[kS]
      pdfArr[kO,kS]  = oPDF[kO] * sPDF[kS]
    endfor


    
    cgPlot, sD, sPDF
    window, 2
    cgPlot, oD, oPDF

    window, 4
    cgPlot, distArr, pdfArr, psym=3


    sepVal = sD
    sepPDF = sD * .0
    for kS=0, nbS-1 do begin

      loc =where( distArr ge binD*kS and distArr lt binD*(kS+1) )
      if loc[0] ne -1 then sepPDF[kS] = total( pdfArr[loc] )
  
    endfor
    sepPDF /= total(sepPDF)


    window, 4
    cgPlot, sepVal, total(sepPDF,/Cumulative), psym=3

;    atv, distARr
stop



 endfor









end