HcalFlexiHardcodeGeometryLoader Class Reference

#include <HcalFlexiHardcodeGeometryLoader.h>

Classes
struct	HBHOCellParameters
struct	HECellParameters
struct	HFCellParameters

Public Member Functions
	HcalFlexiHardcodeGeometryLoader (const edm::ParameterSet &)
CaloSubdetectorGeometry *	load (const HcalTopology &fTopology)

Private Member Functions
void	fillHBHO (CaloSubdetectorGeometry *fGeometry, const std::vector< HBHOCellParameters > &fCells, bool fHB)
void	fillHE (CaloSubdetectorGeometry *fGeometry, const std::vector< HECellParameters > &fCells)
void	fillHF (CaloSubdetectorGeometry *fGeometry, const std::vector< HFCellParameters > &fCells)
std::vector< HBHOCellParameters >	makeHBCells (const HcalTopology &topology)
std::vector< HECellParameters >	makeHECells (const HcalTopology &topology)
std::vector< HECellParameters >	makeHECells_H2 ()
std::vector< HFCellParameters >	makeHFCells ()
std::vector< HBHOCellParameters >	makeHOCells ()

Private Attributes
double	DEGREE2RAD
std::vector< std::vector< int > >	m_segmentation
int	MAX_HCAL_PHI

Detailed Description

Date:

2012/03/22 10:46:31

Revision:

1.3

Author

F.Ratnikov, UMd

Definition at line 18 of file HcalFlexiHardcodeGeometryLoader.h.

Constructor & Destructor Documentation

HcalFlexiHardcodeGeometryLoader::HcalFlexiHardcodeGeometryLoader

(

const edm::ParameterSet &

ps

)

Definition at line 14 of file HcalFlexiHardcodeGeometryLoader.cc.

References gather_cfg::cout, DEGREE2RAD, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, gen::k, M_PI, m_segmentation, MAX_HCAL_PHI, and mergeVDriftHistosByStation::name.

                                                                                          {

  MAX_HCAL_PHI = 72;
  DEGREE2RAD = M_PI / 180.;

  edm::ParameterSet ps0 = ps.getParameter<edm::ParameterSet>("HcalReLabel");
  bool relabel_= ps0.getUntrackedParameter<bool>("RelabelHits",false);
  if (relabel_) {
    edm::ParameterSet ps1 = ps0.getUntrackedParameter<edm::ParameterSet>("RelabelRules");
    m_segmentation.resize(29);
    for (int i=0; i<29; i++) {
      char name[10];
      snprintf(name,10,"Eta%d",i+1);
      if (i>0) {
    m_segmentation[i]=ps1.getUntrackedParameter<std::vector<int> >(name,m_segmentation[i-1]);
      } else {
    m_segmentation[i]=ps1.getUntrackedParameter<std::vector<int> >(name);
      }
#ifdef DebugLog
      std::cout << name;
      for (unsigned int k=0; k<m_segmentation[i].size(); ++k) {
    std::cout << " [" << k << "] " << m_segmentation[i][k];
      }
      std::cout << std::endl;
#endif
    }
  }

}

Member Function Documentation

void HcalFlexiHardcodeGeometryLoader::fillHBHO ( CaloSubdetectorGeometry *  fGeometry, const std::vector< HBHOCellParameters > &  fCells, bool  fHB  )

private

Definition at line 172 of file HcalFlexiHardcodeGeometryLoader.cc.

