RPCGeometryParsFromDD.cc

Go to the documentation of this file.

#include "Geometry/RPCGeometryBuilder/src/RPCGeometryParsFromDD.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"

#include <CondFormats/GeometryObjects/interface/RecoIdealGeometry.h>
#include <DetectorDescription/Core/interface/DDFilter.h>
#include <DetectorDescription/Core/interface/DDFilteredView.h>
#include <DetectorDescription/Core/interface/DDSolid.h>

#include "Geometry/MuonNumbering/interface/MuonDDDNumbering.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/RPCNumberingScheme.h"

#include "DataFormats/GeometryVector/interface/Basic3DVector.h"

#include "CLHEP/Units/GlobalSystemOfUnits.h"

#include <iostream>
#include <algorithm>

RPCGeometryParsFromDD::RPCGeometryParsFromDD()
{ }

RPCGeometryParsFromDD::~RPCGeometryParsFromDD()
{ }

void RPCGeometryParsFromDD::build(const DDCompactView* cview,
                                  const MuonDDDConstants& muonConstants, RecoIdealGeometry& rgeo )
{
  const std::string attribute = "ReadOutName";
  const std::string value     = "MuonRPCHits";

  // Asking only for the MuonRPC's
  DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0.0)};
  DDFilteredView fview(*cview,filter);

  this->buildGeometry(fview, muonConstants, rgeo);
}

void RPCGeometryParsFromDD::buildGeometry(DDFilteredView& fview,
                                          const MuonDDDConstants& muonConstants, RecoIdealGeometry& rgeo)
{
  for ( bool doSubDets = fview.firstChild(); doSubDets==true; doSubDets = fview.nextSibling() ) {

    // Get the Base Muon Number
    MuonDDDNumbering mdddnum(muonConstants);
    MuonBaseNumber   mbn=mdddnum.geoHistoryToBaseNumber(fview.geoHistory());

    // Get the The Rpc det Id
    RPCNumberingScheme rpcnum(muonConstants);
    const int detid = rpcnum.baseNumberToUnitNumber(mbn);
    RPCDetId rpcid(detid);

    DDValue numbOfStrips("nStrips");

    std::vector<const DDsvalues_type* > specs(fview.specifics());
    int nStrips=0;
    for (auto & spec : specs){
      if (DDfetch( spec, numbOfStrips)){
        nStrips=int(numbOfStrips.doubles()[0]);
      }
    }
    if (nStrips == 0 ) std::cout <<"No strip found!!"<<std::endl;

    const std::vector<double> dpar = fview.logicalPart().solid().parameters();

    const std::string name=fview.logicalPart().name().name();
    const std::vector<std::string> strpars = {name};
    const DDTranslation& tran = fview.translation();

    const DDRotationMatrix& rota = fview.rotation();//.Inverse();
    DD3Vector x, y, z;
    rota.GetComponents(x,y,z);
    std::vector<double> pars;
    if (dpar.size()==3){
      const double width     = dpar[0];
      const double length    = dpar[1];
      const double thickness = dpar[2];
      pars = {width, length, thickness, numbOfStrips.doubles()[0]};
    }
    else{
      pars = {
        dpar[4] /*b/2*/, dpar[8] /*B/2*/, dpar[0] /*h/2*/,
        0.4,
        numbOfStrips.doubles()[0] /*h/2*/
      };
    }

    const std::vector<double> vtra = {tran.x(), tran.y(), tran.z()};
    const std::vector<double> vrot = {x.X(), x.Y(), x.Z(),
                                      y.X(), y.Y(), y.Z(),
                                      z.X(), z.Y(), z.Z()};
    rgeo.insert(rpcid.rawId(),vtra,vrot, pars,strpars);
  }

}