#include <EGEnergyCorrector.h>

Public Member Functions
std::pair< double, double >	CorrectedEnergyWithError (const reco::Photon &p, const reco::VertexCollection &vtxcol, EcalClusterLazyTools &clustertools, const edm::EventSetup &es)

std::pair< double, double >	CorrectedEnergyWithError (const reco::GsfElectron &e, const reco::VertexCollection &vtxcol, EcalClusterLazyTools &clustertools, const edm::EventSetup &es)

std::pair< double, double >	CorrectedEnergyWithErrorV3 (const reco::Photon &p, const reco::VertexCollection &vtxcol, double rho, EcalClusterLazyTools &clustertools, const edm::EventSetup &es, bool applyRescale=false)

std::pair< double, double >	CorrectedEnergyWithErrorV3 (const reco::GsfElectron &e, const reco::VertexCollection &vtxcol, double rho, EcalClusterLazyTools &clustertools, const edm::EventSetup &es)

	EGEnergyCorrector ()

void	Initialize (const edm::EventSetup &iSetup, std::string regweights, bool weightsFromDB=false)

Bool_t	IsInitialized () const

	~EGEnergyCorrector ()

Protected Attributes
EcalClusterLocal	_ecalLocal

Bool_t	fIsInitialized

Bool_t	fOwnsForests

const GBRForest *	fReadereb

const GBRForest *	fReaderebvariance

const GBRForest *	fReaderee

const GBRForest *	fReadereevariance

Float_t *	fVals

Detailed Description

Definition at line 22 of file EGEnergyCorrector.h.

Constructor & Destructor Documentation

EGEnergyCorrector::EGEnergyCorrector ( )

Definition at line 16 of file EGEnergyCorrector.cc.

                                      :
 fReadereb(nullptr),
 fReaderebvariance(nullptr),
 fReaderee(nullptr),
 fReadereevariance(nullptr),
 fIsInitialized(kFALSE),
 fOwnsForests(kFALSE),
 fVals(nullptr)
 {
   // Constructor.
 }

EGEnergyCorrector::~EGEnergyCorrector ( )

Definition at line 30 of file EGEnergyCorrector.cc.

References fOwnsForests, fReadereb, fReaderebvariance, fReaderee, fReadereevariance, and fVals.

 {
 
   if (fVals) delete [] fVals;
 
   if (fOwnsForests) {
     if (fReadereb) delete fReadereb;
     if (fReaderebvariance) delete fReaderebvariance;
     if (fReaderee) delete fReaderee;
     if (fReadereevariance) delete fReadereevariance;
   }
 
 }

Member Function Documentation

std::pair< double, double > EGEnergyCorrector::CorrectedEnergyWithError	(	const reco::Photon &	p,
		const reco::VertexCollection &	vtxcol,
		EcalClusterLazyTools &	clustertools,
		const edm::EventSetup &	es
	)

Definition at line 90 of file EGEnergyCorrector.cc.

References _ecalLocal, Abs(), b, reco::deltaPhi(), reco::Photon::e5x5(), EcalBarrel, reco::CaloCluster::energy(), reco::CaloCluster::eta(), fReadereb, fReaderebvariance, fReaderee, fReadereevariance, fVals, GBRForest::GetResponse(), reco::Photon::hadronicOverEm(), reco::CaloCluster::hitsAndFractions(), edm::Ptr< T >::isNonnull(), EcalClusterLocal::localCoordsEB(), cmsBatch::log, reco::CaloCluster::phi(), reco::Photon::r9(), matplotRender::reader, alignCSCRings::s, mathSSE::sqrt(), and reco::Photon::superCluster().

