#include <EnergyResolutionVsLumi.h>

Classes
struct	DegradationAtEta

Public Member Functions
double	calcampDropPhotoDetector (double eta)

double	calcampDropTotal (double eta)

double	calcampDropTransparency (double eta)

double	calcLightCollectionEfficiencyWeighted (DetId id, double z)

double	calcLightCollectionEfficiencyWeighted2 (double eta, double z, double mu_ind=-1.0)

double	calcmuEM (double eta)

double	calcmuHD (double eta)

double	calcmuTot (double eta)

void	calcmuTot ()

double	calcnoiseADC (double eta)

double	calcnoiseIncreaseADC (double eta)

double	calcresolutitonConstantTerm (double eta)

DegradationAtEta	CalculateDegradation (double eta)

	EnergyResolutionVsLumi ()

	EnergyResolutionVsLumi (double lumi, double instlumi)

double	getmuTot (double eta, int ix, int iy)

double	Resolution (double eta, double ene)

void	setInstLumi (double x)

void	setLumi (double x)

void	setLumies (double x, double y)

	~EnergyResolutionVsLumi ()=default

Private Attributes
double	m_instlumi

double	m_lumi

double	mu_eta [10085]

double	vpt_eta [10085]

Detailed Description

Definition at line 18 of file EnergyResolutionVsLumi.h.

Constructor & Destructor Documentation

EnergyResolutionVsLumi::EnergyResolutionVsLumi ( )

Definition at line 9 of file EnergyResolutionVsLumi.cc.

References m_instlumi, and m_lumi.

                                                {
   m_lumi = 0;
   m_instlumi = 0;
 }

EnergyResolutionVsLumi::EnergyResolutionVsLumi	(	double	lumi,
		double	instlumi
	)

inline

Definition at line 22 of file EnergyResolutionVsLumi.h.

References calcmuTot(), dqmdumpme::lumi, m_instlumi, and m_lumi.

                                                        {
     m_lumi = lumi;
     m_instlumi = instlumi;
     calcmuTot();
   };

EnergyResolutionVsLumi::~EnergyResolutionVsLumi ( )

default

Member Function Documentation

double EnergyResolutionVsLumi::calcampDropPhotoDetector ( double eta )

Definition at line 148 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::AgingVPT(), m_instlumi, m_lumi, isotrackApplyRegressor::model, and mps_fire::result.

Referenced by calcampDropTotal(), calcLightCollectionEfficiencyWeighted2(), and calcmuTot().

                                                                   {
   double instLumi = m_instlumi;
   double totLumi = m_lumi;
   EvolutionECAL model;
   double result = model.AgingVPT(instLumi, totLumi, eta);
   return result;
 }

double EnergyResolutionVsLumi::calcampDropTotal ( double eta )

Definition at line 156 of file EnergyResolutionVsLumi.cc.

References calcampDropPhotoDetector(), calcampDropTransparency(), and mps_fire::result.

Referenced by EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios().

                                                           {
   double tra = this->calcampDropTransparency(eta);
   double pho = this->calcampDropPhotoDetector(eta);
   double result = tra * pho;
   return result;
 }

double EnergyResolutionVsLumi::calcampDropTransparency ( double eta )

Definition at line 140 of file EnergyResolutionVsLumi.cc.

References calcmuEM(), calcmuHD(), EvolutionECAL::DegradationMeanEM50GeV(), isotrackApplyRegressor::model, and mps_fire::result.

Referenced by calcampDropTotal().

                                                                  {
   double muEM = this->calcmuEM(eta);
   double muHD = this->calcmuHD(eta);
   EvolutionECAL model;
   double result = model.DegradationMeanEM50GeV(muEM + muHD);
   return result;
 }

double EnergyResolutionVsLumi::calcLightCollectionEfficiencyWeighted	(	DetId	id,
		double	z
	)

Definition at line 90 of file EnergyResolutionVsLumi.cc.

References funct::abs(), EcalBarrel, EcalEndcap, EBDetId::ieta(), EEDetId::ix(), EEDetId::IX_MAX, EEDetId::iy(), EvolutionECAL::LightCollectionEfficiencyWeighted(), EBDetId::MAX_IETA, isotrackApplyRegressor::model, mu_eta, mps_fire::result, findQualityFiles::v, vpt_eta, and z.

Referenced by ECalSD::curve_LY().

                                                                                        {
   double v = 1.0;
   double muTot = 0;
   if (id.subdetId() == EcalBarrel) {
     EBDetId ebId(id);
     int ieta = std::abs(ebId.ieta());
     muTot = mu_eta[ieta];
 
   } else if (id.subdetId() == EcalEndcap) {
     EEDetId eeId(id);
     int ix = eeId.ix();
     int iy = eeId.iy();
 
     muTot = mu_eta[EBDetId::MAX_IETA + ix + iy * (EEDetId::IX_MAX)];
     v = vpt_eta[EBDetId::MAX_IETA + ix + iy * (EEDetId::IX_MAX)];
   } else {
     muTot = 0;
   }
   double zcor = z;
   EvolutionECAL model;
   if (z < 0.02) {
     zcor = 0.02;
   } else if (z > 0.98) {
     zcor = 0.98;
   }
 
   double result = model.LightCollectionEfficiencyWeighted(zcor, muTot) * v;
 
   return result;
 }

