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 "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"

class HcalDDDSimConstants {

public:

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

  HcalCellType::HcalCell    cell(int det, int zside, int depth, int etaR, 
                 int iphi) const;
  std::vector<std::pair<double,double> > getConstHBHE(const int type) const;
  int                       getDepthEta16(int i) const {return ((i==1) ? depthEta16[1] : depthEta16[0]);}
  int                       getDepthEta29(int i) const {return ((i==1) ? depthEta29[1] : depthEta29[0]);}
  std::pair<int,double>     getDetEta(double eta, int depth);
  int                       getEta(int det,int lay, double hetaR);
  std::pair<int,int>        getEtaDepth(int det, int etaR, int phi, 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;}
  unsigned int              findLayer(int layer, const std::vector<HcalParameters::LayerItem>& layerGroup) const;
  const std::vector<double> &  getGparHF() const {return hpar->gparHF;}
  double                    getLayer0Wt(int det, int phi, int zside) const;
  int                       getMaxDepth(const int type) const {return maxDepth[type];}
  std::pair<int,int>        getModHalfHBHE(const int type) const;
  std::pair<double,double>  getPhiCons(int det, int ieta);
  const std::vector<HcalDetId> & getIdHF2QIE() const {return idHF2QIE;}
  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_;}
  int                       maxHFDepth(int ieta, int iphi) const;
  unsigned int              numberOfCells(HcalSubdetector) const;
  int                       phiNumber(int phi, int unit) const;
  void                      printTiles() const;
  int                       unitPhi(int det, int etaR) const;
  int                       unitPhi(double dphi) const; 
       
private:
  void                      initialize();
  double                    deltaEta(int det, int eta, int depth) const;
  double                    getEta(int det, int etaR, int zside, int depth=1) const;
  double                    getEta(double r, double z) const;
  int                       getShift(HcalSubdetector subdet, int depth) const;
  double                    getGain (HcalSubdetector subdet, int depth) const;
  void                      printTileHB(int eta, int depth) const;
  void                      printTileHE(int eta, int depth) const;
  unsigned int              layerGroupSize( unsigned int eta ) const;
  unsigned int              layerGroup( unsigned int eta, unsigned int i ) const;
  
  const HcalParameters* hpar;

  static const int    nDepthMax=9;
  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
  bool                isBH_;         // if HE is BH
  std::vector<HcalDetId> idHF2QIE;   // DetId's of HF modules with 2 QIE's
};

#endif