EcalTrigTowerConstituentsMap Class Reference

#include <EcalTrigTowerConstituentsMap.h>

Classes
struct	MapItem

Public Member Functions
void	assign (const DetId &cell, const EcalTrigTowerDetId &tower)
	set the association between a DetId and a tower More...
std::vector< DetId >	constituentsOf (const EcalTrigTowerDetId &id) const
	Get the constituent detids for this tower id. More...
	EcalTrigTowerConstituentsMap ()
EcalTrigTowerDetId	towerOf (const DetId &id) const
	Get the tower id for this det id (or null if not known) More...

Private Types
typedef boost::multi_index_container< MapItem, boost::multi_index::indexed_by< boost::multi_index::ordered_unique< boost::multi_index::member< MapItem, DetId,&MapItem::cell > >, boost::multi_index::ordered_non_unique< boost::multi_index::member< MapItem, EcalTrigTowerDetId,&MapItem::tower > > > >	EcalTowerMap
typedef EcalTowerMap::nth_index< 0 >::type	EcalTowerMap_by_DetId
typedef EcalTowerMap::nth_index< 1 >::type	EcalTowerMap_by_towerDetId

Private Member Functions
DetId	changeEEDetIdQuadrantAndZ (const DetId &fromid, const int &toQuadrant, const int &tozside) const
int	changeTowerQuadrant (int phiTower, int fromQuadrant, int toQuadrant) const
DetId	wrapEcalTrigTowerDetId (const DetId &id) const
	Wrap a generic EcalTrigTowerDetId to the equivalent one in z+ Quadrant 1 (from 0 < phi < pi/2) More...
DetId	wrapEEDetId (const DetId &id) const
	Wrap a generic EEDetId to the equivalent one in z+ Quadrant 1 (from 0 < phi < pi/2) More...

Private Attributes
EcalTowerMap	m_items

Detailed Description

Author

P.Meridiani

Definition at line 19 of file EcalTrigTowerConstituentsMap.h.

Member Typedef Documentation

typedef boost::multi_index_container< MapItem, boost::multi_index::indexed_by< boost::multi_index::ordered_unique< boost::multi_index::member < MapItem,DetId,&MapItem::cell > >, boost::multi_index::ordered_non_unique< boost::multi_index::member < MapItem,EcalTrigTowerDetId,&MapItem::tower> > > > EcalTrigTowerConstituentsMap::EcalTowerMap

private

Definition at line 57 of file EcalTrigTowerConstituentsMap.h.

typedef EcalTowerMap::nth_index<0>::type EcalTrigTowerConstituentsMap::EcalTowerMap_by_DetId

private

Definition at line 59 of file EcalTrigTowerConstituentsMap.h.

typedef EcalTowerMap::nth_index<1>::type EcalTrigTowerConstituentsMap::EcalTowerMap_by_towerDetId

private

Definition at line 60 of file EcalTrigTowerConstituentsMap.h.

Constructor & Destructor Documentation

EcalTrigTowerConstituentsMap::EcalTrigTowerConstituentsMap

(

)

Definition at line 9 of file EcalTrigTowerConstituentsMap.cc.

{
}

Member Function Documentation

void EcalTrigTowerConstituentsMap::assign	(	const DetId &	cell,
		const EcalTrigTowerDetId &	tower
	)

set the association between a DetId and a tower

Definition at line 159 of file EcalTrigTowerConstituentsMap.cc.

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

Referenced by EcalTrigTowerConstituentsMapBuilder::parseTextMap().

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

DetId EcalTrigTowerConstituentsMap::changeEEDetIdQuadrantAndZ ( const DetId &  fromid, const int &  toQuadrant, const int &  tozside  ) const

private

Definition at line 102 of file EcalTrigTowerConstituentsMap.cc.

References DetId::det(), DetId::Ecal, EcalEndcap, EEDetId::iquadrant(), EEDetId::ix(), EEDetId::iy(), DetId::subdetId(), and EEDetId::XYMODE.

