HcalDDDSimConstants.h

Go to the documentation of this file.

#ifndef Geometry_HcalTowerAlgo_HcalDDDSimConstants_h
#define Geometry_HcalTowerAlgo_HcalDDDSimConstants_h

#include<string>
#include<vector>
#include<iostream>
#include<iomanip>

#include "CondFormats/GeometryObjects/interface/HcalParameters.h"
#include "Geometry/HcalCommonData/interface/HcalCellType.h"
#include "Geometry/HcalCommonData/interface/HcalLayerDepthMap.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"

class HcalDDDSimConstants {

public:

  HcalDDDSimConstants(const HcalParameters* hp);
  ~HcalDDDSimConstants();

  HcalCellType::HcalCell    cell(const int det, const int zside, 
                 const int depth, const int etaR, 
                 const int iphi) const;
  int                       findDepth(const int det, const int eta, 
                      const int phi, const int zside,
                      const int lay) const;
  unsigned int              findLayer(int layer, const std::vector<HcalParameters::LayerItem>& layerGroup) const;
  std::vector<std::pair<double,double> > getConstHBHE(const int type) const;
  int                       getDepthEta16(const int det, const int phi, 
                      const int zside) const;
  int                       getDepthEta16M(const int det) const;
  int                       getDepthEta29(const int phi, int zside, int i) const;
  int                       getDepthEta29M(const int i, const bool planOne) const;
  std::pair<int,double>     getDetEta(const double eta, const int depth);
  int                       getEta(const int det, const int lay, const double hetaR);
  std::pair<int,int>        getEtaDepth(const int det, int etaR, int phi,
                    int zside, int depth, int lay);
  double                    getEtaHO(double& etaR, double& x, double& y, 
                     double& z) const;
  std::pair<int,int>        getiEtaRange(const int i) const {return std::pair<int,int>(hpar->etaMin[i],hpar->etaMax[i]);}
  const std::vector<double> &  getEtaTableHF() const {return hpar->etaTableHF;}
  const std::vector<double> &  getGparHF() const {return hpar->gparHF;}
  const std::vector<HcalDetId> & getIdHF2QIE() const {return idHF2QIE;}
  double                    getLayer0Wt(const int det, const int phi, 
                    const int zside) const;
  int                       getFrontLayer(const int det, const int eta) const;
  int                       getLayerFront(const int det, const int eta, 
                      const int phi, const int zside,
                      const int depth) const;
  int                       getLayerBack(const int det, const int eta, 
                     const int phi, const int zside,
                     const int depth) const;
  int                       getLayerMax(const int eta, const int depth) const;
  int                       getMaxDepth(const int type) const {return maxDepth[type];}
  int                       getMaxDepth(const int det, const int eta, 
                    const int phi, const int zside,
                    bool partialOnly) const;
  std::pair<int,int>        getMaxDepthDet(const int i) const {return ((i==1) ? depthMaxDf_ : depthMaxSp_);}
  int                       getMinDepth(const int det, const int eta, 
                    const int phi, const int zside,
                    bool partialOnly) const;
  std::pair<int,int>        getModHalfHBHE(const int type) const;
  std::pair<double,double>  getPhiCons(const int det, const int ieta) const;
  std::vector<std::pair<int,double> > getPhis(const int subdet, const int ieta) const;
  const std::vector<double> &  getPhiTableHF() const {return hpar->phitable;}
  const std::vector<double> &  getRTableHF()   const {return hpar->rTable;}
  std::vector<HcalCellType> HcalCellTypes()    const;
  std::vector<HcalCellType> HcalCellTypes(HcalSubdetector, int ieta=-1,
                      int depth=-1) const;
  bool                      isBH() const {return isBH_;}
  const HcalLayerDepthMap*  ldMap() const {return &ldmap_;}
  int                       maxHFDepth(const int ieta, const int iphi) const;
  unsigned int              numberOfCells(HcalSubdetector) const;
  int                       phiNumber(const int phi, const int unit) const;
  void                      printTiles() const;
  int                       unitPhi(const int det, const int etaR) const;
  int                       unitPhi(const double dphi) const; 
       
private:
  void                      initialize();
  double                    deltaEta(const int det, const int eta, 
                     const int depth) const;
  double                    getEta(const int det, const int etaR,
                   const int zside, const int depth=1) const;
  double                    getEta(const double r, const double z) const;
  int                       getShift(const HcalSubdetector subdet, 
                     const int depth) const;
  double                    getGain (const HcalSubdetector subdet, 
                     const int depth) const;
  void                      printTileHB(const int eta, const int phi,
                    const int zside, const int depth) const;
  void                      printTileHE(const int eta, const int phi, 
                    const int zside, const int depth) const;
  unsigned int              layerGroupSize(const int eta) const;
  unsigned int              layerGroup(const int eta, const int i) const;
  unsigned int              layerGroup(const int det, const int eta,
                       const int phi, const int zside,
                       const int i) const;

  const HcalParameters* hpar;
  HcalLayerDepthMap     ldmap_;

  static const int    nDepthMax=9;
  static const int    maxLayer_=18;
  static const int    maxLayerHB_=16;

  std::vector<int>    maxDepth; // Maximum depths in HB/HE/HF/HO
  int                 nEta;     // Number of bins in eta for HB and HE
  int                 nR;       // Number of bins in r
  int                 nPhiF;    // Number of bins in phitable
  std::vector<int>    depths[nDepthMax];   // Maximum layer number for a depth 
  int                 nDepth;   // Number of bins in depth0
  int                 nzHB, nmodHB;     // Number of halves and modules in HB
  int                 nzHE, nmodHE;     // Number of halves and modules in HE
  double              etaHO[4], rminHO; // eta in HO ring boundaries
  double              zVcal;    // Z-position  of the front of HF
  double              dzVcal;   // Half length of the HF
  double              dlShort;  // Diference of length between long and short
  int                 depthEta16[2]; // depth index of ieta=16 for HBmax, HEMin
  int                 depthEta29[2]; // maximum depth index for ieta=29
  int                 layFHB[2]; // first layers in HB (normal, special 16)
  int                 layFHE[3]; // first layers in HE (normal, special 16, 18)
  bool                isBH_;         // if HE is BH
  std::vector<HcalDetId> idHF2QIE;   // DetId's of HF modules with 2 QIE's
  std::pair<int,int>  depthMaxDf_, depthMaxSp_; // (subdet,maximum depth) default,special
};

#endif