EffPurFromHistos Class Reference

#include <EffPurFromHistos.h>

Public Types
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement

Public Member Functions
void	compute (DQMStore::IBooker &ibook)
FlavourHistograms< double > &	discriminatorCutEfficScan () const
FlavourHistograms< double > &	discriminatorNoCutEffic () const
bool	doCTagPlots (bool Ctag)
	EffPurFromHistos (const std::string &ext, TH1F *h_d, TH1F *h_u, TH1F *h_s, TH1F *h_c, TH1F *h_b, TH1F *h_g, TH1F *h_ni, TH1F *h_dus, TH1F *h_dusg, TH1F *h_pu, const std::string &label, unsigned int mc, int nBin=100, double startO=0.005, double endO=1.005)
	EffPurFromHistos (const FlavourHistograms< double > &dDiscriminatorFC, const std::string &label, unsigned int mc, DQMStore::IBooker &ibook, int nBin=100, double startO=0.005, double endO=1.005)
void	epsPlot (const std::string &name)
TH1F *	getEffFlavVsBEff_b ()
TH1F *	getEffFlavVsBEff_c ()
TH1F *	getEffFlavVsBEff_d ()
TH1F *	getEffFlavVsBEff_dus ()
TH1F *	getEffFlavVsBEff_dusg ()
TH1F *	getEffFlavVsBEff_g ()
TH1F *	getEffFlavVsBEff_ni ()
TH1F *	getEffFlavVsBEff_pu ()
TH1F *	getEffFlavVsBEff_s ()
TH1F *	getEffFlavVsBEff_u ()
void	plot (TPad *theCanvas=0)
void	plot (const std::string &name, const std::string &ext)
void	psPlot (const std::string &name)
	~EffPurFromHistos ()

Private Member Functions
void	check ()

Private Attributes
std::unique_ptr< FlavourHistograms< double > >	discrCutEfficScan
std::unique_ptr< FlavourHistograms< double > >	discrNoCutEffic
bool	doCTagPlots_
MonitorElement *	EffFlavVsXEff_b
MonitorElement *	EffFlavVsXEff_c
MonitorElement *	EffFlavVsXEff_d
MonitorElement *	EffFlavVsXEff_dus
MonitorElement *	EffFlavVsXEff_dusg
MonitorElement *	EffFlavVsXEff_g
MonitorElement *	EffFlavVsXEff_ni
MonitorElement *	EffFlavVsXEff_pu
MonitorElement *	EffFlavVsXEff_s
MonitorElement *	EffFlavVsXEff_u
TH1F *	effVersusDiscr_b
TH1F *	effVersusDiscr_c
TH1F *	effVersusDiscr_d
TH1F *	effVersusDiscr_dus
TH1F *	effVersusDiscr_dusg
TH1F *	effVersusDiscr_g
TH1F *	effVersusDiscr_ni
TH1F *	effVersusDiscr_pu
TH1F *	effVersusDiscr_s
TH1F *	effVersusDiscr_u
double	endOutput
bool	fromDiscriminatorDistr
std::string	histoExtension
std::string	label_
unsigned int	mcPlots_
int	nBinOutput
double	startOutput

Detailed Description

Definition at line 13 of file EffPurFromHistos.h.

Member Typedef Documentation

typedef dqm::legacy::DQMStore EffPurFromHistos::DQMStore

Definition at line 15 of file EffPurFromHistos.h.

typedef dqm::legacy::MonitorElement EffPurFromHistos::MonitorElement

Definition at line 16 of file EffPurFromHistos.h.

Constructor & Destructor Documentation

EffPurFromHistos::EffPurFromHistos	(	const std::string &	ext,
		TH1F *	h_d,
		TH1F *	h_u,
		TH1F *	h_s,
		TH1F *	h_c,
		TH1F *	h_b,
		TH1F *	h_g,
		TH1F *	h_ni,
		TH1F *	h_dus,
		TH1F *	h_dusg,
		TH1F *	h_pu,
		const std::string &	label,
		unsigned int	mc,
		int	nBin = 100,
		double	startO = 0.005,
		double	endO = 1.005
	)

Definition at line 16 of file EffPurFromHistos.cc.

References check().

    : fromDiscriminatorDistr(false),
      mcPlots_(mc),
      doCTagPlots_(false),
      label_(label),
      histoExtension(ext),
      effVersusDiscr_d(h_d),
      effVersusDiscr_u(h_u),
      effVersusDiscr_s(h_s),
      effVersusDiscr_c(h_c),
      effVersusDiscr_b(h_b),
      effVersusDiscr_g(h_g),
      effVersusDiscr_ni(h_ni),
      effVersusDiscr_dus(h_dus),
      effVersusDiscr_dusg(h_dusg),
      effVersusDiscr_pu(h_pu),
      nBinOutput(nBin),
      startOutput(startO),
      endOutput(endO) {
  // consistency check
  check();
}

EffPurFromHistos::EffPurFromHistos	(	const FlavourHistograms< double > &	dDiscriminatorFC,
		const std::string &	label,
		unsigned int	mc,
		DQMStore::IBooker &	ibook,
		int	nBin = 100,
		double	startO = 0.005,
		double	endO = 1.005
	)

Definition at line 54 of file EffPurFromHistos.cc.

