egamma Namespace Reference

Namespaces
	conv

Classes
struct	ElectronMomentum

Enumerations
enum	HoeMode { HoeMode::SingleTower = 0, HoeMode::TowersBehindCluster = 1 }

Functions
float	classBasedElectronEnergy (reco::GsfElectron const &, reco::BeamSpot const &, EcalClusterFunctionBaseClass const &crackCorrectionFunction)
double	classBasedElectronEnergyUncertainty (reco::GsfElectron const &)
reco::GsfElectron::Classification	classifyElectron (reco::GsfElectron const &)
ElectronMomentum	correctElectronMomentum (reco::GsfElectron const &, TrajectoryStateOnSurface const &)
double	depth1HcalESum (std::vector< CaloTowerDetId > const &towers, CaloTowerCollection const &)
double	depth2HcalESum (std::vector< CaloTowerDetId > const &towers, CaloTowerCollection const &)
double	electronEnergyUncertainty (reco::GsfElectron::Classification c, double eta, double brem, double energy)
std::pair< reco::TrackRef, float >	getClosestCtfToGsf (reco::GsfTrackRef const &, edm::Handle< reco::TrackCollection > const &ctfTracksH, edm::soa::EtaPhiTableView trackEtaPhiTable)
template<typename T >
uint32_t	getRandomSeedFromObj (const edm::Event &iEvent, const T &obj, size_t nrObjs, size_t objNr)
uint32_t	getRandomSeedFromSC (const edm::Event &iEvent, const reco::SuperClusterRef scRef)
bool	hasActiveHcal (std::vector< CaloTowerDetId > const &towers, CaloTowerConstituentsMap const &towerMap, HcalChannelQuality const &hcalQuality, HcalTopology const &hcalTopology)
bool	isBetterElectron (reco::GsfElectron const &, reco::GsfElectron const &)
bool	isInnermostElectron (reco::GsfElectron const &, reco::GsfElectron const &)
int	sharedDets (reco::GsfTrackRef const &, reco::GsfTrackRef const &)
float	sharedEnergy (reco::CaloCluster const &clu1, reco::CaloCluster const &clu2, EcalRecHitCollection const &barrelRecHits, EcalRecHitCollection const &endcapRecHits)
float	sharedEnergy (reco::SuperClusterRef const &sc1, reco::SuperClusterRef const &sc2, EcalRecHitCollection const &barrelRecHits, EcalRecHitCollection const &endcapRecHits)
int	sharedHits (reco::GsfTrackRef const &, reco::GsfTrackRef const &)
double	simpleElectronEnergyUncertainty (reco::GsfElectron const &)
CaloTowerDetId	towerOf (reco::CaloCluster const &cluster, CaloTowerConstituentsMap const &towerMap)
std::vector< CaloTowerDetId >	towersOf (reco::SuperCluster const &sc, CaloTowerConstituentsMap const &towerMap, HoeMode mode=HoeMode::SingleTower)

Enumeration Type Documentation

HoeMode

enum egamma::HoeMode

strong

Enumerator
SingleTower
TowersBehindCluster

Definition at line 16 of file EgammaHadTower.h.

{ SingleTower = 0, TowersBehindCluster = 1 };

Function Documentation

classBasedElectronEnergy()

float egamma::classBasedElectronEnergy	(	reco::GsfElectron const &	electron,
		reco::BeamSpot const &	bs,
		EcalClusterFunctionBaseClass const &	crackCorrectionFunction
	)