Referenced by EGEnergyAnalyzer::analyze(), and IsInitialized().

                                                                                                                                                                                    {
 
   const SuperClusterRef s = p.superCluster();
   const CaloClusterPtr b = s->seed(); //seed  basic cluster
 
   //highest energy basic cluster excluding seed basic cluster
   CaloClusterPtr b2;
   Double_t ebcmax = -99.;
   for (reco::CaloCluster_iterator bit = s->clustersBegin(); bit!=s->clustersEnd(); ++bit) {
     const CaloClusterPtr bc = *bit;
     if (bc->energy() > ebcmax && bc !=b) {
       b2 = bc;
       ebcmax = bc->energy();
     }
   }
 
   //lowest energy basic cluster excluding seed (for pileup mitigation)
   CaloClusterPtr bclast;
   Double_t ebcmin = 1e6;
   for (reco::CaloCluster_iterator bit = s->clustersBegin(); bit!=s->clustersEnd(); ++bit) {
     const CaloClusterPtr bc = *bit;
     if (bc->energy() < ebcmin && bc !=b) {
       bclast = bc;
       ebcmin = bc->energy();
     }
   }
 
   //2nd lowest energy basic cluster excluding seed (for pileup mitigation)
   CaloClusterPtr bclast2;
   ebcmin = 1e6;
   for (reco::CaloCluster_iterator bit = s->clustersBegin(); bit!=s->clustersEnd(); ++bit) {
     const CaloClusterPtr bc = *bit;
     if (bc->energy() < ebcmin && bc !=b && bc!=bclast) {
       bclast2 = bc;
       ebcmin = bc->energy();
     }
   }
 
   Bool_t isbarrel =  b->hitsAndFractions().at(0).first.subdetId()==EcalBarrel;
   Bool_t hasbc2 = b2.isNonnull() && b2->energy()>0.;
   Bool_t hasbclast = bclast.isNonnull() && bclast->energy()>0.;
   Bool_t hasbclast2 = bclast2.isNonnull() && bclast2->energy()>0.;
 
 
   if (isbarrel) {
 
     //basic supercluster variables
     fVals[0]  = s->rawEnergy();
     fVals[1]  = p.r9();
     fVals[2]  = s->eta();
     fVals[3]  = s->phi();
     fVals[4]  = p.e5x5()/s->rawEnergy();
     fVals[5] = p.hadronicOverEm();
     fVals[6] = s->etaWidth();
     fVals[7] = s->phiWidth();
 
 
     //seed basic cluster variables
     double bemax = clustertools.eMax(*b);
     double be2nd = clustertools.e2nd(*b);
     double betop = clustertools.eTop(*b);
     double bebottom = clustertools.eBottom(*b);
     double beleft = clustertools.eLeft(*b);
     double beright = clustertools.eRight(*b);
 
 
     fVals[8] = b->eta()-s->eta();
     fVals[9] = reco::deltaPhi(b->phi(),s->phi());
     fVals[10] = b->energy()/s->rawEnergy();
     fVals[11] = clustertools.e3x3(*b)/b->energy();
     fVals[12] = clustertools.e5x5(*b)/b->energy();
     fVals[13] = sqrt(clustertools.localCovariances(*b)[0]); //sigietaieta
     fVals[14] = sqrt(clustertools.localCovariances(*b)[2]); //sigiphiiphi
     fVals[15] = clustertools.localCovariances(*b)[1];       //sigietaiphi
     fVals[16] = bemax/b->energy();                       //crystal energy ratio gap variables
     fVals[17] = log(be2nd/bemax);
     fVals[18] = log(betop/bemax);
     fVals[19] = log(bebottom/bemax);
     fVals[20] = log(beleft/bemax);
     fVals[21] = log(beright/bemax);
     fVals[22] = (betop-bebottom)/(betop+bebottom);
     fVals[23] = (beleft-beright)/(beleft+beright);
 
 
     double bc2emax = hasbc2 ? clustertools.eMax(*b2) : 0.;
     double bc2e2nd = hasbc2 ? clustertools.e2nd(*b2) : 0.;
     double bc2etop = hasbc2 ? clustertools.eTop(*b2) : 0.;
     double bc2ebottom = hasbc2 ? clustertools.eBottom(*b2) : 0.;
     double bc2eleft = hasbc2 ? clustertools.eLeft(*b2) : 0.;
     double bc2eright = hasbc2 ? clustertools.eRight(*b2) : 0.;
 
     fVals[24] = hasbc2 ? (b2->eta()-s->eta()) : 0.;
     fVals[25] = hasbc2 ? reco::deltaPhi(b2->phi(),s->phi()) : 0.;
     fVals[26] = hasbc2 ? b2->energy()/s->rawEnergy() : 0.;
     fVals[27] = hasbc2 ? clustertools.e3x3(*b2)/b2->energy() : 0.;
     fVals[28] = hasbc2 ? clustertools.e5x5(*b2)/b2->energy() : 0.;
     fVals[29] = hasbc2 ? sqrt(clustertools.localCovariances(*b2)[0]) : 0.;
     fVals[30] = hasbc2 ? sqrt(clustertools.localCovariances(*b2)[2]) : 0.;
     fVals[31] = hasbc2 ? clustertools.localCovariances(*b)[1] : 0.;
     fVals[32] = hasbc2 ? bc2emax/b2->energy() : 0.;
     fVals[33] = hasbc2 ? log(bc2e2nd/bc2emax) : 0.;
     fVals[34] = hasbc2 ? log(bc2etop/bc2emax) : 0.;
     fVals[35] = hasbc2 ? log(bc2ebottom/bc2emax) : 0.;
     fVals[36] = hasbc2 ? log(bc2eleft/bc2emax) : 0.;
     fVals[37] = hasbc2 ? log(bc2eright/bc2emax) : 0.;
     fVals[38] = hasbc2 ? (bc2etop-bc2ebottom)/(bc2etop+bc2ebottom) : 0.;
     fVals[39] = hasbc2 ? (bc2eleft-bc2eright)/(bc2eleft+bc2eright) : 0.;
 
     fVals[40] = hasbclast ? (bclast->eta()-s->eta()) : 0.;
     fVals[41] = hasbclast ? reco::deltaPhi(bclast->phi(),s->phi()) : 0.;
     fVals[42] = hasbclast ? bclast->energy()/s->rawEnergy() : 0.;
     fVals[43] = hasbclast ? clustertools.e3x3(*bclast)/bclast->energy() : 0.;
     fVals[44] = hasbclast ? clustertools.e5x5(*bclast)/bclast->energy() : 0.;
     fVals[45] = hasbclast ? sqrt(clustertools.localCovariances(*bclast)[0]) : 0.;
     fVals[46] = hasbclast ? sqrt(clustertools.localCovariances(*bclast)[2]) : 0.;
     fVals[47] = hasbclast ? clustertools.localCovariances(*bclast)[1] : 0.;
 
     fVals[48] = hasbclast2 ? (bclast2->eta()-s->eta()) : 0.;
     fVals[49] = hasbclast2 ? reco::deltaPhi(bclast2->phi(),s->phi()) : 0.;
     fVals[50] = hasbclast2 ? bclast2->energy()/s->rawEnergy() : 0.;
     fVals[51] = hasbclast2 ? clustertools.e3x3(*bclast2)/bclast2->energy() : 0.;
     fVals[52] = hasbclast2 ? clustertools.e5x5(*bclast2)/bclast2->energy() : 0.;
     fVals[53] = hasbclast2 ? sqrt(clustertools.localCovariances(*bclast2)[0]) : 0.;
     fVals[54] = hasbclast2 ? sqrt(clustertools.localCovariances(*bclast2)[2]) : 0.;
     fVals[55] = hasbclast2 ? clustertools.localCovariances(*bclast2)[1] : 0.;
 
 
     //local coordinates and crystal indices
 
 
     //seed cluster
     float betacry, bphicry, bthetatilt, bphitilt;
     int bieta, biphi;
     _ecalLocal.localCoordsEB(*b,es,betacry,bphicry,bieta,biphi,bthetatilt,bphitilt);
 
     fVals[56] = bieta; //crystal ieta
     fVals[57] = biphi; //crystal iphi
     fVals[58] = bieta%5; //submodule boundary eta symmetry
     fVals[59] = biphi%2; //submodule boundary phi symmetry
     fVals[60] = (TMath::Abs(bieta)<=25)*(bieta%25) + (TMath::Abs(bieta)>25)*((bieta-25*TMath::Abs(bieta)/bieta)%20);  //module boundary eta approximate symmetry
     fVals[61] = biphi%20; //module boundary phi symmetry
     fVals[62] = betacry; //local coordinates with respect to closest crystal center at nominal shower depth
     fVals[63] = bphicry;
 
 
     //2nd cluster (meaningful gap corrections for converted photons)
     float bc2etacry, bc2phicry, bc2thetatilt, bc2phitilt;
     int bc2ieta, bc2iphi;
     if (hasbc2) _ecalLocal.localCoordsEB(*b2,es,bc2etacry,bc2phicry,bc2ieta,bc2iphi,bc2thetatilt,bc2phitilt);
 
     fVals[64] = hasbc2 ? bc2ieta : 0.;
     fVals[65] = hasbc2 ? bc2iphi : 0.;
     fVals[66] = hasbc2 ? bc2ieta%5 : 0.;
     fVals[67] = hasbc2 ? bc2iphi%2 : 0.;
     fVals[68] = hasbc2 ? (TMath::Abs(bc2ieta)<=25)*(bc2ieta%25) + (TMath::Abs(bc2ieta)>25)*((bc2ieta-25*TMath::Abs(bc2ieta)/bc2ieta)%20) : 0.;
     fVals[69] = hasbc2 ? bc2iphi%20 : 0.;
     fVals[70] = hasbc2 ? bc2etacry : 0.;
     fVals[71] = hasbc2 ? bc2phicry : 0.;
 
     fVals[72] = vtxcol.size();
 
   }
   else {
     fVals[0]  = s->rawEnergy();
     fVals[1]  = p.r9();
     fVals[2]  = s->eta();
     fVals[3]  = s->phi();
     fVals[4]  = p.e5x5()/s->rawEnergy();
     fVals[5] = s->etaWidth();
     fVals[6] = s->phiWidth();
     fVals[7] = vtxcol.size();
   }
 
 
   const Double_t varscale = 1.253;
   Double_t den;
   const GBRForest *reader;
   const GBRForest *readervar;
   if (isbarrel) {
     den = s->rawEnergy();
     reader = fReadereb;
     readervar = fReaderebvariance;
   }
   else {
     den = s->rawEnergy() + s->preshowerEnergy();
     reader = fReaderee;
     readervar = fReadereevariance;
   }
 
   Double_t ecor = reader->GetResponse(fVals)*den;
   Double_t ecorerr = readervar->GetResponse(fVals)*den*varscale;
 
   //printf("ecor = %5f, ecorerr = %5f\n",ecor,ecorerr);
 
   return std::pair<double,double>(ecor,ecorerr);
 }

