interpolation without quads

This commit is contained in:
2020-01-24 15:32:48 +01:00
parent e4ab4547db
commit 940c8f05cd

View File

@ -5,7 +5,7 @@
//#include "moench03T1ZmqData.h" //#include "moench03T1ZmqData.h"
//#define DOUBLE_SPH //#define DOUBLE_SPH
//#define MANYFILES //#define MANYFILES
#define WRITE_QUAD //#define WRITE_QUAD
#ifdef DOUBLE_SPH #ifdef DOUBLE_SPH
#include "single_photon_hit_double.h" #include "single_photon_hit_double.h"
#endif #endif
@ -177,6 +177,12 @@ int main(int argc, char *argv[]) {
// if (cl.x>210 && cl.x<240 && cl.y>210 && cl.y<240) { // if (cl.x>210 && cl.x<240 && cl.y>210 && cl.y<240) {
// #endif // #endif
if (cl.x<0 || cl.y<0 || cl.x>NC || cl.y>NR) { if (cl.x<0 || cl.y<0 || cl.x>NC || cl.y>NR) {
if (cl.x<-1 || cl.y<-1 || cl.x>NC+1 || cl.y>NR+1) {
cout <<"**************"<< endl;
cout << cl.x << " " << cl.y << " " << sum << endl;
cl.print();
cout <<"**************"<< endl;
}
} else { } else {
#ifndef FF #ifndef FF
// interp->getInterpolatedPosition(cl.x,cl.y, cl.get_cluster(),int_x, int_y); // interp->getInterpolatedPosition(cl.x,cl.y, cl.get_cluster(),int_x, int_y);