--- DarthVader/TrackerLevel2/src/TrkLevel2.cpp	2007/08/31 14:56:51	1.40
+++ DarthVader/TrackerLevel2/src/TrkLevel2.cpp	2007/11/27 15:28:57	1.42
@@ -61,9 +61,10 @@
 
     };
 
-    TrkParams::SetTrackingMode();
-    TrkParams::SetPrecisionFactor();
-    TrkParams::SetStepMin();
+//     TrkParams::SetTrackingMode();
+//     TrkParams::SetPrecisionFactor();
+//     TrkParams::SetStepMin();
+    TrkParams::SetMiniDefault();
     TrkParams::SetPFA();
 
 };
@@ -105,9 +106,10 @@
 	ypu[ip]      = t.ypu[ip];   
     };
 
-    TrkParams::SetTrackingMode();
-    TrkParams::SetPrecisionFactor(); 
-    TrkParams::SetStepMin();   
+//     TrkParams::SetTrackingMode();
+//     TrkParams::SetPrecisionFactor(); 
+//     TrkParams::SetStepMin();   
+    TrkParams::SetMiniDefault();
     TrkParams::SetPFA();
 
 };
@@ -621,9 +623,10 @@
 //     path_.error   = 0;
 //     readb_();
 
-    TrkParams::SetTrackingMode();
-    TrkParams::SetPrecisionFactor();
-    TrkParams::SetStepMin();
+//     TrkParams::SetTrackingMode();
+//     TrkParams::SetPrecisionFactor();
+//     TrkParams::SetStepMin();
+    TrkParams::SetMiniDefault();
 
     TrkParams::Set(path,1);
     TrkParams::Load(1);
@@ -724,7 +727,7 @@
  */
 Bool_t TrkTrack::EvaluateClusterPositions(){
     
-//     cout << "void TrkTrack::GetClusterPositions() "<<endl;
+//     cout << "void TrkTrack::GetClusterositions() "<<endl;
 
     TrkParams::Load( );
     if( !TrkParams::IsLoaded() )return false;
@@ -772,10 +775,11 @@
  * @see EvaluateClusterPositions()
  *
  * The fitting procedure can be varied by changing the tracking mode, 
- * the fit-precision factor and the minimum number of step.
+ * the fit-precision factor, the minimum number of step, etc.
  * @see SetTrackingMode(int)
  * @see SetPrecisionFactor(double)
  * @see SetStepMin(int)
+ * @see SetDeltaB(int,double)
  */
 void TrkTrack::Fit(double pfixed, int& fail, int iprint, int froml1){
 
@@ -878,6 +882,14 @@
     track_.istepmin = istepmin;
 }
 /**
+ * Set deltaB parameters (id=0,1). By default they are set to zero.
+ */
+void TrkTrack::SetDeltaB(int id, double db){
+    if(id!=0 && id!=1)cout << "void TrkTrack::SetDeltaB(int id,double db) -- wrong input parameters: "<<id<<" "<<db<<endl;
+    TrkParams::SetDeltaB(id,db);
+}
+
+/**
  * Returns 1 if the track is inside the magnet cavity
  * Set the minimum number of steps for tracking precision 
  */
@@ -1473,8 +1485,8 @@
 TrkTrack *TrkLevel2::GetStoredTrack(int is){
 
     if(is >= this->ntrk()){
-	cout << "** TrkLevel2 ** Track "<< is << "doen not exits! " << endl; 
-	cout << "                Stored tracks ntrk() = "<< this->ntrk() << endl;
+	cout << "TrkTrack *TrkLevel2::GetStoredTrack(int) >> Track "<< is << "doen not exits! " << endl; 
+	cout << "Stored tracks ntrk() = "<< this->ntrk() << endl;
 	return 0;
     }
     if(!Track){
@@ -1495,8 +1507,8 @@
 TrkSinglet *TrkLevel2::GetSingletX(int is){
 
 	if(is >= this->nclsx()){
-		cout << "** TrkLevel2 ** Singlet "<< is << "doen not exits! " << endl; 
-		cout << "                Stored x-singlets nclsx() = "<< this->nclsx() << endl;
+		cout << "TrkSinglet *TrkLevel2::GetSingletX(int) >> Singlet "<< is << "doen not exits! " << endl; 
+		cout << "Stored x-singlets nclsx() = "<< this->nclsx() << endl;
 		return 0;
 	}
 	if(!SingletX)return 0;
@@ -1515,8 +1527,8 @@
 TrkSinglet *TrkLevel2::GetSingletY(int is){
 
 	if(is >= this->nclsy()){
-		cout << "** TrkLevel2 ** Singlet "<< is << "doen not exits! " << endl; 
-		cout << "                Stored y-singlets nclsy() = "<< this->nclsx() << endl;
+		cout << "TrkSinglet *TrkLevel2::GetSingletY(int) >> Singlet "<< is << "doen not exits! " << endl; 
+		cout << "Stored y-singlets nclsx() = "<< this->nclsx() << endl;
 		return 0;
 	}
 	if(!SingletY)return 0;
@@ -1536,8 +1548,8 @@
 TrkTrack *TrkLevel2::GetTrack(int it){
     
 	if(it >= this->GetNTracks()){
-		cout << "** TrkLevel2 ** Track "<< it << "does not exits! " << endl; 
-		cout << "                Physical tracks GetNTracks() = "<< this->ntrk() << endl;
+		cout << "TrkTrack *TrkLevel2::GetTrack(int) >> Track "<< it << "does not exits! " << endl; 
+		cout << "Physical tracks GetNTracks() = "<< this->ntrk() << endl;
 		return 0;
 	}
 	
@@ -1574,8 +1586,8 @@
 TrkTrack *TrkLevel2::GetTrackImage(int it){
 
     if(it >= this->GetNTracks()){
-	cout << "** TrkLevel2 ** Track "<< it << "does not exits! " << endl; 
-	cout << "                Physical tracks GetNTracks() = "<< this->ntrk() << endl;
+	cout << "TrkTrack *TrkLevel2::GetTrackImage(int) >> Track "<< it << "does not exits! " << endl; 
+	cout << "Physical tracks GetNTracks() = "<< this->ntrk() << endl;
 	return 0;
     }
 	
@@ -1584,7 +1596,7 @@
     TrkTrack *track = (TrkTrack*)sorted->At(it);
 	
     if(!track->HasImage()){
-	cout << "** TrkLevel2 ** Track "<< it << "does not have image! " << endl; 
+	cout << "TrkTrack *TrkLevel2::GetTrackImage(int) >> Track "<< it << "does not have image! " << endl; 
 	return 0;
     }
     if(!Track)return 0;
@@ -1611,9 +1623,10 @@
 //     path_.error   = 0;
 //     readb_();
 
-    TrkParams::SetTrackingMode();
-    TrkParams::SetPrecisionFactor();
-    TrkParams::SetStepMin();
+//     TrkParams::SetTrackingMode();
+//     TrkParams::SetPrecisionFactor();
+//     TrkParams::SetStepMin();
+    TrkParams::SetMiniDefault();
 
     TrkParams::Set(path,1);
     TrkParams::Load(1);