HcalDDDGeometryLoader Class Reference

#include <HcalDDDGeometryLoader.h>

Public Types
typedef CaloSubdetectorGeometry *	ReturnType

Public Member Functions
	HcalDDDGeometryLoader (const HcalDDDRecConstants *hcons)
ReturnType	load (const HcalTopology &topo, DetId::Detector, int)
ReturnType	load (const HcalTopology &topo)
	Load all of HCAL. More...
virtual	~HcalDDDGeometryLoader ()

Private Member Functions
void	fill (HcalSubdetector, HcalDDDGeometry *, CaloSubdetectorGeometry *)
	HcalDDDGeometryLoader ()
void	makeCell (const HcalDetId &, const HcalCellType &, double, double, CaloSubdetectorGeometry *geom) const

Private Attributes
HcalTopology *	dummyTopology_
const HcalDDDRecConstants *	hcalConstants

Detailed Description

Note

The Geometry as loaded from DDD

Author

S. Banerjee

Definition at line 20 of file HcalDDDGeometryLoader.h.

Member Typedef Documentation

typedef CaloSubdetectorGeometry* HcalDDDGeometryLoader::ReturnType

Definition at line 27 of file HcalDDDGeometryLoader.h.

Constructor & Destructor Documentation

HcalDDDGeometryLoader::HcalDDDGeometryLoader ( const HcalDDDRecConstants *  hcons )

explicit

Definition at line 16 of file HcalDDDGeometryLoader.cc.

  : hcalConstants(hcons) { }

HcalDDDGeometryLoader::~HcalDDDGeometryLoader ( )

virtual

Definition at line 19 of file HcalDDDGeometryLoader.cc.

                                              {
}

HcalDDDGeometryLoader::HcalDDDGeometryLoader ( )

private

Member Function Documentation

void HcalDDDGeometryLoader::fill ( HcalSubdetector  subdet, HcalDDDGeometry *  geometryDDD, CaloSubdetectorGeometry *  geom  )

private

helper functions to make all the ids and cells, and put them into the vectors and mpas passed in.

Definition at line 69 of file HcalDDDGeometryLoader.cc.

References HcalDDDRecConstants::HcalCellTypes(), hcalConstants, i, HcalDDDGeometry::insertCell(), relval_steps::k, GetRecoTauVFromDQM_MC_cff::kk, LogDebug, M_PI, makeCell(), convertSQLiteXML::ok, and phi.

Referenced by load().

                                                           {

  // start by making the new HcalDetIds
  std::vector<HcalCellType> hcalCells = hcalConstants->HcalCellTypes(subdet);
  geometryDDD->insertCell(hcalCells);
#ifdef DebugLog
  LogDebug("HCalGeom") << "HcalDDDGeometryLoader::fill gets " 
               << hcalCells.size() << " cells for subdetector " 
               << subdet;
#endif           
  // Make the new HcalDetIds and the cells

  double deg = M_PI/180.;
  std::vector<HcalDetId> hcalIds;
  for (unsigned int i=0; i<hcalCells.size(); i++) {
    int etaRing  = hcalCells[i].etaBin();
    int depthBin = hcalCells[i].depthSegment();
    int phiInc   = 4/hcalCells[i].nPhiModule();
    unsigned int iphi = 1;
    if (hcalCells[i].unitPhi() == 4) iphi = 3;
    double  dphi = (hcalCells[i].phiBinWidth())*deg;
    double   phi =-(hcalCells[i].phiOffset())*deg + 0.5*dphi;
    std::vector<int>missPlus  = hcalCells[i].missingPhiPlus();
    std::vector<int>missMinus = hcalCells[i].missingPhiMinus();
#ifdef DebugLog
    LogDebug("HCalGeom") << "HcalDDDGeometryLoader: Subdet " << subdet
                         << " eta " << etaRing << " depth " << depthBin
                         << " modules " << hcalCells[i].nPhiModule() << " "
                         << phiInc << " phi " << phi/deg << " " << dphi/deg
                         << " Missing " << missPlus.size() << "/"
                         << missMinus.size();
#endif
    for (int k = 0; k < hcalCells[i].nPhiBins(); k++) {
      bool ok = true;
      for (unsigned int kk = 0; kk < missPlus.size(); kk++)
        if (iphi == (unsigned int)(missPlus[kk])) ok = false;
      if (ok) {
#ifdef DebugLog
        LogDebug("HCalGeom") << "HcalDDDGeometryLoader::fill Cell " << i
                             << " eta " << etaRing << " phi " << iphi << "("
                             << phi/deg << ", " << dphi/deg << ") depth "
                             << depthBin;
#endif
        HcalDetId id(subdet, etaRing, iphi, depthBin);
        hcalIds.push_back(id);
        makeCell(id,hcalCells[i],phi,dphi,geom) ;
      }
      if (hcalCells[i].nHalves() > 1) {
        ok = true;
        for (unsigned int kk = 0; kk < missMinus.size(); kk++)
          if (iphi == (unsigned int)(missMinus[kk])) ok = false;
        if (ok) {
#ifdef DebugLog
          LogDebug("HCalGeom") << "HcalDDDGeometryLoader::fill Cell " << i
                               << " eta " << -etaRing << " phi " << iphi <<" ("
                               << phi/deg << ", " << dphi/deg << ") depth "
                               << depthBin;
#endif
          HcalDetId id(subdet, -etaRing, iphi, depthBin);
          hcalIds.push_back(id);
      makeCell(id,hcalCells[i],phi,dphi,geom) ;
        }
      }
      iphi += phiInc;
      phi  += dphi;
    }
  }
  
  edm::LogInfo("HCalGeom") << "Number of HCAL DetIds made for " << subdet
               << " is " << hcalIds.size();
}