References FlavourHistograms< T >::baseNameDescription(), FlavourHistograms< T >::baseNameTitle(), discrCutEfficScan, discrNoCutEffic, effVersusDiscr_b, effVersusDiscr_c, effVersusDiscr_d, effVersusDiscr_dus, effVersusDiscr_dusg, effVersusDiscr_g, effVersusDiscr_ni, effVersusDiscr_pu, effVersusDiscr_s, effVersusDiscr_u, FlavourHistograms< T >::getHistoVector(), histoExtension, label, FlavourHistograms< T >::lowerBound(), mcPlots_, seedmultiplicitymonitor_newtracking_cfi::nBins, FlavourHistograms< T >::nBins(), mathSSE::sqrt(), and FlavourHistograms< T >::upperBound().

    : fromDiscriminatorDistr(true),
      mcPlots_(mc),
      doCTagPlots_(false),
      label_(label),
      nBinOutput(nBin),
      startOutput(startO),
      endOutput(endO) {
  histoExtension = "_" + dDiscriminatorFC.baseNameTitle();

  discrNoCutEffic.reset(new FlavourHistograms<double>("totalEntries" + histoExtension,
                                                      "Total Entries: " + dDiscriminatorFC.baseNameDescription(),
                                                      dDiscriminatorFC.nBins(),
                                                      dDiscriminatorFC.lowerBound(),
                                                      dDiscriminatorFC.upperBound(),
                                                      false,
                                                      true,
                                                      false,
                                                      "b",
                                                      label,
                                                      mcPlots_,
                                                      ibook));

  // conditional discriminator cut for efficiency histos

  discrCutEfficScan.reset(new FlavourHistograms<double>("effVsDiscrCut" + histoExtension,
                                                        "Eff. vs Disc. Cut: " + dDiscriminatorFC.baseNameDescription(),
                                                        dDiscriminatorFC.nBins(),
                                                        dDiscriminatorFC.lowerBound(),
                                                        dDiscriminatorFC.upperBound(),
                                                        false,
                                                        true,
                                                        false,
                                                        "b",
                                                        label,
                                                        mcPlots_,
                                                        ibook));
  discrCutEfficScan->SetMinimum(1E-4);
  if (mcPlots_) {
    if (mcPlots_ > 2) {
      effVersusDiscr_d = discrCutEfficScan->histo_d();
      effVersusDiscr_u = discrCutEfficScan->histo_u();
      effVersusDiscr_s = discrCutEfficScan->histo_s();
      effVersusDiscr_g = discrCutEfficScan->histo_g();
      effVersusDiscr_dus = discrCutEfficScan->histo_dus();
    } else {
      effVersusDiscr_d = nullptr;
      effVersusDiscr_u = nullptr;
      effVersusDiscr_s = nullptr;
      effVersusDiscr_g = nullptr;
      effVersusDiscr_dus = nullptr;
    }
    effVersusDiscr_c = discrCutEfficScan->histo_c();
    effVersusDiscr_b = discrCutEfficScan->histo_b();
    effVersusDiscr_ni = discrCutEfficScan->histo_ni();
    effVersusDiscr_dusg = discrCutEfficScan->histo_dusg();
    effVersusDiscr_pu = discrCutEfficScan->histo_pu();

    if (mcPlots_ > 2) {
      effVersusDiscr_d->SetXTitle("Discriminant");
      effVersusDiscr_d->GetXaxis()->SetTitleOffset(0.75);
      effVersusDiscr_u->SetXTitle("Discriminant");
      effVersusDiscr_u->GetXaxis()->SetTitleOffset(0.75);
      effVersusDiscr_s->SetXTitle("Discriminant");
      effVersusDiscr_s->GetXaxis()->SetTitleOffset(0.75);
      effVersusDiscr_g->SetXTitle("Discriminant");
      effVersusDiscr_g->GetXaxis()->SetTitleOffset(0.75);
      effVersusDiscr_dus->SetXTitle("Discriminant");
      effVersusDiscr_dus->GetXaxis()->SetTitleOffset(0.75);
    }
    effVersusDiscr_c->SetXTitle("Discriminant");
    effVersusDiscr_c->GetXaxis()->SetTitleOffset(0.75);
    effVersusDiscr_b->SetXTitle("Discriminant");
    effVersusDiscr_b->GetXaxis()->SetTitleOffset(0.75);
    effVersusDiscr_ni->SetXTitle("Discriminant");
    effVersusDiscr_ni->GetXaxis()->SetTitleOffset(0.75);
    effVersusDiscr_dusg->SetXTitle("Discriminant");
    effVersusDiscr_dusg->GetXaxis()->SetTitleOffset(0.75);
    effVersusDiscr_pu->SetXTitle("Discriminant");
    effVersusDiscr_pu->GetXaxis()->SetTitleOffset(0.75);
  } else {
    effVersusDiscr_d = nullptr;
    effVersusDiscr_u = nullptr;
    effVersusDiscr_s = nullptr;
    effVersusDiscr_c = nullptr;
    effVersusDiscr_b = nullptr;
    effVersusDiscr_g = nullptr;
    effVersusDiscr_ni = nullptr;
    effVersusDiscr_dus = nullptr;
    effVersusDiscr_dusg = nullptr;
    effVersusDiscr_pu = nullptr;
  }

  // discr. for computation
  vector<TH1F*> discrCfHistos = dDiscriminatorFC.getHistoVector();

  // discr no cut
  vector<TH1F*> discrNoCutHistos = discrNoCutEffic->getHistoVector();

  // discr no cut
  vector<TH1F*> discrCutHistos = discrCutEfficScan->getHistoVector();

  const int& dimHistos = discrCfHistos.size();  // they all have the same size

  // DISCR-CUT LOOP:
  // fill the histos for eff-pur computations by scanning the discriminatorFC histogram

  // better to loop over bins -> discrCut no longer needed
  const int& nBins = dDiscriminatorFC.nBins();

  // loop over flavours
  for (int iFlav = 0; iFlav < dimHistos; iFlav++) {
    if (discrCfHistos[iFlav] == nullptr)
      continue;
    discrNoCutHistos[iFlav]->SetXTitle("Discriminant");
    discrNoCutHistos[iFlav]->GetXaxis()->SetTitleOffset(0.75);

    // In Root histos, bin counting starts at 1 to nBins.
    // bin 0 is the underflow, and nBins+1 is the overflow.
    const double& nJetsFlav = discrCfHistos[iFlav]->GetEntries();
    double sum = discrCfHistos[iFlav]->GetBinContent(nBins + 1);  //+1 to get the overflow.

    for (int iDiscr = nBins; iDiscr > 0; --iDiscr) {
      // fill all jets into NoCut histo
      discrNoCutHistos[iFlav]->SetBinContent(iDiscr, nJetsFlav);
      discrNoCutHistos[iFlav]->SetBinError(iDiscr, sqrt(nJetsFlav));
      sum += discrCfHistos[iFlav]->GetBinContent(iDiscr);
      discrCutHistos[iFlav]->SetBinContent(iDiscr, sum);
      discrCutHistos[iFlav]->SetBinError(iDiscr, sqrt(sum));
    }
  }

  // divide to get efficiency vs. discriminator cut from absolute numbers
  discrCutEfficScan->divide(*discrNoCutEffic);  // does: histos including discriminator cut / flat histo
  discrCutEfficScan->setEfficiencyFlag();
}

