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

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::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(), csvLumiCalc::lumi, m_instlumi, m_lumi, and ~EnergyResolutionVsLumi().

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

EnergyResolutionVsLumi::~EnergyResolutionVsLumi ( )

virtual

Definition at line 14 of file EnergyResolutionVsLumi.cc.

Referenced by EnergyResolutionVsLumi().

14 {}

Member Function Documentation

double EnergyResolutionVsLumi::calcampDropPhotoDetector ( double eta )

Definition at line 150 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::AgingVPT(), muonGEMDigis_cfi::instLumi, m_instlumi, m_lumi, and mps_fire::result.

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

                                                                   {
   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 158 of file EnergyResolutionVsLumi.cc.

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

Referenced by EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios(), and setLumies().

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

double EnergyResolutionVsLumi::calcampDropTransparency ( double eta )

Definition at line 142 of file EnergyResolutionVsLumi.cc.

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

Referenced by calcampDropTotal(), and setLumies().

                                                                  {
   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 92 of file EnergyResolutionVsLumi.cc.

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

Referenced by ECalSD::curve_LY(), and setLumies().

                                                                                        {
   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 123 of file EnergyResolutionVsLumi.cc.

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

Referenced by setLumies().

                                                                                                          {
   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 48 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::InducedAbsorptionEM(), muonGEMDigis_cfi::instLumi, m_instlumi, and mps_fire::result.

Referenced by calcampDropTransparency(), calcresolutitonConstantTerm(), and setLumies().

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

double EnergyResolutionVsLumi::calcmuHD ( double eta )

Definition at line 55 of file EnergyResolutionVsLumi.cc.

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

Referenced by calcampDropTransparency(), calcresolutitonConstantTerm(), and setLumies().

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

double EnergyResolutionVsLumi::calcmuTot ( double eta )

Definition at line 132 of file EnergyResolutionVsLumi.cc.

References EvolutionECAL::InducedAbsorptionEM(), EvolutionECAL::InducedAbsorptionHadronic(), muonGEMDigis_cfi::instLumi, m_instlumi, m_lumi, 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 62 of file EnergyResolutionVsLumi.cc.

References funct::abs(), EBDetId::approxEta(), calcampDropPhotoDetector(), PVValHelper::eta, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, cmsBatch::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 173 of file EnergyResolutionVsLumi.cc.

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

Referenced by setLumies().

                                                       {
   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 165 of file EnergyResolutionVsLumi.cc.

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

Referenced by EcalTrivialConditionRetriever::produceEcalPedestals(), and setLumies().

                                                               {
   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 189 of file EnergyResolutionVsLumi.cc.

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

Referenced by EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration(), and setLumies().

                                                                      {
   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 16 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
	)

Referenced by setLumies().

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

Definition at line 197 of file EnergyResolutionVsLumi.cc.

References EnergyResolutionVsLumi::DegradationAtEta::ampDropTotal, patCaloMETCorrections_cff::C, CalculateDegradation(), edmIntegrityCheck::d, N, EnergyResolutionVsLumi::DegradationAtEta::noiseIncreaseADC, EnergyResolutionVsLumi::DegradationAtEta::resolutitonConstantTerm, S(), 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);
 }