plot_resolution.C File Reference

#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/LoadFromFile.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH2.h"
#include "TLegend.h"

Go to the source code of this file.

Functions
void	Compare (const std::vector< TH1 * > &hs)
void	reco_minus_true_panels (std::map< std::string, std::map< std::string, Spectrum * >> &plots, bool log)
void	plot_2d (std::map< std::string, std::map< std::string, Spectrum * >> &plots)
void	plot_2d_vars (std::map< std::string, std::map< std::string, Spectrum * >> &plots)
void	plot (std::string label, std::map< std::string, std::map< std::string, Spectrum * >> &plots, bool log)
void	plot_diff (std::map< std::string, std::map< std::string, Spectrum * >> &plots, bool log)
void	plot_resolution ()

Variables
const std::vector< std::string >	ress = {"2"}
const std::vector< int >	cols = {kRed}

Function Documentation

void Compare

(

const std::vector< TH1 * > &

hs

)

Definition at line 29 of file plot_resolution.C.

References make_syst_table_plots::h, and maxy.

Referenced by genie::flux::operator<<(), genie::AxialFormFactor::operator==(), genie::ELFormFactors::operator==(), genie::QELFormFactors::operator==(), genie::DISStructureFunc::operator==(), plot(), plot_diff(), and reco_minus_true_panels().

{
  double maxy = 0;
  TH1* bigh = 0;
  for(TH1* h: hs){
    h->Scale(1/h->Integral());
    if(h->GetBinContent(h->GetMaximumBin()) > maxy){
      maxy = h->GetBinContent(h->GetMaximumBin());
      bigh = h;
    }
  }

  bigh->Draw("hist");
  for(TH1* h: hs) h->Draw("hist same");
}

void plot	(	std::string	label,
		std::map< std::string, std::map< std::string, Spectrum * >> &	plots,
		bool	log
	)

Definition at line 184 of file plot_resolution.C.

References cols, Compare(), e, make_syst_table_plots::h, he, hh, calib::j, MECModelEnuComparisons::leg, plots, ana::Spectrum::POT(), ress, string, and ana::Spectrum::ToTH1().

Referenced by plot_resolution().

{
  new TCanvas;

  TLegend* leg = new TLegend(.6, .6, .85, .85);
  leg->SetFillStyle(0);

  std::vector<TH1*> hs;

  for(unsigned int j = 0; j < ress.size(); ++j){
    const std::string res = ress[j];
    Spectrum* h = plots[res][label+"h"];
    TH1* hh = h->ToTH1(h->POT());
    hh->SetLineColor(cols[j]);
    hh->Scale(1/hh->Integral());
    hh->Draw(j ? "hist same" : "hist");
    hs.push_back(hh);
    if(!log){
      if(label == "dist") hh->GetXaxis()->SetRangeUser(0, 40);
      if(label == "distnorm") hh->GetXaxis()->SetRangeUser(0, 5);
    }
    if(label == "distnorm") hh->GetYaxis()->SetTitle("Events / cm^{2}");

    leg->AddEntry(hh, ("HoughVertex "+res).c_str(), "l");
  }

  Spectrum* e = plots[ress[0]][label+"e"];
  TH1* he = e->ToTH1(e->POT());
  if(!log){
    if(label == "dist") he->GetXaxis()->SetRangeUser(0, 40);
    if(label == "distnorm") he->GetXaxis()->SetRangeUser(0, 5);
  }
  if(label == "distnorm") he->GetYaxis()->SetTitle("Events / cm^{2}");
  he->Scale(1/he->Integral());
  he->Draw("hist same");
  hs.push_back(he);
  leg->AddEntry(he, "ElasticArms", "l");

  Compare(hs);

  leg->Draw();

  if(log){
    gPad->SetLogy();

    gPad->Print(("plots/tot_"+label+"_log.png").c_str());
    gPad->Print(("plots/tot_"+label+"_log.eps").c_str());
  }
  else{
    gPad->Print(("plots/tot_"+label+".png").c_str());
    gPad->Print(("plots/tot_"+label+".eps").c_str());
  }
}

void plot_2d

(

std::map< std::string, std::map< std::string, Spectrum * >> &

plots

)

