CaDoubletSearchWindowMap.h

Go to the documentation of this file.

/* Copyright (C) 2025 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergey Gorbunov [committer] */

 
#pragma once  // include this header only once per compilation unit
 
#include "CaSearchWindowMap.h"
 
#include <boost/serialization/array.hpp>
#include <boost/serialization/string.hpp>
 
#include <cmath>
 
namespace cbm::algo::ca
{
  class DoubletSearchWindowMap {
   public:
    DoubletSearchWindowMap() = default;

    ~DoubletSearchWindowMap() = default;

    void Init(float targetX, float targetY, float targetZ, float refZ, float xMin, float xMax, int nBinsX, float yMin,
              float yMax, int nBinsY)
    {
      fTargetX = targetX;
      fTargetY = targetY;
      fTargetZ = targetZ;
      fRefZ    = refZ;
      fMap.SetRange(xMin, xMax, nBinsX, yMin, yMax, nBinsY);
    }

    SearchWindowMap::SearchWindow GetSearchWindow(float x, float y, float z, float dx, float dy) const
    {
      // hit :  x, y, z
      // track : straight line at (x, y, z)
      // tx = tx = (x - fTargetX) / (z - fTargetZ)
      // ty = (y - fTargetY) / (z - fTargetZ)
      // intersection with the plane at zRef
      // c     = (fRefZ - z) / (z - fTargetZ);
      // xr = x + tx * (zRef-z) = x * (1. + c) - fTargetX * c
      // yr = y + ty * (zRef-z) = y * (1. + c) - fTargetY * c
      // dxr^2 = dx^2 * (1. + c)^2
      // dyr^2 = dy^2 * (1. + c)^2
 
      float c     = (fRefZ - z) / (z - fTargetZ);
      float cp1   = (1.f + c);
      float xe    = x + (x - fTargetX) * c;
      float ye    = y + (y - fTargetY) * c;
      float dxe   = std::fabs(dx * cp1);
      float dye   = std::fabs(dy * cp1);
      auto window = fMap.GetSearchWindow(xe, ye);
      window.xMin -= dxe;
      window.xMax += dxe;
      window.yMin -= dye;
      window.yMax += dye;
      return window;
    };
 
    SearchWindowMap& GetMap() { return fMap; }

    std::string ToString() const;
 
    // TODO: make the members private
   public:
    SearchWindowMap fMap;  
    float fTargetX{0.f};   
    float fTargetY{0.f};   
    float fTargetZ{0.f};   
    float fRefZ{0.f};      

    friend class boost::serialization::access;
    template<class Archive>
    void serialize(Archive& ar, const unsigned int)
    {
      ar& fMap;
      ar& fTargetX;
      ar& fTargetY;
      ar& fTargetZ;
      ar& fRefZ;
    }
  };
}  // namespace cbm::algo::ca