EcalTools Class Reference

#include <EcalTools.h>

Static Public Member Functions
static bool	deadNeighbour (const DetId &id, const EcalChannelStatus &chs, int chStatusThreshold, int dx, int dy)
	return true if the channel at offsets dx,dy is dead More...
static bool	isHGCalDet (DetId::Detector thedet)
	identify HGCal cells More...
static bool	isNextToBoundary (const DetId &id)
	true if near a crack or ecal border More...
static bool	isNextToDead (const DetId &id, const EcalNextToDeadChannel &es)
	true if the channel is near a dead one (in the 3x3) More...
static bool	isNextToDeadFromNeighbours (const DetId &id, const EcalChannelStatus &chs, int chStatusThreshold)
	same as isNextToDead, but will use information from the neighbour More...
static float	swissCross (const DetId &id, const EcalRecHitCollection &recHits, float recHitThreshold, bool avoidIeta85=true)
	the good old 1-e4/e1. Ignore hits below recHitThreshold More...

Static Private Member Functions
static float	recHitApproxEt (const DetId id, const EcalRecHitCollection &recHits)
static float	recHitE (const DetId id, const EcalRecHitCollection &recHits)
static float	recHitE (const DetId id, const EcalRecHitCollection &recHits, int dEta, int dPhi)

Detailed Description

Definition at line 18 of file EcalTools.h.

Member Function Documentation

deadNeighbour()

bool EcalTools::deadNeighbour ( const DetId &  id, const EcalChannelStatus &  chs, int  chStatusThreshold, int  dx, int  dy  )

static

return true if the channel at offsets dx,dy is dead

Definition at line 87 of file EcalTools.cc.

References PVValHelper::dx, PVValHelper::dy, EcalBarrel, EcalEndcap, EcalCondObjectContainer< T >::end(), EcalCondObjectContainer< T >::find(), EBDetId::offsetBy(), and EEDetId::offsetBy().

Referenced by isNextToDeadFromNeighbours().

                                                                                                                  {
  DetId nid;
  if (id.subdetId() == EcalBarrel)
    nid = EBDetId::offsetBy(id, dx, dy);
  else if (id.subdetId() == EcalEndcap)
    nid = EEDetId::offsetBy(id, dx, dy);

  if (!nid)
    return false;

  EcalChannelStatus::const_iterator chIt = chs.find(nid);
  uint16_t dbStatus = 0;
  if (chIt != chs.end()) {
    // note that
    dbStatus = chIt->getStatusCode();
  } else {
    edm::LogError("EcalDBError") << "No channel status found for xtal " << id.rawId()
                                 << "! something wrong with EcalChannelStatus in your DB? ";
  }

  return (dbStatus >= chStatusThreshold);
}

isHGCalDet()

static bool EcalTools::isHGCalDet ( DetId::Detector  thedet )

inlinestatic

identify HGCal cells

Definition at line 49 of file EcalTools.h.

References DetId::Forward, DetId::Hcal, DetId::HGCalEE, DetId::HGCalHSc, and DetId::HGCalHSi.

Referenced by GsfElectronAlgo::createElectron(), EcalRegressionData::fill(), GEDPhotonProducer::fillPhotonCollection(), ElectronSeedProducer::filterClusters(), and EGRegressionModifierV2::modifyObject().

                                                      {
    return (thedet == DetId::Forward || thedet == DetId::Hcal || thedet == DetId::HGCalEE ||
            thedet == DetId::HGCalHSi || thedet == DetId::HGCalHSc);
  }

isNextToBoundary()

static bool EcalTools::isNextToBoundary ( const DetId &  id )

static

true if near a crack or ecal border

isNextToDead()

bool EcalTools::isNextToDead ( const DetId &  id, const EcalNextToDeadChannel &  es  )

static

true if the channel is near a dead one (in the 3x3)

This function will use the ChannelStatus record to determine if the channel is next to a dead one, using bit 12 of the channel status word

Definition at line 54 of file EcalTools.cc.

References EcalCondObjectContainer< T >::end(), and EcalCondObjectContainer< T >::find().

Referenced by InterestingDetIdFromSuperClusterProducer::produce(), and InterestingDetIdCollectionProducer::produce().

                                                                              {
  EcalNextToDeadChannel::const_iterator chIt = dch.find(id);

  if (chIt != dch.end()) {
    return *chIt;
  } else {
    edm::LogError("EcalDBError") << "No NextToDead  status found for xtal " << id.rawId();
  }

  return false;
}

isNextToDeadFromNeighbours()

bool EcalTools::isNextToDeadFromNeighbours ( const DetId &  id, const EcalChannelStatus &  chs, int  chStatusThreshold  )

static

same as isNextToDead, but will use information from the neighbour

