Add new tracker calibration class for TRD system to be used to adjust

error parametrization for reconstructed hits. To be checked by class
CbmTrackerInputQaTrd

macro/run/run_qa.C

@@ -177,6 +177,7 @@ void run_qa(TString dataTra = "data/sis100_muon_jpsi_test", TString dataRaw = "d
     //run->AddTask(new CbmTrdDigitizerPRFQa()); //works put currently doesn't do anything
     //run->AddTask(new CbmTrdHitRateFastQa());  //opens lots of windows
     run->AddTask(new CbmTrdHitProducerQa());
+    run->AddTask(new CbmTrdCalibTracker());
     run->AddTask(new CbmTrackerInputQaTrd());  // Tracker requirements to TRD
   }
   // ------------------------------------------------------------------------

reco/detectors/trd/CMakeLists.txt

@@ -75,6 +75,7 @@ qa/CbmTrdQa.cxx
 qa/CbmTrdRecoQa.cxx
 qa/CbmTrdTracksPidQa.cxx
 
+qa/CbmTrdCalibTracker.cxx
 qa/data/CbmTrdHitMC.cxx
 
 pid/CbmTrdElectronsTrainAnn.cxx

reco/detectors/trd/CbmTrdRecoLinkDef.h

@@ -27,6 +27,7 @@
 #pragma link C++ class CbmTrdQa + ;
 #pragma link C++ class CbmTrdRecoQa + ;
 #pragma link C++ class CbmTrdTracksPidQa + ;
+#pragma link C++ class CbmTrdCalibTracker + ;
 #pragma link C++ class CbmTrdHitMC + ;
 
 #pragma link C++ class std::vector<std::pair<unsigned long, unsigned long>> + ;

reco/detectors/trd/qa/CbmTrdCalibTracker.h