double EnergyResolutionVsLumi::calcLightCollectionEfficiencyWeighted2	(	double	eta,
		double	z,
		double	mu_ind = `-1.0`
	)

Definition at line 121 of file EnergyResolutionVsLumi.cc.

References calcampDropPhotoDetector(), calcmuTot(), EvolutionECAL::LightCollectionEfficiencyWeighted(), isotrackApplyRegressor::model, mps_fire::result, and findQualityFiles::v.

                                                                                                          {
   if (mu_ind < 0)
     mu_ind = this->calcmuTot(eta);
   double v = this->calcampDropPhotoDetector(eta);
   EvolutionECAL model;
   double result = model.LightCollectionEfficiencyWeighted(z, mu_ind) * v;
   return result;
 }

double EnergyResolutionVsLumi::calcmuEM ( double eta )

Definition at line 46 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::InducedAbsorptionEM(), m_instlumi, isotrackApplyRegressor::model, and mps_fire::result.

Referenced by calcampDropTransparency(), and calcresolutitonConstantTerm().

                                                   {
   double instLumi = m_instlumi;
   EvolutionECAL model;
   double result = model.InducedAbsorptionEM(instLumi, eta);
   return result;
 }

double EnergyResolutionVsLumi::calcmuHD ( double eta )

Definition at line 53 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::InducedAbsorptionHadronic(), m_lumi, isotrackApplyRegressor::model, and mps_fire::result.

Referenced by calcampDropTransparency(), and calcresolutitonConstantTerm().

                                                   {
   double totLumi = m_lumi;
   EvolutionECAL model;
   double result = model.InducedAbsorptionHadronic(totLumi, eta);
   return result;
 }

double EnergyResolutionVsLumi::calcmuTot ( double eta )

Definition at line 130 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::InducedAbsorptionEM(), EvolutionECAL::InducedAbsorptionHadronic(), m_instlumi, m_lumi, isotrackApplyRegressor::model, and mps_fire::result.

                                                    {
   double totLumi = m_lumi;
   double instLumi = m_instlumi;
   EvolutionECAL model;
   double muEM = model.InducedAbsorptionEM(instLumi, eta);
   double muH = model.InducedAbsorptionHadronic(totLumi, eta);
   double result = muEM + muH;
   return result;
 }

void EnergyResolutionVsLumi::calcmuTot ( )

Definition at line 60 of file EnergyResolutionVsLumi.cc.

References funct::abs(), EBDetId::approxEta(), calcampDropPhotoDetector(), PVValHelper::eta, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, log, EBDetId::MAX_IETA, mu_eta, alignCSCRings::r, mathSSE::sqrt(), funct::tan(), findQualityFiles::v, EEDetId::validDetId(), and vpt_eta.

Referenced by calcLightCollectionEfficiencyWeighted2(), EnergyResolutionVsLumi(), and setLumies().

                                        {
   for (int iEta = 1; iEta <= EBDetId::MAX_IETA; ++iEta) {
     if (iEta == 0)
       continue;
 
     double eta = EBDetId::approxEta(EBDetId(iEta, 1));
     eta = std::abs(eta);
     double r = calcmuTot(eta);
 
     mu_eta[iEta] = r;
     vpt_eta[iEta] = 1.0;
   }
 
   for (int iX = EEDetId::IX_MIN; iX <= EEDetId::IX_MAX; ++iX) {
     for (int iY = EEDetId::IY_MIN; iY <= EEDetId::IY_MAX; ++iY) {
       if (EEDetId::validDetId(iX, iY, 1)) {
         EEDetId eedetidpos(iX, iY, 1);
         double eta = -log(tan(0.5 * atan(sqrt((iX - 50.0) * (iX - 50.0) + (iY - 50.0) * (iY - 50.0)) * 2.98 / 328.)));
         eta = std::abs(eta);
         double r = calcmuTot(eta);
         double v = calcampDropPhotoDetector(eta);
 
         mu_eta[EBDetId::MAX_IETA + iX + iY * (EEDetId::IX_MAX)] = r;
         vpt_eta[EBDetId::MAX_IETA + iX + iY * (EEDetId::IX_MAX)] = v;
         //std::cout<<"eta/mu/vpt"<<eta<<"/"<<r<<"/"<<v<<std::endl;
       }
     }
   }
 }

double EnergyResolutionVsLumi::calcnoiseADC ( double eta )

Definition at line 171 of file EnergyResolutionVsLumi.cc.

