15 const vector<pair<string, Color_t > > int_mode_colors
25 {
"selectedVtxFailure",
kGreen-5}
30 {
"numucc",
"#nu_{#mu} CC"},
31 {
"numubarcc",
"#nu_{#mu} CC"},
35 {
"nueorbar",
"#nu_e + #bar{#nu}_e"},
37 {
"otherbkgd",
"Other Bkgd"},
38 {
"selectedVtxFailure",
"Vtx Fail"}
45 void plot_muonid_dist(
string inputFilename,
string mode_str,
string outputFile =
"muonid_dist.pdf")
48 if (mode_str.tolower() ==
"fhc")
50 else if (mode_str.tolower() ==
"rhc")
53 throw std::runtime_error(
"Invalid beam mode, plesae give 'FHC' or 'RHC'");
55 TFile *
infile = TFile::Open(inputFilename.c_str());
56 TCanvas *
c1 =
new TCanvas(
"c1",
"c1", 1600, 1200);
57 TLegend *
leg =
new TLegend();
58 bool firstplot =
true;
59 for (pair<string, Color_t> int_mode : int_mode_colors)
62 hist->SetName((
"muonid_" + int_mode.first).c_str());
63 hist->SetTitle(
";MuonID;Events / 12.5E20POT");
64 hist->SetLineColor(int_mode.second);
65 leg->AddEntry(hist, int_mode_titles[int_mode.first].c_str());
72 hist->Draw(
"hist same");
75 c1->SaveAs(outputFile.c_str());
std::unique_ptr< T > LoadFrom(TDirectory *dir, const std::string &label)
std::vector< float > Spectrum
const double kAna2020RHCPOT
void plot_muonid_dist(string inputFilename, string mode_str, string outputFile="muonid_dist.pdf")
map< string, string > int_mode_titles_fhc