Index: external/TrackCovariance/TrkUtil.cc
===================================================================
--- external/TrackCovariance/TrkUtil.cc	(revision 4df491e8c4fe6b4befe6b5766353136a3580f85a)
+++ external/TrackCovariance/TrkUtil.cc	(revision 5862d1fa07d0f47f91791bc8ab5ef06774495194)
@@ -233,4 +233,11 @@
 }
 //
+
+//
+void TrkUtil::SetBfield(Double_t Bz)
+{
+	fBz = Bz;
+}
+
 // Setup chamber volume
 void TrkUtil::SetDchBoundaries(Double_t Rmin, Double_t Rmax, Double_t Zmin, Double_t Zmax)
@@ -267,10 +274,10 @@
 		//	<< ", C= " << C << ", z0= " << z0 << ", ct= " << ct << std::endl;
 		//
-		// Track length per unit phase change 
+		// Track length per unit phase change
 		Double_t Scale = sqrt(1.0 + ct*ct) / (2.0*TMath::Abs(C));
 		//
 		// Find intersections with chamber boundaries
 		//
-		Double_t phRin = 0.0;			// phase of inner cylinder 
+		Double_t phRin = 0.0;			// phase of inner cylinder
 		Double_t phRin2= 0.0;			// phase of inner cylinder intersection (2nd branch)
 		Double_t phRhi = 0.0;			// phase of outer cylinder intersection
@@ -287,5 +294,5 @@
 			//std::cout << "Rin intersection: ph = " << ph<<", z= "<<z << std::endl;
 
-			if (z < fZmax && z > fZmin)	phRin = TMath::Abs(ph);	// Intersection inside chamber volume	
+			if (z < fZmax && z > fZmin)	phRin = TMath::Abs(ph);	// Intersection inside chamber volume
 			//
 			// Include second branch of loopers
@@ -301,5 +308,5 @@
 			Double_t ph = 2 * asin(C*sqrt((fRmax*fRmax - D*D) / (1.0 + 2.0*C*D)));
 			Double_t z = z0 + ct*ph / (2.0*C);
-			if (z < fZmax && z > fZmin)	phRhi = TMath::Abs(ph);	// Intersection inside chamber volume	
+			if (z < fZmax && z > fZmin)	phRhi = TMath::Abs(ph);	// Intersection inside chamber volume
 		}
 		//  ... with left wall
@@ -309,5 +316,5 @@
 			Double_t ph = 2.0*C*Zdir;
 			Double_t Rint = sqrt(D*D + (1.0 + 2.0*C*D)*pow(sin(ph / 2), 2) / (C*C));
-			if (Rint < fRmax && Rint > fRmin)	phZmn = TMath::Abs(ph);	// Intersection inside chamber volume	
+			if (Rint < fRmax && Rint > fRmin)	phZmn = TMath::Abs(ph);	// Intersection inside chamber volume
 		}
 		//  ... with right wall
@@ -317,5 +324,5 @@
 			Double_t ph = 2.0*C*Zdir;
 			Double_t Rint = sqrt(D*D + (1.0 + 2.0*C*D)*pow(sin(ph / 2), 2) / (C*C));
-			if (Rint < fRmax && Rint > fRmin)	phZmx = TMath::Abs(ph);	// Intersection inside chamber volume	
+			if (Rint < fRmax && Rint > fRmin)	phZmx = TMath::Abs(ph);	// Intersection inside chamber volume
 		}
 		//
@@ -390,5 +397,5 @@
 //
 //
-Double_t TrkUtil::Nclusters(Double_t begam) 
+Double_t TrkUtil::Nclusters(Double_t begam)
 {
 	Int_t Opt = fGasSel;
@@ -465,5 +472,5 @@
     	{
 		case 0: std::copy(ncl_He_Iso, ncl_He_Iso + Npt, ncl);	// He-Isobutane
-		break;							
+		break;
 		case 1: std::copy(ncl_He, ncl_He + Npt, ncl);		// pure He
 		break;
@@ -493,6 +500,6 @@
 	TVectorD coeff = Xval * y;
 	Double_t interp = coeff[0] + coeff[1] * begam + coeff[2] * begam * begam;
-	//std::cout << "val1= (" <<x(0)<<", "<< y(0) << "), val2= (" 
-	//	<<x(1)<<", "<< y(1) << "), val3= (" 
+	//std::cout << "val1= (" <<x(0)<<", "<< y(0) << "), val2= ("
+	//	<<x(1)<<", "<< y(1) << "), val3= ("
 	//	<<x(2)<<", "<< y(2)
 	//	<< "), result= (" <<begam<<", "<< interp<<")" << std::endl;
Index: external/TrackCovariance/TrkUtil.h
===================================================================
--- external/TrackCovariance/TrkUtil.h	(revision 4df491e8c4fe6b4befe6b5766353136a3580f85a)
+++ external/TrackCovariance/TrkUtil.h	(revision 5862d1fa07d0f47f91791bc8ab5ef06774495194)
@@ -25,5 +25,4 @@
 	// Service routines
 	//
-	void SetBfield(Double_t Bz) { fBz = Bz; }
 	TVectorD XPtoPar(TVector3 x, TVector3 p, Double_t Q);
 	TVector3 ParToP(TVectorD Par);
@@ -54,5 +53,5 @@
 		Double_t c = 2.99792458e8;	// speed of light m/sec
 		//return TMath::C()*1.0e-9;	// Incompatible with root5
-		return c*1.0e-9; 		// Reduced speed of light	
+		return c*1.0e-9; 		// Reduced speed of light
 	}
 	//
@@ -71,5 +70,6 @@
 	// Cluster counting in gas
 	//
-	// Define gas volume (units = meters) 
+	void SetBfield(Double_t Bz);
+	// Define gas volume (units = meters)
 	void SetDchBoundaries(Double_t Rmin, Double_t Rmax, Double_t Zmin, Double_t Zmax);
 	// Gas mixture selection
