12 #include <dqm/analysis/modules/DQMHistAnalysisMiraBelle.h>
15 #include <dqm/analysis/modules/DQMHistAnalysis.h>
39 DQMHistAnalysisMiraBelleModule::DQMHistAnalysisMiraBelleModule()
43 setDescription(
"Modify and analyze the data quality histograms of MiraBelle");
63 mumu_main =
new TCanvas(
"mumu_main",
"mumu_main", 0, 0, 800, 600);
64 mumu_resolution =
new TCanvas(
"mumu_resolution",
"mumu_resolution", 0, 0, 800, 600);
65 mumu_muon_val =
new TCanvas(
"mumu_muon_val",
"mumu_muon_val", 0, 0, 400, 400);
66 dst_mass =
new TCanvas(
"dst_mass",
"dst_mass", 0, 0, 1200, 400);
67 dst_pi_val =
new TCanvas(
"dst_pi_val",
"dst_pi_val", 0, 0, 800, 400);
68 dst_k_val =
new TCanvas(
"dst_k_val",
"dst_k_val", 0, 0, 800, 400);
78 B2DEBUG(20,
"DQMHistAnalysisMiraBelle: initialized.");
83 B2DEBUG(20,
"DQMHistAnalysisMiraBelle: beginRun called.");
88 B2DEBUG(20,
"DQMHistAnalysisOutputMiraBelle: event called.");
95 TH1* hist_npxd =
findHist(
"PhysicsObjectsMiraBelle/hist_npxd");
96 TH1* hist_nsvd =
findHist(
"PhysicsObjectsMiraBelle/hist_nsvd");
97 TH1* hist_ncdc =
findHist(
"PhysicsObjectsMiraBelle/hist_ncdc");
98 TH1* hist_topdig =
findHist(
"PhysicsObjectsMiraBelle/hist_topdig");
99 TH1* hist_DetPhotonARICH =
findHist(
"PhysicsObjectsMiraBelle/hist_DetPhotonARICH");
100 TH1* hist_klmTotalHits =
findHist(
"PhysicsObjectsMiraBelle/hist_klmTotalHits");
101 TH1* hist_klmClusterLayers =
findHist(
"PhysicsObjectsMiraBelle/hist_klmClusterLayers");
102 TH1* hist_dD0 =
findHist(
"PhysicsObjectsMiraBelle/hist_dD0");
103 TH1* hist_dZ0 =
findHist(
"PhysicsObjectsMiraBelle/hist_dZ0");
104 TH1* hist_dPtcms =
findHist(
"PhysicsObjectsMiraBelle/hist_dPtcms");
105 TH1* hist_Pval =
findHist(
"PhysicsObjectsMiraBelle/hist_Pval");
106 TH1* hist_nExtraCDCHits =
findHist(
"PhysicsObjectsMiraBelle/hist_nExtraCDCHits");
107 TH1* hist_nECLClusters =
findHist(
"PhysicsObjectsMiraBelle/hist_nECLClusters");
108 TH1* hist_muid =
findHist(
"PhysicsObjectsMiraBelle/hist_muid");
109 TH1* hist_inv_p =
findHist(
"PhysicsObjectsMiraBelle/hist_inv_p");
110 TH1* hist_ndf =
findHist(
"PhysicsObjectsMiraBelle/hist_ndf");
111 TH1* hist_D0 =
findHist(
"PhysicsObjectsMiraBelle/hist_D0");
112 TH1* hist_Z0 =
findHist(
"PhysicsObjectsMiraBelle/hist_Z0");
113 TH1* hist_theta =
findHist(
"PhysicsObjectsMiraBelle/hist_theta");
114 TH1* hist_Phi0 =
findHist(
"PhysicsObjectsMiraBelle/hist_Phi0");
115 TH1* hist_Pt =
findHist(
"PhysicsObjectsMiraBelle/hist_Pt");
116 TH1* hist_Mom =
findHist(
"PhysicsObjectsMiraBelle/hist_Mom");
117 TH1* hist_klmTotalBarrelHits =
findHist(
"PhysicsObjectsMiraBelle/hist_klmTotalBarrelHits");
118 TH1* hist_klmTotalEndcapHits =
findHist(
"PhysicsObjectsMiraBelle/hist_klmTotalEndcapHits");
119 TH1* hist_dPhicms =
findHist(
"PhysicsObjectsMiraBelle/hist_dPhicms");
132 mumu_main->cd(9); hist_DetPhotonARICH->Draw();
133 mumu_main->cd(10); hist_klmClusterLayers->Draw();
134 mumu_main->cd(11); hist_nExtraCDCHits->Draw();
135 mumu_main->cd(12); hist_nECLClusters->Draw();
154 float mean_npxd = hist_npxd->GetMean();
155 float mean_nsvd = hist_nsvd->GetMean();
156 float mean_ncdc = hist_ncdc->GetMean();
157 float mean_nklmlayer = hist_klmClusterLayers->GetMean();
158 float mean_nklm = hist_klmTotalHits->GetMean();
159 float mean_nbklm = hist_klmTotalBarrelHits->GetMean();
160 float mean_neklm = hist_klmTotalEndcapHits->GetMean();
161 float mean_topdig = hist_topdig->GetMean();
162 float mean_parich = hist_DetPhotonARICH->GetMean();
163 float mean_ncdc_ex = hist_nExtraCDCHits->GetMean();
164 float mean_necl = hist_nECLClusters->GetMean();
165 float mean_muid = hist_muid->GetMean();
166 float mean_d0 = hist_D0->GetMean();
167 float mean_z0 = hist_Z0->GetMean();
168 float mean_pval = hist_Pval->GetMean();
169 float mean_ndf = hist_ndf->GetMean();
170 float dif_ndf_ncdc = mean_ndf - mean_ncdc;
171 float mean_dd0 = hist_dD0->GetMean();
172 float mean_dz0 = hist_dZ0->GetMean();
173 float mean_dpt = hist_dPtcms->GetMean();
174 float rms_dd0 = hist_dD0->GetRMS();
175 float rms_dz0 = hist_dZ0->GetRMS();
176 float rms_dpt = hist_dPtcms->GetRMS();
177 int ntot = hist_nsvd->GetEntries();
178 float neve_mumu = ntot;
179 float goodmu_frac = hist_muid->GetBinContent(20) / (float)ntot;
180 float goodmu_o_badmu = hist_muid->GetBinContent(20) / (float)hist_muid->GetBinContent(1);
181 float pval_more95 = 0;
182 float pval_less05 = 0;
183 for (
int i = 95; i < 100; i++) pval_more95 += hist_Pval->GetBinContent(i + 1);
184 for (
int i = 0; i < 5; i++) pval_less05 += hist_Pval->GetBinContent(i + 1);
185 float pval_frac_0 = pval_less05 / (float)ntot;
186 float pval_frac_1 = pval_more95 / (float)ntot;
187 float nocdc_frac = hist_ncdc->GetBinContent(1) / (float)ntot;
188 float notop_frac = hist_topdig->GetBinContent(1) / (float)ntot;
189 float noarich_frac = hist_DetPhotonARICH->GetBinContent(1) / (float)ntot;
227 TH1* hist_D0_InvM =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_InvM");
228 TH1* hist_delta_m =
findHist(
"PhysicsObjectsMiraBelleDst/hist_delta_m");
229 TH1* hist_D0_softpi_PID_ALL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_ALL_pion");
230 TH1* hist_D0_softpi_PID_SVD_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_SVD_pion");
231 TH1* hist_D0_softpi_PID_CDC_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_CDC_pion");
232 TH1* hist_D0_softpi_PID_TOP_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_TOP_pion");
233 TH1* hist_D0_softpi_PID_ARICH_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_ARICH_pion");
234 TH1* hist_D0_softpi_PID_ECL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_ECL_pion");
235 TH1* hist_D0_softpi_PID_KLM_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_KLM_pion");
236 TH1* hist_D0_pi_PID_ALL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_ALL_pion");
237 TH1* hist_D0_pi_PID_SVD_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_SVD_pion");
238 TH1* hist_D0_pi_PID_CDC_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_CDC_pion");
239 TH1* hist_D0_pi_PID_TOP_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_TOP_pion");
240 TH1* hist_D0_pi_PID_ARICH_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_ARICH_pion");
241 TH1* hist_D0_pi_PID_ECL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_ECL_pion");
242 TH1* hist_D0_pi_PID_KLM_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_KLM_pion");
243 TH1* hist_D0_K_PID_ALL_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_ALL_kaon");
244 TH1* hist_D0_K_PID_SVD_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_SVD_kaon");
245 TH1* hist_D0_K_PID_CDC_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_CDC_kaon");
246 TH1* hist_D0_K_PID_TOP_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_TOP_kaon");
247 TH1* hist_D0_K_PID_ARICH_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_ARICH_kaon");
248 TH1* hist_D0_K_PID_ECL_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_ECL_kaon");
249 TH1* hist_D0_K_PID_KLM_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_D0_K_PID_KLM_kaon");
250 TH1* hist_sideband_D0_softpi_PID_ALL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_ALL_pion");
251 TH1* hist_sideband_D0_softpi_PID_SVD_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_SVD_pion");
252 TH1* hist_sideband_D0_softpi_PID_CDC_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_CDC_pion");
253 TH1* hist_sideband_D0_softpi_PID_TOP_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_TOP_pion");
254 TH1* hist_sideband_D0_softpi_PID_ARICH_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_ARICH_pion");
255 TH1* hist_sideband_D0_softpi_PID_ECL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_ECL_pion");
256 TH1* hist_sideband_D0_softpi_PID_KLM_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_KLM_pion");
257 TH1* hist_sideband_D0_pi_PID_ALL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_ALL_pion");
258 TH1* hist_sideband_D0_pi_PID_SVD_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_SVD_pion");
259 TH1* hist_sideband_D0_pi_PID_CDC_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_CDC_pion");
260 TH1* hist_sideband_D0_pi_PID_TOP_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_TOP_pion");
261 TH1* hist_sideband_D0_pi_PID_ARICH_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_ARICH_pion");
262 TH1* hist_sideband_D0_pi_PID_ECL_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_ECL_pion");
263 TH1* hist_sideband_D0_pi_PID_KLM_pion =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_KLM_pion");
264 TH1* hist_sideband_D0_K_PID_ALL_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_ALL_kaon");
265 TH1* hist_sideband_D0_K_PID_SVD_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_SVD_kaon");
266 TH1* hist_sideband_D0_K_PID_CDC_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_CDC_kaon");
267 TH1* hist_sideband_D0_K_PID_TOP_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_TOP_kaon");
268 TH1* hist_sideband_D0_K_PID_ARICH_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_ARICH_kaon");
269 TH1* hist_sideband_D0_K_PID_ECL_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_ECL_kaon");
270 TH1* hist_sideband_D0_K_PID_KLM_kaon =
findHist(
"PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_KLM_kaon");
271 TH1* hist_D0_pi0_InvM =
findHist(
"PhysicsObjectsMiraBelleDst2/hist_D0_pi0_InvM");
274 TF1* f_InvM =
new TF1(
"f_InvM",
"[0]*TMath::Gaus(x,[1],[2])+[3]*pow(x-[4],2)+[5]", 1.81, 1.95);
275 f_InvM->SetParNames(
"Height",
"#mu",
"#sigma",
"a",
"b",
"c");
276 f_InvM->SetParameters(hist_D0_InvM->GetMaximum(), 1.86, 5e-3, 0., 1.86, 0.);
277 hist_D0_InvM->Fit(f_InvM,
"",
"", 1.81, 1.95);
278 f_InvM->SetLineColor(kRed);
280 TF1* f_delta_m =
new TF1(
"f_delta_m",
"[0]*TMath::Gaus(x,[1],[2])+[3]*pow(x-[4],2)+[5]", 0.14, 0.16);
281 f_delta_m->SetParNames(
"Height",
"#mu",
"#sigma",
"a",
"b",
"c");
282 f_delta_m->SetParameters(hist_delta_m->GetMaximum(), 0.145, 5e-4, 0., 0.145, 0.);
283 hist_delta_m->Fit(f_delta_m,
"",
"", 0.14, 0.16);
284 f_delta_m->SetLineColor(kRed);
286 TF1* f_pi0_InvM =
new TF1(
"f_pi0_InvM",
"[0]*TMath::Gaus(x,[1],[2])+[3]*pow(x-[4],2)+[5]", 0.09, 0.17);
287 f_pi0_InvM->SetParNames(
"Height",
"#mu",
"#sigma",
"a",
"b",
"c");
288 f_pi0_InvM->SetParameters(hist_D0_pi0_InvM->GetMaximum(), 0.13, 5e-3, 0., 0.13, 0.);
289 hist_D0_pi0_InvM->Fit(f_pi0_InvM,
"",
"", 0.09, 0.17);
290 f_pi0_InvM->SetLineColor(kRed);
293 TF1* f_gaus_InvM =
new TF1(
"f_gaus_InvM",
"gaus", 1.81, 1.95);
294 f_gaus_InvM->SetParameters(f_InvM->GetParameter(0), f_InvM->GetParameter(1), f_InvM->GetParameter(2));
295 float InvM_bin_width = hist_D0_InvM->GetXaxis()->GetBinWidth(1);
296 float neve_dst = f_gaus_InvM->Integral(1.81, 1.95) / InvM_bin_width;
298 TF1* f_gaus_pi0_InvM =
new TF1(
"f_gaus_pi0_InvM",
"gaus", 0.09, 0.17);
299 f_gaus_pi0_InvM->SetParameters(f_pi0_InvM->GetParameter(0), f_pi0_InvM->GetParameter(1), f_pi0_InvM->GetParameter(2));
300 float pi0_InvM_bin_width = hist_D0_pi0_InvM->GetXaxis()->GetBinWidth(1);
301 float neve_pi0 = f_gaus_pi0_InvM->Integral(0.09, 0.17) / pi0_InvM_bin_width;
304 hist_D0_softpi_PID_ALL_pion->Sumw2();
305 hist_D0_softpi_PID_SVD_pion->Sumw2();
306 hist_D0_softpi_PID_CDC_pion->Sumw2();
307 hist_D0_softpi_PID_TOP_pion->Sumw2();
308 hist_D0_softpi_PID_ARICH_pion->Sumw2();
309 hist_D0_softpi_PID_ECL_pion->Sumw2();
310 hist_D0_softpi_PID_KLM_pion->Sumw2();
311 hist_D0_pi_PID_ALL_pion->Sumw2();
312 hist_D0_pi_PID_SVD_pion->Sumw2();
313 hist_D0_pi_PID_CDC_pion->Sumw2();
314 hist_D0_pi_PID_TOP_pion->Sumw2();
315 hist_D0_pi_PID_ARICH_pion->Sumw2();
316 hist_D0_pi_PID_ECL_pion->Sumw2();
317 hist_D0_pi_PID_KLM_pion->Sumw2();
318 hist_D0_K_PID_ALL_kaon->Sumw2();
319 hist_D0_K_PID_SVD_kaon->Sumw2();
320 hist_D0_K_PID_CDC_kaon->Sumw2();
321 hist_D0_K_PID_TOP_kaon->Sumw2();
322 hist_D0_K_PID_ARICH_kaon->Sumw2();
323 hist_D0_K_PID_ECL_kaon->Sumw2();
324 hist_D0_K_PID_KLM_kaon->Sumw2();
325 hist_sideband_D0_softpi_PID_ALL_pion->Sumw2();
326 hist_sideband_D0_softpi_PID_SVD_pion->Sumw2();
327 hist_sideband_D0_softpi_PID_CDC_pion->Sumw2();
328 hist_sideband_D0_softpi_PID_TOP_pion->Sumw2();
329 hist_sideband_D0_softpi_PID_ARICH_pion->Sumw2();
330 hist_sideband_D0_softpi_PID_ECL_pion->Sumw2();
331 hist_sideband_D0_softpi_PID_KLM_pion->Sumw2();
332 hist_sideband_D0_pi_PID_ALL_pion->Sumw2();
333 hist_sideband_D0_pi_PID_SVD_pion->Sumw2();
334 hist_sideband_D0_pi_PID_CDC_pion->Sumw2();
335 hist_sideband_D0_pi_PID_TOP_pion->Sumw2();
336 hist_sideband_D0_pi_PID_ARICH_pion->Sumw2();
337 hist_sideband_D0_pi_PID_ECL_pion->Sumw2();
338 hist_sideband_D0_pi_PID_KLM_pion->Sumw2();
339 hist_sideband_D0_K_PID_ALL_kaon->Sumw2();
340 hist_sideband_D0_K_PID_SVD_kaon->Sumw2();
341 hist_sideband_D0_K_PID_CDC_kaon->Sumw2();
342 hist_sideband_D0_K_PID_TOP_kaon->Sumw2();
343 hist_sideband_D0_K_PID_ARICH_kaon->Sumw2();
344 hist_sideband_D0_K_PID_ECL_kaon->Sumw2();
345 hist_sideband_D0_K_PID_KLM_kaon->Sumw2();
348 hist_sideband_D0_softpi_PID_ALL_pion->Scale(
m_scale_dst);
349 hist_sideband_D0_softpi_PID_SVD_pion->Scale(
m_scale_dst);
350 hist_sideband_D0_softpi_PID_CDC_pion->Scale(
m_scale_dst);
351 hist_sideband_D0_softpi_PID_TOP_pion->Scale(
m_scale_dst);
352 hist_sideband_D0_softpi_PID_ARICH_pion->Scale(
m_scale_dst);
353 hist_sideband_D0_softpi_PID_ECL_pion->Scale(
m_scale_dst);
354 hist_sideband_D0_softpi_PID_KLM_pion->Scale(
m_scale_dst);
355 hist_sideband_D0_pi_PID_ALL_pion->Scale(
m_scale_dst);
356 hist_sideband_D0_pi_PID_SVD_pion->Scale(
m_scale_dst);
357 hist_sideband_D0_pi_PID_CDC_pion->Scale(
m_scale_dst);
358 hist_sideband_D0_pi_PID_TOP_pion->Scale(
m_scale_dst);
359 hist_sideband_D0_pi_PID_ARICH_pion->Scale(
m_scale_dst);
360 hist_sideband_D0_pi_PID_ECL_pion->Scale(
m_scale_dst);
361 hist_sideband_D0_pi_PID_KLM_pion->Scale(
m_scale_dst);
362 hist_sideband_D0_K_PID_ALL_kaon->Scale(
m_scale_dst);
363 hist_sideband_D0_K_PID_SVD_kaon->Scale(
m_scale_dst);
364 hist_sideband_D0_K_PID_CDC_kaon->Scale(
m_scale_dst);
365 hist_sideband_D0_K_PID_TOP_kaon->Scale(
m_scale_dst);
366 hist_sideband_D0_K_PID_ARICH_kaon->Scale(
m_scale_dst);
367 hist_sideband_D0_K_PID_ECL_kaon->Scale(
m_scale_dst);
368 hist_sideband_D0_K_PID_KLM_kaon->Scale(
m_scale_dst);
371 hist_D0_softpi_PID_ALL_pion->Add(hist_sideband_D0_softpi_PID_ALL_pion, -1);
372 hist_D0_softpi_PID_SVD_pion->Add(hist_sideband_D0_softpi_PID_SVD_pion, -1);
373 hist_D0_softpi_PID_CDC_pion->Add(hist_sideband_D0_softpi_PID_CDC_pion, -1);
374 hist_D0_softpi_PID_TOP_pion->Add(hist_sideband_D0_softpi_PID_TOP_pion, -1);
375 hist_D0_softpi_PID_ARICH_pion->Add(hist_sideband_D0_softpi_PID_ARICH_pion, -1);
376 hist_D0_softpi_PID_ECL_pion->Add(hist_sideband_D0_softpi_PID_ECL_pion, -1);
377 hist_D0_softpi_PID_KLM_pion->Add(hist_sideband_D0_softpi_PID_KLM_pion, -1);
378 hist_D0_pi_PID_ALL_pion->Add(hist_sideband_D0_pi_PID_ALL_pion, -1);
379 hist_D0_pi_PID_SVD_pion->Add(hist_sideband_D0_pi_PID_SVD_pion, -1);
380 hist_D0_pi_PID_CDC_pion->Add(hist_sideband_D0_pi_PID_CDC_pion, -1);
381 hist_D0_pi_PID_TOP_pion->Add(hist_sideband_D0_pi_PID_TOP_pion, -1);
382 hist_D0_pi_PID_ARICH_pion->Add(hist_sideband_D0_pi_PID_ARICH_pion, -1);
383 hist_D0_pi_PID_ECL_pion->Add(hist_sideband_D0_pi_PID_ECL_pion, -1);
384 hist_D0_pi_PID_KLM_pion->Add(hist_sideband_D0_pi_PID_KLM_pion, -1);
385 hist_D0_K_PID_ALL_kaon->Add(hist_sideband_D0_K_PID_ALL_kaon, -1);
386 hist_D0_K_PID_SVD_kaon->Add(hist_sideband_D0_K_PID_SVD_kaon, -1);
387 hist_D0_K_PID_CDC_kaon->Add(hist_sideband_D0_K_PID_CDC_kaon, -1);
388 hist_D0_K_PID_TOP_kaon->Add(hist_sideband_D0_K_PID_TOP_kaon, -1);
389 hist_D0_K_PID_ARICH_kaon->Add(hist_sideband_D0_K_PID_ARICH_kaon, -1);
390 hist_D0_K_PID_ECL_kaon->Add(hist_sideband_D0_K_PID_ECL_kaon, -1);
391 hist_D0_K_PID_KLM_kaon->Add(hist_sideband_D0_K_PID_KLM_kaon, -1);
397 hist_D0_InvM->Draw();
398 f_InvM->Draw(
"SAME");
400 hist_delta_m->Draw();
401 f_delta_m->Draw(
"SAME");
403 hist_D0_pi0_InvM->Draw();
404 f_pi0_InvM->Draw(
"SAME");
407 dst_pi_val->cd(1); hist_D0_pi_PID_ALL_pion->Draw();
408 dst_pi_val->cd(2); hist_D0_pi_PID_SVD_pion->Draw();
409 dst_pi_val->cd(3); hist_D0_pi_PID_CDC_pion->Draw();
410 dst_pi_val->cd(4); hist_D0_pi_PID_TOP_pion->Draw();
411 dst_pi_val->cd(5); hist_D0_pi_PID_ARICH_pion->Draw();
412 dst_pi_val->cd(6); hist_D0_pi_PID_ECL_pion->Draw();
413 dst_pi_val->cd(7); hist_D0_pi_PID_KLM_pion->Draw();
416 dst_k_val->cd(1); hist_D0_K_PID_ALL_kaon->Draw();
417 dst_k_val->cd(2); hist_D0_K_PID_SVD_kaon->Draw();
418 dst_k_val->cd(3); hist_D0_K_PID_CDC_kaon->Draw();
419 dst_k_val->cd(4); hist_D0_K_PID_TOP_kaon->Draw();
420 dst_k_val->cd(5); hist_D0_K_PID_ARICH_kaon->Draw();
421 dst_k_val->cd(6); hist_D0_K_PID_ECL_kaon->Draw();
422 dst_k_val->cd(7); hist_D0_K_PID_KLM_kaon->Draw();
425 float mean_D0_InvM = f_InvM->GetParameter(1);
426 float width_D0_InvM = f_InvM->GetParameter(2);
427 float mean_delta_m = f_delta_m->GetParameter(1);
428 float width_delta_m = f_delta_m->GetParameter(2);
429 float mean_pi0_InvM = f_pi0_InvM->GetParameter(1);
430 float width_pi0_InvM = f_pi0_InvM->GetParameter(2);
431 float mean_D0_softpi_PID_ALL_pion = hist_D0_softpi_PID_ALL_pion->GetMean();
432 float mean_D0_softpi_PID_SVD_pion = hist_D0_softpi_PID_SVD_pion->GetMean();
433 float mean_D0_softpi_PID_CDC_pion = hist_D0_softpi_PID_CDC_pion->GetMean();
434 float mean_D0_softpi_PID_TOP_pion = hist_D0_softpi_PID_TOP_pion->GetMean();
435 float mean_D0_softpi_PID_ARICH_pion = hist_D0_softpi_PID_ARICH_pion->GetMean();
436 float mean_D0_softpi_PID_ECL_pion = hist_D0_softpi_PID_ECL_pion->GetMean();
437 float mean_D0_softpi_PID_KLM_pion = hist_D0_softpi_PID_KLM_pion->GetMean();
438 float mean_D0_pi_PID_ALL_pion = hist_D0_pi_PID_ALL_pion->GetMean();
439 float mean_D0_pi_PID_SVD_pion = hist_D0_pi_PID_SVD_pion->GetMean();
440 float mean_D0_pi_PID_CDC_pion = hist_D0_pi_PID_CDC_pion->GetMean();
441 float mean_D0_pi_PID_TOP_pion = hist_D0_pi_PID_TOP_pion->GetMean();
442 float mean_D0_pi_PID_ARICH_pion = hist_D0_pi_PID_ARICH_pion->GetMean();
443 float mean_D0_pi_PID_ECL_pion = hist_D0_pi_PID_ECL_pion->GetMean();
444 float mean_D0_pi_PID_KLM_pion = hist_D0_pi_PID_KLM_pion->GetMean();
445 float mean_D0_K_PID_ALL_kaon = hist_D0_K_PID_ALL_kaon->GetMean();
446 float mean_D0_K_PID_SVD_kaon = hist_D0_K_PID_SVD_kaon->GetMean();
447 float mean_D0_K_PID_CDC_kaon = hist_D0_K_PID_CDC_kaon->GetMean();
448 float mean_D0_K_PID_TOP_kaon = hist_D0_K_PID_TOP_kaon->GetMean();
449 float mean_D0_K_PID_ARICH_kaon = hist_D0_K_PID_ARICH_kaon->GetMean();
450 float mean_D0_K_PID_ECL_kaon = hist_D0_K_PID_ECL_kaon->GetMean();
451 float mean_D0_K_PID_KLM_kaon = hist_D0_K_PID_KLM_kaon->GetMean();
485 B2DEBUG(20,
"DQMHistAnalysisMiraBelle : endRun called");
491 B2DEBUG(20,
"terminate called");
496 TIter nextckey(gROOT->GetListOfCanvases());
497 TObject* cobj = NULL;
499 while ((cobj = (TObject*)nextckey())) {
500 if (cobj->IsA()->InheritsFrom(
"TCanvas")) {
501 if (cobj->GetName() == canvas_name)
505 return (TCanvas*)cobj;