EcalTrigTowerConstituentsMap.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOTOPOLOGY_ECALTRIGTOWERCONSTITUENTSMAP_H
#define GEOMETRY_CALOTOPOLOGY_ECALTRIGTOWERCONSTITUENTSMAP_H 1
 
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
 
#include <boost/multi_index_container.hpp>
#include <boost/multi_index/member.hpp>
#include <boost/multi_index/ordered_index.hpp>
#include <boost/tuple/tuple.hpp>
 
#include <vector>
 
class EcalTrigTowerConstituentsMap {
public:
  EcalTrigTowerConstituentsMap();
 
  EcalTrigTowerDetId towerOf(const DetId& id) const;
 
  std::vector<DetId> constituentsOf(const EcalTrigTowerDetId& id) const;
 
  void assign(const DetId& cell, const EcalTrigTowerDetId& tower);
 
private:
  DetId wrapEEDetId(const DetId& id) const;
  DetId wrapEcalTrigTowerDetId(const DetId& id) const;
 
  DetId changeEEDetIdQuadrantAndZ(const DetId& fromid, const int& toQuadrant, const int& tozside) const;
 
  int changeTowerQuadrant(int phiTower, int fromQuadrant, int toQuadrant) const;
 
  struct MapItem {
    MapItem(const DetId& acell, const EcalTrigTowerDetId& atower) : cell(acell), tower(atower) {}
    DetId cell;
    EcalTrigTowerDetId tower;
  };
 
  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;
 
  EcalTowerMap m_items;
};
 
#endif