CbmLitRadLengthQa.h

Go to the documentation of this file.

/* Copyright (C) 2013-2014 GSI/JINR-LIT, Darmstadt/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Andrey Lebedev [committer] */
 
#ifndef CBMLITRADLENGTHQA_H_
#define CBMLITRADLENGTHQA_H_
 
#include "FairTask.h"
#include "base/CbmLitDetectorSetup.h"
 
#include <map>
#include <string>
using std::map;
using std::pair;
using std::string;
class CbmHistManager;
class TClonesArray;
 
class CbmLitRadLengthQa : public FairTask {
 public:
  CbmLitRadLengthQa();
 
  virtual ~CbmLitRadLengthQa();
 
  virtual InitStatus Init();
 
  virtual void Exec(Option_t* opt);
 
  virtual void Finish();
 
  /* Setters */
  void SetOutputDir(const string& outputDir) { fOutputDir = outputDir; }
 
 private:
  void ReadDataBranches();
 
  void CreateHistograms();
 
  void ExecDetector(const string& pathPattern, const string& detName);
 
  void ExecDetector(const string& detName, Int_t (*getStationId)(const TString&));
 
  void FillHistosDetector(const map<Int_t, map<Int_t, Double_t>>& parMap, const string& histName,
                          map<Int_t, map<Int_t, pair<Double_t, Double_t>>>& xyOnTrack);
 
  static Int_t GetMvdStationId(const TString& nodePath);
 
  static Int_t GetStsStationId(const TString& nodePath);
 
  static Int_t GetTrdStationId(const TString& nodePath);
 
  static Int_t GetMuchStationId(const TString& nodePath);
 
  static Int_t GetMuchAbsorberId(const TString& nodePath);
 
  void SaveMaterialBudgetToFile();
 
  void SaveDetectorMaterialBudgetToFile(const string& detName);
 
  // Material properties of silicon
  //static const Double_t SILICON_DENSITY = 2.33; // g*cm^-3
  //static const Double_t SILICON_A = 28.08855; // silicon atomic weight
  //static const Double_t SILICON_Z = 14.0; // silicon atomic number
  static const Double_t SILICON_RAD_LENGTH;  // cm
 
  CbmHistManager* fHM;  // Histogram manager
  string fOutputDir;    // Output directory for results
 
  // Pointers to data arrays
  TClonesArray* fRadLen;  // RadLen array
 
  CbmLitDetectorSetup fDet;  // Detector setup
 
  CbmLitRadLengthQa(const CbmLitRadLengthQa&);
  CbmLitRadLengthQa& operator=(const CbmLitRadLengthQa&);
 
  ClassDef(CbmLitRadLengthQa, 1);
};
 
#endif /* CBMLITRADLENGTHQA_H_ */