EcalElectronicsMapping.h

Go to the documentation of this file.

#ifndef ECALELECTRONICSMAPPING_H
#define ECALELECTRONICSMAPPING_H 1

#include <memory>
#include <iostream>
#include <string>

#include "DataFormats/EcalDetId/interface/EcalDetIdCollections.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.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 <boost/multi_index/mem_fun.hpp>
#include <boost/multi_index/composite_key.hpp>

#include "RecoEcal/EgammaCoreTools/interface/EcalEtaPhiRegion.h"


#include <vector>
#include <map>


class EcalElectronicsMapping {

  
 public:
  
  EcalElectronicsMapping();
  
  EcalElectronicsId getElectronicsId(const DetId& id) const;
  
  EcalTriggerElectronicsId getTriggerElectronicsId(const DetId& id) const;

  DetId getDetId(const EcalElectronicsId& id) const;

  EcalTriggerElectronicsId getTriggerElectronicsId(const EcalElectronicsId& id) const;

  DetId getDetId(const EcalTriggerElectronicsId& id) const;

  EcalElectronicsId getElectronicsId(const EcalTriggerElectronicsId& id) const;

  std::vector<DetId> dccConstituents(int dccId) const;

  std::vector<DetId> dccTowerConstituents(int dccId, int tower) const;

  std::vector<DetId> stripConstituents(int dccId, int tower, int strip) const;

  std::vector<DetId> tccConstituents(int tccId) const;

  std::vector<DetId> ttConstituents(int tccId, int tt) const;

  std::vector<DetId> pseudoStripConstituents(int tccId, int tt, int pseudostrip) const;

  void assign(const DetId& cell, const EcalElectronicsId&, const EcalTriggerElectronicsId& tower);

  std::pair<int, int> getDCCandSC(EcalScDetId id) const;

  std::vector<EcalScDetId> getEcalScDetId(int DCCid, int DCC_Channel,
                      bool ignoreSingleCrystal = true) const;

  int DCCid(const EBDetId& id) const;

  int TCCid(const EBDetId& id) const;

  int iTT(const EcalTrigTowerDetId& id) const;

  int TCCid(const EcalTrigTowerDetId& id) const;

  int DCCid(const EcalTrigTowerDetId& id) const;

  EcalTrigTowerDetId getTrigTowerDetId(int TCCid, int iTT) const;

  EcalSubdetector subdet(int dccid, int mode) const;
  int zside(int dcctcc, int mode) const;

  bool rightTower(int tower) const;

  // methods used for regional unpacking :
  std::vector<int> GetListofFEDs(const EcalEtaPhiRegion region) const ;
  void GetListofFEDs(const EcalEtaPhiRegion region, std::vector<int> & FEDs) const ;
  int GetFED(double eta, double phi) const ;
  int DCCBoundary(int FED) const;

  // methods for retrieving the Laser Monitoring readout number

  int getLMNumber(const DetId& id) const;


 // Geometry of SM in EB :
  static const int kCrystalsInPhi = EBDetId::kCrystalsInPhi ; // per SM
  static const int kTowersInPhi = EBDetId::kTowersInPhi ;  // per SM

 // Geometry of the Trigger Towers :

  static const int kEBTowersInPhi = EcalTrigTowerDetId::kEBTowersInPhi ; // per SM (in the Barrel)
  static const int kEBTowersPerSM = EcalTrigTowerDetId::kEBTowersPerSM ; // per SM (in the Barrel)
  static const int kEBTowersInEta = EcalTrigTowerDetId::kEBTowersInEta; // per SM (in the Barrel)
  static const int kEETowersInEta = EcalTrigTowerDetId::kEETowersInEta; // Endcap
  static const int kEETowersInPhiPerQuadrant = EcalTrigTowerDetId::kEETowersInPhiPerQuadrant; // per Quadrant (EE)
                                                                                                                                                    
  static const int kEETowersInPhiPerTCC = 4;    // each TCC contains 4 towers in phi
  static const int kEETowersInEtaPerInnerTCC = 7;   // each inner TCC contains 7 towers in eta
  static const int kEETowersInEtaPerOuterTCC = 4;   // each outer TCC contains 4 towers in eta
  static const int iEEEtaMinOuter = 18;         // outer TCC : ieta = 18 -> 21
  static const int iEEEtaMinInner = 22;         // inner TCC : ieta = 22 -> 28


