calTree = getCalTree(obs=getObservation(1342205092, useHsa=True, instrument="PACS"))

for i in range(len(tiles)):
  if tiles[i] == referenceTiles[i]:
     referenceFitsFile = dataDir + "smc-mergedScans-tile" + str(tiles[i]) +  "-" + cameraName + ".fits"
     referenceFrames = FitsArchive().load(referenceFitsFile)
     outputFitsFile = dataDir + "smc-mergedScans-tile" + str(tiles[i]) +  "-" + cameraName + "-aligned.fits"
     FitsArchive().save(outputFitsFile, referenceFrames)
  else:
     """
      Load the frames and align them
     """
     referenceFitsFile = dataDir + "smc-mergedScans-tile" + str(referenceTiles[i]) +  "-" + cameraName + "-aligned.fits"
     fitsFile = dataDir + "smc-mergedScans-tile" + str(tiles[i]) +  "-" + cameraName + ".fits"
     referenceFrames = FitsArchive().load(referenceFitsFile)
     frames = FitsArchive().load(fitsFile)
     referenceFrames.setActive("Scanamorphos_NonOverlapMask", True)
     frames.setActive("Scanamorphos_NonOverlapMask", True)
     referenceFrames, frames = scanamorphosAlignFrames(referenceFrames, frames, calTree=calTree, createMask=False)
     
     """
      Write the fits files
     """
     frames.setActive("Scanamorphos_NonOverlapMask", False)
     outputFitsFile = dataDir + "smc-mergedScans-tile" + str(tiles[i]) +  "-" + cameraName + "-aligned.fits"
     FitsArchive().save(outputFitsFile, frames)
     del(referenceFrames, frames)