 		photonVals[y][x] = 0;
 		lambdaVals[y][x] = 0;
 	      }*/
-	      if (x > 393 || y > 411 || (/*(x + y > 700) ||*/ (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 (x > 393 || y > 411 || (/*(x + y > 700) ||*/ (y < 235 || (x < 300 && y < 314) || x < 255)  && !(y < 183 && y > 39 && x > 290 && x < 393) &&
+	      !(x < 170 && x > 60 && y < 98 && y > 40) &&
+	      !(x < 290 && x > 254 && y < 120 && y > 39) &&
+	      !(x < 88 && x > 19 && y < 130 && y > 97)))
 //	      if ((y < 233 || (x < 350 && y < 370) || x < 255))
 	      {
 		photonVals[y][x] = 0;

intens[idx] = 0;*/
intensityfactor = 1e-13;
}
  intensityfactor *= pow(1.2, /*by*/ - 0 * 0.5);
  likelihood likelihooder(myspherer, intensityfactor);

  return likelihooder;

__host__ __device__ likelihood(idealsphere& spherer, float factor) : spherer(spherer), factor(factor) {}

  __host__ __device__ float getIntensity(const int& data) {
  	float lbase = spherer.dLambdas[data] * spherer.lfactor;
	lbase += spherer.ebase;
	if (lbase < 2e-3) lbase = 2e-3;
  	float intensity = spherer(data) * factor + lbase;

	return intensity;
  }
