BTagCalibration.cc File Reference

#include "Core/BTagging/interface/BTagCalibration.h"
#include <iostream>
#include <exception>
#include <algorithm>
#include <sstream>
#include <fstream>

Include dependency graph for BTagCalibration.cc:

Classes
class	BTagCalibrationReader::BTagCalibrationReaderImpl
struct	BTagCalibrationReader::BTagCalibrationReaderImpl::TmpEntry

Functions
std::string	th1ToFormulaLin (const TH1 *hist)
std::string	th1ToFormulaBinTree (const TH1 *hist, int start=0, int end=-1)

Function Documentation

th1ToFormulaBinTree()

std::string th1ToFormulaBinTree	(	const TH1 *	hist,
		int	start = 0,
		int	end = -1
	)

                                                                        {
  if (end == -1) {                      // initialize
    start = 0.;
    end = hist->GetNbinsX()+1;
    TH1* h2 = (TH1*) hist->Clone();
    h2->SetBinContent(start, 0);  // kill underflow
    h2->SetBinContent(end, 0);    // kill overflow
    std::string res = th1ToFormulaBinTree(h2, start, end);
    delete h2;
    return res;
  }
  if (start == end) {                   // leave is reached
    char tmp_buff[20];
    sprintf(tmp_buff, "%g", hist->GetBinContent(start));
    return std::string(tmp_buff);
  }
  if (start == end - 1) {               // no parenthesis for neighbors
    char tmp_buff[70];
    sprintf(tmp_buff,
            "x<%g ? %g:%g",
            hist->GetXaxis()->GetBinUpEdge(start),
            hist->GetBinContent(start),
            hist->GetBinContent(end));
    return std::string(tmp_buff);
  }
 
  // top-down recursion
  std::stringstream buff;
  int mid = (end-start)/2 + start;
  char tmp_buff[25];
  sprintf(tmp_buff,
          "x<%g ? (",
          hist->GetXaxis()->GetBinUpEdge(mid));
  buff << tmp_buff
       << th1ToFormulaBinTree(hist, start, mid)
       << ") : ("
       << th1ToFormulaBinTree(hist, mid+1, end)
       << ")";
  return buff.str();
}

th1ToFormulaLin()

std::string th1ToFormulaLin

(

const TH1 *

hist

)

                                           {
  int nbins = hist->GetNbinsX();
  TAxis const* axis = hist->GetXaxis();
  std::stringstream buff;
  buff << "x<" << axis->GetBinLowEdge(1) << " ? 0. : ";  // default value
  for (int i=1; i<nbins+1; ++i) {
    char tmp_buff[50];
    sprintf(tmp_buff,
            "x<%g ? %g : ",  // %g is the smaller one of %e or %f
            axis->GetBinUpEdge(i),
            hist->GetBinContent(i));
    buff << tmp_buff;
  }
  buff << 0.;  // default value
  return buff.str();
}