99 const std::vector<CDCRecoHit3D>& foundHits,
100 const typename AInBoxAlgorithm::HoughBox& node)
102 TGraph* allHitsGraph =
new TGraph();
103 allHitsGraph->SetLineWidth(2);
104 allHitsGraph->SetLineColor(9);
107 const Vector3D& recoPos3D = recoHit3D.getRecoPos3D();
108 const double R = std::sqrt(recoPos3D.
x() * recoPos3D.
x() + recoPos3D.
y() * recoPos3D.
y());
109 const double Z = recoPos3D.
z();
110 allHitsGraph->SetPoint(allHitsGraph->GetN(),
R, Z);
113 static int nevent(0);
114 TCanvas canv(
"trackCanvas",
"CDC stereo hits in an event", 0, 0, 1600, 1200);
116 allHitsGraph->Draw(
"APL*");
117 allHitsGraph->GetXaxis()->SetLimits(0, 120);
118 allHitsGraph->GetYaxis()->SetRangeUser(-180, 180);
120 TGraph* foundHitsGraph =
new TGraph();
121 foundHitsGraph->SetMarkerStyle(8);
122 foundHitsGraph->SetMarkerColor(2);
125 const Vector3D& recoPos3D = recoHit3D.getRecoPos3D();
126 const double R = std::sqrt(recoPos3D.
x() * recoPos3D.
x() + recoPos3D.
y() * recoPos3D.
y());
127 const double Z = recoPos3D.
z();
128 foundHitsGraph->SetPoint(foundHitsGraph->GetN(),
R, Z);
130 foundHitsGraph->Draw(
"P");
132 const double centerX = (AInBoxAlgorithm::BoxAlgorithm::centerX(node));
133 const double deltaX = (AInBoxAlgorithm::BoxAlgorithm::deltaX(node));
134 const double centerY = (AInBoxAlgorithm::BoxAlgorithm::centerY(node));
135 const double centerZ = (AInBoxAlgorithm::BoxAlgorithm::centerZ(node));
137 TF1* candidateL =
new TF1(
"candL", AInBoxAlgorithm::BoxAlgorithm::debugLine(), 0, 120);
138 TF1* candidateH =
new TF1(
"candH", AInBoxAlgorithm::BoxAlgorithm::debugLine(), 0, 120);
139 TF1* candidateMean =
new TF1(
"candMean", AInBoxAlgorithm::BoxAlgorithm::debugLine(), 0, 120);
141 candidateL->SetParameters(centerX - deltaX, centerY, centerZ - 100.0 * deltaX);
142 candidateH->SetParameters(centerX + deltaX, centerY, centerZ + 100.0 * deltaX);
143 candidateMean->SetParameters(centerX, centerY, centerZ);
145 candidateL->SetLineColor(9);
146 candidateH->SetLineColor(41);
147 candidateMean->SetLineColor(2);
149 candidateL->Draw(
"same");
150 candidateH->Draw(
"same");
151 candidateMean->Draw(
"same");
152 canv.SaveAs(Form(
"CDCRLHits_%i.png", nevent));