Acceptance.cpp

Go to the documentation of this file.

/*
 * Thermal-FIST package
 * 
 * Copyright (c) 2015-2018 Volodymyr Vovchenko
 *
 * GNU General Public License (GPLv3 or later)
 */
#include "HRGEventGenerator/Acceptance.h"

#include <fstream>

namespace thermalfist {

  int Acceptance::ReadAcceptanceFunction(Acceptance::AcceptanceFunction & func, std::string filename)
  {
    double ymin = 0., ymax = 6.;
    double ptmin = 0., ptmax = 2.5;
    func.ys.resize(0);
    func.pts.resize(0);
    func.probs.resize(0);
    std::ifstream fin(filename.c_str());
    if (!fin.is_open()) return 0;
    fin >> func.dy >> func.dpt;
    double ty, tpt, prob;
    func.ys.resize(0);
    func.pts.resize(0);
    func.probs.resize(0);
    while (fin >> ty >> tpt >> prob) {
      if (tpt<ptmin || tpt>ptmax || ty<ymin || ty>ymax) continue;
      func.ys.push_back(ty);
      func.pts.push_back(tpt);
      func.probs.push_back(prob);
    }
    func.setSpline();
    fin.close();
    return 1;
  }

  double Acceptance::AcceptanceFunction::getAcceptance(const double & y, const double & pt) const {
    double ret = sfunc.Eval(y, pt);
    if (ret < 0.) ret = 0.;
    if (ret > 1.) ret = 1.;
    return ret;
  }

} // namespace thermalfist