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

Clustering4D: rename histograms and clarify axis titles: px - seed

parent d6eace47
......@@ -84,24 +84,25 @@ void Clustering4D::initialize() {
clusterTimes = new TH1F("clusterTimes", title.c_str(), 3e6, 0, 3e9);
title = m_detector->getName() + " Cluster multiplicity;clusters;events";
clusterMultiplicity = new TH1F("clusterMultiplicity", title.c_str(), 50, -0.5, 49.5);
title = m_detector->getName() + " pixel - cluster timestamp;ts_{pixel} - ts_{cluster} [ns] (all pixels w/o seed);events";
pxTimeMinusClsTime = new TH1F("pxTimeMinusClsTime", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625);
title =
m_detector->getName() + " pixel - seed pixel timestamp (all pixels w/o seed);ts_{pixel} - ts_ {seed} [ns];events";
pxTimeMinusSeedTime = new TH1F("pxTimeMinusSeedTime", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625);
title = m_detector->getName() +
" pixel - cluster timestamp;ts_{pixel} - ts_{cluster} [ns] (all pixels w/o seed); pixel charge [e];events";
pxTimeMinusClsTime_vs_pxCharge =
new TH2F("pxTimeMinusClsTime_vs_pxCharge", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
" pixel - seed pixel timestamp (all pixels w/o seed);ts_{pixel} - ts_ {seed} [ns]; pixel charge [e];events";
pxTimeMinusSeedTime_vs_pxCharge =
new TH2F("pxTimeMinusSeedTime_vs_pxCharge", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
title = m_detector->getName() +
" pixel - cluster timestamp;ts_{pixel} - ts_{cluster} [ns] (all pixels w/o seed); pixel charge [e];events";
pxTimeMinusClsTime_vs_pxCharge_2px = new TH2F(
"pxTimeMinusClsTime_vs_pxCharge_2px", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
" pixel - seed pixel timestamp (all pixels w/o seed);ts_{pixel} - ts_ {seed} [ns]; pixel charge [e];events";
pxTimeMinusSeedTime_vs_pxCharge_2px = new TH2F(
"pxTimeMinusSeedTime_vs_pxCharge_2px", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
title = m_detector->getName() +
" pixel - cluster timestamp;ts_{pixel} - ts_{cluster} [ns] (all pixels w/o seed); pixel charge [e];events";
pxTimeMinusClsTime_vs_pxCharge_3px = new TH2F(
"pxTimeMinusClsTime_vs_pxCharge_3px", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
" pixel - seed pixel timestamp (all pixels w/o seed);ts_{pixel} - ts_ {seed} [ns]; pixel charge [e];events";
pxTimeMinusSeedTime_vs_pxCharge_3px = new TH2F(
"pxTimeMinusSeedTime_vs_pxCharge_3px", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
title = m_detector->getName() +
" pixel - cluster timestamp;ts_{pixel} - ts_{cluster} [ns] (all pixels w/o seed); pixel charge [e];events";
pxTimeMinusClsTime_vs_pxCharge_4px = new TH2F(
"pxTimeMinusClsTime_vs_pxCharge_4px", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
" pixel - seed pixel timestamp (all pixels w/o seed);ts_{pixel} - ts_ {seed} [ns]; pixel charge [e];events";
pxTimeMinusSeedTime_vs_pxCharge_4px = new TH2F(
"pxTimeMinusSeedTime_vs_pxCharge_4px", title.c_str(), 1000, -99.5 * 1.5625, 900.5 * 1.5625, 256, -0.5, 255.5);
// Get resolution in time of detector and calculate time cut to be applied
LOG(DEBUG) << "Time cut to be applied for " << m_detector->getName() << " is "
......@@ -211,17 +212,17 @@ StatusCode Clustering4D::run(const std::shared_ptr<Clipboard>& clipboard) {
if(px == cluster->getSeedPixel()) {
continue; // don't fill this histogram for seed pixel!
}
pxTimeMinusClsTime->Fill(static_cast<double>(Units::convert(px->timestamp() - cluster->timestamp(), "ns")));
pxTimeMinusClsTime_vs_pxCharge->Fill(
pxTimeMinusSeedTime->Fill(static_cast<double>(Units::convert(px->timestamp() - cluster->timestamp(), "ns")));
pxTimeMinusSeedTime_vs_pxCharge->Fill(
static_cast<double>(Units::convert(px->timestamp() - cluster->timestamp(), "ns")), px->charge());
if(cluster->size() == 2) {
pxTimeMinusClsTime_vs_pxCharge_2px->Fill(
pxTimeMinusSeedTime_vs_pxCharge_2px->Fill(
static_cast<double>(Units::convert(px->timestamp() - cluster->timestamp(), "ns")), px->charge());
} else if(cluster->size() == 3) {
pxTimeMinusClsTime_vs_pxCharge_3px->Fill(
pxTimeMinusSeedTime_vs_pxCharge_3px->Fill(
static_cast<double>(Units::convert(px->timestamp() - cluster->timestamp(), "ns")), px->charge());
} else if(cluster->size() == 4) {
pxTimeMinusClsTime_vs_pxCharge_4px->Fill(
pxTimeMinusSeedTime_vs_pxCharge_4px->Fill(
static_cast<double>(Units::convert(px->timestamp() - cluster->timestamp(), "ns")), px->charge());
}
}
......
......@@ -52,11 +52,11 @@ namespace corryvreckan {
TH2F* clusterPositionLocal;
TH1F* clusterTimes;
TH1F* clusterMultiplicity;
TH1F* pxTimeMinusClsTime;
TH2F* pxTimeMinusClsTime_vs_pxCharge;
TH2F* pxTimeMinusClsTime_vs_pxCharge_2px;
TH2F* pxTimeMinusClsTime_vs_pxCharge_3px;
TH2F* pxTimeMinusClsTime_vs_pxCharge_4px;
TH1F* pxTimeMinusSeedTime;
TH2F* pxTimeMinusSeedTime_vs_pxCharge;
TH2F* pxTimeMinusSeedTime_vs_pxCharge_2px;
TH2F* pxTimeMinusSeedTime_vs_pxCharge_3px;
TH2F* pxTimeMinusSeedTime_vs_pxCharge_4px;
double time_cut_;
int neighbor_radius_row_;
......
Supports Markdown
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