10 #include <dqm/analysis/modules/DQMHistAnalysisV0.h>
35 addParam(
"OverlayPath", m_OverlayPath,
"Path to CAD drawings", std::string(
""));
40 void DQMHistAnalysisV0Module::initialize()
42 B2INFO(
"DQMHistAnalysisV0: initialized.");
45 for (
int i = 0; i < 32; i++) {
46 m_c_xvsy[i] =
new TCanvas(Form(
"V0Object/c_xvsy[%i]", i), Form(
"c_xvsy[%i]", i), 800, 800);
49 m_c_xvsz =
new TCanvas(
"V0Object/c_xvsz",
"c_xvsz", 1500, 400);
51 auto* m_fh =
new TFile(Form(
"%s/v0cad.root", m_OverlayPath.c_str()));
53 contLevelXY.resize(32);
54 for (
int i = 0; i < 32; i++) {
55 contLevelXY[i] =
new TList();
57 if (m_fh->cd(Form(
"h_%dc", i))) {
58 for (
int j = 0; j < 500; j++) {
59 auto curv = gDirectory->Get(Form(
"Graph_%d", j));
61 contLevelXY[i]->AddLast(curv);
66 contLevelXZ =
new TList();
68 if (m_fh->cd(
"h_xzc")) {
69 for (
int j = 0;
true; j++) {
70 auto curv = gDirectory->Get(Form(
"Graph_%d", j));
72 contLevelXZ->AddLast(curv);
81 void DQMHistAnalysisV0Module::event()
85 gStyle->SetOptStat(0);
86 gStyle->SetPalette(kViridis, 0, 0.7);
87 for (
int i = 0; i < 32; i++) {
88 TH2* h = (TH2*) findHist(Form(
"V0Objects/xvsy[%i]", i));
90 if (h) h->Draw(
"COLZ");
92 TList* c = contLevelXY[i];
93 if (c && c->GetSize() > 0) {
94 auto* curv = (TGraph*)c->First();
95 for (
int j = 0; j < c->GetSize(); j++) {
99 curv = (TGraph*)c->After(curv);
103 m_c_xvsy[i]->Modified();
104 m_c_xvsy[i]->Update();
107 TH2* hxz = (TH2*) findHist(
"V0Objects/xvsz");
113 if (hxz) hxz->Draw(
"COLZ");
115 TList* c = contLevelXZ;
116 if (c && c->GetSize() > 0) {
117 auto* curv = (TGraph*)c->First();
118 for (
int j = 0; j < c->GetSize(); j++) {
122 curv = (TGraph*)c->After(curv);
127 m_c_xvsz->Modified();