FirstRung/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)
***************************/
//==== 20191129 ====
// add a function to output the number of runs
//==================
 
const int run_begin = 59163;
 
const int run_end = 62823;
 
const int max_len = 10000;
bool m_debug = false;
 
void draw_plots(){
//  const TString old_path = "/afs/ihep.ac.cn/bes3/offline/CalibConst/dedx/6.6.3/";
  const TString old_path = "./";
  const TString file_name = Form("DedxConst_%d-%d_705.root", run_begin, run_end);
  const TString str_eps = Form("figs/gain_vs_run_%d-%d.eps", run_begin, run_end);
  const TString str_eps2 = Form("figs/resol_vs_run_%d-%d.eps", run_begin, run_end);
  const TString str_eps3 = Form("figs/rel_resol_vs_run_%d-%d.eps", run_begin, run_end);
 
  gStyle->SetOptStat(0);
  Double_t x[max_len], y[max_len], yy[max_len], yyy[max_len];
  // gain constant vs. run 
  TFile o_file(old_path + file_name);
  TTree *o_tree_run = (TTree*)o_file.Get("runcalib");
  Int_t o_run, n_run;
  Double_t o_gain, o_resol, o_wireg, n_mean, n_wireg;  
  o_tree_run->SetBranchAddress("runno", &o_run);
  o_tree_run->SetBranchAddress("rungain", &o_gain);
  o_tree_run->SetBranchAddress("runresol", &o_resol);
  Int_t m_len = o_tree_run->GetEntries();
  Int_t m_count(0), p_run(-1);
  for(int i=0; i<m_len; i++){
    o_tree_run->GetEntry(i);
    x[i] = o_run;
    y[i] = o_gain;
    if(o_run != p_run){
      p_run = o_run;
      m_count ++;
    }
    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(str_eps);
 
  c1->Clear();
  for(int i=0; i<m_len; i++){
    o_tree_run->GetEntry(i);
    yy[i] = o_resol;
    if(y[i]) yyy[i] = 100.*(yy[i]/y[i]/550.);
    if(m_debug) cout << "x " << x[i] << "  y " << yy[i] << endl;
  }
  gr = new TGraph(m_len, x, yy);
  gr->SetTitle("resolution vs. run");
  gr->Draw("a*");
  c1->SaveAs(str_eps2);
 
  c1->Clear();
  gr = new TGraph(m_len, x, yyy);
  gr->SetTitle("relative resolution vs. run (%)");
  gr->Draw("a*");
  c1->SaveAs(str_eps3);
 
  cout << "There are total " << m_count << " runs (not entries)." << endl; 
}