for i in range(len(tiles)):
  tile = tiles[i]
  if i == 0:
    fitsFile = dataDir + "virgo-tile" + str(tile) +  "-" + cameraName + ".fits"
    map = FitsArchive().load(fitsFile)
  else:
    fitsFile = dataDir + "virgo-tile" + str(tile) +  "-" + cameraName + ".fits"
    tileMap = FitsArchive().load(fitsFile)
    map = PhotHelper.addTile(map, tileMap, tileStartCol[tile], tileStartRow[tile], tileColSize[tile], tileRowSize[tile], nCols, nRows)

# Calculate the background offset and the signal limits
from herschel.pacs.spg.phot import PhotHelper
minimumRelativeCoverage = 0.1
boxSize = 40
backgroundOffset = PhotHelper.getBackgroundOffset(map, minimumRelativeCoverage, boxSize)
signalLim = PhotHelper.getSignalLimits(map, minimumRelativeCoverage)
print backgroundOffset
print signalLim

for i in range(len(tiles)):
  tile = tiles[i]
  if i == 0:
    fitsFile = dataDir + "virgo-tile" + str(tile) +  "Ext-" + cameraName + ".fits"
    mapExt = FitsArchive().load(fitsFile)
  else:
    fitsFile = dataDir + "virgo-tile" + str(tile) +  "Ext-" + cameraName + ".fits"
    tileMap = FitsArchive().load(fitsFile)
    mapExt = PhotHelper.addTile(mapExt, tileMap, tileStartCol[tile], tileStartRow[tile], tileColSize[tile], tileRowSize[tile], nCols, nRows)

map = PhotHelper.combineMaps(map, mapExt)

Display(map)

fitsFile = dataDir + "virgo-" + cameraName + ".fits"
FitsArchive().save(fitsFile, map)