;plot_ellipticities
;
;PURPOSE:
;    This function will take an a 4x12xN array of filtered ellipticities and
;    elliptiticy angles, histogram them, and plot them.
;
;    
;INPUTS:
;    EllipticityArr:
;      EllipticityArr(0,k,*) = Ellipticities of High Pass Filtered at cutoff k
;      EllipticityArr(1,k,*) = Angles of High Pass Filtered at cutoff k
;      EllipticityArr(2,k,*) = Ellipticities of Low Pass Filtered at cutoff k
;      EllipticityArr(3,k,*) = Angles of Low Pass Filtered at cutoff k
;   
;    PlotStyle:
;      0 = Plot to X Windwow
;      1 = Plot a PNG file
;      
function hist_ellipticities,EllipticityArr,PlotStyle

;******************************************************UNFILTERED VALUES
E=EllipticityArr(0,0,*)
Theta=EllipticityArr(1,0,*)

;******************************************************THETA HISTOGRAM
ThetaMin=-!pi/2                         
NoBins=20.
ThetaBinWidth=!pi/NoBins
ThetaMax=ThetaMin+(NoBins-1)*ThetaBinWidth

ThetaHist=histogram(Theta,Nbins=NoBins,locations=LocTheta,$
                 reverse_indices=ThetaInd,$
                 max=ThetaMax,min=ThetaMin)
ThetaBinSize=(ThetaMax-ThetaMin)/(NoBins-1)

;**********************************************************E HISTOGRAM
Emin=0                                     ;ELLIPTICITY
NoEBins=20.
EBinWidth=(1-Emin)/NoEBins
Emax=Emin+(NoEBins)*EBinWidth
Ehist=histogram(E,Nbins=NoEBins,locations=LocE,$
                 reverse_indices=EInd,$
                 max=Emax,min=Emin)
Ebinsize=(Emax-Emin)/(NoEBins-1)

EValsArr=fltarr(4,N_elements(EHist))
EValsArr(0,*)=LocE
EValsArr(1,*)=EHist
EValsArr(2,*)=LocTheta
EValsArr(3,*)=ThetaHist

return,EValsArr

end