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