PndFairHitQATask.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

#ifndef PndFairHitQATask_H
#define PndFairHitQATask_H

// Root includes
#include "TVector3.h"
#include "TString.h"
#include "TGeoManager.h"
#include "TH1.h"

// framework includes
#include "PndPersistencyTask.h"
#include "PndMeasurementBoundaries.h"

#include "FairLogger.h"
#include "FairHit.h"

#include <map>

class TClonesArray;

class PndFairHitQATask : public PndPersistencyTask {
 public:
  PndFairHitQATask(std::string hitBranchName, std::string pointBranchName);

  virtual ~PndFairHitQATask();

  virtual void SetParContainers();
  virtual InitStatus Init();

  virtual void Exec(Option_t *opt);
  virtual void FinishTask();
  void SetWorkdir(TString path) { fWorkdir = path; } 

  void SetHitBranch(std::string hitBranch) { fHitBranchName = hitBranch; }
  void SetPointBranch(std::string pointBranch) { fPointBranchName = pointBranch; }
  void SetOutputPlotName(std::string val = "Output.gif") { fOutputPlotName = val; }

 protected:
  // These methods should be overwritten for a realistic test

  virtual void InitHistos();

  virtual void SetAllBoundaries();

  virtual void CalculateResolutions(FairHit *hit);

  virtual void FillHistos();

  virtual void AssignMeasurements();

  // No need to touch these classes
  bool Test();
  void DrawHistos();
  void SaveHistos();
  void PrintDartMeasurement();

  std::string fHitBranchName;
  std::string fPointBranchName;
  std::string fOutputPlotName{"ResolutionHistos.png"};

  TClonesArray *fMcPoints = nullptr;
  TClonesArray *fHits = nullptr;
  Bool_t fProperDifference = false; 
  TVector3 fSpatialDifference;
  double fTimeDifference = -1.0;
  TString fWorkdir;

  PndMeasurementBoundaries fBoundaries;

  std::map<std::string, TH1 *> fHistos;

  ClassDef(PndFairHitQATask, 2);
};

#endif