;DSL coordinates
;
; energy spectrogram
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'], $
                  data_type=['psif'],/energy, $
                  phi=[-135,-45], theta=[-45,45], erange=[25000,500000],suff='_dusk'
;
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'], $
                  data_type=['psif'],/energy, $
                  phi=[45,135], theta=[-45,45], erange=[25000,500000],suff='_dawn'
;
; phi spectrogram
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'], $
                  data_type=['peir'],angle='phi', $
                  phi=[0,360], theta=[-90,90], erange=[1.5e4,2.5e4]
;
; theta spectrogram
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'], $
                  data_type=['peir'],angle='theta', $
                  phi=[0,360], theta=[-90,90], erange=[1.5e4,2.5e4]
;
tplot,'thb_fgs_gsm thb_psif_en_eflux_dusk thb_peir_an_eflux_*'
tlimit,['07-03-23/11:12','07-03-23/11:22']
;
; FAC coordinates
;
; Energy spectrogram
;
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'],$
                  data_type=['psif'], /energy, $
                  pitch=[0,45], suff='_para', $
                  erange=[25000,500000],regrid=[32,16]
;
; Gyro(velocity)phase spectrogram
;
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'],$
                  data_type=['psif'], angle='gyro', $
                  pitch=[45,135], other_dim='ygsm', suff='_perp', $
                  erange=[100000,150000],regrid=[32,16]
;
; Pitch angle spectrogram
;
thm_part_getspec, probe=['b'], trange=['07-03-23/11:10','07-03-23/11:30'],$
                  data_type=['peer'], angle='pa', $
                  erange=[15000,25000],regrid=[32,16]
;
tplot,'thb_fgs_gsm thb_psif_en_eflux_para  thb_psif_an_eflux_gyro_perp thb_peer_an_eflux_pa'
tlimit,['07-03-23/11:12','07-03-23/11:22']
;
; Sun contamination
;
trange=['08-02-22/04:40','08-02-22/05:00']
timespan,'08-02-22/04:00',1,/hours
sc='c'
;
; phi spectrogram
;
thm_load_fgm,probe=sc,data='fgl',coord='dsl gsm'
;
thm_part_getspec, probe=sc,trange=trange, $
                  theta=[0,45], phi=[0,360], start_angle=-180., $
                  erange=[25000,40000], data_type=['psif'], $
                  suffix='_ne',angle=phi
;
thm_part_getspec, probe=sc,trange=trange, $
                  theta=[-45,0], phi=[0,360], start_angle=-180., $
                  erange=[25000,40000], data_type=['psif'], $
                  suffix='_se',angle=phi
;
tplot,'thc_fgl_dsl thc_fgl_gsm thc_psif_an_eflux_phi_ne thc_psif_an_eflux_phi_se'
tlimit,['08-02-22/04:40','08-02-22/05:00']
;
; Sun contamination: remedy
;
edit3dbins,thm_sst_psif(probe=sc, gettime(/c)), bins2mask
; ON: Button1;    OFF: Button2;   QUIT: Button3
bins2mask = where(bins2mask eq 0)
print,bins2mask
thm_part_getspec, probe=sc,trange=trange, $
                  theta=[0,45], phi=[0,360], start_angle=-180., $
                  erange=[25000,40000], data_type=['psif'], $
                  suffix='_ne',angle=phi,bins2mask=bins2mask
;
thm_part_getspec, probe=sc,trange=trange, $
                  theta=[-45,0], phi=[0,360], start_angle=-180., $
                  erange=[25000,40000], data_type=['psif'], $
                  suffix='_se',angle=phi,bins2mask=bins2mask
;
tplot