EcalDQMCommonUtils.h

Go to the documentation of this file.

#ifndef EcalDQMCommonUtils_H
#define EcalDQMCommonUtils_H

#include <algorithm>
#include <cmath>
#include <iomanip>

#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"

#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EcalElectronicsId.h"
#include "DataFormats/EcalDetId/interface/EcalPnDiodeDetId.h"
#include "DataFormats/EcalDetId/interface/EcalScDetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
#include "DataFormats/EcalDetId/interface/EcalTriggerElectronicsId.h"

#include "FWCore/Utilities/interface/Exception.h"

namespace ecaldqm {

  enum SMName {
    kEEm07,
    kEEm08,
    kEEm09,
    kEEm01,
    kEEm02,
    kEEm03,
    kEEm04,
    kEEm05,
    kEEm06,
    kEBm01,
    kEBm02,
    kEBm03,
    kEBm04,
    kEBm05,
    kEBm06,
    kEBm07,
    kEBm08,
    kEBm09,
    kEBm10,
    kEBm11,
    kEBm12,
    kEBm13,
    kEBm14,
    kEBm15,
    kEBm16,
    kEBm17,
    kEBm18,
    kEBp01,
    kEBp02,
    kEBp03,
    kEBp04,
    kEBp05,
    kEBp06,
    kEBp07,
    kEBp08,
    kEBp09,
    kEBp10,
    kEBp11,
    kEBp12,
    kEBp13,
    kEBp14,
    kEBp15,
    kEBp16,
    kEBp17,
    kEBp18,
    kEEp07,
    kEEp08,
    kEEp09,
    kEEp01,
    kEEp02,
    kEEp03,
    kEEp04,
    kEEp05,
    kEEp06,
    kEEmLow = kEEm07,
    kEEmHigh = kEEm06,
    kEEpLow = kEEp07,
    kEEpHigh = kEEp06,
    kEBmLow = kEBm01,
    kEBmHigh = kEBm18,
    kEBpLow = kEBp01,
    kEBpHigh = kEBp18
  };

  static constexpr int nDCC = 54, nEBDCC = 36, nEEDCC = 18, nDCCMEM = 44, nEEDCCMEM = 8;

  // These lines set the number of TriggerTowers in "outer" and "inner" TCCs,
  // where "outer" := closer to the barrel. These constants are used in
  // setting the binning. There are 16 trigger towers per TCC for "outer" TCCs,
  // and 24 per TCC for "inner" TCCs (but the numbering is from 0 to 27, so
  // 28 bins are required).
  static constexpr int nTTOuter = 16, nTTInner = 28;

  static constexpr int nTCC = 108, kEEmTCCLow = 0, kEEmTCCHigh = 35, kEEpTCCLow = 72, kEEpTCCHigh = 107, kEBTCCLow = 36,
                       kEBTCCHigh = 71, nChannels = EBDetId::kSizeForDenseIndexing + EEDetId::kSizeForDenseIndexing,
                       nTowers = EcalTrigTowerDetId::kEBTotalTowers + EcalScDetId::kSizeForDenseIndexing;

  extern std::vector<unsigned> const memDCC;

  extern double const etaBound;

  // returns DCC ID (1 - 54)
  unsigned dccId(DetId const &, EcalElectronicsMapping const *);
  unsigned dccId(EcalElectronicsId const &);

  unsigned memDCCId(unsigned);     // convert from dccId skipping DCCs without MEM
  unsigned memDCCIndex(unsigned);  // reverse conversion

  // returns TCC ID (1 - 108)
  unsigned tccId(DetId const &, EcalElectronicsMapping const *);
  unsigned tccId(EcalElectronicsId const &, EcalElectronicsMapping const *);

  // returns the data tower id - pass only
  unsigned towerId(DetId const &, EcalElectronicsMapping const *);
  unsigned towerId(EcalElectronicsId const &);

  unsigned ttId(DetId const &, EcalElectronicsMapping const *);
  unsigned ttId(EcalElectronicsId const &, EcalElectronicsMapping const *);

  unsigned rtHalf(DetId const &, EcalElectronicsMapping const *);

  std::pair<unsigned, unsigned> innerTCCs(unsigned);
  std::pair<unsigned, unsigned> outerTCCs(unsigned);

  std::vector<DetId> scConstituents(EcalScDetId const &);

  EcalPnDiodeDetId pnForCrystal(DetId const &, char, EcalElectronicsMapping const *);

  unsigned dccId(std::string const &);
  std::string smName(unsigned);

  int zside(DetId const &);

  double eta(EBDetId const &, CaloGeometry const *);
  double eta(EEDetId const &, CaloGeometry const *);
  double phi(EBDetId const &);
  double phi(EEDetId const &);
  double phi(EcalTrigTowerDetId const &);
  double phi(double);

  bool isForward(DetId const &);

  bool isCrystalId(DetId const &);
  bool isSingleChannelId(DetId const &);
  bool isEcalScDetId(DetId const &);
  bool isEndcapTTId(DetId const &);

  unsigned nCrystals(unsigned);
  unsigned nSuperCrystals(unsigned);

  bool ccuExists(unsigned, unsigned);
}  // namespace ecaldqm

#endif