draw_decomp_plots.C File Reference

#include "OscLib/OscCalcDumb.h"
#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/NueVars.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "CAFAna/Cuts/NueCutsSecondAna.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/LoadFromFile.h"
#include "CAFAna/Analysis/Style.h"
#include "CAFAna/Analysis/SALoaders.h"
#include "3FlavorAna/Plotting/NuePlotStyle.h"
#include "3FlavorAna/Decomp/MichelDecomp.h"
#include "3FlavorAna/Decomp/BENDecomp.h"
#include "CAFAna/Decomp/ProportionalDecomp.h"
#include "CAFAna/Decomp/CheatDecomp.h"
#include "3FlavorAna/Decomp/NumuDecomp.h"
#include "CAFAna/Extrap/ModularExtrap.h"
#include "CAFAna/Prediction/PredictionExtrap.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"
#include "CAFAna/Systs/Systs.h"
#include "CAFAna/Core/SystShifts.h"
#include "Utilities/func/MathUtil.h"
#include "TROOT.h"
#include "TCanvas.h"
#include "TH2.h"
#include "TFile.h"
#include "TLegend.h"
#include "TProfile.h"
#include "TRandom3.h"
#include "THStack.h"
#include "TGaxis.h"
#include "TLine.h"
#include "TStyle.h"
#include "TLatex.h"
#include "TGraph.h"
#include "TColor.h"
#include <iostream>
#include <cmath>

Go to the source code of this file.

Functions
void	Preliminary ()
void	PIDDivisions ()
void	PIDBinLabels (std::string pid)
void	FillWithDimmerColor (TH1 *h, float dim=0.65, bool usealpha=false)
std::vector< TH1 * >	NDHists (MichelDecomp *mb, std::string s)
void	PlotComp (std::vector< TH1 * > hs, std::string label)
void	PlotRat (std::vector< TH1 * > hs, std::string label)
void	DrawRatioAnaBin (TH1 *cc, TH1 *nc, TH1 *ne, int idx)
void	PlotRatio (std::vector< TH1 * > ccs, std::vector< TH1 * > ncs, std::vector< TH1 * > nes)
void	PlotStack (std::vector< TH1 * > ccs, std::vector< TH1 * > ncs, std::vector< TH1 * > nes, bool PlotMC=true, bool PlotDecomp=true)
void	PlotNME (TFile *in)
void	MakeDCMPPlots (MichelDecomp *mb)
void	draw_decomp_plots ()

Function Documentation

void draw_decomp_plots

(

)

Definition at line 510 of file draw_decomp_plots.C.

References calc, om::cout, allTimeWatchdog::endl, in, ana::MichelDecomp::LoadFrom(), MakeDCMPPlots(), RunSnowGlobes::mb, ana::MichelDecomp::NueComponent(), PlotNME(), and ana::Spectrum::POT().

{

  TFile *in = TFile::Open("mdcmp_hists.root");

  in->cd("decomp");
  MichelDecomp* mb = MichelDecomp::LoadFrom(gDirectory).release();
  double ndpot = mb->NueComponent().POT();
  std::cout << "We see a total of " << ndpot << " POT" << std::endl;
  osc::OscCalcDumb calc;

  MakeDCMPPlots(mb);
  PlotNME(in);
}

void DrawRatioAnaBin	(	TH1 *	cc,
		TH1 *	nc,
		TH1 *	ne,
		int	idx
	)

Definition at line 232 of file draw_decomp_plots.C.

References MECModelEnuComparisons::i, and ana::UniqueName().

Referenced by PlotRatio().

{
  TH1 *tcc = (TH1F*)cc->Clone(UniqueName().c_str());
  TH1 *tnc = (TH1F*)nc->Clone(UniqueName().c_str());
  TH1 *tne = (TH1F*)ne->Clone(UniqueName().c_str());

  for (int i = 1; i <= 27; i++){
    if ((i-1)/9 != idx || (i-1)%9 == 0){
      tcc->SetBinContent(i,-1);
      tnc->SetBinContent(i,-1);
      tne->SetBinContent(i,-1);
    }
  }
  tcc->Draw("][,hist,same");
  tnc->Draw("][,hist,same");
  tne->Draw("][,hist,same");
}

