HcalDDDRecConstants.h

#ifndef HcalCommonData_HcalDDDRecConstants_h
#define HcalCommonData_HcalDDDRecConstants_h

/** \class HcalDDDRecConstants
 *
 * this class reads the constant section of
 * the hcal-sim-numbering xml-file
 *  
 *  $Date: 2013/12/25 00:06:50 $
 * \author Sunanda Banerjee, SINP <sunanda.banerjee@cern.ch>
 *
 */

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

#include "DetectorDescription/Core/interface/DDsvalues.h"
#include "Geometry/HcalCommonData/interface/HcalCellType.h"
#include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"

class DDCompactView;    
class DDFilteredView;

class HcalDDDRecConstants {

public:

  HcalDDDRecConstants(const DDCompactView& cpv, 
		      const HcalDDDSimConstants& hcons);
  ~HcalDDDRecConstants();

  struct HcalID {
    int    subdet, eta, phi, depth;
    HcalID(int sub=0, int et=0, int fi=0, int d=0) : subdet(sub), eta(et),
						     phi(fi), depth(d) {}
  };
  struct HcalEtaBin {
    int    ieta, nPhi, depthStart;
    double etaMin, etaMax, phi0, dphi;
    std::vector<std::pair<int, int> > layer;
    HcalEtaBin(int eta=0, double et1=0, double et2=0, int nf=0, double fi0=0,
	       double df=0) : ieta(eta), nPhi(nf),depthStart(0), etaMin(et1), 
			      etaMax(et2), phi0(fi0), dphi(df) {}
  };
  struct HcalActiveLength {
    int    ieta, depth;
    double eta, thick;
    HcalActiveLength(int ie=0, int d=0, double et=0, 
		     double t=0) : ieta(ie), depth(d), eta(et), thick(t) {}
  };

  std::vector<std::pair<double,double> > getConstHBHE(const int type) const {
    if      (type == 0) return gconsHB;
    else if (type == 1) return gconsHE;
    else {std::vector<std::pair<double,double> > gcons; return gcons;}
  }
  std::vector<int>          getDepth(const int i) const {return layerGroup[i];}
  std::vector<HcalEtaBin>   getEtaBins(const int itype) const;
  std::pair<int,int>        getEtaRange(const int i) const
    {return std::pair<int,int>(iEtaMin[i],iEtaMax[i]);}
  std::vector<double>       getEtaTable()   const {return etaTable;}
  std::vector<double>       getEtaTableHF() const {return etaTableHF;}
  std::pair<double,double>  getEtaLimit(const int i) const 
    {return std::pair<double,double>(etaTable[i],etaTable[i+1]);}
  HcalID                    getHCID(int subdet, int ieta, int iphi, int lay,
				    int idepth) const;
  int                       getMaxDepth(const int type) const {return maxDepth[type];}
  int                       getNEta() const {return nEta;}
  double                    getPhiBin(const int i) const {return phibin[i];}
  double                    getPhiOff(const int i) const {return phioff[i];}
  std::vector<double>       getPhiOffs()    const {return phioff;}
  std::vector<double>       getPhiTable()   const {return phibin;}
  std::vector<double>       getPhiTableHF() const {return phibinHF;}
  std::vector<HcalActiveLength> getThickActive(const int type) const;
  std::string               getTopoMode() const {return modeTopo;}
  std::vector<HcalCellType> HcalCellTypes(HcalSubdetector) const;
  unsigned int              numberOfCells(HcalSubdetector) const;
  unsigned int              nCells(HcalSubdetector) const;
  unsigned int              nCells() const;
       
private:
  void                      initialize(const DDCompactView& cpv);
  void                      loadSpecPars(const DDFilteredView& fv);
  void                      loadSimConst();
  std::vector<double>       getDDDArray(const char *, const DDsvalues_type &, 
					int &) const;
  std::string               getDDDString(const std::string &, 
					 const DDsvalues_type &) const;

  bool                tobeInitialized;
  static const int    nEtaMax=100;
  const HcalDDDSimConstants & hcons;
  std::string         modeTopo;   // Mode for topology
  std::vector<double> phioff;     // Phi offset for barrel, endcap, forward
  std::vector<int>    etaGroup;   // Eta Grouping
  std::vector<std::pair<int,int> > etaSimValu; // eta ranges at Sim stage
  std::vector<double> etaTable;   // Eta table (HB+HE)
  std::vector<double> etaTableHF; // Eta table (HF)
  std::vector<int>    ietaMap;    // Map Sim level ieta to Rec level ieta
  int                 iEtaMin[4], iEtaMax[4]; // Minimum and maximum eta
  int                 maxDepth[4];// Maximum depth in HB/HE/HF/HO           
  int                 nEta;       // Number of bins in eta for HB and HE
  std::vector<int>    phiGroup;   // Eta Grouping
  std::vector<double> phibin;     // Phi step for all eta bins (HB, HE, HO)
  std::vector<double> phibinHF;   // Phi step for all eta bins (HF)
  std::vector<int>    phiUnitS;   // Phi unit at SIM stage
  std::vector<int>    layerGroup[nEtaMax];
  std::vector<int>    nOff;     // Speical eta bin #'s in barrel and endcap
  std::vector<std::pair<double,double> > gconsHB; // Geometry constatnts HB
  std::vector<std::pair<double,double> > gconsHE; // Geometry constatnts HE
  int                 nModule[2], nHalves[2];     // Modules, Halves for HB/HE
  enum { kHOSizePreLS1 = 2160, kHFSizePreLS1 = 1728 } ;
};

#endif