Turn off gradient bfield testing for now as the MC is broken.

KFTrack · Mar 14, 2024 · b6e53d9 · b6e53d9
1 parent 157aefe
commit b6e53d9
Show file tree

Hide file tree

Showing 3 changed files with 17 additions and 18 deletions.
diff --git a/Fit/DomainWall.hh b/Fit/DomainWall.hh
@@ -61,15 +61,14 @@ namespace KinKal {
     if(tdir == TimeDir::forwards) {
       prevwt_ += fstate.wData();
       fstate.append(dpfwd_,tdir);
-      // fstate.pData().covariance() = ROOT::Math::Similarity(dpdpdb_,fstate.pData().covariance());
+      // fstate.pData().covariance() = ROOT::Math::Similarity(dpdpdb_,fstate.pData().covariance());  Not tested TODO
       nextwt_ += fstate.wData();
     } else {
       nextwt_ += fstate.wData();
       fstate.append(dpfwd_,tdir);
       // fstate.pData().covariance() = ROOT::Math::SimilarityT(dpdpdb_,fstate.pData().covariance());
       prevwt_ += fstate.wData();
     }
-    // rotate the covariance matrix for the change in reference BField.  TODO
   }
 
   template<class KTRAJ> void DomainWall<KTRAJ>::updateState(MetaIterConfig const& miconfig,bool first) {
@@ -96,13 +95,11 @@ namespace KinKal {
     if( tdir == TimeDir::forwards){
       newpiece.range() = TimeRange(next_->begin(),std::max(ptraj.range().end(),next_->end()));
       newpiece.params() = Parameters(nextwt_);
-//      newpiece.setBNom(time(),next_->bnom());
       newpiece.resetBNom(next_->bnom());
       ptraj.append(newpiece);
     } else {
       newpiece.range() = TimeRange(std::min(ptraj.range().begin(),prev_->begin()),prev_->end());
       newpiece.params() = Parameters(prevwt_);
-//      newpiece.setBNom(time(),prev_->bnom());
       newpiece.resetBNom(prev_->bnom());
       ptraj.prepend(newpiece);
     }

diff --git a/Tests/CentralHelixFit_unit.cc b/Tests/CentralHelixFit_unit.cc
@@ -6,14 +6,15 @@ int main(int argc, char **argv) {
     cout << "Testing gradient field, correction 2" << endl;
     std::vector<std::string> arguments;
     arguments.push_back(argv[0]);
-    arguments.push_back("--Bgrad");
-    arguments.push_back("-0.036"); // mu2e-like field gradient
+//    arguments.push_back("--Bgrad");
+//    arguments.push_back("-0.036"); // mu2e-like field gradient
     arguments.push_back("--Schedule");
-    arguments.push_back("seedfit.txt");
-    arguments.push_back("--extend");
-    arguments.push_back("driftextend.txt");
-    arguments.push_back("--ttree");
-    arguments.push_back("1");
+    arguments.push_back("MCAmbigFixedField.txt");
+//    arguments.push_back("seedfit.txt");
+//    arguments.push_back("--extend");
+//    arguments.push_back("driftextend.txt");
+//    arguments.push_back("--ttree");
+//    arguments.push_back("1");
     std::vector<char*> myargv;
     for (const auto& arg : arguments)
       myargv.push_back((char*)arg.data());

diff --git a/Tests/LoopHelixFit_unit.cc b/Tests/LoopHelixFit_unit.cc
@@ -6,14 +6,15 @@ int main(int argc, char **argv) {
     cout << "Testing gradient field, correction 2" << endl;
     std::vector<std::string> arguments;
     arguments.push_back(argv[0]);
-    arguments.push_back("--Bgrad");
-    arguments.push_back("-0.036"); // mu2e-like field gradient
+//    arguments.push_back("--Bgrad");
+//    arguments.push_back("-0.036"); // mu2e-like field gradient
     arguments.push_back("--Schedule");
-    arguments.push_back("seedfit.txt");
-    arguments.push_back("--extend");
-    arguments.push_back("driftextend.txt");
-    arguments.push_back("--ttree");
-    arguments.push_back("1");
+    arguments.push_back("MCAmbigFixedField.txt");
+//    arguments.push_back("seedfit.txt");
+//    arguments.push_back("--extend");
+//    arguments.push_back("driftextend.txt");
+//    arguments.push_back("--ttree");
+//    arguments.push_back("1");
     std::vector<char*> myargv;
     for (const auto& arg : arguments)
       myargv.push_back((char*)arg.data());