References funct::cos(), gather_cfg::cout, DEGREE2RAD, HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::depth, HcalDetId::depth(), HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::dphi, HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::eta, PV3DBase< T, PVType, FrameType >::eta(), HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::etaMax, HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::etaMin, CaloCellGeometry::getParmPtr(), HcalBarrel, HcalOuter, HcalDetId::ieta(), HcalDetId::iphi(), MAX_HCAL_PHI, CaloSubdetectorGeometry::newCell(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::phiFirst, HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::phiStep, HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::rMax, HcalFlexiHardcodeGeometryLoader::HBHOCellParameters::rMin, funct::sin(), x, detailsBasic3DVector::y, detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by load().

                                                                                                                                                                         {

  for (size_t iCell = 0; iCell < fCells.size(); ++iCell) {
    const HcalFlexiHardcodeGeometryLoader::HBHOCellParameters& param = fCells[iCell];
    for (int iPhi = param.phiFirst; iPhi <= MAX_HCAL_PHI; iPhi += param.phiStep) {
      for (int iside = -1; iside <= 1; iside += 2) { // both detector sides are identical
    HcalDetId hid (fHB ? HcalBarrel : HcalOuter, param.eta*iside, iPhi, param.depth);
    float phiCenter = ((iPhi-1)*360./MAX_HCAL_PHI + 0.5*param.dphi) * DEGREE2RAD; // middle of the cell
    float etaCenter = 0.5*(param.etaMin + param.etaMax);
    float x = param.rMin* cos (phiCenter);
    float y = param.rMin* sin (phiCenter);
    float z = iside * param.rMin * sinh(etaCenter);
    // make cell geometry
    GlobalPoint refPoint (x,y,z); // center of the cell's face
    std::vector<CCGFloat> cellParams;
    cellParams.reserve (5);
    cellParams.push_back (0.5 * (param.etaMax - param.etaMin)); // deta_half
    cellParams.push_back (0.5 * param.dphi * DEGREE2RAD);  // dphi_half
    cellParams.push_back (0.5 * (param.rMax - param.rMin) * cosh (etaCenter)); // dr_half
    cellParams.push_back ( fabs( refPoint.eta() ) ) ;
    cellParams.push_back ( fabs( refPoint.z() ) ) ;
#ifdef DebugLog
    std::cout << "HcalFlexiHardcodeGeometryLoader::fillHBHO-> " << hid << hid.ieta() << '/' << hid.iphi() << '/' << hid.depth() << refPoint << '/' << cellParams [0] << '/' << cellParams [1] << '/' << cellParams [2] << std::endl;
#endif
    fGeometry->newCell(refPoint,  refPoint,  refPoint, 
               CaloCellGeometry::getParmPtr(cellParams, 
                            fGeometry->parMgr(), 
                            fGeometry->parVecVec()),
               hid ) ;
      }
    }
  }
}

void HcalFlexiHardcodeGeometryLoader::fillHE ( CaloSubdetectorGeometry *  fGeometry, const std::vector< HECellParameters > &  fCells  )

private

Definition at line 397 of file HcalFlexiHardcodeGeometryLoader.cc.

References funct::cos(), gather_cfg::cout, DEGREE2RAD, HcalFlexiHardcodeGeometryLoader::HECellParameters::depth, HcalFlexiHardcodeGeometryLoader::HECellParameters::dphi, HcalFlexiHardcodeGeometryLoader::HECellParameters::eta, PV3DBase< T, PVType, FrameType >::eta(), HcalFlexiHardcodeGeometryLoader::HECellParameters::etaMax, HcalFlexiHardcodeGeometryLoader::HECellParameters::etaMin, CaloCellGeometry::getParmPtr(), HcalEndcap, MAX_HCAL_PHI, CaloSubdetectorGeometry::newCell(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), perp(), HcalFlexiHardcodeGeometryLoader::HECellParameters::phiFirst, HcalFlexiHardcodeGeometryLoader::HECellParameters::phiStep, funct::sin(), x, detailsBasic3DVector::y, detailsBasic3DVector::z, PV3DBase< T, PVType, FrameType >::z(), HcalFlexiHardcodeGeometryLoader::HECellParameters::zMax, and HcalFlexiHardcodeGeometryLoader::HECellParameters::zMin.

Referenced by load().

                                                                                                                                                           {

  for (size_t iCell = 0; iCell < fCells.size(); ++iCell) {
    const HcalFlexiHardcodeGeometryLoader::HECellParameters& param = fCells[iCell];
    for (int iPhi = param.phiFirst; iPhi <= MAX_HCAL_PHI; iPhi += param.phiStep) {
      for (int iside = -1; iside <= 1; iside += 2) { // both detector sides are identical
    HcalDetId hid (HcalEndcap, param.eta*iside, iPhi, param.depth);
    float phiCenter = ((iPhi-1)*360./MAX_HCAL_PHI + 0.5*param.dphi) * DEGREE2RAD; // middle of the cell
    float etaCenter = 0.5 * (param.etaMin + param.etaMax);

    float perp = param.zMin / sinh (etaCenter);
    float x = perp * cos (phiCenter);
    float y = perp * sin (phiCenter);
    float z = iside * param.zMin;
    // make cell geometry
    GlobalPoint refPoint (x,y,z); // center of the cell's face
    std::vector<CCGFloat> cellParams;
    cellParams.reserve (5);
    cellParams.push_back (0.5 * (param.etaMax - param.etaMin)); //deta_half
    cellParams.push_back (0.5 * param.dphi * DEGREE2RAD);  // dphi_half
    cellParams.push_back (-0.5 * (param.zMax - param.zMin) / tanh (etaCenter)); // dz_half, "-" means edges in Z
    cellParams.push_back ( fabs( refPoint.eta() ) ) ;
    cellParams.push_back ( fabs( refPoint.z() ) ) ;
#ifdef DebugLog
    std::cout << "HcalFlexiHardcodeGeometryLoader::fillHE-> " << hid << refPoint << '/' << cellParams [0] << '/' << cellParams [1] << '/' << cellParams [2] << std::endl;
#endif
    fGeometry->newCell(refPoint,  refPoint,  refPoint, 
               CaloCellGeometry::getParmPtr(cellParams, 
                            fGeometry->parMgr(), 
                            fGeometry->parVecVec()),
               hid ) ;
      }
    }
  }
}

void HcalFlexiHardcodeGeometryLoader::fillHF ( CaloSubdetectorGeometry *  fGeometry, const std::vector< HFCellParameters > &  fCells  )

private

Definition at line 433 of file HcalFlexiHardcodeGeometryLoader.cc.

References funct::cos(), gather_cfg::cout, DEGREE2RAD, HcalFlexiHardcodeGeometryLoader::HFCellParameters::depth, HcalFlexiHardcodeGeometryLoader::HFCellParameters::dphi, HcalFlexiHardcodeGeometryLoader::HFCellParameters::eta, PV3DBase< T, PVType, FrameType >::eta(), CaloCellGeometry::getParmPtr(), HcalForward, SurfaceOrientation::inner, MAX_HCAL_PHI, CaloSubdetectorGeometry::newCell(), SurfaceOrientation::outer, CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), perp(), HcalFlexiHardcodeGeometryLoader::HFCellParameters::phiFirst, HcalFlexiHardcodeGeometryLoader::HFCellParameters::phiStep, HcalFlexiHardcodeGeometryLoader::HFCellParameters::rMax, HcalFlexiHardcodeGeometryLoader::HFCellParameters::rMin, funct::sin(), x, detailsBasic3DVector::y, detailsBasic3DVector::z, PV3DBase< T, PVType, FrameType >::z(), HcalFlexiHardcodeGeometryLoader::HFCellParameters::zMax, and HcalFlexiHardcodeGeometryLoader::HFCellParameters::zMin.

