--- DarthVader/TrackerLevel2/src/TrkLevel2.cpp	2008/11/20 15:06:26	1.48
+++ DarthVader/TrackerLevel2/src/TrkLevel2.cpp	2009/02/03 10:14:27	1.52
@@ -400,6 +400,23 @@
     return (last_plane-first_plane+1);
 }
 /**
+ * Returns the track "lever-arm" on the x+y view, defined as the distance (in planes) between 
+ * the upper and lower x,y (couple) measurements (the maximum value of lever-arm is 6). 
+ */
+Int_t TrkTrack::GetLeverArmXY(){
+    int first_plane = -1;
+    int last_plane  = -1;
+    for(Int_t ip=0; ip<6; ip++){
+	if( XGood(ip)*YGood(ip) && first_plane == -1 )first_plane = ip;
+	if( XGood(ip)*YGood(ip) && first_plane != -1 )last_plane = ip;
+    }
+    if( first_plane == -1 || last_plane == -1){
+	cout<< "Int_t TrkTrack::GetLeverArmXY() -- XGood(ip)*YGood(ip) always false ??? "<<endl; 
+	return 0;
+    }
+    return (last_plane-first_plane+1);
+}
+/**
  * Returns the reduced chi-square of track x-projection
  */
 Float_t  TrkTrack::GetChi2X(){
@@ -610,7 +627,7 @@
 		   4.52043,
 		   4.29926};
     int index;
-    float fact;
+    float fact=0.;
     for(int i=0; i<6; i++) {
 	index = int((fabs(axv[i])+1.)/2.);
 	if(index>10) index=10;
@@ -724,17 +741,17 @@
 	    if(       XGood(i) && !YGood(i) ){
 		track.xm_a[i] = xm[i] - omega * segment;
 		track.ym_a[i] = ym[i] + segment;
-		track.zm_a[i] = zm[i] + beta * segment;
+//	    track.zm_a[i] = zm[i] + beta * segment;//not used yet
 		track.xm_b[i] = xm[i] + omega * segment;
 		track.ym_b[i] = ym[i] - segment;
-		track.zm_b[i] = zm[i] - beta * segment;
+//	    track.zm_b[i] = zm[i] - beta * segment;//not used yet
 	    }else if( !XGood(i) && YGood(i) ){
 		track.xm_a[i] = xm[i] + segment;
 		track.ym_a[i] = ym[i] + omega * segment;
-		track.zm_a[i] = zm[i] - gamma * segment;
+//	    track.zm_a[i] = zm[i] - gamma * segment;//not used yet
 		track.xm_b[i] = xm[i] - segment;
 		track.ym_b[i] = ym[i] - omega * segment;
-		track.zm_b[i] = zm[i] + gamma * segment;
+//	    track.zm_b[i] = zm[i] + gamma * segment;//not used yet
 	    }
 	}
 	
@@ -1259,7 +1276,7 @@
 
 };
 /**
- * \brief Give the average spatial residual 
+ * \brief Give the anerage spatial residual 
  */
 Float_t TrkTrack::GetResidual_av(int ip, int iv){
     //
@@ -1831,7 +1848,8 @@
 
     if(!Track)return 0;
 
-    TRefArray *sorted = new TRefArray();
+    //    TRefArray *sorted = new TRefArray();
+    TRefArray *sorted = NULL;
 	
     TClonesArray &t  = *Track;
 //    TClonesArray &ts = *PhysicalTrack;
@@ -1869,6 +1887,7 @@
 	
 //	    cout << "i** "<< ((TrkTrack *)t[indi])->image << " " << nfiti <<" "<<chi2i<<endl;
 	};
+	if(!sorted)sorted = new TRefArray( TProcessID::GetProcessWithUID(t[indi]));
 	sorted->Add( (TrkTrack*)t[indi] );	
 		
 	m[indi] = 0;