HGCalRadiationMap.h

Go to the documentation of this file.

#ifndef simcalorimetry_hgcalsimalgos_hgcalradiationmap
#define simcalorimetry_hgcalsimalgos_hgcalradiationmap

#include "DataFormats/ForwardDetId/interface/HGCScintillatorDetId.h"
#include "Geometry/HGCalGeometry/interface/HGCalGeometry.h"
#include "Geometry/HcalTowerAlgo/interface/HcalGeometry.h"
#include "vdt/vdtMath.h"
#include <string>

typedef std::array<double, 8> radiiVec;

class HGCalRadiationMap {
public:
  struct DoseParameters {
    DoseParameters() : a_(0.), b_(0.), c_(0.), d_(0.), e_(0.), f_(0.), g_(0.), h_(0.), i_(0.), j_(0.) {}
    double a_, b_, c_, d_, e_, f_, g_, h_, i_, j_;
  };

  HGCalRadiationMap(){};
  ~HGCalRadiationMap(){};

  typedef std::map<std::pair<int, int>, DoseParameters> doseParametersMap;

  void setGeometry(const CaloSubdetectorGeometry *);
  void setDoseMap(const std::string &, const unsigned int);

  double getDoseValue(const int, const int, const radiiVec &, bool logVal = false);
  double getFluenceValue(const int, const int, const radiiVec &, bool logVal = false);

  const unsigned int &algo() { return algo_; }
  const HGCalGeometry *geom() { return hgcalGeom_; }
  const HGCalTopology *topo() { return hgcalTopology_; }
  const HGCalDDDConstants *ddd() { return hgcalDDD_; }

  inline const doseParametersMap &getDoseMap() { return doseMap_; }

private:
  doseParametersMap readDosePars(const std::string &);

  unsigned int algo_;
  const HGCalGeometry *hgcalGeom_;
  const HGCalTopology *hgcalTopology_;
  const HGCalDDDConstants *hgcalDDD_;
  doseParametersMap doseMap_;
  //conversion from grey to krad
  const double grayToKrad_ = 0.1;
};

#endif