EffPurFromHistos::~EffPurFromHistos

(

)

Definition at line 197 of file EffPurFromHistos.cc.

{}

Member Function Documentation

void EffPurFromHistos::check ( )

private

Definition at line 361 of file EffPurFromHistos.cc.

References effVersusDiscr_b, effVersusDiscr_c, effVersusDiscr_d, effVersusDiscr_dus, effVersusDiscr_dusg, effVersusDiscr_g, effVersusDiscr_ni, effVersusDiscr_pu, effVersusDiscr_s, effVersusDiscr_u, Exception, and mcPlots_.

Referenced by doCTagPlots(), and EffPurFromHistos().

                             {
  // number of bins

  int nBins_d = 0;
  int nBins_u = 0;
  int nBins_s = 0;
  int nBins_g = 0;
  int nBins_dus = 0;
  if (mcPlots_ > 2) {
    nBins_d = effVersusDiscr_d->GetNbinsX();
    nBins_u = effVersusDiscr_u->GetNbinsX();
    nBins_s = effVersusDiscr_s->GetNbinsX();
    nBins_g = effVersusDiscr_g->GetNbinsX();
    nBins_dus = effVersusDiscr_dus->GetNbinsX();
  }
  const int& nBins_c = effVersusDiscr_c->GetNbinsX();
  const int& nBins_b = effVersusDiscr_b->GetNbinsX();
  const int& nBins_ni = effVersusDiscr_ni->GetNbinsX();
  const int& nBins_dusg = effVersusDiscr_dusg->GetNbinsX();
  const int& nBins_pu = effVersusDiscr_pu->GetNbinsX();

  const bool& lNBins = ((nBins_d == nBins_u && nBins_d == nBins_s && nBins_d == nBins_c && nBins_d == nBins_b &&
                         nBins_d == nBins_g && nBins_d == nBins_ni && nBins_d == nBins_dus && nBins_d == nBins_dusg) ||
                        (nBins_c == nBins_b && nBins_c == nBins_dusg && nBins_c == nBins_ni && nBins_c == nBins_pu));

  if (!lNBins) {
    throw cms::Exception("Configuration") << "Input histograms do not all have the same number of bins!\n";
  }

  // start
  float sBin_d = 0;
  float sBin_u = 0;
  float sBin_s = 0;
  float sBin_g = 0;
  float sBin_dus = 0;
  if (mcPlots_ > 2) {
    sBin_d = effVersusDiscr_d->GetBinCenter(1);
    sBin_u = effVersusDiscr_u->GetBinCenter(1);
    sBin_s = effVersusDiscr_s->GetBinCenter(1);
    sBin_g = effVersusDiscr_g->GetBinCenter(1);
    sBin_dus = effVersusDiscr_dus->GetBinCenter(1);
  }
  const float& sBin_c = effVersusDiscr_c->GetBinCenter(1);
  const float& sBin_b = effVersusDiscr_b->GetBinCenter(1);
  const float& sBin_ni = effVersusDiscr_ni->GetBinCenter(1);
  const float& sBin_dusg = effVersusDiscr_dusg->GetBinCenter(1);
  const float& sBin_pu = effVersusDiscr_pu->GetBinCenter(1);

  const bool& lSBin = ((sBin_d == sBin_u && sBin_d == sBin_s && sBin_d == sBin_c && sBin_d == sBin_b &&
                        sBin_d == sBin_g && sBin_d == sBin_ni && sBin_d == sBin_dus && sBin_d == sBin_dusg) ||
                       (sBin_c == sBin_b && sBin_c == sBin_dusg && sBin_c == sBin_ni && sBin_c == sBin_pu));

  if (!lSBin) {
    throw cms::Exception("Configuration")
        << "EffPurFromHistos::check() : Input histograms do not all have the same start bin!\n";
  }

  // end
  float eBin_d = 0;
  float eBin_u = 0;
  float eBin_s = 0;
  float eBin_g = 0;
  float eBin_dus = 0;
  if (mcPlots_ > 2) {
    eBin_d = effVersusDiscr_d->GetBinCenter(nBins_d - 1);
    eBin_u = effVersusDiscr_u->GetBinCenter(nBins_d - 1);
    eBin_s = effVersusDiscr_s->GetBinCenter(nBins_d - 1);
    eBin_g = effVersusDiscr_g->GetBinCenter(nBins_d - 1);
    eBin_dus = effVersusDiscr_dus->GetBinCenter(nBins_d - 1);
  }
  const float& eBin_c = effVersusDiscr_c->GetBinCenter(nBins_d - 1);
  const float& eBin_b = effVersusDiscr_b->GetBinCenter(nBins_d - 1);
  const float& eBin_ni = effVersusDiscr_ni->GetBinCenter(nBins_d - 1);
  const float& eBin_dusg = effVersusDiscr_dusg->GetBinCenter(nBins_d - 1);
  const float& eBin_pu = effVersusDiscr_pu->GetBinCenter(nBins_d - 1);

  const bool& lEBin = ((eBin_d == eBin_u && eBin_d == eBin_s && eBin_d == eBin_c && eBin_d == eBin_b &&
                        eBin_d == eBin_g && eBin_d == eBin_ni && eBin_d == eBin_dus && eBin_d == eBin_dusg) ||
                       (eBin_c == eBin_b && eBin_c == eBin_dusg && eBin_c == eBin_ni && eBin_c == eBin_pu));

  if (!lEBin) {
    throw cms::Exception("Configuration")
        << "EffPurFromHistos::check() : Input histograms do not all have the same end bin!\n";
  }
}

