 intens[idx] = 0;*/
 intensityfactor = 1e-13;
 }
-  intensityfactor *= pow(1.0075, myspherer.tid - 64 * 0.5);
+//  intensityfactor *= pow(1.0075, myspherer.tid - 64 * 0.5);
   likelihood likelihooder(myspherer, intensityfactor);
 
   return likelihooder;

//      logLs.write(logLsVals.data(), H5::PredType::NATIVE_FLOAT, memSpace, expectedLambdaSpace);
      //expectedLambdaFile.flush(H5F_SCOPE_LOCAL);
      //hIntensity2.assign(dIntensity2.begin(), dIntensity2.end());
	int maxR = maxidx / grid.y / block.y / grid.x / block.x;
	int maxX = maxidx % block.x;
	int maxY = (maxidx / (grid.x * block.x)) % block.y;
	int maxI = (maxidx / (block.x * grid.x * block.y)) % grid.y;
	int maxI2 = (maxidx / block.x) % grid.x;

      printf("%d %d %lf %g %g %g %d %d %d %d %g %d %d\n", img, psum, lsum, minval, maxval, hIntensity[0], maxR, maxX, maxY, cudaGetLastError(), maxint, maxI, maxI2);
      fflush(stdout);
    }
}

//    if (rawx < 233 - extrax * 3) rawx += /*29*/ -18 + extrax * 3;
//    if (rawy < 226 - extray * 2) rawy += /*3*/ -12 + extray * 2;*/

    int x = rawx - offsetx;
    int y = rawy - offsety;
    if (x < 0) x += FX;
    if (y < 0) y += FY;

    return dPattern[y * FX + x];
  }
};

