HcalDDDSimConstants Class Reference

#include <HcalDDDSimConstants.h>

Public Member Functions
HcalCellType::HcalCell	cell (int det, int zside, int depth, int etaR, int iphi) 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 (int i) const
int	getDepthEta29 (int i) const
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
const std::vector< double > &	getEtaTableHF () const
const std::vector< double > &	getGparHF () const
const std::vector< HcalDetId > &	getIdHF2QIE () const
std::pair< int, int >	getiEtaRange (const int i) const
const std::vector< double > &	getLayer0Wt () const
int	getMaxDepth (const int type) const
std::pair< int, int >	getModHalfHBHE (const int type) const
std::pair< double, double >	getPhiCons (int det, int ieta)
const std::vector< double > &	getPhiTableHF () const
const std::vector< double > &	getRTableHF () const
std::vector< HcalCellType >	HcalCellTypes () const
std::vector< HcalCellType >	HcalCellTypes (HcalSubdetector, int ieta=-1, int depth=-1) const
	HcalDDDSimConstants (const HcalParameters *hp)
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
	~HcalDDDSimConstants ()

Private Member Functions
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
double	getGain (HcalSubdetector subdet, int depth) const
int	getShift (HcalSubdetector subdet, int depth) const
void	initialize ()
unsigned int	layerGroup (unsigned int eta, unsigned int i) const
unsigned int	layerGroupSize (unsigned int eta) const
void	printTileHB (int eta, int depth) const
void	printTileHE (int eta, int depth) const

Private Attributes
int	depthEta16 [2]
int	depthEta29 [2]
std::vector< int >	depths [nDepthMax]
double	dlShort
double	dzVcal
double	etaHO [4]
const HcalParameters *	hpar
std::vector< HcalDetId >	idHF2QIE
std::vector< int >	maxDepth
int	nDepth
int	nEta
int	nmodHB
int	nmodHE
int	nPhiF
int	nR
int	nzHB
int	nzHE
double	rminHO
double	zVcal

Static Private Attributes
static const int	nDepthMax =9

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 23 of file HcalDDDSimConstants.h.

Constructor & Destructor Documentation

HcalDDDSimConstants::HcalDDDSimConstants

(

const HcalParameters *

hp

)

Definition at line 11 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), and initialize().

                                                                 : hpar(hp) {

#ifdef DebugLog
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::HcalDDDSimConstants (const HcalParameter* hp) constructor";
#endif

  initialize();
#ifdef DebugLog
  std::vector<HcalCellType> cellTypes = HcalCellTypes();
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size()
                << " cells of type HCal (All)";
#endif
}

HcalDDDSimConstants::~HcalDDDSimConstants

(

)

Definition at line 26 of file HcalDDDSimConstants.cc.

                                          { 
#ifdef DebugLog
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants::destructed!!!";
#endif
}

Member Function Documentation

HcalCellType::HcalCell HcalDDDSimConstants::cell	(	int	det,
		int	zside,
		int	depth,
		int	etaR,
		int	iphi
	)		const

Definition at line 32 of file HcalDDDSimConstants.cc.

