Commit 57f5e28c authored by Jens Kroeger's avatar Jens Kroeger
Browse files

Clustering4D: removed outdated debugging histograms

parent a3cf5339
......@@ -29,28 +29,6 @@ void Clustering4D::initialise() {
clusterPositionGlobal = new TH2F("clusterPositionGlobal", title.c_str(), 400, -10., 10., 400, -10., 10.);
title = ";cluster timestamp [ns]; # events";
clusterTimes = new TH1F("clusterTimes", title.c_str(), 3e6, 0, 3e9);
// only temporary histograms for debugging
hDistXClusterPixel =
new TH1D("hDistXClusterPixel", "hDistXClusterPixel; cluster centre x - pixel pos x [um]; # events", 550, -50, 500);
hDistYClusterPixel =
new TH1D("hDistYClusterPixel", "hDistYClusterPixel; cluster centre y - pixel pos y [um]; # events", 550, -50, 500);
hDistXClusterPixel_1px = new TH1D(
"hDistXClusterPixel_1px", "hDistXClusterPixel_1px; cluster centre x - pixel pos x [um]; # events", 550, -50, 500);
hDistYClusterPixel_1px = new TH1D(
"hDistYClusterPixel_1px", "hDistYClusterPixel_1px; cluster centre y - pixel pos y [um]; # events", 550, -50, 500);
hDistXClusterPixel_2px = new TH1D(
"hDistXClusterPixel_2px", "hDistXClusterPixel_2px; cluster centre x - pixel pos x [um]; # events", 550, -50, 500);
hDistYClusterPixel_2px = new TH1D(
"hDistYClusterPixel_2px", "hDistYClusterPixel_2px; cluster centre y - pixel pos y [um]; # events", 550, -50, 500);
hDistXClusterPixel_3px = new TH1D(
"hDistXClusterPixel_3px", "hDistXClusterPixel_3px; cluster centre x - pixel pos x [um]; # events", 550, -50, 500);
hDistYClusterPixel_3px = new TH1D(
"hDistYClusterPixel_3px", "hDistYClusterPixel_3px; cluster centre y - pixel pos y [um]; # events", 550, -50, 500);
hDistXClusterPixel_npx = new TH1D(
"hDistXClusterPixel_npx", "hDistXClusterPixel_npx; cluster centre x - pixel pos x [um]; # events", 550, -50, 500);
hDistYClusterPixel_npx = new TH1D(
"hDistYClusterPixel_npx", "hDistYClusterPixel_npx; cluster centre y - pixel pos y [um]; # events", 550, -50, 500);
}
// Sort function for pixels from low to high times
......@@ -262,37 +240,4 @@ void Clustering4D::calculateClusterCentre(Cluster* cluster) {
cluster->setDetectorID(detectorID);
cluster->setClusterCentre(positionGlobal);
cluster->setClusterCentreLocal(positionLocal);
// for debugging: loop over pixels of cluster again and histogram them:
for(auto& pixel : (*cluster->pixels())) {
auto pixelPosLocal = m_detector->getLocalPosition(pixel->column(), pixel->row());
hDistXClusterPixel->Fill(static_cast<double>(Units::convert(abs(cluster->local().x() - pixelPosLocal.x()), "um")));
hDistYClusterPixel->Fill(static_cast<double>(Units::convert(abs(cluster->local().y() - pixelPosLocal.y()), "um")));
if(cluster->size() == 1) {
hDistXClusterPixel_1px->Fill(
static_cast<double>(Units::convert(abs(cluster->local().x() - pixelPosLocal.x()), "um")));
hDistYClusterPixel_1px->Fill(
static_cast<double>(Units::convert(abs(cluster->local().y() - pixelPosLocal.y()), "um")));
}
if(cluster->size() == 2) {
hDistXClusterPixel_2px->Fill(
static_cast<double>(Units::convert(abs(cluster->local().x() - pixelPosLocal.x()), "um")));
hDistYClusterPixel_2px->Fill(
static_cast<double>(Units::convert(abs(cluster->local().y() - pixelPosLocal.y()), "um")));
}
if(cluster->size() == 3) {
hDistXClusterPixel_3px->Fill(
static_cast<double>(Units::convert(abs(cluster->local().x() - pixelPosLocal.x()), "um")));
hDistYClusterPixel_3px->Fill(
static_cast<double>(Units::convert(abs(cluster->local().y() - pixelPosLocal.y()), "um")));
}
if(cluster->size() > 1) {
hDistXClusterPixel_npx->Fill(
static_cast<double>(Units::convert(abs(cluster->local().x() - pixelPosLocal.x()), "um")));
hDistYClusterPixel_npx->Fill(
static_cast<double>(Units::convert(abs(cluster->local().y() - pixelPosLocal.y()), "um")));
}
} // end for
}
......@@ -39,18 +39,6 @@ namespace corryvreckan {
TH2F* clusterPositionGlobal;
TH1F* clusterTimes;
// temp histogram for debuggin:
TH1D* hDistXClusterPixel;
TH1D* hDistYClusterPixel;
TH1D* hDistXClusterPixel_1px;
TH1D* hDistYClusterPixel_1px;
TH1D* hDistXClusterPixel_2px;
TH1D* hDistYClusterPixel_2px;
TH1D* hDistXClusterPixel_3px;
TH1D* hDistYClusterPixel_3px;
TH1D* hDistXClusterPixel_npx;
TH1D* hDistYClusterPixel_npx;
double timingCut;
int neighbourRadiusRow;
int neighbourRadiusCol;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment