HGCalRadiationMap Class Reference

parses a txt file with dose/fluence parameters and provides functions for noise, etc. More...

#include <HGCalRadiationMap.h>

Inheritance diagram for HGCalRadiationMap:

Classes
struct	DoseParameters

Public Types
typedef std::map< std::pair < int, int >, DoseParameters >	doseParametersMap

Public Member Functions
const unsigned int &	algo ()
double	computeRadius (const HGCScintillatorDetId &)
const HGCalDDDConstants *	ddd ()
const HGCalGeometry *	geom ()
const doseParametersMap &	getDoseMap ()
double	getDoseValue (const int, const int, const double, bool logVal=false)
double	getFluenceValue (const int, const int, const double, bool logVal=false)
	HGCalRadiationMap ()
void	setDoseMap (const std::string &, const unsigned int)
void	setFluenceScaleFactor (double val)
void	setGeometry (const CaloSubdetectorGeometry *)
const HGCalTopology *	topo ()
	~HGCalRadiationMap ()

Private Member Functions
doseParametersMap	readDosePars (const std::string &)

Private Attributes
unsigned int	algo_
doseParametersMap	doseMap_
double	fluenceSFlog10_
const double	grayToKrad_ = 0.1
const HGCalDDDConstants *	hgcalDDD_
const HGCalGeometry *	hgcalGeom_
const HGCalTopology *	hgcalTopology_

Detailed Description

parses a txt file with dose/fluence parameters and provides functions for noise, etc.

Definition at line 14 of file HGCalRadiationMap.h.

Member Typedef Documentation

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

Definition at line 23 of file HGCalRadiationMap.h.

Constructor & Destructor Documentation

HGCalRadiationMap::HGCalRadiationMap

(

)

Definition at line 6 of file HGCalRadiationMap.cc.

: fluenceSFlog10_(0.) {}

HGCalRadiationMap::~HGCalRadiationMap ( )

inline

Definition at line 23 of file HGCalRadiationMap.h.

{};

Member Function Documentation

const unsigned int& HGCalRadiationMap::algo ( )

inline

Definition at line 35 of file HGCalRadiationMap.h.

References algo_.

Referenced by setDoseMap().

{ return algo_; }

double HGCalRadiationMap::computeRadius

(

const HGCScintillatorDetId &

cellId

)

Definition at line 22 of file HGCalRadiationMap.cc.

