#include <CaloTowerConstituentsMap.h>

Classes
struct	MapItem

Public Member Functions
void	assign (const DetId &cell, const CaloTowerDetId &tower)
	set the association between a DetId and a tower More...

	CaloTowerConstituentsMap ()=delete

	CaloTowerConstituentsMap (const HcalTopology hcaltopo, const CaloTowerTopology cttopo)

std::vector< DetId >	constituentsOf (const CaloTowerDetId &id) const
	Get the constituent detids for this tower id ( not yet implemented ) More...

void	sort ()
	done adding to the association More...

CaloTowerDetId	towerOf (const DetId &id) const
	Get the tower id for this det id (or null if not known) More...

void	useStandardEB (bool use=true)
	add standard (hardcoded) EB items? More...

void	useStandardHB (bool use=true)
	add standard (hardcoded) HB items? More...

void	useStandardHE (bool use=true)
	add standard (hardcoded) HE items? More...

void	useStandardHF (bool use=true)
	add standard (hardcoded) HF items? More...

void	useStandardHO (bool use=true)
	add standard (hardcoded) HO items? More...

	~CaloTowerConstituentsMap ()

Private Attributes
const CaloTowerTopology *	m_cttopo

const HcalTopology *	m_hcaltopo

edm::SortedCollection< MapItem >	m_items

std::atomic< std::multimap< CaloTowerDetId, DetId > * >	m_reverseItems

bool	standardEB_

bool	standardHB_

bool	standardHE_

bool	standardHF_

bool	standardHO_

Detailed Description

Author: J. Mans - Minnesota

Definition at line 18 of file CaloTowerConstituentsMap.h.

Constructor & Destructor Documentation

CaloTowerConstituentsMap::CaloTowerConstituentsMap ( )

delete

CaloTowerConstituentsMap::~CaloTowerConstituentsMap ( )

Definition at line 12 of file CaloTowerConstituentsMap.cc.

References m_reverseItems.

                                                     {
   delete m_reverseItems.load();
   m_reverseItems = nullptr;
 }

CaloTowerConstituentsMap::CaloTowerConstituentsMap	(	const HcalTopology *	hcaltopo,
		const CaloTowerTopology *	cttopo
	)

Definition at line 16 of file CaloTowerConstituentsMap.cc.

     : m_hcaltopo(hcaltopo),
       m_cttopo(cttopo),
       standardHB_(false),
       standardHE_(false),
       standardHF_(false),
       standardHO_(false),
       standardEB_(false),
       m_reverseItems(nullptr) {}

Member Function Documentation

void CaloTowerConstituentsMap::assign	(	const DetId &	cell,
		const CaloTowerDetId &	tower
	)

set the association between a DetId and a tower

Definition at line 53 of file CaloTowerConstituentsMap.cc.

References TauDecayModes::dec, Exception, m_items, and DetId::rawId().

Referenced by CaloTowerConstituentsMapBuilder::assignEEtoHE(), and CaloTowerConstituentsMapBuilder::parseTextMap().

                                                                                     {
   if (m_items.find(cell) != m_items.end()) {
     throw cms::Exception("CaloTowers") << "Cell with id " << std::hex << cell.rawId() << std::dec
                                        << " is already mapped to a CaloTower " << m_items.find(cell)->tower
                                        << std::endl;
   }
   m_items.emplace_back(MapItem(cell, tower));
 }

std::vector< DetId > CaloTowerConstituentsMap::constituentsOf ( const CaloTowerDetId & id ) const

Get the constituent detids for this tower id ( not yet implemented )

copy from the items map

Definition at line 69 of file CaloTowerConstituentsMap.cc.