 // DCC values :
  static const int MAX_DCCID = EcalElectronicsId::MAX_DCCID ; //To be updated with correct and final number
  static const int MIN_DCCID = EcalElectronicsId::MIN_DCCID;
  static const int MIN_DCCID_EEM = EcalElectronicsId::MIN_DCCID_EEM ;
  static const int MAX_DCCID_EEM = EcalElectronicsId::MAX_DCCID_EEM ;
  static const int MIN_DCCID_EBM = EcalElectronicsId::MIN_DCCID_EBM ;
  static const int MAX_DCCID_EBM = EcalElectronicsId::MAX_DCCID_EBM ;
  static const int MIN_DCCID_EBP = EcalElectronicsId::MIN_DCCID_EBP ;
  static const int MAX_DCCID_EBP = EcalElectronicsId::MAX_DCCID_EBP ;
  static const int MIN_DCCID_EEP = EcalElectronicsId::MIN_DCCID_EEP ;
  static const int MAX_DCCID_EEP = EcalElectronicsId::MAX_DCCID_EEP ;

  static const int DCCID_PHI0_EBM = EcalElectronicsId::DCCID_PHI0_EBM ;   // contains phi = 0 deg.
  static const int DCCID_PHI0_EBP = EcalElectronicsId::DCCID_PHI0_EBP ;

 // TCC values :
  static const int MAX_TCCID = EcalTriggerElectronicsId::MAX_TCCID ; //To be updated with correct and final number
  static const int MIN_TCCID = EcalTriggerElectronicsId::MIN_TCCID ;
  static const int MIN_TCCID_EEM = EcalTriggerElectronicsId::MIN_TCCID_EEM ;
  static const int MAX_TCCID_EEM = EcalTriggerElectronicsId::MAX_TCCID_EEM ;
  static const int MIN_TCCID_EBM = EcalTriggerElectronicsId::MIN_TCCID_EBM ;
  static const int MAX_TCCID_EBM = EcalTriggerElectronicsId::MAX_TCCID_EBM ;
  static const int MIN_TCCID_EBP = EcalTriggerElectronicsId::MIN_TCCID_EBP ;
  static const int MAX_TCCID_EBP = EcalTriggerElectronicsId::MAX_TCCID_EBP ;
  static const int MIN_TCCID_EEP = EcalTriggerElectronicsId::MIN_TCCID_EEP ;
  static const int MAX_TCCID_EEP = EcalTriggerElectronicsId::MAX_TCCID_EEP ;

  static const int TCCID_PHI0_EBM = EcalTriggerElectronicsId::TCCID_PHI0_EBM ;
  static const int TCCID_PHI0_EBP = EcalTriggerElectronicsId::TCCID_PHI0_EBP ;

  static const int TCCID_PHI0_EEM_IN = EcalTriggerElectronicsId::TCCID_PHI0_EEM_IN ;
  static const int TCCID_PHI0_EEM_OUT = EcalTriggerElectronicsId::TCCID_PHI0_EEM_OUT ;
  static const int TCCID_PHI0_EEP_IN = EcalTriggerElectronicsId::TCCID_PHI0_EEP_IN ;
  static const int TCCID_PHI0_EEP_OUT = EcalTriggerElectronicsId::TCCID_PHI0_EEP_OUT ;

  static const int kTCCinPhi = 18;      // Number of TCC "sectors" in phi

 // LaserMonitoring readout numbers :
  static const int MIN_LM_EEM = 73; // corresponds to MIN_DCCID_EEM
  static const int MIN_LM_EBM = 1;  // corresponds to MIN_DCCID_EBM
  static const int MIN_LM_EBP = 37; // corresponds to MIN_DCCID_EBP
  static const int MIN_LM_EEP = 83; // corresponds to MIN_DCCID_EEP
  static const int MAX_LM =     92; // Total number of LaserModules


 private:

  static const int DCCMODE = 0;
  static const int TCCMODE = 1;