HcalDDDGeometryLoader::ReturnType HcalDDDGeometryLoader::load	(	const HcalTopology &	topo,
		DetId::Detector	det,
		int	subdet
	)

Definition at line 23 of file HcalDDDGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), prof2calltree::count, fill(), relativeConstraints::geom, HcalBarrel, hcalConstants, HcalEndcap, HcalForward, HcalOuter, HcalDDDRecConstants::numberOfCells(), and CaloSubdetectorGeometry::parMgr().

Referenced by HcalDDDGeometryEP::produceIdeal().

                                                                                   {

  HcalSubdetector  hsub        = static_cast<HcalSubdetector>(subdet);


  HcalDDDGeometry* gDDD ( new HcalDDDGeometry(topo) );
  ReturnType geom ( gDDD );

  if ( geom->cornersMgr() == 0 ) {
     const unsigned int count (hcalConstants->numberOfCells(HcalBarrel ) +
                   hcalConstants->numberOfCells(HcalEndcap ) +
                   hcalConstants->numberOfCells(HcalForward) +
                   hcalConstants->numberOfCells(HcalOuter  ) );
     geom->allocateCorners( count ) ;
  }

  //  if( geom->cornersMgr() == 0 )  geom->allocateCorners( 2592 ) ;

  if ( geom->parMgr()     == 0 ) geom->allocatePar( 500, 3 ) ;

  fill (hsub, gDDD, geom );
  return geom ;
}

HcalDDDGeometryLoader::ReturnType HcalDDDGeometryLoader::load

(

const HcalTopology &

topo

)

Load all of HCAL.

Definition at line 48 of file HcalDDDGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), prof2calltree::count, fill(), relativeConstraints::geom, HcalBarrel, hcalConstants, HcalEndcap, HcalForward, HcalOuter, HcalDDDRecConstants::numberOfCells(), and CaloSubdetectorGeometry::parMgr().

                                                    {

  HcalDDDGeometry* gDDD ( new HcalDDDGeometry(topo) );
  ReturnType geom ( gDDD );

  if( geom->cornersMgr() == 0 ) {
    const unsigned int count (hcalConstants->numberOfCells(HcalBarrel ) +
                  hcalConstants->numberOfCells(HcalEndcap ) +
                  hcalConstants->numberOfCells(HcalForward) +
                  hcalConstants->numberOfCells(HcalOuter  ) );
    geom->allocateCorners( count ) ;
  }
  if( geom->parMgr()     == 0 ) geom->allocatePar( 500, 3 ) ;
  
  fill(HcalBarrel,  gDDD, geom); 
  fill(HcalEndcap,  gDDD, geom); 
  fill(HcalForward, gDDD, geom); 
  fill(HcalOuter,   gDDD, geom);
  return geom ;
}

void HcalDDDGeometryLoader::makeCell ( const HcalDetId &  detId, const HcalCellType &  hcalCell, double  phi, double  dphi, CaloSubdetectorGeometry *  geom  ) const

private

Definition at line 144 of file HcalDDDGeometryLoader.cc.

