HFRecoEcalCandidateAlgo Class Reference

#include <HFRecoEcalCandidateAlgo.h>

Public Member Functions
	HFRecoEcalCandidateAlgo (bool correct, double e9e25Cut, double intercept2DCut, double intercept2DSlope, const std::vector< double > &e1e9Cut, const std::vector< double > &eCOREe9Cut, const std::vector< double > &eSeLCut, const reco::HFValueStruct hfvv)
void	produce (const edm::Handle< reco::SuperClusterCollection > &SuperClusters, const reco::HFEMClusterShapeAssociationCollection &AssocShapes, reco::RecoEcalCandidateCollection &RecoECand, int nvtx) const

Private Member Functions
reco::RecoEcalCandidate	correctEPosition (const reco::SuperCluster &original, const reco::HFEMClusterShape &shape, int nvtx) const

Private Attributes
const bool	m_correct
const double	m_e1e9Cuthi
const double	m_e1e9Cutlo
const double	m_e9e25Cut
const double	m_eCOREe9Cuthi
const double	m_eCOREe9Cutlo
const int	m_era
const double	m_eSeLCuthi
const double	m_eSeLCutlo
const reco::HFValueStruct	m_hfvv
const double	m_intercept2DCut
const double	m_intercept2DSlope

Detailed Description

Author

K. Klapoetke – Minnesota

Definition at line 23 of file HFRecoEcalCandidateAlgo.h.

Constructor & Destructor Documentation

HFRecoEcalCandidateAlgo()

HFRecoEcalCandidateAlgo::HFRecoEcalCandidateAlgo	(	bool	correct,
		double	e9e25Cut,
		double	intercept2DCut,
		double	intercept2DSlope,
		const std::vector< double > &	e1e9Cut,
		const std::vector< double > &	eCOREe9Cut,
		const std::vector< double > &	eSeLCut,
		const reco::HFValueStruct	hfvv
	)

Definition at line 22 of file HFRecoEcalCandidateAlgo.cc.

    : m_correct(correct),
      m_e9e25Cut(e9e25Cut),
      m_intercept2DCut(intercept2DCut),
      m_intercept2DSlope(intercept2DSlope),
      m_e1e9Cuthi(e1e9Cut[1]),
      m_eCOREe9Cuthi(eCOREe9Cut[1]),
      m_eSeLCuthi(eSeLCut[1]),
      m_e1e9Cutlo(e1e9Cut[0]),
      m_eCOREe9Cutlo(eCOREe9Cut[0]),
      m_eSeLCutlo(eSeLCut[0]),
      m_era(4),
      m_hfvv(hfvv) {}

Member Function Documentation

correctEPosition()

RecoEcalCandidate HFRecoEcalCandidateAlgo::correctEPosition ( const reco::SuperCluster &  original, const reco::HFEMClusterShape &  shape, int  nvtx  ) const

private

