diff --git a/src/modules/Clustering4D/Clustering4D.cpp b/src/modules/Clustering4D/Clustering4D.cpp index 31bd6067dbde0d3a2e7e7a51a3765c4ca3318cf3..cfe5123cfabcbc485fd73cd5b44260e8021e2f61 100644 --- a/src/modules/Clustering4D/Clustering4D.cpp +++ b/src/modules/Clustering4D/Clustering4D.cpp @@ -86,23 +86,29 @@ void Clustering4D::initialize() { title = m_detector->getName() + " Cluster multiplicity;clusters;events"; clusterMultiplicity = new TH1F("clusterMultiplicity", title.c_str(), 50, -0.5, 49.5); title = m_detector->getName() + " Cluster Uncertainty x;cluster uncertainty x [um];events"; - clusterUncertaintyX = new TH1F("clusterUncertaintyX", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().X(), "um"))); + clusterUncertaintyX = new TH1F( + "clusterUncertaintyX", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().X(), "um"))); title = m_detector->getName() + " Cluster Uncertainty y;cluster uncertainty y [um];events"; - clusterUncertaintyY = new TH1F("clusterUncertaintyY", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().Y(), "um"))); + clusterUncertaintyY = new TH1F( + "clusterUncertaintyY", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().Y(), "um"))); title = m_detector->getName() + " Cluster Uncertainty x vs position;column [px];row [px];<cluster uncertainty x> [um]"; - clusterUncertaintyXvsXY = new TProfile2D("clusterUncertaintyXvsXY", title.c_str(),m_detector->nPixels().X(), - -0.5, - m_detector->nPixels().X() - 0.5, - m_detector->nPixels().Y(), - -0.5, - m_detector->nPixels().Y() - 0.5); + clusterUncertaintyXvsXY = new TProfile2D("clusterUncertaintyXvsXY", + title.c_str(), + m_detector->nPixels().X(), + -0.5, + m_detector->nPixels().X() - 0.5, + m_detector->nPixels().Y(), + -0.5, + m_detector->nPixels().Y() - 0.5); title = m_detector->getName() + " Cluster Uncertainty Y vs position;column [px];row [px];<cluster uncertainty y> [um]"; - clusterUncertaintyYvsXY = new TProfile2D("clusterUncertaintyYvsXY", title.c_str(),m_detector->nPixels().X(), - -0.5, - m_detector->nPixels().X() - 0.5, - m_detector->nPixels().Y(), - -0.5, - m_detector->nPixels().Y() - 0.5); + clusterUncertaintyYvsXY = new TProfile2D("clusterUncertaintyYvsXY", + title.c_str(), + m_detector->nPixels().X(), + -0.5, + m_detector->nPixels().X() - 0.5, + m_detector->nPixels().Y(), + -0.5, + m_detector->nPixels().Y() - 0.5); >>>>>>> cluster_error_plots title = @@ -229,8 +235,10 @@ StatusCode Clustering4D::run(const std::shared_ptr<Clipboard>& clipboard) { clusterTimes->Fill(static_cast<double>(Units::convert(cluster->timestamp(), "ns"))); clusterUncertaintyX->Fill(static_cast<double>(Units::convert(cluster->errorX(), "um"))); clusterUncertaintyY->Fill(static_cast<double>(Units::convert(cluster->errorY(), "um"))); - clusterUncertaintyXvsXY->Fill(cluster->column(), cluster->row(), static_cast<double>(Units::convert(cluster->errorX(), "um"))); - clusterUncertaintyYvsXY->Fill(cluster->column(), cluster->row(), static_cast<double>(Units::convert(cluster->errorY(), "um"))); + clusterUncertaintyXvsXY->Fill( + cluster->column(), cluster->row(), static_cast<double>(Units::convert(cluster->errorX(), "um"))); + clusterUncertaintyYvsXY->Fill( + cluster->column(), cluster->row(), static_cast<double>(Units::convert(cluster->errorY(), "um"))); // to check that cluster timestamp = earliest pixel timestamp if(cluster->size() > 1) { diff --git a/src/modules/ClusteringSpatial/ClusteringSpatial.cpp b/src/modules/ClusteringSpatial/ClusteringSpatial.cpp index ed8e5211fa028aece0c2263ebedfc91a62c6c332..69a32e3daeb4765b62bdd8fcedf824b9b8eea73d 100644 --- a/src/modules/ClusteringSpatial/ClusteringSpatial.cpp +++ b/src/modules/ClusteringSpatial/ClusteringSpatial.cpp @@ -64,9 +64,11 @@ void ClusteringSpatial::initialize() { title = m_detector->getName() + " Cluster multiplicity;clusters;events"; clusterMultiplicity = new TH1F("clusterMultiplicity", title.c_str(), 50, -0.5, 49.5); title = m_detector->getName() + " Cluster Uncertainty x;cluster uncertainty x [um];events"; - clusterUncertaintyX = new TH1F("clusterUncertaintyX", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().X(), "um"))); + clusterUncertaintyX = new TH1F( + "clusterUncertaintyX", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().X(), "um"))); title = m_detector->getName() + " Cluster Uncertainty y;cluster uncertainty y [um];events"; - clusterUncertaintyY = new TH1F("clusterUncertaintyY", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().Y(), "um"))); + clusterUncertaintyY = new TH1F( + "clusterUncertaintyY", title.c_str(), 100, 0, static_cast<double>(Units::convert(m_detector->getPitch().Y(), "um"))); } StatusCode ClusteringSpatial::run(const std::shared_ptr<Clipboard>& clipboard) {