References deltaEta(), HcalCellType::HcalCell::deta, HcalCellType::HcalCell::dphi, HcalParameters::drHB, HcalCellType::HcalCell::drz, HcalParameters::dzHE, HcalCellType::HcalCell::eta, eta, HcalParameters::etaMax, HcalParameters::etaMin, HcalCellType::HcalCell::flagrz, getEta(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, maxDepth, nEta, HcalParameters::noff, nR, HcalCellType::HcalCell::ok, convertSQLiteXML::ok, phi, HcalCellType::HcalCell::phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rHB, HcalParameters::rTable, HcalCellType::HcalCell::rz, tmp, unitPhi(), and HcalParameters::zHE.

Referenced by HcalCellTypes(), and HcalTestAnalysis::update().

                                         {

  double etaMn = hpar->etaMin[0];
  double etaMx = hpar->etaMax[0];
  if (idet==static_cast<int>(HcalEndcap)) {
    etaMn = hpar->etaMin[1]; etaMx = hpar->etaMax[1];
  } else if (idet==static_cast<int>(HcalForward)) {
    etaMn = hpar->etaMin[2]; etaMx = hpar->etaMax[2];
  }
  double eta = 0, deta = 0, phi = 0, dphi = 0, rz = 0, drz = 0;
  bool   ok = false, flagrz = true;
  if ((idet==static_cast<int>(HcalBarrel)||idet==static_cast<int>(HcalEndcap)||
       idet==static_cast<int>(HcalOuter)||idet==static_cast<int>(HcalForward))
      && etaR >=etaMn && etaR <= etaMx && depth > 0)    ok = true;
  if (idet == static_cast<int>(HcalEndcap) && depth>(int)(hpar->zHE.size()))ok=false;
  else if (idet == static_cast<int>(HcalBarrel) && depth > 17)              ok=false;
  else if (idet == static_cast<int>(HcalOuter) && depth != 4)               ok=false;
  else if (idet == static_cast<int>(HcalForward) && depth > maxDepth[2])    ok=false;
  if (ok) {
    eta  = getEta(idet, etaR, zside, depth);
    deta = deltaEta(idet, etaR, depth);
    double fibin, fioff;
    if      (idet == static_cast<int>(HcalBarrel)||
         idet == static_cast<int>(HcalOuter)) {
      fioff = hpar->phioff[0];
      fibin = hpar->phibin[etaR-1];
    } else if (idet == static_cast<int>(HcalEndcap)) {
      fioff = hpar->phioff[1];
      fibin = hpar->phibin[etaR-1];
    } else {
      fioff = hpar->phioff[2];
      fibin = hpar->phitable[etaR-hpar->etaMin[2]];
      if (unitPhi(fibin) > 2) fioff = hpar->phioff[4];
    }
    phi  = fioff + (iphi - 0.5)*fibin;
    dphi = 0.5*fibin;
    if (idet == static_cast<int>(HcalForward)) {
      int ir = nR + hpar->etaMin[2] - etaR - 1;
      if (ir > 0 && ir < nR) {
    rz     = 0.5*(hpar->rTable[ir]+hpar->rTable[ir-1]);
    drz    = 0.5*(hpar->rTable[ir]-hpar->rTable[ir-1]);
      } else {
    ok     = false;
#ifdef DebugLog
    edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: wrong eta " << etaR 
                 << " ("  << ir << "/" << nR << ") Detector "
                 << idet;
#endif
      }
    } else if (etaR <= nEta) {
      int laymin(depth), laymax(depth);
      if (idet == static_cast<int>(HcalOuter)) {
    laymin = (etaR > hpar->noff[2]) ? ((int)(hpar->zHE.size())) : ((int)(hpar->zHE.size()))-1;
    laymax = ((int)(hpar->zHE.size()));
      }
      double d1=0, d2=0;
      if (idet == static_cast<int>(HcalEndcap)) {
    flagrz = false;
    d1     = hpar->zHE[laymin-1] - hpar->dzHE[laymin-1];
    d2     = hpar->zHE[laymax-1] + hpar->dzHE[laymax-1];
      } else {
    d1     = hpar->rHB[laymin-1] - hpar->drHB[laymin-1];
    d2     = hpar->rHB[laymax-1] + hpar->drHB[laymax-1];
      }
      rz     = 0.5*(d2+d1);
      drz    = 0.5*(d2-d1);
    } else {
      ok = false;
#ifdef DebugLog
      edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: wrong depth " << depth
                   << " or etaR " << etaR << " for detector " 
                   << idet;
#endif
    }
  } else {
    ok = false;
#ifdef DebugLog
    edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: wrong depth " << depth
                 << " det " << idet;
#endif
  }
  HcalCellType::HcalCell tmp(ok,eta,deta,phi,dphi,rz,drz,flagrz);

#ifdef DebugLog
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: det/side/depth/etaR/phi "
               << idet  << "/" << zside << "/" << depth << "/" 
               << etaR << "/" << iphi << " Cell Flag " << tmp.ok 
               << " "  << tmp.eta << " " << tmp.deta << " phi " 
               << tmp.phi << " " << tmp.dphi << " r(z) " << tmp.rz
               << " "  << tmp.drz << " " << tmp.flagrz;
#endif
  return tmp;
}

double HcalDDDSimConstants::deltaEta ( int  det, int  eta, int  depth  ) const

private

Definition at line 627 of file HcalDDDSimConstants.cc.

References dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, tmp, z, and zVcal.

Referenced by cell().

                                                                       {

  double tmp = 0;
  if (det == static_cast<int>(HcalForward)) {
    int ir = nR + hpar->etaMin[2] - etaR - 1;
    if (ir > 0 && ir < nR) {
      double z = zVcal;
      if (depth%2 != 1) z += dlShort;
      tmp = 0.5*(getEta(hpar->rTable[ir-1],z)-getEta(hpar->rTable[ir],z));
    }
  } else {
    if (etaR > 0 && etaR < nEta) {
      if (etaR == hpar->noff[1]-1 && depth > 2) {
    tmp = 0.5*(hpar->etaTable[etaR+1]-hpar->etaTable[etaR-1]);
      } else if (det == static_cast<int>(HcalOuter)) {
    if (etaR == hpar->noff[2]) {
      tmp = 0.5*(etaHO[0]-hpar->etaTable[etaR-1]);
    } else if (etaR == hpar->noff[2]+1) {
      tmp = 0.5*(hpar->etaTable[etaR]-etaHO[1]);
    } else if (etaR == hpar->noff[3]) {
      tmp = 0.5*(etaHO[2]-hpar->etaTable[etaR-1]);
    } else if (etaR == hpar->noff[3]+1) {
      tmp = 0.5*(hpar->etaTable[etaR]-etaHO[3]);
    } else {
      tmp = 0.5*(hpar->etaTable[etaR]-hpar->etaTable[etaR-1]);
    }
      } else {
    tmp = 0.5*(hpar->etaTable[etaR]-hpar->etaTable[etaR-1]);
      }
    } 
  }
#ifdef DebugLog
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::deltaEta " << etaR << " " 
               << depth << " ==> " << tmp;
#endif
  return tmp;
}

unsigned int HcalDDDSimConstants::findLayer	(	int	layer,
		const std::vector< HcalParameters::LayerItem > &	layerGroup
	)		const

Definition at line 247 of file HcalDDDSimConstants.cc.

References i.

                                                                                                                 {
  
  unsigned int id = layerGroup.size();
  for (unsigned int i = 0; i < layerGroup.size(); i++) {
    if (layer == (int)(layerGroup[i].layer)) {
      id = i;
      break;
    }
  }
  return id;
}

std::vector< std::pair< double, double > > HcalDDDSimConstants::getConstHBHE

(

const int

type

)

const

Definition at line 128 of file HcalDDDSimConstants.cc.

References HcalParameters::drHB, HcalParameters::dzHE, hpar, i, HcalParameters::rHB, and HcalParameters::zHE.

                                                                                         {

  std::vector<std::pair<double,double> > gcons;
  if (type == 0) {
    for (unsigned int i=0; i<hpar->rHB.size(); ++i) {
      gcons.push_back(std::pair<double,double>(hpar->rHB[i],hpar->drHB[i]));
    }
  } else {
    for (unsigned int i=0; i<hpar->zHE.size(); ++i) {
      gcons.push_back(std::pair<double,double>(hpar->zHE[i],hpar->dzHE[i]));
    }
  }
  return gcons;
}

int HcalDDDSimConstants::getDepthEta16 ( int  i ) const

inline

Definition at line 33 of file HcalDDDSimConstants.h.

References depthEta16.

Referenced by HcalDDDRecConstants::getDepthEta16(), HcalDDDRecConstants::getEtaBins(), HcalDDDRecConstants::getHCID(), HcalDDDRecConstants::getMaxDepth(), and HcalDDDRecConstants::initialize().

{return ((i==1) ? depthEta16[1] : depthEta16[0]);}

int HcalDDDSimConstants::getDepthEta29 ( int  i ) const

inline

Definition at line 34 of file HcalDDDSimConstants.h.

References depthEta29.

Referenced by HcalDDDRecConstants::getEtaBins(), and HcalDDDRecConstants::getMaxDepth().

{return ((i==1) ? depthEta29[1] : depthEta29[0]);}

std::pair< int, double > HcalDDDSimConstants::getDetEta	(	double	eta,
		int	depth
	)

Definition at line 144 of file HcalDDDSimConstants.cc.

References eta, HcalParameters::etaMax, HLT_FULL_cff::etaMin, MuonErrorMatrixValues_cff::etaRange, HcalParameters::etaTable, create_public_lumi_plots::exp, HcalBarrel, HcalEndcap, HcalForward, hpar, i, nEta, funct::tan(), theta(), and zVcal.

Referenced by HcalNumberingFromDDD::unitID().

                                                                        {

  int    hsubdet(0), ieta(0);
  double etaR(0);
  double heta = fabs(eta);
  for (int i = 0; i < nEta; i++)
    if (heta > hpar->etaTable[i]) ieta = i + 1;
  if (heta <= hpar->etaRange[1]) {
    if ((ieta <= hpar->etaMin[1] && depth==3) || ieta > hpar->etaMax[0]) {
      hsubdet = static_cast<int>(HcalEndcap);
    } else {
      hsubdet = static_cast<int>(HcalBarrel);
    }
    etaR    = eta;
  } else {
    hsubdet = static_cast<int>(HcalForward);
    double theta = 2.*atan(exp(-heta));
    double hR    = zVcal*tan(theta);
    etaR    = (eta >= 0. ? hR : -hR);
  }
  return std::pair<int,double>(hsubdet,etaR);
}

int HcalDDDSimConstants::getEta	(	int	det,
		int	lay,
		double	hetaR
	)

Definition at line 167 of file HcalDDDSimConstants.cc.

References etaHO, HcalParameters::etaMax, HcalParameters::etaMin, HLT_FULL_cff::etaMin, HcalParameters::etaTable, HcalBarrel, HcalEndcap, HcalForward, hpar, i, nEta, HcalParameters::noff, and nR.

Referenced by cell(), deltaEta(), getEta(), and HcalNumberingFromDDD::unitID().

                                                              {

  int    ieta(0);
  if (det == static_cast<int>(HcalForward)) { // Forward HCal
    ieta    = hpar->etaMax[2];
    for (int i = nR-1; i > 0; i--)
      if (hetaR < hpar->rTable[i]) ieta = hpar->etaMin[2] + nR - i - 1;
  } else { // Barrel or Endcap
    ieta  = 1;
    for (int i = 0; i < nEta-1; i++)
      if (hetaR > hpar->etaTable[i]) ieta = i + 1;
    if (det == static_cast<int>(HcalBarrel)) {
      if (ieta > hpar->etaMax[0])  ieta = hpar->etaMax[0];
      if (lay == 18) {
    if (hetaR > etaHO[1] && ieta == hpar->noff[2]) ieta++;
      }
    } else if (det == static_cast<int>(HcalEndcap)) {
      if (ieta <= hpar->etaMin[1]) ieta = hpar->etaMin[1];
    }
  }
  return ieta;
}

double HcalDDDSimConstants::getEta ( int  det, int  etaR, int  zside, int  depth = 1  ) const

private

Definition at line 665 of file HcalDDDSimConstants.cc.

References dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, tmp, z, and zVcal.

                                    {

  double tmp = 0;
  if (det == static_cast<int>(HcalForward)) {
    int ir = nR + hpar->etaMin[2] - etaR - 1;
    if (ir > 0 && ir < nR) {
      double z = zVcal;
      if (depth%2 != 1) z += dlShort;
      tmp = 0.5*(getEta(hpar->rTable[ir-1],z)+getEta(hpar->rTable[ir],z));
    }
  } else {
    if (etaR > 0 && etaR < nEta) {
      if (etaR == hpar->noff[1]-1 && depth > 2) {
    tmp = 0.5*(hpar->etaTable[etaR+1]+hpar->etaTable[etaR-1]);
      } else if (det == static_cast<int>(HcalOuter)) {
    if (etaR == hpar->noff[2]) {
      tmp = 0.5*(etaHO[0]+hpar->etaTable[etaR-1]);
    } else if (etaR == hpar->noff[2]+1) {
      tmp = 0.5*(hpar->etaTable[etaR]+etaHO[1]);
    } else if (etaR == hpar->noff[3]) {
      tmp = 0.5*(etaHO[2]+hpar->etaTable[etaR-1]);
    } else if (etaR == hpar->noff[3]+1) {
      tmp = 0.5*(hpar->etaTable[etaR]+etaHO[3]);
    } else {
      tmp = 0.5*(hpar->etaTable[etaR]+hpar->etaTable[etaR-1]);
    }
      } else {
    tmp = 0.5*(hpar->etaTable[etaR]+hpar->etaTable[etaR-1]);
      }
    }
  } 
  if (zside == 0) tmp = -tmp;
#ifdef DebugLog
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::getEta " << etaR << " " 
               << zside << " " << depth << " ==> " << tmp;
#endif
  return tmp;
}

double HcalDDDSimConstants::getEta ( double  r, double  z  ) const

private

Definition at line 705 of file HcalDDDSimConstants.cc.

References dqm-mbProfile::log, funct::tan(), and tmp.

                                                           {

  double tmp = 0;
  if (z != 0) tmp = -log(tan(0.5*atan(r/z)));
#ifdef DebugLog
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::getEta " << r << " " << z 
               << " ==> " << tmp;
#endif
  return tmp;
}

std::pair< int, int > HcalDDDSimConstants::getEtaDepth	(	int	det,
		int	etaR,
		int	phi,
		int	depth,
		int	lay
	)

Definition at line 190 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, depthEta16, depthEta29, HcalParameters::etaMin, HcalBarrel, HcalForward, HcalOuter, hpar, layerGroup(), maxDepth, HcalParameters::noff, HcalParameters::phibin, and HcalParameters::phioff.

Referenced by HcalNumberingFromDDD::unitID().

                                                {

  //Modify the depth index
  if (det == static_cast<int>(HcalForward)) { // Forward HCal
  } else if (det == static_cast<int>(HcalOuter)) {
    depth = 4;
  } else {
    if (lay >= 0) {
      depth= layerGroup( etaR-1, lay-1 );
      if (etaR == hpar->noff[0] && lay > 1) {
    int   kphi   = phi + int((hpar->phioff[3]+0.1)/hpar->phibin[etaR-1]);
    kphi         = (kphi-1)%4 + 1;
    if (kphi == 2 || kphi == 3) depth = layerGroup( etaR-1, lay-2 );
      }
    } else if (det == static_cast<int>(HcalBarrel)) {
      if (depth>maxDepth[0]) depth = maxDepth[0];
    }
    if (etaR >= hpar->noff[1] && depth > depthEta29[0]) {
      etaR -= depthEta29[1];
    } else if (etaR == hpar->etaMin[1]) {
      if (det == static_cast<int>(HcalBarrel)) {
    if (depth > depthEta16[0]) depth = depthEta16[0];
      } else {
    if (depth < depthEta16[1]) depth = depthEta16[1];
      }
    }
  }
  return std::pair<int,int>(etaR,depth);
}

double HcalDDDSimConstants::getEtaHO	(	double &	etaR,
		double &	x,
		double &	y,
		double &	z
	)		const

Definition at line 221 of file HcalDDDSimConstants.cc.

References eta, HcalParameters::etaTable, hpar, alignCSCRings::r, rminHO, mathSSE::sqrt(), and HcalParameters::zHO.

Referenced by HcalNumberingFromDDD::unitID().

                                      {

  if (hpar->zHO.size() > 4) {
    double eta  = fabs(etaR);
    double r    = std::sqrt(x*x+y*y);
    if (r > rminHO) {
      double zz = fabs(z);
      if (zz > hpar->zHO[3]) {
    if (eta <= hpar->etaTable[10]) eta = hpar->etaTable[10]+0.001;
      } else if (zz > hpar->zHO[1]) {
    if (eta <= hpar->etaTable[4])  eta = hpar->etaTable[4]+0.001;
      }
    }
    eta = (z >= 0. ? eta : -eta);
#ifdef DebugLog
    edm::LogInfo ("HCalGeom") << "R " << r << " Z " << z << " eta " << etaR 
                  << ":" << eta;
    if (eta != etaR) edm::LogInfo ("HCalGeom") << "**** Check *****";
#endif
    return eta;
  } else {
    return etaR;
  }
}

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

inline

Definition at line 42 of file HcalDDDSimConstants.h.

References HcalParameters::etaTableHF, and hpar.

{return hpar->etaTableHF;}

double HcalDDDSimConstants::getGain ( HcalSubdetector  subdet, int  depth  ) const

private

Definition at line 736 of file HcalDDDSimConstants.cc.

References HcalParameters::HBGains, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEGains, HcalParameters::HFGains, and hpar.

Referenced by HcalCellTypes().

                                                                           {

  double gain;
  switch(subdet) {
  case HcalEndcap:
    gain = hpar->HEGains[0];
    break;
  case HcalForward:
    gain = hpar->HFGains[(depth-1)%2];
    break;
  case HcalOuter:
    gain = hpar->HBGains[3];
    break;
  default:
    gain = hpar->HBGains[0];
    break;
  }
  return gain;
}

const std::vector<double>& HcalDDDSimConstants::getGparHF ( ) const

inline

Definition at line 44 of file HcalDDDSimConstants.h.

References HcalParameters::gparHF, and hpar.

Referenced by HFFibre::initRun(), HFShower::initRun(), HFShowerParam::initRun(), HFShowerLibrary::initRun(), and HCalSD::update().

{return hpar->gparHF;}

const std::vector<HcalDetId>& HcalDDDSimConstants::getIdHF2QIE ( ) const

inline

Definition at line 49 of file HcalDDDSimConstants.h.

References idHF2QIE.

Referenced by HcalDDDRecConstants::getHFCellParameters().

{return idHF2QIE;}

std::pair<int,int> HcalDDDSimConstants::getiEtaRange ( const int  i ) const

inline

Definition at line 41 of file HcalDDDSimConstants.h.

References HcalParameters::etaMax, HcalParameters::etaMin, hpar, and i.

{return std::pair<int,int>(hpar->etaMin[i],hpar->etaMax[i]);}

const std::vector<double>& HcalDDDSimConstants::getLayer0Wt ( ) const

inline

Definition at line 45 of file HcalDDDSimConstants.h.

References hpar, and HcalParameters::Layer0Wt.

Referenced by HCalSD::update().

{return hpar->Layer0Wt;}

int HcalDDDSimConstants::getMaxDepth ( const int  type ) const

inline

Definition at line 46 of file HcalDDDSimConstants.h.

References maxDepth.

Referenced by HCalSD::update().

{return maxDepth[type];}

std::pair< int, int > HcalDDDSimConstants::getModHalfHBHE

(

const int

type

)

const

Definition at line 259 of file HcalDDDSimConstants.cc.

References nmodHB, nmodHE, nzHB, and nzHE.

                                                                         {

  if (type == 0) {
    return std::pair<int,int>(nmodHB,nzHB);
  } else {
    return std::pair<int,int>(nmodHE,nzHE);
  }
}

std::pair< double, double > HcalDDDSimConstants::getPhiCons	(	int	det,
		int	ieta
	)

Definition at line 268 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, HcalBarrel, HcalForward, hpar, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, and unitPhi().

Referenced by HcalNumberingFromDDD::unitID().

                                                                        {

  double fioff(0), fibin(0);
  if (det == static_cast<int>(HcalForward)) { // Forward HCal
    fioff   = hpar->phioff[2];
    fibin   = hpar->phitable[ieta-hpar->etaMin[2]];
    if  (unitPhi(fibin) > 2) {   // HF double-phi  
      fioff = hpar->phioff[4];
    }
  } else { // Barrel or Endcap
    if (det == static_cast<int>(HcalBarrel)) {
      fioff   = hpar->phioff[0];
    } else {
      fioff   = hpar->phioff[1];
    }
    fibin = hpar->phibin[ieta-1];
  }
  return std::pair<double,double>(fioff,fibin);
}

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

inline

Definition at line 50 of file HcalDDDSimConstants.h.

References hpar, and HcalParameters::phitable.

Referenced by HcalDDDRecConstants::getHFCellParameters(), and HFShowerLibrary::initRun().

{return hpar->phitable;}

const std::vector<double>& HcalDDDSimConstants::getRTableHF ( ) const

inline

Definition at line 51 of file HcalDDDSimConstants.h.

References hpar, and HcalParameters::rTable.

Referenced by HcalDDDRecConstants::getHFCellParameters(), HFShowerPMT::initRun(), HFFibre::initRun(), HFShowerFibreBundle::initRun(), and HFShowerLibrary::initRun().

{return hpar->rTable;}

int HcalDDDSimConstants::getShift ( HcalSubdetector  subdet, int  depth  ) const

private

Definition at line 716 of file HcalDDDSimConstants.cc.

References HcalParameters::HBShift, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEShift, HcalParameters::HFShift, hpar, and edm::shift.

Referenced by HcalCellTypes().

                                                                         {

  int shift;
  switch(subdet) {
  case HcalEndcap:
    shift = hpar->HEShift[0];
    break;
  case HcalForward:
    shift = hpar->HFShift[(depth-1)%2];
    break;
  case HcalOuter:
    shift = hpar->HBShift[3];
    break;
  default:
    shift = hpar->HBShift[0];
    break;
  }
  return shift;
}

std::vector< HcalCellType > HcalDDDSimConstants::HcalCellTypes

(

)

const

Definition at line 288 of file HcalDDDSimConstants.cc.

References HcalBarrel, HcalEndcap, HcalForward, HcalOuter, and i.

Referenced by HcalDDDRecConstants::HcalCellTypes(), HcalDDDSimConstants(), and numberOfCells().

                                                                {

  std::vector<HcalCellType> cellTypes = HcalCellTypes(HcalBarrel);
#ifdef DebugLog
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size()
            << " cells of type HCal Barrel";
  for (unsigned int i=0; i<cellTypes.size(); i++)
    edm::LogInfo ("HCalGeom") << "Cell " << i << " " << cellTypes[i];
#endif

  std::vector<HcalCellType> hoCells   = HcalCellTypes(HcalOuter);
#ifdef DebugLog
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants: " << hoCells.size()
            << " cells of type HCal Outer";
  for (unsigned int i=0; i<hoCells.size(); i++)
    edm::LogInfo ("HCalGeom") << "Cell " << i << " " << hoCells[i];
#endif
  cellTypes.insert(cellTypes.end(), hoCells.begin(), hoCells.end());

  std::vector<HcalCellType> heCells   = HcalCellTypes(HcalEndcap);
#ifdef DebugLog
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants: " << heCells.size()
            << " cells of type HCal Endcap";
  for (unsigned int i=0; i<heCells.size(); i++)
    edm::LogInfo ("HCalGeom") << "Cell " << i << " " << heCells[i];
#endif
  cellTypes.insert(cellTypes.end(), heCells.begin(), heCells.end());

  std::vector<HcalCellType> hfCells   = HcalCellTypes(HcalForward);
#ifdef DebugLog
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants: " << hfCells.size()
            << " cells of type HCal Forward";
  for (unsigned int i=0; i<hfCells.size(); i++)
    edm::LogInfo ("HCalGeom") << "Cell " << i << " " << hfCells[i];
#endif
  cellTypes.insert(cellTypes.end(), hfCells.begin(), hfCells.end());

  return cellTypes;
}

std::vector< HcalCellType > HcalDDDSimConstants::HcalCellTypes	(	HcalSubdetector	subdet,
		int	ieta = -1,
		int	depth = -1
	)		const

Definition at line 328 of file HcalDDDSimConstants.cc.

References cell(), HLT_FULL_cff::depth, dlShort, dzVcal, eta, HcalParameters::etaMax, HcalParameters::etaMin, getGain(), getShift(), HcalEndcap, HcalForward, HcalOuter, hpar, idHF2QIE, relval_2017::k, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::l, maxDepth, nmodHB, nmodHE, HcalParameters::noff, HcalCellType::nPhiBins(), HcalCellType::nPhiModule(), nzHB, nzHE, HcalCellType::HcalCell::ok, phi, HcalCellType::setMissingPhi(), edm::shift, unitPhi(), HcalCellType::unitPhi(), units(), and ecaldqm::zside().

                                                             {

  std::vector<HcalCellType> cellTypes;
  if (subdet == HcalForward) {
    if (dzVcal < 0) return cellTypes;
  }

  int    dmin, dmax, indx, nz, nmod;
  double hsize = 0;
  switch(subdet) {
  case HcalEndcap:
    dmin = 1; dmax = 19; indx = 1; nz = nzHE; nmod = nmodHE;
    break;
  case HcalForward:
    dmin = 1; dmax = maxDepth[2]; indx = 2; nz = 2; nmod = 18; 
    break;
  case HcalOuter:
    dmin = 4; dmax = 4; indx = 0; nz = nzHB; nmod = nmodHB;
    break;
  default:
    dmin = 1; dmax = 17; indx = 0; nz = nzHB; nmod = nmodHB;
    break;
  }
  if (depthl > 0) dmin = dmax = depthl;
  int ietamin = (ieta>0) ? ieta : hpar->etaMin[indx];
  int ietamax = (ieta>0) ? ieta : hpar->etaMax[indx];
  int phi = 1, zside  = 1;

  // Get the Cells 
  int subdet0 = static_cast<int>(subdet);
  for (int depth=dmin; depth<=dmax; depth++) {
    int    shift = getShift(subdet, depth);
    double gain  = getGain (subdet, depth);
    if (subdet == HcalForward) {
      if (depth%2 == 1) hsize = dzVcal;
      else              hsize = dzVcal-0.5*dlShort;
    }
    for (int eta=ietamin; eta<= ietamax; eta++) {
      HcalCellType::HcalCell temp1 = cell(subdet0,zside,depth,eta,phi);
      if (temp1.ok) {
    int units = unitPhi (subdet0, eta);
    HcalCellType temp2(subdet, eta, phi, depth, temp1,
               shift, gain, nz, nmod, hsize, units);
    if (subdet == HcalOuter) {
      if (eta == hpar->noff[4]) {
        std::vector<int> missPlus, missMinus;
        int kk = 7;
        for (int miss=0; miss<hpar->noff[5]; miss++) {
          missPlus.push_back(hpar->noff[kk]);
          kk++;
        }
        for (int miss=0; miss<hpar->noff[6]; miss++) {
          missMinus.push_back(hpar->noff[kk]);
          kk++;
        }
        temp2.setMissingPhi(missPlus, missMinus);
      }
    } else if (subdet == HcalForward) {
      if (idHF2QIE.size() > 0 && depth > 2) {
        int phiInc = 4/temp2.nPhiModule();
        int iphi   = 1;
        if (temp2.unitPhi() == 4) iphi = 3;
        std::vector<int> missPlus, missMinus;
        for (int k=0; k<temp2.nPhiBins(); ++k) {
          bool okp(true), okn(true);
          for (unsigned int l=0; l<idHF2QIE.size(); ++l) {
        if ((eta == idHF2QIE[l].ietaAbs()) && 
            (iphi == idHF2QIE[l].iphi())) {
          if (idHF2QIE[l].zside() == 1) okp = false;
          else                          okn = false;
        }
          }
          if (okp) missPlus.push_back(iphi);
          if (okn) missMinus.push_back(iphi);
          iphi += phiInc;
        }
        temp2.setMissingPhi(missPlus, missMinus);
      }
    }
    cellTypes.push_back(temp2);
      }
    }
  }
  return cellTypes;
}

void HcalDDDSimConstants::initialize ( void  )

private

Definition at line 498 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, depthEta16, depthEta29, depths, dlShort, HcalParameters::dzVcal, dzVcal, etaHO, HcalParameters::etaMax, HLT_FULL_cff::etaMax, HcalParameters::etaMin, HcalParameters::etaTable, HcalParameters::gparHF, HcalForward, hpar, i, idHF2QIE, reco::if(), relval_2017::k, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::l, layerGroup(), layerGroupSize(), HcalParameters::maxDepth, maxDepth, HcalParameters::modHB, HcalParameters::modHE, nEta, nmodHB, nmodHE, HcalParameters::noff, nPhiF, nR, nzHB, nzHE, HcalParameters::rHB, HcalParameters::rHO, rminHO, HcalParameters::rTable, HcalParameters::zHO, and zVcal.

Referenced by HcalDDDSimConstants().

                                           {

  nEta      = hpar->etaTable.size();
  nR        = hpar->rTable.size();
  nPhiF     = nR - 1;

#ifdef DebugLog
  for (int i=0; i<nEta-1; ++i) {
    std::cout << "HcalDDDSimConstants:Read LayerGroup" << i << ":";
    for (unsigned int k=0; k<layerGroupSize( i ); k++) 
      std::cout << " [" << k << "] = " << layerGroup( i, k );
    std::cout << std::endl;
  }
#endif

  // Geometry parameters for HF
  dlShort   = hpar->gparHF[0];
  zVcal     = hpar->gparHF[4];
  dzVcal    = hpar->dzVcal;
#ifdef DebugLog
  std::cout << "HcalDDDSimConstants: dlShort " << dlShort << " zVcal " << zVcal
        << " and dzVcal " << dzVcal << std::endl;
#endif

  //Transform some of the parameters
  maxDepth = hpar->maxDepth;
  maxDepth[0] = maxDepth[1] = 0;
  for (int i=0; i<nEta-1; ++i) {
    unsigned int imx = layerGroupSize( i );
    int laymax = (imx > 0) ? layerGroup( i, imx-1 ) : 0;
    if (i < hpar->etaMax[0]) {
      int laymax0 = (imx > 16) ? layerGroup( i, 16 ) : laymax;
      if (i+1 == hpar->etaMax[0] && laymax0 > 2) laymax0 = 2;
      if (maxDepth[0] < laymax0) maxDepth[0] = laymax0;
    }
    if (i >= hpar->etaMin[1]-1 && i < hpar->etaMax[1]) {
      if (maxDepth[1] < laymax) maxDepth[1] = laymax;
    }
  }
#ifdef DebugLog
  for (int i=0; i<4; ++i)
    std::cout << "Detector Type [" << i << "] iEta " << hpar->etaMin[i] << ":" 
          << hpar->etaMax[i] << " MaxDepth " << maxDepth[i] << std::endl;
#endif

  int maxdepth = (maxDepth[1]>maxDepth[0]) ? maxDepth[1] : maxDepth[0];
  for (int i=0; i<maxdepth; ++i) {
    for (int k=0; k<nEta-1; ++k) {
      int layermx = ((k+1 < hpar->etaMin[1]) && i < maxDepth[0]) ? 17 : (int)layerGroupSize( k );
      int ll      = layermx;
      for (int l=layermx-1; l >= 0; --l) {
    if ((int)layerGroup( k, l ) == i+1) {
      ll = l+1; break;
    }
      }
      depths[i].push_back(ll);
    }

#ifdef DebugLog
    std::cout << "Depth " << i << " with " << depths[i].size() << " etas:";
    for (int k=0; k<nEta-1; ++k) std::cout << " " << depths[i][k];
    std::cout << std::endl;
#endif
  }

  nzHB   = hpar->modHB[0];
  nmodHB = hpar->modHB[1];
  nzHE   = hpar->modHE[0];
  nmodHE = hpar->modHE[1];
#ifdef DebugLog
  std::cout << "HcalDDDSimConstants:: " << nzHB << ":" << nmodHB
        << " barrel and " << nzHE << ":" << nmodHE
        << " endcap half-sectors" << std::endl;
#endif

  if (hpar->rHB.size() > 17 && hpar->zHO.size() > 4) {
    rminHO   = hpar->rHO[0];
    for (int k=0; k<4; ++k) etaHO[k] = hpar->rHO[k+1];
  } else {
    rminHO   =-1.0;
    etaHO[0] = hpar->etaTable[4];
    etaHO[1] = hpar->etaTable[4];
    etaHO[2] = hpar->etaTable[10];
    etaHO[3] = hpar->etaTable[10];
  }
#ifdef DebugLog
  std::cout << "HO Eta boundaries " << etaHO[0] << " " << etaHO[1]
        << " " << etaHO[2] << " " << etaHO[3] << std::endl;
  std::cout << "HO Parameters " << rminHO << " " << hpar->zHO.size();
  for (int i=0; i<4; ++i) std::cout << " eta[" << i << "] = " << etaHO[i];
  for (unsigned int i=0; i<hpar->zHO.size(); ++i) 
    std::cout << " zHO[" << i << "] = " << hpar->zHO[i];
  std::cout << std::endl;
#endif

  int noffsize = 7 + hpar->noff[5] + hpar->noff[6];
  if ((int)(hpar->noff.size()) > (noffsize+3)) {
    depthEta16[0] = hpar->noff[noffsize];
    depthEta16[1] = hpar->noff[noffsize+1];
    depthEta29[0] = hpar->noff[noffsize+2];
    depthEta29[1] = hpar->noff[noffsize+3];
  } else {
    depthEta16[0] = 2;
    depthEta16[1] = 3;
    depthEta29[0] = 2;
    depthEta29[1] = 1;
  }
#ifdef DebugLog
  std::cout << "Depth index at ieta = 16 for HB (max) " << depthEta16[0] 
        << " HE (min) " << depthEta16[1] << "; max depth for itea = 29 : ("
        << depthEta29[0] << ":" << depthEta29[1] << ")" << std::endl;
#endif
  
  if ((int)(hpar->noff.size()) > (noffsize+4)) {
    int npair = hpar->noff[noffsize+4];
    int kk    = noffsize+4;
    for (int k=0; k<npair; ++k) {
      idHF2QIE.push_back(HcalDetId(HcalForward,hpar->noff[kk+1],hpar->noff[kk+2],1));
      kk += 2;
    }
  }
#ifdef DebugLog
  std::cout << idHF2QIE.size() << " detector channels having 2 QIE cards:" 
        << std::endl;
  for (unsigned int k=0; k<idHF2QIE.size(); ++k) 
    std::cout << " [" << k << "] " << idHF2QIE[k] << std::endl;
#endif
}

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

private

Definition at line 864 of file HcalDDDSimConstants.cc.

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

Referenced by getEtaDepth(), and initialize().

                                               {
  unsigned int k = 0;
  for( auto const & it :  hpar->layerGroupEtaSim ) {
    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 HcalDDDSimConstants::layerGroupSize ( unsigned int  eta ) const

private

Definition at line 852 of file HcalDDDSimConstants.cc.

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

Referenced by initialize().

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

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

Definition at line 415 of file HcalDDDSimConstants.cc.

References idHF2QIE, relval_2017::k, maxDepth, and convertSQLiteXML::ok.

Referenced by HcalDDDRecConstants::maxHFDepth(), FastHFShowerLibrary::modifyDepth(), and HCalSD::modifyDepth().

                                                           {
  int mxdepth = maxDepth[2];
  if (idHF2QIE.size() > 0) {
    bool ok(false);
    for (unsigned int k=0; k<idHF2QIE.size(); ++k) {
      if ((eta == idHF2QIE[k].ieta()) && (iphi == idHF2QIE[k].iphi())) {
    ok = true; break;
      }
    }
    if (!ok) mxdepth = 2;
  }
  return mxdepth;
}

unsigned int HcalDDDSimConstants::numberOfCells

(

HcalSubdetector

subdet

)

const

Definition at line 429 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), i, ecaldqm::binning::nPhiBins, and pileupDistInMC::num.

Referenced by HcalDDDRecConstants::nCells(), and HcalDDDRecConstants::numberOfCells().

                                                                           {

  unsigned int num = 0;
  std::vector<HcalCellType> cellTypes = HcalCellTypes(subdet);
  for (unsigned int i=0; i<cellTypes.size(); i++) {
    int ncur = (cellTypes[i].nHalves() > 1) ? 2*cellTypes[i].nPhiBins() :
      cellTypes[i].nPhiBins();
    num += (unsigned int)(ncur);
    num -= (unsigned int)(cellTypes[i].nPhiMissingBins());
  }
#ifdef DebugLog
  edm::LogInfo ("HCalGeom") << "HcalDDDSimConstants:numberOfCells " 
            << cellTypes.size()  << " " << num 
            << " for subdetector " << subdet;
#endif
  return num;
}

int HcalDDDSimConstants::phiNumber	(	int	phi,
		int	unit
	)		const

Definition at line 447 of file HcalDDDSimConstants.cc.

References phi.

Referenced by HcalDDDRecConstants::getHCID(), and HcalNumberingFromDDD::unitID().

                                                           {

  int iphi_skip = phi;
  if      (units==2) iphi_skip  = (phi-1)*2+1;
  else if (units==4) iphi_skip  = (phi-1)*4-1;
  if (iphi_skip < 0) iphi_skip += 72;
  return iphi_skip;
}

void HcalDDDSimConstants::printTileHB ( int  eta, int  depth  ) const

private

Definition at line 756 of file HcalDDDSimConstants.cc.

References funct::cos(), gather_cfg::cout, depths, HcalParameters::dxHB, HcalParameters::dyHB, HcalParameters::etaTable, create_public_lumi_plots::exp, hpar, relval_2017::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHB, min(), HcalParameters::rhoxHB, funct::sin(), SiStripMonitorClusterAlca_cfi::zmax, and SiStripMonitorClusterAlca_cfi::zmin.

Referenced by printTiles().

                                                              {
  std::cout << "HcalDDDSimConstants::printTileHB for eta " << eta << " and depth " << depth << "\n";
  
  double etaL   = hpar->etaTable.at(eta-1);
  double thetaL = 2.*atan(exp(-etaL));
  double etaH   = hpar->etaTable.at(eta);
  double thetaH = 2.*atan(exp(-etaH));
  int    layL=0, layH=0;
  if (depth == 1) {
    layH = depths[0][eta-1];
  } else {
    layL = depths[0][eta-1];
    layH = depths[1][eta-1];
  }
  std::cout << "\ntileHB:: eta|depth " << eta << "|" << depth << " theta " << thetaH/CLHEP::deg << ":" << thetaL/CLHEP::deg << " Layer " << layL << ":" << layH-1 << "\n";
  for (int lay=layL; lay<layH; ++lay) {
    std::vector<double> area(2,0);
    int kk=0;
    for (unsigned int k=0; k<hpar->layHB.size(); ++k) {
      if (lay == hpar->layHB[k]) {
    double zmin = hpar->rhoxHB[k]*std::cos(thetaL)/std::sin(thetaL);
    double zmax = hpar->rhoxHB[k]*std::cos(thetaH)/std::sin(thetaH);
    double dz   = (std::min(zmax,hpar->dxHB[k]) - zmin);
    if (dz > 0) {
      area[kk] = dz*hpar->dyHB[k];
      kk++;
    }
      }
    }
    if (area[0] > 0) std::cout << std::setw(2) << lay << " Area " << std::setw(8) << area[0] << " " << std::setw(8) << area[1] << "\n";
  }
}

void HcalDDDSimConstants::printTileHE ( int  eta, int  depth  ) const

private

Definition at line 789 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, depths, HcalParameters::dx1HE, HcalParameters::dyHE, eta, HcalParameters::etaTable, create_public_lumi_plots::exp, hpar, relval_2017::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHE, bookConverter::max, min(), HcalParameters::phibin, HcalParameters::rhoxHE, funct::tan(), and HcalParameters::zxHE.

Referenced by printTiles().

                                                              {

  double etaL   = hpar->etaTable[eta-1];
  double thetaL = 2.*atan(exp(-etaL));
  double etaH   = hpar->etaTable[eta];
  double thetaH = 2.*atan(exp(-etaH));
  int    layL=0, layH=0;
  if (eta == 16) {
    layH = depths[2][eta-1];
  } else if (depth == 1) {
    layH = depths[0][eta-1];
  } else if (depth == 2) {
    layL = depths[0][eta-1];
    layH = depths[1][eta-1];
  } else {
    layL = depths[1][eta-1];
    layH = depths[2][eta-1];
  }
  double phib  = hpar->phibin[eta-1];
  int nphi = 2;
  if (phib > 6*CLHEP::deg) nphi = 1;
  std::cout << "\ntileHE:: Eta/depth " << eta << "|" << depth << " theta " << thetaH/CLHEP::deg << ":" << thetaL/CLHEP::deg << " Layer " << layL << ":" << layH-1 << " phi " << nphi << "\n";
  for (int lay=layL; lay<layH; ++lay) {
    std::vector<double> area(4,0);
    int kk=0;
    for (unsigned int k=0; k<hpar->layHE.size(); ++k) {
      if (lay == hpar->layHE[k]) {
    double rmin = hpar->zxHE[k]*std::tan(thetaH);
    double rmax = hpar->zxHE[k]*std::tan(thetaL);
    if ((lay != 0 || eta == 18) && 
        (lay != 1 || (eta == 18 && hpar->rhoxHE[k]-hpar->dyHE[k] > 1000) ||
         (eta != 18 && hpar->rhoxHE[k]-hpar->dyHE[k] < 1000)) &&
        rmin+30 < hpar->rhoxHE[k]+hpar->dyHE[k] && 
        rmax > hpar->rhoxHE[k]-hpar->dyHE[k]) {
      rmin = std::max(rmin,hpar->rhoxHE[k]-hpar->dyHE[k]);
      rmax = std::min(rmax,hpar->rhoxHE[k]+hpar->dyHE[k]);
      double dx1 = rmin*std::tan(phib);
      double dx2 = rmax*std::tan(phib);
      double ar1=0, ar2=0;
      if (nphi == 1) {
        ar1 = 0.5*(rmax-rmin)*(dx1+dx2-4.*hpar->dx1HE[k]);
      } else {
        ar1 = 0.5*(rmax-rmin)*(dx1+dx2-2.*hpar->dx1HE[k]);
        ar2 = 0.5*(rmax-rmin)*((rmax+rmin)*tan(10.*CLHEP::deg)-4*hpar->dx1HE[k])-ar1;
      }
      area[kk]   = ar1;
      area[kk+2] = ar2;
      kk++;
    }
      }
    }
    if (area[0] > 0 && area[1] > 0) {
      int lay0 = lay-1;
      if (eta == 18) lay0++;
      if (nphi == 1) {
    std::cout << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " " << std::setw(8) << area[1] << "\n";
      } else {
    std::cout << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " " << std::setw(8) << area[1] << ":" << std::setw(8) << area[2] << " " << std::setw(8) << area[3] << "\n";
      }
    }
  }
}

void HcalDDDSimConstants::printTiles

(

)

const

Definition at line 456 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, HLT_FULL_cff::depth, depthEta16, depths, eta, HcalParameters::etaMax, HLT_FULL_cff::etaMax, HcalParameters::etaMin, hpar, maxDepth, printTileHB(), and printTileHE().

                                           {
 
  std::cout << "Tile Information for HB from " << hpar->etaMin[0] << " to " << hpar->etaMax[0] << "\n\n";
  for (int eta=hpar->etaMin[0]; eta<= hpar->etaMax[0]; eta++) {
    int dmax = 1;
    if (depths[0][eta-1] < 17) dmax = 2;
    if (maxDepth[0] > 2) {
      if (eta == hpar->etaMin[1]) dmax = depthEta16[0];
      else                        dmax = maxDepth[0];
    }
    for (int depth=1; depth<=dmax; depth++) 
      printTileHB(eta, depth);
  }

  std::cout << "\nTile Information for HE from " << hpar->etaMin[1] << " to " << hpar->etaMax[1] << "\n\n";
  for (int eta=hpar->etaMin[1]; eta<= hpar->etaMax[1]; eta++) {
    int dmin=1, dmax=3;
    if (eta == hpar->etaMin[1]) {
      dmax = dmin = depthEta16[1];
    } else if (depths[0][eta-1] > 18) {
      dmax = 1;
    } else if (depths[1][eta-1] > 18) {
      dmax = 2;
    }
    for (int depth=dmin; depth<=dmax; depth++)
      printTileHE(eta, depth);
  }
}

int HcalDDDSimConstants::unitPhi	(	int	det,
		int	etaR
	)		const

Definition at line 485 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, HcalForward, hpar, HcalParameters::phibin, and HcalParameters::phitable.

Referenced by cell(), HcalDDDRecConstants::getHCID(), getPhiCons(), HcalCellTypes(), HcalDDDRecConstants::HcalCellTypes(), HcalDDDRecConstants::initialize(), and HcalNumberingFromDDD::unitID().

                                                        {

  double dphi = (det == static_cast<int>(HcalForward)) ? hpar->phitable[etaR-hpar->etaMin[2]] : hpar->phibin[etaR-1];
  return unitPhi(dphi);
}

int HcalDDDSimConstants::unitPhi

(

double

dphi

)

const

Definition at line 491 of file HcalDDDSimConstants.cc.

References M_PI, and units().

                                                  {

  const double fiveDegInRad = 2*M_PI/72;
  int units = int(dphi/fiveDegInRad+0.5);
  return units;
}

Member Data Documentation

int HcalDDDSimConstants::depthEta16[2]

private

Definition at line 89 of file HcalDDDSimConstants.h.

Referenced by getDepthEta16(), getEtaDepth(), initialize(), and printTiles().

int HcalDDDSimConstants::depthEta29[2]

private

Definition at line 90 of file HcalDDDSimConstants.h.

Referenced by getDepthEta29(), getEtaDepth(), and initialize().

std::vector<int> HcalDDDSimConstants::depths[nDepthMax]

private

Definition at line 81 of file HcalDDDSimConstants.h.

Referenced by initialize(), printTileHB(), printTileHE(), and printTiles().

double HcalDDDSimConstants::dlShort

private

Definition at line 88 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getEta(), HcalCellTypes(), and initialize().

double HcalDDDSimConstants::dzVcal

private

Definition at line 87 of file HcalDDDSimConstants.h.

Referenced by HcalCellTypes(), and initialize().

double HcalDDDSimConstants::etaHO[4]

private

Definition at line 85 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getEta(), and initialize().

const HcalParameters* HcalDDDSimConstants::hpar

private

Definition at line 74 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getConstHBHE(), getDetEta(), getEta(), getEtaDepth(), getEtaHO(), getEtaTableHF(), getGain(), getGparHF(), getiEtaRange(), getLayer0Wt(), getPhiCons(), getPhiTableHF(), getRTableHF(), getShift(), HcalCellTypes(), initialize(), layerGroup(), layerGroupSize(), printTileHB(), printTileHE(), printTiles(), and unitPhi().

std::vector<HcalDetId> HcalDDDSimConstants::idHF2QIE

private

Definition at line 91 of file HcalDDDSimConstants.h.

Referenced by getIdHF2QIE(), HcalCellTypes(), initialize(), and maxHFDepth().

std::vector<int> HcalDDDSimConstants::maxDepth

private

Definition at line 77 of file HcalDDDSimConstants.h.

Referenced by cell(), getEtaDepth(), getMaxDepth(), HcalCellTypes(), initialize(), maxHFDepth(), and printTiles().

int HcalDDDSimConstants::nDepth

private

Definition at line 82 of file HcalDDDSimConstants.h.

const int HcalDDDSimConstants::nDepthMax =9

staticprivate

Definition at line 76 of file HcalDDDSimConstants.h.

int HcalDDDSimConstants::nEta

private

Definition at line 78 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getDetEta(), getEta(), and initialize().

int HcalDDDSimConstants::nmodHB

private

Definition at line 83 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), HcalCellTypes(), and initialize().

int HcalDDDSimConstants::nmodHE

private

Definition at line 84 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), HcalCellTypes(), and initialize().

int HcalDDDSimConstants::nPhiF

private

Definition at line 80 of file HcalDDDSimConstants.h.

Referenced by initialize().

int HcalDDDSimConstants::nR

private

Definition at line 79 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getEta(), and initialize().

int HcalDDDSimConstants::nzHB

private

Definition at line 83 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), HcalCellTypes(), and initialize().

int HcalDDDSimConstants::nzHE

private

Definition at line 84 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), HcalCellTypes(), and initialize().

double HcalDDDSimConstants::rminHO

private

Definition at line 85 of file HcalDDDSimConstants.h.

Referenced by getEtaHO(), and initialize().

double HcalDDDSimConstants::zVcal

private

Definition at line 86 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getDetEta(), getEta(), and initialize().