{
  if (! (fromid.det() == DetId::Ecal && fromid.subdetId() == EcalEndcap) )
    return EEDetId(0);
  
  EEDetId myId(fromid);
  int dQuadrant = toQuadrant-myId.iquadrant(); 
  switch(dQuadrant%4){
  case 0:
    return DetId(EEDetId(myId.ix(),myId.iy(),tozside,EEDetId::XYMODE));
    break;
  case 1:
    /* adjacent tower: they are symetric*/
    return DetId(EEDetId(101-myId.ix(),myId.iy(),tozside,EEDetId::XYMODE));
    break;
  case 2:
    /* opposite quadrant: they are identical*/
    return DetId(EEDetId(101-myId.ix(),101-myId.iy(),tozside,EEDetId::XYMODE));   
    break;
  case 3:
    /* adjacent tower: they are symetric*/
    return DetId(EEDetId(myId.ix(),101-myId.iy(),tozside,EEDetId::XYMODE));   
    break;
  default:
    /*should never be reached*/
    edm::LogError("EcalTrigTowerConstituentsMapError") << "This should never be reached. Profound error!";    
  }
  return EEDetId(0);
}

int EcalTrigTowerConstituentsMap::changeTowerQuadrant ( int  phiTower, int  fromQuadrant, int  toQuadrant  ) const

private

Definition at line 132 of file EcalTrigTowerConstituentsMap.cc.

References EcalTrigTowerDetId::kEETowersInPhiPerQuadrant.

                                                                                                         {
  int newPhiTower = phiTower;
  int dQuadrant = toQuadrant-fromQuadrant; 
  
  switch(dQuadrant%4){
  case 0:
    newPhiTower = phiTower;
    break;
  case 1:
    /* adjacent tower: they are symetric*/
    newPhiTower = EcalTrigTowerDetId::kEETowersInPhiPerQuadrant * 2 - phiTower + 1;
    break;
  case 2:
    /* opposite quadrant: they are identical*/
    newPhiTower = phiTower + EcalTrigTowerDetId::kEETowersInPhiPerQuadrant * 2;
    break;
  case 3:
    /* adjacent tower: they are symetric*/
    newPhiTower = EcalTrigTowerDetId::kEETowersInPhiPerQuadrant * 4 - phiTower + 1;
    break;
  default:
    /*should never be reached*/
    edm::LogError("EcalTrigTowerConstituentsMapError") << "This should never be reached. Profound error!";    
  }
  return newPhiTower;
}

std::vector< DetId > EcalTrigTowerConstituentsMap::constituentsOf

(

const EcalTrigTowerDetId &

id

)

const

Get the constituent detids for this tower id.

Definition at line 167 of file EcalTrigTowerConstituentsMap.cc.

References MillePedeFileConverter_cfg::e, DetId::Ecal, EcalBarrel, EcalEndcap, EcalTriggerTower, EBDetId::ETAPHIMODE, mps_monitormerge::items, m_items, AlCaHLTBitMon_ParallelJobs::p, and HcalDigiParam_cfi::zside.

