HcalDDDRecConstants Class Reference

#include <HcalDDDRecConstants.h>

Classes
struct	HcalActiveLength
struct	HcalEtaBin
struct	HcalID
struct	HFCellParameters

Public Member Functions
std::vector< std::pair< double, double > >	getConstHBHE (const int type) const
const std::vector< int > &	getDepth (const unsigned int i) const
int	getDepthEta16 (int i) const
std::vector< HcalEtaBin >	getEtaBins (const int itype) const
std::pair< double, double >	getEtaLimit (const int i) const
std::pair< double, double >	getEtaPhi (int subdet, int ieta, int iphi) const
std::pair< int, int >	getEtaRange (const int i) const
const std::vector< double > &	getEtaTable () const
const std::vector< double > &	getEtaTableHF () const
HcalID	getHCID (int subdet, int ieta, int iphi, int lay, int idepth) const
std::vector< HFCellParameters >	getHFCellParameters () const
int	getMaxDepth (const int type) const
int	getMaxDepth (const int itype, const int ieta) const
int	getNEta () const
int	getNoff (const int i) const
double	getPhiBin (const int i) const
double	getPhiOff (const int i) const
const std::vector< double > &	getPhiOffs () const
const std::vector< double > &	getPhiTable () const
const std::vector< double > &	getPhiTableHF () const
double	getRZ (int subdet, int ieta, int depth) const
std::vector< HcalActiveLength >	getThickActive (const int type) const
int	getTopoMode () const
int	getTriggerMode () const
std::vector< HcalCellType >	HcalCellTypes (HcalSubdetector) const
	HcalDDDRecConstants (const HcalParameters *hp, const HcalDDDSimConstants &hc)
int	maxHFDepth (int ieta, int iphi) const
unsigned int	nCells (HcalSubdetector) const
unsigned int	nCells () const
unsigned int	numberOfCells (HcalSubdetector) const
	~HcalDDDRecConstants ()

Private Member Functions
void	initialize (void)
unsigned int	layerGroup (unsigned int eta, unsigned int i) const
unsigned int	layerGroupSize (unsigned int eta) const

Private Attributes
std::vector< std::pair< int, int > >	etaSimValu
std::vector< double >	etaTable
std::vector< std::pair< double, double > >	gconsHB
std::vector< std::pair< double, double > >	gconsHE
const HcalDDDSimConstants &	hcons
const HcalParameters *	hpar
std::vector< int >	ietaMap
std::vector< int >	iEtaMax
std::vector< int >	iEtaMin
std::vector< int >	maxDepth
int	nHalves [2]
int	nModule [2]
std::vector< double >	phibin
std::vector< int >	phiUnitS

Detailed Description

this class reads the constant section of the hcal-sim-numbering xml-file

Author

Sunanda Banerjee, SINP sunan.nosp@m.da.b.nosp@m.anerj.nosp@m.ee@c.nosp@m.ern.c.nosp@m.h

Definition at line 22 of file HcalDDDRecConstants.h.

Constructor & Destructor Documentation

HcalDDDRecConstants::HcalDDDRecConstants	(	const HcalParameters *	hp,
		const HcalDDDSimConstants &	hc
	)

Definition at line 13 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, and initialize().

                                                    : 
  hpar(hp), hcons(hc) {

#ifdef DebugLog
  std::cout << "HcalDDDRecConstants::HcalDDDRecConstants (const HcalParameters* hp) constructor" << std::endl;
#endif
  initialize();
}

HcalDDDRecConstants::~HcalDDDRecConstants

(

)

Definition at line 23 of file HcalDDDRecConstants.cc.

References gather_cfg::cout.

                                          { 
#ifdef DebugLog
  std::cout << "HcalDDDRecConstants::destructed!!!" << std::endl;
#endif
}

Member Function Documentation

std::vector<std::pair<double,double> > HcalDDDRecConstants::getConstHBHE ( const int  type ) const

inline

Definition at line 57 of file HcalDDDRecConstants.h.

References gconsHB, and gconsHE.

Referenced by HcalFlexiHardcodeGeometryLoader::makeHBCells(), and HcalFlexiHardcodeGeometryLoader::makeHECells().

                                                                        {
    if      (type == 0) return gconsHB;
    else if (type == 1) return gconsHE;
    else {std::vector<std::pair<double,double> > gcons; return gcons;}
  }

const std::vector< int > & HcalDDDRecConstants::getDepth

(

const unsigned int

i

)

const

Definition at line 570 of file HcalDDDRecConstants.cc.

References hpar, plotBeamSpotDB::last, and HcalParameters::layerGroupEtaRec.

Referenced by HcalTopology::HcalTopology().

                                                                               {
  std::vector<HcalParameters::LayerItem>::const_iterator last = hpar->layerGroupEtaRec.begin();
  for( std::vector<HcalParameters::LayerItem>::const_iterator it = hpar->layerGroupEtaRec.begin(); it != hpar->layerGroupEtaRec.end(); ++it ) {
    if( it->layer == i + 1 )
      return it->layerGroup;
    if( it->layer > i + 1 )
      return last->layerGroup;
    last = it;
  }
  return last->layerGroup;
}

int HcalDDDRecConstants::getDepthEta16 ( int  i ) const

inline

Definition at line 63 of file HcalDDDRecConstants.h.

References HcalDDDSimConstants::getDepthEta16(), and hcons.

Referenced by HcalTopology::depthBinInformation(), and HcalTopology::validRaw().

{return hcons.getDepthEta16(i);}

std::vector< HcalDDDRecConstants::HcalEtaBin > HcalDDDRecConstants::getEtaBins

(

const int

itype

)

const