std::pair< double, double > EGEnergyCorrector::CorrectedEnergyWithError	(	const reco::GsfElectron &	e,
		const reco::VertexCollection &	vtxcol,
		EcalClusterLazyTools &	clustertools,
		const edm::EventSetup &	es
	)

Definition at line 289 of file EGEnergyCorrector.cc.

References _ecalLocal, Abs(), b, reco::deltaPhi(), EcalBarrel, reco::GsfElectron::ecalDrivenSeed(), reco::CaloCluster::energy(), reco::CaloCluster::eta(), fReadereb, fReaderebvariance, fReaderee, fReadereevariance, fVals, GBRForest::GetResponse(), reco::GsfElectron::hcalOverEcal(), reco::CaloCluster::hitsAndFractions(), edm::Ptr< T >::isNonnull(), EcalClusterLocal::localCoordsEB(), cmsBatch::log, reco::CaloCluster::phi(), matplotRender::reader, alignCSCRings::s, mathSSE::sqrt(), and reco::GsfElectron::superCluster().

                                                                                                                                                                                         {
 
   //apply v2 regression to electrons
   //mostly duplicated from photon function above //TODO, make common underlying function
 
   //protection, this doesn't work properly on non-egamma-seeded electrons
   if (!e.ecalDrivenSeed()) return std::pair<double,double>(0.,0.);
 
 
   const SuperClusterRef s = e.superCluster();
   const CaloClusterPtr b = s->seed();
 
   CaloClusterPtr b2;
   Double_t ebcmax = -99.;
   for (reco::CaloCluster_iterator bit = s->clustersBegin(); bit!=s->clustersEnd(); ++bit) {
     const CaloClusterPtr bc = *bit;
     if (bc->energy() > ebcmax && bc !=b) {
       b2 = bc;
       ebcmax = bc->energy();
     }
   }
 
   CaloClusterPtr bclast;
   Double_t ebcmin = 1e6;
   for (reco::CaloCluster_iterator bit = s->clustersBegin(); bit!=s->clustersEnd(); ++bit) {
     const CaloClusterPtr bc = *bit;
     if (bc->energy() < ebcmin && bc !=b) {
       bclast = bc;
       ebcmin = bc->energy();
     }
   }
 
   CaloClusterPtr bclast2;
   ebcmin = 1e6;
   for (reco::CaloCluster_iterator bit = s->clustersBegin(); bit!=s->clustersEnd(); ++bit) {
     const CaloClusterPtr bc = *bit;
     if (bc->energy() < ebcmin && bc !=b && bc!=bclast) {
       bclast2 = bc;
       ebcmin = bc->energy();
     }
   }
 
   Bool_t isbarrel =  b->hitsAndFractions().at(0).first.subdetId()==EcalBarrel;
   Bool_t hasbc2 = b2.isNonnull() && b2->energy()>0.;
   Bool_t hasbclast = bclast.isNonnull() && bclast->energy()>0.;
   Bool_t hasbclast2 = bclast2.isNonnull() && bclast2->energy()>0.;
 
 
   if (isbarrel) {
     fVals[0]  = s->rawEnergy();
     fVals[1]  = clustertools.e3x3(*b)/s->rawEnergy(); //r9
     fVals[2]  = s->eta();
     fVals[3]  = s->phi();
     fVals[4]  = clustertools.e5x5(*b)/s->rawEnergy();
     fVals[5] = e.hcalOverEcal();
     fVals[6] = s->etaWidth();
     fVals[7] = s->phiWidth();
 
 
     double bemax = clustertools.eMax(*b);
     double be2nd = clustertools.e2nd(*b);
     double betop = clustertools.eTop(*b);
     double bebottom = clustertools.eBottom(*b);
     double beleft = clustertools.eLeft(*b);
     double beright = clustertools.eRight(*b);
 
 
     fVals[8] = b->eta()-s->eta();
     fVals[9] = reco::deltaPhi(b->phi(),s->phi());
     fVals[10] = b->energy()/s->rawEnergy();
     fVals[11] = clustertools.e3x3(*b)/b->energy();
     fVals[12] = clustertools.e5x5(*b)/b->energy();
     fVals[13] = sqrt(clustertools.localCovariances(*b)[0]);
     fVals[14] = sqrt(clustertools.localCovariances(*b)[2]);
     fVals[15] = clustertools.localCovariances(*b)[1];
     fVals[16] = bemax/b->energy();
     fVals[17] = log(be2nd/bemax);
     fVals[18] = log(betop/bemax);
     fVals[19] = log(bebottom/bemax);
     fVals[20] = log(beleft/bemax);
     fVals[21] = log(beright/bemax);
     fVals[22] = (betop-bebottom)/(betop+bebottom);
     fVals[23] = (beleft-beright)/(beleft+beright);
 
 
     double bc2emax = hasbc2 ? clustertools.eMax(*b2) : 0.;
     double bc2e2nd = hasbc2 ? clustertools.e2nd(*b2) : 0.;
     double bc2etop = hasbc2 ? clustertools.eTop(*b2) : 0.;
     double bc2ebottom = hasbc2 ? clustertools.eBottom(*b2) : 0.;
     double bc2eleft = hasbc2 ? clustertools.eLeft(*b2) : 0.;
     double bc2eright = hasbc2 ? clustertools.eRight(*b2) : 0.;
 
     fVals[24] = hasbc2 ? (b2->eta()-s->eta()) : 0.;
     fVals[25] = hasbc2 ? reco::deltaPhi(b2->phi(),s->phi()) : 0.;
     fVals[26] = hasbc2 ? b2->energy()/s->rawEnergy() : 0.;
     fVals[27] = hasbc2 ? clustertools.e3x3(*b2)/b2->energy() : 0.;
     fVals[28] = hasbc2 ? clustertools.e5x5(*b2)/b2->energy() : 0.;
     fVals[29] = hasbc2 ? sqrt(clustertools.localCovariances(*b2)[0]) : 0.;
     fVals[30] = hasbc2 ? sqrt(clustertools.localCovariances(*b2)[2]) : 0.;
     fVals[31] = hasbc2 ? clustertools.localCovariances(*b)[1] : 0.;
     fVals[32] = hasbc2 ? bc2emax/b2->energy() : 0.;
     fVals[33] = hasbc2 ? log(bc2e2nd/bc2emax) : 0.;
     fVals[34] = hasbc2 ? log(bc2etop/bc2emax) : 0.;
     fVals[35] = hasbc2 ? log(bc2ebottom/bc2emax) : 0.;
     fVals[36] = hasbc2 ? log(bc2eleft/bc2emax) : 0.;
     fVals[37] = hasbc2 ? log(bc2eright/bc2emax) : 0.;
     fVals[38] = hasbc2 ? (bc2etop-bc2ebottom)/(bc2etop+bc2ebottom) : 0.;
     fVals[39] = hasbc2 ? (bc2eleft-bc2eright)/(bc2eleft+bc2eright) : 0.;
 
     fVals[40] = hasbclast ? (bclast->eta()-s->eta()) : 0.;
     fVals[41] = hasbclast ? reco::deltaPhi(bclast->phi(),s->phi()) : 0.;
     fVals[42] = hasbclast ? bclast->energy()/s->rawEnergy() : 0.;
     fVals[43] = hasbclast ? clustertools.e3x3(*bclast)/bclast->energy() : 0.;
     fVals[44] = hasbclast ? clustertools.e5x5(*bclast)/bclast->energy() : 0.;
     fVals[45] = hasbclast ? sqrt(clustertools.localCovariances(*bclast)[0]) : 0.;
     fVals[46] = hasbclast ? sqrt(clustertools.localCovariances(*bclast)[2]) : 0.;
     fVals[47] = hasbclast ? clustertools.localCovariances(*bclast)[1] : 0.;
 
     fVals[48] = hasbclast2 ? (bclast2->eta()-s->eta()) : 0.;
     fVals[49] = hasbclast2 ? reco::deltaPhi(bclast2->phi(),s->phi()) : 0.;
     fVals[50] = hasbclast2 ? bclast2->energy()/s->rawEnergy() : 0.;
     fVals[51] = hasbclast2 ? clustertools.e3x3(*bclast2)/bclast2->energy() : 0.;
     fVals[52] = hasbclast2 ? clustertools.e5x5(*bclast2)/bclast2->energy() : 0.;
     fVals[53] = hasbclast2 ? sqrt(clustertools.localCovariances(*bclast2)[0]) : 0.;
     fVals[54] = hasbclast2 ? sqrt(clustertools.localCovariances(*bclast2)[2]) : 0.;
     fVals[55] = hasbclast2 ? clustertools.localCovariances(*bclast2)[1] : 0.;
 
 
     float betacry, bphicry, bthetatilt, bphitilt;
     int bieta, biphi;
     _ecalLocal.localCoordsEB(*b,es,betacry,bphicry,bieta,biphi,bthetatilt,bphitilt);
 
     fVals[56] = bieta;
     fVals[57] = biphi;
     fVals[58] = bieta%5;
     fVals[59] = biphi%2;
     fVals[60] = (TMath::Abs(bieta)<=25)*(bieta%25) + (TMath::Abs(bieta)>25)*((bieta-25*TMath::Abs(bieta)/bieta)%20);
     fVals[61] = biphi%20;
     fVals[62] = betacry;
     fVals[63] = bphicry;
 
     float bc2etacry, bc2phicry, bc2thetatilt, bc2phitilt;
     int bc2ieta, bc2iphi;
     if (hasbc2) _ecalLocal.localCoordsEB(*b2,es,bc2etacry,bc2phicry,bc2ieta,bc2iphi,bc2thetatilt,bc2phitilt);
 
     fVals[64] = hasbc2 ? bc2ieta : 0.;
     fVals[65] = hasbc2 ? bc2iphi : 0.;
     fVals[66] = hasbc2 ? bc2ieta%5 : 0.;
     fVals[67] = hasbc2 ? bc2iphi%2 : 0.;
     fVals[68] = hasbc2 ? (TMath::Abs(bc2ieta)<=25)*(bc2ieta%25) + (TMath::Abs(bc2ieta)>25)*((bc2ieta-25*TMath::Abs(bc2ieta)/bc2ieta)%20) : 0.;
     fVals[69] = hasbc2 ? bc2iphi%20 : 0.;
     fVals[70] = hasbc2 ? bc2etacry : 0.;
     fVals[71] = hasbc2 ? bc2phicry : 0.;
 
     fVals[72] = vtxcol.size();
 
   }
   else {
     fVals[0]  = s->rawEnergy();
     fVals[1]  = clustertools.e3x3(*b)/s->rawEnergy(); //r9
     fVals[2]  = s->eta();
     fVals[3]  = s->phi();
     fVals[4]  = clustertools.e5x5(*b)/s->rawEnergy();
     fVals[5] = s->etaWidth();
     fVals[6] = s->phiWidth();
     fVals[7] = vtxcol.size();
   }
 
 
   const Double_t varscale = 1.253;
   Double_t den;
   const GBRForest *reader;
   const GBRForest *readervar;
   if (isbarrel) {
     den = s->rawEnergy();
     reader = fReadereb;
     readervar = fReaderebvariance;
   }
   else {
     den = s->rawEnergy() + s->preshowerEnergy();
     reader = fReaderee;
     readervar = fReadereevariance;
   }
 
   Double_t ecor = reader->GetResponse(fVals)*den;
   Double_t ecorerr = readervar->GetResponse(fVals)*den*varscale;
 
   //printf("ecor = %5f, ecorerr = %5f\n",ecor,ecorerr);
 
   return std::pair<double,double>(ecor,ecorerr);
 
 }

