AnalysisDUT: adding normalised associated cluster size plot

src/modules/AnalysisDUT/AnalysisDUT.cpp

@@ -87,6 +87,7 @@ void AnalysisDUT::initialise() {
     clusterTotAssoc = new TH1F("clusterTotAssociated", "clusterTotAssociated", 10000, 0, 10000);
     clusterTotAssocNorm = new TH1F("clusterTotAssociatedNormalized", "clusterTotAssociatedNormalized", 10000, 0, 10000);
     clusterSizeAssoc = new TH1F("clusterSizeAssociated", "clusterSizeAssociated", 30, 0, 30);
+    clusterSizeAssocNorm = new TH1F("clusterSizeAssociatedNormalized", "clusterSizeAssociatedNormalized", 30, 0, 30);
 
     // In-pixel studies:
     auto pitch_x = static_cast<double>(Units::convert(m_detector->pitch().X(), "um"));
@@ -336,6 +337,7 @@ StatusCode AnalysisDUT::run(std::shared_ptr<Clipboard> clipboard) {
                 residualsTimeVsSignal->Fill(tdistance, cluster->tot());
 
                 clusterSizeAssoc->Fill(static_cast<double>(cluster->size()));
+                clusterSizeAssocNorm->Fill(static_cast<double>(cluster->size()));
 
                 // Fill in-pixel plots: (all as function of track position within pixel cell)
                 if(is_within_roi) {
@@ -382,7 +384,10 @@ StatusCode AnalysisDUT::run(std::shared_ptr<Clipboard> clipboard) {
             hPixelEfficiencyMap->Fill(xmod, ymod, has_associated_cluster);
         }
     }
-
     // Return value telling analysis to keep running
     return StatusCode::Success;
 }
+
+void AnalysisDUT::finalise() {
+    clusterSizeAssocNorm->Scale(1 / clusterSizeAssoc->Integral());
+}

src/modules/AnalysisDUT/AnalysisDUT.h

@@ -21,6 +21,7 @@ namespace corryvreckan {
         // Functions
         void initialise();
         StatusCode run(std::shared_ptr<Clipboard> clipboard);
+        void finalise();
 
     private:
         std::shared_ptr<Detector> m_detector;
@@ -40,6 +41,7 @@ namespace corryvreckan {
 
         TH1F *clusterTotAssoc, *clusterTotAssocNorm;
         TH1F* clusterSizeAssoc;
+        TH1F* clusterSizeAssocNorm;
 
         TProfile2D *rmsxvsxmym, *rmsyvsxmym, *rmsxyvsxmym;
         TProfile2D *qvsxmym, *qMoyalvsxmym, *pxqvsxmym;