KF: add setter for target position

external/InstallKFParticle.cmake

@@ -6,7 +6,7 @@ set(KFPARTICLE_DESTDIR "${CMAKE_BINARY_DIR}/external/KFPARTICLE-prefix")
 # GIT_TAG is a hash for KFParticle tag cbm/v1.1-1 
 download_project_if_needed(PROJECT         kfparticle_source
                            GIT_REPOSITORY  ${KFPARTICLE_SRC_URL}
-                           GIT_TAG         "1eddb796b01d91ac7faf020e3b616556c28da4d0"
+                           GIT_TAG         "1e1c5dee0e2bd98a3df677beb88e497fbef92504"
                            SOURCE_DIR      ${CMAKE_CURRENT_SOURCE_DIR}/KFParticle
                            TEST_FILE       CMakeLists.txt
                           )

reco/KF/CbmKF.cxx

@@ -278,6 +278,10 @@ void CbmKF::GetTargetInfo()
   target.y = global[1];
   target.z = global[2];
 
+  fTargetXYZ[0] = target.x;
+  fTargetXYZ[1] = target.y;
+  fTargetXYZ[2] = target.z;
+
   if (fVerbose) {
     cout << "KALMAN FILTER : === READ TARGET MATERIAL ===" << endl;
     cout << " found targed \"" << targetPath << "\" at ( " << target.x << " " << target.y << " " << target.z << " ) "
@@ -310,6 +314,8 @@ void CbmKF::GetTargetInfo()
   LOG(info) << "Target info: " << target.KFInfo();
 }
 
+std::array<double, 3> CbmKF::GetTargetPosition() { return fTargetXYZ; }
+
 void CbmKF::FindTargetNode(TString& targetPath, TGeoNode*& targetNode)
 {
   if (!targetNode) {  // init at the top of the tree

reco/KF/CbmKF.h

@@ -91,6 +91,8 @@ public:
 
   int GetNMvdStations() const { return CbmKF::Instance()->vMvdMaterial.size(); }
 
+  std::array<double, 3> GetTargetPosition();
+
 private:
   static CbmKF* fInstance;
 
@@ -108,8 +110,12 @@ private:
   CbmKFMaterial* ReadPassive(FairGeoNode* node);
   void GetTargetInfo();
 
+
   void FindTargetNode(TString& targetPath, TGeoNode*& targetNode);
 
+  std::array<double, 3> fTargetXYZ;
+
+
 private:
   CbmKF(const CbmKF&);
   void operator=(const CbmKF&);

reco/KF/CbmKFParticleFinder.cxx

@@ -579,6 +579,9 @@ void CbmKFParticleFinder::SetSuperEventAnalysis()
   fTopoReconstructor->SetMixedEventAnalysis();
 }
 
+
+void CbmKFParticleFinder::SetTarget(const std::array<float, 3>& target) { fTopoReconstructor->SetTarget(target); }
+
 KFParticleFinder* CbmKFParticleFinder::GetKFParticleFinder() { return fTopoReconstructor->GetKFParticleFinder(); }
 void CbmKFParticleFinder::SetMaxDistanceBetweenParticlesCut(float cut)
 {

reco/KF/CbmKFParticleFinder.h

@@ -50,12 +50,15 @@ public:
 
   void SetPIDInformation(CbmKFParticleFinderPID* pid) { fPID = pid; }
 
+
   // set cuts
   void SetPrimaryProbCut(float prob);
 
   // Set SE analysis
   void SetSuperEventAnalysis();
 
+  void SetTarget(const std::array<float, 3>& target);
+
   //KF Particle Finder cuts
   void SetMaxDistanceBetweenParticlesCut(float cut);
   void SetLCut(float cut);