PhysicsTowerOrganizer Class Reference

#include <HBHEIsolatedNoiseAlgos.h>

Classes
struct	towercmp

Public Member Functions
void	findNeighbors (const CaloTowerDetId &id, std::set< const PhysicsTower *> &neighbors) const
void	findNeighbors (const PhysicsTower *twr, std::set< const PhysicsTower *> &neighbors) const
void	findNeighbors (int ieta, int iphi, std::set< const PhysicsTower *> &neighbors) const
const PhysicsTower *	findTower (const CaloTowerDetId &id) const
const PhysicsTower *	findTower (int ieta, int iphi) const
	PhysicsTowerOrganizer (const edm::Handle< HBHERecHitCollection > &hbhehitcoll_h, const edm::Handle< EcalRecHitCollection > &ebhitcoll_h, const edm::Handle< EcalRecHitCollection > &eehitcoll_h, const edm::Handle< std::vector< reco::TrackExtrapolation > > &trackextrapcoll_h, const ObjectValidatorAbs &objectvalidator, const CaloTowerConstituentsMap &ctcm, const CaloGeometry &geo)
virtual	~PhysicsTowerOrganizer ()

Private Member Functions
PhysicsTower *	findTower (const CaloTowerDetId &id)
PhysicsTower *	findTower (int ieta, int iphi)
void	insert_ (CaloTowerDetId &id, const HBHERecHit *hit)
void	insert_ (CaloTowerDetId &id, const EcalRecHit *hit)
void	insert_ (CaloTowerDetId &id, const reco::Track *hit)

Private Attributes
std::set< PhysicsTower, towercmp >	towers_

Detailed Description

Definition at line 171 of file HBHEIsolatedNoiseAlgos.h.

Constructor & Destructor Documentation

PhysicsTowerOrganizer()

PhysicsTowerOrganizer::PhysicsTowerOrganizer	(	const edm::Handle< HBHERecHitCollection > &	hbhehitcoll_h,
		const edm::Handle< EcalRecHitCollection > &	ebhitcoll_h,
		const edm::Handle< EcalRecHitCollection > &	eehitcoll_h,
		const edm::Handle< std::vector< reco::TrackExtrapolation > > &	trackextrapcoll_h,
		const ObjectValidatorAbs &	objectvalidator,
		const CaloTowerConstituentsMap &	ctcm,
		const CaloGeometry &	geo
	)

Definition at line 150 of file HBHEIsolatedNoiseAlgos.cc.

References edm::SortedCollection< T, SORT >::begin(), DetId::Ecal, EcalBarrel, EcalEndcap, edm::SortedCollection< T, SORT >::end(), CaloSubdetectorGeometry::getClosestCell(), CaloGeometry::getSubdetectorGeometry(), hit::id, insert_(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, point, reco::TrackExtrapolation::positions(), CaloTowerConstituentsMap::towerOf(), reco::TrackExtrapolation::track(), HLT_2024v14_cff::track, ObjectValidatorAbs::validHit(), and ObjectValidatorAbs::validTrack().

                             {
  // get some geometries
  const CaloSubdetectorGeometry* gEB = geo.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
  const CaloSubdetectorGeometry* gEE = geo.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);

  // do the HCAL hits
  for (HBHERecHitCollection::const_iterator it = hbhehitcoll_h->begin(); it != hbhehitcoll_h->end(); ++it) {
    const HBHERecHit* hit = &(*it);

    // check that the hit is valid
    if (!objectvalidator.validHit(*hit))
      continue;

    // add the hit to the organizer
    CaloTowerDetId tid = ctcm.towerOf(hit->id());
    insert_(tid, hit);
  }

  // do the EB hits
  for (EcalRecHitCollection::const_iterator it = ebhitcoll_h->begin(); it != ebhitcoll_h->end(); ++it) {
    const EcalRecHit* hit = &(*it);

    if (!objectvalidator.validHit(*hit))
      continue;
    CaloTowerDetId tid = ctcm.towerOf(hit->id());
    insert_(tid, hit);
  }

  // do the EE hits
  for (EcalRecHitCollection::const_iterator it = eehitcoll_h->begin(); it != eehitcoll_h->end(); ++it) {
    const EcalRecHit* hit = &(*it);

    if (!objectvalidator.validHit(*hit))
      continue;
    CaloTowerDetId tid = ctcm.towerOf(hit->id());
    insert_(tid, hit);
  }

  // do the tracks
  for (std::vector<reco::TrackExtrapolation>::const_iterator it = trackextrapcoll_h->begin();
       it != trackextrapcoll_h->end();
       ++it) {
    const reco::TrackExtrapolation* extrap = &(*it);
    const reco::Track* track = &(*(extrap->track()));

    // validate track
    if (!objectvalidator.validTrack(*track))
      continue;

    // get the point
    if (extrap->positions().empty())
      continue;
    const GlobalPoint point(
        extrap->positions().front().x(), extrap->positions().front().y(), extrap->positions().front().z());

    if (std::fabs(point.eta()) < 1.479) {
      EBDetId cell = gEB->getClosestCell(point);
      CaloTowerDetId tid = ctcm.towerOf(cell);
      insert_(tid, track);
    } else {
      EEDetId cell = gEE->getClosestCell(point);
      CaloTowerDetId tid = ctcm.towerOf(cell);
      insert_(tid, track);
    }
  }

  return;
}