Definition at line 30 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, HcalDDDRecConstants::HcalEtaBin::depthStart, etaTable, plotBeamSpotDB::first, HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta29(), hcons, hpar, i, iEtaMax, iEtaMin, relval_2017::k, cmsLHEtoEOSManager::l, HcalDDDRecConstants::HcalEtaBin::layer, layerGroup(), layerGroupSize(), gen::n, nModule, HcalParameters::noff, phibin, HcalParameters::phioff, and edm::second().

Referenced by getThickActive(), HcalCellTypes(), HcalTopology::HcalTopology(), HcalFlexiHardcodeGeometryLoader::makeHBCells(), HcalFlexiHardcodeGeometryLoader::makeHECells(), and nCells().

                                                     {

  std::vector<HcalDDDRecConstants::HcalEtaBin> bins;
  unsigned int type  = (itype == 0) ? 0 : 1;
  unsigned int lymax = (type == 0) ? 17 : 19;
  for (int ieta = iEtaMin[type]; ieta <= iEtaMax[type]; ++ieta) {
    int nfi = (int)((20.001*nModule[itype]*CLHEP::deg)/phibin[ieta-1]);
    HcalDDDRecConstants::HcalEtaBin etabin = HcalDDDRecConstants::HcalEtaBin(ieta, etaTable[ieta-1], etaTable[ieta], nfi, hpar->phioff[type], phibin[ieta-1]);
    int n = (ieta == iEtaMax[type]) ? 0 : 1;
    HcalDDDRecConstants::HcalEtaBin etabin0= HcalDDDRecConstants::HcalEtaBin(ieta, etaTable[ieta-1], etaTable[ieta+n], nfi, hpar->phioff[type], phibin[ieta-1]);
    etabin0.depthStart = hcons.getDepthEta29(0)+1;
    int dstart = -1;
    if (layerGroupSize(ieta-1) > 0) {
      int lmin(0), lmax(0);
      int dep = layerGroup(ieta-1, 0);
      if (type == 1 && ieta == iEtaMin[1]) dep = hcons.getDepthEta16(1);
      unsigned lymx0 = (layerGroupSize(ieta-1) > lymax) ? lymax : layerGroupSize(ieta-1);
      for (unsigned int l=0; l<lymx0; ++l) {
    if ((int)layerGroup( ieta-1, l) == dep) {
      if (lmin == 0) lmin = l + 1;
      lmax = l + 1;
    } else if ((int)layerGroup( ieta-1, l ) > dep) {
      if (dstart < 0) dstart = dep;
      if (type == 1 && ieta+1 == hpar->noff[1] && dep > hcons.getDepthEta29(0)) {
        etabin0.layer.push_back(std::pair<int,int>(lmin,lmax));
      } else {
        etabin.layer.push_back(std::pair<int,int>(lmin,lmax));
      }
      lmin = (l + 1);
      lmax = l;
      dep  = layerGroup(ieta-1, l);
    }
    if (type == 0 && ieta == iEtaMax[type] && dep > hcons.getDepthEta16(0)) break;
    if (type == 1 && ieta == hpar->noff[1] && dep > hcons.getDepthEta29(0)){
      lmax = lymx0;
      break;
    }
      }
      if (lmax >= lmin) {
    if (ieta+1 == hpar->noff[1]) {
      etabin0.layer.push_back(std::pair<int,int>(lmin,lmax));
      bins.push_back(etabin0);
    } else if (ieta == hpar->noff[1]) {
    } else {
      etabin.layer.push_back(std::pair<int,int>(lmin,lmax));
      if (dstart < 0) dstart = dep;
    }
      }
    }
    etabin.depthStart = dstart;
    bins.push_back(etabin);
  }
#ifdef DebugLog
  std::cout << "Prepares " << bins.size() << " eta bins for type " << type 
        << std::endl;
  for (unsigned int i=0; i<bins.size(); ++i) {
    std::cout << "Bin[" << i << "]: Eta = (" << bins[i].ieta << ":"
          << bins[i].etaMin << ":" << bins[i].etaMax << ") Phi = (" 
          << bins[i].nPhi << ":" << bins[i].phi0 << ":" << bins[i].dphi 
          << ") and " << bins[i].layer.size() << " depths (start) "
          << bins[i].depthStart << " :";
    for (unsigned int k=0; k<bins[i].layer.size(); ++k)
      std::cout << " [" << k << "] " << bins[i].layer[k].first << ":"
        << bins[i].layer[k].second;
    std::cout << std::endl;
  }
#endif
  return bins;
}

std::pair<double,double> HcalDDDRecConstants::getEtaLimit ( const int  i ) const

inline

Definition at line 70 of file HcalDDDRecConstants.h.

References etaTable, and i.

    {return std::pair<double,double>(etaTable[i],etaTable[i+1]);}

std::pair< double, double > HcalDDDRecConstants::getEtaPhi	(	int	subdet,
		int	ieta,
		int	iphi
	)		const

Definition at line 101 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, eta, etaTable, HcalParameters::etaTableHF, HcalBarrel, HcalEndcap, HcalOuter, hpar, iEtaMax, iEtaMin, M_PI, phi, phibin, HcalParameters::phioff, HcalParameters::phitable, and csvLumiCalc::unit.