void FillWithDimmerColor	(	TH1 *	h,
		float	dim = 0.65,
		bool	usealpha = false
	)

Definition at line 110 of file draw_decomp_plots.C.

References plot_validation_datamc_2018::color, and R.

Referenced by PlotNME(), and PlotStack().

{
  if ( usealpha ){
    h->SetFillColorAlpha(h->GetLineColor(),dim);
    return;
  }

  TColor *color = gROOT->GetColor(h->GetLineColor());
  float R,G,B,hR,hG,hB,hHue,hSat,hVal;

  color->GetRGB(hR,hG,hB);
  color->RGB2HSV(hR,hG,hB,hHue,hSat,hVal);
  color->HSV2RGB(hHue,dim*hSat,hVal,R,G,B);

  h->SetFillColor(color->GetColor(R,G,B));
}

void MakeDCMPPlots

(

MichelDecomp *

mb

)

Definition at line 498 of file draw_decomp_plots.C.

References NDHists(), and PlotRatio().

Referenced by draw_decomp_plots().

{
  std::vector<TH1*> ccs = NDHists(mb, "cc");
  std::vector<TH1*> ncs = NDHists(mb, "nc");
  std::vector<TH1*> nes = NDHists(mb, "ne");

  PlotRatio(ccs,ncs,nes);
  //PlotStack(ccs,ncs,nes,true,true);  // MC + Decomp
  //PlotStack(ccs,ncs,nes,true,false); // MC Only
  //PlotStack(ccs,ncs,nes,false,true); // Decomp Only
}

std::vector<TH1*> NDHists	(	MichelDecomp *	mb,
		std::string	s
	)

Definition at line 129 of file draw_decomp_plots.C.

References ana::MichelDecomp::AntiNueComponent(), ana::MichelDecomp::AntiNumuComponent(), om::cout, allTimeWatchdog::endl, MECModelEnuComparisons::i, kGreen, kRed, cet::sqlite::max(), ana::MichelDecomp::MC_AntiNueComponent(), ana::MichelDecomp::MC_AntiNumuComponent(), ana::MichelDecomp::MC_NCTotalComponent(), ana::MichelDecomp::MC_NueComponent(), ana::MichelDecomp::MC_NumuComponent(), ana::MichelDecomp::NCTotalComponent(), ana::MichelDecomp::NueComponent(), ana::MichelDecomp::NumuComponent(), ana::Spectrum::ToTH1(), and registry_explorer::v.

Referenced by MakeDCMPPlots().

{
  std::vector<TH1*> v;
  double ndpot = 8.09E20;
  double max = 0;
  std::cout << ndpot << std::endl;
  if (s == "cc") {
      v.push_back((TH1*)(mb->MC_NumuComponent() +
                         mb->MC_AntiNumuComponent()).ToTH1(ndpot));
      v.push_back((TH1*)(mb->NumuComponent() +
                         mb->AntiNumuComponent()).ToTH1(ndpot));
  }
  if (s == "nc") {
    v.push_back((TH1*)mb->MC_NCTotalComponent().ToTH1(ndpot));
    v.push_back((TH1*)mb->NCTotalComponent().ToTH1(ndpot));
  }
  if (s == "ne") {
    v.push_back((TH1*)(mb->MC_NueComponent() +
                       mb->MC_AntiNueComponent()).ToTH1(ndpot));
    v.push_back((TH1*)(mb->NueComponent() +
                       mb->AntiNueComponent()).ToTH1(ndpot));
  }
  std::cout << v.size() << std::endl;
  v.back()->SetLineWidth(2);
  std::cout << v.back()->Integral() << std::endl;
  if (v.back()->GetMaximum() > max) max = v.back()->GetMaximum();

  v[0]->SetMaximum(1.2*max);

  v[0]->SetLineColor(kRed);
  v[1]->SetLineColor(kGreen+2);

  v[0]->Scale(0.001); // Scale to a manageable y-axis range
  v[1]->Scale(0.001);

  for (int i = 0; i < 2; i++){
    v[i]->GetXaxis()->SetTitle("Reconstructed Energy [GeV]");
    v[i]->GetXaxis()->CenterTitle();
    v[i]->GetYaxis()->SetTitle("10^{3} Events / 8.09e20");
    v[i]->GetYaxis()->CenterTitle();
    v[i]->SetTitle("");
  }
  return v;
}