~PhysicsTowerOrganizer()

virtual PhysicsTowerOrganizer::~PhysicsTowerOrganizer ( )

inlinevirtual

Definition at line 185 of file HBHEIsolatedNoiseAlgos.h.

{}

Member Function Documentation

findNeighbors() [1/3]

void PhysicsTowerOrganizer::findNeighbors	(	const CaloTowerDetId &	id,
		std::set< const PhysicsTower *> &	neighbors
	)		const

Definition at line 275 of file HBHEIsolatedNoiseAlgos.cc.

References findTower(), EcalPhiSymFlatTableProducers_cfi::id, hcalRecHitTable_cff::ieta, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), CaloTowerDetId::ietaAbs(), hcalRecHitTable_cff::iphi, CaloTowerDetId::iphi(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, CaloTowerDetId::zside(), and ecaldqm::zside().

Referenced by findNeighbors(), HBHEHitMapOrganizer::getHPDNeighbors(), and HBHEHitMapOrganizer::HBHEHitMapOrganizer().

                                                                                        {
  // correct for the merging of the |ieta|=28-29 towers
  CaloTowerDetId id(tempid);
  if (tempid.ietaAbs() == 29)
    id = CaloTowerDetId((tempid.ietaAbs() - 1) * tempid.zside(), tempid.iphi());

  std::vector<CaloTowerDetId> ids;
  // get the neighbor with higher iphi
  if (id.ietaAbs() <= 20) {
    if (id.iphi() == 72)
      ids.push_back(CaloTowerDetId(id.ieta(), 1));
    else
      ids.push_back(CaloTowerDetId(id.ieta(), id.iphi() + 1));
  } else {
    if (id.iphi() == 71)
      ids.push_back(CaloTowerDetId(id.ieta(), 1));
    else
      ids.push_back(CaloTowerDetId(id.ieta(), id.iphi() + 2));
  }

  // get the neighbor with the lower iphi
  if (id.ietaAbs() <= 20) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId(id.ieta(), 72));
    else
      ids.push_back(CaloTowerDetId(id.ieta(), id.iphi() - 1));
  } else {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId(id.ieta(), 71));
    else
      ids.push_back(CaloTowerDetId(id.ieta(), id.iphi() - 2));
  }

  // get the neighbor with the higher ietaAbs
  if (id.ietaAbs() == 20 && (id.iphi() % 2) == 0)
    ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), id.iphi() - 1));
  else
    ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), id.iphi()));

  // get the neighbor(s) with the lower ietaAbs
  if (id.ietaAbs() == 21) {
    ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi()));
    ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi() + 1));
  } else if (id.ietaAbs() == 1) {
    ids.push_back(CaloTowerDetId(-id.ieta(), id.iphi()));
  } else {
    ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi()));
  }

  // get the neighbor with higher ieta and higher iphi
  if (id.ietaAbs() <= 19 || (id.ietaAbs() == 20 && (id.iphi() % 2) == 0)) {
    if (id.iphi() == 72)
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), 1));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), id.iphi() + 1));
  } else if (id.ietaAbs() >= 21) {
    if (id.iphi() == 71)
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), 1));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), id.iphi() + 2));
  }

  // get the neighbor with higher ieta and lower iphi
  if (id.ietaAbs() <= 19) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), 72));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), id.iphi() - 1));
  } else if (id.ietaAbs() >= 21 || (id.ietaAbs() == 20 && (id.iphi() % 2) == 1)) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), 71));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() + 1) * id.zside(), id.iphi() - 2));
  }

  // get the neighbor with lower ieta and higher iphi
  if (id.ietaAbs() == 1) {
    if (id.iphi() == 72)
      ids.push_back(CaloTowerDetId(-id.ieta(), 1));
    else
      ids.push_back(CaloTowerDetId(-id.ieta(), id.iphi() + 1));
  } else if (id.ietaAbs() <= 20) {
    if (id.iphi() == 72)
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), 1));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi() + 1));
  } else if (id.ietaAbs() >= 21) {
    if (id.iphi() == 71)
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), 1));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi() + 2));
  }

  // get the neighbor with lower ieta and lower iphi
  if (id.ietaAbs() == 1) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId(-id.ieta(), 72));
    else
      ids.push_back(CaloTowerDetId(-id.ieta(), id.iphi() - 1));
  } else if (id.ietaAbs() <= 20) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), 72));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi() - 1));
  } else if (id.ietaAbs() >= 22) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), 71));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi() - 2));
  } else if (id.ietaAbs() == 21) {
    if (id.iphi() == 1)
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), 72));
    else
      ids.push_back(CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi() - 1));
  }

  // clear neighbors
  neighbors.clear();

  // find the neighbors and add them to the eponymous set
  for (std::vector<CaloTowerDetId>::const_iterator it = ids.begin(); it != ids.end(); ++it) {
    const PhysicsTower* twr = findTower(*it);
    if (twr)
      neighbors.insert(twr);
  }

  return;
}

