modified: ../src/components/TrackParamClusterInit.cpp

	modified:   ../src/components/TrackParamTruthInit.cpp
	modified:   ../src/components/TrackParamVertexClusterInit.cpp
	modified:   ../src/components/TrackerSourcesLinker.cpp

JugReco/src/components/TrackParamClusterInit.cpp

@@ -105,8 +105,14 @@ namespace Jug::Reco {
         // add both charges to the track candidate...
         init_trk_params->push_back({pSurface, params,  1});
 
-        params(Acts::eBoundQOverP) = -1/p;
-        init_trk_params->push_back({pSurface, params, -1});
+        Acts::BoundVector  params2;
+        params2(Acts::eBoundLoc0)   = 0.0 * mm ;
+        params2(Acts::eBoundLoc1)   = 0.0 * mm ;
+        params2(Acts::eBoundPhi)    = momentum.Phi();
+        params2(Acts::eBoundTheta)  = momentum.Theta();
+        params2(Acts::eBoundQOverP) = -1/p;
+        params2(Acts::eBoundTime)   = 0 * ns;
+        init_trk_params->push_back({pSurface, params2, -1});
 
         // acts v1.2.0:
         //init_trk_params->emplace_back(Acts::Vector4(0 * mm, 0 * mm, 0 * mm, 0),

JugReco/src/components/TrackParamTruthInit.cpp

@@ -108,6 +108,7 @@ namespace Jug::Reco {
         auto pSurface = Acts::Surface::makeShared<Acts::PerigeeSurface>(
             Acts::Vector3{part.vsx() * mm, part.vsy() * mm, part.vsz() * mm});
 
+        params(Acts::eBoundQOverP) = charge/p;
         init_trk_params->push_back({pSurface, params, charge});
         // std::make_optional(std::move(cov))
 

JugReco/src/components/TrackParamVertexClusterInit.cpp

@@ -121,8 +121,14 @@ namespace Jug::Reco {
           // add both charges to the track candidate...
           init_trk_params->push_back({pSurface, params, 1});
 
-          params(Acts::eBoundQOverP) = -1 / p_cluster;
-          init_trk_params->push_back({pSurface, params, -1});
+          Acts::BoundVector params2;
+          params2(Acts::eBoundLoc0)   = 0.0 * mm;
+          params2(Acts::eBoundLoc1)   = 0.0 * mm;
+          params2(Acts::eBoundPhi)    = momentum.Phi();
+          params2(Acts::eBoundTheta)  = momentum.Theta();
+          params2(Acts::eBoundQOverP) = -1 / p_cluster;
+          params2(Acts::eBoundTime)   = 0 * ns;
+          init_trk_params->push_back({pSurface, params2, -1});
         }
         // init_trk_params->emplace_back(Acts::Vector4(0 * mm, 0 * mm, 0 * mm, 0),
         //                              Acts::Vector3(c.x() * p / len, c.y() * p / len, c.z() * p / len), p, 1,

JugReco/src/components/TrackerSourcesLinker.cpp

@@ -124,7 +124,7 @@ namespace Jug::Reco {
         debug() << (*hits).size() << " hits " << endmsg;
         for (const auto& ahit : *hits) {
 
-          allHits->push_back(ahit);
+          allHits->push_back(ahit.clone());
 
           Acts::SymMatrix2 cov = Acts::SymMatrix2::Zero();
           cov(0, 0) = ahit.covsym_xx() * Acts::UnitConstants::mm * ahit.covsym_xx() * Acts::UnitConstants::mm;