void PIDBinLabels

(

std::string

pid

)

Definition at line 81 of file draw_decomp_plots.C.

References ana::assert(), ana::bins, ana::Binning::Edges(), genie::utils::style::Format(), MECModelEnuComparisons::i, ana::kCVNNLBinning, ana::kLEMNLBinning, ana::kLIDNLBinning, pot, and submit_syst::y.

{
  gPad->Update();
  const double y = gPad->GetUymax()*0.94;
  const ana::Binning * bins =0;
  if (pid=="LID") bins = & ana::kLIDNLBinning;
  if (pid=="LEM") bins = & ana::kLEMNLBinning;
  if (pid=="CVN") bins = & ana::kCVNNLBinning;
  assert (bins);

  double firstwidth;

  for (int i=0; i<3; i++){
    TLatex* ltx = new TLatex(3+6*i,y,TString::Format("%g < %s < %g",
                                                     bins->Edges()[i],
                                                     pid.c_str(),
                                                     bins->Edges()[i+1]));
    ltx->SetTextSize(0.035);
    ltx->SetTextAlign(22);
    ltx->Draw();
    if (i==0) firstwidth=ltx->GetXsize();
  }
  TLatex* pot = new TLatex(3-firstwidth/2,gPad->GetUymax()*0.84,"3.72#times10^{20} POT");
  pot->SetTextSize(0.035);
  pot->SetTextAlign(12);
  pot->Draw();
}

void PIDDivisions

(

)

Definition at line 66 of file draw_decomp_plots.C.

References submit_hadd::l, and central_limit::l1.

{
  TGraph * l1 = new TGraph;
  l1->SetPoint(0, 6, 0);
  l1->SetPoint(1, 6, 50000);
  TGraph * l2 = new TGraph;
  l2->SetPoint(0, 12, 0);
  l2->SetPoint(1, 12, 50000);
  for (TGraph* l: {l1,l2}){
    l->SetLineWidth(1);
    l->Draw("l same");
  }
}

void PlotComp	(	std::vector< TH1 * >	hs,
		std::string	label
	)

Definition at line 174 of file draw_decomp_plots.C.

References allTimeWatchdog::can, plot_validation_datamc::Clone(), h1, h2, h3, and MECModelEnuComparisons::leg.

{
  TH1* h1 = (TH1*)hs[1]->Clone("h1");
  TH1* h2 = (TH1*)hs[2]->Clone("h2");
  TH1* h3 = (TH1*)hs[3]->Clone("h3");

  h1->Divide(hs[0]);
  h2->Divide(hs[0]);
  h3->Divide(hs[0]);

  h1->GetYaxis()->SetTitle("Scale Factor");

  gStyle->SetPadLeftMargin(0.15);

  NuePlotStyle style(label.c_str());
  TCanvas* can = style.MakeSACanvas(hs);

  h1->Draw("same,hist");
  h2->Draw("same,hist");
  h3->Draw("same,hist");

  TLegend *leg = new TLegend(0.28,0.68,0.48,0.88);
  leg->AddEntry(h1, "Prop / MC", "l");
  leg->AddEntry(h2, "Michel / MC", "l");
  leg->AddEntry(h3, "Combo / MC", "l");
  leg->Draw();

  can->SaveAs(("MDCMP_rats_"+label+".png").c_str());
}

void PlotNME

(

TFile *

in

)

Definition at line 432 of file draw_decomp_plots.C.

