5         TFile *
f = TFile::Open(inFile);
 
   28         TH1F *
nc=(TH1F *)f->Get(
"nc");
 
   35         TH1F *
hvzpos=(TH1F *)f->Get(
"hvzpos");
 
   36         TH2F *
hvpos=(TH2F *)f->Get(
"hvpos");
 
   46         double eff1=n_etac_4c->GetBinContent(1)/n_events->GetBinContent(1);
 
   47         std::cout<<
"Efficiency (4C-fit) = "<<eff1<<std::endl;
 
   48         double eff2=n_etac_vtx->GetBinContent(1)/n_events->GetBinContent(1);
 
   49         std::cout<<
"Efficiency (vertex fit) = "<<eff2<<std::endl;
 
   51         double eff3=n_etac_vtx_2->GetBinContent(1)/n_events->GetBinContent(1);
 
   52         std::cout<<
"Efficiency (vertex fit, prefit selection) = "<<eff3<<std::endl;
 
   56         std::cout<<
"Max efficiency="<<nc->Integral(5,20)/nc->GetEntries()<<std::endl;
 
   58         TCanvas *
c1=
new TCanvas(
"c1",
"N charged",600,600);
 
   60         if (saveHistos) c1->SaveAs(
"n_charged_stt.png");
 
   62         TCanvas *
c2=
new TCanvas(
"c2",
"No cuts",600,600);
 
   65         h_mphi_nocuts->SetTitleSize(18);
 
   66         h_mphi_nocuts->GetXaxis()->SetTitleSize(0.05);
 
   67         h_mphi_nocuts->GetXaxis()->SetLabelSize(0.06);
 
   68         h_mphi_nocuts->GetXaxis()->SetTitleSize(0.05);
 
   69         h_mphi_nocuts->GetYaxis()->SetLabelSize(0.06);
 
   70         h_mphi_nocuts->Draw();
 
   71         TLine *l1=
new TLine(1.12,0,1.12,35000);
 
   77         h_etac_nocut->SetTitleSize(18);
 
   78         h_etac_nocut->GetXaxis()->SetTitleSize(0.05);
 
   79         h_etac_nocut->GetXaxis()->SetLabelSize(0.06);
 
   80         h_etac_nocut->GetXaxis()->SetTitleSize(0.05);
 
   81         h_etac_nocut->GetYaxis()->SetLabelSize(0.06);
 
   83         if (saveHistos) c2->SaveAs(
"m_nocuts_stt.png");
 
   85         TCanvas *
c3=
new TCanvas(
"c3",
"MC PID",600,600);
 
   88         h_mphi_pid->SetTitleSize(18);
 
   89         h_mphi_pid->GetXaxis()->SetTitleSize(0.05);
 
   90         h_mphi_pid->GetXaxis()->SetLabelSize(0.06);
 
   91         h_mphi_pid->GetXaxis()->SetTitleSize(0.05);
 
   92         h_mphi_pid->GetYaxis()->SetLabelSize(0.06);
 
   95         TLine *l2=
new TLine(1.12,0,1.12,35000);
 
  101         h_etac_pid->SetTitleSize(18);
 
  102         h_etac_pid->GetXaxis()->SetTitleSize(0.05);
 
  103         h_etac_pid->GetXaxis()->SetLabelSize(0.06);
 
  104         h_etac_pid->GetXaxis()->SetTitleSize(0.05);
 
  105         h_etac_pid->GetYaxis()->SetLabelSize(0.06);
 
  107         if (saveHistos) c3->SaveAs(
"m_pid_stt.png");
 
  127         TCanvas *
c6=
new TCanvas(
"c6",
"Vertex position",600,600);
 
  133         if (saveHistos) c6->SaveAs(
"vertex_pos_stt.png");
 
  143         TCanvas *
c8=
new TCanvas(
"c8",
"Vertex resolution",600,600);
 
  148         TF1 *f1_vtxx= 
new TF1(
"f1_vtxx",
"gaus",-1.,1.);
 
  149         hvtxresX->Fit(f1_vtxx,
"R",
"",-1,1.);
 
  150         double mean_x=f1_vtxx->GetParameter(1);
 
  151         double sigma_x=f1_vtxx->GetParameter(2);
 
  152         TF1 *f2_vtxx= 
new TF1(
"f2_vtxx",
"gaus",-1.,1.);
 
  153         hvtxresX->Fit(f2_vtxx,
"R",
"",mean_x-1.6*sigma_x,mean_x+1.6*sigma_x);
 
  154         double mean_x2=f2_vtxx->GetParameter(1);
 
  155         double sigma_x2=f2_vtxx->GetParameter(2);
 
  156         std::cout<<
"!!!!!!!!!! vertex x resolution = "<<sigma_x2<<std::endl;
 
  160         TF1 *f1_vtxy= 
new TF1(
"f1_vtxy",
"gaus",-1.,1.);
 
  161         hvtxresY->Fit(f1_vtxy,
"R",
"",-1,1.);
 
  162         double mean_y=f1_vtxy->GetParameter(1);
 
  163         double sigma_y=f1_vtxy->GetParameter(2);
 
  164         TF1 *f2_vtxy= 
new TF1(
"f2_vtxy",
"gaus",-1.,1.);
 
  165         hvtxresY->Fit(f2_vtxy,
"R",
"",mean_y-1.6*sigma_y,mean_y+1.6*sigma_y);
 
  166         double mean_y2=f2_vtxy->GetParameter(1);
 
  167         double sigma_y2=f2_vtxy->GetParameter(2);
 
  168         std::cout<<
"!!!!!!!!!! vertex y resolution = "<<sigma_y2<<std::endl;
 
  172         TF1 *f1_vtxz= 
