obsids  = [1342205092]

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()
      # Use the old pointing correction   
      pp = obs.auxiliary.pointing
      orbitEphem = obs.auxiliary.orbitEphemeris
      horizons = obs.auxiliary.horizons
      sets = pp.getSets()
      setsArray = sets.keySet().toArray()
      for set in setsArray:
         pp[set]["filterQuat"].data =  pp[set]["simpleCorrFilterQuat"].data
      scans = photAddInstantPointing(scans, pp, calTree=calTree, orbitEphem=orbitEphem, horizonsProduct=horizons)
      # Create some of the masks
      scans = scanamorphosRemoveTurnarounds(scans, limit=50.0)
      scans = scanamorphosMaskLongTermGlitches(scans, stepAfter=20, galactic=True, calTree=calTree)
      scans = scanamorphosBaselinePreprocessing(scans)
      # Save the frames into disk
      if(camera == "blue"):
         if(scans.meta["blue"].value == "blue2"):
            cameraName = "green"
         else: 
            cameraName = "blue" 
      else: 
         cameraName = "red" 
      outputFitsFile = dataDir + "smc" + str(obsid) +  "-" + cameraName + ".fits"
      FitsArchive().save(outputFitsFile, scans)
      del pp, orbitEphem, horizons, sets, setsArray, set