CbmTrdHitProducerQa.h

Go to the documentation of this file.

/* Copyright (C) 2006-2014 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Matus Kalisky, Denis Bertini [committer], Florian Uhlig, Dominik Smith */
 
// -----------------------------------------------------------------------
// -----                     CbmTrdHitProducerQa                     -----
// -----               Created 13/12/05  by M. Kalisky               -----
// -----------------------------------------------------------------------
 
 
#ifndef CBMTRDHITPRODUCERQA
#define CBMTRDHITPRODUCERQA 1
 
#include "FairTask.h"
 
#include <TFolder.h>
 
class CbmDigiManager;
class CbmMCDataArray;
class CbmQaCanvas;
class TClonesArray;
class TH1F;
 
class CbmTrdHitProducerQa : public FairTask {
 public:
  /* Defaul constructor */
  CbmTrdHitProducerQa();
 
  /* Standard constructor */
  CbmTrdHitProducerQa(const char* name, const char* title = "FairTask");
 
  /* Destructor */
  virtual ~CbmTrdHitProducerQa();
 
  /* Initialisation */
  InitStatus Init();
 
  /* Execution */
  virtual void Exec(Option_t* option);
 
  /* Finish at the end of each event */
  virtual void Finish();
 
  /* Set the momentum cuts */
  void SetMomentumCuts(float CutLower, float CutHigher)
  {
    fMomCutLower = CutLower;
    fMomCutUpper = CutHigher;
  }
 
  /* Set number of TRD stations */
  void SetNumberStations(int nStations) { fNoTrdStations = nStations; }
 
  void SetLayersPerStations(int nLayers) { fNoTrdPerStation = nLayers; }
 
 private:
  TFolder fOutFolder;             
  TFolder* histFolder = nullptr;  
 
  CbmDigiManager* fDigiMan = nullptr;
 
  /* Data branches*/
  TClonesArray* fTrdHitCollection = nullptr;
  CbmMCDataArray* fMCTrackArray   = nullptr;
  CbmMCDataArray* fTrdPoints      = nullptr;
 
  int fNoTrdStations = 4;
 
  int fNoTrdPerStation = 1;
 
  /* Write test histograms */
  void WriteHistograms();
 
  /* Test histograms*/
  std::vector<TH1F*> fvhHitPullX;
  std::vector<TH1F*> fvhHitPullY;
  std::vector<TH1F*> fvhHitPullT;
  std::vector<TH1F*> fvhHitResX;
  std::vector<TH1F*> fvhHitResY;
  std::vector<TH1F*> fvhHitResT;
 
  std::vector<TH1F*> fvhedEcut;
  std::vector<TH1F*> fvhedEall;
  std::vector<TH1F*> fvhpidEcut;
  std::vector<TH1F*> fvhpidEall;
 
  /* Canvases */
  std::vector<CbmQaCanvas*> fvdECanvas;
  std::vector<CbmQaCanvas*> fvPullCanvas;
 
  /* Momentum cuts for energy distributions */
  float fMomCutLower = 1.;
  float fMomCutUpper = 7.;
 
  ClassDef(CbmTrdHitProducerQa, 3)
};
 
#endif