HBHEHitMapOrganizer Class Reference

#include <HBHEIsolatedNoiseAlgos.h>

Public Member Functions
void	getDiHits (std::vector< HBHEHitMap > &v, double energy) const
void	getHPDs (std::vector< HBHEHitMap > &v, double energy) const
void	getMonoHits (std::vector< HBHEHitMap > &v, double energy) const
void	getRBXs (std::vector< HBHEHitMap > &v, double energy) const
	HBHEHitMapOrganizer (const edm::Handle< HBHERecHitCollection > &hbhehitcoll_h, const ObjectValidatorAbs &objvalidator, const PhysicsTowerOrganizer &pto, const HcalFrontEndMap *hfemap)
virtual	~HBHEHitMapOrganizer ()

Private Member Functions
void	getHPDNeighbors (const HBHERecHit *hit, std::vector< const HBHERecHit *> &neighbors, const PhysicsTowerOrganizer &pto)

Private Attributes
std::vector< HBHEHitMap >	dihits_
const HcalFrontEndMap *	hfemap_
std::map< int, HBHEHitMap >	hpds_
std::vector< HBHEHitMap >	monohits_
std::map< int, HBHEHitMap >	rbxs_

Detailed Description

Definition at line 330 of file HBHEIsolatedNoiseAlgos.h.

Constructor & Destructor Documentation

HBHEHitMapOrganizer()

HBHEHitMapOrganizer::HBHEHitMapOrganizer	(	const edm::Handle< HBHERecHitCollection > &	hbhehitcoll_h,
		const ObjectValidatorAbs &	objvalidator,
		const PhysicsTowerOrganizer &	pto,
		const HcalFrontEndMap *	hfemap
	)

Definition at line 743 of file HBHEIsolatedNoiseAlgos.cc.

References edm::SortedCollection< T, SORT >::begin(), dihits_, edm::SortedCollection< T, SORT >::end(), PhysicsTowerOrganizer::findNeighbors(), PhysicsTowerOrganizer::findTower(), getHPDNeighbors(), hfemap_, hpds_, hit::id, hcalRecHitTable_cff::ieta, HBHEHitMap::insert(), hcalRecHitTable_cff::iphi, HcalFrontEndMap::lookupRBXIndex(), HcalFrontEndMap::lookupRMIndex(), monohits_, rbxs_, l1tHGCalTowerProducer_cfi::tower, and ObjectValidatorAbs::validHit().

    : hfemap_(hfemap) {
  // loop over the hits
  for (HBHERecHitCollection::const_iterator it = hbhehitcoll_h->begin(); it != hbhehitcoll_h->end(); ++it) {
    const HBHERecHit* hit = &(*it);
    if (!objvalidator.validHit(*hit))
      continue;

    // get the Physics Tower and the neighbors
    const PhysicsTower* tower = pto.findTower(hit->id().ieta(), hit->id().iphi());

    std::set<const PhysicsTower*> neighbors;
    pto.findNeighbors(hit->id().ieta(), hit->id().iphi(), neighbors);

    // organize the RBXs
    int rbxidnum = hfemap_->lookupRBXIndex(hit->id());
    rbxs_[rbxidnum].insert(hit, tower, neighbors);

    // organize the HPDs
    int hpdidnum = hfemap_->lookupRMIndex(hit->id());
    hpds_[hpdidnum].insert(hit, tower, neighbors);

    // organize the dihits
    std::vector<const HBHERecHit*> hpdneighbors;
    getHPDNeighbors(hit, hpdneighbors, pto);

    if (hpdneighbors.size() == 1) {
      std::vector<const HBHERecHit*> hpdneighborsneighbors;
      getHPDNeighbors(hpdneighbors[0], hpdneighborsneighbors, pto);

      if (hpdneighborsneighbors.size() == 1 && hpdneighborsneighbors[0] == hit &&
          hit->energy() > hpdneighbors[0]->energy()) {
        // we've found two hits who are neighbors in the same HPD, but who have no other
        // neighbors (in the same HPD) in common.  In order not to double-count, we
        // require that the first hit has more energy

        const PhysicsTower* tower2 = pto.findTower(hpdneighbors[0]->id().ieta(), hpdneighbors[0]->id().iphi());
        std::set<const PhysicsTower*> neighbors2;
        pto.findNeighbors(hpdneighbors[0]->id().ieta(), hpdneighbors[0]->id().iphi(), neighbors2);

        HBHEHitMap dihit;
        dihit.insert(hit, tower, neighbors);
        dihit.insert(hpdneighbors[0], tower2, neighbors2);
        dihits_.push_back(dihit);
      }
    } else if (hpdneighbors.empty()) {
      // organize the monohits
      HBHEHitMap monohit;
      monohit.insert(hit, tower, neighbors);
      monohits_.push_back(monohit);
    }

  }  // finished looping over HBHERecHits
  return;
}

~HBHEHitMapOrganizer()

virtual HBHEHitMapOrganizer::~HBHEHitMapOrganizer ( )

inlinevirtual

Definition at line 337 of file HBHEIsolatedNoiseAlgos.h.

{}

Member Function Documentation

getDiHits()

