NoKickCuts.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.                  *
 **************************************************************************/
 
#include <TObject.h>
 
#pragma once
 
 
namespace Belle2 {

 
 
  class NoKickCuts: public TObject {
 
  public:
    NoKickCuts() {}

    explicit NoKickCuts(const std::string& m_fileName)
    {
      if (m_fileName.size() != 0) FillCuts(m_fileName); //else cuts not applied
    }

    enum EMinMax {
      c_Min,
      c_Max
    };

    enum EParameters {
      c_Omega,
      c_D0,
      c_Phi0,
      c_Z0,
      c_Tanlambda
    };

    enum ECutName {
      c_Norm,
      c_Pow,
      c_Bkg
    };

    std::vector<double> cutSelector(double sintheta, double momentum, int layer1, int layer2, EParameters par);
 
    /*  This methods select the cuts as function parameters from TH3F contained
    *  in cutfile in function of track parameter and layer only. EMinMax represent
    *  if you want the maximum or minimum cut, instead ECutName are the three
    *  parameters of the theta-p parametrization of the cuts (see the enum type)
    * input: (first layer, second layer, track parameter, min/max value, fit parameter)
    * output: (value of the fit parameter of a specific cut)
    */
    double getCut(int layer1, int layer2, EParameters par, EMinMax m, ECutName cut);

    void FillCuts(const std::string& m_fileName);
 
  private:
    std::vector<std::vector <std::vector<std::vector <double>>>> m_cutNorm; 
    std::vector<std::vector <std::vector<std::vector <double>>>> m_cutPow; 
    std::vector<std::vector <std::vector<std::vector <double>>>> m_cutBkg; 
 
    int m_nbinpar = 5; 
    int m_nbinlay = 7; 

    ClassDef(NoKickCuts, 1);
  };

}