Looks at neighbours in 3x3 and returns true if one of them has chStatus above chStatusThreshold. Use sparingly, slow. Normally isNextToDead() should be used instead.

Definition at line 66 of file EcalTools.cc.

References deadNeighbour().

Referenced by EcalNextToDeadChannelESProducer::setupNextToDeadChannels().

                                                                                                               {
  if (deadNeighbour(id, chs, chStatusThreshold, 1, 0))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, -1, 0))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, 0, 1))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, 0, -1))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, 1, -1))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, 1, 1))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, -1, 1))
    return true;
  if (deadNeighbour(id, chs, chStatusThreshold, -1, -1))
    return true;

  return false;
}

recHitApproxEt()

float EcalTools::recHitApproxEt ( const DetId  id, const EcalRecHitCollection &  recHits  )

staticprivate

Definition at line 134 of file EcalTools.cc.

References EBDetId::approxEta(), EcalBarrel, recHitE(), and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by swissCross().

                                                                                   {
  // for the time being works only for the barrel
  if (id.subdetId() == EcalBarrel) {
    return recHitE(id, recHits) / cosh(EBDetId::approxEta(id));
  }
  return 0;
}

recHitE() [1/2]

float EcalTools::recHitE ( const DetId  id, const EcalRecHitCollection &  recHits  )

staticprivate

Definition at line 123 of file EcalTools.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by recHitApproxEt(), recHitE(), and swissCross().

                                                                            {
  if (id == DetId(0)) {
    return 0;
  } else {
    EcalRecHitCollection::const_iterator it = recHits.find(id);
    if (it != recHits.end())
      return (*it).energy();
  }
  return 0;
}

recHitE() [2/2]

float EcalTools::recHitE ( const DetId  id, const EcalRecHitCollection &  recHits, int  dEta, int  dPhi  )

staticprivate

Definition at line 110 of file EcalTools.cc.

References EcalBarrel, EcalEndcap, EBDetId::offsetBy(), EEDetId::offsetBy(), recHitE(), and FastTrackerRecHitMaskProducer_cfi::recHits.

                                                                                            {
  // in the barrel:   di = dEta   dj = dPhi
  // in the endcap:   di = dX     dj = dY

  DetId nid;
  if (id.subdetId() == EcalBarrel)
    nid = EBDetId::offsetBy(id, di, dj);
  else if (id.subdetId() == EcalEndcap)
    nid = EEDetId::offsetBy(id, di, dj);

  return (nid == DetId(0) ? 0 : recHitE(nid, recHits));
}

swissCross()

float EcalTools::swissCross ( const DetId &  id, const EcalRecHitCollection &  recHits, float  recHitThreshold, bool  avoidIeta85 = true  )

static

the good old 1-e4/e1. Ignore hits below recHitThreshold

Definition at line 9 of file EcalTools.cc.

References funct::abs(), StorageManager_cfg::e1, EcalBarrel, EcalEndcap, EBDetId::ieta(), recHitApproxEt(), recHitE(), FastTrackerRecHitMaskProducer_cfi::recHits, OccupancyTask_cfi::recHitThreshold, and photons_cff::s4.

Referenced by HiSpikeCleaner::produce(), and ecaldqm::RecoSummaryTask::runOnRecHits().

                                              {
  // compute swissCross
  if (id.subdetId() == EcalBarrel) {
    EBDetId ebId(id);
    // avoid recHits at |eta|=85 where one side of the neighbours is missing
    // (may improve considering also eta module borders, but no
    // evidence for the time being that there the performance is
    // different)
    if (abs(ebId.ieta()) == 85 && avoidIeta85)
      return 0;
    // select recHits with Et above recHitThreshold
    if (recHitApproxEt(id, recHits) < recHitThreshold)
      return 0;
    float s4 = 0;
    float e1 = recHitE(id, recHits);
    // protect against nan (if 0 threshold is given above)
    if (e1 == 0)
      return 0;
    s4 += recHitE(id, recHits, 1, 0);
    s4 += recHitE(id, recHits, -1, 0);
    s4 += recHitE(id, recHits, 0, 1);
    s4 += recHitE(id, recHits, 0, -1);
    return 1 - s4 / e1;
  } else if (id.subdetId() == EcalEndcap) {
    EEDetId eeId(id);
    // select recHits with E above recHitThreshold
    float e1 = recHitE(id, recHits);
    if (e1 < recHitThreshold)
      return 0;
    float s4 = 0;
    // protect against nan (if 0 threshold is given above)
    if (e1 == 0)
      return 0;
    s4 += recHitE(id, recHits, 1, 0);
    s4 += recHitE(id, recHits, -1, 0);
    s4 += recHitE(id, recHits, 0, 1);
    s4 += recHitE(id, recHits, 0, -1);
    return 1 - s4 / e1;
  }
  return 0;
}