RPCGeometryParsFromDD Class Reference

#include <RPCGeometryParsFromDD.h>

Public Member Functions
void	build (const DDCompactView *cview, const MuonDDDConstants &muonConstants, RecoIdealGeometry &rgeo)
	RPCGeometryParsFromDD ()
	~RPCGeometryParsFromDD ()

Private Member Functions
void	buildGeometry (DDFilteredView &fview, const MuonDDDConstants &muonConstants, RecoIdealGeometry &rgeo)

Detailed Description

Build the RPCGeometry ftom the DDD description

Author

Port of: MuDDDRPCBuilder, MuonRPCGeometryBuilder (ORCA)

M. Maggi - INFN Bari

Definition at line 22 of file RPCGeometryParsFromDD.h.

Constructor & Destructor Documentation

RPCGeometryParsFromDD::RPCGeometryParsFromDD

(

)

Implementation of the RPC Geometry Builder from DDD

Author

Port of: MuDDDRPCBuilder (ORCA)

M. Maggi - INFN Bari

Definition at line 25 of file RPCGeometryParsFromDD.cc.

{ }

RPCGeometryParsFromDD::~RPCGeometryParsFromDD

(

)

Definition at line 28 of file RPCGeometryParsFromDD.cc.

{ }

Member Function Documentation

void RPCGeometryParsFromDD::build	(	const DDCompactView *	cview,
		const MuonDDDConstants &	muonConstants,
		RecoIdealGeometry &	rgeo
	)

Definition at line 31 of file RPCGeometryParsFromDD.cc.

References buildGeometry(), ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by RPCRecoIdealDBLoader::beginRun().

{
  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  )

private

Definition at line 44 of file RPCGeometryParsFromDD.cc.

References RPCNumberingScheme::baseNumberToUnitNumber(), gather_cfg::cout, DDfetch(), DDValue::doubles(), DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), RecoIdealGeometry::insert(), createfilelist::int, DDFilteredView::logicalPart(), DDName::name(), dataset::name, DDBase< N, C >::name(), DDFilteredView::nextSibling(), DDSolid::parameters(), DetId::rawId(), DDFilteredView::rotation(), DDLogicalPart::solid(), DDFilteredView::specifics(), AlCaHLTBitMon_QueryRunRegistry::string, DDFilteredView::translation(), ApeEstimator_cff::width, x, y, and z.

Referenced by build().

{
  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 is=specs.begin();is!=specs.end(); is++){
      if (DDfetch( *is, 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};
    DDTranslation tran = fview.translation();

    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);
  }

}