Definition at line 113 of file plot_resolution.C.

References plot_validation_datamc::c, make_syst_table_plots::h, plots, ana::Spectrum::POT(), ress, string, and ana::Spectrum::ToTH2().

Referenced by plot_resolution().

{

  for(std::string res: ress){
    TCanvas* c = new TCanvas;//("c", "c", 800, 800);
    c->Divide(2, 2);

    c->cd(1);
    Spectrum* xd2d = plots[res]["xd2d"];
    TH2* h = xd2d->ToTH2(xd2d->POT());
    h->GetYaxis()->SetTitle("ElasticArms Reco - True X (cm)");
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(2);
    Spectrum* yd2d = plots[res]["yd2d"];
    h = yd2d->ToTH2(yd2d->POT());
    h->GetYaxis()->SetTitle("ElasticArms Reco - True Y (cm)");
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(3);
    Spectrum* zd2d = plots[res]["zd2d"];
    h = zd2d->ToTH2(zd2d->POT());
    h->GetYaxis()->SetTitle("ElasticArms Reco - True Z (cm)");
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(0);
    gPad->Print(("plots/vtx_compare"+res+".png").c_str());
    gPad->Print(("plots/vtx_compare"+res+".eps").c_str());
  }
}

void plot_2d_vars

(

std::map< std::string, std::map< std::string, Spectrum * >> &

plots

)

Definition at line 147 of file plot_resolution.C.

References plot_validation_datamc::c, make_syst_table_plots::h, plots, ana::Spectrum::POT(), ress, string, and ana::Spectrum::ToTH2().

Referenced by plot_resolution().

{

  for(std::string res: ress){
    TCanvas* c = new TCanvas;//("c", "c", 800, 800);
    c->Divide(2, 2);

    c->cd(1);
    Spectrum* disth2d = plots[res]["disth2d"];
    TH2* h = disth2d->ToTH2(disth2d->POT());
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(2);
    Spectrum* diste2d = plots[res]["diste2d"];
    h = diste2d->ToTH2(diste2d->POT());
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(3);
    Spectrum* calh2d = plots[res]["calh2d"];
    h = calh2d->ToTH2(calh2d->POT());
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(4);
    Spectrum* cale2d = plots[res]["cale2d"];
    h = cale2d->ToTH2(cale2d->POT());
    h->Draw("colz");
    gPad->SetLogz();

    c->cd(0);
    gPad->Print(("plots/vtx_distvar"+res+".png").c_str());
    gPad->Print(("plots/vtx_distvar"+res+".eps").c_str());
  }
}

void plot_diff	(	std::map< std::string, std::map< std::string, Spectrum * >> &	plots,
		bool	log
	)

Definition at line 238 of file plot_resolution.C.

References cols, Compare(), calib::j, MECModelEnuComparisons::leg, plots, ana::Spectrum::POT(), ress, string, and ana::Spectrum::ToTH1().

Referenced by plot_resolution().

{
  new TCanvas;

  TLegend* leg = new TLegend(.6, .6, .85, .85);
  leg->SetFillStyle(0);

  std::vector<TH1*> hs;

  for(unsigned int j = 0; j < ress.size(); ++j){
    const std::string res = ress[j];
    Spectrum* distdiffh = plots[res]["distdiff"];
    TH1* hdistdiffh = distdiffh->ToTH1(distdiffh->POT());
    hdistdiffh->SetLineColor(cols[j]);
    hdistdiffh->Scale(1/hdistdiffh->Integral());
    hdistdiffh->Draw(j ? "hist same" : "hist");
    hs.push_back(hdistdiffh);
    if(!log) hdistdiffh->GetXaxis()->SetRangeUser(-20, +20);

    leg->AddEntry(hdistdiffh, ("HoughVertex "+res).c_str(), "l");
  }

  Compare(hs);

  leg->Draw();

  if(log){
    gPad->SetLogy();

    gPad->Print("plots/dist_diff_log.png");
    gPad->Print("plots/dist_diff_log.eps");
  }
  else{
    gPad->Print("plots/dist_diff.png");
    gPad->Print("plots/dist_diff.eps");
  }
}

void plot_resolution

(

)