Referenced by HGCalSimHitValidation::analyzeHits().

                                                                   {
  int ietaAbs = (ieta > 0) ? ieta : -ieta;
  const double fiveDegInRad = 2*M_PI/72;
  double eta(0), phi(0);
  if ((subdet == static_cast<int>(HcalBarrel)) || 
      (subdet == static_cast<int>(HcalEndcap)) ||
      (subdet == static_cast<int>(HcalOuter))) {  // Use Eta Table
    int unit    = (int)(phibin[ietaAbs-1]/fiveDegInRad+0.5);
    int kphi    = (unit == 2) ? ((iphi-1)/2 + 1) : iphi;
    double foff = (ietaAbs <= iEtaMax[0]) ? hpar->phioff[0] : hpar->phioff[1];
    eta         = 0.5*(etaTable[ietaAbs-1]+etaTable[ietaAbs]);
    phi         = foff + (kphi-0.5)*phibin[ietaAbs-1];
  } else {
    ietaAbs    -= iEtaMin[2];
    int unit    = (int)(hpar->phitable[ietaAbs-1]/fiveDegInRad+0.5);
    int kphi    = (unit == 4) ? ((iphi-3)/4 + 1) : ((iphi-1)/2 + 1);
    double foff = (unit > 2) ? hpar->phioff[4] : hpar->phioff[2];
    eta         = 0.5*(hpar->etaTableHF[ietaAbs-1]+hpar->etaTableHF[ietaAbs]);
    phi         = foff + (kphi-0.5)*hpar->phitable[ietaAbs-1];
  }
  if (ieta < 0)   eta  = -eta;
  if (phi > M_PI) phi -= (2*M_PI);
#ifdef DebugLog
  std::cout << "getEtaPhi: subdet|ieta|iphi " << subdet << "|" << ieta << "|"
        << iphi << " eta|phi " << eta << "|" << phi << std::endl;
#endif
  return std::pair<double,double>(eta,phi);
}

std::pair<int,int> HcalDDDRecConstants::getEtaRange ( const int  i ) const

inline

Definition at line 66 of file HcalDDDRecConstants.h.

References i, iEtaMax, and iEtaMin.

Referenced by SimHitsValidationHcal::getLimits(), and HcalTopology::HcalTopology().

    {return std::pair<int,int>(iEtaMin[i],iEtaMax[i]);}

const std::vector<double>& HcalDDDRecConstants::getEtaTable ( ) const

inline

Definition at line 68 of file HcalDDDRecConstants.h.

References etaTable.

Referenced by HcalTopology::HcalTopology(), and CaloTowerHardcodeGeometryLoader::load().

{return etaTable;}

const std::vector<double>& HcalDDDRecConstants::getEtaTableHF ( ) const

inline

Definition at line 69 of file HcalDDDRecConstants.h.

References HcalParameters::etaTableHF, and hpar.

Referenced by HcalTopology::HcalTopology(), and CaloTowerHardcodeGeometryLoader::load().

{return hpar->etaTableHF;}

HcalDDDRecConstants::HcalID HcalDDDRecConstants::getHCID	(	int	subdet,
		int	ieta,
		int	iphi,
		int	lay,
		int	idepth
	)		const

Definition at line 131 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, HLT_FULL_cff::depth, eta, HcalDDDSimConstants::getDepthEta16(), HcalBarrel, HcalEndcap, HcalOuter, hcons, hpar, ietaMap, iEtaMin, layerGroup(), HcalParameters::noff, phi, phibin, HcalParameters::phigroup, HcalDDDSimConstants::phiNumber(), HcalParameters::phioff, phiUnitS, csvLumiCalc::unit, and HcalDDDSimConstants::unitPhi().

Referenced by SimHitsValidationHcal::analyze(), HGCalSimHitValidation::analyze(), and HcalHitRelabeller::relabel().

                                   {

  int    eta(ieta), phi(iphi), depth(idepth);
  if ((subdet == static_cast<int>(HcalOuter)) ||
      ((subdet == static_cast<int>(HcalBarrel)) && (lay > 17))) {
    subdet= static_cast<int>(HcalOuter);
    depth = 4;
  } else if (subdet == static_cast<int>(HcalBarrel) || 
         subdet == static_cast<int>(HcalEndcap)) {
    eta      = ietaMap[ieta-1];
    int unit = phiUnitS[ieta-1];
    int phi0 = (iphi-1)/(hpar->phigroup[eta-1]);
    if (unit == 2) {
      phi0   = (iphi+1)/2;
      phi0   = (phi0-1)/(hpar->phigroup[eta-1]);
    } else if (unit == 4) {
      phi0   = (iphi+5)/4;
      phi0   = (phi0-1)/(hpar->phigroup[eta-1]);
    }
    ++phi0;
    unit     = hcons.unitPhi(phibin[eta-1]);
    phi      = hcons.phiNumber(phi0,unit);
    depth    = layerGroup( eta-1, lay-1 );
    if (eta == iEtaMin[1]) {
      if (subdet == static_cast<int>(HcalBarrel)) {
    if (depth > hcons.getDepthEta16(0)) depth = hcons.getDepthEta16(0);
      } else {
    if (depth < hcons.getDepthEta16(1)) depth = hcons.getDepthEta16(1);
      }
    } else if (eta == hpar->noff[0] && lay > 1) {
      int   kphi   = phi + int((hpar->phioff[3]+0.1)/phibin[eta-1]);
      kphi         = (kphi-1)%4 + 1;
      if (kphi == 2 || kphi == 3) depth = layerGroup( eta-1, lay-2 );
    } else if (eta == hpar->noff[1] && depth > 2) {
       eta = hpar->noff[1]-1;
    }
  } 
#ifdef DebugLog
  std::cout << "getHCID: input " << subdet << ":" << ieta << ":" << iphi
        << ":" << idepth << ":" << lay << " output " << eta << ":" << phi
        << ":" << depth << std::endl;
#endif
  return HcalDDDRecConstants::HcalID(subdet,eta,phi,depth);
}

std::vector< HcalDDDRecConstants::HFCellParameters > HcalDDDRecConstants::getHFCellParameters

(

)

const

Definition at line 178 of file HcalDDDRecConstants.cc.

References funct::abs(), gather_cfg::cout, HLT_FULL_cff::depth, HcalDDDSimConstants::getIdHF2QIE(), HcalDDDSimConstants::getPhiTableHF(), HcalDDDSimConstants::getRTableHF(), hcons, iEtaMin, relval_2017::k, maxDepth, HLT_FULL_cff::nEta, and HLT_FULL_cff::nPhi.

