#include <CaloTowerHardcodeGeometryLoader.h>

Public Member Functions
std::auto_ptr < CaloSubdetectorGeometry >	load (const CaloTowerTopology limits, const HcalTopology hcaltopo, const HcalDDDRecConstants *hcons)

Private Member Functions
void	makeCell (uint32_t din, CaloSubdetectorGeometry *geom) const

Private Attributes
const HcalTopology *	m_hcaltopo

const HcalDDDRecConstants *	m_hcons

const CaloTowerTopology *	m_limits

std::vector< double >	theHBHEEtaBounds

std::vector< double >	theHFEtaBounds

Detailed Description

Author: J. Mans - Minnesota

Definition at line 15 of file CaloTowerHardcodeGeometryLoader.h.

Member Function Documentation

std::auto_ptr< CaloSubdetectorGeometry > CaloTowerHardcodeGeometryLoader::load	(	const CaloTowerTopology *	limits,
		const HcalTopology *	hcaltopo,
		const HcalDDDRecConstants *	hcons
	)

Definition at line 9 of file CaloTowerHardcodeGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), cropTnPTrees::din, relativeConstraints::geom, HcalDDDRecConstants::getEtaTable(), HcalDDDRecConstants::getEtaTableHF(), m_hcaltopo, m_hcons, m_limits, makeCell(), CaloTowerGeometry::numberOfCellsForCorners(), CaloTowerGeometry::numberOfParametersPerShape(), CaloTowerGeometry::numberOfShapes(), CaloSubdetectorGeometry::parMgr(), CaloTowerTopology::sizeForDenseIndexing(), theHBHEEtaBounds, and theHFEtaBounds.

Referenced by CaloTowerHardcodeGeometryEP::produce().

                                                                                                                                                                           {
 
   m_limits = limits;
   m_hcaltopo = hcaltopo;
   m_hcons = hcons;
 
   //get eta limits from hcal rec constants
   theHFEtaBounds   = m_hcons->getEtaTableHF();
   theHBHEEtaBounds = m_hcons->getEtaTable();
 
   CaloTowerGeometry* geom=new CaloTowerGeometry(m_limits);
 
   if( 0 == geom->cornersMgr() ) geom->allocateCorners ( 
      geom->numberOfCellsForCorners() ) ;
   if( 0 == geom->parMgr() ) geom->allocatePar (
      geom->numberOfParametersPerShape()*geom->numberOfShapes(),
      geom->numberOfParametersPerShape() ) ;
 
   // simple loop
   for (uint32_t din = 0; din < m_limits->sizeForDenseIndexing(); ++din) {
     makeCell(din, geom);
   }
   edm::LogInfo("Geometry") << "CaloTowersHardcodeGeometry made " << m_limits->sizeForDenseIndexing() << " towers.";
   return std::auto_ptr<CaloSubdetectorGeometry>(geom); 
 }

void CaloTowerHardcodeGeometryLoader::makeCell	(	uint32_t	din,
		CaloSubdetectorGeometry *	geom
	)		const

private

Definition at line 36 of file CaloTowerHardcodeGeometryLoader.cc.

References funct::abs(), CaloTowerTopology::convertCTtoHcal(), funct::cos(), CaloTowerTopology::detIdFromDenseIndex(), eta, CaloTowerTopology::firstHERing(), HcalTopology::firstHFRing(), CaloCellGeometry::getParmPtr(), CaloTowerDetId::ieta(), CaloTowerTopology::lastHERing(), m_hcaltopo, m_limits, M_PI, CaloSubdetectorGeometry::newCell(), HcalTopology::nPhiBins(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), phi, point, alignCSCRings::r, jetcorrextractor::sign(), funct::sin(), theHBHEEtaBounds, theHFEtaBounds, x, y, and z.

Referenced by load().

                                                              {
   const double EBradius = 143.0; // cm
   const double HOradius = 406.0+1.0;
   const double EEz = 320.0; // rough (cm)
   const double HEz = 568.0; // back (cm)
   const double HFz = 1100.0;
   const double HFthick = 165;
   // Tower 17 is the last EB tower
 
   //use CT topology to get the DetId for this dense index
   CaloTowerDetId id = m_limits->detIdFromDenseIndex(din);
   int ieta = id.ieta();
   int iphi = id.iphi();
   
   //use CT topology to get proper ieta for hcal
   int etaRing=m_limits->convertCTtoHcal(abs(ieta));
   int sign=(ieta>0)?(1):(-1);
   double eta1, eta2;
   if (abs(ieta)>m_limits->lastHERing()) {
     eta1 = theHFEtaBounds[etaRing-m_hcaltopo->firstHFRing()];
     eta2 = theHFEtaBounds[etaRing-m_hcaltopo->firstHFRing()+1];
   } else {
     eta1 = theHBHEEtaBounds[etaRing-1];
     eta2 = theHBHEEtaBounds[etaRing];
   }
   double eta = 0.5*(eta1+eta2);
   double deta = (eta2-eta1);  
 
   // in radians
   double dphi_nominal = 2.0*M_PI / m_hcaltopo->nPhiBins(1); // always the same
   double dphi_half = M_PI / m_hcaltopo->nPhiBins(etaRing); // half-width
   
   double phi_low = dphi_nominal*(iphi-1); // low-edge boundaries are constant...
   double phi = phi_low+dphi_half;
 
   double x,y,z,thickness;
   bool alongZ=true;
   if (abs(ieta)>m_limits->lastHERing()) { // forward
     z=HFz;
     double r=z/sinh(eta);
     x=r * cos(phi);
     y=r * sin(phi);
     thickness=HFthick/tanh(eta);
   } else if (abs(ieta)>m_limits->firstHERing()+1) { // EE-containing
     z=EEz;
     double r=z/sinh(eta);
     x=r * cos(phi);
     y=r * sin(phi);
     thickness=(HEz-EEz)/tanh(eta);
   } else { // EB-containing
     x=EBradius * cos(phi);
     y=EBradius * sin(phi);
     alongZ=false;
     z=EBradius * sinh(eta);
     thickness=(HOradius-EBradius) * cosh(eta);
   }
 
   z*=sign;
   GlobalPoint point(x,y,z);
 
   const double mysign ( !alongZ ? 1 : -1 ) ;
   std::vector<CCGFloat> hh ;
   hh.reserve(5) ;
   hh.push_back( deta/2 ) ;
   hh.push_back( dphi_half ) ;
   hh.push_back( mysign*thickness/2. ) ;
 
   hh.push_back( fabs( eta ) ) ;
   hh.push_back( fabs( z ) ) ;
 
   geom->newCell( point, point, point,
          CaloCellGeometry::getParmPtr( hh, 
                            geom->parMgr(), 
                            geom->parVecVec() ),
          id ) ;
 }