Definition at line 276 of file plot_resolution.C.

References plot_validation_datamc::fname, plot(), plot_2d(), plot_2d_vars(), plot_diff(), plots, reco_minus_true_panels(), runNovaSAM::release, ress, and string.

{
  std::map<std::string, std::map<std::string, Spectrum*>> plots;
  
  // Set up to be able to compare scans over HoughVertex parameter
  // SigmaZ
  // Can be easily modified to scan over other parameters

  for(std::string res: ress){//{"0.5", "1", "2"}){
    const std::string fname = "sigma."+res+".nonswap.root";
    
    for(std::string name: {"xdh", "ydh", "zdh", "xde", "yde", "zde",
          "xd2d", "yd2d", "zd2d", "disth", "diste", "disth2d", "diste2d", "calh2d", "cale2d", "distnormh", "distnorme", "distdiff"}){
      plots[res][name] = LoadFromFile<Spectrum>(fname, name).release();
    }
  }

  reco_minus_true_panels(plots, false);
  reco_minus_true_panels(plots, true);

  plot_2d(plots);
  plot_2d_vars(plots);

  plot("dist", plots, false);
  plot("dist", plots, true);

  plot("distnorm", plots, false);
  plot("distnorm", plots, true);

  plot_diff(plots, false);
  plot_diff(plots, true);
}

void reco_minus_true_panels	(	std::map< std::string, std::map< std::string, Spectrum * >> &	plots,
		bool	log
	)

Definition at line 45 of file plot_resolution.C.

References plot_validation_datamc::c, cols, Compare(), e, make_syst_table_plots::h, he, hh, MECModelEnuComparisons::i, calib::j, label, PandAna.Demos.pi0_spectra::labels, MECModelEnuComparisons::leg, plots, ana::Spectrum::POT(), ress, string, and ana::Spectrum::ToTH1().

Referenced by plot_resolution().

{
  TCanvas* c = new TCanvas;
  c->Divide(2, 2);

  const std::vector<std::string> labels = {"xd", "yd", "zd"};

  TLegend* leg = new TLegend(.1, .2, .9, .8);

  for(int i = 0; i < 3; ++i){
    std::string label = labels[i];

    c->cd(i+1);

    std::vector<TH1*> hs;

    for(unsigned int j = 0; j < ress.size(); ++j){
      const std::string res = ress[j];

      Spectrum* h = plots[res][label+"h"];
      TH1* hh = h->ToTH1(h->POT());
      hh->Scale(1/hh->Integral());
      hh->SetLineColor(cols[j]);
      hs.push_back(hh);
      hh->Draw(j ? "hist same" : "hist");
      if(!log) hh->GetXaxis()->SetRangeUser(-20, +20);

      if(i == 0){
        leg->AddEntry(hh, ("HoughVertex "+res).c_str(), "l");
      }
    }

    Spectrum* e = plots[ress[0]][label+"e"];
    TH1* he = e->ToTH1(e->POT());
    if(!log) he->GetXaxis()->SetRangeUser(-20, +20);
    hs.push_back(he);
    he->Scale(1/he->Integral());
    he->Draw("hist same");

    Compare(hs);

    if(i == 0){
      leg->AddEntry(he, "ElasticArms", "l");
    }
  }

  c->cd(4);
  leg->Draw();

  if(log){
    c->cd(1);
    gPad->SetLogy();
    c->cd(2);
    gPad->SetLogy();
    c->cd(3);
    gPad->SetLogy();
    c->cd(0);

    gPad->Print("plots/vtx_res_log.png");
    gPad->Print("plots/vtx_res_log.eps");
  }
  else{
    c->cd(0);
    gPad->Print("plots/vtx_res.png");
    gPad->Print("plots/vtx_res.eps");
  }
}

Variable Documentation

const std::vector<int> cols = {kRed}

Definition at line 26 of file plot_resolution.C.

Referenced by plot(), plot_diff(), and reco_minus_true_panels().

const std::vector<std::string> ress = {"2"}

Definition at line 25 of file plot_resolution.C.

Referenced by plot(), plot_2d(), plot_2d_vars(), plot_diff(), plot_resolution(), and reco_minus_true_panels().