PndGemFindHitsQA.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                   PndGemFindHitsQA header file          -----
// -----                  Created 02.06.2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------

#ifndef PNDGEMFINDHITSQA_H
#define PNDGEMFINDHITSQA_H

#include "TH1F.h"
#include "TH2F.h"
#include "TH3F.h"

#include "FairMCPoint.h"
#include "FairTask.h"
#include "FairRootFileSink.h"

#include "PndGemHit.h"
#include "PndGemDigiPar.h"

#include <vector>

class PndGemFindHitsQA : public FairTask {

 public:
  PndGemFindHitsQA();

  PndGemFindHitsQA(Int_t iVerbose);

  virtual ~PndGemFindHitsQA();

  virtual void Exec(Option_t *opt);

  void SetVerbose(const Int_t &verbose) { fVerbose = verbose; };

  void SetPointEffDist(Double_t dn) { fPointEffDist = dn; };

 private:
  PndGemDigiPar *fDigiPar;

  TClonesArray *fMCPointArray;

  TClonesArray *fGemHitArray;

  TClonesArray *fGemData[10];
  Int_t fGemDataPointer[1000];
  Int_t fGemPointNumber;

  Int_t fNofEvents; 

  TList *fHistoList;

  Int_t fHistPlaneDivs;
  Double_t fHistWidth[4][2];

  Double_t fPointEffDist;

  // XY distance from point to hit
  // vs station(up to 4),sensor(2 per station),position(divide into 4x4 regions)
  TH2F *fhPointToHit[4][2][4][4];
  // distance of point to the closest hit
  //  TH3F* fhPointClosest     [4][2];
  TH2F *fhPointNof[4][2];
  TH2F *fhPointReco[4][2];
  TH2F *fhPointRecoEff[4][2];
  TH1F *fhPointRadNof[4][2];
  TH1F *fhPointRadReco[4][2];
  TH1F *fhPointRadRecoEff[4][2];

  TH2F *fhPointMatch[4][2];
  TH2F *fhPointMatchEff[4][2];
  TH1F *fhPointRadMatch[4][2];
  TH1F *fhPointRadMatchEff[4][2];

  TH2F *fhHitNof[4][2];
  TH2F *fhHitFake[4][2];
  TH2F *fhHitFakeProb[4][2];
  TH1F *fhHitRadNof[4][2];
  TH1F *fhHitRadFake[4][2];
  TH1F *fhHitRadFakeProb[4][2];

  TH2F *fhHitMultipleRate[4][2];

  TH1F *fhCloseHits[4][2];

  TH2F *fhTrueMatchDiXYPerSt[4][2];
  TH1F *fhTrueMatchDistPerSt[4][2];

  TH2F *fhTrueMatchDiXY;
  TH1F *fhTrueMatchDist;
  TH1F *fhTrueMatchValue;
  TH2F *fhTrueMatchDistValue;

  TH1F *fhTrueMatchNofPerHit;
  TH1F *fhTrueMatchNofPerPoint;

  void CreateHistos();
  void DivideHistos(TH1 *hist1, TH1 *hist2, TH1 *hist3);

  Int_t GetPointVector(Int_t arrayId, Int_t entryId, std::vector<Int_t> &pointVector, Bool_t printInfo = kFALSE);

  virtual void SetParContainers();

  virtual void Finish();

  virtual InitStatus Init();

  virtual InitStatus ReInit();

  ClassDef(PndGemFindHitsQA, 1);
};

#endif