#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 create_public_lumi_plots::year, and years_.

Referenced by HBHEDarkening::LumiYearComp::operator()().

                                                                                                                                                               :
     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 46 of file HBHEDarkening.h.

References degradation().

46 {}

Member Function Documentation

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

Definition at line 88 of file HBHEDarkening.cc.

References funct::abs(), degradationYear(), getYearForLumi(), ieta_shift_, AlCaHLTBitMon_QueryRunRegistry::string, create_public_lumi_plots::year, and years_.

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

                                                                        {
     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.size());
     
     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 72 of file HBHEDarkening.cc.

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

Referenced by degradation(), and get_ieta_shift().

                                                                                                  {
     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 = year.sumlumi_ - intlumi;
     
     //calculate degradation
     return std::exp(-(intlumiToUse*doseToUse)/decayConst);
 }

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

private

Definition at line 47 of file HBHEDarkening.cc.

References dosemaps_.

Referenced by degradationYear(), and get_ieta_shift().

                                                              {
     //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 50 of file HBHEDarkening.h.

References degradationYear(), dose(), getYearForLumi(), ieta_shift_, readDoseMap(), AlCaHLTBitMon_QueryRunRegistry::string, and create_public_lumi_plots::year.

Referenced by HBHERecalibration::setup().

50 { return ieta_shift_; }

HBHEDarkening::ieta_shift_

int ieta_shift_

Definition: HBHEDarkening.h:62

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

private

Definition at line 63 of file HBHEDarkening.cc.

References Exception, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, or, and years_.

Referenced by degradation(), and get_ieta_shift().

                                                            {
     //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 29 of file HBHEDarkening.cc.

References Exception, mps_splice::line, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by get_ieta_shift(), and 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;
 }