--- PamelaDigitizer/DigitizeTOF.cxx	2009/08/11 14:56:30	1.6
+++ PamelaDigitizer/DigitizeTOF.cxx	2014/01/17 13:08:40	1.8
@@ -189,7 +189,7 @@
   } // NTHTOF < 200
 // ======  ADC ======
   for(Int_t i=0; i<48; i++){
-    if (QevePmt_pC[i] < 800.)  ADCtof[i]= (Int_t)(ADC_pC0A + ADC_pC1A*QevePmt_pC[i] + ADC_pC2A*pow(QevePmt_pC[i],2) + ADC_pC3A*pow(QevePmt_pC[i],3));
+    if (QevePmt_pC[i] <= 800.)  ADCtof[i]= (Int_t)(ADC_pC0A + ADC_pC1A*QevePmt_pC[i] + ADC_pC2A*pow(QevePmt_pC[i],2) + ADC_pC3A*pow(QevePmt_pC[i],3));
     if (QevePmt_pC[i] > 800.)  ADCtof[i]= (Int_t)(ADC_pC0B + ADC_pC1B*QevePmt_pC[i] + ADC_pC2B*pow(QevePmt_pC[i],2) + ADC_pC3B*pow(QevePmt_pC[i],3));
     if (QevePmt_pC[i] > 2485.) ADCtof[i]= (Int_t)(1758. + 0.54*QevePmt_pC[i]);  //assuming a fictional 0.54 ch/pC above ADCsat
     if (ADCtof[i]>ADCsat) ADCtof[i]=ADCsat;
@@ -203,12 +203,21 @@
 
 // ======  build  TDC coincidence  ======
 
+  //
+  for(Int_t i=0; i<48; i++) {
+    if((tdcpmt[i] - c1_S[i]) > 1e-7) {
+      tdcpmt[i] = 0.;
+      tdc[i] = 0.;
+    }
+  }// cycle to introduce a window for tdc
+
+
   Float_t t_coinc = 0;
-  Int_t ilast = 100;
+  //  Int_t ilast = 100;
   for (Int_t ii=0; ii<48;ii++)
     if (tdc[ii] > t_coinc) {
       t_coinc = tdc[ii];
-      ilast = ii;
+      //      ilast = ii;
     }
   
   //     cout<<ilast<<" "<<t_coinc<<endl;
@@ -373,7 +382,7 @@
   UChar_t crcTof=0x00;
   UChar_t *pc=&crcTof, *pc2;
   pc2=pTof;
-  for (Int_t jp=0; jp < 23; jp++){
+  for (Int_t jp=0; jp < 22; jp++){ // cecilia: fixed 23->22
     //crcTof = crc8(...)
     Crc8Tof(pc2++,pc);
     //    printf("%2i --- %x\n",jp,crcTof);