diff --git a/common-tools/clas-tracking/src/main/java/org/jlab/clas/tracking/utilities/MatrixOps.java b/common-tools/clas-tracking/src/main/java/org/jlab/clas/tracking/utilities/MatrixOps.java
index c1ce8072ab..afdf17834d 100644
--- a/common-tools/clas-tracking/src/main/java/org/jlab/clas/tracking/utilities/MatrixOps.java
+++ b/common-tools/clas-tracking/src/main/java/org/jlab/clas/tracking/utilities/MatrixOps.java
@@ -128,7 +128,7 @@ public double[][] MatrixTranspose(Object obj1) {
 
         for(int i = 0; i < r; i++) {
             for (int j = 0; j < c; j++) {
-                result[i][j] = arr1[j][i];
+                result[j][i] = arr1[i][j];
             }
         }
         arr1 = null;
diff --git a/etc/bankdefs/hipo4/dc.json b/etc/bankdefs/hipo4/dc.json
index 4558300006..54ce6acbb6 100644
--- a/etc/bankdefs/hipo4/dc.json
+++ b/etc/bankdefs/hipo4/dc.json
@@ -529,36 +529,47 @@
     "name": "TimeBasedTrkg::TBCovMat",
     "group": 20600,
     "item" : 37,
-    "info": "reconstructed track covariance matrix",
+    "info": "reconstructed track covariance matrix in lab frame",
     "entries": [
         {"name":"id",     "type":"S", "info":"id of the track"},
-        {"name":"C11",    "type":"F", "info":"C11 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C12",    "type":"F", "info":"C12 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C13",    "type":"F", "info":"C13 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C14",    "type":"F", "info":"C14 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C15",    "type":"F", "info":"C15 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C21",    "type":"F", "info":"C21 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C22",    "type":"F", "info":"C22 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C23",    "type":"F", "info":"C23 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C24",    "type":"F", "info":"C24 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C25",    "type":"F", "info":"C25 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C31",    "type":"F", "info":"C31 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C32",    "type":"F", "info":"C32 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C33",    "type":"F", "info":"C33 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C34",    "type":"F", "info":"C34 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C35",    "type":"F", "info":"C35 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C41",    "type":"F", "info":"C41 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C42",    "type":"F", "info":"C42 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C43",    "type":"F", "info":"C43 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C44",    "type":"F", "info":"C44 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C45",    "type":"F", "info":"C45 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C51",    "type":"F", "info":"C51 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C52",    "type":"F", "info":"C52 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C53",    "type":"F", "info":"C53 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C54",    "type":"F", "info":"C54 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C55",    "type":"F", "info":"C55 covariance matrix element at last superlayer used in the fit"}
+        {"name":"C11",    "type":"F", "info":"var(x) at vertex in lab frame"},
+        {"name":"C12",    "type":"F", "info":"cov(x, y) at vertex in lab frame"},
+        {"name":"C13",    "type":"F", "info":"cov(x, z) at vertex in lab frame"}, 
+        {"name":"C14",    "type":"F", "info":"cov(x, p_theta) at vertex in lab frame"},
+        {"name":"C15",    "type":"F", "info":"cov(x, p_phi) at vertex in lab frame"},
+        {"name":"C16",    "type":"F", "info":"cov(x, p) at vertex in lab frame"},        
+        {"name":"C21",    "type":"F", "info":"cov(y, x) at vertex in lab frame"},
+        {"name":"C22",    "type":"F", "info":"var(y) at vertex in lab frame"},
+        {"name":"C23",    "type":"F", "info":"cov(y, z) at vertex in lab frame"}, 
+        {"name":"C24",    "type":"F", "info":"cov(y, p_theta) at vertex in lab frame"},
+        {"name":"C25",    "type":"F", "info":"cov(y, p_phi) at vertex in lab frame"},
+        {"name":"C26",    "type":"F", "info":"cov(y, p) at vertex in lab frame"},        
+        {"name":"C31",    "type":"F", "info":"cov(z, x) at vertex in lab frame"},
+        {"name":"C32",    "type":"F", "info":"cov(z, y) at vertex in lab frame"},
+        {"name":"C33",    "type":"F", "info":"var(z) at vertex in lab frame"}, 
+        {"name":"C34",    "type":"F", "info":"cov(z, p_theta) at vertex in lab frame"},
+        {"name":"C35",    "type":"F", "info":"cov(z, p_phi) at vertex in lab frame"},
+        {"name":"C36",    "type":"F", "info":"cov(z, p) at vertex in lab frame"},        
+        {"name":"C41",    "type":"F", "info":"cov(p_theta, x) at vertex in lab frame"},
+        {"name":"C42",    "type":"F", "info":"cov(p_theta, y) at vertex in lab frame"},
+        {"name":"C43",    "type":"F", "info":"cov(p_theta, z) at vertex in lab frame"}, 
+        {"name":"C44",    "type":"F", "info":"var(p_theta) at vertex in lab frame"},
+        {"name":"C45",    "type":"F", "info":"cov(p_theta, p_phi) at vertex in lab frame"},
+        {"name":"C46",    "type":"F", "info":"cov(p_theta, p) at vertex in lab frame"},        
+        {"name":"C51",    "type":"F", "info":"cov(p_phi, x) at vertex in lab frame"},
+        {"name":"C52",    "type":"F", "info":"cov(p_phi, y) at vertex in lab frame"},
+        {"name":"C53",    "type":"F", "info":"cov(p_phi, z) at vertex in lab frame"}, 
+        {"name":"C54",    "type":"F", "info":"cov(p_phi, p_theta) at vertex in lab frame"},
+        {"name":"C55",    "type":"F", "info":"var(p_phi) at vertex in lab frame"},
+        {"name":"C56",    "type":"F", "info":"cov(p_phi, p) at vertex in lab frame"},        
+        {"name":"C61",    "type":"F", "info":"cov(p, x) at vertex in lab frame"},
+        {"name":"C62",    "type":"F", "info":"cov(p, y) at vertex in lab frame"},
+        {"name":"C63",    "type":"F", "info":"cov(p, z) at vertex in lab frame"}, 
+        {"name":"C64",    "type":"F", "info":"cov(p, p_theta) at vertex in lab frame"},
+        {"name":"C65",    "type":"F", "info":"cov(p, p_phi) at vertex in lab frame"},
+        {"name":"C66",    "type":"F", "info":"var(p) at vertex in lab frame"}
     ]
-  },
+  },  
   {
     "name": "TimeBasedTrkg::Trajectory",
     "group": 20600,
diff --git a/etc/bankdefs/hipo4/dcnn.json b/etc/bankdefs/hipo4/dcnn.json
index 74276a745d..145f3cba6b 100644
--- a/etc/bankdefs/hipo4/dcnn.json
+++ b/etc/bankdefs/hipo4/dcnn.json
@@ -390,36 +390,47 @@
     "name": "TimeBasedTrkg::AICovMat",
     "group": 20600,
     "item" : 77,
-    "info": "reconstructed track covariance matrix",
+    "info": "reconstructed track covariance matrix in lab frame",
     "entries": [
         {"name":"id",     "type":"S", "info":"id of the track"},
-        {"name":"C11",    "type":"F", "info":"C11 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C12",    "type":"F", "info":"C12 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C13",    "type":"F", "info":"C13 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C14",    "type":"F", "info":"C14 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C15",    "type":"F", "info":"C15 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C21",    "type":"F", "info":"C21 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C22",    "type":"F", "info":"C22 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C23",    "type":"F", "info":"C23 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C24",    "type":"F", "info":"C24 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C25",    "type":"F", "info":"C25 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C31",    "type":"F", "info":"C31 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C32",    "type":"F", "info":"C32 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C33",    "type":"F", "info":"C33 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C34",    "type":"F", "info":"C34 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C35",    "type":"F", "info":"C35 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C41",    "type":"F", "info":"C41 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C42",    "type":"F", "info":"C42 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C43",    "type":"F", "info":"C43 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C44",    "type":"F", "info":"C44 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C45",    "type":"F", "info":"C45 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C51",    "type":"F", "info":"C51 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C52",    "type":"F", "info":"C52 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C53",    "type":"F", "info":"C53 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C54",    "type":"F", "info":"C54 covariance matrix element at last superlayer used in the fit"},
-        {"name":"C55",    "type":"F", "info":"C55 covariance matrix element at last superlayer used in the fit"}
+        {"name":"C11",    "type":"F", "info":"var(x) at vertex in lab frame"},
+        {"name":"C12",    "type":"F", "info":"cov(x, y) at vertex in lab frame"},
+        {"name":"C13",    "type":"F", "info":"cov(x, z) at vertex in lab frame"}, 
+        {"name":"C14",    "type":"F", "info":"cov(x, p_theta) at vertex in lab frame"},
+        {"name":"C15",    "type":"F", "info":"cov(x, p_phi) at vertex in lab frame"},
+        {"name":"C16",    "type":"F", "info":"cov(x, p) at vertex in lab frame"},        
+        {"name":"C21",    "type":"F", "info":"cov(y, x) at vertex in lab frame"},
+        {"name":"C22",    "type":"F", "info":"var(y) at vertex in lab frame"},
+        {"name":"C23",    "type":"F", "info":"cov(y, z) at vertex in lab frame"}, 
+        {"name":"C24",    "type":"F", "info":"cov(y, p_theta) at vertex in lab frame"},
+        {"name":"C25",    "type":"F", "info":"cov(y, p_phi) at vertex in lab frame"},
+        {"name":"C26",    "type":"F", "info":"cov(y, p) at vertex in lab frame"},        
+        {"name":"C31",    "type":"F", "info":"cov(z, x) at vertex in lab frame"},
+        {"name":"C32",    "type":"F", "info":"cov(z, y) at vertex in lab frame"},
+        {"name":"C33",    "type":"F", "info":"var(z) at vertex in lab frame"}, 
+        {"name":"C34",    "type":"F", "info":"cov(z, p_theta) at vertex in lab frame"},
+        {"name":"C35",    "type":"F", "info":"cov(z, p_phi) at vertex in lab frame"},
+        {"name":"C36",    "type":"F", "info":"cov(z, p) at vertex in lab frame"},        
+        {"name":"C41",    "type":"F", "info":"cov(p_theta, x) at vertex in lab frame"},
+        {"name":"C42",    "type":"F", "info":"cov(p_theta, y) at vertex in lab frame"},
+        {"name":"C43",    "type":"F", "info":"cov(p_theta, z) at vertex in lab frame"}, 
+        {"name":"C44",    "type":"F", "info":"var(p_theta) at vertex in lab frame"},
+        {"name":"C45",    "type":"F", "info":"cov(p_theta, p_phi) at vertex in lab frame"},
+        {"name":"C46",    "type":"F", "info":"cov(p_theta, p) at vertex in lab frame"},        
+        {"name":"C51",    "type":"F", "info":"cov(p_phi, x) at vertex in lab frame"},
+        {"name":"C52",    "type":"F", "info":"cov(p_phi, y) at vertex in lab frame"},
+        {"name":"C53",    "type":"F", "info":"cov(p_phi, z) at vertex in lab frame"}, 
+        {"name":"C54",    "type":"F", "info":"cov(p_phi, p_theta) at vertex in lab frame"},
+        {"name":"C55",    "type":"F", "info":"var(p_phi) at vertex in lab frame"},
+        {"name":"C56",    "type":"F", "info":"cov(p_phi, p) at vertex in lab frame"},        
+        {"name":"C61",    "type":"F", "info":"cov(p, x) at vertex in lab frame"},
+        {"name":"C62",    "type":"F", "info":"cov(p, y) at vertex in lab frame"},
+        {"name":"C63",    "type":"F", "info":"cov(p, z) at vertex in lab frame"}, 
+        {"name":"C64",    "type":"F", "info":"cov(p, p_theta) at vertex in lab frame"},
+        {"name":"C65",    "type":"F", "info":"cov(p, p_phi) at vertex in lab frame"},
+        {"name":"C66",    "type":"F", "info":"var(p) at vertex in lab frame"}
     ]
-  },
+  },  
   {
     "name": "TimeBasedTrkg::AITrajectory",
     "group": 20600,
diff --git a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/Banks.java b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/Banks.java
index aac1368be3..b677dd8f90 100644
--- a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/Banks.java
+++ b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/Banks.java
@@ -136,7 +136,7 @@ public String getTrajBank() {
 
     public String getCovmatBank() {
         return this.getOutputBank("CovMat");
-    }
+    }        
     
     public String getRecEventBank() {
         return this.getRecBank("Event");
diff --git a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java
index 529be90740..01e734df15 100644
--- a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java
+++ b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java
@@ -446,50 +446,69 @@ public DataBank fillHBTracksBank(DataEvent event, List<Track> candlist) {
     public DataBank fillHBTrajectoryBank(DataEvent event, List<Track> candlist) {
         return this.fillTrajectoryBank(event, candlist);
     }
-    
+        
     /**
      *
      * @param event hipo event
      * @param candlist tracks
-     * @return covariance matrix from HB fits to be used for starting TB tracking
+     * @return covariance matrix for momentum and vertex in lab frame
      */
-    private DataBank fillTrackCovMatBank(DataEvent event, List<Track> candlist) {
+    private DataBank fillTrackCovMatLabBank(DataEvent event, List<Track> candlist) {
 
         DataBank bank = event.createBank(bankNames.getCovmatBank(), candlist.size());
 
         for (int i = 0; i < candlist.size(); i++) {
             bank.setShort("id", i, (short) candlist.get(i).get_Id());
-            if(candlist.get(i).get_CovMat()!=null) {
-                bank.setFloat("C11", i, (float) candlist.get(i).get_CovMat().get(0, 0));
-                bank.setFloat("C12", i, (float) candlist.get(i).get_CovMat().get(0, 1));
-                bank.setFloat("C13", i, (float) candlist.get(i).get_CovMat().get(0, 2));
-                bank.setFloat("C14", i, (float) candlist.get(i).get_CovMat().get(0, 3));
-                bank.setFloat("C15", i, (float) candlist.get(i).get_CovMat().get(0, 4));
-                bank.setFloat("C21", i, (float) candlist.get(i).get_CovMat().get(1, 0));
-                bank.setFloat("C22", i, (float) candlist.get(i).get_CovMat().get(1, 1));
-                bank.setFloat("C23", i, (float) candlist.get(i).get_CovMat().get(1, 2));
-                bank.setFloat("C24", i, (float) candlist.get(i).get_CovMat().get(1, 3));
-                bank.setFloat("C25", i, (float) candlist.get(i).get_CovMat().get(1, 4));
-                bank.setFloat("C31", i, (float) candlist.get(i).get_CovMat().get(2, 0));
-                bank.setFloat("C32", i, (float) candlist.get(i).get_CovMat().get(2, 1));
-                bank.setFloat("C33", i, (float) candlist.get(i).get_CovMat().get(2, 2));
-                bank.setFloat("C34", i, (float) candlist.get(i).get_CovMat().get(2, 3));
-                bank.setFloat("C35", i, (float) candlist.get(i).get_CovMat().get(2, 4));
-                bank.setFloat("C41", i, (float) candlist.get(i).get_CovMat().get(3, 0));
-                bank.setFloat("C42", i, (float) candlist.get(i).get_CovMat().get(3, 1));
-                bank.setFloat("C43", i, (float) candlist.get(i).get_CovMat().get(3, 2));
-                bank.setFloat("C44", i, (float) candlist.get(i).get_CovMat().get(3, 3));
-                bank.setFloat("C45", i, (float) candlist.get(i).get_CovMat().get(3, 4));
-                bank.setFloat("C51", i, (float) candlist.get(i).get_CovMat().get(4, 0));
-                bank.setFloat("C52", i, (float) candlist.get(i).get_CovMat().get(4, 1));
-                bank.setFloat("C53", i, (float) candlist.get(i).get_CovMat().get(4, 2));
-                bank.setFloat("C54", i, (float) candlist.get(i).get_CovMat().get(4, 3));
-                bank.setFloat("C55", i, (float) candlist.get(i).get_CovMat().get(4, 4));
+            if(candlist.get(i).get_CMInLab()!=null) {
+                double[][] CM = candlist.get(i).get_CMInLab();
+                bank.setFloat("C11", i, (float) CM[0][0]);
+                bank.setFloat("C12", i, (float) CM[0][1]);
+                bank.setFloat("C13", i, (float) CM[0][2]);
+                bank.setFloat("C14", i, (float) CM[0][3]);
+                bank.setFloat("C15", i, (float) CM[0][4]);
+                bank.setFloat("C16", i, (float) CM[0][5]);
+                
+                bank.setFloat("C21", i, (float) CM[1][0]);
+                bank.setFloat("C22", i, (float) CM[1][1]);
+                bank.setFloat("C23", i, (float) CM[1][2]);
+                bank.setFloat("C24", i, (float) CM[1][3]);
+                bank.setFloat("C25", i, (float) CM[1][4]);
+                bank.setFloat("C26", i, (float) CM[1][5]);
+                
+                bank.setFloat("C31", i, (float) CM[2][0]);
+                bank.setFloat("C32", i, (float) CM[2][1]);
+                bank.setFloat("C33", i, (float) CM[2][2]);
+                bank.setFloat("C34", i, (float) CM[2][3]);
+                bank.setFloat("C35", i, (float) CM[2][4]);
+                bank.setFloat("C36", i, (float) CM[2][5]);
+                
+                bank.setFloat("C41", i, (float) CM[3][0]);
+                bank.setFloat("C42", i, (float) CM[3][1]);
+                bank.setFloat("C43", i, (float) CM[3][2]);
+                bank.setFloat("C44", i, (float) CM[3][3]);
+                bank.setFloat("C45", i, (float) CM[3][4]);
+                bank.setFloat("C46", i, (float) CM[3][5]);
+                
+                bank.setFloat("C51", i, (float) CM[4][0]);
+                bank.setFloat("C52", i, (float) CM[4][1]);
+                bank.setFloat("C53", i, (float) CM[4][2]);
+                bank.setFloat("C54", i, (float) CM[4][3]);
+                bank.setFloat("C55", i, (float) CM[4][4]);
+                bank.setFloat("C56", i, (float) CM[4][5]);
+                
+                bank.setFloat("C61", i, (float) CM[5][0]);
+                bank.setFloat("C62", i, (float) CM[5][1]);
+                bank.setFloat("C63", i, (float) CM[5][2]);
+                bank.setFloat("C64", i, (float) CM[5][3]);
+                bank.setFloat("C65", i, (float) CM[5][4]);
+                bank.setFloat("C66", i, (float) CM[5][5]);
+               
             }
         }
         //bank.show();
         return bank;
     }
+    
      /**
      *
      * @param event the EvioEvent
@@ -939,7 +958,7 @@ public void fillAllTBBanks(DataEvent event, List<FittedHit> fhits, List<FittedCl
                     this.fillTBCrossesBank(event, crosses),
                     this.fillTBTracksBank(event, trkcands),
                     this.fillTrajectoryBank(event, trkcands),
-                    this.fillTrackCovMatBank(event, trkcands)
+                    this.fillTrackCovMatLabBank(event, trkcands)
                     );
         }
         if (crosses != null && trkcands == null) {
diff --git a/reconstruction/dc/src/main/java/org/jlab/rec/dc/track/Track.java b/reconstruction/dc/src/main/java/org/jlab/rec/dc/track/Track.java
index d9c8d75825..96b07460fe 100644
--- a/reconstruction/dc/src/main/java/org/jlab/rec/dc/track/Track.java
+++ b/reconstruction/dc/src/main/java/org/jlab/rec/dc/track/Track.java
@@ -11,6 +11,8 @@
 import org.jlab.rec.dc.segment.Segment;
 import org.jlab.rec.dc.trajectory.StateVec;
 import org.jlab.rec.dc.trajectory.Trajectory;
+import org.jlab.clas.tracking.utilities.MatrixOps;
+import org.jlab.clas.tracking.utilities.MatrixOps.Libr;
 
 /**
  * A class representing track candidates in the DC.  A track has a trajectory represented by an ensemble of geometrical state vectors along its path, 
@@ -42,6 +44,7 @@ public void setFinalStateVec(StateVec finalStateVec) {
     private int _Q;
     private double _P;
     private Matrix _CovMat;
+    private double[][] _CMInLab = new double[6][6];
 
     private Point3D _Region3CrossPoint;
     private Point3D _Region3CrossDir;
@@ -396,6 +399,90 @@ public void set_CovMat(Matrix _CovMat) {
         this._CovMat = _CovMat;
     }
     
+    /**
+     * 
+     * @return covariance matrix in lab frame
+     */
+    public double[][] get_CMInLab() {
+        return _CMInLab;
+    }
+    
+     /**
+     * 
+     * @param _CMInLab covariance matrix in lab frame
+     */
+    public void set_CMInLab(double[][] _CMInLab) {
+        this._CMInLab = _CMInLab;
+    }
+    
+    public void transCMToGlobal(){
+        MatrixOps mo = new MatrixOps(Libr.EJML);
+        
+        int sector = this.getSector();
+        
+        int charge = this.get_Q();
+        double p = this.get_P();                
+        
+        Vector3D momGlobal = this.get_pAtOrig();
+        Cross C = this.get(this.size() - 1);
+        Point3D momLocal = C.getCoordsInTiltedSector(momGlobal.x(), momGlobal.y(), momGlobal.z());
+        if(momLocal.z() == 0) set_CMInLab(new double[6][6]);
+        double tx = momLocal.x()/momLocal.z();
+        double ty = momLocal.y()/momLocal.z();
+        double theta = Math.atan(Math.sqrt(tx*tx + ty*ty));
+        double phi = Math.atan2(ty, tx);
+        
+        //jm: Jacobi matrix for transformation from (x, y, tx, ty, Q) to (x, y, z, θ, φ, p) in the tilted sector frame
+        double dpdQ = -charge * p * p;
+        double dphidtx = -ty / (tx*tx + ty*ty);
+        double dphidty = tx / (tx*tx + ty*ty);
+        double dthetadtx = tx/(1 + tx*tx + ty*ty)/Math.sqrt(tx*tx + ty*ty);
+        double dthetadty = ty/(1 + tx*tx + ty*ty)/Math.sqrt(tx*tx + ty*ty);
+        double[][] jm = {{1, 0, 0, 0, 0}, {0, 1, 0, 0, 0}, {0, 0, 0, 0, 0},
+                         {0, 0, dthetadtx, dthetadty, 0}, {0, 0, dphidtx, dphidty, 0}, {0, 0, 0, 0, dpdQ}
+                        };
+        double[][] jmT = mo.MatrixTranspose(jm);
+        
+        //jt: Jacobi matrix for transformation from (x, y, z, θ, φ, p) in the tilted sector frame to (x’, y’, z’, θ’, φ’, p’) in the sector frame (rotate 250 around y)
+        double numerator11 = Constants.SIN25 * Math.cos(theta) * Math.cos(phi) + Constants.COS25 * Math.sin(theta);
+        double numerator12 = -Constants.SIN25 * Math.sin(theta) * Math.sin(phi);
+        double denominator1 = Math.sqrt(1 - Math.pow(Constants.SIN25*Math.sin(theta)*Math.cos(phi) - Constants.COS25*Math.cos(theta), 2));         
+        double numerator21 = Constants.SIN25 * Math.sin(phi);
+        double numerator22 = Constants.COS25 * Math.pow(Math.sin(theta),2) + Constants.SIN25 * Math.sin(theta) * Math.cos(theta) * Math.cos(phi);
+        double denoninator2 = Math.pow(Math.sin(theta) * Math.sin(phi), 2) + Math.pow(Constants.COS25 * Math.sin(theta) * Math.cos(phi) +  Constants.SIN25 * Math.cos(theta),2);
+        double[][] jt = {{Constants.COS25, 0, Constants.SIN25, 0, 0, 0}, {0, 1, 0, 0, 0, 0}, {-Constants.SIN25, 0, Constants.COS25, 0, 0, 0}, 
+                         {0, 0, 0, numerator11/denominator1, numerator12/denominator1, 0}, {0, 0, 0, numerator21/denoninator2, numerator22/denoninator2, 0}, {0, 0, 0, 0, 0, 1}
+                        };
+        double[][] jtT = mo.MatrixTranspose(jt);
+        
+        //js: Jacobi matrix for transformation from  (x’, y’, z’, θ’, φ’, p’)  to (x’’, y’’, z’’, θ’’, φ’’, p’’)  in the lab frame (rotate around z’)
+        double[][] js = {{Constants.COSSECTOR60[sector - 1], -Constants.SINSECTOR60[sector - 1], 0, 0, 0, 0}, {Constants.SINSECTOR60[sector - 1], Constants.COSSECTOR60[sector - 1], 0, 0, 0, 0}, {0, 0, 1, 0, 0, 0}, 
+                         {0, 0, 0, 1, 0, 0}, {0, 0, 0, 0, 1, 0}, {0, 0, 0, 0, 0, 1}
+                        };
+        double[][] jsT = mo.MatrixTranspose(js);
+        
+        Matrix CM = this.get_CovMat();
+        double[][] CMOrig = {
+                                {CM.get(0, 0), CM.get(0, 1), CM.get(0, 2), CM.get(0, 3), CM.get(0, 4)},
+                                {CM.get(1, 0), CM.get(1, 1), CM.get(1, 2), CM.get(1, 3), CM.get(1, 4)},
+                                {CM.get(2, 0), CM.get(2, 1), CM.get(2, 2), CM.get(2, 3), CM.get(2, 4)},
+                                {CM.get(3, 0), CM.get(3, 1), CM.get(3, 2), CM.get(3, 3), CM.get(3, 4)},
+                                {CM.get(4, 0), CM.get(4, 1), CM.get(4, 2), CM.get(4, 3), CM.get(4, 4)}
+            
+                            };
+                
+        double[][] transCM1 = mo.MatrixMultiplication(jm, CMOrig);
+        double[][] transCM2 = mo.MatrixMultiplication(transCM1, jmT);
+        
+        double[][] transCM3 = mo.MatrixMultiplication(jt, transCM2);
+        double[][] transCM4 = mo.MatrixMultiplication(transCM3, jtT);
+        
+        double[][] transCM5 = mo.MatrixMultiplication(js, transCM4);
+        double[][] transCM6 = mo.MatrixMultiplication(transCM5, jsT);
+                      
+        set_CMInLab(transCM6);
+     }
+                
     /**
      * 
      * @param fitConvergenceStatus fit convergence status 0 if OK, 1 if the fit exits before converging
diff --git a/reconstruction/dc/src/main/java/org/jlab/service/dc/DCTBEngine.java b/reconstruction/dc/src/main/java/org/jlab/service/dc/DCTBEngine.java
index 93abb7cf23..6a09909377 100644
--- a/reconstruction/dc/src/main/java/org/jlab/service/dc/DCTBEngine.java
+++ b/reconstruction/dc/src/main/java/org/jlab/service/dc/DCTBEngine.java
@@ -263,8 +263,9 @@ public boolean processDataEvent(DataEvent event) {
                     }
                                         
                     // get CovMat at vertex
-                    Point3D VTCS = crosses.get(0).getCoordsInTiltedSector(TrackArray1.get_Vtx0().x(), TrackArray1.get_Vtx0().y(), TrackArray1.get_Vtx0().z());
-                    TrackArray1.set_CovMat(kFZRef.propagateToVtx(crosses.get(0).get_Sector(), VTCS.z()));
+                    Point3D VTCS = TrackArray1.get(TrackArray1.size()-1).getCoordsInTiltedSector(TrackArray1.get_Vtx0().x(), TrackArray1.get_Vtx0().y(), TrackArray1.get_Vtx0().z());
+                    TrackArray1.set_CovMat(kFZRef.propagateToVtx(TrackArray1.get(TrackArray1.size()-1).get_Sector(), VTCS.z()));                    
+                    TrackArray1.transCMToGlobal();
                     
                     double deltaPathToVtx =  kFZRef.getDeltaPathToVtx(TrackArray1.get(TrackArray1.size()-1).get_Sector(), VTCS.z());                                
                     List<org.jlab.rec.dc.trajectory.StateVec> kfStateVecsAlongTrajectory = setKFStateVecsAlongTrajectory(kFZRef, deltaPathToVtx);
@@ -306,9 +307,10 @@ public boolean processDataEvent(DataEvent event) {
                         continue;
                     }
 
-                    // get CovMat at vertex
-                    Point3D VTCS = crosses.get(0).getCoordsInTiltedSector(TrackArray1.get_Vtx0().x(), TrackArray1.get_Vtx0().y(), TrackArray1.get_Vtx0().z());
-                    TrackArray1.set_CovMat(kFZRef.propagateToVtx(crosses.get(0).get_Sector(), VTCS.z()));
+                    // get CovMat at vertex                    
+                    Point3D VTCS = TrackArray1.get(TrackArray1.size()-1).getCoordsInTiltedSector(TrackArray1.get_Vtx0().x(), TrackArray1.get_Vtx0().y(), TrackArray1.get_Vtx0().z());
+                    TrackArray1.set_CovMat(kFZRef.propagateToVtx(TrackArray1.get(TrackArray1.size()-1).get_Sector(), VTCS.z()));                    
+                    TrackArray1.transCMToGlobal();
                     
                     double deltaPathToVtx =  kFZRef.getDeltaPathToVtx(TrackArray1.get(TrackArray1.size()-1).get_Sector(), VTCS.z());                                
                     List<org.jlab.rec.dc.trajectory.StateVec> kfStateVecsAlongTrajectory = setKFStateVecsAlongTrajectory(kFZRef, deltaPathToVtx);