Index: PndEmcXtal.h
===================================================================
--- PndEmcXtal.h	(revision 7277)
+++ PndEmcXtal.h	(working copy)
@@ -52,6 +52,7 @@
   const TVector3& centre() const;
   const TVector3& frontCentre() const;
   const TVector3& normalToFrontFace() const ;
+  const TVector3& axisVector() const;
   const TGeoTrap& geometry() const;
   const TGeoRotation& rotation() const;
 
@@ -66,6 +67,7 @@
   TVector3 fCentre;
   TVector3 fFrontCentre;
   TVector3 fNormalToFrontFace;
+  TVector3 fAxis;
 
   ClassDef(PndEmcXtal,1)
 }; 
Index: PndEmcXtal.cxx
===================================================================
--- PndEmcXtal.cxx	(revision 7277)
+++ PndEmcXtal.cxx	(working copy)
@@ -70,12 +70,14 @@
 
         fFrontCentre = pos + axis_vector;
 		fNormalToFrontFace = -1.0 * normal_vector;
+        fAxis = -1.0 * axis_vector.Unit();
     } else {
         //cout << "Crystal orientation in module " << id->Index() / 100000000
         //    << " -dz" << endl;
 
         fFrontCentre = pos - axis_vector;
 		fNormalToFrontFace = +1.0 * normal_vector;
+        fAxis = +1.0 * axis_vector.Unit();
     }
 }
 
@@ -137,6 +139,12 @@
   return fNormalToFrontFace;
 }
 
+const TVector3&
+PndEmcXtal::axisVector() const
+{
+  return fAxis;
+}
+
 const TGeoTrap&
 PndEmcXtal::geometry() const
 { 
Index: ../EmcData/PndEmcDigi.cxx
===================================================================
--- ../EmcData/PndEmcDigi.cxx	(revision 7277)
+++ ../EmcData/PndEmcDigi.cxx	(working copy)
@@ -117,7 +117,7 @@
 	//cout << "depth" << endl;
 	
 	TVector3 pos = xtal->frontCentre();
-	TVector3 norm=xtal->normalToFrontFace();
+	TVector3 norm = xtal->axisVector();
 	norm*=fPositionDepth;
 	pos += norm;