EcalRegionCabling.h

Go to the documentation of this file.

#ifndef EcalRegionCabling_H
#define EcalRegionCabling_H

#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/EcalMapping/interface/ESElectronicsMapper.h"

#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

class EcalRegionCabling {
 public:
  EcalRegionCabling(edm::ParameterSet & conf, const EcalElectronicsMapping * m)
  : mapping_(m)
  {
    const edm::ParameterSet esMap = conf.getParameter<edm::ParameterSet>("esMapping");
    es_mapping_ = new ESElectronicsMapper(esMap);
  }
  
  ~EcalRegionCabling(){
    // this pointer is own by this object.
    delete es_mapping_;
  }
  const EcalElectronicsMapping * mapping() const  { return mapping_;}
  const ESElectronicsMapper * es_mapping() const { return es_mapping_;}

  static uint32_t maxElementIndex() {return (FEDNumbering::MAXECALFEDID - FEDNumbering::MINECALFEDID +1);}
  static uint32_t maxESElementIndex() { return (FEDNumbering::MAXPreShowerFEDID - FEDNumbering::MINPreShowerFEDID +1);}

  static uint32_t elementIndex(const int FEDindex) {
    //do a test for the time being
    if (FEDindex > FEDNumbering::MAXECALFEDID || FEDindex < FEDNumbering::MINECALFEDID) {
      edm::LogError("IncorrectMapping")<<"FEDindex: "<< FEDindex
                    <<" is not between: "<<(int) FEDNumbering::MINECALFEDID
                    <<" and "<<(int)FEDNumbering::MAXECALFEDID;
      return 0;}
    uint32_t eI = FEDindex - FEDNumbering::MINECALFEDID;
    return eI; }

  static uint32_t esElementIndex(const int FEDindex) {
    //do a test for the time being
    if (FEDindex > FEDNumbering::MAXPreShowerFEDID || FEDindex < FEDNumbering::MINPreShowerFEDID) {
      edm::LogError("IncorrectMapping")<<"FEDindex: "<< FEDindex
                    <<" is not between: "<<(int) FEDNumbering::MINPreShowerFEDID
                    <<" and "<<(int)FEDNumbering::MAXPreShowerFEDID;
      return 0;}
    uint32_t eI = FEDindex - FEDNumbering::MINPreShowerFEDID;
    return eI; }

  static int fedIndex(const uint32_t index){ 
    int fI = index+FEDNumbering::MINECALFEDID; 
    return fI;}
    
  static int esFedIndex(const uint32_t index){ 
    int fI = index+FEDNumbering::MINPreShowerFEDID; 
    return fI;}

    
  uint32_t elementIndex(const double eta, const double phi) const{
    int FEDindex = mapping()->GetFED(eta,phi);
    return elementIndex(FEDindex); }

 private:
  const EcalElectronicsMapping * mapping_;
  const ESElectronicsMapper * es_mapping_;
};

#endif