#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 24 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 15 of file EnergyResolutionVsLumi.cc.

Referenced by EnergyResolutionVsLumi().

16 {

17 }

Member Function Documentation

double EnergyResolutionVsLumi::calcampDropPhotoDetector ( double eta )

Definition at line 159 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 168 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 150 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 100 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 130 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 52 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 60 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 139 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 68 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 185 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 176 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 203 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 19 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 212 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);
 
 }