CheckRung/check/draw_plots.cxx

Go to the documentation of this file.

/************* draw some plots for validate the dE/dx calibration
Gain (Bhabha) plots:
-- gain constant vs. run #
-- mean dE/dx vs run# (with new constants applied)
-- gain constant vs. wire #
-- mean dE/dx vs wire# (with new constants applied)
-- distribution of new wire gain - old wire gain (to see how stable they are)
  (one 1-D plot of difference for all wires)
***************************/
 
TString old_path = "../../RungEangle/";
const int max_len = 10000;
bool m_debug = false;
 
void draw_plots(){
  gStyle->SetOptStat(0);
  Double_t x[max_len], y[max_len];
  // gain constant vs. run 
  TFile o_file(old_path + "rungain.root");
  TTree *o_tree_run = (TTree*)o_file.Get("runcalib");
  TTree *o_tree_wire = (TTree*)o_file.Get("wiregcalib");
  Int_t o_run, n_run;
  Double_t o_gain, o_wireg, n_mean, n_wireg;  
  o_tree_run->SetBranchAddress("runno", &o_run);
  o_tree_run->SetBranchAddress("rungain", &o_gain);
  Int_t m_len = o_tree_run->GetEntries();
  for(int i=0; i<m_len; i++){
    o_tree_run->GetEntry(i);
    x[i] = o_run;
    y[i] = o_gain;
    if(m_debug) cout << "x " << x[i] << "  y " << y[i] << endl;
  }
  TGraph *gr = new TGraph(m_len, x, y);
  gr->SetTitle("gain constant vs. run");
  TCanvas *c1 = new TCanvas("c1", "c1", 10, 10, 800, 642);
  gr->Draw("a*");
  c1->SaveAs("figs/gain_vs_run.eps");
 
  // mean vs run
  TFile n_file("../rungain.root");
  TTree *n_tree_run = (TTree*)n_file.Get("runcalib");
  n_tree_run->SetBranchAddress("runno", &n_run);
  n_tree_run->SetBranchAddress("runmean", &n_mean);
  if(m_len != n_tree_run->GetEntries()){
    cout << "the num of enetries are different!" << endl;
    cout << "old " << m_len << "  new " << n_tree_run->GetEntries() << endl;
    return 1;
  }
  for(int i=0; i<m_len; i++){
    n_tree_run->GetEntry(i);
    x[i] = n_run;
    y[i] = n_mean;
  }
  gr = new TGraph(m_len, x, y);
  gr->SetTitle("mean dE/dx vs. run (with new constants)");
  c1->Clear();
  gr->Draw("a*");
  c1->SaveAs("figs/mean_vs_run.eps");
 
  // gain constant vs wire
  TFile o_file("../../WireGain/wiregain.root");
  TTree *o_tree_wire = (TTree*)o_file.Get("wiregcalib");
  o_tree_wire->SetBranchAddress("wireg", y);
  o_tree_wire->GetEntry(0);
  m_len = 6796;
  for(int i=0; i<m_len; i++){
    x[i] = i;
    if(m_debug) cout << "i " << i << "   y[i] " << y[i] << endl;
  }
  gr = new TGraph(m_len, x, y);
  gr->SetTitle("gain constant vs. wire");
  c1->Clear();
  gr->Draw("a*");
  c1->SaveAs("figs/wireg_vs_wire.eps");
 
  // mean dE/dx vs wire
  TFile n_file("wiregain.root");
  TTree *n_tree_wire = (TTree*)n_file.Get("wiregcalib");
  n_tree_wire->SetBranchAddress("fitmean", y);
  n_tree_wire->GetEntry(0);
  m_len = 6796;
  for(int i=0; i<m_len; i++){
    x[i] = i;
    if(y[i]<480) y[i]=y[i-1];
    if(m_debug) cout << "i " << i << "   y[i] " << y[i] << endl;
  }
  gr = new TGraph(m_len, x, y);
  gr->SetTitle("dE/dx vs. wire with new constants");
  c1->Clear();
  gr->Draw("a*");
  c1->SaveAs("figs/mean_vs_wire.eps");
 
  n_tree_wire->SetBranchAddress("fitsigma", y);
  n_tree_wire->GetEntry(0);
  m_len = 6796;
  for(int i=0; i<m_len; i++){
    x[i] = i;
    if(m_debug) cout << "i " << i << "   y[i] " << y[i] << endl;
  }
  gr = new TGraph(m_len, x, y);
  gr->SetTitle("dE/dx sigma vs. wire with new constants");
  c1->Clear();
  gr->Draw("a*");
  c1->SaveAs("figs/sigma_vs_wire.eps");
 
 
// Hit-level Rad. Bhabha plots
// -- "surface plot" of ave. per hit dE/dx vs. dca & entrance angle
//   (can show with and without correction applied)
//   -- Since 2-D surfaces can be hard to see clearly, one can also make 1-D plots:
//   cut in slices of dca and make several 1-D plots vs. angle, 
//        cut in slices of angle and plot 1-D vs. dca 
  TFile o_file("../../DocaEangle/docaeangle.root");
  TH2F *o_hist = (TH2F*)o_file.Get("dEdxVsDocaEAng");
  o_hist->SetTitle("dE/dx vs. doca and enta, without new constants");
  c1->Clear();
//  double m_max_dedx = o_hist->GetBinContent(o_hist->GetMaximumBin());
//  double m_min_dedx = o_hist->GetBinContent(o_hist->GetMinimumBin());
//  double m_col_ratio = (550.-m_min_dedx)/(m_max_dedx-m_min_dedx);
//  const UInt_t Number = 3;
//  Double_t Red[Number]    = { 0.0, 0.00, 0.00};
//  Double_t Green[Number]  = { 0.0, 52./255, 0.00};
//  Double_t Blue[Number]   = { 0.0, 114./255, 0.00};
//  Double_t Length[Number] = { 0.00, m_col_ratio, 1.00};
//  Int_t nb=200;
//  TColor::CreateGradientColorTable(Number,Length,Red,Green,Blue,nb);
  gStyle->SetPalette(1,0);
  o_hist->Draw("colz");
  c1->SaveAs("figs/dedx_vs_doca_and_enta_without_cor.eps");
 
  TFile n_file("docaeangle.root");
  TH2F *n_hist = (TH2F*)n_file.Get("dEdxVsDocaEAng");
  n_hist->SetTitle("dE/dx vs. doca and enta, with new constants");
  c1->Clear();
//  m_max_dedx = n_hist->GetBinContent(n_hist->GetMaximumBin());
//  m_min_dedx = n_hist->GetBinContent(n_hist->GetMinimumBin());
//  m_col_ratio = (550.-m_min_dedx)/(m_max_dedx-m_min_dedx);
//  Length[1] = m_col_ratio;
//  TColor::CreateGradientColorTable(Number,Length,Red,Green,Blue,nb);
  n_hist->Draw("colz");
  c1->SaveAs("figs/dedx_vs_doca_and_enta_with_cor.eps");
 
  TFile o_file("dEdxVsCostheta.root");
  TCanvas *on_can = (TCanvas*)o_file.Get("fitmean");
  on_can->SetTitle("dE/dx vs. cos, with/without correction");
  on_can->Draw();
  on_can->SaveAs("figs/dedx_vs_cos_with_without_cor.eps");
 
  on_can = (TCanvas*)o_file.Get("fitsigma");
  on_can->SetTitle("dE/dx sigma vs. cos, with/without correction");
  on_can->Draw();
  on_can->SaveAs("figs/dedx_sigma_vs_cos_with_without_cor.eps");
 
  TFile o_file("dEdxVsCosNP.root");
  on_can = (TCanvas*)o_file.Get("fitmean");
  on_can->SetTitle("dE/dx vs. cos, electron & positron");
  on_can->Draw();
  on_can->SaveAs("figs/dedx_vs_cos_NP.eps");
 
  on_can = (TCanvas*)o_file.Get("fitsigma");
  on_can->SetTitle("dE/dx sigma vs. cos, electron & positron");
  on_can->Draw();
  on_can->SaveAs("figs/dedx_sigma_vs_cos_NP.eps");
 
}