void EffPurFromHistos::compute

(

DQMStore::IBooker &

ibook

)

Definition at line 447 of file EffPurFromHistos.cc.

References HistoProviderDQM::book1D(), doCTagPlots_, EffFlavVsXEff_b, EffFlavVsXEff_c, EffFlavVsXEff_d, EffFlavVsXEff_dus, EffFlavVsXEff_dusg, EffFlavVsXEff_g, EffFlavVsXEff_ni, EffFlavVsXEff_pu, EffFlavVsXEff_s, EffFlavVsXEff_u, effVersusDiscr_b, effVersusDiscr_c, effVersusDiscr_d, effVersusDiscr_dus, effVersusDiscr_dusg, effVersusDiscr_g, effVersusDiscr_ni, effVersusDiscr_pu, effVersusDiscr_s, effVersusDiscr_u, endOutput, dqm::impl::MonitorElement::Fill(), RecoBTag::findBinClosestYValue(), dqm::legacy::MonitorElement::getTH1F(), histoExtension, label_, mcPlots_, nBinOutput, anotherprimaryvertexanalyzer_cfi::nBinX, dqm::impl::MonitorElement::setEfficiencyFlag(), startOutput, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                     {
  if (!mcPlots_) {
    EffFlavVsXEff_d = nullptr;
    EffFlavVsXEff_u = nullptr;
    EffFlavVsXEff_s = nullptr;
    EffFlavVsXEff_c = nullptr;
    EffFlavVsXEff_b = nullptr;
    EffFlavVsXEff_g = nullptr;
    EffFlavVsXEff_ni = nullptr;
    EffFlavVsXEff_dus = nullptr;
    EffFlavVsXEff_dusg = nullptr;
    EffFlavVsXEff_pu = nullptr;
    return;
  }

  // to have shorter names ......
  const std::string& hE = histoExtension;
  std::string hX = "";
  TString Title = "";
  if (!doCTagPlots_) {
    hX = "FlavEffVsBEff_";
    Title = "b";
  } else {
    hX = "FlavEffVsCEff_";
    Title = "c";
  }

  // create histograms from base name and extension as given from user
  // BINNING MUST BE IDENTICAL FOR ALL OF THEM!!
  HistoProviderDQM prov("Btag", label_, ibook);
  if (mcPlots_ > 2) {
    EffFlavVsXEff_d = (prov.book1D(hX + "D" + hE, hX + "D" + hE, nBinOutput, startOutput, endOutput));
    EffFlavVsXEff_d->setEfficiencyFlag();
    EffFlavVsXEff_u = (prov.book1D(hX + "U" + hE, hX + "U" + hE, nBinOutput, startOutput, endOutput));
    EffFlavVsXEff_u->setEfficiencyFlag();
    EffFlavVsXEff_s = (prov.book1D(hX + "S" + hE, hX + "S" + hE, nBinOutput, startOutput, endOutput));
    EffFlavVsXEff_s->setEfficiencyFlag();
    EffFlavVsXEff_g = (prov.book1D(hX + "G" + hE, hX + "G" + hE, nBinOutput, startOutput, endOutput));
    EffFlavVsXEff_g->setEfficiencyFlag();
    EffFlavVsXEff_dus = (prov.book1D(hX + "DUS" + hE, hX + "DUS" + hE, nBinOutput, startOutput, endOutput));
    EffFlavVsXEff_dus->setEfficiencyFlag();
  } else {
    EffFlavVsXEff_d = nullptr;
    EffFlavVsXEff_u = nullptr;
    EffFlavVsXEff_s = nullptr;
    EffFlavVsXEff_g = nullptr;
    EffFlavVsXEff_dus = nullptr;
  }
  EffFlavVsXEff_c = (prov.book1D(hX + "C" + hE, hX + "C" + hE, nBinOutput, startOutput, endOutput));
  EffFlavVsXEff_c->setEfficiencyFlag();
  EffFlavVsXEff_b = (prov.book1D(hX + "B" + hE, hX + "B" + hE, nBinOutput, startOutput, endOutput));
  EffFlavVsXEff_b->setEfficiencyFlag();
  EffFlavVsXEff_ni = (prov.book1D(hX + "NI" + hE, hX + "NI" + hE, nBinOutput, startOutput, endOutput));
  EffFlavVsXEff_ni->setEfficiencyFlag();
  EffFlavVsXEff_dusg = (prov.book1D(hX + "DUSG" + hE, hX + "DUSG" + hE, nBinOutput, startOutput, endOutput));
  EffFlavVsXEff_dusg->setEfficiencyFlag();
  EffFlavVsXEff_pu = (prov.book1D(hX + "PU" + hE, hX + "PU" + hE, nBinOutput, startOutput, endOutput));
  EffFlavVsXEff_pu->setEfficiencyFlag();

  if (mcPlots_ > 2) {
    EffFlavVsXEff_d->getTH1F()->SetXTitle(Title + "-jet efficiency");
    EffFlavVsXEff_d->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
    EffFlavVsXEff_d->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_d->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_u->getTH1F()->SetXTitle(Title + "-jet efficiency");
    EffFlavVsXEff_u->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
    EffFlavVsXEff_u->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_u->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_s->getTH1F()->SetXTitle(Title + "-jet efficiency");
    EffFlavVsXEff_s->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
    EffFlavVsXEff_s->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_s->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_g->getTH1F()->SetXTitle(Title + "-jet efficiency");
    EffFlavVsXEff_g->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
    EffFlavVsXEff_g->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_g->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_dus->getTH1F()->SetXTitle(Title + "-jet efficiency");
    EffFlavVsXEff_dus->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
    EffFlavVsXEff_dus->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
    EffFlavVsXEff_dus->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
  }
  EffFlavVsXEff_c->getTH1F()->SetXTitle(Title + "-jet efficiency");
  EffFlavVsXEff_c->getTH1F()->SetYTitle("c-jet efficiency");
  EffFlavVsXEff_c->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_c->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_b->getTH1F()->SetXTitle(Title + "-jet efficiency");
  EffFlavVsXEff_b->getTH1F()->SetYTitle("b-jet efficiency");
  EffFlavVsXEff_b->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_b->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_ni->getTH1F()->SetXTitle(Title + "-jet efficiency");
  EffFlavVsXEff_ni->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
  EffFlavVsXEff_ni->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_ni->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_dusg->getTH1F()->SetXTitle(Title + "-jet efficiency");
  EffFlavVsXEff_dusg->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
  EffFlavVsXEff_dusg->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_dusg->getTH1F()->GetYaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_pu->getTH1F()->SetXTitle(Title + "-jet efficiency");
  EffFlavVsXEff_pu->getTH1F()->SetYTitle("non " + Title + "-jet efficiency");
  EffFlavVsXEff_pu->getTH1F()->GetXaxis()->SetTitleOffset(0.75);
  EffFlavVsXEff_pu->getTH1F()->GetYaxis()->SetTitleOffset(0.75);

  // loop over eff. vs. discriminator cut b-histo and look in which bin the closest entry is;
  // use fact that eff decreases monotonously

  // any of the histos to be created can be taken here:
  MonitorElement* EffFlavVsXEff = EffFlavVsXEff_b;

  const int& nBinX = EffFlavVsXEff->getTH1F()->GetNbinsX();

  for (int iBinX = 1; iBinX <= nBinX; iBinX++) {  // loop over the bins on the x-axis of the histograms to be filled

    const float& effXBinWidth = EffFlavVsXEff->getTH1F()->GetBinWidth(iBinX);
    const float& effXMid = EffFlavVsXEff->getTH1F()->GetBinCenter(iBinX);  // middle of b-efficiency bin
    const float& effXLeft = effXMid - 0.5 * effXBinWidth;                  // left edge of bin
    const float& effXRight = effXMid + 0.5 * effXBinWidth;                 // right edge of bin
    // find the corresponding bin in the efficiency versus discriminator cut histo: closest one in efficiency

    int binClosest = -1;
    if (!doCTagPlots_) {
      binClosest = findBinClosestYValue(effVersusDiscr_b, effXMid, effXLeft, effXRight);
    } else {
      binClosest = findBinClosestYValue(effVersusDiscr_c, effXMid, effXLeft, effXRight);
    }

    const bool& binFound = (binClosest > 0);
    //
    if (binFound) {
      // fill the histos
      if (mcPlots_ > 2) {
        EffFlavVsXEff_d->Fill(effXMid, effVersusDiscr_d->GetBinContent(binClosest));
        EffFlavVsXEff_u->Fill(effXMid, effVersusDiscr_u->GetBinContent(binClosest));
        EffFlavVsXEff_s->Fill(effXMid, effVersusDiscr_s->GetBinContent(binClosest));
        EffFlavVsXEff_g->Fill(effXMid, effVersusDiscr_g->GetBinContent(binClosest));
        EffFlavVsXEff_dus->Fill(effXMid, effVersusDiscr_dus->GetBinContent(binClosest));
      }
      EffFlavVsXEff_c->Fill(effXMid, effVersusDiscr_c->GetBinContent(binClosest));
      EffFlavVsXEff_b->Fill(effXMid, effVersusDiscr_b->GetBinContent(binClosest));
      EffFlavVsXEff_ni->Fill(effXMid, effVersusDiscr_ni->GetBinContent(binClosest));
      EffFlavVsXEff_dusg->Fill(effXMid, effVersusDiscr_dusg->GetBinContent(binClosest));
      EffFlavVsXEff_pu->Fill(effXMid, effVersusDiscr_pu->GetBinContent(binClosest));

      if (mcPlots_ > 2) {
        EffFlavVsXEff_d->getTH1F()->SetBinError(iBinX, effVersusDiscr_d->GetBinError(binClosest));
        EffFlavVsXEff_u->getTH1F()->SetBinError(iBinX, effVersusDiscr_u->GetBinError(binClosest));
        EffFlavVsXEff_s->getTH1F()->SetBinError(iBinX, effVersusDiscr_s->GetBinError(binClosest));
        EffFlavVsXEff_g->getTH1F()->SetBinError(iBinX, effVersusDiscr_g->GetBinError(binClosest));
        EffFlavVsXEff_dus->getTH1F()->SetBinError(iBinX, effVersusDiscr_dus->GetBinError(binClosest));
      }
      EffFlavVsXEff_c->getTH1F()->SetBinError(iBinX, effVersusDiscr_c->GetBinError(binClosest));
      EffFlavVsXEff_b->getTH1F()->SetBinError(iBinX, effVersusDiscr_b->GetBinError(binClosest));
      EffFlavVsXEff_ni->getTH1F()->SetBinError(iBinX, effVersusDiscr_ni->GetBinError(binClosest));
      EffFlavVsXEff_dusg->getTH1F()->SetBinError(iBinX, effVersusDiscr_dusg->GetBinError(binClosest));
      EffFlavVsXEff_pu->getTH1F()->SetBinError(iBinX, effVersusDiscr_pu->GetBinError(binClosest));
    }
  }
}

