CaloTowerConstituentsMap Class Reference

#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 56 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 71 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(), HcalTopology::depthBinInformation(), ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, spr::find(), CaloTowerTopology::firstHERing(), CaloTowerTopology::firstHFRing(), HcalTopology::getMergePositionFlag(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, mps_monitormerge::items, CaloTowerTopology::lastHBRing(), CaloTowerTopology::lastHERing(), CaloTowerTopology::lastHFRing(), CaloTowerTopology::lastHORing(), m_cttopo, m_hcaltopo, m_items, m_reverseItems, EBDetId::MAX_IETA, HcalTopology::mergedDepthDetId(), sd, standardEB_, standardHB_, standardHE_, standardHF_, standardHO_, HcalDetId::zside(), and ecaldqm::zside().

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

                                                                                        {
  std::vector<DetId> 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 63 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 28 of file CaloTowerConstituentsMap.cc.

References CaloTowerTopology::convertHcaltoCT(), DetId::Ecal, EcalBarrel, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, HcalDetId::ietaAbs(), reco::if(), 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 EgammaIsoHcalDetIdCollectionProducer::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;
}

void CaloTowerConstituentsMap::useStandardEB

(

bool

use = true

)

add standard (hardcoded) EB items?

Definition at line 205 of file CaloTowerConstituentsMap.cc.

References standardEB_.

                                                     {
  standardEB_=use;
}

void CaloTowerConstituentsMap::useStandardHB

(

bool

use = true

)

add standard (hardcoded) HB items?

Definition at line 193 of file CaloTowerConstituentsMap.cc.

References standardHB_.

                                                     {
  standardHB_=use;
}

void CaloTowerConstituentsMap::useStandardHE

(

bool

use = true

)

add standard (hardcoded) HE items?

Definition at line 196 of file CaloTowerConstituentsMap.cc.

References standardHE_.

                                                     {
  standardHE_=use;
}

void CaloTowerConstituentsMap::useStandardHF

(

bool

use = true

)

add standard (hardcoded) HF items?

Definition at line 202 of file CaloTowerConstituentsMap.cc.

References standardHF_.

                                                     {
  standardHF_=use;
}

void CaloTowerConstituentsMap::useStandardHO

(

bool

use = true

)

add standard (hardcoded) HO items?

Definition at line 199 of file CaloTowerConstituentsMap.cc.

References standardHO_.

                                                     {
  standardHO_=use;
}

Member Data Documentation

const CaloTowerTopology* CaloTowerConstituentsMap::m_cttopo

private

Definition at line 49 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), and towerOf().

const HcalTopology* CaloTowerConstituentsMap::m_hcaltopo

private

Definition at line 48 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf().

edm::SortedCollection<MapItem> CaloTowerConstituentsMap::m_items

private

Definition at line 65 of file CaloTowerConstituentsMap.h.

Referenced by assign(), constituentsOf(), sort(), and towerOf().

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

mutableprivate

Definition at line 66 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), and ~CaloTowerConstituentsMap().

bool CaloTowerConstituentsMap::standardEB_

private

Definition at line 55 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), towerOf(), and useStandardEB().

bool CaloTowerConstituentsMap::standardHB_

private

Definition at line 51 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), towerOf(), and useStandardHB().

bool CaloTowerConstituentsMap::standardHE_

private

Definition at line 52 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), towerOf(), and useStandardHE().

bool CaloTowerConstituentsMap::standardHF_

private

Definition at line 53 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), towerOf(), and useStandardHF().

bool CaloTowerConstituentsMap::standardHO_

private

Definition at line 54 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf(), towerOf(), and useStandardHO().