TIBLayerBuilder Class Reference

#include <TIBLayerBuilder.h>

Public Member Functions
TIBLayer *	build (const GeometricDet *aTIBLayer, const TrackerGeometry *theGeomDetGeometry)
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 15 of file TIBLayerBuilder.h.

Constructor & Destructor Documentation

TIBLayerBuilder::TIBLayerBuilder ( )

inline

Definition at line 17 of file TIBLayerBuilder.h.

{};

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(), i, and align::TIBLayer.

Referenced by GeometricSearchTrackerBuilder::build().

{
  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 37 of file TIBLayerBuilder.cc.

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

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