Commit e3f89e20 authored by Jens Kroeger's avatar Jens Kroeger
Browse files

Merge branch 'clusterMultiplicity' into 'master'

Adding plot for cluster multiplicity (clusters per plane per event)

See merge request !158
parents 1cc74b56 20ff4964
Pipeline #1074007 passed with stages
in 17 minutes and 20 seconds
......@@ -28,6 +28,8 @@ 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);
title = m_detector->name() + " Cluster multiplicity;clusters;events";
clusterMultiplicity = new TH1F("clusterMultiplicity", title.c_str(), 50, 0, 50);
}
// Sort function for pixels from low to high times
......@@ -116,6 +118,8 @@ StatusCode Clustering4D::run(std::shared_ptr<Clipboard> clipboard) {
deviceClusters->push_back(cluster);
}
clusterMultiplicity->Fill(static_cast<double>(deviceClusters->size()));
// Put the clusters on the clipboard
if(deviceClusters->size() > 0) {
clipboard->put(m_detector->name(), "clusters", reinterpret_cast<Objects*>(deviceClusters));
......
......@@ -38,6 +38,7 @@ namespace corryvreckan {
TH1F* clusterCharge;
TH2F* clusterPositionGlobal;
TH1F* clusterTimes;
TH1F* clusterMultiplicity;
double timingCut;
int neighbour_radius_row;
......
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