Referenced by HcalFlexiHardcodeGeometryLoader::makeHFCells().

                                               {

  std::vector<HcalDDDRecConstants::HFCellParameters> cells;
  unsigned int nEta = hcons.getPhiTableHF().size();
  if (maxDepth[2] > 0) {
    for (unsigned int k=0; k<nEta; ++k) {
      int ieta = iEtaMin[2] + k;
      int dphi = (int)(0.001 + hcons.getPhiTableHF()[k]/(5.0*CLHEP::deg));
      int iphi = (dphi == 4) ? 3 : 1;
      int nphi = 72/dphi;
      double rMin = hcons.getRTableHF()[nEta-k-1]/CLHEP::cm;
      double rMax = hcons.getRTableHF()[nEta-k]/CLHEP::cm;
      HcalDDDRecConstants::HFCellParameters cell1( ieta,1,iphi,dphi,nphi,rMin,rMax);
      cells.push_back(cell1);
      HcalDDDRecConstants::HFCellParameters cell2(-ieta,1,iphi,dphi,nphi,rMin,rMax);
      cells.push_back(cell2);
    }
  }
  if (maxDepth[2] > 2) {
    if (hcons.getIdHF2QIE().size() > 0) {
      for (unsigned int k=0; k<hcons.getIdHF2QIE().size(); ++k) {
    int ieta = hcons.getIdHF2QIE()[k].ieta();
    int ind  = std::abs(ieta) - iEtaMin[2];
    int dphi = (int)(0.001 + hcons.getPhiTableHF()[ind]/(5.0*CLHEP::deg));
    int iphi = hcons.getIdHF2QIE()[k].iphi();
    double rMin = hcons.getRTableHF()[nEta-ind-1]/CLHEP::cm;
    double rMax = hcons.getRTableHF()[nEta-ind]/CLHEP::cm;
    HcalDDDRecConstants::HFCellParameters cell1( ieta,3,iphi,dphi,1,rMin,rMax);
    cells.push_back(cell1);
      }
    } else {
      for (unsigned int k=0; k<nEta; ++k) {
    int ieta = iEtaMin[2] + k;
    int dphi = (int)(0.001 + hcons.getPhiTableHF()[k]/(5.0*CLHEP::deg));
    int iphi = (dphi == 4) ? 3 : 1;
    int nphi = 72/dphi;
    double rMin = hcons.getRTableHF()[nEta-k-1]/CLHEP::cm;
    double rMax = hcons.getRTableHF()[nEta-k]/CLHEP::cm;
    HcalDDDRecConstants::HFCellParameters cell1( ieta,3,iphi,dphi,nphi,rMin,rMax);
    cells.push_back(cell1);
    HcalDDDRecConstants::HFCellParameters cell2(-ieta,3,iphi,dphi,nphi,rMin,rMax);
    cells.push_back(cell2);
      }
    }
  }
#ifdef DebugLog
  std::cout << "HcalDDDRecConstants returns " << cells.size() 
        << " HF cell parameters" << std::endl;
  for (unsigned int k=0; k<cells.size(); ++k)
    std::cout << "Cell[" << k <<"] : (" << cells[k].ieta <<", "<< cells[k].depth
          << ", " << cells[k].firstPhi << ", " << cells[k].stepPhi << ", "
          << cells[k].nPhi << ", " << cells[k].rMin << ", "
          << cells[k].rMax << ")" << std::endl;
#endif
  return cells;
}

int HcalDDDRecConstants::getMaxDepth ( const int  type ) const

inline

Definition at line 75 of file HcalDDDRecConstants.h.

References maxDepth.

Referenced by HGCalDigiClient::beginRun(), HGCalRecHitsClient::beginRun(), HGCalSimHitsClient::beginRun(), HcalSimHitsClient::beginRun(), SimHitsValidationHcal::bookHistograms(), HcalTopology::depthBinInformation(), HGCalDigiValidation::dqmBeginRun(), HGCalRecHitValidation::dqmBeginRun(), HGCalSimHitValidation::dqmBeginRun(), HcalTopology::HcalTopology(), HcalTopologyIdealEP::produce(), and HcalTopology::validRaw().

{return maxDepth[type];}

int HcalDDDRecConstants::getMaxDepth	(	const int	itype,
		const int	ieta
	)		const

Definition at line 235 of file HcalDDDRecConstants.cc.

References HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta29(), hcons, hpar, iEtaMax, layerGroup(), layerGroupSize(), and HcalParameters::noff.

                                                                           {

  int lmax(0);
  unsigned int type  = (itype == 0) ? 0 : 1;
  unsigned int lymax = (type == 0) ? 17 : 19;
  if (layerGroupSize(ieta-1) > 0) {
    if (layerGroupSize(ieta-1) < lymax) lymax = layerGroupSize(ieta-1);
    lmax = (int)(layerGroup(ieta-1, lymax-1));
    if (type == 0 && ieta == iEtaMax[type]) lmax = hcons.getDepthEta16(0);
    if (type == 1 && ieta >= hpar->noff[1]) lmax = hcons.getDepthEta29(0);
  }
  return lmax;
}

int HcalDDDRecConstants::getNEta ( ) const

inline

Definition at line 77 of file HcalDDDRecConstants.h.

References HcalParameters::etagroup, and hpar.

Referenced by HcalTopology::HcalTopology().

{return hpar->etagroup.size();}

int HcalDDDRecConstants::getNoff ( const int  i ) const

inline

Definition at line 78 of file HcalDDDRecConstants.h.

References hpar, i, and HcalParameters::noff.

Referenced by HcalTopology::validRaw().

{return hpar->noff[i];}

double HcalDDDRecConstants::getPhiBin ( const int  i ) const

inline

Definition at line 79 of file HcalDDDRecConstants.h.

References i, and phibin.

{return phibin[i];}

double HcalDDDRecConstants::getPhiOff ( const int  i ) const

