--- tof/flight/ToFReprocessing/ToFreproc.cpp 2012/04/12 12:04:20 1.6 +++ tof/flight/ToFReprocessing/ToFreproc.cpp 2014/01/17 15:16:24 1.7 @@ -255,12 +255,14 @@ if ( !(iev%100000) ) printf(" ATIME %u re %u \n",(Int_t)tabs,(UInt_t)iev); //================================================================== - Float_t dedx_corr_m[2000][48],dedx_corr[48]; + Float_t dedx_corr_m[2000][48];//,dedx_corr[48]; Double_t mtime[2000],t1,t2,tm; - Float_t yhelp1,yhelp2,slope,inter,thelp1,thelp2; + // Float_t yhelp1,yhelp2,slope,inter,thelp1,thelp2; + Float_t yhelp1,yhelp2; Float_t xmean1,xwidth1; - Int_t ical,ii,j,jj; + Int_t ical=0; + Int_t ii,j,jj; if (iev==0) { @@ -323,8 +325,8 @@ //== interpolate betwen time limits //================================================================== - thelp1 = mtime[ical]; - thelp2 = mtime[ical+1]; + // thelp1 = mtime[ical]; + // thelp2 = mtime[ical+1]; for (ii=0; ii<48;ii++) { yhelp1 = fabs(dedx_corr_m[ical][ii]); @@ -333,9 +335,9 @@ yhelp2 = fabs(dedx_corr_m[ical+1][ii]); // yhelp2 = 6.; if ( yhelp2 < 0.1 ) yhelp2 = 4.; - slope = (yhelp2-yhelp1)/(thelp2-thelp1); - inter = yhelp1 - slope*thelp1; - dedx_corr[ii] = slope*tabs + inter; + // slope = (yhelp2-yhelp1)/(thelp2-thelp1); + // inter = yhelp1 - slope*thelp1; + // dedx_corr[ii] = slope*tabs + inter; // if (ii==0) cout<