NoKickCuts.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Valerio Bertacchi                                        *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#include <TObject.h>
 
#pragma once
 
 
namespace Belle2 {
  class NoKickCuts: public TObject {
 
  public:
    NoKickCuts() {}
 
    explicit NoKickCuts(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 cointained
    *  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(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);
  };
}