"""
  Get the calibration tree
"""
calTree = getCalTree(obs=getObservation(obsids[0], useHsa=True, instrument="PACS"))


"""
  Merge all the aligned frames from the red camera
"""
camera = "red"
cameraName = "red"

fitsFile = dataDir + str(obsids[0]) +  "-" + str(obsids[1])  +  "-" + cameraName + "-aligned.fits"
mergedFrames = FitsArchive().load(fitsFile)

for i in range(len(obsids)/2 - 1):
   fitsFile = dataDir + str(obsids[2*(i+1)]) +  "-" + str(obsids[2*(i+1) + 1])  +  "-" + cameraName + "-aligned.fits"
   mergedFrames.join(FitsArchive().load(fitsFile))

"""
  Create the wcs
"""
wcs = wcs4map(mergedFrames, outputPixelsize=outputPixelSize, crota=crota, calTree=calTree)

"""
  Create the red map including all data
"""
finalMap, mi = photProject(mergedFrames, outputPixelsize=outputPixelSize, pixfrac=pixfrac, calTree=calTree, wcs=wcs)

outputFitsFile = dataDir + "field" + field + "Ext" + "-" + cameraName + ".fits"
FitsArchive().save(outputFitsFile, finalMap)

"""
  Create the red map without the borders
"""
mergedFrames.setActive("Scanamorphos_NonOverlapMask", True)

finalMap, mi = photProject(mergedFrames, outputPixelsize=outputPixelSize, pixfrac=pixfrac, calTree=calTree, wcs=wcs)

outputFitsFile = dataDir + "field" + field + "-" + cameraName + ".fits"
FitsArchive().save(outputFitsFile, finalMap)


"""
  Merge all the aligned frames from the blue camera
"""
camera = "blue"
cameraName = "blue"
try:
  fitsFile = dataDir + str(obsids[0]) +  "-" + str(obsids[1])  +  "-" + cameraName + "-aligned.fits"
  mergedFrames = FitsArchive().load(fitsFile)
except:
  if cameraName == "blue": 
    cameraName = "green"
  fitsFile = dataDir + str(obsids[0]) +  "-" + str(obsids[1])  +  "-" + cameraName + "-aligned.fits"
  mergedFrames = FitsArchive().load(fitsFile)

for i in range(len(obsids)/2 - 1):
   fitsFile = dataDir + str(obsids[2*(i+1)]) +  "-" + str(obsids[2*(i+1) + 1])  +  "-" + cameraName + "-aligned.fits"
   mergedFrames.join(FitsArchive().load(fitsFile))

"""
  Create the blue map including all data
"""
finalMap, mi = photProject(mergedFrames, outputPixelsize=outputPixelSize, pixfrac=pixfrac, calTree=calTree, wcs=wcs)

outputFitsFile = dataDir + "field" + field + "Ext" + "-" + cameraName + ".fits"
FitsArchive().save(outputFitsFile, finalMap)

"""
  Create the blue map without the borders
"""
mergedFrames.setActive("Scanamorphos_NonOverlapMask", True)

finalMap, mi = photProject(mergedFrames, outputPixelsize=outputPixelSize, pixfrac=pixfrac, calTree=calTree, wcs=wcs)

outputFitsFile = dataDir + "field" + field + "-" + cameraName + ".fits"
FitsArchive().save(outputFitsFile, finalMap)