References allTimeWatchdog::can, cc(), FillWithDimmerColor(), MECModelEnuComparisons::i, ana::kBeamNueBackgroundColor, ana::kNCBackgroundColor, ana::kNumuBackgroundColor, MECModelEnuComparisons::leg, nc, and ana::UniqueName().

Referenced by draw_decomp_plots().

{
  TCanvas *can = new TCanvas();

  TH2F *da2D = (TH2F*)in->Get("MichelDecomp_Plots/templates/DataTemplate");
  TH2F *cc2D = (TH2F*)in->Get("MichelDecomp_Plots/templates/NumuTemplate");
  TH2F *nc2D = (TH2F*)in->Get("MichelDecomp_Plots/templates/NCTemplate");
  TH2F *ne2D = (TH2F*)in->Get("MichelDecomp_Plots/templates/NueTemplate");

  TH1D *da = (TH1D*)da2D->ProjectionY();
  TH1D *cc = (TH1D*)cc2D->ProjectionY();
  TH1D *nc = (TH1D*)nc2D->ProjectionY();
  TH1D *ne = (TH1D*)ne2D->ProjectionY();

  da->Scale(0.001);
  cc->Scale(0.001);
  ne->Scale(0.001);
  nc->Scale(0.001);

  THStack *st = new THStack(UniqueName().c_str(),"");

  ne->SetLineColor(kBeamNueBackgroundColor);
  ne->SetLineWidth(2);
  ne->SetFillStyle(1001);
  FillWithDimmerColor(ne,0.65,true);
  st->Add(ne);

  nc->SetLineColor(kNCBackgroundColor);
  nc->SetLineWidth(2);
  nc->SetFillStyle(1001);
  FillWithDimmerColor(nc,0.25,true);
  st->Add(nc);

  cc->SetLineColor(kNumuBackgroundColor);
  cc->SetLineWidth(2);
  cc->SetFillStyle(1001);
  FillWithDimmerColor(cc,0.10,true);
  st->Add(cc);


  da->GetXaxis()->SetTitle("Number of Michel Electrons");
  da->GetXaxis()->CenterTitle();
  da->GetYaxis()->SetTitle("10^{3} Events / 8.09#times10^{20}");
  da->GetYaxis()->CenterTitle();
  da->SetMarkerStyle(8);
  da->GetYaxis()->SetRangeUser(0,27);

  for (int i = 1; i <= 3; i++)
    da->SetBinError(i,da->GetBinContent(i)/1000);

  da->Draw("E1");
  st->Draw("same");

  TLegend *leg = new TLegend(0.5,0.54,0.8,0.78);
  leg->AddEntry(da, "ND Data","lep");
  leg->AddEntry(nc, "NC", "f");
  leg->AddEntry(cc, "#nu_{#mu} CC", "f");
  leg->AddEntry(ne, "Beam #nu_{e}", "f");
  leg->Draw();

  can->SaveAs("Nue2017_NME.pdf");
  can->SaveAs("Nue2017_NME.png");
  can->SaveAs("Nue2017_NME.eps");
  can->SaveAs("Nue2017_NME.gif");
}

void PlotRat	(	std::vector< TH1 * >	hs,
		std::string	label
	)

Definition at line 204 of file draw_decomp_plots.C.

References allTimeWatchdog::can, make_mec_shifts_plots::GetMaximum(), MECModelEnuComparisons::i, makeTrainCVSamples::int, MECModelEnuComparisons::leg, and cet::sqlite::max().

{
  double max = 0;
  for (int i = 0; i < (int)hs.size(); i++)
    if (hs[i]->GetMaximum() > max) max = hs[i]->GetMaximum();

  hs[0]->GetYaxis()->SetRangeUser(0,1.2*max);

  gStyle->SetPadLeftMargin(0.15);

  NuePlotStyle style(label.c_str());
  TCanvas* can = style.MakeSACanvas(hs);

  hs[0]->Draw("hist,same");
  hs[1]->Draw("same,hist");
  hs[2]->Draw("same,hist");
  hs[3]->Draw("same,hist");

  TLegend *leg = new TLegend(0.4,0.68,0.6,0.88);
  leg->AddEntry(hs[0], "MC", "l");
  leg->AddEntry(hs[1], "Prop", "l");
  leg->AddEntry(hs[2], "Michel", "l");
  leg->AddEntry(hs[3], "Combo", "l");
  leg->Draw();

  can->SaveAs(("MDCMP_hist_"+label+".png").c_str());
}

