obsids = [# virgo1
          1342212284, 1342212285,
          1342212301, 1342212302,
          1342211397, 1342211398,
          1342198147, 1342198148,
          1342210964, 1342210983,
          # virgo2
          1342199776, 1342199777,
          1342211415, 1342211416,
          1342201124, 1342201125,
          1342187265, 1342187266,
          # virgo3
          1342212357, 1342212358,
          1342212378, 1342212379,
          1342201139, 1342201140,
          1342198186, 1342198187,
          # virgo4
          1342212433, 1342212434,
          1342224141, 1342224032,
          1342199758, 1342199759,
          1342212407, 1342212408]

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

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

mergedFrames = None
for i in range(len(obsids)/2):
   fitsFile = dataDir + str(obsids[2*i]) +  "-" + str(obsids[2*i + 1])  +  "-" + cameraName + "-aligned.fits"
   f = FitsArchive().load(fitsFile)
   f = PhotHelper.selectMapRegion(f, maskedMap) 
   if f != None:
     if mergedFrames == None:
       mergedFrames = f
     else:
       mergedFrames.join(f)

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

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

Display(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 + "virgo-tile" + str(tile) + "-" + cameraName + ".fits"
FitsArchive().save(outputFitsFile, finalMap)


"""
  Merge all the aligned frames from the blue camera
"""
cameraName = "green"

mergedFrames = None
for i in range(len(obsids)/2):
   fitsFile = dataDir + str(obsids[2*i]) +  "-" + str(obsids[2*i + 1])  +  "-" + cameraName + "-aligned.fits"
   f = FitsArchive().load(fitsFile)
   f = PhotHelper.selectMapRegion(f, maskedMap) 
   if f != None:
     if mergedFrames == None:
       mergedFrames = f
     else:
       mergedFrames.join(f)

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

outputFitsFile = dataDir + "virgo-tile" + str(tile) + "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 + "virgo-tile" + str(tile) + "-" + cameraName + ".fits"
FitsArchive().save(outputFitsFile, finalMap)