CaPolynomialSearchWindow.cxx

Go to the documentation of this file.

/* Copyright (C) 2022 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergey Gorbunov, Sergei Zharko [committer] */
 
#include "CaPolynomialSearchWindow.h"
 
#include <cassert>
#include <iomanip>
#include <sstream>
 
using cbm::algo::ca::PolynomialSearchWindow;
 
// ---------------------------------------------------------------------------------------------------------------------
//
PolynomialSearchWindow::PolynomialSearchWindow(int stationID, int trackGrID)
  : fStationID(stationID)
  , fTrackGroupID(trackGrID)
{
  assert(stationID > -1);
  assert(trackGrID > -1);
 
  // Case for constant windows (TEMPORARY: we should add selection of different windows):
  static_assert(kNpars == 1);
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDxMaxVsX0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDxMaxVsX0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDxMinVsX0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDxMinVsX0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDxMaxVsY0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDxMaxVsY0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDxMinVsY0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDxMinVsY0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDyMaxVsX0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDyMaxVsX0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDyMinVsX0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDyMinVsX0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDyMaxVsY0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDyMaxVsY0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
void PolynomialSearchWindow::SetParamDyMinVsY0(int id, float val)
{
  assert(id > -1 && id < kNpars);
  fvParams[kDyMinVsY0 * kNpars + id] = val;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
std::string PolynomialSearchWindow::ToString() const
{
  using std::setw;
  std::stringstream msg;
  msg << "----- CA hits search window: \n";
  msg << "\tstation ID:      " << fStationID << '\n';
  msg << "\ttracks group ID: " << fTrackGroupID << '\n';
  msg << "\tparameters:\n";
  msg << "\t\t" << setw(6) << "No." << ' ';
  msg << setw(12) << "dx_max(x0)" << ' ';
  msg << setw(12) << "dx_min(x0)" << ' ';
  msg << setw(12) << "dx_max(y0)" << ' ';
  msg << setw(12) << "dx_min(y0)" << ' ';
  msg << setw(12) << "dy_max(x0)" << ' ';
  msg << setw(12) << "dy_min(x0)" << ' ';
  msg << setw(12) << "dy_max(y0)" << ' ';
  msg << setw(12) << "dy_min(y0)" << '\n';
  for (int iPar = 0; iPar < kNpars; ++iPar) {
    msg << "\t\t" << setw(6) << iPar << ' ';
    for (int iDep = 0; iDep < kNdeps; ++iDep) {
      msg << setw(12) << fvParams[iDep * kNpars + iPar] << ' ';
    }
  }
  return msg.str();
}