Referenced by EcalRecHitsValidation::analyze(), EcalDetIdToBeRecoveredProducer::produce(), ecaldqm::TimingClient::producePlots(), EcalRecHitWorkerRecover::run(), ecaldqm::TrigPrimTask::runOnEmulTPs(), ecaldqm::TrigPrimTask::runOnRealTPs(), EcalDeadCellTriggerPrimitiveFilter::setEvtRecHitstatus(), and ecaldqm::DQWorkerClient::towerAverage_().

                                                                                                {
  std::vector<DetId> items;
  
  if ( id.det() == DetId::Ecal && id.subdetId() == EcalTriggerTower && id.subDet() == EcalBarrel ) 
    {
      //--------------------
      // Ecal Barrel
      //--------------------
      // trigger towers are 5x5 crystals in the barrel
      int etaxtalMin=(id.ietaAbs()-1)*5+1;
      int phixtalMin=((id.iphi()-1)*5+11)%360;
      if(phixtalMin<=0)phixtalMin+=360;
      int etaxtalMax=id.ietaAbs()*5;
      int phixtalMax=((id.iphi())*5+10)%360;
      if(phixtalMax<=0) phixtalMax+=360;
      for(int e=etaxtalMin;e<=etaxtalMax;e++) 
    for(int p=phixtalMin;p<=phixtalMax;p++) 
      items.emplace_back(DetId(EBDetId(id.zside()*e,p,EBDetId::ETAPHIMODE)));
    }
  else if (id.det() == DetId::Ecal && id.subdetId() == EcalTriggerTower && id.subDet() == EcalEndcap) 
    {
      //--------------------
      // Ecal Endcap
      //--------------------
      //DetId myId=wrapEcalTrigTowerDetId(id);  
      EcalTowerMap_by_towerDetId::const_iterator lb,ub;
      //boost::tuples::tie(lb,ub)=get<1>(m_items).equal_range(myId);
      boost::tuples::tie(lb,ub)=boost::get<1>(m_items).equal_range(id);
      while (lb!=ub)
    {
      //EEDetId mappedId((*lb).cell);
      //items.emplace_back(changeEEDetIdQuadrantAndZ(mappedId,id.iquadrant(),id.zside()));
      items.emplace_back((*lb).cell);
      ++lb;
    }
    }
  
  return items;
}

EcalTrigTowerDetId EcalTrigTowerConstituentsMap::towerOf

(

const DetId &

id

)

const

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

Definition at line 13 of file EcalTrigTowerConstituentsMap.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, mps_fire::i, m_items, EBDetId::tower(), and EEDetId::zside().

Referenced by EcalRecHitWorkerRecover::estimateEnergy(), EcalTBReadout::findTTlist(), EcalDeadCellTriggerPrimitiveFilter::getChannelStatusMaps(), EcalDeadCellDeltaRFilter::getChannelStatusMaps(), EcalSelectiveReadout::getCrystalInterest(), PFRecoTauDiscriminationAgainstElectronDeadECAL::loopXtals(), EcalTPSkimmer::produce(), EcalDetIdToBeRecoveredProducer::produce(), ecaldqm::OccupancyClient::producePlots(), EcalTBReadout::readOut(), EcalReadoutTools::readOutUnitOf(), PFEcalBarrelRecHitCreator::readOutUnitOf(), EcalSelectiveReadoutValidation::readOutUnitOf(), EcalRecHitWorkerRecover::run(), EcalSelectiveReadout::runSelectiveReadout0(), EcalSelectiveReadoutValidation::setTtEtSums(), and EcalSelectiveReadoutSuppressor::setTtFlags().

                                                                              {
  if (id.det() == DetId::Ecal && id.subdetId() == EcalBarrel)
    {
      //--------------------
      // Ecal Barrel
      //--------------------
      EBDetId myId(id);
      return myId.tower();
    }
  else if (id.det() == DetId::Ecal && id.subdetId() == EcalEndcap)
    {
      //--------------------
      // Ecal Endcap
      //--------------------
      EEDetId originalId(id);
      // DetId wrappedId=wrapEEDetId(id);
      DetId wrappedId(originalId);
      EcalTowerMap::const_iterator i=m_items.find(wrappedId);
      if (i!=m_items.end()) 
    {
      int etaTower=i->tower.ietaAbs();
      int phiTower=i->tower.iphi();
      //trigger tower <-> crystal maping read 
      //..........from file and done only for 1 quadrant
      //move from quadrant 1 to the actual one:
      // phiTower = changeTowerQuadrant(phiTower, 1, originalId.iquadrant());
      // std::cout << originalId.zside() <<  " " << etaTower << " " << phiTower << std::endl;
      return EcalTrigTowerDetId(originalId.zside(),EcalEndcap,etaTower,phiTower);
    }
    }
  return EcalTrigTowerDetId(0);
}

