k2d = '/home/bcarry/work/mining/k2/' fC = 0 lC = 7 ccd = 21 - 2 pix = 2200 * 1024. totArea = ccd*pix fov=160. for kC=fC, lC do begin readcol, k2d+'K2Campaign'+strtrim(string(kC,format='(I)'),2)+'targets.csv', $ epic, ra, dec, mag, prog, delimiter=',', $ format='(L,F,F,F,A)', /Silent nbFOV = n_elements(dec) valid = where( ra ne 0 and dec ne 0, nbFOV) epic=epic[valid] ra=ra[valid] dec=dec[valid] mag=mag[valid] muRa = median(ra) muDec= median(dec) area = fov*nbFov frac = 100.*area/totArea ;---ecliptic ecL = fltarr(nbFOV) ecB = fltarr(nbFOV) for k=0, nbFov-1 do begin aa=framecoord_eq2ec( ra[k], dec[k] ) ecL[k]=aa[0] ecB[k]=aa[1] endfor ; stop ; cgwindow ; cgPlot, ra[valid], dec[valid], /iso, psym=3, /AddCmd ; cgwindow ; cgPlot, ecL, ecB, /iso, psym=3 binSize=30./3600. cgHistoPlot, ecB, binSize=binSize, location=bArr, histData=filling miss = where( filling eq 0, nbMiss ) len = max(bArr) - min(bArr) fillFac = 100-100.*nbMiss*binSize/len meanEnc = mean( filling ) stdEnc = stddev( filling ) print, ' C'+strtrim(string(kC,format='(I)'),2)+' & '+string(nbFOV,format='(I5)')+' & '+$ string(frac,format='(F6.2)')+' & '+$ string(fillFac,format='(F6.2)')+' & '+$ string(meanEnc,format='(F4.1)')+' & '+$ string(stdEnc,format='(F4.1)')+' \\ ' ; binSize=300./3600. ; cgHistoPlot, ecB, binSize=binSize, location=bArr, histData=filling ;----histo ; cgHistoPlot, mag, binSize=0.5, /oProbability, histData=histo, location=kArr ; cumul=total(histo,/Cumulative) ; ; cutK=0.8 ; low = where( cumul/max(cumul) le 0.8 ) ; magK = max( kArr[low] ) ; print, magK endfor end