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