FlavourHistograms<double>& EffPurFromHistos::discriminatorCutEfficScan ( ) const

inline

Definition at line 67 of file EffPurFromHistos.h.

References discrCutEfficScan.

{ return *discrCutEfficScan; }

FlavourHistograms<double>& EffPurFromHistos::discriminatorNoCutEffic ( ) const

inline

Definition at line 66 of file EffPurFromHistos.h.

References discrNoCutEffic.

{ return *discrNoCutEffic; }

bool EffPurFromHistos::doCTagPlots ( bool  Ctag )

inline

Definition at line 69 of file EffPurFromHistos.h.

References check(), and doCTagPlots_.

                              {
    doCTagPlots_ = Ctag;
    return doCTagPlots_;
  }

void EffPurFromHistos::epsPlot

(

const std::string &

name

)

Definition at line 199 of file EffPurFromHistos.cc.

References discrCutEfficScan, discrNoCutEffic, fromDiscriminatorDistr, and plot().

Referenced by getEffFlavVsBEff_pu().

                                                    {
  if (fromDiscriminatorDistr) {
    discrNoCutEffic->epsPlot(name);
    discrCutEfficScan->epsPlot(name);
  }
  plot(name, ".eps");
}

TH1F* EffPurFromHistos::getEffFlavVsBEff_b ( )