std::pair< double, double > EGEnergyCorrector::CorrectedEnergyWithErrorV3	(	const reco::Photon &	p,
		const reco::VertexCollection &	vtxcol,
		double	rho,
		EcalClusterLazyTools &	clustertools,
		const edm::EventSetup &	es,
		bool	applyRescale = `false`
	)

Definition at line 483 of file EGEnergyCorrector.cc.

References _ecalLocal, Abs(), b, reco::deltaPhi(), reco::Photon::e5x5(), EcalBarrel, reco::CaloCluster::energy(), reco::CaloCluster::eta(), fReadereb, fReaderebvariance, fReaderee, fReadereevariance, fVals, GBRForest::GetResponse(), reco::Photon::hadTowOverEm(), reco::CaloCluster::hitsAndFractions(), EcalClusterLocal::localCoordsEB(), Min(), reco::CaloCluster::phi(), reco::Photon::r9(), matplotRender::reader, rho, alignCSCRings::s, mathSSE::sqrt(), and reco::Photon::superCluster().

Referenced by IsInitialized().

                                                                                                                                                                                                                     {
 
   const SuperClusterRef s = p.superCluster();
   const CaloClusterPtr b = s->seed(); //seed  basic cluster
 
   Bool_t isbarrel =  b->hitsAndFractions().at(0).first.subdetId()==EcalBarrel;
 
   //basic supercluster variables
   fVals[0]  = s->rawEnergy();
   fVals[1]  = s->eta();
   fVals[2]  = s->phi();
   fVals[3]  = p.r9();
   fVals[4]  = p.e5x5()/s->rawEnergy();
   fVals[5] = s->etaWidth();
   fVals[6] = s->phiWidth();
   fVals[7] = s->clustersSize();
   fVals[8] = p.hadTowOverEm();
   fVals[9] = rho;
   fVals[10] = vtxcol.size();
 
   //seed basic cluster variables
   double bemax = clustertools.eMax(*b);
   double be2nd = clustertools.e2nd(*b);
   double betop = clustertools.eTop(*b);
   double bebottom = clustertools.eBottom(*b);
   double beleft = clustertools.eLeft(*b);
   double beright = clustertools.eRight(*b);
 
   double be2x5max = clustertools.e2x5Max(*b);
   double be2x5top = clustertools.e2x5Top(*b);
   double be2x5bottom = clustertools.e2x5Bottom(*b);
   double be2x5left = clustertools.e2x5Left(*b);
   double be2x5right = clustertools.e2x5Right(*b);
 
   fVals[11] = b->eta()-s->eta();
   fVals[12] = reco::deltaPhi(b->phi(),s->phi());
   fVals[13] = b->energy()/s->rawEnergy();
   fVals[14] = clustertools.e3x3(*b)/b->energy();
   fVals[15] = clustertools.e5x5(*b)/b->energy();
   fVals[16] = sqrt(clustertools.localCovariances(*b)[0]); //sigietaieta
   fVals[17] = sqrt(clustertools.localCovariances(*b)[2]); //sigiphiiphi
   fVals[18] = clustertools.localCovariances(*b)[1];       //sigietaiphi
   fVals[19] = bemax/b->energy();                       //crystal energy ratio gap variables
   fVals[20] = be2nd/b->energy();
   fVals[21] = betop/b->energy();
   fVals[22] = bebottom/b->energy();
   fVals[23] = beleft/b->energy();
   fVals[24] = beright/b->energy();
   fVals[25] = be2x5max/b->energy();                       //crystal energy ratio gap variables
   fVals[26] = be2x5top/b->energy();
   fVals[27] = be2x5bottom/b->energy();
   fVals[28] = be2x5left/b->energy();
   fVals[29] = be2x5right/b->energy();
 
   if (isbarrel) {
     //local coordinates and crystal indices (barrel only)
 
     //seed cluster
     float betacry, bphicry, bthetatilt, bphitilt;
     int bieta, biphi;
     _ecalLocal.localCoordsEB(*b,es,betacry,bphicry,bieta,biphi,bthetatilt,bphitilt);
 
     fVals[30] = bieta; //crystal ieta
     fVals[31] = biphi; //crystal iphi
     fVals[32] = bieta%5; //submodule boundary eta symmetry
     fVals[33] = biphi%2; //submodule boundary phi symmetry
     fVals[34] = (TMath::Abs(bieta)<=25)*(bieta%25) + (TMath::Abs(bieta)>25)*((bieta-25*TMath::Abs(bieta)/bieta)%20);  //module boundary eta approximate symmetry
     fVals[35] = biphi%20; //module boundary phi symmetry
     fVals[36] = betacry; //local coordinates with respect to closest crystal center at nominal shower depth
     fVals[37] = bphicry;
 
   }
   else {
     //preshower energy ratio (endcap only)
     fVals[30]  = s->preshowerEnergy()/s->rawEnergy();
   }
 
 //   if (isbarrel) {
 //     for (int i=0; i<38; ++i) printf("%i: %5f\n",i,fVals[i]);
 //   }
 //   else for (int i=0; i<31; ++i) printf("%i: %5f\n",i,fVals[i]);
 
   Double_t den;
   const GBRForest *reader;
   const GBRForest *readervar;
   if (isbarrel) {
     den = s->rawEnergy();
     reader = fReadereb;
     readervar = fReaderebvariance;
   }
   else {
     den = s->rawEnergy() + s->preshowerEnergy();
     reader = fReaderee;
     readervar = fReadereevariance;
   }
 
   Double_t ecor = reader->GetResponse(fVals)*den;
 
 
   //apply shower shape rescaling - for Monte Carlo only, and only for calculation of energy uncertainty
   if (applyRescale) {
     if (isbarrel) {
       fVals[3] = 1.0045*p.r9() +0.001; //r9
       fVals[5] = 1.04302*s->etaWidth() - 0.000618; //etawidth
       fVals[6] = 1.00002*s->phiWidth() - 0.000371;  //phiwidth
       fVals[14] = fVals[3]*s->rawEnergy()/b->energy();  //compute consistent e3x3/eseed after r9 rescaling
       if (fVals[15]<=1.0)  // rescale e5x5/eseed only if value is <=1.0, don't allow scaled values to exceed 1.0
         fVals[15] = TMath::Min(1.0,1.0022*p.e5x5()/b->energy());
 
       fVals[4] = fVals[15]*b->energy()/s->rawEnergy(); // compute consistent e5x5()/rawEnergy() after e5x5/eseed resacling
 
       fVals[16] = 0.891832*sqrt(clustertools.localCovariances(*b)[0]) + 0.0009133; //sigietaieta
       fVals[17] = 0.993*sqrt(clustertools.localCovariances(*b)[2]); //sigiphiiphi
 
       fVals[19] = 1.012*bemax/b->energy();                       //crystal energy ratio gap variables
       fVals[20] = 1.0*be2nd/b->energy();
       fVals[21] = 0.94*betop/b->energy();
       fVals[22] = 0.94*bebottom/b->energy();
       fVals[23] = 0.94*beleft/b->energy();
       fVals[24] = 0.94*beright/b->energy();
       fVals[25] = 1.006*be2x5max/b->energy();                       //crystal energy ratio gap variables
       fVals[26] = 1.09*be2x5top/b->energy();
       fVals[27] = 1.09*be2x5bottom/b->energy();
       fVals[28] = 1.09*be2x5left/b->energy();
       fVals[29] = 1.09*be2x5right/b->energy();
 
     }
     else {
       fVals[3] = 1.0086*p.r9() -0.0007;           //r9
       fVals[4] = TMath::Min(1.0,1.0022*p.e5x5()/s->rawEnergy());  //e5x5/rawenergy
       fVals[5] = 0.903254*s->etaWidth() +0.001346;  //etawidth
       fVals[6] = 0.99992*s->phiWidth() +  4.8e-07;   //phiwidth
       fVals[13] = TMath::Min(1.0,1.0022*b->energy()/s->rawEnergy());  //eseed/rawenergy (practically equivalent to e5x5)
 
 
       fVals[14] = fVals[3]*s->rawEnergy()/b->energy(); //compute consistent e3x3/eseed after r9 rescaling
 
       fVals[16] = 0.9947*sqrt(clustertools.localCovariances(*b)[0]) + 0.00003; //sigietaieta
 
       fVals[19] = 1.005*bemax/b->energy();                       //crystal energy ratio gap variables
       fVals[20] = 1.02*be2nd/b->energy();
       fVals[21] = 0.96*betop/b->energy();
       fVals[22] = 0.96*bebottom/b->energy();
       fVals[23] = 0.96*beleft/b->energy();
       fVals[24] = 0.96*beright/b->energy();
       fVals[25] = 1.0075*be2x5max/b->energy();                       //crystal energy ratio gap variables
       fVals[26] = 1.13*be2x5top/b->energy();
       fVals[27] = 1.13*be2x5bottom/b->energy();
       fVals[28] = 1.13*be2x5left/b->energy();
       fVals[29] = 1.13*be2x5right/b->energy();
     }
 
   }
 
   Double_t ecorerr = readervar->GetResponse(fVals)*den;
 
   //printf("ecor = %5f, ecorerr = %5f\n",ecor,ecorerr);
 
   return std::pair<double,double>(ecor,ecorerr);
 }

