PndMvdQATask.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                 PndMvdQATask header file             -----
// -----                  Created 20/3/2019  by T.Stockmanns               -----
// -------------------------------------------------------------------------

#ifndef PndMvdQATask_H
#define PndMvdQATask_H

//#pragma once

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

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

#include <array>

class PndMvdQATask : public PndPersistencyTask {
 public:
  PndMvdQATask();

  virtual ~PndMvdQATask();

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

  virtual void Exec(Option_t *opt);
  virtual void FinishTask();
  void SetWorkdir(TString path) { fWorkdir = path; } // only needed to give cdash the path to the stored pictures

 protected:
  void CalculateDifference(PndSdsHit *hit);
  std::array<Double_t, 6> DoubleGaussFit(TH1 *);
  bool TestPixel();
  bool TestStrip();
  void DrawPixelHistos();
  void DrawStripHistos();
  void SetPixelBoundaries();
  void SetStripBoundaries();

  void PrintDartMeasurement();

 private:
  PndGeoHandling *fGeoH;
  TClonesArray *fMcPoints;
  TClonesArray *fStripHits;
  TClonesArray *fPixelHits;
  Bool_t fProperDifference; 
  TVector3 fSpatialDifference;
  TString fWorkdir;
  double fTimeDifference;
  double fEnergyDifference;

  PndMeasurementBoundaries fBoundaries;

  std::array<double, 6> fPixelSpatialRes;
  std::vector<std::pair<double, double>> fPixelBoundaries;
  std::array<double, 6> fStripSpatialRes;
  std::vector<std::pair<double, double>> fStripBoundaries;

  TH1 *fHStripSingle;
  TH1 *fHStripDouble;
  TH1 *fHStripAll;
  TH1 *fHStripClusterSize;
  TH1 *fHStripTime;
  TH1 *fHStripEnergy;

  TH1 *fHPixelSingle;
  TH1 *fHPixelDouble;
  TH1 *fHPixelAll;
  TH1 *fHPixelClusterSize;
  TH1 *fHPixelTime;
  TH1 *fHPixelEnergy;

  ClassDef(PndMvdQATask, 1);
};

#endif