  //  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 EcalElectronicsId& aelid, const EcalTriggerElectronicsId& atrelid) : cell(acell),elid(aelid),trelid(atrelid) { }
    DetId cell;
    EcalElectronicsId elid;
    EcalTriggerElectronicsId trelid;
    int dccId() const { return elid.dccId(); }
    int towerId() const {return elid.towerId(); }
    int stripId() const {return elid.stripId(); }
    int tccId() const { return trelid.tccId(); }
    int ttId() const { return trelid.ttId(); }
    int pseudoStripId() const { return trelid.pseudoStripId(); }
  };


  //hashed indexes to be preferred to ordered (faster for lookup, here we are not interested in ordering...)
  typedef boost::multi_index::multi_index_container<
    MapItem,
    boost::multi_index::indexed_by<
    // hashed_unique< member < MapItem,DetId,&MapItem::cell > >,
    // hashed_unique< member < MapItem,EcalElectronicsId,&MapItem::elid > >,
    // hashed_unique< member < MapItem,EcalTriggerElectronicsId,&MapItem::trelid > >
     boost::multi_index::ordered_unique< boost::multi_index::member < MapItem,DetId,&MapItem::cell > >,
     boost::multi_index::ordered_unique< boost::multi_index::member < MapItem,EcalElectronicsId,&MapItem::elid > >,
     boost::multi_index::ordered_unique< boost::multi_index::member < MapItem,EcalTriggerElectronicsId,&MapItem::trelid > >,
     boost::multi_index::ordered_non_unique< boost::multi_index::const_mem_fun < MapItem, int, &MapItem::dccId > >,
     boost::multi_index::ordered_non_unique<
    boost::multi_index::composite_key<
        MapItem,
        boost::multi_index::const_mem_fun < MapItem, int, &MapItem::dccId > ,
        boost::multi_index::const_mem_fun < MapItem, int, &MapItem::towerId  >  
    > >,
     boost::multi_index::ordered_non_unique<
        boost::multi_index::composite_key<
                MapItem,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::dccId > ,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::towerId  >,
        boost::multi_index::const_mem_fun < MapItem, int, &MapItem::stripId >
    > >,
     boost::multi_index::ordered_non_unique< boost::multi_index::const_mem_fun < MapItem, int, &MapItem::tccId > >,
     boost::multi_index::ordered_non_unique<
        boost::multi_index::composite_key<
                MapItem,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::tccId > ,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::ttId  >
        > >,
     boost::multi_index::ordered_non_unique<
        boost::multi_index::composite_key<
                MapItem,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::tccId > ,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::ttId  >,
                boost::multi_index::const_mem_fun < MapItem, int, &MapItem::pseudoStripId >
        > >
    >
    > EcalElectronicsMap;
  
   typedef EcalElectronicsMap::nth_index<0>::type EcalElectronicsMap_by_DetId;
   typedef EcalElectronicsMap::nth_index<1>::type EcalElectronicsMap_by_ElectronicsId;
   typedef EcalElectronicsMap::nth_index<2>::type EcalElectronicsMap_by_TriggerElectronicsId;

   typedef EcalElectronicsMap::nth_index<3>::type EcalElectronicsMap_by_DccId;
   typedef EcalElectronicsMap::nth_index<4>::type EcalElectronicsMap_by_DccId_and_TowerId;
   typedef EcalElectronicsMap::nth_index<5>::type EcalElectronicsMap_by_DccId_TowerId_and_StripId;
  
   typedef EcalElectronicsMap::nth_index<6>::type EcalElectronicsMap_by_TccId;
   typedef EcalElectronicsMap::nth_index<7>::type EcalElectronicsMap_by_TccId_and_TtId;
   typedef EcalElectronicsMap::nth_index<8>::type EcalElectronicsMap_by_TccId_TtId_and_PseudostripId;


  //Needed only in the EE (contains only first quadrant object) 
  EcalElectronicsMap m_items;


  // Maps between DCC and LaserMonitoring readout numbers (take care that EB DCCs and two EE DCCs
  // actually correspond to two LMs. The map contain only the first one).
  // The maps are filled in the constructor of EcalElectronicsMapping.

  std::map<int, int> LaserMonitoringMap_EB;
  std::map<int, int> LaserMonitoringMap_EE;

};

#endif