std::pair< double, double > EGEnergyCorrector::CorrectedEnergyWithErrorV3	(	const reco::GsfElectron &	e,
		const reco::VertexCollection &	vtxcol,
		double	rho,
		EcalClusterLazyTools &	clustertools,
		const edm::EventSetup &	es
	)

Definition at line 645 of file EGEnergyCorrector.cc.

References _ecalLocal, Abs(), b, reco::deltaPhi(), EcalBarrel, reco::CaloCluster::energy(), reco::CaloCluster::eta(), fReadereb, fReaderebvariance, fReaderee, fReadereevariance, fVals, GBRForest::GetResponse(), reco::GsfElectron::hcalOverEcalBc(), reco::CaloCluster::hitsAndFractions(), EcalClusterLocal::localCoordsEB(), reco::CaloCluster::phi(), matplotRender::reader, rho, alignCSCRings::s, mathSSE::sqrt(), and reco::GsfElectron::superCluster().

                                                                                                                                                                                                       {
 
   const SuperClusterRef s = e.superCluster();
   const CaloClusterPtr b = s->seed(); //seed  basic cluster
 
   Bool_t isbarrel =  b->hitsAndFractions().at(0).first.subdetId()==EcalBarrel;
 
   //basic supercluster variables
   fVals[0]  = s->rawEnergy();
   fVals[1]  = s->eta();
   fVals[2]  = s->phi();
   fVals[3]  = clustertools.e3x3(*b)/s->rawEnergy(); //r9
   fVals[4]  = clustertools.e5x5(*b)/s->rawEnergy();
   fVals[5] = s->etaWidth();
   fVals[6] = s->phiWidth();
   fVals[7] = s->clustersSize();
   fVals[8] = e.hcalOverEcalBc();
   fVals[9] = rho;
   fVals[10] = vtxcol.size();
 
   //seed basic cluster variables
   double bemax = clustertools.eMax(*b);
   double be2nd = clustertools.e2nd(*b);
   double betop = clustertools.eTop(*b);
   double bebottom = clustertools.eBottom(*b);
   double beleft = clustertools.eLeft(*b);
   double beright = clustertools.eRight(*b);
 
   double be2x5max = clustertools.e2x5Max(*b);
   double be2x5top = clustertools.e2x5Top(*b);
   double be2x5bottom = clustertools.e2x5Bottom(*b);
   double be2x5left = clustertools.e2x5Left(*b);
   double be2x5right = clustertools.e2x5Right(*b);
 
   fVals[11] = b->eta()-s->eta();
   fVals[12] = reco::deltaPhi(b->phi(),s->phi());
   fVals[13] = b->energy()/s->rawEnergy();
   fVals[14] = clustertools.e3x3(*b)/b->energy();
   fVals[15] = clustertools.e5x5(*b)/b->energy();
   fVals[16] = sqrt(clustertools.localCovariances(*b)[0]); //sigietaieta
   fVals[17] = sqrt(clustertools.localCovariances(*b)[2]); //sigiphiiphi
   fVals[18] = clustertools.localCovariances(*b)[1];       //sigietaiphi
   fVals[19] = bemax/b->energy();                       //crystal energy ratio gap variables
   fVals[20] = be2nd/b->energy();
   fVals[21] = betop/b->energy();
   fVals[22] = bebottom/b->energy();
   fVals[23] = beleft/b->energy();
   fVals[24] = beright/b->energy();
   fVals[25] = be2x5max/b->energy();                       //crystal energy ratio gap variables
   fVals[26] = be2x5top/b->energy();
   fVals[27] = be2x5bottom/b->energy();
   fVals[28] = be2x5left/b->energy();
   fVals[29] = be2x5right/b->energy();
 
   if (isbarrel) {
     //local coordinates and crystal indices (barrel only)
 
     //seed cluster
     float betacry, bphicry, bthetatilt, bphitilt;
     int bieta, biphi;
     _ecalLocal.localCoordsEB(*b,es,betacry,bphicry,bieta,biphi,bthetatilt,bphitilt);
 
     fVals[30] = bieta; //crystal ieta
     fVals[31] = biphi; //crystal iphi
     fVals[32] = bieta%5; //submodule boundary eta symmetry
     fVals[33] = biphi%2; //submodule boundary phi symmetry
     fVals[34] = (TMath::Abs(bieta)<=25)*(bieta%25) + (TMath::Abs(bieta)>25)*((bieta-25*TMath::Abs(bieta)/bieta)%20);  //module boundary eta approximate symmetry
     fVals[35] = biphi%20; //module boundary phi symmetry
     fVals[36] = betacry; //local coordinates with respect to closest crystal center at nominal shower depth
     fVals[37] = bphicry;
 
   }
   else {
     //preshower energy ratio (endcap only)
     fVals[30]  = s->preshowerEnergy()/s->rawEnergy();
   }
 
   Double_t den;
   const GBRForest *reader;
   const GBRForest *readervar;
   if (isbarrel) {
     den = s->rawEnergy();
     reader = fReadereb;
     readervar = fReaderebvariance;
   }
   else {
     den = s->rawEnergy() + s->preshowerEnergy();
     reader = fReaderee;
     readervar = fReadereevariance;
   }
 
   Double_t ecor = reader->GetResponse(fVals)*den;
   Double_t ecorerr = readervar->GetResponse(fVals)*den;
 
   //printf("ecor = %5f, ecorerr = %5f\n",ecor,ecorerr);
 
   return std::pair<double,double>(ecor,ecorerr);
 }