inline

Definition at line 80 of file HcalDDDRecConstants.h.

References hpar, i, and HcalParameters::phioff.

{return hpar->phioff[i];}

const std::vector<double>& HcalDDDRecConstants::getPhiOffs ( ) const

inline

Definition at line 81 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::phioff.

Referenced by HcalTopology::HcalTopology().

{return hpar->phioff;}

const std::vector<double>& HcalDDDRecConstants::getPhiTable ( ) const

inline

Definition at line 82 of file HcalDDDRecConstants.h.

References phibin.

Referenced by HcalTopology::HcalTopology().

{return phibin;}

const std::vector<double>& HcalDDDRecConstants::getPhiTableHF ( ) const

inline

Definition at line 83 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::phitable.

Referenced by HcalTopology::HcalTopology().

{return hpar->phitable;}

double HcalDDDRecConstants::getRZ	(	int	subdet,
		int	ieta,
		int	depth
	)		const

Definition at line 249 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, HLT_FULL_cff::etaMax, plotBeamSpotDB::first, gconsHB, gconsHE, HcalBarrel, relval_2017::k, layerGroup(), and layerGroupSize().

Referenced by HGCalSimHitValidation::analyzeHits().

                                                                       {

  int ietaAbs = (ieta > 0) ? ieta : -ieta;
  double rz(0);
#ifdef DebugLog
  int    lay(0);
#endif
  if (ietaAbs < hpar->etaMax[1]) {
    for (unsigned int k=0; k< layerGroupSize( ietaAbs-1 ); ++k) {
      if (depth == (int)layerGroup( ietaAbs-1, k )) {
    rz = ((subdet == static_cast<int>(HcalBarrel)) ? (gconsHB[k].first) :
          (gconsHE[k].first));
    if (rz > 10.) {
#ifdef DebugLog
      lay = k;
#endif
      break;
    }
      }
    }
  }
#ifdef DebugLog
  std::cout << "getRZ: subdet|ieta|depth " << subdet << "|" << ieta << "|"
        << depth << " lay|rz " << lay << "|" << rz << std::endl;
#endif
  return rz;
}

std::vector< HcalDDDRecConstants::HcalActiveLength > HcalDDDRecConstants::getThickActive

(

const int

type

)

const

Definition at line 278 of file HcalDDDRecConstants.cc.

References funct::cos(), gather_cfg::cout, HcalDDDRecConstants::HcalActiveLength::depth, HLT_FULL_cff::depth, eta, HcalDDDRecConstants::HcalActiveLength::eta, create_public_lumi_plots::exp, plotBeamSpotDB::first, gconsHB, gconsHE, getEtaBins(), i, HcalDDDRecConstants::HcalActiveLength::ieta, j, relval_2017::k, pileupReCalc_HLTpaths::scale, edm::second(), funct::sin(), theta(), and HcalDDDRecConstants::HcalActiveLength::thick.

Referenced by HcalRaddamMuon::beginRun(), and HcalHBHEMuonAnalyzer::beginRun().

                                                        {

  std::vector<HcalDDDRecConstants::HcalActiveLength> actives;
  std::vector<HcalDDDRecConstants::HcalEtaBin> bins = getEtaBins(type);
#ifdef DebugLog
  unsigned int kount(0);
#endif
  for (unsigned int k=0; k<bins.size(); ++k) {
    int    ieta  = bins[k].ieta;
    double eta   = 0.5*(bins[k].etaMin+bins[k].etaMax);
    double theta = 2*atan(exp(-eta));
    double scale = 1.0/((type == 0) ? sin(theta) : cos(theta));
    int    depth = bins[k].depthStart;
    for (unsigned int i = 0; i < bins[k].layer.size(); ++i) {
      double thick(0);
      for (int j = bins[k].layer[i].first; j <= bins[k].layer[i].second; ++j) {
    if (type == 0 || j > 1) 
      thick += ((type == 0) ? gconsHB[j-1].second : gconsHE[j-1].second);
      }
      thick *= (2.*scale);
      HcalDDDRecConstants::HcalActiveLength active(ieta,depth,eta,thick);
      actives.push_back(active);
      ++depth;
#ifdef DebugLog
      kount++;
      std::cout << "getThickActive: [" << kount << "] eta:" << active.ieta 
        << ":" << active.eta << " depth " << active.depth << " thick " 
        << active.thick << std::endl;
#endif
    }
  }
  return actives;
}

int HcalDDDRecConstants::getTopoMode ( ) const

inline

Definition at line 86 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::topologyMode.

Referenced by HcalTopology::HcalTopology(), and HcalTopologyIdealEP::produce().

{return ((hpar->topologyMode)&0xFF);}

int HcalDDDRecConstants::getTriggerMode ( ) const

inline

Definition at line 87 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::topologyMode.

Referenced by HcalTopology::HcalTopology().

{return (((hpar->topologyMode)>>8)&0xFF);}

std::vector< HcalCellType > HcalDDDRecConstants::HcalCellTypes

(

HcalSubdetector

subdet

)

const

Definition at line 313 of file HcalDDDRecConstants.cc.

References newFWLiteAna::bin, KineDebug3::count(), gather_cfg::cout, HLT_FULL_cff::depth, HLT_FULL_cff::etaMax, HLT_FULL_cff::etaMin, etaSimValu, plotBeamSpotDB::first, getEtaBins(), HcalBarrel, HcalDDDSimConstants::HcalCellTypes(), HcalEndcap, hcons, hpar, HLT_FULL_cff::nPhi, HcalParameters::phioff, groupFilesInBlocks::temp, csvLumiCalc::unit, and HcalDDDSimConstants::unitPhi().

