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

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

      spherer.rfactor = 0.005;
      spherer.roffset = -10;
      spherer.baseoffsetx = NX / 2 - 10 - 21 - 0.5; // good val -10
      spherer.baseoffsety = NY / 2 + 10 - 11 - 0.5; // good val +10
      dPhotons.assign(photonVals.data(), photonVals.data() + NY * NX);
      dLambdas.assign(lambdaVals.data(), lambdaVals.data() + NY * NX);

	{
		maxval = hIntensity[k];
		maxint = hIntensity2[k];
		maxidx = k;
	}
	if (hIntensity[k] < minval) minval = hIntensity[k];
}
