pro Darkcurrent_Diff, Date = Date, RootDir = RootDir, ReducedDir = ReducedDir, $
	              ObjectName = ObjectName, NReads = NReads, $
		      OverRide = OverRide, CheckFiles = CheckFiles
 Common FitsCube, FCB, unit

 If not(keyword_set(Date)) then Date='07Apr26/'
 If not(keyword_set(RootDir)) then RootDir='/nfs/slac/g/ki/ki09/lances/'
 If not(keyword_set(ReducedDir)) then ReducedDir = Rootdir+'Reduced/'+Date
 cd, RootDir+Date

  If not(keyword_set(ObjectName)) then ObjectName = 'Dark'
  If not(keyword_set(NReads)) then NReads = 15
  If not(keyword_set(ThisFilter)) then ThisFilter = 'blank'
  If not(keyword_set(OverRide)) then OverRide = 0
  If not(keyword_set(CheckFiles)) then CheckFiles = 0

  ;Get all of the Dark Files for this Number of Reads
  Dark_Files = file_search(ObjectName+'*'+strtrim(NReads,2)+'_Reads*.fits')
  Dummy = 0
  Valid_Indices = 0

  ;Get a few kewyords
  Fits_Read,  Dark_Files(0), NoData, FitsHeader_Dark, /header_only
  Naxis1     = long(sxpar(FitsHeader_Dark, 'NAXIS1'))
  Naxis2     = long(sxpar(FitsHeader_Dark, 'NAXIS2'))
  TotalReads = sxpar(FitsHeader_Dark ,'NAXIS3')
  ArrSize = long(Naxis1*Naxis2)

  ;Only use the images that are certain to have Blank as the filter
  for i = 0, n_elements(Dark_Files) - 1 do begin
     Fits_Read, Dark_Files(i), NoData, FitsHeader, /header_only
     Filter = sxpar(FitsHeader ,'FILTER')
     print, i, ' ' , Filter
     If stregex(Filter, 'Blank' , /boolean) and $
        not stregex(Dark_Files(i), '42_44', /boolean) then begin
          valid_indices = [valid_indices, i]
          If CheckFiles then begin
             Fits_Read, Dark_Files(i), Data, Header
             Tvscl, congrid(Data(*,*,0), 512,512)
             print, Dark_Files(i)
          endif
     EndIf
     
   endfor

  ;IPC in 29, 28 : (2324:2328, 2022:2026) 
  Im = Return_Cubes_Diff(Dark_Files(valid_indices(10)), $
			 Dark_Files(valid_indices(11)), $
			 ArrSize, Naxis1, Naxis2, 1, 1 )

  ;Remove 60 Hz noise if necessary
  Remove60Hz = 0
  if Remove60Hz then begin
     medvector=median(im(128:Naxis2-1,*),dimension=1)
     stop
     for row = 0, Naxis2-1 do begin
       im(*,row) = im(*, row) - medvector(row)
     endfor
  endif
 
  NoBins = 100
  CountMin = -100.
  CountMax = 100.

  DarkHist = histogram(Im,Nbins=NoBins,locations=LocCounts,$
                 reverse_indices=Count_Ind,$
                 max=CountMax,min=CountMin)

  DarkBinSize=(CountMax-CountMin)/(NoBins-1)

  plot, LocCounts+DarkBinSize/2, DarkHist, psym=10,$
        xtitle='Counts',$
        ytitle='Number of Pixels'
  Fits_Open, ReducedDir+'Diff.fits', FCB, /write 
  Fits_Write, FCB, Im 
  stop

end