References funct::abs(), m_lumi, isotrackApplyRegressor::model, EvolutionECAL::NoiseFactorFE(), and mps_fire::result.

                                                       {
   double totLumi = m_lumi;
   double Nadc = 1.1;
   double result = 1.0;
   EvolutionECAL model;
   if (std::abs(eta) < 1.497) {
     Nadc = 1.1;
     result = model.NoiseFactorFE(totLumi, eta) * Nadc;
   } else {
     Nadc = 2.0;
     result = Nadc;
     // endcaps no increase in ADC
   }
   return result;
 }

double EnergyResolutionVsLumi::calcnoiseIncreaseADC ( double eta )

Definition at line 163 of file EnergyResolutionVsLumi.cc.

References m_lumi, isotrackApplyRegressor::model, EvolutionECAL::NoiseFactorFE(), and mps_fire::result.

Referenced by EcalTrivialConditionRetriever::produceEcalPedestals().

                                                               {
   double totLumi = m_lumi;
   EvolutionECAL model;
   double result = model.NoiseFactorFE(totLumi, eta);
   return result;
   // noise increase in ADC
 }

double EnergyResolutionVsLumi::calcresolutitonConstantTerm ( double eta )

Definition at line 187 of file EnergyResolutionVsLumi.cc.

References calcmuEM(), calcmuHD(), isotrackApplyRegressor::model, EvolutionECAL::ResolutionConstantTermEM50GeV(), and mps_fire::result.

Referenced by EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration().

                                                                      {
   double muEM = this->calcmuEM(eta);
   double muHD = this->calcmuHD(eta);
   EvolutionECAL model;
   double result = model.ResolutionConstantTermEM50GeV(muEM + muHD);
   return result;
 }

EnergyResolutionVsLumi::DegradationAtEta EnergyResolutionVsLumi::CalculateDegradation ( double eta )

Definition at line 14 of file EnergyResolutionVsLumi.cc.

Referenced by Resolution().

                                                                                               {
   DegradationAtEta result;
 
   result.eta = eta;
   double totLumi = m_lumi;
   double instLumi = m_instlumi;
 
   EvolutionECAL model;
 
   // Index of induced absorption due to EM damages in PWO4
   result.muEM = model.InducedAbsorptionEM(instLumi, eta);
 
   // Index of induced absorption due to hadron damages in PWO4
   result.muHD = model.InducedAbsorptionHadronic(totLumi, eta);
 
   // Average degradation of amplitude due to transparency change
   result.ampDropTransparency = model.DegradationMeanEM50GeV(result.muEM + result.muHD);
 
   // Average degradation of amplitude due to photo-detector aging
   result.ampDropPhotoDetector = model.AgingVPT(instLumi, totLumi, eta);
 
   result.ampDropTotal = result.ampDropTransparency * result.ampDropPhotoDetector;
 
   // Noise increase in ADC counts due to photo-detector and front-end
   result.noiseIncreaseADC = model.NoiseFactorFE(totLumi, eta);
 
   // Resolution degradation due to LY non-uniformity caused by transparency loss
   result.resolutitonConstantTerm = model.ResolutionConstantTermEM50GeV(result.muEM + result.muHD);
 
   return result;
 }

double EnergyResolutionVsLumi::getmuTot	(	double	eta,
		int	ix,
		int	iy
	)

double EnergyResolutionVsLumi::Resolution	(	double	eta,
		double	ene
	)

Definition at line 195 of file EnergyResolutionVsLumi.cc.

References EnergyResolutionVsLumi::DegradationAtEta::ampDropTotal, gen::C, CalculateDegradation(), ztail::d, N, EnergyResolutionVsLumi::DegradationAtEta::noiseIncreaseADC, EnergyResolutionVsLumi::DegradationAtEta::resolutitonConstantTerm, and mathSSE::sqrt().

                                                                 {
   // Initial parameters for ECAL resolution
   double S;
   double Nadc;
   double adc2GeV;
   double C;
   if (eta < 1.497) {
     S = 0.028;  // CMS note 2006/148 (EB test beam)
     Nadc = 1.1;
     adc2GeV = 0.039;
     C = 0.003;
   } else {
     S = 0.052;  //  CMS DN 2009/002
     Nadc = 2.0;
     adc2GeV = 0.069;
     C = 0.004;
   }
 
   DegradationAtEta d = CalculateDegradation(eta);
 
   // adjust resolution parameters
   S /= sqrt(d.ampDropTotal);
   Nadc *= d.noiseIncreaseADC;
   adc2GeV /= d.ampDropTotal;
   double N = Nadc * adc2GeV * 3.0;  // 3x3 noise in GeV
   C = sqrt(C * C + d.resolutitonConstantTerm * d.resolutitonConstantTerm);
 
   return sqrt(S * S / ene + N * N / ene / ene + C * C);
 }