ElectronRegressionValueMapProducer.cc File Reference

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/Common/interface/View.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/PatCandidates/interface/Electron.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
#include "RecoEgamma/EgammaTools/interface/EcalClusterLocal.h"
#include "TVector2.h"
#include <memory>
#include <vector>

Go to the source code of this file.

Classes
class	ElectronRegressionValueMapProducer

Functions
template<typename LazyTools >
void	calculateValues (EcalClusterLazyToolsBase *tools_tocast, const edm::Ptr< reco::GsfElectron > &iEle, const edm::EventSetup &iSetup, std::vector< float > &vsigmaIEtaIPhi, std::vector< float > &veMax, std::vector< float > &ve2nd, std::vector< float > &veTop, std::vector< float > &veBottom, std::vector< float > &veLeft, std::vector< float > &veRight, std::vector< float > &vclusterMaxDR, std::vector< float > &vclusterMaxDRDPhi, std::vector< float > &vclusterMaxDRDEta, std::vector< float > &vclusterMaxDRRawEnergy, std::vector< float > &vclusterRawEnergy0, std::vector< float > &vclusterRawEnergy1, std::vector< float > &vclusterRawEnergy2, std::vector< float > &vclusterDPhiToSeed0, std::vector< float > &vclusterDPhiToSeed1, std::vector< float > &vclusterDPhiToSeed2, std::vector< float > &vclusterDEtaToSeed0, std::vector< float > &vclusterDEtaToSeed1, std::vector< float > &vclusterDEtaToSeed2, std::vector< int > &veleIPhi, std::vector< int > &veleIEta, std::vector< float > &veleCryPhi, std::vector< float > &veleCryEta)
static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker < edm::ParameterSetDescriptionFiller < ElectronRegressionValueMapProducer > >	s_filler__LINE__ ("ElectronRegressionValueMapProducer")
static const edm::MakerPluginFactory::PMaker < edm::WorkerMaker < ElectronRegressionValueMapProducer > >	s_maker__LINE__ ("ElectronRegressionValueMapProducer")

Function Documentation

template<typename LazyTools >

void calculateValues ( EcalClusterLazyToolsBase *  tools_tocast, const edm::Ptr< reco::GsfElectron > &  iEle, const edm::EventSetup &  iSetup, std::vector< float > &  vsigmaIEtaIPhi, std::vector< float > &  veMax, std::vector< float > &  ve2nd, std::vector< float > &  veTop, std::vector< float > &  veBottom, std::vector< float > &  veLeft, std::vector< float > &  veRight, std::vector< float > &  vclusterMaxDR, std::vector< float > &  vclusterMaxDRDPhi, std::vector< float > &  vclusterMaxDRDEta, std::vector< float > &  vclusterMaxDRRawEnergy, std::vector< float > &  vclusterRawEnergy0, std::vector< float > &  vclusterRawEnergy1, std::vector< float > &  vclusterRawEnergy2, std::vector< float > &  vclusterDPhiToSeed0, std::vector< float > &  vclusterDPhiToSeed1, std::vector< float > &  vclusterDPhiToSeed2, std::vector< float > &  vclusterDEtaToSeed0, std::vector< float > &  vclusterDEtaToSeed1, std::vector< float > &  vclusterDEtaToSeed2, std::vector< int > &  veleIPhi, std::vector< int > &  veleIEta, std::vector< float > &  veleCryPhi, std::vector< float > &  veleCryEta  )

inline

Definition at line 145 of file ElectronRegressionValueMapProducer.cc.

