NoKickCutsEvalModule.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.                *
 **************************************************************************/
 
#ifndef NOKICKCUTSEVAL_H
#define NOKICKCUTSEVAL_H
 
#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;
 
    void terminate() override;
 
    // enum EParameters {
    //   omega,
    //   d0,
    //   phi0,
    //   z0,
    //   tanlambda
    // };
 
    double deltaParEval(hitXP hit1, hitXP hit2, NoKickCuts::EParameters par, bool is0 = false);
 
 
    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]",
      ""
    };
  };
}
#endif