#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)

void	calcmuTot ()

double	calcmuTot (double eta)

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)

virtual	~EnergyResolutionVsLumi ()

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() [1/2]

EnergyResolutionVsLumi::EnergyResolutionVsLumi ( )

Definition at line 9 of file EnergyResolutionVsLumi.cc.

                                                {
   m_lumi = 0;
   m_instlumi = 0;
 }

References m_instlumi, and m_lumi.

◆ EnergyResolutionVsLumi() [2/2]

EnergyResolutionVsLumi::EnergyResolutionVsLumi	(	double	lumi,
		double	instlumi
	)

inline

Definition at line 22 of file EnergyResolutionVsLumi.h.

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

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

◆ ~EnergyResolutionVsLumi()

EnergyResolutionVsLumi::~EnergyResolutionVsLumi ( )

virtual

Definition at line 14 of file EnergyResolutionVsLumi.cc.

14 {}

Member Function Documentation

◆ calcampDropPhotoDetector()

double EnergyResolutionVsLumi::calcampDropPhotoDetector ( double eta )

Definition at line 150 of file EnergyResolutionVsLumi.cc.

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

References PVValHelper::eta, muonGEMDigis_cfi::instLumi, m_instlumi, m_lumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

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

◆ calcampDropTotal()

double EnergyResolutionVsLumi::calcampDropTotal ( double eta )

Definition at line 158 of file EnergyResolutionVsLumi.cc.

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

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

Referenced by EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios().

◆ calcampDropTransparency()

double EnergyResolutionVsLumi::calcampDropTransparency ( double eta )

Definition at line 142 of file EnergyResolutionVsLumi.cc.

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

References calcmuEM(), calcmuHD(), ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

Referenced by calcampDropTotal().

◆ calcLightCollectionEfficiencyWeighted()

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

Definition at line 92 of file EnergyResolutionVsLumi.cc.

                                                                                        {
   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;
 }

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

Referenced by ECalSD::curve_LY().

◆ calcLightCollectionEfficiencyWeighted2()

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

Definition at line 123 of file EnergyResolutionVsLumi.cc.

                                                                                                          {
   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;
 }

References calcampDropPhotoDetector(), calcmuTot(), ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, mps_fire::result, findQualityFiles::v, and z.

◆ calcmuEM()

double EnergyResolutionVsLumi::calcmuEM ( double eta )

Definition at line 48 of file EnergyResolutionVsLumi.cc.

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

References PVValHelper::eta, muonGEMDigis_cfi::instLumi, m_instlumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

Referenced by calcampDropTransparency(), and calcresolutitonConstantTerm().

◆ calcmuHD()

double EnergyResolutionVsLumi::calcmuHD ( double eta )

Definition at line 55 of file EnergyResolutionVsLumi.cc.

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

References PVValHelper::eta, m_lumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

Referenced by calcampDropTransparency(), and calcresolutitonConstantTerm().

◆ calcmuTot() [1/2]

void EnergyResolutionVsLumi::calcmuTot ( )

Definition at line 62 of file EnergyResolutionVsLumi.cc.

                                        {
   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;
       }
     }
   }
 }

References funct::abs(), EBDetId::approxEta(), calcampDropPhotoDetector(), PVValHelper::eta, L1TowerCalibrationProducer_cfi::iEta, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, dqm-mbProfile::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().

◆ calcmuTot() [2/2]

double EnergyResolutionVsLumi::calcmuTot ( double eta )

Definition at line 132 of file EnergyResolutionVsLumi.cc.

                                                    {
   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;
 }

References PVValHelper::eta, muonGEMDigis_cfi::instLumi, m_instlumi, m_lumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

◆ calcnoiseADC()

double EnergyResolutionVsLumi::calcnoiseADC ( double eta )

Definition at line 173 of file EnergyResolutionVsLumi.cc.

                                                       {
   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;
 }

References funct::abs(), PVValHelper::eta, m_lumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

◆ calcnoiseIncreaseADC()

double EnergyResolutionVsLumi::calcnoiseIncreaseADC ( double eta )

Definition at line 165 of file EnergyResolutionVsLumi.cc.

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

References PVValHelper::eta, m_lumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

