19 #include "TIterator.h"
37 std::cout <<
"Opening file: " << inFile << std::endl;
38 TFile inputFile(inFile.c_str());
39 TList* layers = inputFile.GetListOfKeys();
40 std::cout <<
"Layers to print: " << std::endl;
46 TCanvas* canvas1 =
new TCanvas();
47 TCanvas* canvas2 =
new TCanvas();
49 while ((obj =
next())) {
51 TDirectory* dir = inputFile.GetDirectory(obj->GetName());
52 TH2F* r_z = (TH2F*)dir->Get(
"r_z");
53 TH2F* r_z_assigned = (TH2F*)dir->Get(
"r_z_assigned");
54 TH2F* x_y = (TH2F*)dir->Get(
"x_y");
55 TH2F* x_y_assigned = (TH2F*)dir->Get(
"x_y_assigned");
56 if (entry == 17) entry = 20;
57 if (entry == 10 || entry == 50) entry++;
58 if (r_z && r_z_assigned) {
61 r_z->SetMarkerColor(TColor::GetColorDark(entry));
62 r_z->SetMarkerStyle(6);
63 r_z->GetXaxis()->SetTitle(
"z");
64 r_z->GetYaxis()->SetTitle(
"r");
67 r_z_assigned->SetStats(0);
68 r_z_assigned->SetMarkerColor(TColor::GetColorBright(entry));
69 r_z_assigned->SetMarkerStyle(6);
70 r_z_assigned->GetXaxis()->SetTitle(
"z");
71 r_z_assigned->GetYaxis()->SetTitle(
"r");
72 r_z_assigned->Draw(
"same");
75 r_z_assigned->SetDirectory(0);
77 if (x_y && x_y_assigned) {
80 x_y->SetMarkerColor(TColor::GetColorDark(entry));
81 x_y->SetMarkerStyle(6);
82 x_y->GetXaxis()->SetTitle(
"x");
83 x_y->GetYaxis()->SetTitle(
"y");
86 x_y_assigned->SetStats(0);
87 x_y_assigned->SetMarkerColor(TColor::GetColorBright(entry));
88 x_y_assigned->SetMarkerStyle(6);
89 x_y_assigned->GetXaxis()->SetTitle(
"x");
90 x_y_assigned->GetYaxis()->SetTitle(
"y");
91 x_y_assigned->Draw(
"same");
94 x_y_assigned->SetDirectory(0);
101 std::cout <<
"Opening file: " << infile_geoRZ << std::endl;
102 TFile inputFileRZ(infile_geoRZ.c_str());
104 TH2F* geo_rz = (TH2F*)inputFileRZ.Get(histName_geoRZ.c_str());
107 geo_rz->Draw(
"same");
109 std::cout <<
"Can not access histogram with name: " << histName_geoRZ
111 geo_rz->SetDirectory(0);
115 std::cout <<
"Opening file: " << infile_geoXY << std::endl;
116 TFile inputFileXY(infile_geoXY.c_str());
118 TH2F* geo_xy = (TH2F*)inputFileXY.Get(histName_geoXY.c_str());
121 geo_xy->Draw(
"same");
123 std::cout <<
"Can not access histogram with name: " << histName_geoXY
125 geo_xy->SetDirectory(0);