References reco::deltaPhi(), reco::deltaR(), edm::detail::isnan(), EcalClusterLocal::localCoordsEB(), EcalClusterLocal::localCoordsEE(), bookConverter::max, cmsHarvester::sep, and mathSSE::sqrt().

                                                          {
  LazyTools* tools = static_cast<LazyTools*>(tools_tocast);
  
  const auto& the_sc  = iEle->superCluster();
  const auto& theseed = the_sc->seed();
  
  std::vector<float> vCov = tools->localCovariances( *theseed );
  
  const float eMax = tools->eMax( *theseed );
  const float e2nd = tools->e2nd( *theseed );
  const float eTop = tools->eTop( *theseed );
  const float eLeft = tools->eLeft( *theseed );
  const float eRight = tools->eRight( *theseed );
  const float eBottom = tools->eBottom( *theseed );
  
  float dummy;
  int iPhi;
  int iEta;
  float cryPhi;
  float cryEta;
  EcalClusterLocal _ecalLocal;
  if (iEle->isEB()) 
    _ecalLocal.localCoordsEB(*theseed, iSetup, cryEta, cryPhi, iEta, iPhi, dummy, dummy);
  else 
    _ecalLocal.localCoordsEE(*theseed, iSetup, cryEta, cryPhi, iEta, iPhi, dummy, dummy);
  
  double see = (isnan(vCov[0]) ? 0. : sqrt(vCov[0]));
  double spp = (isnan(vCov[2]) ? 0. : sqrt(vCov[2]));
  double sep;    
  if (see*spp > 0)
    sep = vCov[1] / (see * spp);
  else if (vCov[1] > 0)
    sep = 1.0;
  else
    sep = -1.0;
  
  vsigmaIEtaIPhi.push_back(sep);
  veMax.push_back(eMax);
  ve2nd.push_back(e2nd);
  veTop.push_back(eTop);
  veBottom.push_back(eBottom);
  veLeft.push_back(eLeft);
  veRight.push_back(eRight);
  veleIPhi.push_back(iPhi);
  veleIEta.push_back(iEta);
  veleCryPhi.push_back(cryPhi);
  veleCryEta.push_back(cryEta);
  
    // loop over all clusters that aren't the seed
  auto clusend = the_sc->clustersEnd();
  int numberOfClusters =  the_sc->clusters().size();
  
  std::vector<float> _clusterRawEnergy;
  _clusterRawEnergy.resize(std::max(3, numberOfClusters), 0);
  std::vector<float> _clusterDEtaToSeed;
  _clusterDEtaToSeed.resize(std::max(3, numberOfClusters), 0);
  std::vector<float> _clusterDPhiToSeed;
  _clusterDPhiToSeed.resize(std::max(3, numberOfClusters), 0);
  float _clusterMaxDR     = 999.;
  float _clusterMaxDRDPhi = 999.;
  float _clusterMaxDRDEta = 999.;
  float _clusterMaxDRRawEnergy = 0.;
  
  size_t iclus = 0;
  float maxDR = 0;
  edm::Ptr<reco::CaloCluster> pclus;
  for( auto clus = the_sc->clustersBegin(); clus != clusend; ++clus ) {
    pclus = *clus;
    
    if(theseed == pclus ) 
      continue;
    _clusterRawEnergy.push_back(pclus->energy());
    _clusterDPhiToSeed.push_back(reco::deltaPhi(pclus->phi(),theseed->phi()));
    _clusterDEtaToSeed.push_back(pclus->eta() - theseed->eta());
    
    // find cluster with max dR
    if(reco::deltaR(*pclus, *theseed) > maxDR) {
      maxDR = reco::deltaR(*pclus, *theseed);
      _clusterMaxDR = maxDR;
      _clusterMaxDRDPhi = _clusterDPhiToSeed[iclus];
      _clusterMaxDRDEta = _clusterDEtaToSeed[iclus];
      _clusterMaxDRRawEnergy = _clusterRawEnergy[iclus];
    }
    
    ++iclus;
  }
  
  vclusterMaxDR.push_back(_clusterMaxDR);
  vclusterMaxDRDPhi.push_back(_clusterMaxDRDPhi);
  vclusterMaxDRDEta.push_back(_clusterMaxDRDEta);
  vclusterMaxDRRawEnergy.push_back(_clusterMaxDRRawEnergy);
  vclusterRawEnergy0.push_back(_clusterRawEnergy[0]); 
  vclusterRawEnergy1.push_back(_clusterRawEnergy[1]); 
  vclusterRawEnergy2.push_back(_clusterRawEnergy[2]); 
  vclusterDPhiToSeed0.push_back(_clusterDPhiToSeed[0]);
  vclusterDPhiToSeed1.push_back(_clusterDPhiToSeed[1]);
  vclusterDPhiToSeed2.push_back(_clusterDPhiToSeed[2]);
  vclusterDEtaToSeed0.push_back(_clusterDEtaToSeed[0]);
  vclusterDEtaToSeed1.push_back(_clusterDEtaToSeed[1]);
  vclusterDEtaToSeed2.push_back(_clusterDEtaToSeed[2]);
}

static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker<edm::ParameterSetDescriptionFiller< ElectronRegressionValueMapProducer > > s_filler__LINE__ ( "ElectronRegressionValueMapProducer"  )

static

static const edm::MakerPluginFactory ::PMaker< edm::WorkerMaker< ElectronRegressionValueMapProducer > > s_maker__LINE__ ( "ElectronRegressionValueMapProducer"  )

static