GEMGeometryBuilderFromDDD Class Reference

#include <GEMGeometryBuilderFromDDD.h>

Public Member Functions
GEMGeometry *	build (const DDCompactView *cview, const MuonDDDConstants &muonConstants)
	GEMGeometryBuilderFromDDD ()
	~GEMGeometryBuilderFromDDD ()

Private Types
typedef ReferenceCountingPointer< BoundPlane >	RCPBoundPlane

Private Member Functions
RCPBoundPlane	boundPlane (const DDFilteredView &fv, Bounds *bounds, bool isOddChamber) const
GEMChamber *	buildChamber (DDFilteredView &fv, GEMDetId detId) const
GEMEtaPartition *	buildEtaPartition (DDFilteredView &fv, GEMDetId detId) const
GEMGeometry *	buildGeometry (DDFilteredView &fview, const MuonDDDConstants &muonConstants)
GEMSuperChamber *	buildSuperChamber (DDFilteredView &fv, GEMDetId detId) const

Private Attributes
std::map< GEMDetId, std::vector< GEMDetId > >	chids

Detailed Description

Build the GEMGeometry ftom the DDD description

Author

M. Maggi - INFN Bari

Definition at line 25 of file GEMGeometryBuilderFromDDD.h.

Member Typedef Documentation

typedef ReferenceCountingPointer<BoundPlane> GEMGeometryBuilderFromDDD::RCPBoundPlane

private

Definition at line 39 of file GEMGeometryBuilderFromDDD.h.

Constructor & Destructor Documentation

GEMGeometryBuilderFromDDD::GEMGeometryBuilderFromDDD

(

)

Implementation of the GEM Geometry Builder from DDD

Author

Port of: MuDDDGEMBuilder (ORCA)

M. Maggi - INFN Bari

Definition at line 28 of file GEMGeometryBuilderFromDDD.cc.

{ }

GEMGeometryBuilderFromDDD::~GEMGeometryBuilderFromDDD

(

)

Definition at line 31 of file GEMGeometryBuilderFromDDD.cc.

{ }

Member Function Documentation

GEMGeometryBuilderFromDDD::RCPBoundPlane GEMGeometryBuilderFromDDD::boundPlane ( const DDFilteredView &  fv, Bounds *  bounds, bool  isOddChamber  ) const

private

Definition at line 238 of file GEMGeometryBuilderFromDDD.cc.

References TkRotation< T >::rotateAxes(), idealTransformation::rotation, DDFilteredView::rotation(), DDFilteredView::translation(), x, y, and z.

Referenced by buildChamber(), buildEtaPartition(), and buildSuperChamber().

                                                               {
  // 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);
  // LogDebug("GEMGeometryBuilderFromDDD") << "translation: "<< fv.translation() << std::endl;
  // LogDebug("GEMGeometryBuilderFromDDD") << "rotation   : "<< fv.rotation() << std::endl;
  // LogDebug("GEMGeometryBuilderFromDDD") << "INVERSE rotation manually: \n"
  //        << 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()));
  
  //Change of axes for the forward
  Basic3DVector<float> newX(1.,0.,0.);
  Basic3DVector<float> newY(0.,0.,1.);
  Basic3DVector<float> newZ(0.,1.,0.);
  // Odd chambers are inverted in gem.xml
  if (isOddChamber) newY *= -1;
  
  rotResult.rotateAxes(newX, newY, newZ);

  return RCPBoundPlane( new BoundPlane( posResult, rotResult, bounds));
}

GEMGeometry * GEMGeometryBuilderFromDDD::build	(	const DDCompactView *	cview,
		const MuonDDDConstants &	muonConstants
	)

Definition at line 34 of file GEMGeometryBuilderFromDDD.cc.

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

Referenced by GEMGeometryESModule::produce().

