Misc. fixes for the pvertexer.

Author: shahor02

Detectors/Vertexing/include/DetectorsVertexing/PVertexer.h

@@ -37,7 +37,7 @@
 #include <TFile.h>
 #include <TStopwatch.h>
 
-// #define _PV_DEBUG_TREE_ // if enabled, produce dbscan and vertex comparison dump
+#define _PV_DEBUG_TREE_ // if enabled, produce dbscan and vertex comparison dump
 
 namespace o2
 {
@@ -263,7 +263,6 @@ void PVertexer::createTracksPool(const TR& tracks, gsl::span<const o2d::GlobalTr
 
   // check all containers
   float vtxErr2 = 0.5 * (mMeanVertex.getSigmaX2() + mMeanVertex.getSigmaY2()) + mPVParams->meanVertexExtraErrSelection * mPVParams->meanVertexExtraErrSelection;
-  o2d::DCA dca;
 
   for (uint32_t i = 0; i < ntGlo; i++) {
     int id = sortedTrackID[i];

Detectors/Vertexing/include/DetectorsVertexing/PVertexerParams.h

@@ -40,7 +40,7 @@ struct PVertexerParams : public o2::conf::ConfigurableParamHelper<PVertexerParam
   long maxTimeMSPerCluster = 10000; ///< max allowed time per TZCluster processing, ms
 
   // track selection
-  float meanVertexExtraErrSelection = 0.01; ///< extra error to meanvertex sigma used when selecting tracks
+  float meanVertexExtraErrSelection = 0.02; ///< extra error to meanvertex sigma used when selecting tracks
   float dcaTolerance = 1.3; ///< consider tracks within this abs DCA to mean vertex
   float pullIniCut = 9;     ///< cut on pull (n^2 sigma) on dca to mean vertex
   float maxTimeErrorMUS = 10.0; ///< max time error in ms of the track to account

Detectors/Vertexing/src/PVertexer.cxx

@@ -1051,13 +1051,14 @@ bool PVertexer::relateTrackToMeanVertex(o2::track::TrackParCov& trc, float vtxEr
 {
   o2d::DCA dca;
   auto z = trc.getZAt(0., mBz);
-  if (z < -999.) {
+  if (z < -998.) {
     z = mMeanVertex.getZ();
   }
   mMeanVertex.setMeanXYVertexAtZ(mMeanVertexSeed, z);
-  return o2::base::Propagator::Instance()->propagateToDCA(mMeanVertex, trc, mBz, 2.0f,
-                                                          o2::base::Propagator::MatCorrType::USEMatCorrLUT, &dca, nullptr, 0, mPVParams->dcaTolerance) &&
-         (dca.getY() * dca.getY() / (dca.getSigmaY2() + vtxErr2) < mPVParams->pullIniCut);
+  if (!o2::base::Propagator::Instance()->propagateToDCA(mMeanVertex, trc, mBz, 2.0f, o2::base::Propagator::MatCorrType::USEMatCorrLUT, &dca, nullptr, 0, mPVParams->dcaTolerance)) {
+    return false;
+  }
+  return dca.getY() * dca.getY() / (dca.getSigmaY2() + vtxErr2) < mPVParams->pullIniCut;
 }
 
 //______________________________________________