Referenced by load().

                                                                                                                                                           {

  for (size_t iCell = 0; iCell < fCells.size(); ++iCell) {
    const HcalFlexiHardcodeGeometryLoader::HFCellParameters& param = fCells[iCell];
    for (int iPhi = param.phiFirst; iPhi <= MAX_HCAL_PHI; iPhi += param.phiStep) {
      for (int iside = -1; iside <= 1; iside += 2) { // both detector sides are identical
    HcalDetId hid (HcalForward, param.eta*iside, iPhi, param.depth);
    float phiCenter = ((iPhi-1)*360./MAX_HCAL_PHI + 0.5*param.dphi) * DEGREE2RAD; // middle of the cell
    GlobalPoint inner (param.rMin, 0, param.zMin);
    GlobalPoint outer (param.rMax, 0, param.zMin);
    float iEta = inner.eta();
    float oEta = outer.eta();
    float etaCenter = 0.5 * ( iEta + oEta );
    
    float perp = param.zMin / sinh (etaCenter);
    float x = perp * cos (phiCenter);
    float y = perp * sin (phiCenter);
    float z = iside * param.zMin;
    // make cell geometry
    GlobalPoint refPoint (x,y,z); // center of the cell's face
    std::vector<CCGFloat> cellParams;
    cellParams.reserve (5);
    cellParams.push_back (0.5 * ( iEta - oEta )); // deta_half
    cellParams.push_back (0.5 * param.dphi * DEGREE2RAD);  // dphi_half
    cellParams.push_back (0.5 * (param.zMax - param.zMin)); // dz_half
    cellParams.push_back ( fabs( refPoint.eta()));
    cellParams.push_back ( fabs( refPoint.z() ) ) ;
#ifdef DebugLog
    std::cout << "HcalFlexiHardcodeGeometryLoader::fillHF-> " << hid << refPoint << '/' << cellParams [0] << '/' << cellParams [1] << '/' << cellParams [2] << std::endl;
#endif  
    fGeometry->newCell(refPoint,  refPoint,  refPoint, 
               CaloCellGeometry::getParmPtr(cellParams, 
                            fGeometry->parMgr(), 
                            fGeometry->parVecVec()),
               hid ) ;
      }
    }
  }
}