Referenced by EcalTrivialConditionRetriever::produceEcalPedestals().

◆ calcresolutitonConstantTerm()

double EnergyResolutionVsLumi::calcresolutitonConstantTerm ( double eta )

Definition at line 189 of file EnergyResolutionVsLumi.cc.

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

References calcmuEM(), calcmuHD(), ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

Referenced by EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration().

◆ CalculateDegradation()

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

Definition at line 16 of file EnergyResolutionVsLumi.cc.

                                                                                               {
   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;
 }

References PVValHelper::eta, muonGEMDigis_cfi::instLumi, m_instlumi, m_lumi, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and mps_fire::result.

Referenced by Resolution().

◆ getmuTot()

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

◆ Resolution()

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

Definition at line 197 of file EnergyResolutionVsLumi.cc.

                                                                 {
   // 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);
 }

References gen::C, CalculateDegradation(), ztail::d, PVValHelper::eta, N, S(), and mathSSE::sqrt().

◆ setInstLumi()

void EnergyResolutionVsLumi::setInstLumi ( double x )

inline

Definition at line 46 of file EnergyResolutionVsLumi.h.

46 { m_instlumi = x; };

References m_instlumi, and x.

Referenced by EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration(), EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios(), and EcalTrivialConditionRetriever::produceEcalPedestals().

◆ setLumi()

void EnergyResolutionVsLumi::setLumi ( double x )

inline

Definition at line 45 of file EnergyResolutionVsLumi.h.

45 { m_lumi = x; };

References m_lumi, and x.

Referenced by EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration(), EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios(), and EcalTrivialConditionRetriever::produceEcalPedestals().

◆ setLumies()

void EnergyResolutionVsLumi::setLumies	(	double	x,
		double	y
	)

inline

Definition at line 47 of file EnergyResolutionVsLumi.h.

47 { m_lumi = x, m_instlumi = y, calcmuTot(); };

References calcmuTot(), m_instlumi, m_lumi, x, and y.

Referenced by ECalSD::ECalSD().

Member Data Documentation

◆ m_instlumi

double EnergyResolutionVsLumi::m_instlumi

private

Definition at line 67 of file EnergyResolutionVsLumi.h.

Referenced by calcampDropPhotoDetector(), calcmuEM(), calcmuTot(), CalculateDegradation(), EnergyResolutionVsLumi(), setInstLumi(), and setLumies().

◆ m_lumi

double EnergyResolutionVsLumi::m_lumi

private

Definition at line 66 of file EnergyResolutionVsLumi.h.

Referenced by calcampDropPhotoDetector(), calcmuHD(), calcmuTot(), calcnoiseADC(), calcnoiseIncreaseADC(), CalculateDegradation(), EnergyResolutionVsLumi(), setLumi(), and setLumies().

◆ mu_eta

double EnergyResolutionVsLumi::mu_eta[10085]

private

Definition at line 68 of file EnergyResolutionVsLumi.h.

Referenced by calcLightCollectionEfficiencyWeighted(), and calcmuTot().

◆ vpt_eta

double EnergyResolutionVsLumi::vpt_eta[10085]

private

Definition at line 69 of file EnergyResolutionVsLumi.h.

Referenced by calcLightCollectionEfficiencyWeighted(), and calcmuTot().

Classes

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EnergyResolutionVsLumi() [1/2]

◆ EnergyResolutionVsLumi() [2/2]

◆ ~EnergyResolutionVsLumi()

Member Function Documentation

◆ calcampDropPhotoDetector()

◆ calcampDropTotal()

◆ calcampDropTransparency()

◆ calcLightCollectionEfficiencyWeighted()

◆ calcLightCollectionEfficiencyWeighted2()

◆ calcmuEM()

◆ calcmuHD()

◆ calcmuTot() [1/2]

◆ calcmuTot() [2/2]

◆ calcnoiseADC()

◆ calcnoiseIncreaseADC()

◆ calcresolutitonConstantTerm()

◆ CalculateDegradation()

◆ getmuTot()

◆ Resolution()

◆ setInstLumi()

◆ setLumi()

◆ setLumies()

Member Data Documentation

◆ m_instlumi

◆ m_lumi

◆ mu_eta

◆ vpt_eta