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;

  static EcalTrigTowerDetId barrelTowerOf(const DetId& id);

  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