FastTimeTopology.h

Go to the documentation of this file.

#ifndef Geometry_CaloTopology_FastTimeTopology_h
#define Geometry_CaloTopology_FastTimeTopology_h 1
 
#include "DataFormats/ForwardDetId/interface/ForwardSubdetector.h"
#include "DataFormats/ForwardDetId/interface/FastTimeDetId.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/HGCalCommonData/interface/FastTimeDDDConstants.h"
#include <vector>
#include <iostream>
 
class FastTimeTopology : public CaloSubdetectorTopology {
public:
  FastTimeTopology(const FastTimeDDDConstants& hdcons, ForwardSubdetector subdet, int type);
 
  ~FastTimeTopology() override {}
 
  DetId goNorth(const DetId& id) const override { return offsetBy(id, 0, +1); }
  std::vector<DetId> north(const DetId& id) const override {
    DetId nextId = goNorth(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }
 
  DetId goSouth(const DetId& id) const override { return offsetBy(id, 0, -1); }
  std::vector<DetId> south(const DetId& id) const override {
    DetId nextId = goSouth(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }
 
  DetId goEast(const DetId& id) const override { return offsetBy(id, +1, 0); }
  std::vector<DetId> east(const DetId& id) const override {
    DetId nextId = goEast(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }
 
  DetId goWest(const DetId& id) const override { return offsetBy(id, -1, 0); }
  std::vector<DetId> west(const DetId& id) const override {
    DetId nextId = goWest(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }
 
  std::vector<DetId> up(const DetId& id) const override {
    std::vector<DetId> vNeighborsDetId;
    return vNeighborsDetId;
  }
 
  std::vector<DetId> down(const DetId& id) const override {
    std::vector<DetId> vNeighborsDetId;
    return vNeighborsDetId;
  }
 
  uint32_t detId2denseId(const DetId& id) const override;
  DetId denseId2detId(uint32_t denseId) const override;
  virtual uint32_t detId2denseGeomId(const DetId& id) const;
 
  bool valid(const DetId& id) const override;
  bool validHashIndex(uint32_t ix) const { return (ix < kSizeForDenseIndexing); }
 
  unsigned int totalModules() const { return kSizeForDenseIndexing; }
  unsigned int totalGeomModules() const { return (unsigned int)(2 * kHGeomHalf_); }
  int numberCells() const { return kHGeomHalf_; }
 
  const FastTimeDDDConstants& dddConstants() const { return hdcons_; }
 
  DetId offsetBy(const DetId startId, int nrStepsX, int nrStepsY) const;
  DetId switchZSide(const DetId startId) const;
 
  struct DecodedDetId {
    DecodedDetId() : iPhi(0), iEtaZ(0), zside(0), iType(0), subdet(0) {}
    int iPhi, iEtaZ, zside, iType, subdet;
  };
 
  DecodedDetId geomDenseId2decId(const uint32_t& hi) const;
  DecodedDetId decode(const DetId& id) const;
  DetId encode(const DecodedDetId& id_) const;
 
  ForwardSubdetector subDetector() const { return subdet_; }
  int detectorType() const { return type_; }
 
private:
  DetId changeXY(const DetId& id, int nrStepsX, int nrStepsY) const;
 
  const FastTimeDDDConstants& hdcons_;
  ForwardSubdetector subdet_;
  int type_, nEtaZ_, nPhi_, kHGhalf_, kHGeomHalf_;
  unsigned int kSizeForDenseIndexing;
};
 
#endif