Definition at line 295 of file ElectronEnergyCorrector.cc.

                                                                                                    {
  auto elClass = electron.classification();
 
  // new corrections from N. Chanon et al., taken from EcalClusterCorrectionObjectSpecific.cc
  float corr = 1.;
  float corr2 = 1.;
  float energy = electron.superCluster()->energy();
 
  //int subdet = electron.superCluster()->seed()->hitsAndFractions()[0].first.subdetId();
 
  if (electron.isEB()) {
    float cetacorr = fEta(electron.superCluster()->rawEnergy(), electron.superCluster()->eta(), 0) /
                     electron.superCluster()->rawEnergy();
    energy = electron.superCluster()->rawEnergy() * cetacorr;  //previously in CMSSW
    //energy = superCluster.rawEnergy()*fEta(e5x5, superCluster.seed()->eta(), 0)/e5x5;
  } else if (electron.isEE()) {
    energy = electron.superCluster()->rawEnergy() + electron.superCluster()->preshowerEnergy();
  } else {
    edm::LogWarning("ElectronEnergyCorrector::classBasedParameterizationEnergy")
        << "nor barrel neither endcap electron !";
  }
 
  corr = fBremEta(electron.superCluster()->phiWidth() / electron.superCluster()->etaWidth(),
                  electron.superCluster()->eta(),
                  0,
                  elClass);
 
  float et = energy * TMath::Sin(2 * TMath::ATan(TMath::Exp(-electron.superCluster()->eta()))) / corr;
 
  if (electron.isEB()) {
    corr2 = corr * fEt(et, 0, elClass);
  } else if (electron.isEE()) {
    corr2 = corr * fEnergy(energy / corr, 1, elClass);
  } else {
    edm::LogWarning("ElectronEnergyCorrector::classBasedParameterizationEnergy")
        << "nor barrel neither endcap electron !";
  }
 
  float newEnergy = energy / corr2;
 
  // cracks
  double crackcor = 1.;
  for (reco::CaloCluster_iterator cIt = electron.superCluster()->clustersBegin();
       cIt != electron.superCluster()->clustersEnd();
       ++cIt) {
    const reco::CaloClusterPtr cc = *cIt;
    crackcor *= (electron.superCluster()->rawEnergy() + cc->energy() * (crackCorrectionFunction.getValue(*cc) - 1.)) /
                electron.superCluster()->rawEnergy();
  }
  newEnergy *= crackcor;
 
  return newEnergy;

References alignCSCRings::corr, metsig::electron, HCALHighEnergyHPDFilter_cfi::energy, EgHLTOffHistBins_cfi::et, hcaldqm::quantity::fEnergy, and EcalClusterFunctionBaseClass::getValue().

Referenced by GsfElectronAlgo::createElectron().

classBasedElectronEnergyUncertainty()

double egamma::classBasedElectronEnergyUncertainty

(

reco::GsfElectron const &

electron

)

Definition at line 270 of file ElectronEnergyCorrector.cc.

                                                                                  {
  double ecalEnergy = electron.correctedEcalEnergy();
  double eleEta = electron.superCluster()->eta();
  double brem = electron.superCluster()->etaWidth() / electron.superCluster()->phiWidth();
  return egamma::electronEnergyUncertainty(electron.classification(), eleEta, brem, ecalEnergy);

References metsig::electron, and electronEnergyUncertainty().

Referenced by GsfElectronAlgo::createElectron().

classifyElectron()

reco::GsfElectron::Classification egamma::classifyElectron

(

reco::GsfElectron const &

)

Referenced by GsfElectronAlgo::createElectron().

correctElectronMomentum()

egamma::ElectronMomentum egamma::correctElectronMomentum	(	reco::GsfElectron const &	electron,
		TrajectoryStateOnSurface const &	vtxTsos
	)

Definition at line 22 of file ElectronMomentumCorrector.cc.

                                                                                                  {
  int elClass = electron.classification();
 
  //=======================================================================================
  // cluster energy
  //=======================================================================================
 
  float scEnergy = electron.correctedEcalEnergy();
  float errorEnergy = electron.correctedEcalEnergyError();
 
  //=======================================================================================
  // track  momentum
  //=======================================================================================
 
  // basic values
  float trackMomentum = electron.trackMomentumAtVtx().R();
  //float errorTrackMomentum = 999. ;
 
  // tracker momentum scale corrections (Mykhailo Dalchenko)
  double scale = 1.;
  if (electron.isEB()) {
    if (elClass == 0) {
      scale = 1. / (0.00104 * sqrt(trackMomentum) + 1);
    }
    if (elClass == 1) {
      scale = 1. / (0.0017 * sqrt(trackMomentum) + 0.9986);
    }
    if (elClass == 3) {
      scale = 1. / (1.004 - 0.00021 * trackMomentum);
    }
    if (elClass == 4) {
      scale = 0.995;
    }
  } else if (electron.isEE()) {
    if (elClass == 3) {
      scale = 1. / (1.01432 - 0.00201872 * trackMomentum + 0.0000142621 * trackMomentum * trackMomentum);
    }
    if (elClass == 4) {
      scale = 1. / (0.996859 - 0.000345347 * trackMomentum);
    }
  }
  if (scale < 0.)
    scale = 1.;  // CC added protection
  trackMomentum = trackMomentum * scale;
 
  // error (must be done after trackMomentum rescaling)
  MultiGaussianState1D qpState(MultiGaussianStateTransform::multiState1D(vtxTsos, 0));
  GaussianSumUtilities1D qpUtils(qpState);
  float errorTrackMomentum = trackMomentum * trackMomentum * sqrt(qpUtils.mode().variance());
 
  //=======================================================================================
  // combination
  //=======================================================================================
 
  float finalMomentum = electron.p4().t();  // initial
  float finalMomentumError = 999.;
 
  // first check for large errors
 
  if (errorTrackMomentum / trackMomentum > 0.5 && errorEnergy / scEnergy <= 0.5) {
    finalMomentum = scEnergy;
    finalMomentumError = errorEnergy;
  } else if (errorTrackMomentum / trackMomentum <= 0.5 && errorEnergy / scEnergy > 0.5) {
    finalMomentum = trackMomentum;
    finalMomentumError = errorTrackMomentum;
  } else if (errorTrackMomentum / trackMomentum > 0.5 && errorEnergy / scEnergy > 0.5) {
    if (errorTrackMomentum / trackMomentum < errorEnergy / scEnergy) {
      finalMomentum = trackMomentum;
      finalMomentumError = errorTrackMomentum;
    } else {
      finalMomentum = scEnergy;
      finalMomentumError = errorEnergy;
    }
  }
 
  // then apply the combination algorithm
  else {
    // calculate E/p and corresponding error
    float eOverP = scEnergy / trackMomentum;
    float errorEOverP = sqrt((errorEnergy / trackMomentum) * (errorEnergy / trackMomentum) +
                             (scEnergy * errorTrackMomentum / trackMomentum / trackMomentum) *
                                 (scEnergy * errorTrackMomentum / trackMomentum / trackMomentum));
 
    bool eleIsNotInCombination = false;
    if ((eOverP > 1 + 2.5 * errorEOverP) || (eOverP < 1 - 2.5 * errorEOverP) || (eOverP < 0.8) || (eOverP > 1.3)) {
      eleIsNotInCombination = true;
    }
    if (eleIsNotInCombination) {
      if (eOverP > 1) {
        finalMomentum = scEnergy;
        finalMomentumError = errorEnergy;
      } else {
        if (elClass == reco::GsfElectron::GOLDEN) {
          finalMomentum = scEnergy;
          finalMomentumError = errorEnergy;
        }
        if (elClass == reco::GsfElectron::BIGBREM) {
          if (scEnergy < 36) {
            finalMomentum = trackMomentum;
            finalMomentumError = errorTrackMomentum;
          } else {
            finalMomentum = scEnergy;
            finalMomentumError = errorEnergy;
          }
        }
        if (elClass == reco::GsfElectron::BADTRACK) {
          finalMomentum = scEnergy;
          finalMomentumError = errorEnergy;
        }
        if (elClass == reco::GsfElectron::SHOWERING) {
          if (scEnergy < 30) {
            finalMomentum = trackMomentum;
            finalMomentumError = errorTrackMomentum;
          } else {
            finalMomentum = scEnergy;
            finalMomentumError = errorEnergy;
          }
        }
        if (elClass == reco::GsfElectron::GAP) {
          if (scEnergy < 60) {
            finalMomentum = trackMomentum;
            finalMomentumError = errorTrackMomentum;
          } else {
            finalMomentum = scEnergy;
            finalMomentumError = errorEnergy;
          }
        }
      }
    }
 
    else {
      // combination
      finalMomentum = (scEnergy / errorEnergy / errorEnergy + trackMomentum / errorTrackMomentum / errorTrackMomentum) /
                      (1 / errorEnergy / errorEnergy + 1 / errorTrackMomentum / errorTrackMomentum);
      float finalMomentumVariance = 1 / (1 / errorEnergy / errorEnergy + 1 / errorTrackMomentum / errorTrackMomentum);
      finalMomentumError = sqrt(finalMomentumVariance);
    }
  }
 
  //=======================================================================================
  // final set
  //=======================================================================================
 
  auto const& oldMomentum = electron.p4();
 
  return {{oldMomentum.x() * finalMomentum / oldMomentum.t(),
           oldMomentum.y() * finalMomentum / oldMomentum.t(),
           oldMomentum.z() * finalMomentum / oldMomentum.t(),
           finalMomentum},
          errorTrackMomentum,
          finalMomentumError};

References reco::GsfElectron::BADTRACK, reco::GsfElectron::BIGBREM, metsig::electron, EgHLTOffHistBins_cfi::eOverP, reco::GsfElectron::GAP, reco::GsfElectron::GOLDEN, GaussianSumUtilities1D::mode(), MultiGaussianStateTransform::multiState1D(), HLT_FULL_cff::scale, reco::GsfElectron::SHOWERING, mathSSE::sqrt(), reco::btau::trackMomentum, and SingleGaussianState1D::variance().

Referenced by GsfElectronAlgo::createElectron().

depth1HcalESum()

double egamma::depth1HcalESum	(	std::vector< CaloTowerDetId > const &	towers,
		CaloTowerCollection const &	towerCollection
	)

Definition at line 73 of file EgammaHadTower.cc.

                                                                                                                 {
  double esum = 0.;
  for (auto const& tower : towerCollection) {
    if (std::find(towers.begin(), towers.end(), tower.id()) != towers.end()) {
      esum += tower.ietaAbs() < 18 || tower.ietaAbs() > 29 ? tower.hadEnergy() : tower.hadEnergyHeInnerLayer();
    }
  }
  return esum;
}

References spr::find(), hgcalTowerProducer_cfi::tower, HLT_FULL_cff::towerCollection, and HLT_FULL_cff::towers.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth1BehindClusters(), SuperClusterImporter::importToBlock(), and EgammaHLTBcHcalIsolationProducersRegional::produce().

depth2HcalESum()

double egamma::depth2HcalESum	(	std::vector< CaloTowerDetId > const &	towers,
		CaloTowerCollection const &	towerCollection
	)

Definition at line 83 of file EgammaHadTower.cc.

                                                                                                                 {
  double esum = 0.;
  for (auto const& tower : towerCollection) {
    if (std::find(towers.begin(), towers.end(), tower.id()) != towers.end()) {
      esum += tower.hadEnergyHeOuterLayer();
    }
  }
  return esum;
}

References spr::find(), hgcalTowerProducer_cfi::tower, HLT_FULL_cff::towerCollection, and HLT_FULL_cff::towers.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth2BehindClusters(), SuperClusterImporter::importToBlock(), and EgammaHLTBcHcalIsolationProducersRegional::produce().

electronEnergyUncertainty()

double egamma::electronEnergyUncertainty	(	reco::GsfElectron::Classification	c,
		double	eta,
		double	brem,
		double	energy
	)

Definition at line 402 of file EnergyUncertaintyElectronSpecific.cc.

                                                                                                            {
  if (c == GsfElectron::GOLDEN)
    return computeEnergyUncertaintyGolden(eta, brem, energy);
  if (c == GsfElectron::BIGBREM)
    return computeEnergyUncertaintyBigbrem(eta, brem, energy);
  if (c == GsfElectron::SHOWERING)
    return computeEnergyUncertaintyShowering(eta, brem, energy);
  if (c == GsfElectron::BADTRACK)
    return computeEnergyUncertaintyBadTrack(eta, brem, energy);
  if (c == GsfElectron::GAP)
    return computeEnergyUncertaintyCracks(eta, brem, energy);
  throw cms::Exception("GsfElectronAlgo|InternalError") << "unknown classification";
}

References HltBtagPostValidation_cff::c, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, and Exception.

Referenced by classBasedElectronEnergyUncertainty().

getClosestCtfToGsf()

std::pair< TrackRef, float > egamma::getClosestCtfToGsf	(	reco::GsfTrackRef const &	gsfTrackRef,
		edm::Handle< reco::TrackCollection > const &	ctfTracksH,
		edm::soa::EtaPhiTableView	trackEtaPhiTable
	)

Definition at line 16 of file GsfElectronTools.cc.

                                                                                  {
    float maxFracShared = 0;
    TrackRef ctfTrackRef = TrackRef();
    const TrackCollection* ctfTrackCollection = ctfTracksH.product();
 
    float gsfEta = gsfTrackRef->eta();
    float gsfPhi = gsfTrackRef->phi();
    const HitPattern& gsfHitPattern = gsfTrackRef->hitPattern();
 
    constexpr float dR2 = 0.3 * 0.3;
 
    unsigned int counter = 0;
    for (auto ctfTkIter = ctfTrackCollection->begin(); ctfTkIter != ctfTrackCollection->end(); ctfTkIter++, counter++) {
      // dont want to look at every single track in the event!
      using namespace edm::soa::col;
      if (reco::deltaR2(gsfEta, gsfPhi, trackTable.get<Eta>(counter), trackTable.get<Phi>(counter)) > dR2)
        continue;
 
      unsigned int shared = 0;
      int gsfHitCounter = 0;
      int numGsfInnerHits = 0;
      int numCtfInnerHits = 0;
      // get the CTF Track Hit Pattern
      const HitPattern& ctfHitPattern = ctfTkIter->hitPattern();
 
      for (auto elHitsIt = gsfTrackRef->recHitsBegin(); elHitsIt != gsfTrackRef->recHitsEnd();
           elHitsIt++, gsfHitCounter++) {
        if (!((**elHitsIt).isValid()))  //count only valid Hits
        {
          continue;
        }
 
        // look only in the pixels/TIB/TID
        uint32_t gsfHit = gsfHitPattern.getHitPattern(HitPattern::TRACK_HITS, gsfHitCounter);
        if (!(HitPattern::pixelHitFilter(gsfHit) || HitPattern::stripTIBHitFilter(gsfHit) ||
              HitPattern::stripTIDHitFilter(gsfHit))) {
          continue;
        }
 
        numGsfInnerHits++;
 
        int ctfHitsCounter = 0;
        numCtfInnerHits = 0;
        for (auto ctfHitsIt = ctfTkIter->recHitsBegin(); ctfHitsIt != ctfTkIter->recHitsEnd();
             ctfHitsIt++, ctfHitsCounter++) {
          if (!((**ctfHitsIt).isValid()))  //count only valid Hits!
          {
            continue;
          }
 
          uint32_t ctfHit = ctfHitPattern.getHitPattern(HitPattern::TRACK_HITS, ctfHitsCounter);
          if (!(HitPattern::pixelHitFilter(ctfHit) || HitPattern::stripTIBHitFilter(ctfHit) ||
                HitPattern::stripTIDHitFilter(ctfHit))) {
            continue;
          }
 
          numCtfInnerHits++;
 
          if ((**elHitsIt).sharesInput(&(**ctfHitsIt), TrackingRecHit::all)) {
            shared++;
            break;
          }
 
        }  //ctfHits iterator
 
      }  //gsfHits iterator
 
      if ((numGsfInnerHits == 0) || (numCtfInnerHits == 0)) {
        continue;
      }
 
      if (static_cast<float>(shared) / std::min(numGsfInnerHits, numCtfInnerHits) > maxFracShared) {
        maxFracShared = static_cast<float>(shared) / std::min(numGsfInnerHits, numCtfInnerHits);
        ctfTrackRef = TrackRef(ctfTracksH, counter);
      }
 
    }  //ctfTrack iterator
 
    return {ctfTrackRef, maxFracShared};
  }

References TrackingRecHit::all, reco::deltaR2(), reco::HitPattern::getHitPattern(), reco::HitPattern::hitPattern, min(), reco::HitPattern::pixelHitFilter(), edm::Handle< T >::product(), reco::HitPattern::stripTIBHitFilter(), reco::HitPattern::stripTIDHitFilter(), and reco::HitPattern::TRACK_HITS.

Referenced by GsfElectronCoreEcalDrivenProducer::produce(), LowPtGsfElectronCoreProducer::produce(), and GEDGsfElectronCoreProducer::produceElectronCore().

getRandomSeedFromObj()

template<typename T >

uint32_t egamma::getRandomSeedFromObj	(	const edm::Event &	iEvent,
		const T &	obj,
		size_t	nrObjs,
		size_t	objNr
	)

Definition at line 21 of file EgammaRandomSeeds.h.

                                                                                                   {
    std::seed_seq seeder = {int(iEvent.id().event()),
                            int(iEvent.id().luminosityBlock()),
                            int(iEvent.id().run()),
                            int(nrObjs),
                            int(std::numeric_limits<int>::max() * obj.phi() / M_PI) & 0xFFF,
                            int(objNr)};
    uint32_t seed = 0, tries = 10;
    do {
      seeder.generate(&seed, &seed + 1);
      tries++;
    } while (seed == 0 && tries < 10);
    return seed ? seed : iEvent.id().event() + 10000 * objNr;
  }

References iEvent, createfilelist::int, M_PI, SiStripPI::max, getGTfromDQMFile::obj, and fileCollector::seed.

Referenced by CalibratedPhotonProducerT< T >::setSemiDetRandomSeed(), and CalibratedElectronProducerT< T >::setSemiDetRandomSeed().

getRandomSeedFromSC()

uint32_t egamma::getRandomSeedFromSC	(	const edm::Event &	iEvent,
		const reco::SuperClusterRef	scRef
	)

Definition at line 5 of file EgammaRandomSeeds.cc.

                                                                                            {
  const int offset = 0;  //for future expansion
  std::seed_seq seeder = {int(iEvent.id().event()),
                          int(iEvent.id().luminosityBlock()),
                          int(iEvent.id().run()),
                          int(scRef->seed()->seed().rawId()),
                          int(scRef->seed()->hitsAndFractions().size()),
                          offset};
  uint32_t seed = 0, tries = 10;
  do {
    seeder.generate(&seed, &seed + 1);
    tries++;
  } while (seed == 0 && tries < 10);
  return seed ? seed : iEvent.id().event() + 10000 * scRef.key();
}

References iEvent, createfilelist::int, edm::Ref< C, T, F >::key(), hltrates_dqm_sourceclient-live_cfg::offset, and fileCollector::seed.

Referenced by CalibratedPhotonProducerT< T >::setSemiDetRandomSeed(), and CalibratedElectronProducerT< T >::setSemiDetRandomSeed().

hasActiveHcal()

bool egamma::hasActiveHcal	(	std::vector< CaloTowerDetId > const &	towers,
		CaloTowerConstituentsMap const &	towerMap,
		HcalChannelQuality const &	hcalQuality,
		HcalTopology const &	hcalTopology
	)

Definition at line 93 of file EgammaHadTower.cc.

                                                             {
  bool active = false;
  int statusMask = ((1 << HcalChannelStatus::HcalCellOff) | (1 << HcalChannelStatus::HcalCellMask) |
                    (1 << HcalChannelStatus::HcalCellDead));
#ifdef EDM_ML_DEBUG
  std::cout << "DEBUG: hasActiveHcal called with " << towers.size() << " detids. First tower detid ieta "
            << towers.front().ieta() << " iphi " << towers.front().iphi() << std::endl;
#endif
  for (auto towerid : towers) {
    unsigned int ngood = 0, nbad = 0;
    for (DetId id : towerMap.constituentsOf(towerid)) {
      if (id.det() != DetId::Hcal) {
        continue;
      }
      HcalDetId hid(id);
      if (hid.subdet() != HcalBarrel && hid.subdet() != HcalEndcap)
        continue;
#ifdef EDM_ML_DEBUG
      std::cout << "EgammaHadTower DetId " << std::hex << id.rawId() << "  hid.rawId  " << hid.rawId() << std::dec
                << "   sub " << hid.subdet() << "   ieta " << hid.ieta() << "   iphi " << hid.iphi() << "   depth "
                << hid.depth() << std::endl;
#endif
      // Sunanda's fix for 2017 Plan1
      // and removed protection
      int status =
          hcalQuality.getValues((DetId)(hcalTopology.idFront(HcalDetId(id))), /*throwOnFail=*/true)->getValue();
 
#ifdef EDM_ML_DEBUG
      std::cout << "channels status = " << std::hex << status << std::dec << "  int value = " << status << std::endl;
#endif
 
      if (status & statusMask) {
#ifdef EDM_ML_DEBUG
        std::cout << "          BAD!" << std::endl;
#endif
        nbad++;
      } else {
        ngood++;
      }
    }
#ifdef EDM_ML_DEBUG
    std::cout << "    overall ngood " << ngood << " nbad " << nbad << "\n";
#endif
    if (nbad == 0 || (ngood > 0 && nbad < ngood)) {
      active = true;
    }
  }
  return active;
}

References CaloTowerConstituentsMap::constituentsOf(), gather_cfg::cout, TauDecayModes::dec, HcalDetId::depth(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalChannelStatus::HcalCellDead, HcalChannelStatus::HcalCellMask, HcalChannelStatus::HcalCellOff, HcalEndcap, HcalTopology::idFront(), HcalDetId::ieta(), HcalDetId::iphi(), DetId::rawId(), mps_update::status, HcalDetId::subdet(), and HLT_FULL_cff::towers.

Referenced by GEDPhotonProducer::fillPhotonCollection(), and ElectronHcalHelper::hasActiveHcal().

isBetterElectron()

bool egamma::isBetterElectron	(	reco::GsfElectron const &	e1,
		reco::GsfElectron const &	e2
	)

Definition at line 16 of file EgAmbiguityTools.cc.

                                                                              {
    return (std::abs(e1.eSuperClusterOverP() - 1) < std::abs(e2.eSuperClusterOverP() - 1));
  }

References funct::abs(), StorageManager_cfg::e1, and reco::GsfElectron::eSuperClusterOverP().

Referenced by isInnermostElectron(), and GsfElectronProducer::setAmbiguityData().

isInnermostElectron()

bool egamma::isInnermostElectron	(	reco::GsfElectron const &	e1,
		reco::GsfElectron const &	e2
	)

Definition at line 20 of file EgAmbiguityTools.cc.

                                                                                 {
    // retreive first valid hit
    int gsfHitCounter1 = 0;
    for (auto const& elHit : e1.gsfTrack()->recHits()) {
      if (elHit->isValid())
        break;
      gsfHitCounter1++;
    }
 
    int gsfHitCounter2 = 0;
    for (auto const& elHit : e2.gsfTrack()->recHits()) {
      if (elHit->isValid())
        break;
      gsfHitCounter2++;
    }
 
    uint32_t gsfHit1 = e1.gsfTrack()->hitPattern().getHitPattern(HitPattern::TRACK_HITS, gsfHitCounter1);
    uint32_t gsfHit2 = e2.gsfTrack()->hitPattern().getHitPattern(HitPattern::TRACK_HITS, gsfHitCounter2);
 
    if (HitPattern::getSubStructure(gsfHit1) != HitPattern::getSubStructure(gsfHit2)) {
      return (HitPattern::getSubStructure(gsfHit1) < HitPattern::getSubStructure(gsfHit2));
    } else if (HitPattern::getLayer(gsfHit1) != HitPattern::getLayer(gsfHit2)) {
      return (HitPattern::getLayer(gsfHit1) < HitPattern::getLayer(gsfHit2));
    } else {
      return isBetterElectron(e1, e2);
    }
  }

References StorageManager_cfg::e1, reco::GsfElectron::gsfTrack(), and isBetterElectron().

Referenced by GsfElectronProducer::setAmbiguityData().

sharedDets()

int egamma::sharedDets	(	reco::GsfTrackRef const &	gsfTrackRef1,
		reco::GsfTrackRef const &	gsfTrackRef2
	)

Definition at line 95 of file EgAmbiguityTools.cc.

                                                                                   {
    //get the Hit Pattern for the gsfTracks
    const HitPattern& gsfHitPattern1 = gsfTrackRef1->hitPattern();
    const HitPattern& gsfHitPattern2 = gsfTrackRef2->hitPattern();
 
    unsigned int shared = 0;
 
    int gsfHitCounter1 = 0;
    for (trackingRecHit_iterator elHitsIt1 = gsfTrackRef1->recHitsBegin(); elHitsIt1 != gsfTrackRef1->recHitsEnd();
         elHitsIt1++, gsfHitCounter1++) {
      if (!((**elHitsIt1).isValid())) {
        //count only valid Hits
        continue;
      }
      //if (gsfHitCounter1>1) continue; // to test only the first hit of the track 1
      uint32_t gsfHit = gsfHitPattern1.getHitPattern(HitPattern::TRACK_HITS, gsfHitCounter1);
      if (!(HitPattern::pixelHitFilter(gsfHit) || HitPattern::stripTIBHitFilter(gsfHit) ||
            HitPattern::stripTOBHitFilter(gsfHit) || HitPattern::stripTECHitFilter(gsfHit) ||
            HitPattern::stripTIDHitFilter(gsfHit))) {
        continue;
      }
 
      int gsfHitsCounter2 = 0;
      for (trackingRecHit_iterator gsfHitsIt2 = gsfTrackRef2->recHitsBegin(); gsfHitsIt2 != gsfTrackRef2->recHitsEnd();
           gsfHitsIt2++, gsfHitsCounter2++) {
        if (!((**gsfHitsIt2).isValid())) {
          //count only valid Hits!
          continue;
        }
 
        uint32_t gsfHit2 = gsfHitPattern2.getHitPattern(HitPattern::TRACK_HITS, gsfHitsCounter2);
        if (!(HitPattern::pixelHitFilter(gsfHit2) || HitPattern::stripTIBHitFilter(gsfHit2) ||
              HitPattern::stripTOBHitFilter(gsfHit2) || HitPattern::stripTECHitFilter(gsfHit2) ||
              HitPattern::stripTIDHitFilter(gsfHit2)))
          continue;
        if ((**elHitsIt1).geographicalId() == (**gsfHitsIt2).geographicalId())
          shared++;
      }  //gsfHits2 iterator
    }    //gsfHits1 iterator
 
    //std::cout << "[sharedHits] number of shared dets " << shared << std::endl;
    //return shared/min(gsfTrackRef1->numberOfValidHits(),gsfTrackRef2->numberOfValidHits());
    return shared;
  }

References reco::HitPattern::getHitPattern().

sharedEnergy() [1/2]

float egamma::sharedEnergy	(	reco::CaloCluster const &	clu1,
		reco::CaloCluster const &	clu2,
		EcalRecHitCollection const &	barrelRecHits,
		EcalRecHitCollection const &	endcapRecHits
	)

Definition at line 140 of file EgAmbiguityTools.cc.

                                                                {
    double fractionShared = 0;
 
    for (auto const& h1 : clu1.hitsAndFractions()) {
      for (auto const& h2 : clu2.hitsAndFractions()) {
        if (h1.first != h2.first)
          continue;
 
        // here we have common Xtal id
        EcalRecHitCollection::const_iterator itt;
        if (h1.first.subdetId() == EcalBarrel) {
          if ((itt = barrelRecHits.find(h1.first)) != barrelRecHits.end())
            fractionShared += itt->energy();
        } else if (h1.first.subdetId() == EcalEndcap) {
          if ((itt = endcapRecHits.find(h1.first)) != endcapRecHits.end())
            fractionShared += itt->energy();
        }
      }
    }
 
    //std::cout << "[sharedEnergy] shared energy /min(energy1,energy2) " << fractionShared << std::endl;
    return fractionShared;
  }

References HLT_FULL_cff::barrelRecHits, EcalBarrel, EcalEndcap, HLT_FULL_cff::endcapRecHits, HLT_FULL_cff::fractionShared, and reco::CaloCluster::hitsAndFractions().

Referenced by RealisticHitToClusterAssociator::findAndMergeInvisibleClusters(), GsfElectronProducer::setAmbiguityData(), and sharedEnergy().

sharedEnergy() [2/2]

float egamma::sharedEnergy	(	reco::SuperClusterRef const &	sc1,
		reco::SuperClusterRef const &	sc2,
		EcalRecHitCollection const &	barrelRecHits,
		EcalRecHitCollection const &	endcapRecHits
	)

Definition at line 167 of file EgAmbiguityTools.cc.

                                                                {
    double energyShared = 0;
    for (CaloCluster_iterator icl1 = sc1->clustersBegin(); icl1 != sc1->clustersEnd(); icl1++) {
      for (CaloCluster_iterator icl2 = sc2->clustersBegin(); icl2 != sc2->clustersEnd(); icl2++) {
        energyShared += sharedEnergy(**icl1, **icl2, barrelRecHits, endcapRecHits);
      }
    }
    return energyShared;
  }

References HLT_FULL_cff::barrelRecHits, HLT_FULL_cff::endcapRecHits, and sharedEnergy().

sharedHits()

int egamma::sharedHits	(	reco::GsfTrackRef const &	gsfTrackRef1,
		reco::GsfTrackRef const &	gsfTrackRef2
	)

Definition at line 48 of file EgAmbiguityTools.cc.

                                                                                   {
    //get the Hit Pattern for the gsfTracks
    HitPattern const& gsfHitPattern1 = gsfTrackRef1->hitPattern();
    HitPattern const& gsfHitPattern2 = gsfTrackRef2->hitPattern();
 
    unsigned int shared = 0;
 
    int gsfHitCounter1 = 0;
    for (trackingRecHit_iterator elHitsIt1 = gsfTrackRef1->recHitsBegin(); elHitsIt1 != gsfTrackRef1->recHitsEnd();
         elHitsIt1++, gsfHitCounter1++) {
      if (!(*elHitsIt1)->isValid()) {
        //count only valid Hits
        continue;
      }
      //if (gsfHitCounter1>1) continue; // test only the first hit of the track 1
      uint32_t gsfHit = gsfHitPattern1.getHitPattern(HitPattern::TRACK_HITS, gsfHitCounter1);
      if (!(HitPattern::pixelHitFilter(gsfHit) || HitPattern::stripTIBHitFilter(gsfHit) ||
            HitPattern::stripTOBHitFilter(gsfHit) || HitPattern::stripTECHitFilter(gsfHit) ||
            HitPattern::stripTIDHitFilter(gsfHit))) {
        continue;
      }
 
      int gsfHitsCounter2 = 0;
      for (trackingRecHit_iterator gsfHitsIt2 = gsfTrackRef2->recHitsBegin(); gsfHitsIt2 != gsfTrackRef2->recHitsEnd();
           gsfHitsIt2++, gsfHitsCounter2++) {
        if (!(**gsfHitsIt2).isValid()) {
          //count only valid Hits
          continue;
        }
        uint32_t gsfHit2 = gsfHitPattern2.getHitPattern(HitPattern::TRACK_HITS, gsfHitsCounter2);
        if (!(HitPattern::pixelHitFilter(gsfHit2) || HitPattern::stripTIBHitFilter(gsfHit2) ||
              HitPattern::stripTOBHitFilter(gsfHit2) || HitPattern::stripTECHitFilter(gsfHit2) ||
              HitPattern::stripTIDHitFilter(gsfHit2))) {
          continue;
        }
        if ((*elHitsIt1)->sharesInput(&(**gsfHitsIt2), TrackingRecHit::some)) {
          //if (comp.equals(&(**elHitsIt1),&(**gsfHitsIt2))) {
          shared++;
        }
      }  //gsfHits2 iterator
    }    //gsfHits1 iterator
 
    //std::cout << "[sharedHits] number of shared hits " << shared << std::endl;
    return shared;
  }

References reco::HitPattern::getHitPattern(), and TrackingRecHit::some.

simpleElectronEnergyUncertainty()

double egamma::simpleElectronEnergyUncertainty

(

reco::GsfElectron const &

electron

)

Definition at line 277 of file ElectronEnergyCorrector.cc.

                                                                              {
  double error = 999.;
  double ecalEnergy = electron.correctedEcalEnergy();
 
  if (electron.isEB()) {
    float parEB[3] = {5.24e-02, 2.01e-01, 1.00e-02};
    error = ecalEnergy * energyError(ecalEnergy, parEB);
  } else if (electron.isEE()) {
    float parEE[3] = {1.46e-01, 9.21e-01, 1.94e-03};
    error = ecalEnergy * energyError(ecalEnergy, parEE);
  } else {
    edm::LogWarning("ElectronEnergyCorrector::simpleParameterizationUncertainty")
        << "nor barrel neither endcap electron !";
  }
 
  return error;

References metsig::electron, and relativeConstraints::error.

Referenced by GsfElectronAlgo::createElectron().

towerOf()

CaloTowerDetId egamma::towerOf	(	reco::CaloCluster const &	cluster,
		CaloTowerConstituentsMap const &	towerMap
	)

Definition at line 10 of file EgammaHadTower.cc.

                                                                                                       {
  DetId detid = cluster.seed();
  if (detid.det() != DetId::Ecal) {
    // Basic clusters of hybrid super-cluster do not have the seed set; take the first DetId instead
    // Should be checked . The single Tower Mode should be favoured until fixed
    detid = cluster.hitsAndFractions()[0].first;
    if (detid.det() != DetId::Ecal) {
      CaloTowerDetId tower;
      return tower;
    }
  }
  CaloTowerDetId id(towerMap.towerOf(detid));
  return id;
}

References DetId::det(), DetId::Ecal, reco::CaloCluster::hitsAndFractions(), triggerObjects_cff::id, reco::CaloCluster::seed(), hgcalTowerProducer_cfi::tower, and CaloTowerConstituentsMap::towerOf().

Referenced by ecaldqm::TrigPrimTask::runOnDigis(), and towersOf().

towersOf()

std::vector< CaloTowerDetId > egamma::towersOf	(	reco::SuperCluster const &	sc,
		CaloTowerConstituentsMap const &	towerMap,
		HoeMode	mode = HoeMode::SingleTower
	)

Definition at line 25 of file EgammaHadTower.cc.

                                                           {
  constexpr unsigned int nMaxClusters = 4;
 
  std::vector<CaloTowerDetId> towers;
  std::vector<reco::CaloClusterPtr> orderedClusters;
 
  // in this mode, check only the tower behind the seed
  if (mode == HoeMode::SingleTower) {
    towers.push_back(towerOf(*sc.seed(), towerMap));
  }
 
  // in this mode check the towers behind each basic cluster
  if (mode == HoeMode::TowersBehindCluster) {
    // Loop on the basic clusters
    for (auto it = sc.clustersBegin(); it != sc.clustersEnd(); ++it) {
      orderedClusters.push_back(*it);
    }
    std::sort(orderedClusters.begin(), orderedClusters.end(), [](auto& c1, auto& c2) { return (*c1 > *c2); });
    unsigned nclusters = orderedClusters.size();
    for (unsigned iclus = 0; iclus < nclusters && iclus < nMaxClusters; ++iclus) {
      // Get the tower
      CaloTowerDetId id = towerOf(*(orderedClusters[iclus]), towerMap);
#ifdef EDM_ML_DEBUG
      std::cout << "CaloTowerId " << id << std::endl;
#endif
      if (std::find(towers.begin(), towers.end(), id) == towers.end()) {
        towers.push_back(id);
      }
    }
  }
  //  if(towers.size() > 4) {
  //    std::cout << " NTOWERS " << towers.size() << " ";
  //    for(unsigned i=0; i<towers.size() ; ++i) {
  //      std::cout << towers[i] << " ";
  //    }
  //    std::cout <<  std::endl;
  //    for ( unsigned iclus=0 ; iclus < orderedClusters.size(); ++iclus) {
  //      // Get the tower
  //      CaloTowerDetId id = towerOf(*(orderedClusters[iclus]));
  //      std::cout << " Pos " << orderedClusters[iclus]->position() << " " << orderedClusters[iclus]->energy() << " " << id ;
  //    }
  //    std::cout << std::endl;
  //  }
  return towers;
}

References alignmentValidation::c1, reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), gather_cfg::cout, spr::find(), triggerObjects_cff::id, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, reco::SuperCluster::seed(), towerOf(), and HLT_FULL_cff::towers.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hasActiveHcal(), ElectronHcalHelper::hcalTowersBehindClusters(), SuperClusterImporter::importToBlock(), and EgammaHLTBcHcalIsolationProducersRegional::produce().