;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 plot_ellipticities,EllipticityArr,PlotStyle

if PlotStyle eq 1 then begin
  set_plot,'z'
  erase
  device, set_font='Courier'
  device,set_resolution=[800,600]
  !p.charsize=1
  !p.charthick=1.2
  !x.thick=2
  !y.thick=2
  !p.thick=4
  !p.noerase=0
  loadct,39

  white='FFFFFF'x
  black='000000'x
  blue='FF0000'x
  red='0000FF'x
endif else begin
  set_plot,'x'
endelse

!p.multi=[0,1,1]

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

plot,LocTheta+ThetaBinSize/2,ThetaHist,psym=10,xtitle='Theta',$
     ytitle='Number of Events',xrange=[-1.6,1.6],$
     background=white,color=black


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

plot,LocE,Ehist,psym=10,xtitle='e',$
     ytitle='Number of Events',$
     background=white,color=black

stop
return,1

end