pro feplot,fefile,dcfile,prn
@fe55setup

  if n_elements(prn) eq 0 then prn=0

  print,'Reading in dark frame...'
  fzrd,dcfile,d,hd,eng
  print,'Reading in Fe55 exposure...'
  fzrd,fefile,f,hd,eng
  tit=strupcase(fefile)+'   '+string(hd.date)

  print,'Plotting DC and FE histograms...'
  dcplots,f,d,tit,fe,hfx,hfy,prn
  if prn eq 0 then begin
    print,'--> Please press any key to continue'
    waite=get_kbrd(1)
    erase
  endif

  lowlim=replicate(lower,1076,2)
  uplim=replicate(upper,1076,2)

  print,'Fitting Fe55 line...'
  fe55fit,fe,lowlim,uplim,tit,prn
  print,'Done !...'

end