DTGeometryBuilderFromDDD Class Reference

#include <DTGeometryBuilderFromDDD.h>

Public Member Functions
void	build (boost::shared_ptr< DTGeometry > theGeometry, const DDCompactView *cview, const MuonDDDConstants &muonConstants)
	DTGeometryBuilderFromDDD ()
	Constructor. More...
virtual	~DTGeometryBuilderFromDDD ()
	Destructor. More...

Private Types
typedef ReferenceCountingPointer < Plane >	RCPPlane

Private Member Functions
DTChamber *	buildChamber (DDFilteredView &fv, const std::string &type, const MuonDDDConstants &muonConstants) const
	create the chamber More...
void	buildGeometry (boost::shared_ptr< DTGeometry > theGeometry, DDFilteredView &fv, const MuonDDDConstants &muonConstants) const
DTLayer *	buildLayer (DDFilteredView &fv, DTSuperLayer *sl, const std::string &type, const MuonDDDConstants &muonConstants) const
	create the layer More...
DTSuperLayer *	buildSuperLayer (DDFilteredView &fv, DTChamber *chamber, const std::string &type, const MuonDDDConstants &muonConstants) const
	create the SL More...
std::vector< double >	extractParameters (DDFilteredView &fv) const
	get parameter also for boolean solid. More...
RCPPlane	plane (const DDFilteredView &fv, Bounds *bounds) const

Detailed Description

Build the DTGeometry from the DDD description.

Author

N. Amapane - CERN.

Port of: MuBarDDDGeomBuilder, MuBarDetBuilder (ORCA) by S. Lacaprara, M. Case

Definition at line 25 of file DTGeometryBuilderFromDDD.h.

Member Typedef Documentation

typedef ReferenceCountingPointer<Plane> DTGeometryBuilderFromDDD::RCPPlane

private

Definition at line 59 of file DTGeometryBuilderFromDDD.h.

Constructor & Destructor Documentation

DTGeometryBuilderFromDDD::DTGeometryBuilderFromDDD

(

)

Constructor.

Definition at line 31 of file DTGeometryBuilderFromDDD.cc.

{}

DTGeometryBuilderFromDDD::~DTGeometryBuilderFromDDD ( )

virtual

Destructor.

Definition at line 33 of file DTGeometryBuilderFromDDD.cc.

{}

Member Function Documentation

void DTGeometryBuilderFromDDD::build	(	boost::shared_ptr< DTGeometry >	theGeometry,
		const DDCompactView *	cview,
		const MuonDDDConstants &	muonConstants
	)

Definition at line 36 of file DTGeometryBuilderFromDDD.cc.

References DDFilteredView::addFilter(), AND, alcazmumu_cfi::filter, matches, DDSpecificsFilter::setCriteria(), AlCaHLTBitMon_QueryRunRegistry::string, and relativeConstraints::value.

Referenced by AlignmentMonitorAsAnalyzer::analyze(), AlignmentProducer::createGeometries_(), DTGeometryESModule::geometryCallback_(), MuonAlignmentInputMethod::idealDTGeometry(), MisalignedMuonESProducer::produce(), and MuonAlignmentOutputXML::write().

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

  std::string attribute = "MuStructure"; 
  std::string value     = "MuonBarrelDT";
  DDValue val(attribute, value, 0.0);

  // Asking only for the Muon DTs
  DDSpecificsFilter filter;
  filter.setCriteria(val,  // name & value of a variable 
             DDCompOp::matches,
             DDLogOp::AND, 
             true, // compare strings otherwise doubles
             true  // use merged-specifics or simple-specifics
             );
  DDFilteredView fview(*cview);
  fview.addFilter(filter);
  buildGeometry(theGeometry, fview, muonConstants);
}

DTChamber * DTGeometryBuilderFromDDD::buildChamber ( DDFilteredView &  fv, const std::string &  type, const MuonDDDConstants &  muonConstants  ) const

