CbmStsFindTracksQa.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                  CbmStsFindTracksQa header file               -----
// -----                  Created 11/01/06  by V. Friese               -----
// -------------------------------------------------------------------------
 
 
#ifndef CBMSTSTRACKFINDERQA_H
#define CBMSTSTRACKFINDERQA_H 1
 
#include "FairTask.h"
#include "TStopwatch.h"
#include "TVector3.h"
 
class TH1F;
class CbmEvent;
class CbmMCDataArray;
class CbmStsSetup;
 
 
class CbmStsFindTracksQa : public FairTask {
 
public:
  CbmStsFindTracksQa(Int_t iVerbose = 1);
 
 
  CbmStsFindTracksQa(Int_t minHits, Double_t quota, Int_t iVerbose);
 
 
  virtual ~CbmStsFindTracksQa();
 
 
  virtual void SetParContainers();
 
 
  virtual InitStatus Init();
 
 
  virtual InitStatus ReInit();
 
 
  virtual void Exec(Option_t* opt);
 
 
private:
  virtual void Finish();
 
  InitStatus GetGeometry();
 
 
  void GetTargetPosition();
 
 
  void CreateHistos();
 
 
  void Reset();
 
 
  void FillHitMap(CbmEvent* event);
 
 
  void
  FillMatchMap(CbmEvent* event, Int_t& nRec, Int_t& nGhosts, Int_t& nClones);
 
 
  void DivideHistos(TH1* histo1, TH1* histo2, TH1* histo3);
 
 
  void ProcessEvent(CbmEvent* event = NULL);
 
 
  std::map<Int_t, std::map<Int_t, Int_t>> fHitMap;
 
 
  std::map<Int_t, Int_t> fMatchMap;
 
 
  std::map<Int_t, Double_t> fQualiMap;
 
 
  TClonesArray* fEvents;       
  CbmMCDataArray* fMCTracks;   
  CbmMCDataArray* fStsPoints;  
  TClonesArray* fStsHits;      
  TClonesArray* fStsHitMatch;  
  TClonesArray* fStsTracks;    
  TClonesArray* fMatches;      
 
  Bool_t fLegacy;
 
 
  TVector3 fTargetPos;  // Target centre position
  CbmStsSetup* fSetup;  // STS setup interface
  Int_t fNStations;     // Number of STS stations
 
 
  Int_t
    fMinStations;  // Minimal number of stations with hits for considered MCTrack
  Double_t fQuota;  // True/all hits for track to be considered reconstructed
 
 
  TH1F *fhMomAccAll, *fhMomRecAll, *fhMomEffAll;     // eff. vs. p, all
  TH1F *fhMomAccPrim, *fhMomRecPrim, *fhMomEffPrim;  // eff. vs. p, vertex
  TH1F *fhMomAccSec, *fhMomRecSec, *fhMomEffSec;     // eff. vs. p, non-vertex
  TH1F *fhNpAccAll, *fhNpRecAll, *fhNpEffAll;        // eff. vs. np, all
  TH1F *fhNpAccPrim, *fhNpRecPrim, *fhNpEffPrim;     // eff. vs. np, vertex
  TH1F *fhNpAccSec, *fhNpRecSec, *fhNpEffSec;        // eff. vs. np, non-vertex
  TH1F *fhZAccSec, *fhZRecSec, *fhZEffSec;           // eff. vs. z, non-vertex
  TH1F *fhNhClones, *fhNhGhosts;  // # hits of clones and ghosts
 
 
  TList* fHistoList;
 
 
  Int_t fNAccAll, fNAccPrim, fNAccRef, fNAccSec;
  Int_t fNRecAll, fNRecPrim, fNRecRef, fNRecSec;
  Int_t fNGhosts, fNClones;
  Int_t fNEvents;       
  Int_t fNEventsFailed; 
  Double_t fTime;       
  TStopwatch fTimer;
 
  CbmStsFindTracksQa(const CbmStsFindTracksQa&);
  CbmStsFindTracksQa operator=(const CbmStsFindTracksQa&);
 
  ClassDef(CbmStsFindTracksQa, 2);
};
 
 
#endif