From e17451ca7a022b98f21d17d60e977c9b6e2fc27d Mon Sep 17 00:00:00 2001
From: mmusich <marco.musich@cern.ch>
Date: Wed, 13 Mar 2024 10:31:59 +0100
Subject: [PATCH] fix OfflinePixel3DPrimaryVertices after PR #43846

---
 .../OfflinePixel3DPrimaryVertices_cfi.py      | 26 ++++++++++++++-----
 1 file changed, 20 insertions(+), 6 deletions(-)

diff --git a/RecoVertex/PrimaryVertexProducer/python/OfflinePixel3DPrimaryVertices_cfi.py b/RecoVertex/PrimaryVertexProducer/python/OfflinePixel3DPrimaryVertices_cfi.py
index b6ad7e2ed8eec..437d77bb286ef 100644
--- a/RecoVertex/PrimaryVertexProducer/python/OfflinePixel3DPrimaryVertices_cfi.py
+++ b/RecoVertex/PrimaryVertexProducer/python/OfflinePixel3DPrimaryVertices_cfi.py
@@ -12,12 +12,26 @@
         maxEta = 100.,
     ),
 
-    TkClusParameters = dict(
-        TkDAClusParameters = dict(
-            dzCutOff = 4.0,
-            Tmin = 2.4,
-            vertexSize = 0.01,
-            uniquetrkweight = 0.9
+    TkClusParameters = cms.PSet(
+        algorithm   = cms.string("DA_vect"),
+        TkDAClusParameters = cms.PSet(
+            coolingFactor = cms.double(0.6),  # moderate annealing speed
+            zrange = cms.double(4.),          # consider only clusters within 4 sigma*sqrt(T) of a track
+            delta_highT = cms.double(1.e-2),  # convergence requirement at high T
+            delta_lowT = cms.double(1.e-3),   # convergence requirement at low T
+            convergence_mode = cms.int32(0),  # 0 = two steps, 1 = dynamic with sqrt(T)
+            Tmin = cms.double(2.4),           # end of vertex splitting, (adjusted for Pixel vertices)
+            Tpurge = cms.double(2.0),         # cleaning
+            Tstop = cms.double(0.5),          # end of annealing
+            vertexSize = cms.double(0.01),    # added in quadrature to track-z resolution (adjusted for Pixel vertices)
+            d0CutOff = cms.double(4.),        # downweight high IP tracks (adjusted for Pixel vertices)
+            dzCutOff = cms.double(3.),        # outlier rejection after freeze-out (T<Tmin)
+            zmerge = cms.double(1e-2),        # merge intermediat clusters separated by less than zmerge
+            uniquetrkweight = cms.double(0.9),# require at least two tracks with this weight at T=Tpurge
+            uniquetrkminp = cms.double(0.0),  # minimal a priori track weight for counting unique tracks
+            runInBlocks = cms.bool(False),    # activate the DA running in blocks of z sorted tracks
+            block_size = cms.uint32(10000),   # block size in tracks
+            overlap_frac = cms.double(0.0)    # overlap between consecutive blocks (blocks_size*overlap_frac)
         )
     ),