;+ ;FUNCTION: tp_4d(dat,ENERGY=en,ERANGE=er,EBINS=ebins,ANGLE=an,ARANGE=ar,BINS=bins,MASS=ms,m_int=mi,q=q,mincnt=mincnt) ;INPUT: ; dat: structure, 4d data structure filled by themis routines mvn_sta_c6.pro, mvn_sta_d0.pro, etc. ;KEYWORDS ; ENERGY: fltarr(2), optional, min,max energy range for integration ; ERANGE: fltarr(2), optional, min,max energy bin numbers for integration ; EBINS: bytarr(na), optional, energy bins array for integration ; 0,1=exclude,include, ; na = dat.nenergy ; ANGLE: fltarr(2,2), optional, angle range for integration ; theta min,max (0,0),(1,0) -90<theta<90 ; phi min,max (0,1),(1,1) 0<phi<360 ; ARANGE: fltarr(2), optional, min,max angle bin numbers for integration ; BINS: bytarr(nb), optional, angle bins array for integration ; 0,1=exclude,include, ; nb = dat.ntheta ; BINS: bytarr(na,nb), optional, energy/angle bins array for integration ; 0,1=exclude,include ;PURPOSE: ; Returns the perp temperature of a beam from APID C8 in units of eV ;NOTES: ; Function normally called by "get_4dt" to ; generate time series data for "tplot.pro". ; ;CREATED BY: ; J.McFadden 2014-03-12 ;LAST MODIFICATION: ;- function tp_4d,dat2,ENERGY=en,ERANGE=er,EBINS=ebins,ANGLE=an,ARANGE=ar,BINS=bins,MASS=ms,m_int=mi,q=q,mincnt=mincnt if dat2.valid eq 0 then begin print,'Invalid Data' return, 0 endif if dat2.apid ne 'C8' and dat2.apid ne 'c8' then begin print,'Invalid Data: Data must be Maven APID C8' return, 0 endif dat = conv_units(dat2,"counts") ; initially use counts nth = dat.ndef n_e = dat.nenergy n_m = dat.nmass mass_amu = dat.mass_arr data = dat.data energy = dat.energy if keyword_set(en) then begin ind = where(energy lt en[0] or energy gt en[1],count) if count ne 0 then data[ind]=1.e-20 endif if keyword_set(mincnt) then if total(data) lt mincnt then return,0 charge=dat.charge if keyword_set(q) then charge=q energy=(dat.energy+charge*dat.sc_pot/abs(charge))>0. ; energy/charge analyzer, require positive energy v = energy^.5 sth = sin(dat.theta/!radeg) v0 = total(v*sth*data,2)/(total(data,2)>1.e-20) tp = total((v*sth - v0#replicate(1.,nth))^2*data)/(total(data)>1.e-20) return, tp ; eV end