void EGEnergyCorrector::Initialize	(	const edm::EventSetup &	iSetup,
		std::string	regweights,
		bool	weightsFromDB = `false`
	)

Definition at line 45 of file EGEnergyCorrector.cc.

References fIsInitialized, fOwnsForests, fReadereb, fReaderebvariance, fReaderee, fReadereevariance, fVals, edm::EventSetup::get(), and edm::ESHandle< T >::product().

Referenced by EGEnergyAnalyzer::analyze().

                                                                                                         {
     fIsInitialized = kTRUE;
 
     if (fVals) delete [] fVals;
     if (fOwnsForests) {
       if (fReadereb) delete fReadereb;
       if (fReaderebvariance) delete fReaderebvariance;
       if (fReaderee) delete fReaderee;
       if (fReadereevariance) delete fReadereevariance;
     }
 
     fVals = new Float_t[73];
 
     if (weightsFromDB) { //weights from event setup
 
       edm::ESHandle<GBRForest> readereb;
       edm::ESHandle<GBRForest> readerebvar;
       edm::ESHandle<GBRForest> readeree;
       edm::ESHandle<GBRForest> readereevar;
 
       iSetup.get<GBRWrapperRcd>().get(regweights + "_EBCorrection", readereb);
       iSetup.get<GBRWrapperRcd>().get(regweights + "_EBUncertainty", readerebvar);
       iSetup.get<GBRWrapperRcd>().get(regweights + "_EECorrection", readeree);
       iSetup.get<GBRWrapperRcd>().get(regweights + "_EEUncertainty", readereevar);
 
       fReadereb = readereb.product();
       fReaderebvariance = readerebvar.product();
       fReaderee = readeree.product();
       fReadereevariance = readereevar.product();
 
     }
     else { //weights from root file
       fOwnsForests = kTRUE;
 
       TFile *fgbr = TFile::Open(regweights.c_str(),"READ");
       fReadereb = (GBRForest*)fgbr->Get("EBCorrection");
       fReaderebvariance = (GBRForest*)fgbr->Get("EBUncertainty");
       fReaderee = (GBRForest*)fgbr->Get("EECorrection");
       fReadereevariance = (GBRForest*)fgbr->Get("EEUncertainty");
       fgbr->Close();
     }
 
 }