References CaloTowerTopology::convertCTtoHcal(), gather_cfg::cout, HcalDetId::crystal_ieta_high(), HcalDetId::crystal_ieta_low(), HcalDetId::crystal_iphi_high(), HcalDetId::crystal_iphi_low(), TauDecayModes::dec, HcalTopology::depthBinInformation(), ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, spr::find(), CaloTowerTopology::firstHERing(), CaloTowerTopology::firstHFRing(), HcalTopology::getMergePositionFlag(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, mps_monitormerge::items, dqmiolumiharvest::j, CaloTowerTopology::lastHBRing(), CaloTowerTopology::lastHERing(), CaloTowerTopology::lastHFRing(), CaloTowerTopology::lastHORing(), m_cttopo, m_hcaltopo, m_items, m_reverseItems, EBDetId::MAX_IETA, HcalTopology::mergedDepthDetId(), FastTimerService_cff::range, sd, standardEB_, standardHB_, standardHE_, standardHF_, standardHO_, HcalDetId::zside(), and ecaldqm::zside().

Referenced by spr::eECALmatrix(), CaloTowersCreationAlgo::hadShwrPos(), EgammaHadTower::hasActiveHcal(), CaloTowersCreationAlgo::makeEcalBadChs(), and CaloTowersCreationAlgo::makeHcalDropChMap().

                                                                                         {
 #ifdef EDM_ML_DEBUG
   std::cout << "Get constituent of " << std::hex << id.rawId() << std::dec << " ID " << id << " ieta " << id.ieta()
             << std::endl;
 #endif
   std::vector<DetId> items;
   if (id.ieta() == 0)
     return items;
 
   // build reverse map if needed
   if (!m_reverseItems.load(std::memory_order_acquire)) {
     std::unique_ptr<std::multimap<CaloTowerDetId, DetId>> ptr{new std::multimap<CaloTowerDetId, DetId>};
     for (auto m_item : m_items)
       ptr->insert(std::pair<CaloTowerDetId, DetId>(m_item.tower, m_item.cell));
     std::multimap<CaloTowerDetId, DetId>* expected = nullptr;
     if (m_reverseItems.compare_exchange_strong(expected, ptr.get(), std::memory_order_acq_rel)) {
       ptr.release();
     }
   }
 
   std::multimap<CaloTowerDetId, DetId>::const_iterator j;
   auto range = (*m_reverseItems.load(std::memory_order_acquire)).equal_range(id);
   for (j = range.first; j != range.second; j++)
     items.emplace_back(j->second);
 
   // dealing with topo dependency...
   //use cttopo when dealing with calotower detids
   int nd, sd;
   int hcal_ieta = m_cttopo->convertCTtoHcal(id.ietaAbs());
 
   if (standardHB_) {
     if (id.ietaAbs() <= m_cttopo->lastHBRing()) {
       m_hcaltopo->depthBinInformation(HcalBarrel, hcal_ieta, id.iphi(), id.zside(), nd, sd);
       for (int i = 0; i < nd; i++) {
         if (m_hcaltopo->getMergePositionFlag()) {
           HcalDetId hid =
               m_hcaltopo->mergedDepthDetId(HcalDetId(HcalBarrel, hcal_ieta * id.zside(), id.iphi(), i + sd));
           if (std::find(items.begin(), items.end(), hid) == items.end()) {
             items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
             std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
           }
         } else {
           HcalDetId hid(HcalBarrel, hcal_ieta * id.zside(), id.iphi(), i + sd);
           items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
           std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
         }
       }
     }
   }
   if (standardHO_) {
     if (id.ietaAbs() <= m_cttopo->lastHORing()) {
       m_hcaltopo->depthBinInformation(HcalOuter, hcal_ieta, id.iphi(), id.zside(), nd, sd);
       for (int i = 0; i < nd; i++) {
         HcalDetId hid(HcalOuter, hcal_ieta * id.zside(), id.iphi(), i + sd);
         items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
         std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
       }
     }
   }
   if (standardHE_) {
     if (id.ietaAbs() >= m_cttopo->firstHERing() && id.ietaAbs() <= m_cttopo->lastHERing()) {
       m_hcaltopo->depthBinInformation(HcalEndcap, hcal_ieta, id.iphi(), id.zside(), nd, sd);
       for (int i = 0; i < nd; i++) {
         if (m_hcaltopo->getMergePositionFlag()) {
           HcalDetId hid =
               m_hcaltopo->mergedDepthDetId(HcalDetId(HcalEndcap, hcal_ieta * id.zside(), id.iphi(), i + sd));
           if (std::find(items.begin(), items.end(), hid) == items.end()) {
             items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
             std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
           }
         } else {
           HcalDetId hid(HcalEndcap, hcal_ieta * id.zside(), id.iphi(), i + sd);
           items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
           std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
         }
       }
     }
   }
   if (standardHF_) {
     if (id.ietaAbs() >= m_cttopo->firstHFRing() && id.ietaAbs() <= m_cttopo->lastHFRing()) {
       m_hcaltopo->depthBinInformation(HcalForward, hcal_ieta, id.iphi(), id.zside(), nd, sd);
       for (int i = 0; i < nd; i++) {
         HcalDetId hid(HcalForward, hcal_ieta * id.zside(), id.iphi(), i + sd);
         items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
         std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
       }
       // special handling for first HF tower
       if (id.ietaAbs() == m_cttopo->firstHFRing()) {
         int hcal_ieta2 = hcal_ieta - 1;
         m_hcaltopo->depthBinInformation(HcalForward, hcal_ieta2, id.iphi(), id.zside(), nd, sd);
         for (int i = 0; i < nd; i++) {
           HcalDetId hid(HcalForward, hcal_ieta2 * id.zside(), id.iphi(), i + sd);
           items.emplace_back(hid);
 #ifdef EDM_ML_DEBUG
           std::cout << id << " Depth " << i << ":" << i + sd << " " << hid << "\n";
 #endif
         }
       }
     }
   }
   if (standardEB_ && hcal_ieta <= EBDetId::MAX_IETA / 5) {
     HcalDetId hid(HcalBarrel, hcal_ieta * id.zside(), id.iphi(), 1);  // for the limits
     int etaMin, etaMax;
     if (hid.zside() == -1) {
       etaMin = hid.crystal_ieta_high();
       etaMax = hid.crystal_ieta_low();
     } else {
       etaMin = hid.crystal_ieta_low();
       etaMax = hid.crystal_ieta_high();
     }
     for (int ie = etaMin; ie <= etaMax; ie++)
       for (int ip = hid.crystal_iphi_low(); ip <= hid.crystal_iphi_high(); ip++)
         items.emplace_back(EBDetId(ie, ip));
   }
   return items;
 }

void CaloTowerConstituentsMap::sort ( )

done adding to the association

Definition at line 62 of file CaloTowerConstituentsMap.cc.

References m_items.

                                     {
   m_items.sort();
 
   //  for (auto const & it : m_items)
   //    std::cout << std::hex << it.cell.rawId() << " " << it.tower.rawId() << std::dec << std::endl;
 }

CaloTowerDetId CaloTowerConstituentsMap::towerOf ( const DetId & id ) const

Get the tower id for this det id (or null if not known)

Definition at line 26 of file CaloTowerConstituentsMap.cc.

References CaloTowerTopology::convertHcaltoCT(), DetId::Ecal, EcalBarrel, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, HcalDetId::ietaAbs(), HcalDetId::iphi(), m_cttopo, m_items, DetId::null(), standardEB_, standardHB_, standardHE_, standardHF_, standardHO_, HcalDetId::subdet(), EBDetId::tower_ieta(), EBDetId::tower_iphi(), and HcalDetId::zside().

Referenced by EGHcalRecHitSelector::addDetIds(), CaloTowersCreationAlgo::assignHitEcal(), CaloTowersCreationAlgo::assignHitHcal(), spr::eECALmatrix(), CaloTowersCreationAlgo::makeHcalDropChMap(), PhysicsTowerOrganizer::PhysicsTowerOrganizer(), CaloTowersCreationAlgo::rescale(), and EgammaHadTower::towerOf().

                                                                       {
   CaloTowerDetId tid;  // null to start with
 
   edm::SortedCollection<MapItem>::const_iterator i = m_items.find(id);
   if (i != m_items.end())
     tid = i->tower;
 
   //use hcaltopo when dealing with hcal detids
   if (tid.null()) {
     if (id.det() == DetId::Hcal) {
       HcalDetId hid(id);
       if ((hid.subdet() == HcalBarrel && standardHB_) || (hid.subdet() == HcalEndcap && standardHE_) ||
           (hid.subdet() == HcalOuter && standardHO_) || (hid.subdet() == HcalForward && standardHF_)) {
         tid = CaloTowerDetId(m_cttopo->convertHcaltoCT(hid.ietaAbs(), hid.subdet()) * hid.zside(), hid.iphi());
       }
     } else if (id.det() == DetId::Ecal) {
       EcalSubdetector esd = (EcalSubdetector)id.subdetId();
       if (esd == EcalBarrel && standardEB_) {
         EBDetId ebid(id);
         tid = CaloTowerDetId(ebid.tower_ieta(), ebid.tower_iphi());
       }
     }
   }
 
   return tid;
 }