calTree = getCalTree(obs=getObservation(1342192699, useHsa=True, instrument="PACS")) # Load the reference fits file and save the aligned version referenceFitsFile = dataDir + "smc-bridge-part1-mergedScans-" + cameraName + ".fits" referenceFrames = FitsArchive().load(referenceFitsFile) outputFitsFile = dataDir + "smc-bridge-part1-mergedScans-" + cameraName + "-aligned.fits" FitsArchive().save(outputFitsFile, referenceFrames) # Load the frames that should be aligned fitsFile = dataDir + "smc-bridge-part2-mergedScans-" + cameraName + ".fits" frames = FitsArchive().load(fitsFile) # Activate the non overlap masks and align the frames referenceFrames.setActive("Scanamorphos_NonOverlapMask", True) frames.setActive("Scanamorphos_NonOverlapMask", True) referenceFrames, frames = scanamorphosAlignFrames(referenceFrames, frames, calTree=calTree, createMask=False) # Deactivate the masks and save the aligned frames referenceFrames.setActive("Scanamorphos_NonOverlapMask", False) frames.setActive("Scanamorphos_NonOverlapMask", False) outputFitsFile = dataDir + "smc-bridge-part2-mergedScans-" + cameraName + "-aligned.fits" FitsArchive().save(outputFitsFile, frames)