inline

Definition at line 53 of file EffPurFromHistos.h.

References EffFlavVsXEff_b, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_b->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_c ( )

inline

Definition at line 52 of file EffPurFromHistos.h.

References EffFlavVsXEff_c, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_c->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_d ( )

inline

Definition at line 49 of file EffPurFromHistos.h.

References EffFlavVsXEff_d, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_d->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_dus ( )

inline

Definition at line 56 of file EffPurFromHistos.h.

References EffFlavVsXEff_dus, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_dus->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_dusg ( )

inline

Definition at line 57 of file EffPurFromHistos.h.

References EffFlavVsXEff_dusg, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_dusg->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_g ( )

inline

Definition at line 54 of file EffPurFromHistos.h.

References EffFlavVsXEff_g, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_g->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_ni ( )

inline

Definition at line 55 of file EffPurFromHistos.h.

References EffFlavVsXEff_ni, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_ni->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_pu ( )

inline

Definition at line 58 of file EffPurFromHistos.h.

References EffFlavVsXEff_pu, epsPlot(), dqm::legacy::MonitorElement::getTH1F(), Skims_PA_cff::name, plot(), psPlot(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_pu->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_s ( )

inline

Definition at line 51 of file EffPurFromHistos.h.

References EffFlavVsXEff_s, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_s->getTH1F(); };

TH1F* EffPurFromHistos::getEffFlavVsBEff_u ( )

inline

Definition at line 50 of file EffPurFromHistos.h.

References EffFlavVsXEff_u, and dqm::legacy::MonitorElement::getTH1F().

Referenced by BTagDifferentialPlot::fillHisto().

{ return EffFlavVsXEff_u->getTH1F(); };

void EffPurFromHistos::plot

(

TPad *

theCanvas = 0

)

Definition at line 225 of file EffPurFromHistos.cc.

References doCTagPlots_, MillePedeFileConverter_cfg::e, EffFlavVsXEff_b, EffFlavVsXEff_c, EffFlavVsXEff_d, EffFlavVsXEff_dus, EffFlavVsXEff_dusg, EffFlavVsXEff_g, EffFlavVsXEff_ni, EffFlavVsXEff_pu, EffFlavVsXEff_s, EffFlavVsXEff_u, dqm::legacy::MonitorElement::getTH1F(), mcPlots_, and plotscripts::setTDRStyle().