{
  std::string attribute = "MuStructure"; // "ReadOutName"; // could come from .orcarc
  std::string value     = "MuonEndCapGEM"; // "MuonGEMHits"; // could come from .orcarc

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

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

GEMChamber * GEMGeometryBuilderFromDDD::buildChamber ( DDFilteredView &  fv, GEMDetId  detId  ) const

private

Definition at line 173 of file GEMGeometryBuilderFromDDD.cc.

References boundPlane(), relativeConstraints::chamber, GEMDetId::chamber(), GEMDetId::chamberId(), allConversions_cfi::dz, LogDebug, DDFilteredView::logicalPart(), DDName::name(), DDBase< N, C >::name(), DDSolid::parameters(), DDLogicalPart::solid(), DDBooleanSolid::solidA(), and DDBooleanSolid::solidB().

Referenced by buildGeometry().

                                                                                            {
  LogDebug("GEMGeometryBuilderFromDDD") << "buildChamber "<<fv.logicalPart().name().name() <<" "<< detId <<std::endl;
  
  DDBooleanSolid solid = (DDBooleanSolid)(fv.logicalPart().solid());
  std::vector<double> dpar = solid.solidA().parameters();
  
  double dy = dpar[0]/cm;//length is along local Y
  double dz = dpar[3]/cm;// thickness is long local Z
  double dx1= dpar[4]/cm;// bottom width is along local X
  double dx2= dpar[8]/cm;// top width is along local X
  dpar = solid.solidB().parameters();
  dz += dpar[3]/cm;// chamber thickness
  
  bool isOdd = detId.chamber()%2;
  RCPBoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(dx1,dx2,dy,dz), isOdd ));
  LogDebug("GEMGeometryBuilderFromDDD") << "size "<< dx1 << " " << dx2 << " " << dy << " " << dz <<std::endl;
  
  GEMChamber* chamber = new GEMChamber(detId.chamberId(), surf);
  return chamber;
}

GEMEtaPartition * GEMGeometryBuilderFromDDD::buildEtaPartition ( DDFilteredView &  fv, GEMDetId  detId  ) const

private

Definition at line 194 of file GEMGeometryBuilderFromDDD.cc.

References boundPlane(), GEMDetId::chamber(), DDfetch(), DDValue::doubles(), GeomDetEnumerators::GEM, LogDebug, DDFilteredView::logicalPart(), DDName::name(), dataset::name, DDBase< N, C >::name(), DDSolid::parameters(), DDLogicalPart::solid(), DDFilteredView::specifics(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by buildGeometry().

                                                                                                      {
  LogDebug("GEMGeometryBuilderFromDDD") << "buildEtaPartition "<<fv.logicalPart().name().name() <<" "<< detId <<std::endl;
  
  // EtaPartition specific parameter (nstrips and npads) 
  DDValue numbOfStrips("nStrips");
  DDValue numbOfPads("nPads");
  std::vector<const DDsvalues_type* > specs(fv.specifics());
  std::vector<const DDsvalues_type* >::iterator is = specs.begin();
  double nStrips = 0., nPads = 0.;
  for (;is != specs.end(); is++){
    if (DDfetch( *is, numbOfStrips)) nStrips = numbOfStrips.doubles()[0];
    if (DDfetch( *is, numbOfPads))   nPads = numbOfPads.doubles()[0];
  }
  LogDebug("GEMGeometryBuilderFromDDD") 
    << ((nStrips == 0. ) ? ("No nStrips found!!") : ("Number of strips: " + boost::lexical_cast<std::string>(nStrips))); 
  LogDebug("GEMGeometryBuilderFromDDD") 
    << ((nPads == 0. ) ? ("No nPads found!!") : ("Number of pads: " + boost::lexical_cast<std::string>(nPads)));
  
  // EtaPartition specific parameter (size) 
  std::vector<double> dpar = fv.logicalPart().solid().parameters();

  double be = dpar[4]/cm; // half bottom edge
  double te = dpar[8]/cm; // half top edge
  double ap = dpar[0]/cm; // half apothem
  double ti = 0.4/cm;     // half thickness
  
  std::vector<float> pars;
  pars.push_back(be); 
  pars.push_back(te); 
  pars.push_back(ap); 
  pars.push_back(nStrips);
  pars.push_back(nPads);
  
  bool isOdd = detId.chamber()%2;
  RCPBoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(be, te, ap, ti), isOdd ));
  std::string name = fv.logicalPart().name().name();
  GEMEtaPartitionSpecs* e_p_specs = new GEMEtaPartitionSpecs(GeomDetEnumerators::GEM, name, pars);
  
  LogDebug("GEMGeometryBuilderFromDDD") << "size "<< be << " " << te << " " << ap << " " << ti <<std::endl;
  GEMEtaPartition* etaPartition = new GEMEtaPartition(detId, surf, e_p_specs);
  return etaPartition;
}