Referenced by HcalDDDGeometryLoader::fill(), and numberOfCells().

                                                               {

  if (subdet == HcalBarrel || subdet == HcalEndcap) {
    std::vector<HcalCellType> cells;
    int isub   = (subdet == HcalBarrel) ? 0 : 1;
    std::vector<HcalDDDRecConstants::HcalEtaBin> etabins = getEtaBins(isub);
    for (unsigned int bin=0; bin<etabins.size(); ++bin) {
      std::vector<HcalCellType> temp;
      std::vector<int>          count;
      std::vector<double>       dmin, dmax;
      for (unsigned int il=0; il<etabins[bin].layer.size(); ++il) {
    HcalCellType cell(subdet, 0, 0, 0, HcalCellType::HcalCell());
    temp.push_back(cell);
    count.push_back(0);
    dmin.push_back(0);
    dmax.push_back(0);
      }
      int ieta = etabins[bin].ieta;
      for (int keta=etaSimValu[ieta-1].first; keta<=etaSimValu[ieta-1].second;
       ++keta) {
    std::vector<HcalCellType> cells = hcons.HcalCellTypes(subdet,keta,-1);
    for (unsigned int ic=0; ic<cells.size(); ++ic) {
      for (unsigned int il=0; il<etabins[bin].layer.size(); ++il) {
        if (cells[ic].depthSegment() >= etabins[bin].layer[il].first &&
        cells[ic].depthSegment() <= etabins[bin].layer[il].second) {
          if (count[il] == 0) {
        temp[il] = cells[ic];
        dmin[il] = cells[ic].depthMin();
        dmax[il] = cells[ic].depthMax();
          }
          ++count[il];
          if (cells[ic].depthMin() < dmin[il]) dmin[il] = cells[ic].depthMin();
          if (cells[ic].depthMax() > dmax[il]) dmax[il] = cells[ic].depthMax();
          break;
        }
      }
    }
      }
      int unit = hcons.unitPhi(etabins[bin].dphi);
      for (unsigned int il=0; il<etabins[bin].layer.size(); ++il) {
    int depth = etabins[bin].depthStart + (int)(il);
    temp[il].setEta(ieta,etabins[bin].etaMin,etabins[bin].etaMax);
    temp[il].setPhi(etabins[bin].nPhi,unit,etabins[bin].dphi/CLHEP::deg,
            hpar->phioff[isub]/CLHEP::deg);
    temp[il].setDepth(depth,dmin[il],dmax[il]);
    cells.push_back(temp[il]);
      }
    }
#ifdef DebugLog
    std::cout << "HcalDDDRecConstants: found " << cells.size() << " cells for sub-detector type " << isub << std::endl;
    for (unsigned int ic=0; ic<cells.size(); ++ic)
      std::cout << "Cell[" << ic << "] " << cells[ic] << std::endl;
#endif
    return cells;
  } else {
    return hcons.HcalCellTypes(subdet,-1,-1);
  }
}

void HcalDDDRecConstants::initialize ( void  )

private

Definition at line 417 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, HcalParameters::drHB, HcalParameters::dzHE, HcalParameters::etagroup, HcalParameters::etaMax, HLT_FULL_cff::etaMax, HcalParameters::etaMin, HLT_FULL_cff::etaMin, etaSimValu, HcalParameters::etaTable, etaTable, HcalParameters::etaTableHF, Exception, plotBeamSpotDB::first, citk::for(), gconsHB, gconsHE, HcalDDDSimConstants::getDepthEta16(), hcons, hpar, i, ietaMap, iEtaMax, iEtaMin, relval_2017::k, layerGroup(), layerGroupSize(), HcalParameters::maxDepth, maxDepth, HcalParameters::modHB, HcalParameters::modHE, HLT_FULL_cff::nEta, nHalves, nModule, HcalParameters::phibin, phibin, HcalParameters::phigroup, HcalParameters::phitable, phiUnitS, HcalParameters::rHB, edm::second(), csvLumiCalc::unit, HcalDDDSimConstants::unitPhi(), and HcalParameters::zHE.