Referenced by epsPlot(), getEffFlavVsBEff_pu(), plot(), and psPlot().

                                              {
  //fixme:
  bool btppNI = false;
  bool btppColour = true;

  //   if ( !btppTitle ) gStyle->SetOptTitle( 0 );
  setTDRStyle()->cd();

  if (plotCanvas)
    plotCanvas->cd();

  gPad->UseCurrentStyle();
  gPad->SetFillColor(0);
  gPad->SetLogy(1);
  gPad->SetGridx(1);
  gPad->SetGridy(1);

  int col_c;
  int col_g;
  int col_dus;
  int col_ni;

  int mStyle_c;
  int mStyle_g;
  int mStyle_dus;
  int mStyle_ni;

  // marker size(same for all)
  float mSize = gPad->GetWh() * gPad->GetHNDC() / 500.;  //1.2;

  if (btppColour) {
    col_c = 6;
    col_g = 3;    // g in green
    col_dus = 4;  // uds in blue
    col_ni = 5;   // ni in ??
    mStyle_c = 20;
    mStyle_g = 20;
    mStyle_dus = 20;
    mStyle_ni = 20;
  } else {
    col_c = 1;
    col_g = 1;
    col_dus = 1;
    col_ni = 1;
    mStyle_c = 22;
    mStyle_g = 29;
    mStyle_dus = 20;
    mStyle_ni = 27;
  }

  TString Title = "";
  if (!doCTagPlots_) {
    Title = "b";
  } else {
    Title = "c";
  }

  // for the moment: plot c,dus,g
  if (mcPlots_ > 2) {
    EffFlavVsXEff_dus->getTH1F()->GetXaxis()->SetTitle(Title + "-jet efficiency");
    EffFlavVsXEff_dus->getTH1F()->GetYaxis()->SetTitle("non " + Title + "-jet efficiency");
    EffFlavVsXEff_dus->getTH1F()->GetYaxis()->SetTitleOffset(0.25);
    EffFlavVsXEff_dus->getTH1F()->SetMaximum(1.1);
    EffFlavVsXEff_dus->getTH1F()->SetMinimum(1.e-5);
    EffFlavVsXEff_dus->getTH1F()->SetMarkerColor(col_dus);
    EffFlavVsXEff_dus->getTH1F()->SetLineColor(col_dus);
    EffFlavVsXEff_dus->getTH1F()->SetMarkerSize(mSize);
    EffFlavVsXEff_dus->getTH1F()->SetMarkerStyle(mStyle_dus);
    EffFlavVsXEff_dus->getTH1F()->SetStats(false);
    EffFlavVsXEff_dus->getTH1F()->Draw("pe");

    EffFlavVsXEff_g->getTH1F()->SetMarkerColor(col_g);
    EffFlavVsXEff_g->getTH1F()->SetLineColor(col_g);
    EffFlavVsXEff_g->getTH1F()->SetMarkerSize(mSize);
    EffFlavVsXEff_g->getTH1F()->SetMarkerStyle(mStyle_g);
    EffFlavVsXEff_g->getTH1F()->SetStats(false);
    EffFlavVsXEff_g->getTH1F()->Draw("peSame");
  }
  EffFlavVsXEff_c->getTH1F()->SetMarkerColor(col_c);
  EffFlavVsXEff_c->getTH1F()->SetLineColor(col_c);
  EffFlavVsXEff_c->getTH1F()->SetMarkerSize(mSize);
  EffFlavVsXEff_c->getTH1F()->SetMarkerStyle(mStyle_c);
  EffFlavVsXEff_c->getTH1F()->SetStats(false);
  EffFlavVsXEff_c->getTH1F()->Draw("peSame");

  if (mcPlots_ > 2) {
    EffFlavVsXEff_d->getTH1F()->SetMinimum(0.01);
    EffFlavVsXEff_u->getTH1F()->SetMinimum(0.01);
    EffFlavVsXEff_s->getTH1F()->SetMinimum(0.01);
    EffFlavVsXEff_g->getTH1F()->SetMinimum(0.01);
    EffFlavVsXEff_dus->getTH1F()->SetMinimum(0.01);
  }
  EffFlavVsXEff_c->getTH1F()->SetMinimum(0.01);
  EffFlavVsXEff_b->getTH1F()->SetMinimum(0.01);
  EffFlavVsXEff_ni->getTH1F()->SetMinimum(0.01);
  EffFlavVsXEff_dusg->getTH1F()->SetMinimum(0.01);
  EffFlavVsXEff_pu->getTH1F()->SetMinimum(0.01);

  // plot separately u,d and s
  //  EffFlavVsXEff_d ->GetXaxis()->SetTitle( Title + "-jet efficiency" );
  //  EffFlavVsXEff_d ->GetYaxis()->SetTitle( "non " + Title + "-jet efficiency" );
  //  EffFlavVsXEff_d ->GetYaxis()->SetTitleOffset( 1.25 );
  //  EffFlavVsXEff_d ->SetMaximum    ( 1.1 );
  //  EffFlavVsXEff_d ->SetMinimum    ( 1.e-5 );
  //  EffFlavVsXEff_d ->SetMarkerColor( col_dus );
  //  EffFlavVsXEff_d ->SetLineColor  ( col_dus );
  //  EffFlavVsXEff_d ->SetMarkerSize ( mSize );
  //  EffFlavVsXEff_d ->SetMarkerStyle( mStyle_dus );
  //  EffFlavVsXEff_d ->SetStats    ( false );
  //  EffFlavVsXEff_d ->Draw("pe");
  //
  //  EffFlavVsXEff_u   ->SetMarkerColor( col_g );
  //  EffFlavVsXEff_u   ->SetLineColor  ( col_g );
  //  EffFlavVsXEff_u   ->SetMarkerSize ( mSize );
  //  EffFlavVsXEff_u   ->SetMarkerStyle( mStyle_g );
  //  EffFlavVsXEff_u   ->SetStats    ( false );
  //  EffFlavVsXEff_u   ->Draw("peSame");
  //
  //  EffFlavVsXEff_s   ->SetMarkerColor( col_c );
  //  EffFlavVsXEff_s   ->SetLineColor  ( col_c );
  //  EffFlavVsXEff_s   ->SetMarkerSize ( mSize );
  //  EffFlavVsXEff_s   ->SetMarkerStyle( mStyle_c );
  //  EffFlavVsXEff_s   ->SetStats    ( false );
  //  EffFlavVsXEff_s   ->Draw("peSame");

  // only if asked: NI
  if (btppNI) {
    EffFlavVsXEff_ni->getTH1F()->SetMarkerColor(col_ni);
    EffFlavVsXEff_ni->getTH1F()->SetLineColor(col_ni);
    EffFlavVsXEff_ni->getTH1F()->SetMarkerSize(mSize);
    EffFlavVsXEff_ni->getTH1F()->SetMarkerStyle(mStyle_ni);
    EffFlavVsXEff_ni->getTH1F()->SetStats(false);
    EffFlavVsXEff_ni->getTH1F()->Draw("peSame");
  }
}

