EcalTBHodoscopeGeometryLoaderFromDDD Class Reference

#include <EcalTBHodoscopeGeometryLoaderFromDDD.h>

Public Member Functions
	EcalTBHodoscopeGeometryLoaderFromDDD ()
std::unique_ptr < CaloSubdetectorGeometry >	load (const DDCompactView *cpv)
virtual	~EcalTBHodoscopeGeometryLoaderFromDDD ()

Private Member Functions
std::string	getDDDString (std::string s, DDFilteredView *fv)
DDFilter *	getDDFilter ()
unsigned int	getDetIdForDDDNode (const DDFilteredView &fv)
void	makeGeometry (const DDCompactView *cpv, CaloSubdetectorGeometry *ebg)

Private Attributes
EcalHodoscopeNumberingScheme	_scheme

Detailed Description

Author

P. Meridiani - INFN Roma 1

Definition at line 22 of file EcalTBHodoscopeGeometryLoaderFromDDD.h.

Constructor & Destructor Documentation

EcalTBHodoscopeGeometryLoaderFromDDD::EcalTBHodoscopeGeometryLoaderFromDDD ( )

inline

Definition at line 24 of file EcalTBHodoscopeGeometryLoaderFromDDD.h.

{};

virtual EcalTBHodoscopeGeometryLoaderFromDDD::~EcalTBHodoscopeGeometryLoaderFromDDD ( )

inlinevirtual

Definition at line 26 of file EcalTBHodoscopeGeometryLoaderFromDDD.h.

{};

Member Function Documentation

std::string EcalTBHodoscopeGeometryLoaderFromDDD::getDDDString ( std::string  s, DDFilteredView *  fv  )

private

Referenced by makeGeometry().

DDFilter * EcalTBHodoscopeGeometryLoaderFromDDD::getDDFilter ( )

private

Definition at line 113 of file EcalTBHodoscopeGeometryLoaderFromDDD.cc.

Referenced by makeGeometry().

                                                            {
  return new DDAndFilter<DDSpecificsMatchesValueFilter, DDSpecificsMatchesValueFilter>(
      DDSpecificsMatchesValueFilter{DDValue("SensitiveDetector", "EcalTBH4BeamDetector", 0)},
      DDSpecificsMatchesValueFilter{DDValue("ReadOutName", "EcalTBH4BeamHits", 0)});
}

unsigned int EcalTBHodoscopeGeometryLoaderFromDDD::getDetIdForDDDNode ( const DDFilteredView &  fv )

private

Definition at line 96 of file EcalTBHodoscopeGeometryLoaderFromDDD.cc.

References _scheme, EcalBaseNumber::addLevel(), cms::cuda::assert(), DDFilteredView::geoHistory(), EcalHodoscopeNumberingScheme::getUnitID(), mps_fire::i, and parents.

Referenced by makeGeometry().

                                                                                              {
  // perform some consistency checks
  // get the parents and grandparents of this node
  DDGeoHistory parents = fv.geoHistory();

  assert(parents.size() >= 3);

  EcalBaseNumber baseNumber;
  //baseNumber.setSize(parents.size());

  for (unsigned int i = 1; i <= parents.size(); i++) {
    baseNumber.addLevel(parents[parents.size() - i].logicalPart().name().name(), parents[parents.size() - i].copyno());
  }

  return _scheme.getUnitID(baseNumber);
}

std::unique_ptr< CaloSubdetectorGeometry > EcalTBHodoscopeGeometryLoaderFromDDD::load

(

const DDCompactView *

cpv

)

Definition at line 21 of file EcalTBHodoscopeGeometryLoaderFromDDD.cc.

References gather_cfg::cout, and makeGeometry().

Referenced by EcalTBHodoscopeGeometryEP::produce().

                                                                                                          {
  std::cout << "[EcalTBHodoscopeGeometryLoaderFromDDD]:: start the construction of EcalTBHodoscope" << std::endl;

  std::unique_ptr<CaloSubdetectorGeometry> ebg(new EcalTBHodoscopeGeometry());

  makeGeometry(cpv, ebg.get());

  std::cout << "[EcalTBHodoscopeGeometryLoaderFromDDD]:: Returning EcalTBHodoscopeGeometry" << std::endl;

  return ebg;
}

