;hist_ellipticities_image
;
;PURPOSE:
;    This function will take an a 4xN array of e1s, e2s, ellipticities and
;    elliptiticy angles, histogram them, and plot them.
;
;    
;INPUTS:
;    EllipticityArr:
;      EllipticityArr(0,*) = E1
;      EllipticityArr(1,*) = E2
;      EllipticityArr(2,*) = E
;      EllipticityArr(3,*) = Theta
;   
;    PlotStyle:
;      0 = Plot to X Windwow
;      1 = Plot a PNG file
;      
function hist_ellipticities_image, EllipticityArr, PlotStyle

;******************************************************UNFILTERED VALUES
E=EllipticityArr(2,*)
Theta=EllipticityArr(3,*)

;******************************************************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

stop
end