;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