void PlotRatio	(	std::vector< TH1 * >	ccs,
		std::vector< TH1 * >	ncs,
		std::vector< TH1 * >	nes
	)

Definition at line 250 of file draw_decomp_plots.C.

References file_size_ana::axes, allTimeWatchdog::can, plot_validation_datamc::Clone(), DrawRatioAnaBin(), ana::kBeamNueBackgroundColor, ana::kNCBackgroundColor, ana::kNumuBackgroundColor, MECModelEnuComparisons::leg, ana::Nue2017FourBinAxis(), ana::Nue2017FourBinDivisions(), ana::Nue2017FourBinLabels(), one(), and ana::UniqueName().

Referenced by MakeDCMPPlots().

{
  TH1F* ccratio = (TH1F*) ccs[1]->Clone("ccratio");
  ccratio->Divide(ccs[0]);
  ccratio->SetLineColor(kNumuBackgroundColor);
  ccratio->SetFillStyle(0);
  ccratio->GetYaxis()->SetTitle("Ratio");
  TH1F* ncratio = (TH1F*) ncs[1]->Clone("ncratio");
  ncratio->Divide(ncs[0]);
  ncratio->SetLineColor(kNCBackgroundColor);
  ncratio->SetFillStyle(0);
  TH1F* neratio = (TH1F*) nes[1]->Clone("neratio");
  neratio->Divide(nes[0]);
  neratio->SetLineColor(kBeamNueBackgroundColor);
  neratio->SetFillStyle(0);

  TLine* one = new TLine(0,1,27,1);
  one->SetLineStyle(2);

  TLegend *leg = new TLegend(0.41,0.62,0.61,0.82);
  leg->SetFillStyle(0);
  leg->AddEntry(neratio, "Beam #nu_{e} / MC", "l");
  leg->AddEntry(ccratio, "#nu_{#mu} CC / MC", "l");
  leg->AddEntry(ncratio, "NC / MC", "l");

  TCanvas* can = new TCanvas();
  TH1 *axes = new TH1F(UniqueName().c_str(),"",27,0,27);
  // Explicitly cut out ratio of 15ish in first bin
  axes->GetYaxis()->SetRangeUser(0,2.2);
  axes->GetYaxis()->SetTitle("Decomposition Scale");
  axes->GetYaxis()->CenterTitle();
  axes->GetXaxis()->SetTitle("Reconstructed neutrino energy (GeV)");
  axes->GetXaxis()->CenterTitle();
  //axes->GetYaxis()->SetTitleOffset(1.2);
  axes->Draw();

  Nue2017FourBinLabels(0.98,0.028,kGray+3,false);
  Nue2017FourBinAxis(axes,true,false);
  Nue2017FourBinDivisions(kGray+3,3);

  // The data / MC ratio is about 16 in the first bin.
  // Since there's no events down there, let's just suppress the histogram
  // for that bin.
  DrawRatioAnaBin(ccratio,ncratio,neratio,0);
  DrawRatioAnaBin(ccratio,ncratio,neratio,1);
  DrawRatioAnaBin(ccratio,ncratio,neratio,2);


  //Preliminary(); //uncomment this line if plot is preliminary
  gPad->SetFillStyle(0);
  /*
  ccratio->Draw("][,hist,same");
  ncratio->Draw("][,hist,same");
  neratio->Draw("][,hist,same");
  */
  one->Draw();
  leg->Draw();

  can->SaveAs("Nue2017_MDCMP_ratio.pdf");
  can->SaveAs("Nue2017_MDCMP_ratio.png");
  can->SaveAs("Nue2017_MDCMP_ratio.eps");
  can->SaveAs("Nue2017_MDCMP_ratio.gif");

}

