Skip to content
Snippets Groups Projects
Commit f7f380b7 authored by Jens Kroeger's avatar Jens Kroeger
Browse files

Clustering4D: add histograms for cluster ToT and seed cluster ToT (as no calibration is performed)

parent f6d7ea69
No related branches found
No related tags found
No related merge requests found
......@@ -37,12 +37,16 @@ void Clustering4D::initialise() {
// Cluster plots
std::string title = m_detector->name() + " Cluster size;cluster size;events";
clusterSize = new TH1F("clusterSize", title.c_str(), 100, 0, 100);
title = m_detector->name() + " Cluster seed ToT;cluster seed ToT [lsb];events";
clusterSeedTot = new TH1F("clusterSeedTot", title.c_str(), 256, 0, 256);
title = m_detector->name() + " Cluster seed charge;cluster seed charge [e];events";
clusterSeedCharge = new TH1F("clusterSeedCharge", title.c_str(), 256, 0, 256);
title = m_detector->name() + " Cluster Width - Rows;cluster width [rows];events";
clusterWidthRow = new TH1F("clusterWidthRow", title.c_str(), 25, 0, 25);
title = m_detector->name() + " Cluster Width - Columns;cluster width [columns];events";
clusterWidthColumn = new TH1F("clusterWidthColumn", title.c_str(), 100, 0, 100);
title = m_detector->name() + " Cluster ToT;cluster ToT [lsb];events";
clusterTot = new TH1F("clusterTot", title.c_str(), 1024, 0, 1024);
title = m_detector->name() + " Cluster Charge;cluster charge [e];events";
clusterCharge = new TH1F("clusterCharge", title.c_str(), 5000, 0, 50000);
title = m_detector->name() + " Cluster Position (Global);x [mm];y [mm];events";
......@@ -134,7 +138,9 @@ StatusCode Clustering4D::run(std::shared_ptr<Clipboard> clipboard) {
clusterSize->Fill(static_cast<double>(cluster->size()));
clusterWidthRow->Fill(cluster->rowWidth());
clusterWidthColumn->Fill(cluster->columnWidth());
clusterTot->Fill(cluster->charge()); // no calibration --> same as ToT
clusterCharge->Fill(cluster->charge());
clusterSeedTot->Fill(cluster->getSeedPixel()->charge()); // no calibration --> same as ToT
clusterSeedCharge->Fill(cluster->getSeedPixel()->charge());
clusterPositionGlobal->Fill(cluster->global().x(), cluster->global().y());
clusterTimes->Fill(static_cast<double>(Units::convert(cluster->timestamp(), "ns")));
......
......@@ -42,9 +42,11 @@ namespace corryvreckan {
// Cluster histograms
TH1F* clusterSize;
TH1F* clusterSeedTot;
TH1F* clusterSeedCharge;
TH1F* clusterWidthRow;
TH1F* clusterWidthColumn;
TH1F* clusterTot;
TH1F* clusterCharge;
TH2F* clusterPositionGlobal;
TH1F* clusterTimes;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment