PVertexerHelpers.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#include "DetectorsVertexing/PVertexerHelpers.h"
 
using namespace o2::vertexing;
 
void VertexSeed::print() const
{
  auto terr2 = tMeanAccErr > 0 ? 1. / tMeanAccErr : 0.;
  LOGF(info, "VtxSeed: Scale: %+e ScalePrev: %+e |NScaleIncreased: %d NSlowDecrease: %d| WChi2: %e WSum: %e | TMean: %e TMeanE: %e\n",
       scaleSigma2, scaleSigma2Prev, nScaleIncrease, nScaleSlowConvergence, wghChi2, wghSum, tMeanAcc * terr2, std::sqrt(terr2));
  double dZP, rmsZP, dZN, rmsZN, dTP, rmsTP, dTN, rmsTN;
  double dZ, rmsZ, dT, rmsT;
  PVertex::print();
}
 
int SeedHistoTZ::findPeakBin()
{
  if (nEntries < 2) {
    return -1;
  }
  int maxBin = -1, ib = filledBins.size(), last = ib;
  float maxv = 0.;
  while (ib--) {
    auto bin = filledBins[ib];
    auto v = getBinContent(bin);
    if (v > maxv) {
      maxv = v;
      maxBin = bin;
    } else if (v <= 0.) {                  // bin was emptied
      filledBins[ib] = filledBins[--last]; // move last non-empty bin in place of emptied one
    }
  }
  filledBins.resize(last);
  return maxBin;
}
 
void TrackVF::reportBadTrack(const o2::track::TrackParCov& src, const TimeEst& t_est, GTrackID _gid)
{
  constexpr int MaxRep = 10;
  static int count = 0;
  if (count < MaxRep) {
    std::string msg = fmt::format("Will reject bad track {} {} at t={}", _gid.asString(), src.asString(), t_est.getTimeStamp());
    if (++count == MaxRep) {
      msg += " | further warnings will be suppressed";
    }
    LOG(warn) << msg;
  }
}