void EcalTBHodoscopeGeometryLoaderFromDDD::makeGeometry ( const DDCompactView *  cpv, CaloSubdetectorGeometry *  ebg  )

private

Definition at line 33 of file EcalTBHodoscopeGeometryLoaderFromDDD.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), dqmiodatasetharvest::ctr, ddbox, Exception, alcazmumu_cfi::filter, DDFilteredView::firstChild(), getDDDString(), getDDFilter(), getDetIdForDDDNode(), CaloCellGeometry::getParmPtr(), mps_fire::i, CaloCellGeometry::k_ScaleFromDDDtoGeant, EBDetId::kSizeForDenseIndexing, DDFilteredView::logicalPart(), CaloSubdetectorGeometry::newCell(), DDFilteredView::nextSibling(), DDSolid::parameters(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), MetAnalyzer::pv(), DDFilteredView::rotation(), DDSolid::shape(), DDLogicalPart::solid(), AlCaHLTBitMon_QueryRunRegistry::string, DDFilteredView::translation(), x, y, and z.

Referenced by load().

                                                                                                              {
  if (ebg->cornersMgr() == nullptr)
    ebg->allocateCorners(EBDetId::kSizeForDenseIndexing);
  if (ebg->parMgr() == nullptr)
    ebg->allocatePar(10, 3);

  std::unique_ptr<DDFilter> filter{getDDFilter()};

  DDFilteredView fv(*cpv, *filter);

  bool doSubDets;
  for (doSubDets = fv.firstChild(); doSubDets; doSubDets = fv.nextSibling()) {
#if 0
      std::string answer = getDDDString("ReadOutName",&fv);
      if (answer != "EcalTBH4BeamHits")
        continue;
#endif

    const DDSolid& solid = fv.logicalPart().solid();

    if (solid.shape() != DDSolidShape::ddbox) {
      throw cms::Exception("DDException")
          << std::string(__FILE__)
          << "\n CaloGeometryEcalTBHodoscope::upDate(...): currently only box fiber shapes supported ";
      edm::LogWarning("EcalTBHodoscopeGeometry") << "Wrong shape for sensitive volume!" << solid;
    }

    std::vector<double> pv = solid.parameters();

    // use preshower strip as box in space representation

    // rotate the box and then move it
    DD3Vector x, y, z;
    fv.rotation().GetComponents(x, y, z);
    CLHEP::Hep3Vector hx(x.X(), x.Y(), x.Z());
    CLHEP::Hep3Vector hy(y.X(), y.Y(), y.Z());
    CLHEP::Hep3Vector hz(z.X(), z.Y(), z.Z());
    CLHEP::HepRotation hrot(hx, hy, hz);
    CLHEP::Hep3Vector htran(fv.translation().X(), fv.translation().Y(), fv.translation().Z());

    const HepGeom::Transform3D ht3d(hrot,  // only scale translation
                                    CaloCellGeometry::k_ScaleFromDDDtoGeant * htran);

    const HepGeom::Point3D<float> ctr(ht3d * HepGeom::Point3D<float>(0, 0, 0));

    const GlobalPoint refPoint(ctr.x(), ctr.y(), ctr.z());

    std::vector<CCGFloat> vv;
    vv.reserve(pv.size() + 1);
    for (unsigned int i(0); i != pv.size(); ++i) {
      vv.emplace_back(CaloCellGeometry::k_ScaleFromDDDtoGeant * pv[i]);
    }
    vv.emplace_back(0.);  // tilt=0 here
    const CCGFloat* pP(CaloCellGeometry::getParmPtr(vv, ebg->parMgr(), ebg->parVecVec()));

    const DetId detId(getDetIdForDDDNode(fv));

    //Adding cell to the Geometry

    ebg->newCell(refPoint, refPoint, refPoint, pP, detId);
  }  // loop over all children
}

Member Data Documentation

EcalHodoscopeNumberingScheme EcalTBHodoscopeGeometryLoaderFromDDD::_scheme

private

Definition at line 39 of file EcalTBHodoscopeGeometryLoaderFromDDD.h.

Referenced by getDetIdForDDDNode().