DTGeometryParsFromDD Class Reference

#include <DTGeometryParsFromDD.h>

Public Types
enum	DTDetTag { DTChamberTag, DTSuperLayerTag, DTLayerTag }

Public Member Functions
void	build (const DDCompactView *cview, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig)
void	build (const cms::DDCompactView *cview, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rgeo)
	DTGeometryParsFromDD ()
	Constructor. More...
virtual	~DTGeometryParsFromDD ()
	Destructor. More...

Private Types
typedef std::pair< std::vector< double >, std::vector< double > >	PosRotPair

Private Member Functions
void	buildGeometry (DDFilteredView &fv, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
void	buildGeometry (cms::DDFilteredView &fv, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
std::vector< double >	extractParameters (DDFilteredView &fv) const
	get parameter also for boolean solid. More...
void	insertChamber (DDFilteredView &fv, const std::string &type, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
	create the chamber More...
void	insertChamber (cms::DDFilteredView &fv, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
	create the chamber More...
void	insertLayer (DDFilteredView &fv, const std::string &type, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
	create the layer More...
void	insertLayer (cms::DDFilteredView &fv, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
	create the layer More...
void	insertSuperLayer (DDFilteredView &fv, const std::string &type, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
	create the SL More...
void	insertSuperLayer (cms::DDFilteredView &fv, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rig) const
	create the SL More...
PosRotPair	plane (const DDFilteredView &fv) const
PosRotPair	plane (const cms::DDFilteredView &fv) const

Detailed Description

Build the RPCGeometry from the DDD and DD4hep description

DD4hep part added to the original old file (DD version) made by Stefano Lacaprara (INFN LNL)

Author

: Sergio Lo Meo (sergi.nosp@m.o.lo.nosp@m..meo@.nosp@m.cern.nosp@m..ch) Created: Tue, 26 Jan 2021

Definition at line 29 of file DTGeometryParsFromDD.h.

Member Typedef Documentation

PosRotPair

typedef std::pair<std::vector<double>, std::vector<double> > DTGeometryParsFromDD::PosRotPair

private

Definition at line 68 of file DTGeometryParsFromDD.h.

Member Enumeration Documentation

DTDetTag

enum DTGeometryParsFromDD::DTDetTag

Enumerator
DTChamberTag
DTSuperLayerTag
DTLayerTag

Definition at line 43 of file DTGeometryParsFromDD.h.

{ DTChamberTag, DTSuperLayerTag, DTLayerTag };

Constructor & Destructor Documentation

DTGeometryParsFromDD()

DTGeometryParsFromDD::DTGeometryParsFromDD

(

)

Constructor.

Definition at line 37 of file DTGeometryParsFromDD.cc.

{}

~DTGeometryParsFromDD()

DTGeometryParsFromDD::~DTGeometryParsFromDD ( )

virtual

Destructor.

Definition at line 39 of file DTGeometryParsFromDD.cc.

{}

Member Function Documentation

build() [1/2]

void DTGeometryParsFromDD::build	(	const DDCompactView *	cview,
		const MuonGeometryConstants &	muonConstants,
		RecoIdealGeometry &	rig
	)

Definition at line 42 of file DTGeometryParsFromDD.cc.

References ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by DTRecoIdealDBLoader::beginRun().

                                                         {
  //  cout << "DTGeometryParsFromDD::build" << endl;
  //   static const string t0 = "DTGeometryParsFromDD::build";
  //   TimeMe timer(t0,true);

  std::string attribute = "MuStructure";
  std::string value = "MuonBarrelDT";

  // Asking only for the Muon DTs
  DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0.0)};
  DDFilteredView fview(*cview, filter);
  buildGeometry(fview, muonConstants, rig);
  //cout << "RecoIdealGeometry " << rig.size() << endl;
}

build() [2/2]

void DTGeometryParsFromDD::build	(	const cms::DDCompactView *	cview,
		const MuonGeometryConstants &	muonConstants,
		RecoIdealGeometry &	rgeo
	)

Definition at line 61 of file DTGeometryParsFromDD.cc.

References ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                          {
  const std::string attribute = "MuStructure";
  const std::string value = "MuonBarrelDT";
  const cms::DDFilter filter(attribute, value);
  cms::DDFilteredView fview(*cview, filter);
  buildGeometry(fview, muonConstants, rgeo);
}

buildGeometry() [1/2]

void DTGeometryParsFromDD::buildGeometry ( DDFilteredView &  fv, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

Definition at line 72 of file DTGeometryParsFromDD.cc.

References DDfetch(), DDFilteredView::firstChild(), DDFilteredView::mergedSpecifics(), DDFilteredView::nextSibling(), submitPVValidationJobs::params, DDFilteredView::parent(), and heppy_batch::val.

                                                                       {
  // static const string t0 = "DTGeometryParsFromDD::buildGeometry";
  // TimeMe timer(t0,true);

  edm::LogVerbatim("DTGeometryParsFromDD") << "(0) DTGeometryParsFromDD - DDD ";

  bool doChamber = fv.firstChild();

  // Loop on chambers
  int ChamCounter = 0;
  while (doChamber) {
    ChamCounter++;
    DDValue val("Type");
    const DDsvalues_type params(fv.mergedSpecifics());
    string type;
    if (DDfetch(&params, val))
      type = val.strings()[0];
    // FIXME
    val = DDValue("FEPos");
    string FEPos;
    if (DDfetch(&params, val))
      FEPos = val.strings()[0];
    insertChamber(fv, type, muonConstants, rig);

    // Loop on SLs
    bool doSL = fv.firstChild();
    int SLCounter = 0;
    while (doSL) {
      SLCounter++;
      insertSuperLayer(fv, type, muonConstants, rig);

      bool doL = fv.firstChild();
      int LCounter = 0;
      // Loop on SLs
      while (doL) {
        LCounter++;
        insertLayer(fv, type, muonConstants, rig);

        // fv.parent();
        doL = fv.nextSibling();  // go to next layer
      }                          // layers

      fv.parent();
      doSL = fv.nextSibling();  // go to next SL
    }                           // sls

    fv.parent();
    doChamber = fv.nextSibling();  // go to next chamber
  }                                // chambers
}

buildGeometry() [2/2]

void DTGeometryParsFromDD::buildGeometry ( cms::DDFilteredView &  fv, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

Definition at line 264 of file DTGeometryParsFromDD.cc.

References cms::DDFilteredView::down(), cms::DDFilteredView::firstChild(), cms::DDFilteredView::nextSibling(), cms::DDFilteredView::parent(), and cms::DDFilteredView::sibling().

                                                                       {
  edm::LogVerbatim("DTGeometryParsFromDD") << "(0) DTGeometryParsFromDD - DD4hep ";

  bool doChamber = fv.firstChild();

  while (doChamber) {
    insertChamber(fv, muonConstants, rig);

    bool doSL = fv.nextSibling();
    while (doSL) {
      insertSuperLayer(fv, muonConstants, rig);

      fv.down();
      bool doLayers = fv.sibling();
      while (doLayers) {
        insertLayer(fv, muonConstants, rig);

        doLayers = fv.sibling();
      }

      doSL = fv.nextSibling();
    }

    fv.parent();
    doChamber = fv.firstChild();
  }
}

extractParameters()

vector< double > DTGeometryParsFromDD::extractParameters ( DDFilteredView &  fv ) const

private

get parameter also for boolean solid.

Definition at line 213 of file DTGeometryParsFromDD.cc.

References A, cms::cuda::bs, ddbox, DDFilteredView::logicalPart(), DDSolid::parameters(), DDSolid::shape(), and DDLogicalPart::solid().

                                                                               {
  vector<double> par;
  if (fv.logicalPart().solid().shape() != DDSolidShape::ddbox) {
    DDBooleanSolid bs(fv.logicalPart().solid());
    DDSolid A = bs.solidA();
    while (A.shape() != DDSolidShape::ddbox) {
      DDBooleanSolid bs(A);
      A = bs.solidA();
    }
    par = A.parameters();
  } else {
    par = fv.logicalPart().solid().parameters();
  }
  return par;
}

insertChamber() [1/2]

void DTGeometryParsFromDD::insertChamber ( DDFilteredView &  fv, const std::string &  type, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

create the chamber

insertChamber() [2/2]

void DTGeometryParsFromDD::insertChamber ( cms::DDFilteredView &  fv, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

create the chamber

Definition at line 327 of file DTGeometryParsFromDD.cc.

References DTNumberingScheme::baseNumberToUnitNumber(), MuonGeometryNumbering::geoHistoryToBaseNumber(), cms::DDFilteredView::history(), RecoIdealGeometry::insert(), and cms::DDFilteredView::parameters().

                                                                       {
  MuonGeometryNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.history()));
  DTChamberId detId(rawid);

  vector<double> par_temp = fv.parameters();
  vector<double> par(4);
  par[0] = DTChamberTag;  //DTChamberTag is the ID of a Chamber
  par[1] = par_temp[0] / dd4hep::mm;
  par[2] = par_temp[1] / dd4hep::mm;
  par[3] = par_temp[2] / dd4hep::mm;

  edm::LogVerbatim("DTGeometryParsFromDD")
      << "(1) DD4hep, Chamber DetID " << rawid << " " << par[0] << " " << par[1] << " " << par[2] << " " << par[3];
  PosRotPair posRot(plane(fv));
  rig.insert(rawid, posRot.first, posRot.second, par);
}

insertLayer() [1/2]

void DTGeometryParsFromDD::insertLayer ( DDFilteredView &  fv, const std::string &  type, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

create the layer

Definition at line 173 of file DTGeometryParsFromDD.cc.

References DDFilteredView::copyno(), DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonGeometryNumbering::geoHistoryToBaseNumber(), DTNumberingScheme::getDetId(), RecoIdealGeometry::insert(), DDFilteredView::nextSibling(), DDFilteredView::parent(), and findQualityFiles::size.

                                                                     {
  MuonGeometryNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.getDetId(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
  DTLayerId layId(rawid);
  //cout << "inserting Layer " << layId << endl;

  // Layer specific parameter (size)
  vector<double> par;
  par.emplace_back(DTLayerTag);
  vector<double> size = extractParameters(fv);
  par.insert(par.end(), size.begin(), size.end());

  // Loop on wires
  bool doWire = fv.firstChild();
  int WCounter = 0;
  int firstWire = fv.copyno();
  //float wireLength = par[1]/cm;
  while (doWire) {
    WCounter++;
    doWire = fv.nextSibling();  // next wire
  }
  vector<double> sensSize = extractParameters(fv);
  //int lastWire=fv.copyno();
  par.emplace_back(firstWire);
  par.emplace_back(WCounter);
  par.emplace_back(sensSize[1]);
  fv.parent();

  PosRotPair posRot(plane(fv));

  edm::LogVerbatim("DTGeometryParsFromDD")
      << "(3) DDD, Layer DetID " << rawid << " " << par[0] << " " << par[1] << " " << par[2] << " " << par[3] << " "
      << par[4] << " " << par[5] << " " << par[6];
  rig.insert(layId, posRot.first, posRot.second, par);
}

insertLayer() [2/2]

void DTGeometryParsFromDD::insertLayer ( cms::DDFilteredView &  fv, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

create the layer

Definition at line 369 of file DTGeometryParsFromDD.cc.

References DTNumberingScheme::baseNumberToUnitNumber(), cms::DDFilteredView::checkChild(), cms::DDFilteredView::down(), MuonGeometryNumbering::geoHistoryToBaseNumber(), cms::DDFilteredView::history(), RecoIdealGeometry::insert(), cms::DDFilteredView::parameters(), cms::DDFilteredView::sibling(), cms::DDFilteredView::up(), and cms::DDFilteredView::volume().

                                                                     {
  MuonGeometryNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.history()));
  DTLayerId layId(rawid);

  vector<double> par_temp = fv.parameters();
  vector<double> par(4);
  par[0] = DTLayerTag;  //DTLayerTag is the ID of a Layer
  par[1] = par_temp[0] / dd4hep::mm;
  par[2] = par_temp[1] / dd4hep::mm;
  par[3] = par_temp[2] / dd4hep::mm;

  fv.down();
  bool doWire = fv.sibling();
  int firstWire = fv.volume()->GetNumber();
  auto const& wpar = fv.parameters();
  float wireLength = wpar[1] / dd4hep::mm;

  int WCounter = 0;
  while (doWire) {
    doWire = fv.checkChild();
    WCounter++;
  }

  par.emplace_back(firstWire);
  par.emplace_back(WCounter);
  par.emplace_back(wireLength);

  fv.up();

  PosRotPair posRot(plane(fv));

  edm::LogVerbatim("DTGeometryParsFromDD")
      << "(3) DD4hep, Layer DetID " << rawid << " " << par[0] << " " << par[1] << " " << par[2] << " " << par[3] << " "
      << par[4] << " " << par[5] << " " << par[6];
  rig.insert(layId, posRot.first, posRot.second, par);
}

insertSuperLayer() [1/2]

void DTGeometryParsFromDD::insertSuperLayer ( DDFilteredView &  fv, const std::string &  type, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

create the SL

Definition at line 151 of file DTGeometryParsFromDD.cc.

References DDFilteredView::geoHistory(), MuonGeometryNumbering::geoHistoryToBaseNumber(), DTNumberingScheme::getDetId(), RecoIdealGeometry::insert(), and findQualityFiles::size.

                                                                          {
  MuonGeometryNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.getDetId(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
  DTSuperLayerId slId(rawid);
  //cout << "inserting SuperLayer " << slId << endl;

  // Slayer specific parameter (size)
  vector<double> par;
  par.emplace_back(DTSuperLayerTag);
  vector<double> size = extractParameters(fv);
  par.insert(par.end(), size.begin(), size.end());
  edm::LogVerbatim("DTGeometryParsFromDD")
      << "(2) DDD, Super Layer DetID " << rawid << " " << par[0] << " " << par[1] << " " << par[2] << " " << par[3];
  // Ok this is the slayer position...
  PosRotPair posRot(plane(fv));
  rig.insert(slId, posRot.first, posRot.second, par);
}

insertSuperLayer() [2/2]

void DTGeometryParsFromDD::insertSuperLayer ( cms::DDFilteredView &  fv, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rig  ) const

private

create the SL

Definition at line 348 of file DTGeometryParsFromDD.cc.

References DTNumberingScheme::baseNumberToUnitNumber(), MuonGeometryNumbering::geoHistoryToBaseNumber(), cms::DDFilteredView::history(), RecoIdealGeometry::insert(), and cms::DDFilteredView::parameters().

                                                                          {
  MuonGeometryNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.history()));
  DTSuperLayerId slId(rawid);

  vector<double> par_temp = fv.parameters();
  vector<double> par(4);
  par[0] = DTSuperLayerTag;  //DTSuperLayerTag is the ID of a SuperLayer
  par[1] = par_temp[0] / dd4hep::mm;
  par[2] = par_temp[1] / dd4hep::mm;
  par[3] = par_temp[2] / dd4hep::mm;

  edm::LogVerbatim("DTGeometryParsFromDD")
      << "(2) DD4hep, Super Layer DetID " << rawid << " " << par[0] << " " << par[1] << " " << par[2] << " " << par[3];
  PosRotPair posRot(plane(fv));
  rig.insert(slId, posRot.first, posRot.second, par);
}

plane() [1/2]

DTGeometryParsFromDD::PosRotPair DTGeometryParsFromDD::plane ( const DDFilteredView &  fv ) const

private

Definition at line 229 of file DTGeometryParsFromDD.cc.

References angle_units::operators::convertMmToCm(), idealTransformation::rotation, DDFilteredView::rotation(), DDFilteredView::translation(), and x.

                                                                                         {
  // extract the position
  const DDTranslation& trans(fv.translation());

  std::vector<double> gtran(3);
  gtran[0] = convertMmToCm(trans.x());
  gtran[1] = convertMmToCm(trans.y());
  gtran[2] = convertMmToCm(trans.z());

  edm::LogVerbatim("DTGeometryParsFromDD") << "(4) DDD, Position "
                                           << " " << gtran[0] << " " << gtran[1] << " " << gtran[2];
  // now the rotation
  //     'active' and 'passive' rotations are inverse to each other
  const DDRotationMatrix& rotation = fv.rotation();  //REMOVED .Inverse();
  DD3Vector x, y, z;
  rotation.GetComponents(x, y, z);

  std::vector<double> grmat(9);
  grmat[0] = x.X();
  grmat[1] = x.Y();
  grmat[2] = x.Z();

  grmat[3] = y.X();
  grmat[4] = y.Y();
  grmat[5] = y.Z();

  grmat[6] = z.X();
  grmat[7] = z.Y();
  grmat[8] = z.Z();

  return pair<std::vector<double>, std::vector<double> >(gtran, grmat);
}

plane() [2/2]

DTGeometryParsFromDD::PosRotPair DTGeometryParsFromDD::plane ( const cms::DDFilteredView &  fv ) const

private

Definition at line 294 of file DTGeometryParsFromDD.cc.

References idealTransformation::rotation, cms::DDFilteredView::rotation(), cms::DDFilteredView::trans(), and x.

                                                                                            {
  const Double_t* tr = fv.trans();

  std::vector<double> gtran(3);

  gtran[0] = tr[0] / dd4hep::cm;
  gtran[1] = tr[1] / dd4hep::cm;
  gtran[2] = tr[2] / dd4hep::cm;

  edm::LogVerbatim("DTGeometryParsFromDD") << "(4) DD4hep, Position "
                                           << " " << gtran[0] << " " << gtran[1] << " " << gtran[2];

  DDRotationMatrix rotation = fv.rotation();
  DD3Vector x, y, z;
  rotation.GetComponents(x, y, z);

  std::vector<double> grmat(9);

  grmat[0] = x.X();
  grmat[1] = x.Y();
  grmat[2] = x.Z();

  grmat[3] = y.X();
  grmat[4] = y.Y();
  grmat[5] = y.Z();

  grmat[6] = z.X();
  grmat[7] = z.Y();
  grmat[8] = z.Z();

  return pair<std::vector<double>, std::vector<double> >(gtran, grmat);
}