Referenced by HcalDDDRecConstants().

                                         {

  //Eta grouping
  int nEta  = (int)(hpar->etagroup.size());
  if (nEta != (int)(hpar->phigroup.size())) {
    edm::LogError("HCalGeom") << "HcalDDDRecConstants: sizes of the vectors "
                              << " etaGroup (" << nEta << ") and phiGroup ("
                              << hpar->phigroup.size() << ") do not match";
    throw cms::Exception("DDException") << "HcalDDDRecConstants: inconsistent array sizes" << nEta << ":" << hpar->phigroup.size();
  }

  // First eta table
  iEtaMin     = hpar->etaMin;
  iEtaMax     = hpar->etaMax;
  etaTable.clear(); ietaMap.clear(); etaSimValu.clear();
  int ieta(0), ietaHB(0), ietaHE(0);
  etaTable.push_back(hpar->etaTable[ieta]);
  for (int i=0; i<nEta; ++i) {
    int ef = ieta+1;
    ieta  += (hpar->etagroup[i]);
    if (ieta >= (int)(hpar->etaTable.size())) {
      edm::LogError("HCalGeom") << "Going beyond the array boundary "
                << hpar->etaTable.size() << " at index " << i 
                << " of etaTable from SimConstant";
      throw cms::Exception("DDException") << "Going beyond the array boundary "
                      << hpar->etaTable.size() 
                      << " at index " << i 
                      << " of etaTable from SimConstant";
    } else {
      etaTable.push_back(hpar->etaTable[ieta]);
      etaSimValu.push_back(std::pair<int,int>(ef,ieta));
    }
    for (int k=0; k<(hpar->etagroup[i]); ++k) ietaMap.push_back(i+1);
    if (ieta <= hpar->etaMax[0]) ietaHB = i+1;
    if (ieta <= hpar->etaMin[1]) ietaHE = i+1;
  }
  iEtaMin[1] = ietaHE;
  iEtaMax[0] = ietaHB;

  // Then Phi bins
  ieta = 0;
  phibin.clear(); phiUnitS.clear();
  for (int i=0; i<nEta; ++i) {
    double dphi = (hpar->phigroup[i])*(hpar->phibin[ieta]);
    phibin.push_back(dphi);
    ieta += (hpar->etagroup[i]);
  }
  for (unsigned int i=1; i<hpar->etaTable.size(); ++i) {
    int unit = hcons.unitPhi(hpar->phibin[i-1]);
    phiUnitS.push_back(unit);
  }

#ifdef DebugLog
  std::cout << "Modified eta/deltaphi table for " << nEta << " bins" << std::endl;
  for (int i=0; i<nEta; ++i)
    std::cout << "Eta[" << i << "] = " << etaTable[i] << ":" << etaTable[i+1]
          << ":" << etaSimValu[i].first << ":" << etaSimValu[i].second
          << " PhiBin[" << i << "] = " << phibin[i]/CLHEP::deg <<std::endl;
  std::cout << "PhiUnitS";
  for (unsigned int i=0; i<phiUnitS.size(); ++i)
    std::cout << " [" << i << "] = " << phiUnitS[i];
  std::cout << std::endl;
  std::cout << "EtaTableHF";
  for (unsigned int i=0; i<hpar->etaTableHF.size(); ++i)
    std::cout << " [" << i << "] = " << hpar->etaTableHF[i];
  std::cout << std::endl;
  std::cout << "PhiBinHF";
  for (unsigned int i=0; i<hpar->phitable.size(); ++i)
    std::cout << " [" << i << "] = " << hpar->phitable[i];
  std::cout << std::endl;
#endif

  //Now the depths
  maxDepth    = hpar->maxDepth;
  maxDepth[0] = maxDepth[1] = 0;
  for (int i=0; i<nEta; ++i) {
    unsigned int imx = layerGroupSize(i);
    int laymax = (imx > 0) ? layerGroup(i,imx-1) : 0;
    if (i < iEtaMax[0]) {
      int laymax0 = (imx > 16) ? layerGroup(i,16) : laymax;
      if (i+1 == iEtaMax[0]) laymax0 = hcons.getDepthEta16(0);
#ifdef DebugLog
      std::cout << "HB " << i << " " << imx << " " << laymax << " " << laymax0 << std::endl;
#endif
      if (maxDepth[0] < laymax0) maxDepth[0] = laymax0;
    }
    if (i >= iEtaMin[1]-1 && i < iEtaMax[1]) {
#ifdef DebugLog
      std::cout << "HE " << i << " " << imx << " " << laymax << std::endl;
#endif
      if (maxDepth[1] < laymax) maxDepth[1] = laymax;
    }
  }
#ifdef DebugLog
  for (int i=0; i<4; ++i) 
    std::cout << "Detector Type[" << i << "] iEta " << iEtaMin[i] << ":"
              << iEtaMax[i] << " MaxDepth " << maxDepth[i] << std::endl; 
#endif

  //Now the geometry constants
  nModule[0] = hpar->modHB[0];
  nHalves[0] = hpar->modHB[1];
  for (unsigned int i=0; i<hpar->rHB.size(); ++i) {
    gconsHB.push_back(std::pair<double,double>(hpar->rHB[i]/CLHEP::cm,
                           hpar->drHB[i]/CLHEP::cm));
  }
#ifdef DebugLog
  std::cout << "HB with " << nModule[0] << " modules and " << nHalves[0]
        <<" halves and " << gconsHB.size() << " layers" << std::endl;
  for (unsigned int i=0; i<gconsHB.size(); ++i) 
    std::cout << "rHB[" << i << "] = " << gconsHB[i].first << " +- "
          << gconsHB[i].second << std::endl; 
#endif
  nModule[1] = hpar->modHE[0];
  nHalves[1] = hpar->modHE[1];
  for (unsigned int i=0; i<hpar->zHE.size(); ++i) {
    gconsHE.push_back(std::pair<double,double>(hpar->zHE[i]/CLHEP::cm,
                           hpar->dzHE[i]/CLHEP::cm));
  }
#ifdef DebugLog
  std::cout << "HE with " << nModule[1] << " modules and " << nHalves[1] 
        <<" halves and " << gconsHE.size() << " layers" << std::endl;
  for (unsigned int i=0; i<gconsHE.size(); ++i) 
    std::cout << "zHE[" << i << "] = " << gconsHE[i].first << " +- "
          << gconsHE[i].second << std::endl; 
#endif
}

unsigned int HcalDDDRecConstants::layerGroup ( unsigned int  eta, unsigned int  i  ) const

private

Definition at line 557 of file HcalDDDRecConstants.cc.

References hpar, relval_2017::k, and HcalParameters::layerGroupEtaRec.

Referenced by getEtaBins(), getHCID(), getMaxDepth(), getRZ(), and initialize().

                                               {
  unsigned int k = 0;
  for( auto const & it :  hpar->layerGroupEtaRec ) {
    if( it.layer == eta + 1 )  {
      return it.layerGroup.at( i );
    }
    if( it.layer > eta + 1 ) break;
    k = it.layerGroup.at( i );
  }
  return k;
}

unsigned int HcalDDDRecConstants::layerGroupSize ( unsigned int  eta ) const

private

Definition at line 545 of file HcalDDDRecConstants.cc.

References hpar, relval_2017::k, and HcalParameters::layerGroupEtaRec.

Referenced by getEtaBins(), getMaxDepth(), getRZ(), and initialize().

                                                                         {
  unsigned int k = 0;
  for( auto const & it : hpar->layerGroupEtaRec ) {
    if( it.layer == eta + 1 ) {
      return it.layerGroup.size();
    }
    if( it.layer > eta + 1 ) break;
    k = it.layerGroup.size();
  }
  return k;
}