void HBHEHitMapOrganizer::getDiHits	(	std::vector< HBHEHitMap > &	v,
		double	energy
	)		const

Definition at line 820 of file HBHEIsolatedNoiseAlgos.cc.

References dihits_, hcalRecHitTable_cff::energy, and findQualityFiles::v.

                                                                                 {
  for (std::vector<HBHEHitMap>::const_iterator it = dihits_.begin(); it != dihits_.end(); ++it) {
    if (it->hitEnergy() > energy)
      v.push_back(*it);
  }
  return;
}

getHPDNeighbors()

void HBHEHitMapOrganizer::getHPDNeighbors ( const HBHERecHit *  hit, std::vector< const HBHERecHit *> &  neighbors, const PhysicsTowerOrganizer &  pto  )

private

Definition at line 836 of file HBHEIsolatedNoiseAlgos.cc.

References PhysicsTowerOrganizer::findNeighbors(), PhysicsTowerOrganizer::findTower(), hfemap_, HBHERecHit::id(), hit::id, HcalFrontEndMap::lookupRMIndex(), and groupFilesInBlocks::temp.

Referenced by HBHEHitMapOrganizer().

                                                                            {
  std::set<const PhysicsTower*> temp;
  pto.findNeighbors(hit->id().ieta(), hit->id().iphi(), temp);

  // make sure to include the same tower that the hit is in
  temp.insert(pto.findTower(hit->id().ieta(), hit->id().iphi()));

  // loop over the rechits in the temp neighbors
  for (std::set<const PhysicsTower*>::const_iterator it1 = temp.begin(); it1 != temp.end(); ++it1) {
    for (std::set<const HBHERecHit*>::const_iterator it2 = (*it1)->hcalhits.begin(); it2 != (*it1)->hcalhits.end();
         ++it2) {
      const HBHERecHit* hit2(*it2);
      if (hit != hit2 && hfemap_->lookupRMIndex(hit->id()) == hfemap_->lookupRMIndex(hit2->id())) {
        neighbors.push_back(hit2);
      }
    }
  }
  return;
}

getHPDs()

void HBHEHitMapOrganizer::getHPDs	(	std::vector< HBHEHitMap > &	v,
		double	energy
	)		const

Definition at line 811 of file HBHEIsolatedNoiseAlgos.cc.

References hcalRecHitTable_cff::energy, hpds_, genParticles_cff::map, and findQualityFiles::v.

                                                                               {
  for (std::map<int, HBHEHitMap>::const_iterator it = hpds_.begin(); it != hpds_.end(); ++it) {
    const HBHEHitMap& map = it->second;
    if (map.hitEnergy() > energy)
      v.push_back(map);
  }
  return;
}

getMonoHits()

void HBHEHitMapOrganizer::getMonoHits	(	std::vector< HBHEHitMap > &	v,
		double	energy
	)		const

Definition at line 828 of file HBHEIsolatedNoiseAlgos.cc.

References hcalRecHitTable_cff::energy, monohits_, and findQualityFiles::v.

                                                                                   {
  for (std::vector<HBHEHitMap>::const_iterator it = monohits_.begin(); it != monohits_.end(); ++it) {
    if (it->hitEnergy() > energy)
      v.push_back(*it);
  }
  return;
}

getRBXs()

void HBHEHitMapOrganizer::getRBXs	(	std::vector< HBHEHitMap > &	v,
		double	energy
	)		const

Definition at line 802 of file HBHEIsolatedNoiseAlgos.cc.

References hcalRecHitTable_cff::energy, genParticles_cff::map, rbxs_, and findQualityFiles::v.

                                                                               {
  for (std::map<int, HBHEHitMap>::const_iterator it = rbxs_.begin(); it != rbxs_.end(); ++it) {
    const HBHEHitMap& map = it->second;
    if (map.hitEnergy() > energy)
      v.push_back(map);
  }
  return;
}

Member Data Documentation

dihits_

std::vector<HBHEHitMap> HBHEHitMapOrganizer::dihits_

private

Definition at line 347 of file HBHEIsolatedNoiseAlgos.h.

Referenced by getDiHits(), and HBHEHitMapOrganizer().

hfemap_

const HcalFrontEndMap* HBHEHitMapOrganizer::hfemap_

private

Definition at line 345 of file HBHEIsolatedNoiseAlgos.h.

Referenced by getHPDNeighbors(), and HBHEHitMapOrganizer().

hpds_

std::map<int, HBHEHitMap> HBHEHitMapOrganizer::hpds_

private

Definition at line 346 of file HBHEIsolatedNoiseAlgos.h.

Referenced by getHPDs(), and HBHEHitMapOrganizer().

monohits_

std::vector<HBHEHitMap> HBHEHitMapOrganizer::monohits_

private

Definition at line 347 of file HBHEIsolatedNoiseAlgos.h.

Referenced by getMonoHits(), and HBHEHitMapOrganizer().

rbxs_

std::map<int, HBHEHitMap> HBHEHitMapOrganizer::rbxs_

private

Definition at line 346 of file HBHEIsolatedNoiseAlgos.h.

Referenced by getRBXs(), and HBHEHitMapOrganizer().