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 21 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 30 of file ME0GeometryBuilderFromCondDB.cc.

{ }

Member Function Documentation

ME0Geometry * ME0GeometryBuilderFromCondDB::build

(

const RecoIdealGeometry &

rgeo

)

Definition at line 33 of file ME0GeometryBuilderFromCondDB.cc.

References RecoIdealGeometry::detIds(), geometry, hcalTTPDigis_cfi::id, GeomDetEnumerators::ME0, dataset::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;
    float nstrip = -999.;
    float npad = -999.;
    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;
}