obsids  = [1342195728,
          1342195713,
          1342195712,
          1342195708,
          1342195707,
          1342195684,
          1342195683,
          1342195669,
          1342195668,
          #
          1342187189,
          1342187188,
          # cscans
          1342202086,
          1342202087,
          1342202202,
          1342202203,
          1342202216,
          1342202217,
          1342202224, 
          1342202225,
          1342202243,
          1342202244]

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

for camera in ["red", "blue"]:
   for obsid in obsids:
      # Load the frames from the archive
      obs = getObservation(obsid, useHsa=True, instrument="PACS")
      level1 = PacsContext(obs.level1)
      scans = level1.averaged.getCamera(camera).product.selectAll()
      # Create some of the masks
      scans = scanamorphosRemoveTurnarounds(scans, limit=50.0)
      scans = scanamorphosMaskLongTermGlitches(scans, stepAfter=20, galactic=True, calTree=calTree)
      # Save the frames into disk
      if(camera == "blue"):
         if(scans.meta["blue"].value == "blue2"):
            cameraName = "green"
         else: 
            cameraName = "blue" 
      else: 
         cameraName = "red" 
      outputFitsFile = dataDir + "lmc" + str(obsid) +  "-" + cameraName + ".fits"
      FitsArchive().save(outputFitsFile, scans)