private

create the chamber

SL the definition of length, width, thickness depends on the local reference frame of the Det

Definition at line 115 of file DTGeometryBuilderFromDDD.cc.

References DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), DTNumberingScheme::getDetId(), and create_public_lumi_plots::width.

                                                                                                                   {
  MuonDDDNumbering mdddnum (muonConstants);
  DTNumberingScheme dtnum (muonConstants);
  int rawid = dtnum.getDetId(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
  DTChamberId detId(rawid);  

  // Chamber specific parameter (size) 
  // FIXME: some trouble for boolean solids?
  vector<double> par = extractParameters(fv);

  float width = par[0]/cm;     // r-phi  dimension - different in different chambers
  float length = par[1]/cm;    // z      dimension - constant 125.55 cm
  float thickness = par[2]/cm; // radial thickness - almost constant about 18 cm

  // width is along local X
  // length is along local Y
  // thickness is long local Z

  RCPPlane surf(plane(fv, new RectangularPlaneBounds(width, length, thickness) ));

  DTChamber* chamber = new DTChamber(detId, surf);

  return chamber;
}

void DTGeometryBuilderFromDDD::buildGeometry ( boost::shared_ptr< DTGeometry >  theGeometry, DDFilteredView &  fv, const MuonDDDConstants &  muonConstants  ) const

private

Definition at line 61 of file DTGeometryBuilderFromDDD.cc.

References DDfetch(), DDFilteredView::firstChild(), DDFilteredView::mergedSpecifics(), DDFilteredView::nextSibling(), DDFilteredView::parent(), and DDValue::strings().

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

  //DTGeometry* theGeometry = new DTGeometry;

  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];
    DTChamber* chamber = buildChamber(fv,type, muonConstants);

    // Loop on SLs
    bool doSL = fv.firstChild();
    int SLCounter=0;
    while (doSL) {
      SLCounter++;
      DTSuperLayer* sl = buildSuperLayer(fv, chamber, type, muonConstants);
      theGeometry->add(sl);

      bool doL = fv.firstChild();
      int LCounter=0;
      // Loop on SLs
      while (doL) {
        LCounter++;
        DTLayer* layer = buildLayer(fv, sl, type, muonConstants);
        theGeometry->add(layer);

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

      fv.parent();
      doSL = fv.nextSibling(); // go to next SL
    } // sls
    theGeometry->add(chamber);

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

DTLayer * DTGeometryBuilderFromDDD::buildLayer ( DDFilteredView &  fv, DTSuperLayer *  sl, const std::string &  type, const MuonDDDConstants &  muonConstants  ) const

private

create the layer

Definition at line 173 of file DTGeometryBuilderFromDDD.cc.

References DTSuperLayer::add(), DDFilteredView::copyno(), DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), DTNumberingScheme::getDetId(), DDFilteredView::nextSibling(), ecaldqm::topology(), and create_public_lumi_plots::width.

                                                                       {

  MuonDDDNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.getDetId(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
  DTLayerId layId(rawid);

  // Layer specific parameter (size)
  vector<double> par = extractParameters(fv);
  float width = par[0]/cm;     // r-phi  dimension - changes in different chambers
  float length = par[1]/cm;    // z      dimension - constant 126.8 cm
  float thickness = par[2]/cm; // radial thickness - almost constant about 20 cm

  RCPPlane surf(plane(fv, new RectangularPlaneBounds(width, length, thickness) ));

  // Loop on wires
  bool doWire = fv.firstChild();
  int WCounter=0;
  int firstWire=fv.copyno();
  par = extractParameters(fv);
  float wireLength = par[1]/cm;
  while (doWire) {
    WCounter++;
    doWire = fv.nextSibling(); // next wire
  }
  //int lastWire=fv.copyno();
  DTTopology topology(firstWire, WCounter, wireLength);

  DTLayerType layerType;

  DTLayer* layer = new DTLayer(layId, surf, topology, layerType, sl);

  sl->add(layer);
  return layer;
}

DTSuperLayer * DTGeometryBuilderFromDDD::buildSuperLayer ( DDFilteredView &  fv, DTChamber *  chamber, const std::string &  type, const MuonDDDConstants &  muonConstants  ) const

private

create the SL

Definition at line 142 of file DTGeometryBuilderFromDDD.cc.

References DTChamber::add(), DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), DTNumberingScheme::getDetId(), and create_public_lumi_plots::width.

                                                                         {

  MuonDDDNumbering mdddnum(muonConstants);
  DTNumberingScheme dtnum(muonConstants);
  int rawid = dtnum.getDetId(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
  DTSuperLayerId slId(rawid);

  // Slayer specific parameter (size)
  vector<double> par = extractParameters(fv);

  float width = par[0]/cm;     // r-phi  dimension - changes in different chambers
  float length = par[1]/cm;    // z      dimension - constant 126.8 cm
  float thickness = par[2]/cm; // radial thickness - almost constant about 20 cm

  // Ok this is the slayer position...
  RCPPlane surf(plane(fv, new RectangularPlaneBounds(width, length, thickness) ));

  DTSuperLayer* slayer = new DTSuperLayer(slId, surf, chamber);

  //LocalPoint lpos(10,20,30);
  //GlobalPoint gpos=slayer->toGlobal(lpos);

  // add to the chamber
  chamber->add(slayer);
  return slayer;
}

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

private

get parameter also for boolean solid.

Definition at line 213 of file DTGeometryBuilderFromDDD.cc.

References funct::A, ddbox, DDFilteredView::logicalPart(), DDSolid::parameters(), DDSolid::shape(), DDLogicalPart::solid(), and DDBooleanSolid::solidA().

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

DTGeometryBuilderFromDDD::RCPPlane DTGeometryBuilderFromDDD::plane ( const DDFilteredView &  fv, Bounds *  bounds  ) const

private

Definition at line 230 of file DTGeometryBuilderFromDDD.cc.

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

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

  const Surface::PositionType posResult(float(trans.x()/cm), 
                                        float(trans.y()/cm), 
                                        float(trans.z()/cm));
  // now the rotation
  //  DDRotationMatrix tmp = fv.rotation();
  // === DDD uses 'active' rotations - see CLHEP user guide ===
  //     ORCA uses 'passive' rotation. 
  //     'active' and 'passive' rotations are inverse to each other
  //  DDRotationMatrix tmp = fv.rotation();
  DDRotationMatrix rotation = fv.rotation();//REMOVED .Inverse();
  DD3Vector x, y, z;
  rotation.GetComponents(x,y,z);
//   std::cout << "INVERSE rotation by its own operator: "<< fv.rotation() << std::endl;
//   std::cout << "INVERSE rotation manually: "
//      << x.X() << ", " << x.Y() << ", " << x.Z() << std::endl
//      << y.X() << ", " << y.Y() << ", " << y.Z() << std::endl
//      << z.X() << ", " << z.Y() << ", " << z.Z() << std::endl;

  Surface::RotationType rotResult(float(x.X()),float(x.Y()),float(x.Z()),
                                  float(y.X()),float(y.Y()),float(y.Z()),
                                  float(z.X()),float(z.Y()),float(z.Z())); 

//   std::cout << "rotation by its own operator: "<< tmp << std::endl;
//   DD3Vector tx, ty,tz;
//   tmp.GetComponents(tx, ty, tz);
//   std::cout << "rotation manually: "
//      << tx.X() << ", " << tx.Y() << ", " << tx.Z() << std::endl
//      << ty.X() << ", " << ty.Y() << ", " << ty.Z() << std::endl
//      << tz.X() << ", " << tz.Y() << ", " << tz.Z() << std::endl;

  return RCPPlane( new Plane( posResult, rotResult, bounds));
}