 		photonVals[y][x] = 0;
 		lambdaVals[y][x] = 0;
 	      }*/
-	      if (y < 233 || (x < 300 && y < 314) || x < 255)
-//	      if (y < 233 || (x < 350 && y < 370) || x < 255)
+	      if (x > 393 || y > 411 || ((y < 235 || (x < 300 && y < 314) || x < 255)  && !(y < 183 && y > 39 && x > 312 && x < 393) &&
+	      !(x < 170 && x > 60 && y < 98 && y > 40)))
+//	      if ((y < 233 || (x < 350 && y < 370) || x < 255))
 	      {
 		photonVals[y][x] = 0;
 		lambdaVals[y][x] = 0;

__device__ likelihood getObjects(idealsphere& myspherer, unsigned int tx, unsigned int ty, unsigned int zval, unsigned int bx, unsigned int by, float lsum, float psum)
{
  myspherer.r = exp(myspherer.roffset + (zval) * myspherer.rfactor);
  myspherer.offsetx = (tx) * 4 + myspherer.baseoffsetx;
  myspherer.offsety = (ty * 4 + myspherer.baseoffsety);
  myspherer.lfactor = 0.25 / sqrt(lsum) * (((int) bx)) + 1.0;
//    myspherer.lfactor = 1.0;


  
  __host__ __device__ float operator () (const int& data) {
    if (!dLambdas[data]) return 0;
    int rawx = data % NX;
    int rawy = data / NX;

    if (rawx < 216) rawx += 31;
    if (rawy < 198) rawy += 31;

    float x = rawx - offsetx;
    float y = rawy - offsety;
    
    float q = sqrt(x * x + y * y);
    float val = 3 * (sinpif(2 * q * r)  - 2 * ((float) CUDART_PI_F) * q * r * cospif(2 * q * r));
