#include <HBHEDarkening.h>

Classes
struct	LumiYear

struct	LumiYearComp

Public Member Functions
float	degradation (float intlumi, int ieta, int lay) const

int	get_ieta_shift () const

	HBHEDarkening (int ieta_shift, float drdA, float drdB, const std::map< int, std::vector< std::vector< float >>> &dosemaps, const std::vector< LumiYear > &years)

	~HBHEDarkening ()

Static Public Member Functions
static std::vector < std::vector< float > >	readDoseMap (const std::string &fullpath)

Private Member Functions
float	degradationYear (const LumiYear &year, float intlumi, int ieta, int lay) const

float	dose (int ieta, int lay, int energy) const

std::string	getYearForLumi (float intlumi) const

Private Attributes
std::map< int, std::vector < std::vector< float > > >	dosemaps_

float	drdA_

float	drdB_

int	ieta_shift_

std::vector< LumiYear >	years_

Detailed Description

Definition at line 17 of file HBHEDarkening.h.

Constructor & Destructor Documentation

HBHEDarkening::HBHEDarkening	(	int	ieta_shift,
		float	drdA,
		float	drdB,
		const std::map< int, std::vector< std::vector< float >>> &	dosemaps,
		const std::vector< LumiYear > &	years
	)

Definition at line 16 of file HBHEDarkening.cc.

References years_.

     : ieta_shift_(ieta_shift), drdA_(drdA), drdB_(drdB), dosemaps_(dosemaps), years_(years) {
   //finish initializing years
   std::sort(years_.begin(), years_.end());
   //sum up int lumi
   float sumlumi = 0.0;
   for (auto& year : years_) {
     sumlumi += year.intlumi_;
     year.sumlumi_ = sumlumi;
   }
 }

HBHEDarkening::~HBHEDarkening ( )

inline

Definition at line 45 of file HBHEDarkening.h.

45 {}

Member Function Documentation

float HBHEDarkening::degradation	(	float	intlumi,
		int	ieta,
		int	lay
	)		const

Definition at line 98 of file HBHEDarkening.cc.

References funct::abs(), cms::cuda::assert(), degradationYear(), getYearForLumi(), ieta_shift_, uploadConditions::response, AlCaHLTBitMon_QueryRunRegistry::string, and years_.

Referenced by HcalDigitizer::darkening(), HCalSD::getEnergyDeposit(), and HBHERecalibration::initialize().

                                                                        {
   ieta = abs(ieta);
   //shift ieta tower index to act as array index
   ieta -= ieta_shift_;
   //shift layer index by 1 to act as array index
   lay -= 1;
 
   //accumulate degradation over years
   float response = 1.0;
   std::string yearForLumi = getYearForLumi(intlumi);
   assert(!yearForLumi.empty());
 
   for (const auto& year : years_) {
     response *= degradationYear(year, intlumi, ieta, lay);
     if (year.year_ == yearForLumi)
       break;
   }
 
   return response;
 }

float HBHEDarkening::degradationYear	(	const LumiYear &	year,
		float	intlumi,
		int	ieta,
		int	lay
	)		const

private

Definition at line 80 of file HBHEDarkening.cc.

References dose(), drdA_, drdB_, HBHEDarkening::LumiYear::energy_, funct::exp(), HBHEDarkening::LumiYear::intlumi_, HBHEDarkening::LumiYear::lumirate_, funct::pow(), and HBHEDarkening::LumiYear::sumlumi_.

Referenced by degradation().

                                                                                                  {
   float doseToUse = dose(ieta, lay, year.energy_);
   if (doseToUse == 0.0)
     return 1.0;
 
   //apply dose rate dependence model to the provided year
   //get krad/hr from Mrad/fb-1 and fb-1/hr
   float decayConst = drdA_ * std::pow(1000 * doseToUse * year.lumirate_, drdB_);
 
   //determine if this is a partial year
   float intlumiToUse = year.intlumi_;
   if (intlumi < year.sumlumi_)
     intlumiToUse = intlumi - (year.sumlumi_ - year.intlumi_);
 
   //calculate degradation
   return std::exp(-(intlumiToUse * doseToUse) / decayConst);
 }

float HBHEDarkening::dose	(	int	ieta,
		int	lay,
		int	energy
	)		const

private

Definition at line 51 of file HBHEDarkening.cc.

References dosemaps_.

Referenced by degradationYear().

                                                              {
   //existence check
   const auto dosemapIt = dosemaps_.find(energy);
   if (dosemapIt == dosemaps_.end())
     return 0.0;
 
   //bounds check
   const auto& dosemap = dosemapIt->second;
   if (ieta < 0 or ieta >= int(dosemap.size()))
     return 0.0;
 
   //bounds check
   const auto& doserow = dosemap[ieta];
   if (lay < 0 or lay >= int(doserow.size()))
     return 0.0;
 
   return doserow[lay];
 }

int HBHEDarkening::get_ieta_shift ( ) const

inline

Definition at line 49 of file HBHEDarkening.h.

References ieta_shift_.

Referenced by HBHERecalibration::setup().

49 { return ieta_shift_; }

HBHEDarkening::ieta_shift_

int ieta_shift_

Definition: HBHEDarkening.h:61

std::string HBHEDarkening::getYearForLumi ( float intlumi ) const

private

Definition at line 70 of file HBHEDarkening.cc.

References Exception, cuda_std::lower_bound(), or, and years_.

Referenced by degradation().

                                                            {
   //compare based on sum lumi value
   auto lb = std::lower_bound(years_.begin(), years_.end(), intlumi, LumiYearComp());
   if (lb == years_.end() or lb->sumlumi_ < intlumi) {
     throw cms::Exception("ValueError") << "HBHEDarkening: insufficient LHC run information provided to simulate "
                                        << intlumi << "/fb - check the python config" << std::endl;
   }
   return lb->year_;
 }

std::vector< std::vector< float > > HBHEDarkening::readDoseMap ( const std::string & fullpath )

static

Definition at line 32 of file HBHEDarkening.cc.

References Exception, timingPdfMaker::infile, geometryCSVtoXML::line, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by HBHEDarkeningEP::produce().

                                                                                 {
   std::ifstream infile(fullpath.c_str());
   if (!infile.is_open()) {
     throw cms::Exception("FileNotFound") << "Unable to open '" << fullpath << "'" << std::endl;
   }
   std::string line;
   std::vector<std::vector<float>> result;
   while (getline(infile, line)) {
     //space-separated
     std::stringstream linestream(line);
     std::vector<float> lineresult;
     float doseval;
     while (linestream >> doseval)
       lineresult.push_back(doseval);
     result.push_back(lineresult);
   }
   return result;
 }