draw_separation.cxx

Go to the documentation of this file.

void draw_separation(){
 
  int step = 3;
  if(step == 1){
    TFile *f = new TFile("pion.root");
    TTree *t = (TTree*)f->Get("n103");
    
    TFile *f_save = new TFile("figs/histo.root","recreate");
    TH2F *vec_h[5];
    for(int i=0; i<5; i++){
      vec_h[i] = new TH2F(Form("h%d",i), Form("h%d",i), 100, 0., 1.5, 100, 0., 20.);
    }
    
    for(int j=0; j<5; j++){
      cout << j << " th project" << endl; 
      t->Project(Form("h%d",j), Form("abs(expectid[2]-expectid[%d])/sqrt(sigmaid[2]*sigmaid[2] + sigmaid[%d]*sigmaid[%d]):ptrk", j, j, j));
    }
    
    f_save->cd();
    for(int k=0; k<5; k++) vec_h[k]->Write();
    f_save->Close();
  }
 
  else if(step == 2){
    TFile *f = new TFile("figs/histo.root", "update");
    TProfile *vec_p[5];
    for(int i=0; i<5; i++){
      TH2F *h_temp = (TH2F*)gDirectory->Get(Form("h%d",i));
      vec_p[i] = h_temp->ProfileX();
      vec_p[i]->Write();
      delete h_temp;
    }
    f->Close();
  }
 
  else if(step == 3){
    TFile *f = new TFile("figs/histo.root", "update");
    TProfile *vec_p[5];
    for(int i=0; i<5; i++){
      vec_p[i] = (TProfile*)f->Get(Form("h%d_pfx",i));
    }
    TCanvas *c1 = new TCanvas("c1","separation power", 10, 10, 800, 642);
    vec_p[3]->Draw();
    vec_p[0]->Draw("same");
    vec_p[4]->Draw("same");
    c1->SaveAs("figs/sepa.C");
    f->Close();
  }
 
}