TIBLayerBuilder Class Reference

#include <TIBLayerBuilder.h>

Public Member Functions
TIBLayer *	build (const GeometricDet *aTIBLayer, const TrackerGeometry *theGeomDetGeometry) __attribute__((cold))
void	constructRings (std::vector< const GeometricDet * > &theGeometricRods, std::vector< std::vector< const GeometricDet * > > &innerGeometricDetRings, std::vector< std::vector< const GeometricDet * > > &outerGeometricDetRings)
	TIBLayerBuilder ()

Detailed Description

A concrete builder for TIBLayer

Definition at line 14 of file TIBLayerBuilder.h.

Constructor & Destructor Documentation

TIBLayerBuilder::TIBLayerBuilder ( )

inline

Definition at line 16 of file TIBLayerBuilder.h.

References build(), and constructRings().

{};

Member Function Documentation

TIBLayer * TIBLayerBuilder::build	(	const GeometricDet *	aTIBLayer,
		const TrackerGeometry *	theGeomDetGeometry
	)

Definition at line 10 of file TIBLayerBuilder.cc.

References TIBRingBuilder::build(), GeometricDet::components(), mps_fire::i, and align::TIBLayer.

Referenced by GeometricSearchTrackerBuilder::build(), and TIBLayerBuilder().

                                                                                                         {
  vector<const GeometricDet*> theGeometricRods = aTIBLayer->components();

  vector<vector<const GeometricDet*> > innerGeometricDetRings;
  vector<vector<const GeometricDet*> > outerGeometricDetRings;

  constructRings(theGeometricRods, innerGeometricDetRings, outerGeometricDetRings);

  TIBRingBuilder myRingBuilder;

  vector<const TIBRing*> innerRings;
  vector<const TIBRing*> outerRings;

  for (unsigned int i = 0; i < innerGeometricDetRings.size(); i++) {
    innerRings.push_back(myRingBuilder.build(innerGeometricDetRings[i], theGeomDetGeometry));
    outerRings.push_back(myRingBuilder.build(outerGeometricDetRings[i], theGeomDetGeometry));
  }

  return new TIBLayer(innerRings, outerRings);
}

void TIBLayerBuilder::constructRings	(	std::vector< const GeometricDet * > &	theGeometricRods,
		std::vector< std::vector< const GeometricDet * > > &	innerGeometricDetRings,
		std::vector< std::vector< const GeometricDet * > > &	outerGeometricDetRings
	)

Definition at line 31 of file TIBLayerBuilder.cc.

References mps_fire::i, LogDebug, and edm::max().

Referenced by TIBLayerBuilder().

                                                                                                   {
  double meanPerp = 0;
  for (vector<const GeometricDet*>::const_iterator it = theGeometricRods.begin(); it != theGeometricRods.end(); it++) {
    meanPerp = meanPerp + (*it)->positionBounds().perp();
  }
  meanPerp = meanPerp / theGeometricRods.size();

  vector<const GeometricDet*> theInnerGeometricRods;
  vector<const GeometricDet*> theOuterGeometricRods;

  for (vector<const GeometricDet*>::const_iterator it = theGeometricRods.begin(); it != theGeometricRods.end(); it++) {
    if ((*it)->positionBounds().perp() < meanPerp)
      theInnerGeometricRods.push_back(*it);
    if ((*it)->positionBounds().perp() > meanPerp)
      theOuterGeometricRods.push_back(*it);
  }

  size_t innerLeftRodMaxSize = 0;
  size_t innerRightRodMaxSize = 0;
  size_t outerLeftRodMaxSize = 0;
  size_t outerRightRodMaxSize = 0;

  for (vector<const GeometricDet*>::const_iterator it = theInnerGeometricRods.begin();
       it != theInnerGeometricRods.end();
       it++) {
    if ((*it)->positionBounds().z() < 0)
      innerLeftRodMaxSize = max(innerLeftRodMaxSize, (**it).components().size());
    if ((*it)->positionBounds().z() > 0)
      innerRightRodMaxSize = max(innerRightRodMaxSize, (**it).components().size());
  }

  for (vector<const GeometricDet*>::const_iterator it = theOuterGeometricRods.begin();
       it != theOuterGeometricRods.end();
       it++) {
    if ((*it)->positionBounds().z() < 0)
      outerLeftRodMaxSize = max(outerLeftRodMaxSize, (**it).components().size());
    if ((*it)->positionBounds().z() > 0)
      outerRightRodMaxSize = max(outerRightRodMaxSize, (**it).components().size());
  }

  LogDebug("TkDetLayers") << "innerLeftRodMaxSize: " << innerLeftRodMaxSize;
  LogDebug("TkDetLayers") << "innerRightRodMaxSize: " << innerRightRodMaxSize;

  LogDebug("TkDetLayers") << "outerLeftRodMaxSize: " << outerLeftRodMaxSize;
  LogDebug("TkDetLayers") << "outerRightRodMaxSize: " << outerRightRodMaxSize;

  for (unsigned int i = 0; i < (innerLeftRodMaxSize + innerRightRodMaxSize); i++) {
    innerGeometricDetRings.push_back(vector<const GeometricDet*>());
  }

  for (unsigned int i = 0; i < (outerLeftRodMaxSize + outerRightRodMaxSize); i++) {
    outerGeometricDetRings.push_back(vector<const GeometricDet*>());
  }

  for (unsigned int ringN = 0; ringN < innerLeftRodMaxSize; ringN++) {
    for (vector<const GeometricDet*>::const_iterator it = theInnerGeometricRods.begin();
         it != theInnerGeometricRods.end();
         it++) {
      if ((*it)->positionBounds().z() < 0) {
        if ((**it).components().size() > ringN)
          innerGeometricDetRings[ringN].push_back((**it).components()[ringN]);
      }
    }
  }

  for (unsigned int ringN = 0; ringN < innerRightRodMaxSize; ringN++) {
    for (vector<const GeometricDet*>::const_iterator it = theInnerGeometricRods.begin();
         it != theInnerGeometricRods.end();
         it++) {
      if ((*it)->positionBounds().z() > 0) {
        if ((**it).components().size() > ringN)
          innerGeometricDetRings[innerLeftRodMaxSize + ringN].push_back((**it).components()[ringN]);
      }
    }
  }

  for (unsigned int ringN = 0; ringN < outerLeftRodMaxSize; ringN++) {
    for (vector<const GeometricDet*>::const_iterator it = theOuterGeometricRods.begin();
         it != theOuterGeometricRods.end();
         it++) {
      if ((*it)->positionBounds().z() < 0) {
        if ((**it).components().size() > ringN)
          outerGeometricDetRings[ringN].push_back((**it).components()[ringN]);
      }
    }
  }

  for (unsigned int ringN = 0; ringN < outerRightRodMaxSize; ringN++) {
    for (vector<const GeometricDet*>::const_iterator it = theOuterGeometricRods.begin();
         it != theOuterGeometricRods.end();
         it++) {
      if ((*it)->positionBounds().z() > 0) {
        if ((**it).components().size() > ringN)
          outerGeometricDetRings[outerLeftRodMaxSize + ringN].push_back((**it).components()[ringN]);
      }
    }
  }
}