Definition at line 43 of file HFRecoEcalCandidateAlgo.cc.

                                                                            {
  double corEta = original.eta();
  //piece-wise log energy correction to eta
  double logel = log(shape.eLong3x3() / 100.0);
  double lowC = 0.00911;
  double hiC = 0.0193;
  double logSlope = 0.0146;
  double logIntercept = -0.00988;
  double sgn = (original.eta() > 0) ? (1.0) : (-1.0);
  if (logel <= 1.25)
    corEta += (lowC)*sgn;
  else if ((logel > 1.25) && (logel <= 2))
    corEta += (logIntercept + logSlope * logel) * sgn;
  else if (logel > 2)
    corEta += hiC * sgn;
  //poly-3 cell eta correction to eta (de)
  double p0 = 0.0125;
  double p1 = -0.0475;
  double p2 = 0.0732;
  double p3 = -0.0506;
  double de = sgn * (p0 + (p1 * (shape.CellEta())) + (p2 * (shape.CellEta() * shape.CellEta())) +
                     (p3 * (shape.CellEta() * shape.CellEta() * shape.CellEta())));
  corEta += de;
 
  double corPhi = original.phi();
  //poly-3 cell phi correction to phi (dp)
  p0 = 0.0115;
  p1 = -0.0394;
  p2 = 0.0486;
  p3 = -0.0335;
  double dp = p0 + (p1 * (shape.CellPhi())) + (p2 * (shape.CellPhi() * shape.CellPhi())) +
              (p3 * (shape.CellPhi() * shape.CellPhi() * shape.CellPhi()));
  corPhi += dp;
 
  double corEnergy = original.energy();
  //energy corrections
  //flat correction for using only long fibers
  double energyCorrect = 0.7397;
  corEnergy = corEnergy / energyCorrect;
  //corection based on ieta for pileup and general
  //find ieta
  int ieta = 0;
  double etabounds[30] = {
      2.846, 2.957, 3.132, 3.307, 3.482, 3.657, 3.833, 4.006, 4.184, 4.357, 4.532, 4.709, 4.882, 5.184};
  for (int kk = 0; kk <= 12; kk++) {
    if ((fabs(corEta) < etabounds[kk + 1]) && (fabs(corEta) > etabounds[kk])) {
      ieta = (corEta > 0) ? (kk + 29) : (-kk - 29);
    }
  }
  //check if ieta is in bounds, should be [-41, -29] U [29, 41]
  if (ieta != 0)
    corEnergy = (m_hfvv.PUSlope(ieta) * 1.0 * (nvtx - 1) + m_hfvv.PUIntercept(ieta) * 1.0) * corEnergy * 1.0 *
                m_hfvv.EnCor(ieta);
 
  //re-calculate the energy vector
  double corPx = corEnergy * cos(corPhi) / cosh(corEta);
  double corPy = corEnergy * sin(corPhi) / cosh(corEta);
  double corPz = corEnergy * tanh(corEta);
 
  //make the candidate
  RecoEcalCandidate corCand(0, math::XYZTLorentzVector(corPx, corPy, corPz, corEnergy), math::XYZPoint(0, 0, 0));
 
  return corCand;
}

References reco::HFEMClusterShape::CellEta(), reco::HFEMClusterShape::CellPhi(), funct::cos(), Calorimetry_cff::dp, reco::HFEMClusterShape::eLong3x3(), reco::HFValueStruct::EnCor(), LEDCalibrationChannels::ieta, GetRecoTauVFromDQM_MC_cff::kk, dqm-mbProfile::log, m_hfvv, definitions::original, p1, p2, p3, reco::HFValueStruct::PUIntercept(), reco::HFValueStruct::PUSlope(), FWPFMaths::sgn(), and funct::sin().

Referenced by produce().

produce()

void HFRecoEcalCandidateAlgo::produce	(	const edm::Handle< reco::SuperClusterCollection > &	SuperClusters,
		const reco::HFEMClusterShapeAssociationCollection &	AssocShapes,
		reco::RecoEcalCandidateCollection &	RecoECand,
		int	nvtx
	)		const

Analyze the hits