findNeighbors() [2/3]

void PhysicsTowerOrganizer::findNeighbors	(	const PhysicsTower *	twr,
		std::set< const PhysicsTower *> &	neighbors
	)		const

Definition at line 405 of file HBHEIsolatedNoiseAlgos.cc.

References findNeighbors(), and PhysicsTower::id.

                                                                                                               {
  findNeighbors(twr->id, neighbors);
  return;
}

findNeighbors() [3/3]

void PhysicsTowerOrganizer::findNeighbors	(	int	ieta,
		int	iphi,
		std::set< const PhysicsTower *> &	neighbors
	)		const

Definition at line 410 of file HBHEIsolatedNoiseAlgos.cc.

References findNeighbors(), hcalRecHitTable_cff::ieta, and hcalRecHitTable_cff::iphi.

                                                                                                          {
  findNeighbors(CaloTowerDetId(ieta, iphi), neighbors);
  return;
}

findTower() [1/4]

const PhysicsTower * PhysicsTowerOrganizer::findTower ( const CaloTowerDetId &  id ) const

inline

Definition at line 247 of file HBHEIsolatedNoiseAlgos.cc.

References EcalPhiSymFlatTableProducers_cfi::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), hcalRecHitTable_cff::iphi, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, towers_, and ecaldqm::zside().

Referenced by findNeighbors(), findTower(), HBHEHitMapOrganizer::getHPDNeighbors(), HBHEHitMapOrganizer::HBHEHitMapOrganizer(), and insert_().

                                                                                   {
  // create dummy PhysicsTower
  PhysicsTower dummy;

  // correct for the merging of the |ieta|=28-29 towers
  if (id.ietaAbs() == 29)
    dummy.id = CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi());
  else
    dummy.id = id;

  // search on the dummy
  std::set<PhysicsTower, towercmp>::iterator it = towers_.find(dummy);

  if (it == towers_.end())
    return nullptr;
  return &(*it);
}

findTower() [2/4]

const PhysicsTower * PhysicsTowerOrganizer::findTower ( int  ieta, int  iphi  ) const

inline

Definition at line 265 of file HBHEIsolatedNoiseAlgos.cc.

References findTower(), hcalRecHitTable_cff::ieta, and hcalRecHitTable_cff::iphi.

                                                                             {
  CaloTowerDetId tid(ieta, iphi);
  return findTower(tid);
}

