AnalysisATLASpix: add histograms to investigate left/right/main peak of time...

AnalysisATLASpix: add histograms to investigate left/right/main peak of time residual (also for 1px clusters)

src/modules/AnalysisTimingATLASpix/AnalysisTimingATLASpix.cpp

@@ -44,6 +44,7 @@ AnalysisTimingATLASpix::AnalysisTimingATLASpix(Configuration config, std::shared
         m_clusterSizeCut = m_config.get<size_t>("cluster_size_cut");
     }
     m_highTotCut = m_config.get<int>("high_tot_cut", 40);
+    m_lowTotCut = m_config.get<int>("low_tot_cut", 40);
     m_highChargeCut = m_config.get<double>("high_charge_cut", static_cast<double>(m_highTotCut));
     m_leftTailCut = m_config.get<double>("left_tail_cut", static_cast<double>(Units::convert(-10, "ns")));
     m_rightTailCut = m_config.get<double>("right_tail_cut", static_cast<double>(Units::convert(10, "ns")));
@@ -361,9 +362,24 @@ void AnalysisTimingATLASpix::initialise() {
     hTot_leftTail =
         new TH1F("hTot_leftTail", "ToT (left tail of time residual);seed pixel ToT [lsb]; # events", 2 * 64, -64, 64);
     hTot_rightTail =
-        new TH1F("hTot_rightTail", "ToT (left tail of time residual);seed pixel ToT [lsb]; # events", 2 * 64, -64, 64);
-    hTot_mainPeak =
-        new TH1F("hTot_mainPeak", "ToT (main peak of time residual);seed pixel ToT [lsb]; # events", 2 * 64, -64, 64);
+        new TH1F("hTot_rightTail", "ToT (right tail of time residual);seed pixel ToT [lsb]; # events", 2 * 64, -64, 64);
+    hTot_mainPeak = new TH1F(
+        "hTot_mainPeak", "ToT (main peak of time residual, 1px clusters);seed pixel ToT [lsb]; # events", 2 * 64, -64, 64);
+    hTot_leftTail_1px = new TH1F("hTot_leftTail_1px",
+                                 "ToT (left tail of time residual, 1px clusters);seed pixel ToT [lsb]; # events",
+                                 2 * 64,
+                                 -64,
+                                 64);
+    hTot_rightTail_1px = new TH1F("hTot_rightTail_1px",
+                                  "ToT (right tail of time residual, 1px clusters);seed pixel ToT [lsb]; # events",
+                                  2 * 64,
+                                  -64,
+                                  64);
+    hTot_mainPeak_1px = new TH1F("hTot_mainPeak_1px",
+                                 "ToT (main peak of time residual, 1px clusters);seed pixel ToT [lsb]; # events",
+                                 2 * 64,
+                                 -64,
+                                 64);
     hPixelTimestamp_leftTail = new TH1F("pixelTimestamp_leftTail",
                                         "pixelTimestamp (left tail of time residual); pixel timestamp [ms]; # events",
                                         3e6,
@@ -652,16 +668,25 @@ StatusCode AnalysisTimingATLASpix::run(std::shared_ptr<Clipboard> clipboard) {
                             hTot_leftTail->Fill(cluster->getSeedPixel()->raw());
                             hPixelTimestamp_leftTail->Fill(cluster->getSeedPixel()->timestamp());
                             hClusterSize_leftTail->Fill(static_cast<double>(cluster->size()));
+                            if(cluster->size() == 1) {
+                                hTot_leftTail_1px->Fill(cluster->getSeedPixel()->raw());
+                            }
                         } else if(track->timestamp() - cluster->timestamp() > m_rightTailCut) {
                             hClusterMap_rightTail->Fill(cluster->column(), cluster->row());
                             hTot_rightTail->Fill(cluster->getSeedPixel()->raw());
                             hPixelTimestamp_rightTail->Fill(cluster->getSeedPixel()->timestamp());
                             hClusterSize_rightTail->Fill(static_cast<double>(cluster->size()));
+                            if(cluster->size() == 1) {
+                                hTot_rightTail_1px->Fill(cluster->getSeedPixel()->raw());
+                            }
                         } else {
                             hClusterMap_mainPeak->Fill(cluster->column(), cluster->row());
                             hTot_mainPeak->Fill(cluster->getSeedPixel()->raw());
                             hPixelTimestamp_mainPeak->Fill(cluster->getSeedPixel()->timestamp());
                             hClusterSize_mainPeak->Fill(static_cast<double>(cluster->size()));
+                            if(cluster->size() == 1) {
+                                hTot_mainPeak_1px->Fill(cluster->getSeedPixel()->raw());
+                            }
                         }
                     }
                 }

src/modules/AnalysisTimingATLASpix/AnalysisTimingATLASpix.h

@@ -80,6 +80,7 @@ namespace corryvreckan {
 
         TH2F* hTotVsTime;
         TH2F* hTotVsTime_high;
+        TH2F* hTotVsTime_low;
 
         // Control Plots for "left/right tail" and "main peak" of time correlation
         TH2F* hInPixelMap_leftTail;
@@ -91,6 +92,9 @@ namespace corryvreckan {
         TH1F* hTot_leftTail;
         TH1F* hTot_rightTail;
         TH1F* hTot_mainPeak;
+        TH1F* hTot_leftTail_1px;
+        TH1F* hTot_rightTail_1px;
+        TH1F* hTot_mainPeak_1px;
         TH1F* hPixelTimestamp_leftTail;
         TH1F* hPixelTimestamp_rightTail;
         TH1F* hPixelTimestamp_mainPeak;
@@ -115,6 +119,7 @@ namespace corryvreckan {
         double m_clusterChargeCut;
         size_t m_clusterSizeCut;
         int m_highTotCut;       // for pixel->tot()
+        int m_lowTotCut;        // for pixel->tot()
         double m_highChargeCut; // for cluster->charge()
         double m_leftTailCut;
         double m_rightTailCut;