Phase2EndcapRingBuilder Class Reference

#include <Phase2EndcapRingBuilder.h>

Public Member Functions
Phase2EndcapRing *	build (const GeometricDet *aPhase2EndcapRing, const TrackerGeometry *theGeomDetGeometry, const bool useBrothers=true) __attribute__((cold))
	Phase2EndcapRingBuilder ()

Detailed Description

A concrete builder for Phase2EndcapRing

Definition at line 15 of file Phase2EndcapRingBuilder.h.

Constructor & Destructor Documentation

Phase2EndcapRingBuilder::Phase2EndcapRingBuilder ( )

inline

Definition at line 17 of file Phase2EndcapRingBuilder.h.

{};

Member Function Documentation

Phase2EndcapRing * Phase2EndcapRingBuilder::build	(	const GeometricDet *	aPhase2EndcapRing,
		const TrackerGeometry *	theGeomDetGeometry,
		const bool	useBrothers = true
	)

Definition at line 7 of file Phase2EndcapRingBuilder.cc.

References GeometricDet::components(), TrackerGeometry::idToDet(), LogDebug, and LogTrace.

Referenced by Phase2EndcapLayerBuilder::build(), and Phase2OTBarrelLayerBuilder::build().

{
  vector<const GeometricDet*>  allGeometricDets = aPhase2EndcapRing->components();
  vector<const GeometricDet*>  compGeometricDets;
  LogDebug("TkDetLayers") << "Phase2EndcapRingBuilder with #Modules: " << allGeometricDets.size() << std::endl;

  vector<const GeomDet*> frontGeomDets;
  vector<const GeomDet*> backGeomDets;
  double meanZ = 0;

  if(!useBrothers){
  
    //---- to evaluate meanZ
    for(vector<const GeometricDet*>::const_iterator compGeometricDets=allGeometricDets.begin(); compGeometricDets!=allGeometricDets.end();compGeometricDets++){
        LogTrace("TkDetLayers") << " compGeometricDets->positionBounds().perp() " << (*compGeometricDets)->positionBounds().z() << std::endl;
        meanZ = meanZ + (*compGeometricDets)->positionBounds().z();
    }
    meanZ = meanZ/allGeometricDets.size();
    LogDebug("TkDetLayers") << " meanZ " << meanZ << std::endl;
  
    for(vector<const GeometricDet*>::const_iterator compGeometricDets=allGeometricDets.begin();
        compGeometricDets!=allGeometricDets.end();compGeometricDets++){
      const GeomDet* theGeomDet = theGeomDetGeometry->idToDet( (*compGeometricDets)->geographicalID() );
  
      if( fabs( (*compGeometricDets)->positionBounds().z() ) < fabs(meanZ) )
        frontGeomDets.push_back(theGeomDet);
  
      if( fabs( (*compGeometricDets)->positionBounds().z() ) > fabs(meanZ) )
        backGeomDets.push_back(theGeomDet);      

      if( fabs( (*compGeometricDets)->positionBounds().z() ) == fabs(meanZ) )
        throw DetLayerException("Not possible to assiciate this GeometricDet in front or back");
  
    }
  
    LogDebug("TkDetLayers") << "frontGeomDets.size(): " << frontGeomDets.size() ;
    LogDebug("TkDetLayers") << "backGeomDets.size(): " << backGeomDets.size() ;
  
    return new Phase2EndcapRing(frontGeomDets,backGeomDets);

  } else {

    vector<const GeomDet*> frontGeomDetBrothers;
    vector<const GeomDet*> backGeomDetBrothers;
    vector<const GeometricDet*>  compGeometricDets;
  
    //---- to evaluate meanZ
    double meanZ = 0;
    double meanZBrothers = 0;
    for(vector<const GeometricDet*>::const_iterator it=allGeometricDets.begin(); it!=allGeometricDets.end();it++){
      compGeometricDets = (*it)->components();
      if (compGeometricDets.size() != 2){
        throw DetLayerException("Phase2OTEndcapRing is considered as a stack but does not have two components");
      } else {
        LogTrace("TkDetLayers") << " compGeometricDets[0]->positionBounds().perp() " << compGeometricDets[0]->positionBounds().z() << std::endl;
        LogTrace("TkDetLayers") << " compGeometricDets[1]->positionBounds().perp() " << compGeometricDets[1]->positionBounds().z() << std::endl;
        meanZ = meanZ + compGeometricDets[0]->positionBounds().z();
        meanZBrothers = meanZBrothers + compGeometricDets[1]->positionBounds().z();
      }

    }
    meanZ = meanZ/allGeometricDets.size();
    meanZBrothers = meanZBrothers/allGeometricDets.size();
    LogDebug("TkDetLayers") << " meanZ Lower " << meanZ << std::endl;
    LogDebug("TkDetLayers") << " meanZ Upper " << meanZBrothers << std::endl;
  
    for(vector<const GeometricDet*>::const_iterator it=allGeometricDets.begin();
        it!=allGeometricDets.end();it++){
      compGeometricDets = (*it)->components(); 
      const GeomDet* theGeomDet = theGeomDetGeometry->idToDet( compGeometricDets[0]->geographicalID() );
  
      if( fabs( compGeometricDets[0]->positionBounds().z() ) < fabs(meanZ) )
        frontGeomDets.push_back(theGeomDet);
  
      if( fabs( compGeometricDets[0]->positionBounds().z() ) > fabs(meanZ) )
        backGeomDets.push_back(theGeomDet);      
  
      const GeomDet* theGeomDetBrother = theGeomDetGeometry->idToDet( compGeometricDets[1]->geographicalID() );
  
      if( fabs( compGeometricDets[1]->positionBounds().z() ) < fabs(meanZBrothers) )
        frontGeomDetBrothers.push_back(theGeomDetBrother);
  
      if( fabs( compGeometricDets[1]->positionBounds().z() ) > fabs(meanZBrothers) )
        backGeomDetBrothers.push_back(theGeomDetBrother);

      if( fabs( compGeometricDets[0]->positionBounds().z() ) == fabs(meanZ) || fabs( compGeometricDets[1]->positionBounds().z() ) == fabs(meanZBrothers) )
        throw DetLayerException("Not possible to assiciate components of this GeometricDet in front or back");
    }
  
    LogDebug("TkDetLayers") << "frontGeomDets.size(): " << frontGeomDets.size() ;
    LogDebug("TkDetLayers") << "backGeomDets.size(): " << backGeomDets.size() ;
    LogDebug("TkDetLayers") << "frontGeomDetsBro.size(): " << frontGeomDetBrothers.size() ;
    LogDebug("TkDetLayers") << "backGeomDetsBro.size(): " << backGeomDetBrothers.size() ;
  
    return new Phase2EndcapRing(frontGeomDets,backGeomDets,frontGeomDetBrothers,backGeomDetBrothers);
  }
  return new Phase2EndcapRing(frontGeomDets,backGeomDets);

}