NoKickCutsEvalModule.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/core/Module.h>
#include <tracking/dataobjects/hitXP.h>
#include <tracking/trackFindingVXD/sectorMapTools/NoKickCuts.h>
#include <tracking/trackFindingVXD/sectorMapTools/NoKickRTSel.h>
#include <TFile.h>
#include <math.h>
 
namespace Belle2 {

  class NoKickCutsEvalModule : public Module {
  public:
    NoKickCutsEvalModule();

    void initialize() override;

    void event() override;

    void endRun() override;
 

    // enum EParameters {
    //   omega,
    //   d0,
    //   phi0,
    //   z0,
    //   tanlambda
    // };

    double deltaParEval(const hitXP& hit1, const hitXP& hit2, NoKickCuts::EParameters par, bool is0 = false);
 

    static double cutFunction(int p, double pwidth);

    // double cutFunction(int p, double pwid)
    // {
    //   double out;
    //   double mom = p * pwid + c_pmin;
    //   out = -3.971 * pow(10, -7) / pow(mom, 3.373) + 1;
    //   return out;
    // }
 
  private:
 
    const double c_pmin = 0.025; 
    const double c_pmax = 10.; 
    const double c_tmin = 17.*M_PI / 180.; 
    const double c_tmax = 5. / 6.*M_PI; 
    const int c_nbin = 5000; 
    const int c_nbinp = 200; 
    const int c_nbinpar = 5; 
    const int c_nbinlay = 7;
    const int c_nbint = 3; 
    double c_pwidth = (c_pmax - c_pmin) / (double)c_nbinp; 
    double c_twidth = (c_tmax - c_tmin) / (double)c_nbint; 
    const double c_multLimit = 1; 
    const double c_over = NAN; 
    int m_pCounter = 0; 
    int m_tCounter = 0; 
    int m_globCounter = 0; 
    bool c_validationON; 
    bool c_fitMethod; 
 
    NoKickRTSel m_trackSel; 
    TFile* m_outputFile = nullptr; 
    std::vector<double> m_histoLim; 
    std::vector<std::vector<std::vector<std::vector<std::vector<TH1F*>>>>> m_histo; 

    std::vector<TString> m_namePar = {
      "#omega",
      "d0",
      "#phi0",
      "z0",
      "tan#lambda"
    };

    std::vector<TString> m_unitPar = {
      "[cm^{-1}]",
      "[cm]",
      "[rad]",
      "[cm]",
      ""
    };
  };

}