//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Implementation of resCalc
//      Uses existing GFTracks as reference and calculates
//      residua in XY for TpcClusters.
//      Developped for Test Bench tracking systems using the
//      RK track representation
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Felix Boehmer                     (original author)
//      Maxence Vandenbroucke             (Application of TpcCdcFit2DResCalc to GFTracks)
//      Physik Department E18, TUM
//
//-----------------------------------------------------------

#ifndef TPCREFGFTRKRES_HH
#define TPCREFGFTRKRES_HH

#include <TH1.h>
#include <TH2.h>

#include "TString.h"
#include "AbsRefTrackResCalc.h"
#include "TpcAlignmentManager.h"
class GFRecoHitFactory;

class TpcRefGFTrkResCalc : public AbsRefTrackResCalc {

 public:
  TpcRefGFTrkResCalc();
  virtual ~TpcRefGFTrkResCalc();

  virtual int calc();
  virtual bool init();
 
  void setTpcSPHitBranchName(const TString& name) {fTpcSPHitBranchName=name;}
  void setGFTrackBranchName(const TString& name) {fGFTrackBranchName=name;}
  void setRefGFTrackBranchName(const TString& name) {fRefGFTrackBranchName=name;}
  void setTpcClusterBranchName(const TString& name) {fTpcClusterBranchName=name;}

  void setTrackRepId(const int id){fNRep=id;}
  void setRequireGoodFit(bool opt=true) {fRequireGoodFit = opt;}
  void setUseClusters(bool opt=true){fUseClusters=opt;}
  void writeHists();
  
 private:
  TClonesArray* fRefGFTrackArray;
  TClonesArray* fTpcSPHitArray;
  TClonesArray* fGFTrackArray;
  TClonesArray* fTpcClusterArray;
  TString fRefGFTrackBranchName;
  TString fGFTrackBranchName;
  TString fTpcSPHitBranchName;
  TString fTpcClusterBranchName;
  
  TpcAlignmentManager* fAlMan;

  int fNRep;                       //which rep to use	
  int fDetID;

  bool fRequireGoodFit;
  bool fUseClusters;
  bool fAlignedMatching;
  
  double fTrackdPhiCut;
  double fTrackMomentumCut;
  double fTrackPosCut;

  TH1D* dPhiHist;
  TH2D* dPhiVsPhiHist;
  TH1D* dPhiPosHist;
  TH2D* dPhiVsPhiPosHist;

  TH1D* dPhiTpcTrkPosHist;
  TH1D* dPhiRefTrkPosHist;

  TH1D* rRefHist;
  TH2D* xyRefHist;
  TH1D* momRefHist;
  TH1D* phiRefHist;

  TH1D* rRefCutHist;
  TH2D* xyRefCutHist;
  TH1D* momRefCutHist;
  TH1D* phiRefCutHist;
  TH1D* dPhiRefCutTrkPosHist;
  TH1D* dPhiRefCutTrkPosHist2;


  TH1D* rHist;
  TH1D* momHist;
  TH1D* phiHist;


  TH1D* rHistCut;
  TH1D* momHistCut;
  std::vector<TH1*> hists;
  TString fQAFileName;
  

 public:
  ClassDef(TpcRefGFTrkResCalc,3)
  
};

#endif