draw_plots.cxx File Reference

Go to the source code of this file.

Functions
void	draw_plots ()

Variables
const int	run_begin = 59163
const int	run_end = 62823
const int	max_len = 10000
bool	m_debug = false

Function Documentation

draw_plots()

void draw_plots

(

)

Definition at line 21 of file FirstRung/draw_plots.cxx.

                 {
//  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; 
}

Variable Documentation

m_debug

bool m_debug = false

Definition at line 19 of file FirstRung/draw_plots.cxx.

Referenced by draw_plots().

max_len

const int max_len = 10000

Definition at line 18 of file FirstRung/draw_plots.cxx.

Referenced by draw_plots().

run_begin

const int run_begin = 59163

Definition at line 14 of file FirstRung/draw_plots.cxx.

Referenced by draw_plots().

run_end

const int run_end = 62823

Definition at line 16 of file FirstRung/draw_plots.cxx.

Referenced by draw_plots().