References geom(), HGCalGeometry::getPosition(), funct::pow(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by HGCHEbackDigitizer::runRealisticDigitizer().

                                                                          {
  GlobalPoint global = geom()->getPosition(cellId);
  return std::sqrt(std::pow(global.x(), 2) + std::pow(global.y(), 2));
}

const HGCalDDDConstants* HGCalRadiationMap::ddd ( )

inline

Definition at line 38 of file HGCalRadiationMap.h.

References hgcalDDD_.

{ return hgcalDDD_; }

const HGCalGeometry* HGCalRadiationMap::geom ( )

inline

Definition at line 36 of file HGCalRadiationMap.h.

References hgcalGeom_.

Referenced by computeRadius(), and setGeometry().

{ return hgcalGeom_; }

const doseParametersMap& HGCalRadiationMap::getDoseMap ( )

inline

Definition at line 40 of file HGCalRadiationMap.h.

References doseMap_.

Referenced by HGCalSciNoiseMap::scaleByDose().

{ return doseMap_; }

double HGCalRadiationMap::getDoseValue	(	const int	subdet,
		const int	layer,
		const double	radius,
		bool	logVal = false
	)

Definition at line 28 of file HGCalRadiationMap.cc.

References doseMap_, grayToKrad_, submitPVResolutionJobs::key, log, funct::pow(), alignCSCRings::r, and diffTwoXMLs::r2.

Referenced by HGCalSciNoiseMap::scaleByDose().

                                                                                                          {
  std::pair<int, int> key(subdet, layer);

  if (doseMap_.find(key) == doseMap_.end()) {
    return logVal ? -10. : 0.;
  }

  double r(radius - doseMap_[key].doff_);
  double r2(r * r);
  double r3(r2 * r);
  double r4(r3 * r);

  double cellDoseLog10 =
      doseMap_[key].a_ + doseMap_[key].b_ * r + doseMap_[key].c_ * r2 + doseMap_[key].d_ * r3 + doseMap_[key].e_ * r4;

  return logVal ? cellDoseLog10 * M_LN10 + log(grayToKrad_) : std::pow(10, cellDoseLog10) * grayToKrad_;
}

double HGCalRadiationMap::getFluenceValue	(	const int	subdet,
		const int	layer,
		const double	radius,
		bool	logVal = false
	)

Definition at line 47 of file HGCalRadiationMap.cc.

References doseMap_, fluenceSFlog10_, submitPVResolutionJobs::key, funct::pow(), alignCSCRings::r, and diffTwoXMLs::r2.

Referenced by HGCalSciNoiseMap::scaleByDose().

                                                                                                             {
  std::pair<int, int> key(subdet, layer);

  double r(radius - doseMap_[key].foff_);
  double r2(r * r);
  double r3(r2 * r);
  double r4(r3 * r);

  double cellFluenceLog10 =
      doseMap_[key].f_ + doseMap_[key].g_ * r + doseMap_[key].h_ * r2 + doseMap_[key].i_ * r3 + doseMap_[key].j_ * r4;
  cellFluenceLog10 += fluenceSFlog10_;

  return logVal ? cellFluenceLog10 * M_LN10 : std::pow(10, cellFluenceLog10);
}

std::map< std::pair< int, int >, HGCalRadiationMap::DoseParameters > HGCalRadiationMap::readDosePars ( const std::string &  fullpath )

private

Definition at line 63 of file HGCalRadiationMap.cc.

References HGCalRadiationMap::DoseParameters::a_, HGCalRadiationMap::DoseParameters::b_, HGCalRadiationMap::DoseParameters::c_, HGCalRadiationMap::DoseParameters::d_, HGCalRadiationMap::DoseParameters::doff_, HGCalRadiationMap::DoseParameters::e_, Exception, HGCalRadiationMap::DoseParameters::f_, HGCalRadiationMap::DoseParameters::foff_, personalPlayback::fp, edm::FileInPath::fullPath(), HGCalRadiationMap::DoseParameters::g_, HGCalRadiationMap::DoseParameters::h_, HGCalRadiationMap::DoseParameters::i_, timingPdfMaker::infile, HGCalRadiationMap::DoseParameters::j_, submitPVResolutionJobs::key, phase1PixelTopology::layer, geometryCSVtoXML::line, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by setDoseMap().

                               {
  doseParametersMap result;

  //no dose file means no aging
  if (fullpath.empty())
    return result;

  edm::FileInPath fp(fullpath);
  std::ifstream infile(fp.fullPath());
  if (!infile.is_open()) {
    throw cms::Exception("FileNotFound") << "Unable to open '" << fullpath << "'" << std::endl;
  }
  std::string line;
  while (getline(infile, line)) {
    int subdet;
    int layer;
    DoseParameters dosePars;

    //space-separated
    std::stringstream linestream(line);
    linestream >> subdet >> layer >> dosePars.a_ >> dosePars.b_ >> dosePars.c_ >> dosePars.d_ >> dosePars.e_ >>
        dosePars.doff_ >> dosePars.f_ >> dosePars.g_ >> dosePars.h_ >> dosePars.i_ >> dosePars.j_ >> dosePars.foff_;

    std::pair<int, int> key(subdet, layer);
    result[key] = dosePars;
  }
  return result;
}

void HGCalRadiationMap::setDoseMap	(	const std::string &	fullpath,
		const unsigned int	algo
	)

Definition at line 9 of file HGCalRadiationMap.cc.

References algo(), algo_, doseMap_, and readDosePars().

Referenced by HGCalSciNoiseMap::setDoseMap().

                                                                                     {
  doseMap_ = readDosePars(fullpath);
  algo_ = algo;
}

void HGCalRadiationMap::setFluenceScaleFactor ( double  val )

inline

Definition at line 41 of file HGCalRadiationMap.h.

References fluenceSFlog10_.

Referenced by HGCDigitizerBase::HGCDigitizerBase(), and HGCHEbackDigitizer::HGCHEbackDigitizer().

{ fluenceSFlog10_ = log10(val); }

void HGCalRadiationMap::setGeometry

(

const CaloSubdetectorGeometry *

geom

)

Definition at line 15 of file HGCalRadiationMap.cc.

References HGCalTopology::dddConstants(), geom(), hgcalDDD_, hgcalGeom_, hgcalTopology_, and HGCalGeometry::topology().

Referenced by HGCHEbackDigitizer::runRealisticDigitizer().

                                                                       {
  hgcalGeom_ = static_cast<const HGCalGeometry*>(geom);
  hgcalTopology_ = &(hgcalGeom_->topology());
  hgcalDDD_ = &(hgcalTopology_->dddConstants());
}

const HGCalTopology* HGCalRadiationMap::topo ( )

inline

Definition at line 37 of file HGCalRadiationMap.h.

References hgcalTopology_.

{ return hgcalTopology_; }

Member Data Documentation

unsigned int HGCalRadiationMap::algo_

private

Definition at line 46 of file HGCalRadiationMap.h.

Referenced by algo(), and setDoseMap().

doseParametersMap HGCalRadiationMap::doseMap_

private

Definition at line 50 of file HGCalRadiationMap.h.

Referenced by getDoseMap(), getDoseValue(), getFluenceValue(), and setDoseMap().

double HGCalRadiationMap::fluenceSFlog10_

private

Definition at line 53 of file HGCalRadiationMap.h.

Referenced by getFluenceValue(), and setFluenceScaleFactor().

const double HGCalRadiationMap::grayToKrad_ = 0.1

private

Definition at line 52 of file HGCalRadiationMap.h.

Referenced by getDoseValue().

const HGCalDDDConstants* HGCalRadiationMap::hgcalDDD_

private

Definition at line 49 of file HGCalRadiationMap.h.

Referenced by ddd(), and setGeometry().

const HGCalGeometry* HGCalRadiationMap::hgcalGeom_

private

Definition at line 47 of file HGCalRadiationMap.h.

Referenced by geom(), and setGeometry().

const HGCalTopology* HGCalRadiationMap::hgcalTopology_

private

Definition at line 48 of file HGCalRadiationMap.h.

Referenced by setGeometry(), and topo().