DetId EcalTrigTowerConstituentsMap::wrapEcalTrigTowerDetId ( const DetId &  id ) const

private

Wrap a generic EcalTrigTowerDetId to the equivalent one in z+ Quadrant 1 (from 0 < phi < pi/2)

Definition at line 74 of file EcalTrigTowerConstituentsMap.cc.

References DetId::det(), DetId::Ecal, EcalEndcap, EcalTriggerTower, EcalTrigTowerDetId::ietaAbs(), EcalTrigTowerDetId::iphi(), EcalTrigTowerDetId::iquadrant(), EcalTrigTowerDetId::kEETowersInPhiPerQuadrant, EcalTrigTowerDetId::subDet(), and DetId::subdetId().

{
  EcalTrigTowerDetId etid(id);

  if (! (etid.det() == DetId::Ecal && etid.subdetId() == EcalTriggerTower && etid.subDet() == EcalEndcap) )
    return EcalTrigTowerDetId(0);
  
  switch((etid.iquadrant()-1)%4)
    {
    case 0:
      return DetId(EcalTrigTowerDetId(1,EcalEndcap,etid.ietaAbs(),etid.iphi()));
      break;
    case 1:
      return DetId(EcalTrigTowerDetId(1,EcalEndcap,etid.ietaAbs(),EcalTrigTowerDetId::kEETowersInPhiPerQuadrant * 2 - etid.iphi() + 1));
      break;
    case 2:
      return DetId(EcalTrigTowerDetId(1,EcalEndcap,etid.ietaAbs(),etid.iphi() - EcalTrigTowerDetId::kEETowersInPhiPerQuadrant * 2 ));
      break;
    case 3:
      return DetId(EcalTrigTowerDetId(1,EcalEndcap,etid.ietaAbs(),EcalTrigTowerDetId::kEETowersInPhiPerQuadrant * 4 - etid.iphi() + 1));
      break;
    default:
      /*should never be reached*/
      edm::LogError("EcalTrigTowerConstituentsMapError") << "This should never be reached. Profound error!";    
    }
  return EcalTrigTowerDetId(0);
}

DetId EcalTrigTowerConstituentsMap::wrapEEDetId ( const DetId &  id ) const

private

Wrap a generic EEDetId to the equivalent one in z+ Quadrant 1 (from 0 < phi < pi/2)

Definition at line 46 of file EcalTrigTowerConstituentsMap.cc.

References DetId::det(), DetId::Ecal, EcalEndcap, EEDetId::iquadrant(), EEDetId::ix(), EEDetId::iy(), DetId::subdetId(), and EEDetId::XYMODE.

{
  if (! (eeid.det() == DetId::Ecal && eeid.subdetId() == EcalEndcap) )
    return EEDetId(0);
  
  EEDetId myId(eeid);
  switch((myId.iquadrant()-1)%4)
    {
    case 0:
      return DetId(EEDetId(myId.ix(),myId.iy(),1,EEDetId::XYMODE));
      break;
    case 1:
      return DetId(EEDetId(101-myId.ix(),myId.iy(),1,EEDetId::XYMODE));
      break;
    case 2:
      return DetId(EEDetId(101-myId.ix(),101-myId.iy(),1,EEDetId::XYMODE));   
      break;
    case 3:
      return DetId(EEDetId(myId.ix(),101-myId.iy(),1,EEDetId::XYMODE));   
      break;
    default:
      /*should never be reached*/
      edm::LogError("EcalTrigTowerConstituentsMapError") << "This should never be reached. Profound error!";    
    }
  return EEDetId(0);
}

Member Data Documentation

EcalTowerMap EcalTrigTowerConstituentsMap::m_items

private

Definition at line 62 of file EcalTrigTowerConstituentsMap.h.

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