SVDClusterAbsoluteTimeShifterAlgorithm.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 <calibration/CalibrationAlgorithm.h>
 
#include <TMath.h>
 
#include <optional>
 
namespace Belle2 {
  class SVDClusterAbsoluteTimeShifterAlgorithm : public CalibrationAlgorithm {
  public:

    explicit SVDClusterAbsoluteTimeShifterAlgorithm(const std::string& str);

    virtual ~SVDClusterAbsoluteTimeShifterAlgorithm() override {}

    void setAllowedTimeShift(float value) {m_allowedTimeShift = value;}

    float getAllowedTimeShift() {return m_allowedTimeShift;}

    void setMaximumAllowedShift(const float& value) {m_maximumAllowedShift = value;}

    float getMaximumAllowedShift() {return m_maximumAllowedShift;}

    void setTimeAlgorithm(const std::vector<TString>& lst) {m_timeAlgorithms = lst;}

    std::vector<TString> getTimeAlgorithm() {return m_timeAlgorithms;}

    void setMinEntries(const int& minEntries) {m_minEntries = minEntries;}

    int getMinEntries() {return m_minEntries;}
 
  protected:

    virtual EResult calibrate() override;

    virtual bool isBoundaryRequired(const Calibration::ExpRun& currentRun) override;

    virtual void boundaryFindingSetup(std::vector<Calibration::ExpRun> /*runs*/, int /*iteration = 0*/) override
    {
      m_previousTimeMeanL3V.reset();
    }
 
 
  private:
 
    std::string m_id = ""; 
    std::optional<float> m_previousTimeMeanL3V; 
    float m_allowedTimeShift = 2.; 
    float m_maximumAllowedShift = 15.; 
    int m_minEntries = 100; 
    std::vector<TString> m_timeAlgorithms = {"CoG3", "ELS3", "CoG6"}; 
    int m_innerLayer = 3; 
    int m_outerLayer = 6; 
  };

  inline double AbssingleGaus(const double* x, const double* par)
  {
    double N = std::fabs(par[0]);
    double a = par[1];
    double b = par[2];
    double e = std::fabs(par[3]);
    return N * TMath::Gaus(x[0], a, b, true) + e;
  };

  inline double AbsdoubleGaus(const double* x, const double* par)
  {
    double N = std::fabs(par[0]);
    double f = std::fabs(par[1]);
    double a = par[2];
    double b = par[3];
    double c = par[4];
    double d = par[5];
    double e = std::fabs(par[6]);
    return (N * (f * TMath::Gaus(x[0], a, b) + (1 - f) * TMath::Gaus(x[0], c, d)) /
            (TMath::Sqrt(2. * TMath::Pi()) * ((b - d) * f + d)) + e);
  };

} // namespace Belle2