Definition at line 110 of file HFRecoEcalCandidateAlgo.cc.

                                                      {
  //get super's and cluster shapes and associations
  for (unsigned int i = 0; i < SuperClusters->size(); ++i) {
    const SuperCluster& supClus = (*SuperClusters)[i];
    reco::SuperClusterRef theClusRef = edm::Ref<SuperClusterCollection>(SuperClusters, i);
    const HFEMClusterShapeRef clusShapeRef = AssocShapes.find(theClusRef)->val;
    const HFEMClusterShape& clusShape = *clusShapeRef;
 
    // basic candidate
    double px = supClus.energy() * cos(supClus.phi()) / cosh(supClus.eta());
    double py = supClus.energy() * sin(supClus.phi()) / cosh(supClus.eta());
    double pz = supClus.energy() * tanh(supClus.eta());
    RecoEcalCandidate theCand(0, math::XYZTLorentzVector(px, py, pz, supClus.energy()), math::XYZPoint(0, 0, 0));
 
    // correct it?
    if (m_correct)
      theCand = correctEPosition(supClus, clusShape, nvtx);
 
    double e9e25 = clusShape.eLong3x3() / clusShape.eLong5x5();
    double e1e9 = clusShape.eLong1x1() / clusShape.eLong3x3();
    double eSeL = hf_egamma::eSeLCorrected(clusShape.eShort3x3(), clusShape.eLong3x3(), 4);
    double var2d = (clusShape.eCOREe9() - (eSeL * m_intercept2DSlope));
 
    bool isAcceptable = true;
    isAcceptable = isAcceptable && (e9e25 > m_e9e25Cut);
    isAcceptable = isAcceptable && (var2d > m_intercept2DCut);
    isAcceptable = isAcceptable && ((e1e9 < m_e1e9Cuthi) && (e1e9 > m_e1e9Cutlo));
    isAcceptable = isAcceptable && ((clusShape.eCOREe9() < m_eCOREe9Cuthi) && (clusShape.eCOREe9() > m_eCOREe9Cutlo));
    isAcceptable = isAcceptable && ((clusShape.eSeL() < m_eSeLCuthi) && (clusShape.eSeL() > m_eSeLCutlo));
 
    if (isAcceptable) {
      theCand.setSuperCluster(theClusRef);
      RecoECand.push_back(theCand);
    }
  }
}

References correctEPosition(), funct::cos(), reco::HFEMClusterShape::eCOREe9(), reco::HFEMClusterShape::eLong1x1(), reco::HFEMClusterShape::eLong3x3(), reco::HFEMClusterShape::eLong5x5(), reco::CaloCluster::energy(), reco::HFEMClusterShape::eSeL(), hf_egamma::eSeLCorrected(), reco::HFEMClusterShape::eShort3x3(), reco::CaloCluster::eta(), edm::AssociationMap< Tag >::find(), mps_fire::i, m_correct, m_e1e9Cuthi, m_e1e9Cutlo, m_e9e25Cut, m_eCOREe9Cuthi, m_eCOREe9Cutlo, m_eSeLCuthi, m_eSeLCutlo, m_intercept2DCut, m_intercept2DSlope, reco::CaloCluster::phi(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, funct::sin(), and edm::helpers::KeyVal< K, V >::val.

Referenced by HFRecoEcalCandidateProducer::produce(), and HLTHFRecoEcalCandidateProducer::produce().

Member Data Documentation

m_correct

const bool HFRecoEcalCandidateAlgo::m_correct

private

Definition at line 45 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_e1e9Cuthi

const double HFRecoEcalCandidateAlgo::m_e1e9Cuthi

private

Definition at line 49 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_e1e9Cutlo

const double HFRecoEcalCandidateAlgo::m_e1e9Cutlo

private

Definition at line 52 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_e9e25Cut

const double HFRecoEcalCandidateAlgo::m_e9e25Cut

private

Definition at line 46 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_eCOREe9Cuthi

const double HFRecoEcalCandidateAlgo::m_eCOREe9Cuthi

private

Definition at line 50 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_eCOREe9Cutlo

const double HFRecoEcalCandidateAlgo::m_eCOREe9Cutlo

private

Definition at line 53 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_era

const int HFRecoEcalCandidateAlgo::m_era

private

Definition at line 55 of file HFRecoEcalCandidateAlgo.h.

m_eSeLCuthi

const double HFRecoEcalCandidateAlgo::m_eSeLCuthi

private

Definition at line 51 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_eSeLCutlo

const double HFRecoEcalCandidateAlgo::m_eSeLCutlo

private

Definition at line 54 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_hfvv

const reco::HFValueStruct HFRecoEcalCandidateAlgo::m_hfvv

private

Definition at line 56 of file HFRecoEcalCandidateAlgo.h.

Referenced by correctEPosition().

m_intercept2DCut

const double HFRecoEcalCandidateAlgo::m_intercept2DCut

private

Definition at line 47 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().

m_intercept2DSlope

const double HFRecoEcalCandidateAlgo::m_intercept2DSlope

private

Definition at line 48 of file HFRecoEcalCandidateAlgo.h.

Referenced by produce().