; IDL> .r plot_usmag_xlumin_all
;
in_dir1='/disk/yla/takeda/REU_2026/qsbp_lumin/'
in_dir2='/disk/yla/takeda/REU_2026/qsbp_usmag/'
;-----------------------------
files1=file_search(in_dir1+'lumin_amap_*.genx',co=nn1)
files2=file_search(in_dir2+'usmag_hmi_*.genx',co=nn2)
;--- sanity check ---
fname1=file_break(files1)
fname2=file_break(files2)
for i=0,nn1-1 do print,fname1(i)+' - '+fname2(i)
;
for i=0,nn1-1 do begin
  restgen,file=files1(i),lumin_xrt_i
  restgen,file=files2(i),usmag_hmi_i
;
;  help,lumin_xrt_i,/str
;  help,usmag_hmi_i,/str
;
  if i eq 0 then begin
    lumin_qs=lumin_xrt_i.lumin_qs
    lumin_bps=lumin_xrt_i.lumin_bps
    usmag_qs=alog10(usmag_hmi_i.usmag_qs)
    usmag_bps=alog10(usmag_hmi_i.usmag_bps)
    qs_day=i
    bp_day=bytarr(n_elements(lumin_xrt_i.lumin_bps))+i
    yyyymmdd_arr=strmid(fname1(i),11,8)
  endif else begin
    lumin_qs=[lumin_qs, lumin_xrt_i.lumin_qs]
    lumin_bps=[lumin_bps, lumin_xrt_i.lumin_bps]
    usmag_qs=[usmag_qs, alog10(usmag_hmi_i.usmag_qs)]
    usmag_bps=[usmag_bps, alog10(usmag_hmi_i.usmag_bps)]
    qs_day=[qs_day,i]
    bp_day=[bp_day,bytarr(n_elements(lumin_xrt_i.lumin_bps))+i]
    yyyymmdd_arr=[yyyymmdd_arr,strmid(fname1(i),11,8)]
  endelse
endfor
;
help,lumin_qs,lumin_bps
help,usmag_qs,usmag_bps
;
wdef,0,512,512
; clr=[255,60,90,10,170,215,245]
;
set_plot,'z'
device,set_resolution=[512,512]
loadct,38
tvlct,rr,gg,bb,/get
;
erase
plot,usmag_bps,lumin_bps,/nodata,xrange=[18.5,22.5], xsty=1, $
     yran=[20,28],xtitle='HMI unsigned mag flux [log_Mx]', $
     ytitle='XRT luminosity [log_erg/s]', $
     title='Lx-B for BPs(+) and QS(x) with B_cutoff 10G'
oplot,usmag_bps,lumin_bps,psym=1
oplot,usmag_qs,lumin_qs,psym=7
;
;loadct,0
snap1=tvrd()
device,/close & set_plot,'x'
;
tvscl,snap1
;
;write_png,'./plot_usmag_xlumin_all.png',snap1,rr,gg,bb
;
FIN:
!p.charsize=1.0
end