SVDClusterTimeShifter.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <framework/logging/Logger.h>
 
#include "TObject.h"
#include "TString.h"
#include <map>
 
namespace Belle2 {
  class SVDClusterTimeShifter: public TObject {
  public:
    SVDClusterTimeShifter(const TString& uniqueID = "")
      : m_uniqueID(uniqueID)
    {
      m_svdClusterTimeShift.clear();
    };
 
    Double_t getClusterTimeShift(const TString& alg,
                                 const int& layer, const int& sensor,
                                 const bool& isU, const int& size) const
    {
      if (auto searchAlg = m_svdClusterTimeShift.find(alg); // search for time alg
          searchAlg != m_svdClusterTimeShift.end()) {
        TString sensorType = getSensorType(layer, sensor, isU);
        if (auto searchShift = (searchAlg->second).find(sensorType); // search for shift values
            searchShift != (searchAlg->second).end()) {
          int maxClusters = (searchShift->second).size();
          if (maxClusters == 0)
            return 0.;
          else if (size > maxClusters)
            return (searchShift->second)[maxClusters - 1];
          else
            return (searchShift->second)[size - 1];
        }
      }
      return 0.;    // returns zero if map is empty
    };
 
    void setClusterTimeShift(const TString& alg, const TString& sensorType,
                             const std::vector<Double_t>& shiftValues)
    {
      B2INFO("Shift values for " << alg << " and " << sensorType << " is set.");
      m_svdClusterTimeShift[alg][sensorType] = shiftValues;
    };
 
    TString get_uniqueID() const {return m_uniqueID;}
 
    static std::string name;
 
  protected:
 
    TString getSensorType(const int& layer, const int& sensor, const bool& isU) const
    {
      return TString::Format("L%iS%iS%c", layer, sensor, (isU ? 'U' : 'V'));
    };
 
  private:
 
    TString m_uniqueID;
 
    std::map<TString, std::map<TString, std::vector<Double_t>>> m_svdClusterTimeShift;
 
    ClassDef(SVDClusterTimeShifter, 1); 
  };
 
}