--- PamCut/CollectionActions/Histo2DActions/TofDedxVsBetaHistoAction/TofDedxVsBetaHistoAction.cpp	2012/01/03 08:45:27	1.1.2.2
+++ PamCut/CollectionActions/Histo2DActions/TofDedxVsBetaHistoAction/TofDedxVsBetaHistoAction.cpp	2011/12/20 17:41:49	1.4
@@ -78,64 +78,80 @@
   static unsigned int badLayers;
 
   static int trkSeqNo;
-  static int idx;
   if (_standAlone)
-    trkSeqNo = -1;
+    trkSeqNo = 0;
   else
-    trkSeqNo = event->GetTrack(0)->GetTrkTrack()->GetSeqNo();
-
-  for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
-    if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == trkSeqNo)
-      idx = i;
-  }
+    for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
+      if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo())
+        trkSeqNo = i;
+    }
 
   dEdx = 0.;
   badLayers = 0;
   if ((_layers & S11) == S11) {
-    dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 0, 100);
+    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100);
+    // Patch to exclude S115b
+    ToFTrkVar *tofTrk = event->GetToFLevel2()->GetToFTrkVar(trkSeqNo);
+    bool has8 = false, has9 = false;
+    float dedx9 = 0.;
+    for (Int_t ib = 0; ib < tofTrk->npmtadc; ib++) {
+      if (tofTrk->pmtadc[ib] == 8) {
+        has8 = true;
+      }
+      if (tofTrk->pmtadc[ib] == 9) {
+        has9 = true;
+        dedx9 = tofTrk->dedx[ib];
+      }
+    }
+
+    if (has9) {
+      dEdxLayer -= (dedx9 / ((int) has8 + 1));
+      dEdxLayer *= 2.;
+    }
+    //End of patch
     if (dEdxLayer > 0. && dEdxLayer < 4090)
       dEdx += dEdxLayer;
     else
       badLayers++;
   }
   if ((_layers & S12) == S12) {
-    dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 1, 100);
+    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100);
     if (dEdxLayer > 0. && dEdxLayer < 4090)
       dEdx += dEdxLayer;
     else
       badLayers++;
   }
   if ((_layers & S21) == S21) {
-    dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 2, 100);
+    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100);
     if (dEdxLayer > 0. && dEdxLayer < 4090)
       dEdx += dEdxLayer;
     else
       badLayers++;
   }
   if ((_layers & S22) == S22) {
-    dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 3, 100);
+    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
     if (dEdxLayer > 0. && dEdxLayer < 4090)
       dEdx += dEdxLayer;
     else
       badLayers++;
   }
   if ((_layers & S31) == S31) {
-    dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 4, 100);
+    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 4, 100);
     if (dEdxLayer > 0. && dEdxLayer < 4090)
       dEdx += dEdxLayer;
     else
       badLayers++;
   }
   if ((_layers & S32) == S32) {
-    dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 5, 100);
+    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 5, 100);
     if (dEdxLayer > 0. && dEdxLayer < 4090)
       dEdx += dEdxLayer;
     else
       badLayers++;
   }
   if (badLayers == 0) {
-    dEdx /= _nLayers - badLayers;
-    Fill(event->GetToFLevel2()->GetToFTrkVar(idx)->beta[12], dEdx);
+    dEdx /= _nLayers;
+    Fill(event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12], dEdx);
   }
   else
     _badEvents++;