nus_error_barchart.C

Go to the documentation of this file.

//////////////////////////////////////////////////////////////////////////

#include <string>
#include <map>

#include "TCanvas.h"

#include "CAFAna/Analysis/Plots.h"

const std::string outdir = "/nova/ana/users/gsdavies/nus/Ana01";

void nus_error_barchart()
{

  new TCanvas;
  std::map<std::string, double> sig_error
  {
    {"ND Composition",           7.0 },
    {"Calibration",              5.8 },
    {"Normalization",            4.9 },
    {"ND External Activity",     4.1 },
    {"Beam Flux",                3.4 },
    {"Scintillation Model",      2.4 },
    {"Simulation Statistics",    2.0 },
    {"Neutrino Interaction",     1.6 },
    {"Acceptance",               1.0 },
    {"Osc. Parameters",          0.7 },
  };

  double statErrorSig = 12.8;
  ana::CountingExperimentErrorBarChart(sig_error, statErrorSig, true);
  gPad->Print((outdir + "/nus_ana01_unc_sig.pdf").c_str());
  gPad->Print((outdir + "/nus_ana01_unc_sig.eps").c_str());
  gPad->Print((outdir + "/nus_ana01_unc_sig.png").c_str());

  new TCanvas;
  std::map<std::string, double> bkgd_error
  {
    {"ND Composition",           10.4 },
    {"Calibration",              6.0  },
    {"Normalization",            4.9  },
    {"ND External Activity",     1.7  },
    {"Beam Flux",                3.6  },
    {"Scintillation Model",      1.8  },
    {"Simulation Statistics",    4.8  },
    {"Neutrino Interaction",     4.8  },
    {"Acceptance",               0.6  },
    {"Osc. Parameters",          10.7 },
  };
  double statErrBkgd = 20.9;
  ana::CountingExperimentErrorBarChart(bkgd_error, statErrBkgd, false);
  gPad->Print((outdir + "/nus_ana01_unc_bkgd.pdf").c_str());
  gPad->Print((outdir + "/nus_ana01_unc_bkgd.eps").c_str());
  gPad->Print((outdir + "/nus_ana01_unc_bkgd.png").c_str());

}