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 *topology)
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
edm::SortedCollection< MapItem >	m_items
std::atomic< std::multimap < CaloTowerDetId, DetId > * >	m_reverseItems
const HcalTopology *	m_topology
bool	standardEB_
bool	standardHB_
bool	standardHE_
bool	standardHF_
bool	standardHO_

Detailed Description

Author

J. Mans - Minnesota

Definition at line 17 of file CaloTowerConstituentsMap.h.

Constructor & Destructor Documentation

CaloTowerConstituentsMap::CaloTowerConstituentsMap ( )

delete

CaloTowerConstituentsMap::~CaloTowerConstituentsMap

(

)

Definition at line 9 of file CaloTowerConstituentsMap.cc.

References m_reverseItems.

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

CaloTowerConstituentsMap::CaloTowerConstituentsMap

(

const HcalTopology *

topology

)

Definition at line 13 of file CaloTowerConstituentsMap.cc.

                                                                                :
  m_topology(topology),
  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 54 of file CaloTowerConstituentsMap.cc.

References edm::hlt::Exception, m_items, and DetId::rawId().

Referenced by 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.push_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 HcalDetId::crystal_ieta_high(), HcalDetId::crystal_ieta_low(), HcalDetId::crystal_iphi_high(), HcalDetId::crystal_iphi_low(), HcalTopology::depthBinInformation(), HcalTopology::firstHERing(), HcalTopology::firstHFRing(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, i, j, HcalTopology::lastHBRing(), HcalTopology::lastHERing(), HcalTopology::lastHFRing(), HcalTopology::lastHORing(), m_items, m_reverseItems, m_topology, EBDetId::MAX_IETA, sd, standardEB_, standardHB_, standardHE_, standardHF_, standardHO_, and HcalDetId::zside().

Referenced by 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 i=m_items.begin(); i!=m_items.end(); i++)
          ptr->insert(std::pair<CaloTowerDetId,DetId>(i->tower,i->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.push_back(j->second);

  // dealing with topo dependency...
  int nd, sd;

  if (standardHB_) {
    if (id.ietaAbs()<=m_topology->lastHBRing()) {
      m_topology->depthBinInformation(HcalBarrel,id.ietaAbs(),nd,sd);
      for (int i=0; i<nd; i++)
    items.push_back(HcalDetId(HcalBarrel,id.ieta(),id.iphi(),i+sd));
    }
  }
  if (standardHO_) {
    if (id.ietaAbs()<=m_topology->lastHORing()) {
      m_topology->depthBinInformation(HcalOuter,id.ietaAbs(),nd,sd);
      for (int i=0; i<nd; i++)
    items.push_back(HcalDetId(HcalOuter,id.ieta(),id.iphi(),i+sd));
    }
  }
  if (standardHE_) {
    if (id.ietaAbs()>=m_topology->firstHERing() && id.ietaAbs()<=m_topology->lastHERing()) {
      m_topology->depthBinInformation(HcalEndcap,id.ietaAbs(),nd,sd);
      for (int i=0; i<nd; i++)
    items.push_back(HcalDetId(HcalEndcap,id.ieta(),id.iphi(),i+sd));
    }
  }
  if (standardHF_) {
    if (id.ietaAbs()>m_topology->firstHFRing() && id.ietaAbs()<=m_topology->lastHFRing()) { 
      int ieta=id.ieta();
      m_topology->depthBinInformation(HcalForward,id.ietaAbs(),nd,sd);
      for (int i=0; i<nd; i++)
    items.push_back(HcalDetId(HcalForward,ieta,id.iphi(),i+sd));
       if (id.ietaAbs() == 30) {
       ieta = 29*id.zside();
       m_topology->depthBinInformation(HcalForward,ieta,nd,sd);
       for (int i=0; i<nd; i++)
         items.push_back(HcalDetId(HcalForward,ieta,id.iphi(),i+sd));
       }
    }
  }
  if (standardEB_ && id.ietaAbs()<=EBDetId::MAX_IETA/5) {
    HcalDetId hid(HcalBarrel,id.ieta(),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.push_back(EBDetId(ie,ip));
  }
  return items;
}

void CaloTowerConstituentsMap::sort

(

)

done adding to the association

Definition at line 61 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 24 of file CaloTowerConstituentsMap.cc.

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

Referenced by CaloTowersCreationAlgo::assignHitEcal(), CaloTowersCreationAlgo::assignHitHcal(), FastL1Region::FillEMCrystals(), 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;

  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_) ) {
    if ((hid.subdet()==HcalForward) && hid.ietaAbs()==29)  // special handling for tower 29
      tid=CaloTowerDetId(30*hid.zside(),hid.iphi());
    else 
      tid=CaloTowerDetId(hid.ieta(),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 156 of file CaloTowerConstituentsMap.cc.

References standardEB_.

                                                     {
  standardEB_=use;
}

void CaloTowerConstituentsMap::useStandardHB

(

bool

use = true

)

add standard (hardcoded) HB items?

Definition at line 144 of file CaloTowerConstituentsMap.cc.

References standardHB_.

                                                     {
  standardHB_=use;
}

void CaloTowerConstituentsMap::useStandardHE

(

bool

use = true

)

add standard (hardcoded) HE items?

Definition at line 147 of file CaloTowerConstituentsMap.cc.

References standardHE_.

                                                     {
  standardHE_=use;
}

void CaloTowerConstituentsMap::useStandardHF

(

bool

use = true

)

add standard (hardcoded) HF items?

Definition at line 153 of file CaloTowerConstituentsMap.cc.

References standardHF_.

                                                     {
  standardHF_=use;
}

void CaloTowerConstituentsMap::useStandardHO

(

bool

use = true

)

add standard (hardcoded) HO items?

Definition at line 150 of file CaloTowerConstituentsMap.cc.

References standardHO_.

                                                     {
  standardHO_=use;
}

Member Data Documentation

edm::SortedCollection<MapItem> CaloTowerConstituentsMap::m_items

private

Definition at line 63 of file CaloTowerConstituentsMap.h.

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

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

mutableprivate

Definition at line 64 of file CaloTowerConstituentsMap.h.

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

const HcalTopology* CaloTowerConstituentsMap::m_topology

private

Definition at line 47 of file CaloTowerConstituentsMap.h.

Referenced by constituentsOf().

bool CaloTowerConstituentsMap::standardEB_

private

Definition at line 53 of file CaloTowerConstituentsMap.h.

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

bool CaloTowerConstituentsMap::standardHB_

private

Definition at line 49 of file CaloTowerConstituentsMap.h.

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

bool CaloTowerConstituentsMap::standardHE_

private

Definition at line 50 of file CaloTowerConstituentsMap.h.

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

bool CaloTowerConstituentsMap::standardHF_

private

Definition at line 51 of file CaloTowerConstituentsMap.h.

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

bool CaloTowerConstituentsMap::standardHO_

private

Definition at line 52 of file CaloTowerConstituentsMap.h.

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