int HcalDDDRecConstants::maxHFDepth ( int  ieta, int  iphi  ) const

inline

Definition at line 89 of file HcalDDDRecConstants.h.

References hcons, and HcalDDDSimConstants::maxHFDepth().

Referenced by HcalTopology::validRaw().

{return hcons.maxHFDepth(ieta,iphi);}

unsigned int HcalDDDRecConstants::nCells

(

HcalSubdetector

subdet

)

const

Definition at line 394 of file HcalDDDRecConstants.cc.

References getEtaBins(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, i, kHOSizePreLS1, HLT_FULL_cff::nPhi, and HcalDDDSimConstants::numberOfCells().

                                                                     {

  if (subdet == HcalBarrel || subdet == HcalEndcap) {
    int isub   = (subdet == HcalBarrel) ? 0 : 1;
    std::vector<HcalDDDRecConstants::HcalEtaBin> etabins = getEtaBins(isub);
    unsigned int ncell(0);
    for (unsigned int i=0; i<etabins.size(); ++i) {
      ncell += (((unsigned int)(etabins[i].nPhi))*(etabins[i].layer.size()));
    }
    return ncell;
  } else if (subdet == HcalOuter) {
    return kHOSizePreLS1;
  } else if (subdet == HcalForward) {
    return (unsigned int)(hcons.numberOfCells(subdet));
  } else {
    return 0;
  }
}

unsigned int HcalDDDRecConstants::nCells

(

)

const

Definition at line 413 of file HcalDDDRecConstants.cc.

References HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

                                               {
  return (nCells(HcalBarrel)+nCells(HcalEndcap)+nCells(HcalOuter)+nCells(HcalForward));
}

unsigned int HcalDDDRecConstants::numberOfCells

(

HcalSubdetector

subdet

)

const

Definition at line 372 of file HcalDDDRecConstants.cc.

References HcalBarrel, HcalCellTypes(), HcalEndcap, hcons, i, nHalves, ecaldqm::binning::nPhiBins, pileupDistInMC::num, and HcalDDDSimConstants::numberOfCells().

Referenced by HcalDDDGeometryLoader::load().

                                                                            {

  if (subdet == HcalBarrel || subdet == HcalEndcap) {
    unsigned int num = 0;
    std::vector<HcalCellType> cellTypes = HcalCellTypes(subdet);
    for (unsigned int i=0; i<cellTypes.size(); i++) {
      num += (unsigned int)(cellTypes[i].nPhiBins());
      if (cellTypes[i].nHalves() > 1) 
    num += (unsigned int)(cellTypes[i].nPhiBins());
      num -= (unsigned int)(cellTypes[i].nPhiMissingBins());
    }
#ifdef DebugLog
    edm::LogInfo ("HCalGeom") << "HcalDDDRecConstants:numberOfCells " 
                  << cellTypes.size()  << " " << num 
                  << " for subdetector " << subdet;
#endif
    return num;
  } else {
    return hcons.numberOfCells(subdet);
  }
}

Member Data Documentation

std::vector<std::pair<int,int> > HcalDDDRecConstants::etaSimValu

private

Definition at line 101 of file HcalDDDRecConstants.h.

Referenced by HcalCellTypes(), and initialize().

std::vector<double> HcalDDDRecConstants::etaTable

private

Definition at line 102 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), getEtaLimit(), getEtaPhi(), getEtaTable(), and initialize().

std::vector<std::pair<double,double> > HcalDDDRecConstants::gconsHB

private

Definition at line 108 of file HcalDDDRecConstants.h.

Referenced by getConstHBHE(), getRZ(), getThickActive(), and initialize().

std::vector<std::pair<double,double> > HcalDDDRecConstants::gconsHE

private

Definition at line 109 of file HcalDDDRecConstants.h.

Referenced by getConstHBHE(), getRZ(), getThickActive(), and initialize().

const HcalDDDSimConstants& HcalDDDRecConstants::hcons

private

Definition at line 100 of file HcalDDDRecConstants.h.

Referenced by getDepthEta16(), getEtaBins(), getHCID(), getHFCellParameters(), getMaxDepth(), HcalCellTypes(), initialize(), maxHFDepth(), nCells(), and numberOfCells().

const HcalParameters* HcalDDDRecConstants::hpar

private

Definition at line 99 of file HcalDDDRecConstants.h.

Referenced by getDepth(), getEtaBins(), getEtaPhi(), getEtaTableHF(), getHCID(), getMaxDepth(), getNEta(), getNoff(), getPhiOff(), getPhiOffs(), getPhiTableHF(), getTopoMode(), getTriggerMode(), HcalCellTypes(), initialize(), layerGroup(), and layerGroupSize().

std::vector<int> HcalDDDRecConstants::ietaMap

private

Definition at line 103 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().

std::vector<int> HcalDDDRecConstants::iEtaMax

private

Definition at line 104 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), getEtaPhi(), getEtaRange(), getMaxDepth(), and initialize().

std::vector<int> HcalDDDRecConstants::iEtaMin

private

Definition at line 104 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), getEtaPhi(), getEtaRange(), getHCID(), getHFCellParameters(), and initialize().

std::vector<int> HcalDDDRecConstants::maxDepth

private

Definition at line 105 of file HcalDDDRecConstants.h.

Referenced by getHFCellParameters(), getMaxDepth(), and initialize().

int HcalDDDRecConstants::nHalves[2]

private

Definition at line 110 of file HcalDDDRecConstants.h.

Referenced by initialize(), and numberOfCells().

int HcalDDDRecConstants::nModule[2]

private

Definition at line 110 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), and initialize().

std::vector<double> HcalDDDRecConstants::phibin

private

Definition at line 106 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), getEtaPhi(), getHCID(), getPhiBin(), getPhiTable(), and initialize().

std::vector<int> HcalDDDRecConstants::phiUnitS

private

Definition at line 107 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().