CaloSubdetectorGeometry * HcalFlexiHardcodeGeometryLoader::load

(

const HcalTopology &

fTopology

)

Definition at line 44 of file HcalFlexiHardcodeGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), fillHBHO(), fillHE(), fillHF(), HcalTopologyMode::H2, HcalGeometry::k_NumberOfParametersPerShape, makeHBCells(), makeHECells(), makeHECells_H2(), makeHFCells(), makeHOCells(), HcalTopology::mode(), HcalTopology::ncells(), CaloSubdetectorGeometry::numberOfShapes(), and CaloSubdetectorGeometry::parMgr().

Referenced by HcalHardcodeGeometryEP::produceIdeal().

                                                                                            {
  CaloSubdetectorGeometry* hcalGeometry = new HcalGeometry (fTopology);
  if( 0 == hcalGeometry->cornersMgr() ) hcalGeometry->allocateCorners ( fTopology.ncells() );
  if( 0 == hcalGeometry->parMgr() ) hcalGeometry->allocatePar (hcalGeometry->numberOfShapes(),
                                   HcalGeometry::k_NumberOfParametersPerShape ) ;
  if (fTopology.mode() == HcalTopologyMode::H2) {  // TB geometry
    fillHBHO (hcalGeometry, makeHBCells(fTopology), true);
    fillHBHO (hcalGeometry, makeHOCells(), false);
    fillHE (hcalGeometry, makeHECells_H2());
 } else { // regular geometry
    fillHBHO (hcalGeometry, makeHBCells(fTopology), true);
    fillHBHO (hcalGeometry, makeHOCells(), false);
    fillHF (hcalGeometry, makeHFCells());
    fillHE (hcalGeometry, makeHECells(fTopology));
  }
  return hcalGeometry;
}

std::vector< HcalFlexiHardcodeGeometryLoader::HBHOCellParameters > HcalFlexiHardcodeGeometryLoader::makeHBCells ( const HcalTopology &  topology )

private

Definition at line 64 of file HcalFlexiHardcodeGeometryLoader.cc.

References gather_cfg::cout, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, i, m_segmentation, HcalTopology::mode(), query::result, and HcalTopologyMode::SLHC.

Referenced by load().

                                                                                                                                         {

  const float HBRMIN = 181.1;
  const float HBRMAX = 288.8;
    
  float normalDepths[2] = {HBRMIN, HBRMAX};
  float ring15Depths[3] = {HBRMIN, 258.4, HBRMAX};
  float ring16Depths[3] = {HBRMIN, 190.4, 232.6};
  float layerDepths[18] = {HBRMIN, 188.7, 194.7, 200.7, 206.7, 212.7, 218.7,
               224.7, 230.7, 236.7, 242.7, 249.3, 255.9, 262.5,
               269.1, 275.7, 282.3, HBRMAX};
  float slhcDepths[4]   = {HBRMIN, 214., 239., HBRMAX};
#ifdef DebugLog
  std::cout <<"FlexiGeometryLoader called for "<< topology.mode() << ":" << HcalTopologyMode::SLHC << std::endl;
#endif
  std::vector <HcalFlexiHardcodeGeometryLoader::HBHOCellParameters> result;
  for(int iring = 1; iring <= 16; ++iring) {
    std::vector<float> depths;
    if (topology.mode() != HcalTopologyMode::SLHC) {
      if (iring == 15) {
    for (int i=0; i<3; ++i) depths.push_back(ring15Depths[i]);
      } else if (iring == 16) {
    for (int i=0; i<3; ++i) depths.push_back(ring16Depths[i]);
      } else {
    for (int i=0; i<2; ++i) depths.push_back(normalDepths[i]);
      }
    } else {
      if (m_segmentation.size() >= (unsigned int)(iring)) {
    int depth = m_segmentation[iring-1][0];
    depths.push_back(layerDepths[depth]);
    int layer = 1;
    for (unsigned int i=1; i<m_segmentation[iring-1].size(); ++i) {
      if (depth != m_segmentation[iring-1][i]) {
        depth = m_segmentation[iring-1][i];
        layer = i;
        if (iring != 16 || depth < 3)
          depths.push_back(layerDepths[depth]);
      }
      if (i >= 17) break;
    }
    if (layer <= 17) depths.push_back(HBRMAX);
      } else {
    for (int i=0; i<4; ++i) {
      if (iring != 16 || i < 3) {
        depths.push_back(slhcDepths[i]);
      }
    }
      }
    }
    unsigned int ndepth=depths.size()-1;
    unsigned int startingDepth=1;
    float etaMin=(iring-1)*0.087;
    float etaMax=iring*0.087;
    // topology.depthBinInformation(HcalBarrel, iring, ndepth, startingDepth);
#ifdef DebugLog
    std::cout << "HBRing " << iring << " eta " << etaMin << ":" << etaMax << " depths " << ndepth << ":" << startingDepth;
    for (unsigned int i=0; i<depths.size(); ++i) std::cout << ":" << depths[i];
    std::cout << "\n";
#endif
    for (unsigned int idepth = startingDepth; idepth <= ndepth; ++idepth) {
      float rmin = depths[idepth-1];
      float rmax = depths[idepth];
#ifdef DebugLog
      std::cout << "HB " << idepth << " R " << rmin << ":" << rmax << "\n";
#endif
      result.push_back(HcalFlexiHardcodeGeometryLoader::HBHOCellParameters(iring, (int)idepth, 1, 1, 5, rmin, rmax, etaMin, etaMax));
    }
  }
  return result;
}

