;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