void PlotStack	(	std::vector< TH1 * >	ccs,
		std::vector< TH1 * >	ncs,
		std::vector< TH1 * >	nes,
		bool	PlotMC = true,
		bool	PlotDecomp = true
	)

Definition at line 318 of file draw_decomp_plots.C.

References file_size_ana::axes, allTimeWatchdog::can, plot_validation_datamc::Clone(), e, FillWithDimmerColor(), ana::kBeamNueBackgroundColor, ana::kNCBackgroundColor, ana::kNumuBackgroundColor, MECModelEnuComparisons::leg, cet::sqlite::max(), std::max(), ana::Nue2017FourBinAxis(), ana::Nue2017FourBinDivisions(), ana::Nue2017FourBinLabels(), string, and ana::UniqueName().

{

  TH1D* test = new TH1D (UniqueName().c_str(),"",10,0,1);
  test->SetLineStyle(2);
  test->SetLineColor(kBlack);
  test->SetLineWidth(2);

  TH1F * data = (TH1F*)ccs[1]->Clone(UniqueName().c_str());
  data->Add(ncs[1]);
  data->Add(nes[1]);
  data->SetLineColor(kBlack);
  data->SetMarkerStyle(kFullCircle);
  data->SetMarkerSize(1);



  THStack *combostack = new THStack(UniqueName().c_str(),"");
  nes[1]->SetLineColor(kBeamNueBackgroundColor);
  nes[1]->SetLineWidth(2);
  nes[1]->SetFillStyle(1001);
  FillWithDimmerColor(nes[1],0.65,true);
  combostack->Add(nes[1]);
  ccs[1]->SetLineColor(kNumuBackgroundColor);
  ccs[1]->SetLineWidth(2);
  ccs[1]->SetFillStyle(1001);
  FillWithDimmerColor(ccs[1],0.25, true);
  combostack->Add(ccs[1]);
  ncs[1]->SetLineColor(kNCBackgroundColor);
  ncs[1]->SetLineWidth(2);
  ncs[1]->SetFillStyle(1001);
  FillWithDimmerColor(ncs[1], 0.10, true);
  combostack->Add(ncs[1]);

  THStack *mcstack = new THStack(UniqueName().c_str(),"");
  nes[0]->SetLineColor(kBeamNueBackgroundColor);
  nes[0]->SetLineStyle(2);
  FillWithDimmerColor(nes[0],1e-6,true);
  if (PlotMC && !PlotDecomp){
    nes[0]->SetLineStyle(1);
    FillWithDimmerColor(nes[0],0.65,true);
    nes[0]->SetFillStyle(1001);
  }
  nes[0]->SetLineWidth(2);
  mcstack->Add(nes[0]);
  ccs[0]->SetLineColor(kNumuBackgroundColor);
  ccs[0]->SetLineStyle(2);
  ccs[0]->SetFillStyle(0);
  FillWithDimmerColor(ccs[0],1e-6,true);
  if (PlotMC && !PlotDecomp){
    ccs[0]->SetLineStyle(1);
    FillWithDimmerColor(ccs[0],0.25,true);
    ccs[0]->SetFillStyle(1001);
  }
  ccs[0]->SetLineWidth(2);
  mcstack->Add(ccs[0]);
  ncs[0]->SetLineColor(kNCBackgroundColor);
  ncs[0]->SetLineStyle(2);
  ncs[0]->SetFillStyle(0);
  FillWithDimmerColor(ncs[0],1e-6,true);
  if (PlotMC && !PlotDecomp){
    ncs[0]->SetLineStyle(1);
    FillWithDimmerColor(ncs[0],0.10,true);
    ncs[0]->SetFillStyle(1001);
  }
  ncs[0]->SetLineWidth(2);
  mcstack->Add(ncs[0]);

  TCanvas* can = new TCanvas();

  TH1 *axes = new TH1F(UniqueName().c_str(),"",27,0,27);
  // Explicitly cut out ratio of 15ish in first bin
  double max= std::max(mcstack->GetMaximum(),combostack->GetMaximum());
  axes->GetYaxis()->SetRangeUser(0,1.5*max);
  axes->GetYaxis()->SetTitle("10^{3} Events / 8.09#times10^{20}");
  axes->GetYaxis()->CenterTitle();
  axes->GetXaxis()->SetTitle("Reconstructed neutrino energy (GeV)");
  axes->GetXaxis()->CenterTitle();
  axes->Draw();

  Nue2017FourBinLabels(0.98,0.028,kGray+3,false);
  Nue2017FourBinAxis(axes,true,false);
  Nue2017FourBinDivisions(kGray+3,3);


  can->SetFrameFillColor(0);
  can->SetFrameFillStyle(4000);
  gPad->SetFillStyle(0);

  if (PlotMC) mcstack->Draw("noclear,hist,same");
  if (PlotDecomp) combostack->Draw("noclear,hist,same");
  data->Draw("E1,same");

  TLegend *leg = new TLegend(0.41,0.53,0.61,0.83);
  leg->SetTextSize(0.038);
  leg->AddEntry(data, "ND Data","lep");
  leg->AddEntry(ncs[1], "NC", "f");
  leg->AddEntry(ccs[1], "#nu_{#mu} CC", "f");
  leg->AddEntry(nes[1], "Beam #nu_{e}", "f");
  if (PlotDecomp && PlotMC) leg->AddEntry(test, "Uncorrected MC", "l");
  leg->Draw();

  std::string appen = "";
  if ( PlotMC && !PlotDecomp) appen = "_MCOnly";
  if (!PlotMC &&  PlotDecomp) appen = "_DecompOnly";
  can->SaveAs(("Nue2017_MDCMP_stack"+appen+".pdf").c_str());
  can->SaveAs(("Nue2017_MDCMP_stack"+appen+".png").c_str());
  can->SaveAs(("Nue2017_MDCMP_stack"+appen+".eps").c_str());
  can->SaveAs(("Nue2017_MDCMP_stack"+appen+".gif").c_str());
}

