RawTowerCalibration.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RawTowerCalibration.h

#ifndef CALORECO_RAWTOWERCALIBRATION_H
#define CALORECO_RAWTOWERCALIBRATION_H

#include <fun4all/SubsysReco.h>

#include <phparameter/PHParameters.h>

#include <cmath>
#include <iostream>
#include <string>

class CDBTTree;
class PHCompositeNode;
class RawTowerContainer;
class TowerInfoContainerv1;
class RawTowerGeomContainer;

class RawTowerCalibration : public SubsysReco
{
 public:
  RawTowerCalibration(const std::string &name = "RawTowerCalibration");
  ~RawTowerCalibration() override;

  int InitRun(PHCompositeNode *topNode) override;
  int process_event(PHCompositeNode *topNode) override;
  int End(PHCompositeNode *topNode) override;
  void Detector(const std::string &d)
  {
    m_Detector = d;
    _tower_calib_params.set_name(d);
  }
  void TowerType(const int type)
  {
    _tower_type = type;
  }

  enum enu_calib_algorithm
  {
    kNo_calibration = 0,

    kSimple_linear_calibration = 1,

    kTower_by_tower_calibration = 2,

    // use conditions DB file/wrapper (non-xml) file for most gain tracing correction factors
    kDbfile_tbt_gain_corr = 3
  };
  enum ProcessTowerType
  {
    kRawTowerOnly = 0,
    kTowerInfoOnly = 1,
    kBothTowers = 2
  };

  enu_calib_algorithm get_calib_algorithm() const
  {
    return _calib_algorithm;
  }

  void set_calib_algorithm(enu_calib_algorithm calibAlgorithm)
  {
    _calib_algorithm = calibAlgorithm;
  }

  double get_calib_const_GeV_ADC() const
  {
    return _calib_const_GeV_ADC;
  }

  void set_calib_const_GeV_ADC(double calibConstGeVAdc)
  {
    _calib_const_GeV_ADC = calibConstGeVAdc;
  }

  void set_variable_GeV_ADC(const bool value)
  {
    _GeV_ADC_file = value;
  }

  std::string get_calib_tower_node_prefix() const
  {
    return _calib_tower_node_prefix;
  }

  void set_calib_tower_node_prefix(const std::string &calibTowerNodePrefix)
  {
    _calib_tower_node_prefix = calibTowerNodePrefix;
  }

  double get_pedstal_ADC() const
  {
    return _pedstal_ADC;
  }

  void set_pedstal_ADC(double pedstalAdc)
  {
    _pedstal_ADC = pedstalAdc;
  }

  void set_variable_pedestal(const bool value)
  {
    _pedestal_file = value;
  }

  std::string get_raw_tower_node_prefix() const
  {
    return _raw_tower_node_prefix;
  }

  void set_raw_tower_node_prefix(const std::string &rawTowerNodePrefix)
  {
    _raw_tower_node_prefix = rawTowerNodePrefix;
  }

  void set_zero_suppression_GeV(double)
  {
    std::cout << "RawTowerCalibration::set_zero_suppression_GeV is deprecated!" << std::endl
              << "  See discussion at https://github.com/sPHENIX-Collaboration/coresoftware/pull/867" << std::endl
              << std::endl;
  }

  PHParameters &GetCalibrationParameters()
  {
    return _tower_calib_params;
  }

  void set_CalibrationFileName(const char *inCalFname)
  {
    m_CalibrationFileName = inCalFname;
  }
  void set_UseConditionsDB(const bool setUseCondDB)
  {
    m_UseConditionsDB = setUseCondDB;
  }

  void set_towerinfo(RawTowerCalibration::ProcessTowerType UseTowerInfo)
  {
    m_UseTowerInfo = UseTowerInfo;
  }

 protected:
  void CreateNodes(PHCompositeNode *topNode);


  RawTowerContainer *_calib_towers = nullptr;
  RawTowerContainer *_raw_towers = nullptr;

  TowerInfoContainerv1 *_calib_towerinfos = nullptr;
  TowerInfoContainerv1 *_raw_towerinfos = nullptr;

  RawTowerGeomContainer *rawtowergeom = nullptr;

  enu_calib_algorithm _calib_algorithm;

  std::string m_Detector;
  std::string RawTowerNodeName;
  std::string RawTowerInfoNodeName;
  std::string CaliTowerNodeName;
  std::string CaliTowerInfoNodeName;
  std::string TowerGeomNodeName;

  std::string _calib_tower_node_prefix;
  std::string _raw_tower_node_prefix;

  double _pedstal_ADC = NAN;

  bool _pedestal_file = false;

  double _calib_const_GeV_ADC = NAN;

  bool _GeV_ADC_file = false;

  int _tower_type = -1;

  PHParameters _tower_calib_params;

  std::string m_CalibrationFileName;
  bool m_UseConditionsDB = false;

  CDBTTree *m_CDBTTree = nullptr;
  RawTowerCalibration::ProcessTowerType m_UseTowerInfo = RawTowerCalibration::ProcessTowerType::kBothTowers;  // 0 just produce RawTowers, 1 just produce TowerInfo objects, and 2 produce both
};

#endif