void EffPurFromHistos::plot	(	const std::string &	name,
		const std::string &	ext
	)

Definition at line 209 of file EffPurFromHistos.cc.

References doCTagPlots_, histoExtension, plot(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                       {
  std::string hX = "";
  std::string Title = "";
  if (!doCTagPlots_) {
    hX = "FlavEffVsBEff";
    Title = "b";
  } else {
    hX = "FlavEffVsCEff";
    Title = "c";
  }
  TCanvas tc((hX + histoExtension).c_str(),
             ("Flavour misidentification vs. " + Title + "-tagging efficiency " + histoExtension).c_str());
  plot(&tc);
  tc.Print((name + hX + histoExtension + ext).c_str());
}

void EffPurFromHistos::psPlot

(

const std::string &

name

)

Definition at line 207 of file EffPurFromHistos.cc.

References plot().

Referenced by getEffFlavVsBEff_pu().

{ plot(name, ".ps"); }

Member Data Documentation

std::unique_ptr<FlavourHistograms<double> > EffPurFromHistos::discrCutEfficScan

private

Definition at line 85 of file EffPurFromHistos.h.

Referenced by discriminatorCutEfficScan(), EffPurFromHistos(), and epsPlot().

std::unique_ptr<FlavourHistograms<double> > EffPurFromHistos::discrNoCutEffic

private

Definition at line 85 of file EffPurFromHistos.h.

Referenced by discriminatorNoCutEffic(), EffPurFromHistos(), and epsPlot().

bool EffPurFromHistos::doCTagPlots_

private

Definition at line 80 of file EffPurFromHistos.h.

Referenced by compute(), doCTagPlots(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_b

private

Definition at line 113 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_b(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_c

private

Definition at line 112 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_c(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_d

private

Definition at line 109 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_d(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_dus

private

Definition at line 116 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_dus(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_dusg

private

Definition at line 117 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_dusg(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_g

private

Definition at line 114 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_g(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_ni

private

Definition at line 115 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_ni(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_pu

private

Definition at line 118 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_pu(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_s

private

Definition at line 111 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_s(), and plot().

MonitorElement* EffPurFromHistos::EffFlavVsXEff_u

private

Definition at line 110 of file EffPurFromHistos.h.

Referenced by compute(), getEffFlavVsBEff_u(), and plot().

TH1F* EffPurFromHistos::effVersusDiscr_b

private

Definition at line 95 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_c

private

Definition at line 94 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_d

private

Definition at line 91 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_dus

private

Definition at line 98 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_dusg

private

Definition at line 99 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_g

private

Definition at line 96 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_ni

private

Definition at line 97 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_pu

private

Definition at line 100 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_s

private

Definition at line 93 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

TH1F* EffPurFromHistos::effVersusDiscr_u

private

Definition at line 92 of file EffPurFromHistos.h.

Referenced by check(), compute(), and EffPurFromHistos().

double EffPurFromHistos::endOutput

private

Definition at line 107 of file EffPurFromHistos.h.

Referenced by compute().

bool EffPurFromHistos::fromDiscriminatorDistr

private

Definition at line 77 of file EffPurFromHistos.h.

Referenced by epsPlot().

std::string EffPurFromHistos::histoExtension

private

Definition at line 83 of file EffPurFromHistos.h.

Referenced by compute(), EffPurFromHistos(), and plot().

std::string EffPurFromHistos::label_

private

Definition at line 81 of file EffPurFromHistos.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), compute(), and Mixins._Labelable::setLabel().

unsigned int EffPurFromHistos::mcPlots_

private

Definition at line 79 of file EffPurFromHistos.h.

Referenced by check(), compute(), EffPurFromHistos(), and plot().

int EffPurFromHistos::nBinOutput

private

Definition at line 105 of file EffPurFromHistos.h.

Referenced by compute().

double EffPurFromHistos::startOutput

private

Definition at line 106 of file EffPurFromHistos.h.

Referenced by compute().