std::vector< HcalFlexiHardcodeGeometryLoader::HECellParameters > HcalFlexiHardcodeGeometryLoader::makeHECells ( const HcalTopology &  topology )

private

Definition at line 208 of file HcalFlexiHardcodeGeometryLoader.cc.

References gather_cfg::cout, SiPixelRawToDigiRegional_cfi::deltaPhi, HcalTopology::firstHEDoublePhiRing(), HcalTopology::firstHETripleDepthRing(), i, HcalTopology::lastHERing(), m_segmentation, HcalTopology::mode(), query::result, HcalTopologyMode::SLHC, SiStripMonitorClusterAlca_cfi::zmax, and SiStripMonitorClusterAlca_cfi::zmin.

Referenced by load().

                                                                                                                                      {

  std::vector <HcalFlexiHardcodeGeometryLoader::HECellParameters> result;
  const float HEZMIN = 400.458;
  const float HEZMID = 436.168;
  const float HEZMAX = 549.268;
  float normalDepths[3] = {HEZMIN, HEZMID, HEZMAX};
  float tripleDepths[4] = {HEZMIN, 418.768, HEZMID, HEZMAX};
  float slhcDepths[5]   = {HEZMIN, 418.768, HEZMID, 493., HEZMAX};
  float ring16Depths[2] = {418.768,470.968};
  float ring16slhcDepths[3] = {418.768, 450., 470.968};
  float ring17Depths[2] = {409.698,514.468};
  float ring17slhcDepths[5] = {409.698, 435., 460., 495., 514.468};
  float ring18Depths[3] = {391.883,427.468,540.568};
  float ring18slhcDepths[5] = {391.883, 439.,  467., 504. , 540.568};
  float etaBounds[] = {0.087*15, 0.087*16, 0.087*17, 0.087*18,  0.087*19,
               1.74, 1.83,  1.93, 2.043, 2.172, 2.322, 2.500,
               2.650, 2.868, 3.000};
  float layerDepths[19] = {HEZMIN, 408.718, 416.978, 425.248, 433.508, 441.768,
               450.038,458.298, 466.558, 474.828, 483.088, 491.348,
               499.618,507.878, 516.138, 524.398, 532.668, 540.928,
               HEZMAX};

  // count by ring - 16
  for(int iringm16=0; iringm16 <= 13; ++iringm16) {
    int iring = iringm16 + 16;
    std::vector<float> depths;
    unsigned int startingDepth = 1;
    if (topology.mode() != HcalTopologyMode::SLHC) {
      if (iring == 16)     
    {for (int i=0; i<2; ++i) depths.push_back(ring16Depths[i]); startingDepth = 3;}
      else if (iring == 17) 
    for (int i=0; i<2; ++i) depths.push_back(ring17Depths[i]);
      else if (iring == 18) 
    for (int i=0; i<3; ++i) depths.push_back(ring18Depths[i]);
      else if (iring == topology.lastHERing()) 
    for (int i=0; i<3; ++i) depths.push_back(tripleDepths[i]);
      else if (iring >= topology.firstHETripleDepthRing())
    for (int i=0; i<4; ++i) depths.push_back(tripleDepths[i]);
      else
    for (int i=0; i<3; ++i) depths.push_back(normalDepths[i]);
    } else {
      if (m_segmentation.size() >= (unsigned int)(iring)) {
    int depth = m_segmentation[iring-1][0];
    if (iring == 16)      depths.push_back(ring16Depths[0]);
    else if (iring == 17) depths.push_back(ring17Depths[0]);
    else if (iring == 18) depths.push_back(ring18Depths[0]);
    else                  depths.push_back(layerDepths[depth]);
    int layer = 1;
    float lastDepth = depths[0];
    for (unsigned int i=1; i<m_segmentation[iring-1].size(); ++i) {
      if (depth != m_segmentation[iring-1][i]) {
        depth = m_segmentation[iring-1][i];
        layer = i;
        if (layerDepths[depth] > lastDepth && (iring != 16 || depth > 3)) {
          depths.push_back(layerDepths[depth]);
          lastDepth = layerDepths[depth];
        }
      }
    }
    if (layer <= 17) depths.push_back(HEZMAX);
    if (iring == 16) startingDepth = 3;
      } else {
    if (iring == 16)     {for (int i=0; i<3; ++i) depths.push_back(ring16slhcDepths[i]); startingDepth = 3;}
    else if (iring == 17) for (int i=0; i<5; ++i) depths.push_back(ring17slhcDepths[i]);
    else if (iring == 18) for (int i=0; i<5; ++i) depths.push_back(ring18slhcDepths[i]);
    else                  for (int i=0; i<5; ++i) depths.push_back(slhcDepths[i]);
      }
    }
    float etamin = etaBounds[iringm16];
    float etamax = etaBounds[iringm16+1];
    unsigned int ndepth = depths.size()-1;
    //    topology.depthBinInformation(HcalEndcap, iring, ndepth, startingDepth);
#ifdef DebugLog
    std::cout << "HERing " << iring << " eta " << etamin << ":" << etamax << " depths " << ndepth << ":" << startingDepth;
    for (unsigned int i=0; i<depths.size(); ++i) std::cout << ":" << depths[i];
    std::cout << "\n";
#endif
    for (unsigned int idepth = 0; idepth < ndepth; ++idepth) {
      int depthIndex = (int)(idepth + startingDepth);
      float zmin = depths[idepth];
      float zmax = depths[idepth+1];
      if (depthIndex <= 7) {
#ifdef DebugLog
    std::cout << "HE Depth " << idepth << ":" << depthIndex << " Z " << zmin << ":" << zmax << "\n";
#endif
    int stepPhi = (iring >= topology.firstHEDoublePhiRing() ? 2 : 1);
    int deltaPhi =  (iring >= topology.firstHEDoublePhiRing() ? 10 : 5);
    if (topology.mode() != HcalTopologyMode::SLHC &&
        iring == topology.lastHERing()-1 && idepth == ndepth-1) {
#ifdef DebugLog
      std::cout << "HE iEta " << iring << " Depth " << depthIndex << " Eta " << etamin << ":" << etaBounds[iringm16+2] << std::endl;
#endif
      result.push_back(HcalFlexiHardcodeGeometryLoader::HECellParameters(iring, depthIndex, 1, stepPhi, deltaPhi, zmin, zmax, etamin, etaBounds[iringm16+2]));
    } else {
#ifdef DebugLog
      std::cout << "HE iEta " << iring << " Depth " << depthIndex << " Eta " << etamin << ":" << etamax << std::endl;
#endif
      result.push_back(HcalFlexiHardcodeGeometryLoader::HECellParameters(iring, depthIndex, 1, stepPhi, deltaPhi, zmin, zmax, etamin, etamax));
    }
      }
    }
  }

  return result;
}