+/* Copyright (C) 2022 Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest
+   SPDX-License-Identifier: GPL-3.0-only
+   Authors: Alexandru Bercuci[committer]*/
+
+/// @file CbmTrdCalibTracker.cxx
+/// @author Alexandru Bercuci
+/// @date 02.05.2022
+
+#ifndef CBMTRDCALIBTRACKER_H
+#define CBMTRDCALIBTRACKER_H
+
+#include "CbmQaCanvas.h"
+#include "CbmQaHist.h"
+
+#include <FairTask.h>
+
+#include <TFolder.h>
+#include <TH1D.h>
+#include <TH1F.h>
+#include <TH1I.h>
+#include <TH2F.h>
+#include <TParameter.h>
+#include <TProfile.h>
+#include <TProfile2D.h>
+
+class CbmDigiManager;
+class CbmMCDataManager;
+class CbmMCEventList;
+class CbmMCDataArray;
+class CbmTimeSlice;
+class CbmTrdParSetGeo;
+class CbmTrdParSetDigi;
+
+class TClonesArray;
+
+///
+/// CbmTrdCalibTracker class represents all the CA tracker requirements for the TRD detector.
+/// When this QA test is passed, the tracker must work (at least that is the idea).
+///
+/// The class ensures that the tracker has the correct understanding of the TRD geometry and interfaces
+/// and checks the quality of all tracking-related aspects of the TRD data.
+///
+class CbmTrdCalibTracker : public FairTask {
+
+public:
+  /// Constructor
+  CbmTrdCalibTracker(const char* name = "TrdQaTracker", Int_t verbose = 1);
+
+  /// Destructor
+  ~CbmTrdCalibTracker();
+
+  /// FairTask: Intialisation at begin of run.
+  InitStatus Init() { return ReInit(); }
+
+  /// FairTask: Reinitialisation.
+  InitStatus ReInit();
+
+  /// FairTask:  Intialise parameter containers.
+  void SetParContainers();
+
+  /// FairTask: Action at end of run. For this task and all of the subtasks.
+  void Finish();
+
+  /// TTask: Clear all data structures created by a previous execution of a task.
+  void Clear(Option_t* /*option*/ = "") {}
+
+  /// TTask: Process a timeslice
+  void Exec(Option_t*);
+
+  /// Prepare the Qa output and return it as a reference to an internal TFolder.
+  /// The reference is non-const, because FairSink can not write const objects
+  TFolder& GetQa();
+
+private:
+  /// Check the geometry
+  InitStatus GeometryQa();
+
+  /// Analysis of hit uncertainty (pull) distributions
+  void ResolutionQa();
+
+  /// Free the memory etc.
+  void DeInit();
+
+
+  // Setup
+
+  Bool_t fIsTrdInSetup {false};
+  Bool_t fIsMcPresent {false};
+
+  Int_t fNtrackingStations {0};
+
+  CbmTimeSlice* fTimeSlice {nullptr};
+  CbmTrdParSetGeo* fTrdGeoPar {nullptr};
+  CbmTrdParSetDigi* fTrdDigiPar {nullptr};
+  CbmDigiManager* fDigiManager {nullptr};
+
+  /// MC data
+  CbmMCEventList* fMcEventList {nullptr};  // list of MC events connected to the current time slice
+  CbmMCDataManager* fMcManager {nullptr};
+  CbmMCDataArray* fMcTracks {nullptr};
+  CbmMCDataArray* fMcPoints {nullptr};
+
+  /// Data
+  TClonesArray* fClusters {nullptr};
+  TClonesArray* fHits {nullptr};
+  TClonesArray* fHitMatches {nullptr};
+
+  /// Output
+  TClonesArray* fHitsMc {nullptr};
+
+
+  TFolder fOutFolder {"TrdQaTracker", "TrdQaTracker"};  /// output folder with histos and canvases
+  TFolder* fHistFolder {nullptr};                       /// subfolder for histograms
+
+  TParameter<int> fNevents {"nEvents", 0};  /// number of processed events
+
+  /// Histogram for Residual Distribution
+  CbmQaHist<TH1D> fh1DresidualU {"h1DresidualU", "Trd1D: Residual U;(U_{reco} - U_{MC})(cm)", 100, -.5, .5};
+  CbmQaHist<TH1D> fh1DresidualV {"h1DresidualV", "Trd1D: Residual V;(V_{reco} - V_{MC})(cm)", 100, -10, 10};
+  CbmQaHist<TH1D> fh1DresidualT {"h1DresidualT", "Trd1D: Residual T;(T_{reco} - T_{MC})(ns)", 100, -50, 50};
+
+  CbmQaHist<TH1D> fh2DresidualX {"h2DresidualX", "Trd2D: Residual X;(X_{reco} - X_{MC})(cm)", 100, -0.05, 0.05};
+  CbmQaHist<TH1D> fh2DresidualY {"h2DresidualY", "Trd2D: Residual Y;(Y_{reco} - Y_{MC})(cm)", 1000, -2, 2};
+  CbmQaHist<TH1D> fh2DresidualT {"h2DresidualT", "Trd2D: Residual T;(T_{reco} - T_{MC})(ns)", 100, -100, 100};
+
+  /// Histogram for PULL Distribution
+  CbmQaHist<TH1D> fh1DpullU {"h1DpullU", "Trd1D: Pull U;(U_{reco} - U_{MC}) / #sigmaU_{reco}", 100, -5, 5};
+  CbmQaHist<TH1D> fh1DpullV {"h1DpullV", "Trd1D: Pull V;(V_{reco} - V_{MC}) / #sigmaV_{reco}", 100, -5, 5};
+  CbmQaHist<TH1D> fh1DpullT {"h1DpullT", "Trd1D: Pull T;(T_{reco} - T_{MC}) / #sigmaT_{reco}", 100, -5, 5};
+
+  CbmQaHist<TH1D> fh2DpullX {"h2DpullX", "Trd2D: Pull X;(X_{reco} - X_{MC}) / #sigmaX_{reco}", 100, -5, 5};
+  CbmQaHist<TH1D> fh2DpullY {"h2DpullY", "Trd2D: Pull Y;(Y_{reco} - Y_{MC}) / #sigmaY_{reco}", 100, -5, 5};
+  CbmQaHist<TH1D> fh2DpullT {"h2DpullT", "Trd2D: Pull T;(T_{reco} - T_{MC}) / #sigmaT_{reco}", 100, -5, 5};
+
+  /// List of the above histogramms
+  std::vector<CbmQaHist<TH1D>*> fHistList;
+
+  /// hits purity
+  std::vector<CbmQaHist<TH1I>> fhPointsPerHit;
+
+  /// hits efficiency
+  std::vector<CbmQaHist<TH1I>> fhHitsPerPoint;
+
+  /// hits efficiency
+  std::vector<CbmQaHist<TProfile2D>> fhEfficiencyXY;
+  std::vector<CbmQaHist<TProfile>> fhEfficiencyR;
+
+  /// Canvaces: collection of histogramms
+
+  CbmQaCanvas fCanvResidual {"cResidual", "Residual Distribution", 3 * 500, 2 * 500};
+  CbmQaCanvas fCanvPull {"cPull", "Pull Distribution", 3 * 500, 2 * 500};
+  CbmQaCanvas fCanvEfficiencyXY {"cEfficiencyXY", "Efficiency XY: % reconstructed McPoint", 2 * 500, 2 * 500};
+  CbmQaCanvas fCanvEfficiencyR {"cEfficiencyR", "Efficiency R: % reconstructed McPoint", 2 * 500, 2 * 500};
+  CbmQaCanvas fCanvHitsPerPoint {"cHitsPerMcPoint", "Efficiency: Hits Per McPoint", 2 * 500, 2 * 500};
+  CbmQaCanvas fCanvPointsPerHit {"cMcPointsPerHit", "Purity: McPoints per Hit", 2 * 500, 2 * 500};
+
+private:
+  /// Suppressed copy constructor
+  CbmTrdCalibTracker(const CbmTrdCalibTracker&) = delete;
+
+  /// Suppressed assignment operator
+  CbmTrdCalibTracker& operator=(const CbmTrdCalibTracker&) = delete;
+
+  ClassDef(CbmTrdCalibTracker, 0);
+};
+
+#endif