ME0GeometryBuilderFromCondDB Class Reference

#include <ME0GeometryBuilderFromCondDB.h>

Public Member Functions
ME0Geometry *	build (const RecoIdealGeometry &rgeo)
	ME0GeometryBuilderFromCondDB ()
	~ME0GeometryBuilderFromCondDB ()

Detailed Description

Build the ME0Geometry from the DDD description stored in Condition DB

Author

M. Maggi - INFN Bari

Definition at line 20 of file ME0GeometryBuilderFromCondDB.h.

Constructor & Destructor Documentation

ME0GeometryBuilderFromCondDB::ME0GeometryBuilderFromCondDB

(

)

Implementation of the ME0 Geometry Builder from DDD stored in CondDB

Author

M. Maggi - INFN Bari

Definition at line 27 of file ME0GeometryBuilderFromCondDB.cc.

{}

ME0GeometryBuilderFromCondDB::~ME0GeometryBuilderFromCondDB

(

)

Definition at line 29 of file ME0GeometryBuilderFromCondDB.cc.

{}

Member Function Documentation

ME0Geometry * ME0GeometryBuilderFromCondDB::build

(

const RecoIdealGeometry &

rgeo

)

Definition at line 31 of file ME0GeometryBuilderFromCondDB.cc.

References RecoIdealGeometry::detIds(), geometry, triggerObjects_cff::id, GeomDetEnumerators::ME0, Skims_PA_cff::name, gen::npad, makeMuonMisalignmentScenario::rot, TkRotation< T >::rotateAxes(), RecoIdealGeometry::rotStart(), RecoIdealGeometry::shapeStart(), AlCaHLTBitMon_QueryRunRegistry::string, RecoIdealGeometry::strStart(), and RecoIdealGeometry::tranStart().

Referenced by ME0GeometryESModule::produce().

                                                                              {
  const std::vector<DetId>& detids(rgeo.detIds());
  ME0Geometry* geometry = new ME0Geometry();

  std::string name;
  std::vector<double>::const_iterator tranStart;
  std::vector<double>::const_iterator shapeStart;
  std::vector<double>::const_iterator rotStart;
  std::vector<std::string>::const_iterator strStart;

  for (unsigned int id = 0; id < detids.size(); ++id) {
    ME0DetId me0id(detids[id]);

    tranStart = rgeo.tranStart(id);
    shapeStart = rgeo.shapeStart(id);
    rotStart = rgeo.rotStart(id);
    strStart = rgeo.strStart(id);
    name = *(strStart);

    Surface::PositionType pos(*(tranStart) / cm, *(tranStart + 1) / cm, *(tranStart + 2) / cm);
    // CLHEP way
    Surface::RotationType rot(*(rotStart + 0),
                              *(rotStart + 1),
                              *(rotStart + 2),
                              *(rotStart + 3),
                              *(rotStart + 4),
                              *(rotStart + 5),
                              *(rotStart + 6),
                              *(rotStart + 7),
                              *(rotStart + 8));

    Bounds* bounds = nullptr;
    float be = *(shapeStart + 0) / cm;
    float te = *(shapeStart + 1) / cm;
    float ap = *(shapeStart + 2) / cm;
    float ti = *(shapeStart + 3) / cm;
    float nstrip = *(shapeStart + 4);
    float npad = *(shapeStart + 5);
    //  TrapezoidalPlaneBounds*
    bounds = new TrapezoidalPlaneBounds(be, te, ap, ti);

    std::vector<float> pars;
    pars.emplace_back(be);  //b/2;
    pars.emplace_back(te);  //B/2;
    pars.emplace_back(ap);  //h/2;
    pars.emplace_back(nstrip);
    pars.emplace_back(npad);

    ME0EtaPartitionSpecs* e_p_specs = new ME0EtaPartitionSpecs(GeomDetEnumerators::ME0, name, pars);

    //Change of axes for the forward
    Basic3DVector<float> newX(1., 0., 0.);
    Basic3DVector<float> newY(0., 0., 1.);
    //      if (tran[2] > 0. )
    newY *= -1;
    Basic3DVector<float> newZ(0., 1., 0.);
    rot.rotateAxes(newX, newY, newZ);

    BoundPlane* bp = new BoundPlane(pos, rot, bounds);
    ReferenceCountingPointer<BoundPlane> surf(bp);
    ME0EtaPartition* mep = new ME0EtaPartition(me0id, surf, e_p_specs);
    geometry->add(mep);
  }
  return geometry;
}