std::vector< HcalFlexiHardcodeGeometryLoader::HECellParameters > HcalFlexiHardcodeGeometryLoader::makeHECells_H2 ( )

private

Definition at line 317 of file HcalFlexiHardcodeGeometryLoader.cc.

References i, and query::result.

Referenced by load().

                                                                                                             {

  const float HEZMIN_H2 = 400.715;
  const float HEZMID_H2 = 436.285;
  const float HEZMAX_H2 = 541.885;
    
  HcalFlexiHardcodeGeometryLoader::HECellParameters cells [] = {
    // eta, depth, firstPhi, stepPhi, deltaPhi, zMin, zMax, etaMin, etaMax
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 16, 3, 1, 1, 5, 409.885,   462.685,   1.305, 1.373),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 17, 1, 1, 1, 5, HEZMIN_H2, 427.485,   1.373, 1.444),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 17, 2, 1, 1, 5, 427.485,   506.685,   1.373, 1.444),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 18, 1, 1, 1, 5, HEZMIN_H2, HEZMID_H2, 1.444, 1.521),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 18, 2, 1, 1, 5, HEZMID_H2, 524.285,   1.444, 1.521),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 19, 1, 1, 1, 5, HEZMIN_H2, HEZMID_H2, 1.521, 1.603),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 19, 2, 1, 1, 5, HEZMID_H2, HEZMAX_H2, 1.521, 1.603),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 20, 1, 1, 1, 5, HEZMIN_H2, HEZMID_H2, 1.603, 1.693),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 20, 2, 1, 1, 5, HEZMID_H2, HEZMAX_H2, 1.603, 1.693),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 21, 1, 1, 2, 5, HEZMIN_H2, HEZMID_H2, 1.693, 1.79),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 21, 2, 1, 2, 5, HEZMID_H2, HEZMAX_H2, 1.693, 1.79),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 22, 1, 1, 2,10, HEZMIN_H2, HEZMID_H2, 1.79, 1.88),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 22, 2, 1, 2,10, HEZMID_H2, HEZMAX_H2, 1.79, 1.88),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 23, 1, 1, 2,10, HEZMIN_H2, HEZMID_H2, 1.88, 1.98),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 23, 2, 1, 2,10, HEZMID_H2, HEZMAX_H2, 1.88, 1.98),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 24, 1, 1, 2,10, HEZMIN_H2, 418.685,   1.98, 2.09),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 24, 2, 1, 2,10, 418.685,   HEZMID_H2, 1.98, 2.09),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 24, 3, 1, 2,10, HEZMID_H2, HEZMAX_H2, 1.98, 2.09),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 25, 1, 1, 2,10, HEZMIN_H2, 418.685,   2.09, 2.21),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 25, 2, 1, 2,10, 418.685,   HEZMID_H2, 2.09, 2.21),
    HcalFlexiHardcodeGeometryLoader::HECellParameters ( 25, 3, 1, 2,10, HEZMID_H2, HEZMAX_H2, 2.09, 2.21)
  };
  int nCells = sizeof(cells)/sizeof(HcalFlexiHardcodeGeometryLoader::HECellParameters);
  std::vector <HcalFlexiHardcodeGeometryLoader::HECellParameters> result;
  result.reserve (nCells);
  for (int i = 0; i < nCells; ++i) result.push_back (cells[i]);
  return result;
}