new TF1(
"f1_vtxz",
"gaus",-1.,1.);
 
  173         hvtxresZ->Fit(f1_vtxz,
"R",
"",-1,1.);
 
  174         double mean_z=f1_vtxz->GetParameter(1);
 
  175         double sigma_z=f1_vtxz->GetParameter(2);
 
  176         TF1 *f2_vtxz= 
new TF1(
"f2_vtxz",
"gaus",-1.,1.);
 
  177         hvtxresZ->Fit(f2_vtxz,
"R",
"",mean_z-1.6*sigma_z,mean_z+1.6*sigma_z);
 
  178         double mean_z2=f2_vtxz->GetParameter(1);
 
  179         double sigma_z2=f2_vtxz->GetParameter(2);
 
  180         std::cout<<
"!!!!!!!!!! vertex z resolution = "<<sigma_z2<<std::endl;
 
  183         if (saveHistos) c8->SaveAs(
"vertex_res_stt.png");
 
  185         TCanvas *
c9=
new TCanvas(
"c9",
"m vertex",600,600);
 
  188         h_mphi_vtx->SetTitleSize(18);
 
  189         h_mphi_vtx->GetXaxis()->SetTitleSize(0.05);
 
  190         h_mphi_vtx->GetXaxis()->SetLabelSize(0.06);
 
  191         h_mphi_vtx->GetXaxis()->SetTitleSize(0.05);
 
  192         h_mphi_vtx->GetYaxis()->SetLabelSize(0.06);
 
  195         h_etac_vtx->SetTitleSize(18);
 
  196         h_etac_vtx->GetXaxis()->SetTitleSize(0.05);
 
  197         h_etac_vtx->GetXaxis()->SetLabelSize(0.06);
 
  198         h_etac_vtx->GetXaxis()->SetTitleSize(0.05);
 
  199         h_etac_vtx->GetYaxis()->SetLabelSize(0.06);
 
  201         if (saveHistos) c9->SaveAs(
"m_vtx_stt.png");
 
  204         double mean_phi, range_phi, sigma1_phi, mean1_phi, sigma2_phi, mean_etac, range_etac;
 
  205         double sigma1, mean1, sigma2;
 
  249         TCanvas *
c11=
new TCanvas(
"c11",
"Mass final (Vertex fit)",600,600);
 
  252         h_mphi_final_vtx->SetTitleSize(18);
 
  253         h_mphi_final_vtx->GetXaxis()->SetTitleSize(0.05);
 
  254         h_mphi_final_vtx->GetXaxis()->SetLabelSize(0.06);
 
  255         h_mphi_final_vtx->GetXaxis()->SetTitleSize(0.05);
 
  256         h_mphi_final_vtx->GetYaxis()->SetLabelSize(0.06);
 
  258         h_mphi_final_vtx->Draw();
 
  263         TF1 *f1_phi_vtx = 
new TF1(
"f1_phi_vtx",
"gaus",0.9,1.1);
 
  264         h_mphi_final_vtx->Fit(f1_phi_vtx,
"R",
"",mean_phi-range_phi,mean_phi+range_phi);
 
  266         sigma1_phi=f1_phi_vtx->GetParameter(2);
 
  267         mean1_phi=f1_phi_vtx->GetParameter(1);
 
  269         TF1 *f2_phi_vtx = 
new TF1(
"f2_phi_vtx",
"gaus",0.9,1.1);
 
  270         h_mphi_final_vtx->Fit(f2_phi_vtx,
"R",
"",mean1_phi-1.6*sigma1_phi,mean1_phi+1.6*sigma1_phi);
 
  272         sigma2_phi=f2_phi_vtx->GetParameter(2);
 
  273         std::cout<<
"!!!!!!!!!!!! sigma phi (vertex fit)="<<sigma2_phi<<std::endl;
 
  275         TLine *l3=
new TLine(1.0,0,1.0,7000);
 
  280         TLine *l4=
new TLine(1.04,0,1.04,7000);
 
  287         h_etac_phimass_vtx->SetTitleSize(18);
 
  288         h_etac_phimass_vtx->GetXaxis()->SetTitleSize(0.05);
 
  289         h_etac_phimass_vtx->GetXaxis()->SetLabelSize(0.06);
 
  290         h_etac_phimass_vtx->GetXaxis()->SetTitleSize(0.05);
 
  291         h_etac_phimass_vtx->GetYaxis()->SetLabelSize(0.06);
 
  293         h_etac_phimass_vtx->Draw();
 
  298         TF1 *f1_vtx = 
new TF1(
"f1_vtx",
"gaus",2.8,3.2);
 
  299         h_etac_phimass_vtx->Fit(f1_vtx,
"R",
"",mean_etac-range_etac,mean_etac+range_etac);
 
  301         sigma1=f1_vtx->GetParameter(2);
 
  302         mean1=f1_vtx->GetParameter(1);
 
  304         TF1 *f2_vtx = 
new TF1(
"f2_vtx",
"gaus",2.8,3.2);
 
  305         h_etac_phimass_vtx->Fit(f2_vtx,
"R",
"",mean1-1.6*sigma1,mean1+1.6*sigma1);
 
  307         sigma2=f2_vtx->GetParameter(2);
 
  308         std::cout<<
"!!!!!!!!!!!!!!! sigma eta_c (vertex fit) ="<<sigma2<<std::endl;
 
  310         TLine *l5=
new TLine(2.9,0,2.9,1200);
 
  315         TLine *l6=
new TLine(3.06,0,3.06,1200);
 
  320         if (saveHistos) c11->SaveAs(
"m_final_vtx_stt.png");
 
TH1F * h_mphi_final_vtx_2
TH1F * h_etac_phimass_vtx
TH1F * h_etac_phimass_vtx_2