void Preliminary

(

)

Definition at line 55 of file draw_decomp_plots.C.

Referenced by accum_nue_numu_equivalent(), accum_nue_numu_equivalent_fhc(), accum_nue_numu_equivalent_rhc(), bicount(), calcAlphaBetaEachBin(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), CompareDecompDataMC(), CompareSlices(), demoFitContours(), Draw1DProfile(), DrawCovMx(), FD_Data_PosComp(), FD_DataOnly_PosComp(), joint_fit_2018_contours(), joint_fit_2019_contours(), main(), make_angle_plot(), make_plots(), MakeCanvasForData(), MakeCanvasForSim(), MakeHBar(), MakeNuEnergyPlots(), MakeRawPlot(), MakeResultPlot(), MakeSplitCans(), makeSystTable(), MatrixValidation(), median_contours(), monoprob(), NCPaperPlot(), ND_DataMC_Comp_Systs(), nova_t2k_comparison(), Nus17vsAn01_NuFACT_2D(), NusContoursPlusMINOS(), Plot1DSlice(), Plot2DSlice(), plot_DataMCComp_numu(), plot_joint_fit_2020_contours(), plot_ND_DataMC(), plot_ND_DataMC_energybinning(), plot_nd_spectra_2018(), plot_recoE_numu(), plot_time(), plotContProf(), plotDataPred(), plotNumuWrongSign(), PlotNZ(), ana::PlotStack(), PlotSyst(), PlotSysts(), plotWrongSignFraction(), PlotXY(), PositionComparison(), print_tables(), ana::BENDecomp::SavePlotsKa(), ana::BENDecomp::SavePlotsPi(), sensitivity_plot(), and Tutorial2019FitContours().

{
  TLatex* prelim = new TLatex(.9, .95, "NOvA Preliminary");
  prelim->SetTextColor(kBlue);
  prelim->SetNDC();
  prelim->SetTextSize(2/30.);
  prelim->SetTextAlign(32);
  prelim->Draw();
}