std::vector< HcalFlexiHardcodeGeometryLoader::HFCellParameters > HcalFlexiHardcodeGeometryLoader::makeHFCells ( )

private

Definition at line 355 of file HcalFlexiHardcodeGeometryLoader.cc.

References i, and query::result.

Referenced by load().

                                                                                                          {

  const float HFZMIN1 = 1115.;
  const float HFZMIN2 = 1137.;
  const float HFZMAX = 1280.1;
    
  HcalFlexiHardcodeGeometryLoader::HFCellParameters cells [] = {
    // eta, depth, firstPhi, stepPhi, deltaPhi, zMin, zMax, rMin, rMax
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (29, 1, 1, 2, 10, HFZMIN1, HFZMAX,116.2,130.0),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (29, 2, 1, 2, 10, HFZMIN2, HFZMAX,116.2,130.0),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (30, 1, 1, 2, 10, HFZMIN1, HFZMAX, 97.5,116.2),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (30, 2, 1, 2, 10, HFZMIN2, HFZMAX, 97.5,116.2),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (31, 1, 1, 2, 10, HFZMIN1, HFZMAX, 81.8, 97.5),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (31, 2, 1, 2, 10, HFZMIN2, HFZMAX, 81.8, 97.5),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (32, 1, 1, 2, 10, HFZMIN1, HFZMAX, 68.6, 81.8),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (32, 2, 1, 2, 10, HFZMIN2, HFZMAX, 68.6, 81.8),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (33, 1, 1, 2, 10, HFZMIN1, HFZMAX, 57.6, 68.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (33, 2, 1, 2, 10, HFZMIN2, HFZMAX, 57.6, 68.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (34, 1, 1, 2, 10, HFZMIN1, HFZMAX, 48.3, 57.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (34, 2, 1, 2, 10, HFZMIN2, HFZMAX, 48.3, 57.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (35, 1, 1, 2, 10, HFZMIN1, HFZMAX, 40.6, 48.3),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (35, 2, 1, 2, 10, HFZMIN2, HFZMAX, 40.6, 48.3),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (36, 1, 1, 2, 10, HFZMIN1, HFZMAX, 34.0, 40.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (36, 2, 1, 2, 10, HFZMIN2, HFZMAX, 34.0, 40.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (37, 1, 1, 2, 10, HFZMIN1, HFZMAX, 28.6, 34.0),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (37, 2, 1, 2, 10, HFZMIN2, HFZMAX, 28.6, 34.0),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (38, 1, 1, 2, 10, HFZMIN1, HFZMAX, 24.0, 28.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (38, 2, 1, 2, 10, HFZMIN2, HFZMAX, 24.0, 28.6),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (39, 1, 1, 2, 10, HFZMIN1, HFZMAX, 20.1, 24.0),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (39, 2, 1, 2, 10, HFZMIN2, HFZMAX, 20.1, 24.0),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (40, 1, 3, 4, 20, HFZMIN1, HFZMAX, 16.9, 20.1),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (40, 2, 3, 4, 20, HFZMIN2, HFZMAX, 16.9, 20.1),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (41, 1, 3, 4, 20, HFZMIN1, HFZMAX, 12.5, 16.9),
    HcalFlexiHardcodeGeometryLoader::HFCellParameters (41, 2, 3, 4, 20, HFZMIN2, HFZMAX, 12.5, 16.9)
  };
  int nCells = sizeof(cells)/sizeof(HcalFlexiHardcodeGeometryLoader::HFCellParameters);
  std::vector <HcalFlexiHardcodeGeometryLoader::HFCellParameters> result;
  result.reserve (nCells);
  for (int i = 0; i < nCells; ++i) result.push_back (cells[i]);
  return result;
}

std::vector< HcalFlexiHardcodeGeometryLoader::HBHOCellParameters > HcalFlexiHardcodeGeometryLoader::makeHOCells ( )

private

Definition at line 138 of file HcalFlexiHardcodeGeometryLoader.cc.

References i, and query::result.

Referenced by load().

                                                                                                            {
  const float HORMIN0 = 390.0;
  const float HORMIN1 = 412.6;
  const float HORMAX = 413.6;
  
  HcalFlexiHardcodeGeometryLoader::HBHOCellParameters cells [] = {
    // eta, depth, firstPhi, stepPhi, deltaPhi, rMin, rMax, etaMin, etaMax
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 1, 4, 1, 1, 5, HORMIN0, HORMAX, 0.087*0, 0.087*1),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 2, 4, 1, 1, 5, HORMIN0, HORMAX, 0.087*1, 0.087*2),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 3, 4, 1, 1, 5, HORMIN0, HORMAX, 0.087*2, 0.087*3),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 4, 4, 1, 1, 5, HORMIN0, HORMAX, 0.087*3, 0.3075),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 5, 4, 1, 1, 5, HORMIN1, HORMAX, 0.3395,  0.087*5),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 6, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*5, 0.087*6),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 7, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*6, 0.087*7),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 8, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*7, 0.087*8),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters ( 9, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*8, 0.087*9),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters (10, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*9,  0.8494),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters (11, 4, 1, 1, 5, HORMIN1, HORMAX, 0.873, 0.087*11),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters (12, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*11, 0.087*12),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters (13, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*12, 0.087*13),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters (14, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*13, 0.087*14),
    HcalFlexiHardcodeGeometryLoader::HBHOCellParameters (15, 4, 1, 1, 5, HORMIN1, HORMAX, 0.087*14, 0.087*15)
  };
  int nCells = sizeof(cells)/sizeof(HcalFlexiHardcodeGeometryLoader::HBHOCellParameters);
  std::vector <HcalFlexiHardcodeGeometryLoader::HBHOCellParameters> result;
  result.reserve (nCells);
  for (int i = 0; i < nCells; ++i) result.push_back (cells[i]);
  return result;
}

Member Data Documentation

double HcalFlexiHardcodeGeometryLoader::DEGREE2RAD

private

Definition at line 86 of file HcalFlexiHardcodeGeometryLoader.h.

Referenced by fillHBHO(), fillHE(), fillHF(), and HcalFlexiHardcodeGeometryLoader().

std::vector<std::vector<int> > HcalFlexiHardcodeGeometryLoader::m_segmentation

private

Definition at line 88 of file HcalFlexiHardcodeGeometryLoader.h.

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

int HcalFlexiHardcodeGeometryLoader::MAX_HCAL_PHI

private

Definition at line 85 of file HcalFlexiHardcodeGeometryLoader.h.

Referenced by fillHBHO(), fillHE(), fillHF(), and HcalFlexiHardcodeGeometryLoader().