from herschel.pacs.spg.phot import PhotHelper

for field in fields:
  # Load the maps
  redMap = FitsArchive().load(dataDir + "field" + str(field) + "-red.fits")
  redMapExt = FitsArchive().load(dataDir + "field" + str(field) + "Ext-red.fits")
  blueCamera = "blue"
  try:
    blueMap = FitsArchive().load(dataDir + "field" + str(field) + "-" + blueCamera + ".fits")
    blueMapExt = FitsArchive().load(dataDir + "field" + str(field) + "Ext-" + blueCamera + ".fits")
  except:
    blueCamera = "green"
    blueMap = FitsArchive().load(dataDir + "field" + str(field) + "-" + blueCamera + ".fits")
    blueMapExt = FitsArchive().load(dataDir + "field" + str(field) + "Ext-" + blueCamera + ".fits")
  
  # Combine the maps
  redMapCombined =  PhotHelper.combineMaps(redMap, redMapExt)
  blueMapCombined =  PhotHelper.combineMaps(blueMap, blueMapExt)
  
  # Remove the maps offsets
  PhotHelper.removeOffset(redMapCombined, redOffset)
  PhotHelper.removeOffset(blueMapCombined, blueOffset)
  
  # Save the maps
  FitsArchive().save(dataDir + "field" + str(field) + "-red-corr.fits", redMapCombined)
  FitsArchive().save(dataDir + "field" + str(field) + "-" + blueCamera + "-corr.fits", blueMapCombined)