GEMGeometry * GEMGeometryBuilderFromDDD::buildGeometry ( DDFilteredView &  fview, const MuonDDDConstants &  muonConstants  )

private

Definition at line 47 of file GEMGeometryBuilderFromDDD.cc.

References GEMRing::add(), GEMStation::add(), GEMRegion::add(), GEMChamber::add(), GEMSuperChamber::add(), GEMGeometry::add(), GEMNumberingScheme::baseNumberToUnitNumber(), buildChamber(), buildEtaPartition(), buildSuperChamber(), GEMGeometry::chamber(), DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), geometry, GEMSuperChamber::id(), GEMDetId::layer(), LogDebug, DDFilteredView::logicalPart(), GEMDetId::maxStationId, DDName::name(), dataset::name, DDBase< N, C >::name(), DDFilteredView::nextSibling(), DDFilteredView::parent(), relativeConstraints::ring, GEMStation::setName(), Validation_hcalonly_cfi::sign, relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, and GEMGeometry::superChambers().

Referenced by build().

{
  LogDebug("GEMGeometryBuilderFromDDD") <<"Building the geometry service";
  GEMGeometry* geometry = new GEMGeometry();
  LogDebug("GEMGeometryBuilderFromDDD") << "About to run through the GEM structure\n" 
                    <<" First logical part "
                    <<fv.logicalPart().name().name();
 

  
  bool doSuper = fv.firstChild();
  LogDebug("GEMGeometryBuilderFromDDD") << "doSuperChamber = " << doSuper;
  // loop over superchambers
  while (doSuper){

    // getting chamber id from eta partitions
    fv.firstChild();fv.firstChild();
    MuonDDDNumbering mdddnumCh(muonConstants);
    GEMNumberingScheme gemNumCh(muonConstants);
    int rawidCh = gemNumCh.baseNumberToUnitNumber(mdddnumCh.geoHistoryToBaseNumber(fv.geoHistory()));
    GEMDetId detIdCh = GEMDetId(rawidCh);
    // back to chambers
    fv.parent();fv.parent();

    // currently there is no superchamber in the geometry
    // only 2 chambers are present separated by a gap.
    // making superchamber out of the first chamber layer including the gap between chambers
    if (detIdCh.layer() == 1){// only make superChambers when doing layer 1
      GEMSuperChamber *gemSuperChamber = buildSuperChamber(fv, detIdCh);
      geometry->add(gemSuperChamber);
    }
    GEMChamber *gemChamber = buildChamber(fv, detIdCh);
    
    // loop over chambers
    // only 1 chamber
    bool doChambers = fv.firstChild();
    while (doChambers){
      
    // loop over GEMEtaPartitions
      bool doEtaPart = fv.firstChild();
      while (doEtaPart){

    MuonDDDNumbering mdddnum(muonConstants);
    GEMNumberingScheme gemNum(muonConstants);
    int rawid = gemNum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
    GEMDetId detId = GEMDetId(rawid);

    GEMEtaPartition *etaPart = buildEtaPartition(fv, detId);
    gemChamber->add(etaPart);
    geometry->add(etaPart);
    doEtaPart = fv.nextSibling();
      }
      fv.parent();

      geometry->add(gemChamber);
      
      doChambers = fv.nextSibling();
    }
    fv.parent();

    doSuper = fv.nextSibling();
  }
  
  auto& superChambers(geometry->superChambers());
  // construct the regions, stations and rings. 
  for (int re = -1; re <= 1; re = re+2) {
    GEMRegion* region = new GEMRegion(re);
    for (int st=1; st<=GEMDetId::maxStationId; ++st) {
      GEMStation* station = new GEMStation(re, st);
      std::string sign( re==-1 ? "-" : "");
      std::string name("GE" + sign + std::to_string(st) + "/1");
      station->setName(name);
      for (int ri=1; ri<=1; ++ri) {
    GEMRing* ring = new GEMRing(re, st, ri);
    for (unsigned sch=0; sch<superChambers.size(); ++sch){
      GEMSuperChamber* superChamber = const_cast<GEMSuperChamber*>(superChambers.at(sch));
      const GEMDetId detId(superChamber->id());
      if (detId.region() != re || detId.station() != st || detId.ring() != ri) continue;
      
      auto ch1 = geometry->chamber(GEMDetId(detId.region(),detId.ring(),detId.station(),1,detId.chamber(),0));
      auto ch2 = geometry->chamber(GEMDetId(detId.region(),detId.ring(),detId.station(),2,detId.chamber(),0));
      superChamber->add(const_cast<GEMChamber*>(ch1));
      superChamber->add(const_cast<GEMChamber*>(ch2));
      
      ring->add(superChamber);
      LogDebug("GEMGeometryBuilderFromDDD") << "Adding super chamber " << detId << " to ring: " 
                        << "re " << re << " st " << st << " ri " << ri << std::endl;
    }
    LogDebug("GEMGeometryBuilderFromDDD") << "Adding ring " <<  ri << " to station " << "re " << re << " st " << st << std::endl;
    station->add(const_cast<GEMRing*>(ring));
    geometry->add(const_cast<GEMRing*>(ring));
      }
      LogDebug("GEMGeometryBuilderFromDDD") << "Adding station " << st << " to region " << re << std::endl;
      region->add(const_cast<GEMStation*>(station));
      geometry->add(const_cast<GEMStation*>(station));
    }
    LogDebug("GEMGeometryBuilderFromDDD") << "Adding region " << re << " to the geometry " << std::endl;
    geometry->add(const_cast<GEMRegion*>(region));
  }
  
  return geometry;
}

