 
 
     float factor = (/*0.01f + 0.99f **/ expf(- gr / (2 * sigmag * sigmag)));
-    //factor *= sinc(qx / sigma * 0.5) * sinc(qy / sigma * 0.5);
+    factor *= sinc(qx / sigma * 0.5) * sinc(qy / sigma * 0.5);
 
     //float factor = 1.0f;
 

intens[idx] = 0;*/
intensityfactor = 0;
}
  intensityfactor *= pow(1.1f, -1.f + by * 2.f);
  likelihood likelihooder(myspherer, intensityfactor);

  return likelihooder;

	}
      
//      if (psum - lsum < 2500) continue;
      dim3 grid(48, 1, imgcount);
      dim3 block(22, 22, 1);

      int base = (grid.y * block.y * grid.x * block.x * block.z);
