aug14_track_eff_ptmom_count_barrel.cxx File Reference

#include "TCut.h"
#include "cut_pipijpsi.h"
#include "TTree.h"
#include "TGraph.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH1.h"
#include "TString.h"
#include "TText.h"
#include "TGraphErrors.h"
#include <iostream>

Go to the source code of this file.

Functions
void	f_eff (TTree *m_tree, TCut &cut_all, TCut &cut_four, Double_t *m_eff, Double_t *m_err, Double_t *m_xcen, Double_t *m_xerr)
void	aug14_track_eff_ptmom_count_barrel (TString file_fir_name="../../651_pipijpsi_inclusive.root", TString file_sec_name="../../651_pipijpsi_inclusive.root", TCut m_cut_run="", TCut m_cut_run="", TString tree_name="infmom", Int_t m_bins=7)

Variables
TCut	m_cut_three_yes = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && Npi0<1 && maxene<0.3 && epratio<0.9 && (((eveflag==1||eveflag==6) && missmass<0.25) || eveflag==4)"
TCut	m_cut_four_yes = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && Npi0<1 && maxene<0.3 && epratio<0.9 && eveflag==4"
TCut	m_cut_three_no = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && epratio<1.0 && (eveflag==1||eveflag==6 || eveflag==4)"
TCut	m_cut_four_no = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && epratio<1.0 && eveflag==4"

Function Documentation

aug14_track_eff_ptmom_count_barrel()

void aug14_track_eff_ptmom_count_barrel	(	TString	file_fir_name = "../../651_pipijpsi_inclusive.root",
		TString	file_sec_name = "../../651_pipijpsi_inclusive.root",
		TCut	m_cut_run = "",
		TCut	m_cut_run = "",
		TString	tree_name = "infmom",
		Int_t	m_bins = 7 )

Definition at line 36 of file aug14_track_eff_ptmom_count_barrel.cxx.

                                                                                                                                                                                                                                               {
 
  gROOT->SetStyle("BES");
  TFile *m_file = new TFile(file_fir_name);
  TTree *m_tree = (TTree*)m_file->Get(tree_name);
 
  Double_t m_data_eff[7], m_data_err[7], m_mc_eff[7], m_mc_err[7];
  Double_t m_xcen[7], m_xerr[7];
  f_eff(m_tree, m_cut_three_yes, m_cut_four_yes, m_data_eff, m_data_err, m_xcen, m_xerr);
 
  TFile *m_file2 = new TFile(file_sec_name);
  TTree *m_tree2 = (TTree*)m_file2->Get(tree_name);
  f_eff(m_tree2, m_cut_three_no, m_cut_four_no, m_mc_eff, m_mc_err, m_xcen, m_xerr);
 
  // create the TGraphErrors and draw it
  TCanvas *c1 = new TCanvas("c1","c1",10,10,800,643);
  TGraphErrors *gr_data = new TGraphErrors(m_bins, m_xcen, m_data_eff, m_xerr, m_data_err);
  gr_data->SetTitle("Tracking effciency of #pi^{-}");
  gr_data->GetXaxis()->SetTitle("transverse momentum (GeV)");
  gr_data->GetYaxis()->SetTitle("efficiency");
  gr_data->GetYaxis()->SetRangeUser(0.65,1.0);
  gr_data->SetMarkerColor(2);
  gr_data->SetMarkerStyle(20);
 
  TGraphErrors *gr_mc = new TGraphErrors(m_bins, m_xcen, m_mc_eff, m_xerr, m_mc_err);
  gr_mc->SetMarkerColor(4);
  gr_mc->SetMarkerStyle(21);
  gr_mc->GetXaxis()->SetTitle("transverse momentum (GeV)");
  gr_mc->GetYaxis()->SetTitle("efficiency");
  gr_data->Draw("AP");
  gr_mc->Draw("P");
  c1->Update();
 
  TText m_text;
  m_text.SetTextSize(0.08);
  m_text.SetTextColor(2);
  m_text.DrawTextNDC(0.5,0.4,"bg 2.4%");
 
  m_text.SetTextColor(4);
  m_text.DrawTextNDC(0.5,0.3,"bg 12%");
 
  m_text.SetTextColor(1);
  m_text.DrawTextNDC(0.5,0.2,"651, Barrel");
  m_text.DrawTextNDC(0.5,0.5,"Inclusive MC");
 
  c1->Print("figs/nov30_track_eff_ptmom_count_barrel_651_bg_study.eps");
  c1->Print("figs/nov30_track_eff_ptmom_count_barrel_651_bg_study.pdf");
}

f_eff()

void f_eff	(	TTree *	m_tree,
		TCut &	cut_all,
		TCut &	cut_four,
		Double_t *	m_eff,
		Double_t *	m_err,
		Double_t *	m_xcen,
		Double_t *	m_xerr )

Definition at line 20 of file aug14_track_eff_ptmom_count_barrel.cxx.

                                                                                                                              {
  TH1F *hpionm = new TH1F("hpionm","hpionm",7,0.05,0.4);
  TH1F *hfourm = new TH1F("hfourm","hfourm",7,0.05,0.4);
  m_tree->Draw("tppionratiom*mompionm>>hpionm",cut_all+"fabs(costhe[1])<0.7"); 
  m_tree->Draw("tppionratiom*mompionm>>hfourm",cut_four+"fabs(costhe[1])<0.7"); 
 
  for(Int_t i=1; i<8; i++){
    m_eff[i-1] = hpionm->GetBinContent(i) ? hfourm->GetBinContent(i)/hpionm->GetBinContent(i) : 0;
    m_err[i-1] = m_eff[i-1]*(hpionm->GetBinContent(i)) ? sqrt((m_eff[i-1])*(1-m_eff[i-1])/((hpionm->GetBinContent(i)))) : 0 ;
    m_xcen[i-1] = hfourm->GetBinCenter(i);
    m_xerr[i-1] = hfourm->GetBinWidth(i)/2;
    cout << "pt mom: " << 0.05+(i-1)*(0.4-0.05)/7 << "   eff: " <<m_eff[i-1] << "   err: " << m_err[i-1] << endl; 
  }
}

Referenced by aug14_track_eff_ptmom_count_barrel().

Variable Documentation

m_cut_four_no

TCut m_cut_four_no = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && epratio<1.0 && eveflag==4"

Definition at line 18 of file aug14_track_eff_ptmom_count_barrel.cxx.

Referenced by aug14_track_eff_ptmom_count_barrel().

m_cut_four_yes

TCut m_cut_four_yes = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && Npi0<1 && maxene<0.3 && epratio<0.9 && eveflag==4"

Definition at line 15 of file aug14_track_eff_ptmom_count_barrel.cxx.

Referenced by aug14_track_eff_ptmom_count_barrel().

m_cut_three_no

TCut m_cut_three_no = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && epratio<1.0 && (eveflag==1||eveflag==6 || eveflag==4)"

Definition at line 17 of file aug14_track_eff_ptmom_count_barrel.cxx.

Referenced by aug14_track_eff_ptmom_count_barrel().

m_cut_three_yes

TCut m_cut_three_yes = "invjpsi>3.0 && invjpsi<3.2 && pipidang<0.98 && cmslepm<1.7 && cmslepp<1.7 && cmslepm>1.4 && cmslepp>1.4 && mompionm<0.5 && mompionp<0.5 && Npi0<1 && maxene<0.3 && epratio<0.9 && (((eveflag==1||eveflag==6) && missmass<0.25) || eveflag==4)"

Definition at line 14 of file aug14_track_eff_ptmom_count_barrel.cxx.

Referenced by aug14_track_eff_ptmom_count_barrel().