References funct::abs(), funct::cos(), HcalCellType::depthMax(), HcalCellType::depthMin(), HcalCellType::depthType(), eta, HcalCellType::etaMax(), HcalCellType::etaMin(), create_public_lumi_plots::exp, CaloCellGeometry::getParmPtr(), HcalCellType::halfSize(), HcalBarrel, HcalForward, HcalOuter, AnalysisDataFormats_SUSYBSMObjects::hp, GeomDetEnumerators::isBarrel(), LogDebug, CaloSubdetectorGeometry::newCell(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), point, alignCSCRings::r, diffTwoXMLs::r1, diffTwoXMLs::r2, jetcorrextractor::sign(), funct::sin(), HcalDetId::subdet(), funct::tan(), theta(), x, y, z, and HcalDetId::zside().

Referenced by fill().

                                                      {

  // the two eta boundaries of the cell
  double          eta1   = hcalCell.etaMin();
  double          eta2   = hcalCell.etaMax();
  HcalSubdetector subdet = detId.subdet();
  double          eta    = 0.5*(eta1+eta2) * detId.zside();
  double          deta   = (eta2-eta1);
  double          theta  = 2.0*atan(exp(-eta));

  // barrel vs forward
  bool rzType   = hcalCell.depthType();
  bool isBarrel = (subdet == HcalBarrel || subdet == HcalOuter);

  double          z, r, thickness;
#ifdef DebugLog
  double          r0, r1, r2;
#endif
  if (rzType) {
    r          = hcalCell.depthMin();
    if (isBarrel) {
      z         = r * sinh(eta); // sinh(eta) == 1/tan(theta)
      thickness = (hcalCell.depthMax() - r) * cosh(eta); // cosh(eta) == 1/sin(theta)
#ifdef DebugLog
      r1        = r;
      r2        = hcalCell.depthMax();
      r0        = 0.5*(r1+r2);
#endif
    } else {
      z         = r * sinh(eta2);
      thickness = 2. * hcalCell.halfSize();
      r         = z/sinh(std::abs(eta));
#ifdef DebugLog
      r0        = z/sinh(std::abs(eta));
      r1        = z/sinh(std::abs(eta)+0.5*deta);
      r2        = z/sinh(std::abs(eta)-0.5*deta);
#endif
    }
#ifdef DebugLog
    LogDebug("HCalGeom") << "HcalDDDGeometryLoader::makeCell SubDet " << subdet
             << " eta = " << eta << " theta = " << theta
             << " r = " << r << " thickness = " << thickness
                         << " r0-r2 (" << r0 << ":" << r1 << ":" << r2 << ")";
#endif
  } else {
    z          = hcalCell.depthMin();
    thickness  = hcalCell.depthMax() - z;
    z         *= detId.zside(); // get the sign right.
    r          = z * tan(theta);
    thickness /= std::abs(cos(theta));
#ifdef DebugLog
    r0         = z/sinh(std::abs(eta));
    r1         = z/sinh(std::abs(eta)+0.5*deta);
    r2         = z/sinh(std::abs(eta)-0.5*deta);
    LogDebug("HCalGeom") << "HcalDDDGeometryLoader::makeCell SubDet " << subdet
                         << " eta = " << eta << " theta = " << theta
                         << " z = " << z << " r = " << r << " thickness = "
                         << thickness << " r0-r2 (" << r0 << ":" << r1 << ":"
                         << r2 << ")";    
#endif
  }

  double x = r * cos(phi);
  double y = r * sin(phi);
  GlobalPoint point(x,y,z);

#ifdef DebugLog
  LogDebug("HCalGeom") << "HcalDDDGeometryLoader::makeCell for " << detId
               << " Point " << point << " deta = " << deta 
               << " dphi = " << dphi << " thickness = " << thickness
               << " isBarrel = " << isBarrel << " " << rzType;
#endif

  std::vector<CCGFloat> hp ;
  hp.reserve(3) ;
  
  if (subdet==HcalForward) {
    hp.push_back(deta/2.) ;
    hp.push_back(dphi/2.) ;
    hp.push_back(thickness/2.) ;
  } else { 
    const double sign ( isBarrel ? 1 : -1 ) ;
    hp.push_back(deta/2.) ;
    hp.push_back(dphi/2.) ;
    hp.push_back(sign*thickness/2.) ;
  }
  geom->newCell( point, point, point,
         CaloCellGeometry::getParmPtr( hp, 
                           geom->parMgr(), 
                           geom->parVecVec() ),
         detId ) ;
}

Member Data Documentation

HcalTopology* HcalDDDGeometryLoader::dummyTopology_

private

Definition at line 46 of file HcalDDDGeometryLoader.h.

const HcalDDDRecConstants* HcalDDDGeometryLoader::hcalConstants

private

Definition at line 44 of file HcalDDDGeometryLoader.h.

Referenced by fill(), and load().