CbmModelsQA.h

Go to the documentation of this file.

/*
 *====================================================================
 *
 *  CBM Models
 *  
 *  Authors: V.Vovchenko
 *
 *  e-mail : 
 *
 *====================================================================
 *
 *  KF Particles Finder performance 
 *
 *====================================================================
 */
 
#ifndef _CbmModelsQA_h_
#define _CbmModelsQA_h_
 
//#include "KFMCParticle.h"
//#include "KFParticleMatch.h"
#include "CbmModelBase.h"
#include "FairTask.h"
#include "KFQA/CbmKFPartEfficiencies.h"
#include "KFQA/CbmKFTrErrMCPoints.h"
#include <cstring>
#include <map>
#include <vector>
 
class TClonesArray;
class CbmVertex;
class TDirectory;
class TH1F;
class TH2F;
class KFParticleTopoReconstructor;
 
class CbmModelsQA : public FairTask {
public:
  CbmModelsQA(Int_t iVerbose                  = 1,
              int findParticlesMode           = 1,
              int perf                        = 3,
              KFParticleTopoReconstructor* tr = 0,
              const char* name                = "CbmModelsQA",
              const char* title               = "Cbm Models QA",
              float ekin_                     = 24.08);
  ~CbmModelsQA();
 
  void SetOutFile(TString outname) { outfileName = outname; }
 
  virtual InitStatus ReInit();
  virtual InitStatus Init();
  void Exec(Option_t* option);
  void Finish();
 
 
  void AddHRGAnalysis(int TracksType       = -1,
                      double SystError     = 0.,
                      TString name         = "MC",
                      int EventStats       = 1,
                      Bool_t UseWidth      = false,
                      Bool_t UseStatistics = false,
                      double rad           = 0.);
 
  void AddMultiscatteringAnalysis(int TracksType   = -1,
                                  double SystError = 0.,
                                  TString name     = "MC",
                                  int EventStats   = 1);
 
  void AddImpactParameterAnalysis(int TracksType     = -1,
                                  double SystError   = 0.,
                                  TString name       = "MC",
                                  TString InputTable = "");
 
  void AddInverseSlopeAnalysis(int PDGID         = -211,
                               const char* pname = "pi-",
                               int TracksType    = -1,
                               double SystError  = 0.,
                               TString name      = "MC",
                               int EventStats    = 1);
 
  void AddBoltzmannAnalysis(int PDGID         = -211,
                            const char* pname = "pi-",
                            int TracksType    = -1,
                            double SystError  = 0.,
                            TString name      = "MC",
                            int EventStats    = 1);
 
  void AddBlastWaveAnalysis(int PDGID         = -211,
                            const char* pname = "pi-",
                            int TracksType    = -1,
                            double SystError  = 0.,
                            TString name      = "MC",
                            int EventStats    = 1,
                            double Tlong      = 0.132);
 
private:
  void WriteHistos(TObject* obj);
  void WriteHistosCurFile(TObject* obj);
 
  CbmVertex* fPrimVtx;
 
  // Names of files
  TString outfileName;
 
  TDirectory* histodir;
 
  std::vector<int> vStsHitMatch;
  std::vector<int> vStsPointMatch;
  std::vector<int> vMvdPointMatch;
  std::vector<int> vMCTrackMatch;
 
  int fNEvents;
 
  KFParticleTopoReconstructor* fTopoReconstructor;
 
  float ekin;
 
  ClassDef(CbmModelsQA, 1);
 
private:
  CbmModelsQA(const CbmModelsQA&);
  void operator=(const CbmModelsQA&);
  std::vector<CbmModelBase*> Models;
  TDirectory* histodirmod;
};
 
#endif  // _CbmModelsQA_h_