findTower() [3/4]

PhysicsTower * PhysicsTowerOrganizer::findTower ( const CaloTowerDetId &  id )

private

Definition at line 226 of file HBHEIsolatedNoiseAlgos.cc.

References EcalPhiSymFlatTableProducers_cfi::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), hcalRecHitTable_cff::iphi, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, towers_, and ecaldqm::zside().

                                                                       {
  // create dummy PhysicsTower
  PhysicsTower dummy;

  // correct for the merging of the |ieta|=28-29 towers
  if (id.ietaAbs() == 29)
    dummy.id = CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi());
  else
    dummy.id = id;

  // search on the dummy
  std::set<PhysicsTower, towercmp>::iterator it = towers_.find(dummy);

  if (it == towers_.end())
    return nullptr;

  // for whatever reason, I can't get a non-const out of the find method
  PhysicsTower& twr = const_cast<PhysicsTower&>(*it);
  return &twr;
}

findTower() [4/4]

PhysicsTower * PhysicsTowerOrganizer::findTower ( int  ieta, int  iphi  )

private

Definition at line 270 of file HBHEIsolatedNoiseAlgos.cc.

References findTower(), hcalRecHitTable_cff::ieta, and hcalRecHitTable_cff::iphi.

                                                                 {
  CaloTowerDetId tid(ieta, iphi);
  return findTower(tid);
}

insert_() [1/3]

void PhysicsTowerOrganizer::insert_ ( CaloTowerDetId &  id, const HBHERecHit *  hit  )

private

Definition at line 415 of file HBHEIsolatedNoiseAlgos.cc.

References findTower(), PhysicsTower::hcalhits, EcalPhiSymFlatTableProducers_cfi::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), hcalRecHitTable_cff::iphi, towers_, and ecaldqm::zside().

Referenced by PhysicsTowerOrganizer().

                                                                             {
  PhysicsTower* twr = findTower(id);
  if (twr == nullptr) {
    PhysicsTower dummy;
    if (id.ietaAbs() == 29)
      dummy.id = CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi());
    else
      dummy.id = id;
    dummy.hcalhits.insert(hit);
    towers_.insert(dummy);
  } else {
    twr->hcalhits.insert(hit);
  }
  return;
}

insert_() [2/3]

void PhysicsTowerOrganizer::insert_ ( CaloTowerDetId &  id, const EcalRecHit *  hit  )

private

Definition at line 431 of file HBHEIsolatedNoiseAlgos.cc.

References PhysicsTower::ecalhits, findTower(), EcalPhiSymFlatTableProducers_cfi::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), hcalRecHitTable_cff::iphi, towers_, and ecaldqm::zside().

                                                                             {
  PhysicsTower* twr = findTower(id);
  if (twr == nullptr) {
    PhysicsTower dummy;
    if (id.ietaAbs() == 29)
      dummy.id = CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi());
    else
      dummy.id = id;
    dummy.ecalhits.insert(hit);
    towers_.insert(dummy);
  } else {
    twr->ecalhits.insert(hit);
  }
  return;
}

insert_() [3/3]

void PhysicsTowerOrganizer::insert_ ( CaloTowerDetId &  id, const reco::Track *  hit  )

private

Definition at line 447 of file HBHEIsolatedNoiseAlgos.cc.

References findTower(), EcalPhiSymFlatTableProducers_cfi::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), hcalRecHitTable_cff::iphi, towers_, HLT_2024v14_cff::track, PhysicsTower::tracks, and ecaldqm::zside().

                                                                              {
  PhysicsTower* twr = findTower(id);
  if (twr == nullptr) {
    PhysicsTower dummy;
    if (id.ietaAbs() == 29)
      dummy.id = CaloTowerDetId((id.ietaAbs() - 1) * id.zside(), id.iphi());
    else
      dummy.id = id;
    dummy.tracks.insert(track);
    towers_.insert(dummy);
  } else {
    twr->tracks.insert(track);
  }
  return;
}

Member Data Documentation

towers_

std::set<PhysicsTower, towercmp> PhysicsTowerOrganizer::towers_

private

Definition at line 206 of file HBHEIsolatedNoiseAlgos.h.

Referenced by findTower(), and insert_().