GEMSuperChamber * GEMGeometryBuilderFromDDD::buildSuperChamber ( DDFilteredView &  fv, GEMDetId  detId  ) const

private

Definition at line 150 of file GEMGeometryBuilderFromDDD.cc.

References boundPlane(), GEMDetId::chamber(), allConversions_cfi::dz, LogDebug, DDFilteredView::logicalPart(), DDName::name(), DDBase< N, C >::name(), DDSolid::parameters(), DDLogicalPart::solid(), DDBooleanSolid::solidA(), DDBooleanSolid::solidB(), and GEMDetId::superChamberId().

Referenced by buildGeometry().

                                                                                                      {
  LogDebug("GEMGeometryBuilderFromDDD") << "buildSuperChamber "<<fv.logicalPart().name().name() <<" "<< detId <<std::endl;
  
  DDBooleanSolid solid = (DDBooleanSolid)(fv.logicalPart().solid());
  std::vector<double> dpar = solid.solidA().parameters();
  
  double dy = dpar[0]/cm;//length is along local Y
  double dz = dpar[3]/cm;// thickness is long local Z
  double dx1= dpar[4]/cm;// bottom width is along local X
  double dx2= dpar[8]/cm;// top width is along local X
  dpar = solid.solidB().parameters();
  dz += dpar[3]/cm;// chamber thickness
  dz *=2; // 2 chambers in superchamber
  dz += 2.105;// gap between chambers
  
  bool isOdd = detId.chamber()%2;
  RCPBoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(dx1,dx2,dy,dz), isOdd ));
  LogDebug("GEMGeometryBuilderFromDDD") << "size "<< dx1 << " " << dx2 << " " << dy << " " << dz <<std::endl;
  
  GEMSuperChamber* superChamber = new GEMSuperChamber(detId.superChamberId(), surf);
  return superChamber;
}

Member Data Documentation

std::map<GEMDetId,std::vector<GEMDetId> > GEMGeometryBuilderFromDDD::chids

private

Definition at line 37 of file GEMGeometryBuilderFromDDD.h.